针对非均匀校正后的热红外图像剩余非均匀噪声去除方法
【专利摘要】本发明公开了一种针对非均匀校正后的热红外图像剩余非均匀噪声去除方法。其实施步骤如下：1)利用红外成像系统在不同温度下对黑体成像得到红外图像；2)将得到的黑体图像进行非均匀校正；3)求得校正后黑体图像的像素平均值；4)将校正后的黑体图像中每个像素值与其平均值相减存储结果；5)拍摄外景图像并进行非均匀校正；6)将校正后图像像素值与不同温度下图像平均值做比较找到最接近值对应的相减数组，将像素值与数组中对应位置的值相减得到新的图像。本发明弥补了传统红外图像噪声去除方法不足，提供了一种更有效、更精确去除剩余非均匀空间噪声的方法。本发明能够有效去除红外图像中的剩余非均匀空间噪声，并不损失图像原信息，获得高质量的复原图像。
【专利说明】
针对非均匀校正后的热红外图像剩余非均匀噪声去除方法
技术领域：
[0001] 本发明属于高灵敏度红外成像及应用技术领域，基于一种制冷型红外成像硬件系 统，涉及红外图像非均匀校正后的剩余空间噪声的去除。
【背景技术】：
[0002] 随着红外成像技术的发展，人们对热红外成像系统灵敏度水平的要求越来越高， 获取更高灵敏度红外图像的需求也更加迫切，经过非均匀性校正的红外图像已不能满足高 灵敏度场合的需要，因为即使经过非均匀性校正，红外图像仍然存在剩余非均匀性空间噪 声，在高灵敏度场合，这会极大的影响对信号的分辨能力，必须采取新的措施以进一步提高 红外图像的灵敏度水平。高灵敏度红外成像系统一般采用制冷型红外焦平面器件，由于焦 平面阵列中各探测元响应特性不完全一致，不同的探测光敏元件对相同的红外辐射会产生 不同的输出信号，使得红外图像出现非均匀性噪声。常规的非均匀性校正算法无法完全去 除该噪声，但在一般应用场合该剩余噪声不影响图像质量，所以可见文献中尚无针对非均 匀处理后的红外图像进一步去除非均匀性噪声的方法。本方法基于此问题而提出，基于红 外焦平面硬件特性，针对性的进一步削弱非均匀性噪声，可以在不损失图像信息的前提下 有效抑制剩余非均匀性噪声。
【发明内容】
：
[0003] 发明目的:本发明提供了一种有效、精确去除剩余非均匀空间噪声的全新方法。这 种方法根据像元的响应特性以及在定标图像校正后像元与理想值偏差的固定性出发，创造 性的提出了针对剩余非均匀性空间噪声的去除方法：
[0004] 本发明提出的红外成像系统剩余非均匀空间噪声去除方法，包括以下步骤：
[0005] 1)利用基于MXN规模的焦平面阵列的红外成像系统，分别在A，A+e，A+2e，……B的 温度下对均匀的黑体进行成像，得到Η幅图像……Yh，A为设定的起始温度值，B为设定 的终止温度值，e为温度间隔值；
[0006] 2)选取(：1、〇2、(：3、〇4四个温度下的四幅黑体图像? 1、？2、？3^4作为进行非均匀校正的 标准图像，用这四幅图像分别对步骤1)中所得到的Η幅黑体图像进行非均匀校正，得到Η幅 校正后的黑体图像Xi、X 2……Xh;非均匀性校正所用系数由如下公式求得：
[0009]其中，为取定的定标点，Vi,j( (H)为测得的第i行，第j列个光敏元对应的响应输 出；)对红外焦平面阵列所有光敏元的输出信号求平均得；
[00?0] 3)对步骤2)中非均勾校正后的每幅图像求像素点平均值，计为averi，aver2,...... averH，将计算得到的平均值结果作为每幅黑体图像的标准值：
[0012]其中X(i，j)是图像第i行，第j列元素的像素值；
[0013] 4)将步骤2)中非均匀校正后的每幅图像中的每个像素点值与步骤3)中求得的对 应图像的平均值求差，得到Η个Μ X N的二维数组Di，D2……Dh，每个数组中存储着对应图像求 得的差值；
[0014] 5)拍摄外景图像得到红外图像f，用步骤2)中选取的四幅黑体图像 图像f进行非均匀校正，得到校正后的图像F;
[0015] 6)将校正后图像F中的每个像素点值F(i，j)分别与步骤3)中得到的平均值进行比 较，找到与像素点的值最接近的平均值averx，找到最接近的平均值后，再找出与该平均值 所对应的步骤4)中求出的差值二维矩阵Dx，其中X为1~Η中的值；
[0016] 7)将图像F中的像素点F( i，j)与步骤6)中找到的Dx数组中与该像素位置对应的点 Dx(i，j)值求差，将求得的结果作为新的像素值，得到最终去除空间噪声后的红外图像。
[0017] 本发明相对于现有方法具有如下的优点及效果：
[0018] 1本发明紧接传统非均匀校正方法之后，创造性的进一步降低了剩余非均匀性噪 声，为降低红外图像空间噪声提供了新的手段。
[0019] 2本发明针对传统非均匀校正方法处理后的红外图像的剩余非均匀性噪声，在实 验室以不同像元在不同辐射量于固定定标点校正后对均值的偏差作为补偿，降低了在高灵 敏度运用场合剩余非均匀性噪声对图像的干扰，降低了空间噪声。，
[0020] 3本发明针对非均匀性噪声的固定性得到的偏差补偿数据集，可广泛用于各种场 景对固定图案噪声的消除，提高图像的非均匀性，无场景使用限制。
【附图说明】：
[0021] 图1是本发明的红外成像系统剩余非均匀的空间噪声去除方法框图。
[0022] 图2是本发明中（10,40)位置上的黑体图像像素在各温度下与其均值求得的差值 分布曲线。
[0023]图3(a，b)是本发明实施例中非均匀性校正后经过直方图2/255拉伸后的黑体图像 (a) 与用本文方法降低剩余非均匀空间噪声的同样按照2/255同样比例拉伸后的黑体图像 (b) 的对比。
[0024]图4(a，b)是本发明实施例中待去噪外景图像(a)与经过处理后的外景图像(b)的 对比。
【具体实施方式】：
[0025]下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限 于此。
[0026]如图1，以320X256焦平面为硬件的红外成像系统为例，改进的红外系统剩余非均 匀空间噪声去除方法，包括以下步骤：
[0027] 1)利用红外成像系统，分别在6,8，10，12，14，15，17,20,22,24,25,27摄氏度下对 均匀黑体进行成像，得到12幅黑体图像Y^Ys、……Y12，温度的间隔不等。
[0028] 2)选取10、15、20、25摄氏度四个温度下的四幅黑体图像^、？2、内、？ 4作为进行非均 匀校正的标准图像。用这四幅图像分别对1)中所得到的12幅黑体图像进行非均匀校正，得 到12幅校正后的黑体图像。
[0029] 3)对2)中非均勾校正后的每幅图像求像素点平均值，计为averi = 36108，aver2 = 3737，aver3 = 38853，aver4 = 40313."aver7 = 44174."averi〇 = 49857 ???avemz 52449。将计 算得到的平均值结果作为每幅黑体图像的标准值。
[0030] 4)将2)中非均匀校正后的每幅图像中的每个像素点值与3)中求得的对应图像的 平均值求差。得到12个256 X 320的二维数组D!，D2……D46，每个数组中存储着对应图像求得 的差值，图2所示是各温度下（10,40)位置的像素与各平均值求得的差。
[0031] 5)拍摄外景图像得到红外图像f，用2)中选取的四幅图像对图像f进行非均匀校 正，得到校正后的图像F，如图4(a)所示。
[0032] 6)将校正后图像F中的每个像素点值F(i，j)分别与3)中得到的平均值进行比较， 找到与像素点的值最接近的平均值averx，其中i的值为1~256，j的值为1~320，找到最接 近的平均值后，再找出与该平均值所对应的4)中求出的差值二维矩阵Dx，其中X为1~12中 的值。
[0033] 7)将图像F中的像素点F(i，j)与6)中找到的Dx数组中与该像素位置对应的点Dx(i， j)值求差，将求得的结果作为新的像素值，得到最终去除空间噪声后的红外图像。
[0034]如图3(a)所示，是没有经过本文方法处理的非均匀性校正拉伸后黑体图像，与原 始图像拉伸比例为2/255。其中可以明显的看到存在着剩余非均匀性空间噪声，剩余非均匀 性为0.044052 %。图3 (b)为经过本文方法处理后的图像，可以看到剩余非均匀性噪声降低 为0.039722%。
[0035] 用本文的去除剩余非均匀性空间噪声的方法对外景图像进行降噪，如图4(b)所示 为本实施例的去剩余非均匀性空间噪声结果，将图4(b)与图4(a)中的原始空间噪声比较， 可见，其中的空间固定噪声被大幅降低，大幅提高了图像质量，凸显了弱信号目标。
[0036] 以上所述仅为本发明的优选实施方式，本发明的保护范围并不仅限于上述实施方 式，凡是属于本发明原理的技术方案均属于本发明的保护范围。对于本领域的技术人员而 言，在不脱离本发明的原理的前提下进行的若干改进和润饰，这些改进也应视为本发明的 保护范围。
【主权项】
1. 一种针对非均匀校正后的热红外图像剩余非均匀噪声去除方法，包括如下步骤： 1) 利用基于MXN规模的焦平面阵列的红外成像系统，分别在A，A+e，A+2e，……B的温度 下对均匀的黑体进行成像，得到Η幅图像……Yh，A为设定的起始温度值，B为设定的终 止温度值，e为温度间隔值； 2) 选取(^、(^(^(^四个温度下的四幅黑体图像心心内心作为进行非均匀校正的标准 图像，用这四幅图像分别对步骤1)中所得到的Η幅黑体图像进行非均匀校正，得到Η幅校正 后的黑体图像Xi、X 2……Xh;非均匀性校正所用系数由如下公式求得：其中，Φ k为取定的定标点，Vi, j ( Φ k)为测得的第i行，第j列个光敏元对应的响应输出； 内舍)对红外焦平面阵列所有光敏元的输出信号求平均得； 3) 对步骤2)中非均勾校正后的每幅图像求像素点平均值，计为averi，aver2，...... averH，将计算得到的平均值结果作为毎幅黑体图像的标准倌：其中X(i，j)是图像第i行，第j列元素的像素值； 4) 将步骤2)中非均匀校正后的每幅图像中的每个像素点值与步骤3)中求得的对应图 像的平均值求差，得到Η个MXN的二维数组……Dh，每个数组中存储着对应图像求得的 差值； 5) 拍摄外景图像得到红外图像f，用步骤2)中选取的四幅黑体图像&、^内力对图像€ 进行非均匀校正，得到校正后的图像F; 6) 将校正后图像F中的每个像素点值F(i，j)分别与步骤3)中得到的平均值进行比较， 找到与像素点的值最接近的平均值averx，找到最接近的平均值后，再找出与该平均值所对 应的步骤4)中求出的差值二维矩阵Dx，其中X为1~Η中的值； 7) 将图像F中的像素点F(i，j)与步骤6)中找到的Dx数组中与该像素位置对应的点Dx(i， j)值求差，将求得的结果作为新的像素值，得到最终去除空间噪声后的红外图像。
【文档编号】G06T5/00GK105869129SQ201610235714
【公开日】2016年8月17日
【申请日】2016年4月15日
【发明人】姬弘桢, 王建宇, 李春来, 孙羽, 金健, 陈小文
【申请人】中国科学院上海技术物理研究所