红外图像非均匀性矫正两点法实现系统

资 源 简 介

本项目旨在解决红外焦平面阵列（IRFPA）在制造过程中因工艺差异导致的各像素点响应不一致问题，即非均匀性噪声。本系统采用经典的物理定标两点法，通过搜集红外探测器在低温和高温两个已知标定黑体辐射源下的原始响应数据，建立像素级的线性响应模型。功能涵盖了定标图像的多帧累加平均，以有效滤除读出电路产生的随机噪声并获取稳定的基准均值；随后通过计算每个像素点相对于全局均值的斜率和截距，生成精确的增益校正矩阵和偏置校正矩阵。在校正阶段，系统将上述矩阵应用于待处理的原始红外图像，通过实时的矩阵线性变换消除图像中的固定图形

详 情 说 明

红外图像非均匀性矫正两点法实现系统

项目介绍

本系统针对红外焦平面阵列（IRFPA）中各像素由于制造工艺差异导致的响应不一致问题，实现了基于物理定标的经典两点法非均匀性矫正算法。系统通过模拟物体在不同温度下的热辐射响应，建立像素级的线性校正模型，能够有效消除红外图像中的固定图形噪声（FPN）和条纹，显著提升红外成像质量。该方案不仅包含了核心算法的实现，还集成了完整的仿真环境、量化评估模型以及可视化分析工具。

功能特性

• 高精度模拟仿真：支持模拟真实的探测器非均匀性参数（增益偏离与偏置偏差）以及读取电路的随机高斯噪声。
• 物理定标流程实现：完整模拟了低温黑体与高温黑体两个标定点的辐射数据获取逻辑。
• 多帧降噪处理：通过对原始定标数据进行多帧累加平均，大幅削弱随时间变化的随机噪声对定标精度的干扰。
• 矩阵化高效运算：直接利用MATLAB的矩阵并行化处理能力计算增益和偏置矩阵，实现像素级的实时校正模拟。
• 综合量化评估：内置非均匀性（NU）指标及信噪比（SNR）分析功能，直观展示矫正前后的性能提升。
• 多维度数据可视化：提供原始图像、矫正图像、校正矩阵分布图、灰度直方图以及自动化评估报告。

系统要求

• 运行环境：MATLAB R2016b 或更高版本。
• 必备工具箱：Image Processing Toolbox (用于加载示例场景图像及基础图像处理操作)。

实现逻辑与算法分析

系统主程序的实现过程严格遵循红外定标的工程步骤：

1. 模拟环境构建

首先定义图像的分辨率（256x256）与累加帧数。随后通过正态分布随机数生成模拟探测器的“真实增益分布”和“真实偏置分布”，并设定低温与高温两个定标源的辐射强度值。

2. 原始数据采集与降噪

针对高、低温两个参考源，分别生成多帧（本系统中为50帧）含有随机噪声的原始响应序列。系统通过计算三维矩阵在时间轴上的均值，获取稳定的低温均值帧和高温均值帧。这一步是确保定标参数准确性的核心，能有效降低读出电路噪声的影响。

3. 两点定标参数计算

系统基于线性响应模型 $Y = kX + b$ 求解每个像素的校正系数：

• 全局期望值计算：计算所有像素对于高、低温源的全局平均响应，作为校正后的目标值。
• 增益校正矩阵（Gain Map）：通过计算全局温差期望与每个像素实际温差响应的比值来确定。为了防止分母为零导致的异常，程序加入了极小值补偿处理。
• 偏置校正矩阵（Offset Map）：依据求得的增益系数，结合低温点的参考均值，计算出每个像素的平移量。

4. 模拟探测器输出与校正应用

加载一张标准图像（如cameraman）作为真实场景，将其叠加预设的非均匀性模型和随机噪声，模拟出真实的探测器输出图像。随后将该图像与计算得到的增益、偏置矩阵进行线性组合运算，还原出接近真实的场景信息。

5. 性能评估与结果呈现

• 非均匀性（NU）计算：通过统计均匀测试区域的像素标准差与均值之比，量化评估固定图形噪声的消除程度。
• 信噪比（SNR）分析：通过计算均方误差（MSE），量化图像在矫正过程中的信号质量提升。
• 可视化界面：程序最后生成一个包含六个子图的综合看板，不仅直观对照矫正前后的视觉效果，还展示了校正系数的空间分布和像素灰度分布的演变。

使用方法

1. 启动MATLAB软件。
2. 将包含系统脚本的文件夹设置为当前工作路径。
3. 在命令行窗口输入主程序函数名并回车。
4. 程序将自动生成包含50帧定标数据的模拟环境，完成计算并弹出包含量化报告的可视化窗口。
5. 在控制台窗口可查看具体的非均匀性数值下降比例及信噪比提升分贝（dB）。