一种基于双波段的红外图像盲元补偿方法与流程

1.本发明属于双波段红外图像处理技术领域，具体涉及一种基于双波段的红外图像盲元补偿方法。


背景技术：

2.红外焦平面阵列由于受到制造工艺等的影响，常常会出现盲元，盲元包括固定盲元和随机盲元，其表现形式通常是在图像上形成暗点或者亮点。盲元会影响探测器成像的质量，因此需要对焦平面阵列进行盲元检测和补偿。
3.双波段红外焦平面阵列包括中长波段、中短波段、长短波段等，相对于传统的单波段红外焦平面阵列具有较大的优势，在环境条件极其复杂多变的现代战场上，双波段红外焦平面阵列可以获取更多的目标信息，提高目标的检测和识别机率。但是双波段红外焦平面阵列依然存在盲元的困扰。。
4.此前，国内对盲元校正的研究通常是在单波段层面，基本没有基于双波段的盲元补偿算法相关的研究，基于单波段的盲元补偿算法由于缺失了真实的场景信息，对于盲元的补偿是通过利用邻域的场景信息来对盲元进行估计补偿，对大盲元簇(如7
×
7以上盲元簇)及位于边缘位置的盲元补偿效果不尽如人意。
5.因此如何克服现有技术的不足是目前双波段红外图像处理技术领域亟需解决的问题

技术实现要素：

6.本发明的目的是为了解决现有技术的不足，提供一种基于双波段的红外图像盲元补偿方法，利用双波段背景辐射信息的差异，可以对任一单个波段中的盲元进行更好的补偿，从而改善了现有双波段红外图像中的盲元现象。
7.为实现上述目的，本发明采用的技术方案如下：
8.一种基于双波段的红外图像盲元补偿方法，包括如下步骤：
9.步骤s1，对输入的双波段红外图像xi(i，j)、yi(i，j)分别进行盲元检测；
10.步骤s2，对双波段红外图像进行边缘提取，提取边缘图像ex(i，j)、ey(i，j)；
11.步骤s3，根据步骤s1检测情况及步骤s2提取情况，对双波段红外图像的像元情况进行分类，并根据分类情况对像元进行补偿后输出图像。
12.进一步，优选的是，所述的双波段红外图像包括中长波段、中短波段、短长波段、中中波段图像。
13.进一步，优选的是，所述的基于双波段的红外图像盲元补偿方法，包括如下步骤：
14.步骤s1，对输入的双波段红外图像xi(i，j)、yi(i，j)分别进行盲元检测，设立阈值h，初始盲元表置1，若某一波段像元(i，j)连续h帧被判断为盲元，则对应波段盲元表bpx(i，j)或bpy(i，j)置0；
15.步骤s2，对双波段红外图像进行边缘提取，提取边缘图像ex(i，j)、ey(i，j)，对边缘部分置0，非边缘部分置1；
16.步骤s3，根据步骤s1检测情况及步骤s2提取情况，对双波段红外图像的像元情况
进行分类，并根据分类情况对像元进行补偿后输出图像。
17.0和1只是对盲元及边缘状态的一个标记，0代表是盲元和边缘，1代表不是盲元和边缘，便于补偿时进行相应的识别计算。
18.进一步，优选的是，步骤s2中，使用sobel算子进行边缘提取。
19.进一步，优选的是，步骤s3中，对双波段红外图像的像元情况进行分类，具体为：
20.(1)像元在一个波段为盲元，而在另一个波段非盲元，且此像元未位于边缘；
21.(2)像元在一个波段为盲元，而在另一个波段非盲元，且此像元位于边缘；
22.(3)像元在两个波段都为盲元，且此像元未位于边缘；
23.(4)像元在两个波段都为盲元，且此像元位于边缘。
24.进一步，优选的是，步骤s3中，根据分类情况，采取不同的策略对像元进行补偿，当像元在两个波段都为盲元时，采取单波段的补偿策略对盲元进行补偿，如(1)，(2)所示；当像元仅在单波段为盲元，而另一波段不为盲元时，依据盲元在两个波段中的邻域相似性来作为补偿系数，取用两个波段的信息对盲元进行补偿，如(3)，(4)所示；具体补偿过程如下：
25.(1)像元(i，j)在两个波段都为盲元，且此像元未位于边缘：
[0026][0027]
(2)像元(i，j)在两个波段都为盲元，且此像元位于边缘：
[0028][0029]
(3)像元(i，j)在一个波段为盲元，而在另一个波段非盲元，且此像元未位于边缘：
[0030][0031]
此时，xo(i，j)为盲元补偿后盲元所在波段输出图像中(i，j)位置处的像素值；另一波段不需要补偿；
[0032]
(4)像元(i，j)在一个波段为盲元，而在另一个波段非盲元，且此像元位于边缘：
[0033][0034]
此时，xo(i，j)为盲元补偿后盲元所在波段输出图像中(i，j)位置处的像素值；另一波段不需要补偿；
[0035]
其中，xi(i，j)、yi(i，j)为双波段输入图像中(i，j)位置处的像素值；xo(i，j)、yo(i，j)为盲元补偿后双波段输出图像中(i，j)位置处的像素值，m、n代表(i，j)像素周边m
×
n大小的邻域，ex(i，j)、ey(i，j)为边缘图像中(i，j)位置处像素值；m、n取值范围没有限制，只要小于红外图像大小即可。
[0036]
盲元邻域相似性通过余弦相似度dist(x
imn
，y
imn
)来度量，dist(x
imn
，y
imn
)计算公式如下：
[0037][0038]
当相似度高时，可用于盲元补偿的另一波段信息占比大，当相似度低的时候，另一
波段信息占比少。
[0039]
本发明中使用的盲元补偿方法，发挥了双波段红外焦平面阵列的优势，可以同时利用两个波段的信息来对盲元进行补偿。当像元在一个波段为盲元，另一个波段非盲元时，可以同时使用盲元波段的邻域信息以及非盲元波段的真实场景信息来进行盲元补偿，有助于提升大盲元簇(如7
×
7以上盲元簇)及边缘位置处盲元的补偿效果，相较于单波段盲元补偿算法可以更好的恢复出盲元位置的真实情况。
[0040]
本发明与现有技术相比，其有益效果为：
[0041]
本发明中使用的盲元补偿方法，保留了单波段盲元补偿算法快速便捷的优势，又发挥了双波段红外焦平面阵列的优势，可以同时利用两个波段的信息来对盲元进行补偿。本发明首先对盲元在双波段中位置分布的情况进行了分类，针对性的使用了不同的补偿策略对其进行补偿。当像元在一个波段为盲元，另一个波段非盲元时，可以同时使用盲元波段的邻域信息以及非盲元波段的真实场景信息来进行盲元补偿，有助于提升大盲元簇(如7
×
7以上盲元簇)及边缘位置处盲元的补偿效果，相较于单波段盲元补偿算法可以更好的恢复出盲元位置的真实情况。
附图说明
[0042]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0043]
图1为本发明基于双波段的红外图像盲元补偿方法的流程图；
[0044]
图2为本发明基于双波段的红外图像盲元补偿方法的含盲元的中长波示例图；其中，(a)长波图像；(b)为中波图像；
[0045]
图3为补偿效果图；其中，(a)为采用现有方法；(b)为采用本发明方法。
具体实施方式
[0046]
下面结合实施例对本发明作进一步的详细描述。
[0047]
本领域技术人员将会理解，下列实施例仅用于说明本发明，而不应视为限定本发明的范围。实施例中未注明具体技术或条件者，按照本领域内的文献所描述的技术或条件或者按照产品说明书进行。所用材料或设备未注明生产厂商者，均为可以通过购买获得的常规产品。
[0048]
实施例1
[0049]
一种基于双波段的红外图像盲元补偿方法，包括如下步骤：
[0050]
步骤s1，对输入的双波段红外图像xi(i，j)、yi(i，j)分别进行盲元检测；
[0051]
步骤s2，对双波段红外图像进行边缘提取，提取边缘图像ex(i，j)、ey(i，j)；
[0052]
步骤s3，根据步骤s1检测情况及步骤s2提取情况，对双波段红外图像的像元情况进行分类，并根据分类情况对像元进行补偿后输出图像。
[0053]
实施例2
[0054]
一种基于双波段的红外图像盲元补偿方法，包括如下步骤：
[0055]
步骤s1，对输入的双波段红外图像xi(i，j)、yi(i，j)分别进行盲元检测，设立阈值h，初始盲元表置1，若某一波段像元(i，j)连续h帧被判断为盲元，则对应波段盲元表bpx(i，j)或bpy(i，j)置0；
[0056]
步骤s2，对双波段红外图像进行边缘提取，提取边缘图像ex(i，j)、ey(i，j)，对边缘部分置0，非边缘部分置1；
[0057]
步骤s3，根据步骤s1检测情况及步骤s2提取情况，对双波段红外图像的像元情况进行分类，并根据分类情况对像元进行补偿后输出图像。
[0058]
所述的双波段红外图像包括中长波段、中短波段、短长波段、中中波段图像。
[0059]
实施例3
[0060]
一种基于双波段的红外图像盲元补偿方法，包括如下步骤：
[0061]
步骤s1，对输入的双波段红外图像xi(i，j)、yi(i，j)分别进行盲元检测，设立阈值h，初始盲元表置1，若某一波段像元(i，j)连续h帧被判断为盲元，则对应波段盲元表bpx(i，j)或bpy(i，j)置0；
[0062]
步骤s2，对双波段红外图像进行边缘提取，提取边缘图像ex(i，j)、ey(i，j)，对边缘部分置0，非边缘部分置1；
[0063]
步骤s3，根据步骤s1检测情况及步骤s2提取情况，对双波段红外图像的像元情况进行分类，并根据分类情况对像元进行补偿后输出图像。
[0064]
所述的双波段红外图像包括中长波段、中短波段、短长波段、中中波段图像。
[0065]
步骤s2中，使用sobel算子进行边缘提取。
[0066]
步骤s3中，对双波段红外图像的像元情况进行分类，具体为：
[0067]
(1)像元在一个波段为盲元，而在另一个波段非盲元，且此像元未位于边缘；
[0068]
(2)像元在一个波段为盲元，而在另一个波段非盲元，且此像元位于边缘；
[0069]
(3)像元在两个波段都为盲元，且此像元未位于边缘；
[0070]
(4)像元在两个波段都为盲元，且此像元位于边缘。
[0071]
实施例4
[0072]
一种基于双波段的红外图像盲元补偿方法，包括如下步骤：
[0073]
步骤s1，对输入的双波段红外图像xi(i，j)、yi(i，j)分别进行盲元检测，设立阈值h，初始盲元表置1，若某一波段像元(i，j)连续h帧被判断为盲元，则对应波段盲元表bpx(i，j)或bpy(i，j)置0：
[0074]
步骤s2，对双波段红外图像进行边缘提取，提取边缘图像ex(i，j)、ey(i，j)，对边缘部分置0，非边缘部分置1；
[0075]
步骤s3，根据步骤s1检测情况及步骤s2提取情况，对双波段红外图像的像元情况进行分类，并根据分类情况对像元进行补偿后输出图像。
[0076]
所述的双波段红外图像包括中长波段、中短波段、短长波段、中中波段图像。
[0077]
步骤s2中，使用sobel算子进行边缘提取。
[0078]
步骤s3中，对双波段红外图像的像元情况进行分类，具体为：
[0079]
(1)像元在一个波段为盲元，而在另一个波段非盲元，且此像元未位于边缘；
[0080]
(2)像元在一个波段为盲元，而在另一个波段非盲元，且此像元位于边缘；
[0081]
(3)像元在两个波段都为盲元，且此像元未位于边缘；
[0082]
(4)像元在两个波段都为盲元，且此像元位于边缘。
[0083]
步骤s3中，根据分类情况，采取不同的策略对像元进行补偿，当像元在两个波段都为盲元时，采取单波段的补偿策略对盲元进行补偿，如(1)，(2)所示；当像元仅在单波段为盲元，而另一波段不为盲元时，依据盲元在两个波段中的邻域相似性来作为补偿系数，取用两个波段的信息对盲元进行补偿，如(3)，(4)所示；具体补偿过程如下：
[0084]
(1)像元(i，j)在两个波段都为盲元，且此像元未位于边缘：
[0085][0086]
(2)像元(i，j)在两个波段都为盲元，且此像元位于边缘：
[0087][0088]
(3)像元(i，j)在一个波段为盲元，而在另一个波段非盲元，且此像元未位于边缘：
[0089][0090]
(4)像元(i，j)在一个波段为盲元，而在另一个波段非盲元，且此像元位于边缘：
[0091][0092]
其中，xi(i，j)、yi(i，j)为双波段输入图像中(i，j)位置处的像素值；xo(i，j)、yo(i，j)为盲元补偿后双波段输出图像中(i，j)位置处的像素值，m、n代表(i，j)像素周边m
×
n大小的邻域，ex(i，j)、ey(i，j)为边缘图像中(i，j)位置处像素值；
[0093]
盲元邻域相似性通过余弦相似度dist(x
imn
，y
imn
)来度量，dist(x
imn
，y
imn
)计算公式如下：
[0094][0095]
应用实例
[0096]
结合图1，以384
×
288像素中/长双波段红外图像为例，所述方法包括：
[0097]
步骤s1，输入384
×
288像素中/长双波段红外图像，设阈值h为10，分别对其进行中波盲元检测和长波盲元检测，若某一波段像元(i，j)连续10帧被判断为盲元，则对应波段盲元表bpx(i，j)或bpy(i，j)置0；
[0098]
步骤s2，使用sobel算子，分别对中波图像和长波图像进行边缘提取，并将提取的边缘图像中的边缘部分置0，非边缘部分置1；
[0099]
步骤s3，根据步骤s1检测情况及步骤s2提取情况，对双波段红外图像的像元情况进行分类，并根据分类情况对像元进行补偿后输出图像。
[0100]
其中m＝3，n＝3。
[0101]
步骤s3中，对双波段红外图像的像元情况进行分类，具体为：
[0102]
(1)像元在一个波段为盲元，而在另一个波段非盲元，且此像元未位于边缘；
[0103]
(2)像元在一个波段为盲元，而在另一个波段非盲元，且此像元位于边缘；
[0104]
(3)像元在两个波段都为盲元，且此像元未位于边缘；
[0105]
(4)像元在两个波段都为盲元，且此像元位于边缘。
[0106]
步骤s3中，根据分类情况，采取不同的策略对像元进行补偿，当像元在两个波段都为盲元时，采取单波段的补偿策略对盲元进行补偿，如(1)，(2)所示；当像元仅在单波段为盲元，而另一波段不为盲元时，依据盲元在两个波段中的邻域相似性来作为补偿系数，取用两个波段的信息对盲元进行补偿，如(3)，(4)所示；具体补偿过程如下：
[0107]
(1)像元(i，j)在两个波段都为盲元，且此像元未位于边缘：
[0108][0109]
(2)像元(i，j)在两个波段都为盲元，且此像元位于边缘：
[0110][0111]
(3)像元(i，j)在一个波段为盲元，而在另一个波段非盲元，且此像元未位于边缘：
[0112][0113]
(4)像元(i，j)在一个波段为盲元，而在另一个波段非盲元，且此像元位于边缘：
[0114][0115]
其中，xi(i，j)、yi(i，j)为双波段输入图像中(i，j)位置处的像素值；xo(i，j)、yo(i，j)为盲元补偿后双波段输出图像中(i，j)位置处的像素值，m、n代表(i，j)像素周边m
×
n大小的邻域，ex(i，j)、ey(i，j)为边缘图像中(i，j)位置处像素值；
[0116]
盲元邻域相似性通过余弦相似度dist(x
imn
，y
imn
)来度量，dist(x
imn
，y
imn
)计算公式如下：
[0117][0118]
实例图像如图2所示，其中，中波图像无盲元，长波图像中存在零散的小盲元及一个7
×
10的大盲元簇。
[0119]
效果图如图3所示。本实例的补偿效果如图3(b)所示，相较于通常的单波段补偿算法结果如图3(a)所示，本发明对于零散小盲元的补偿效果和单波段效果相似，但对于7
×
10的大盲元簇的补偿效果明显优于单波段补偿算法。
[0120]
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。