[发明专利]基于局部中值直方图的自适应红外图像去条纹算法

说明书

本发明公开了一种基于局部中值直方图的自适应红外图像去条纹算法，通过计算滑动窗口内场景复杂度，找出原始红外图像中每一列与相邻列场景变化最小的窗口。计算该窗口所有列的累积直方图，得到窗口中心列的中值直方图，根据中值直方图计算出窗口中心列的各个像素校正后的灰度值，用原始灰度值均值减去校正后的灰度值均值就可以得到该列的条纹值。将整列都减去条纹值即为该列最终输出的灰度值。本发明通过滑动窗口找到图像中场景变化不明显的区域，校正效果不受场景影响，可以消除传统中值直方图均衡去条纹算法中图像场景变化剧烈时引入新的噪声的问题。

技术领域

本发明属于红外图像非均匀性校正领域，具体涉及一种基于局部中值直方图的自适应红外图像去条纹算法。

背景技术

目前，红外图像已被广泛应用于工业、医学和军事等领域来进行低可视度下的侦测。在理想情况下，红外成像系统对均匀辐射的红外光，获取的数字图像上每个像素点的灰度值应该完全一样。但实际上，受限于固态电子的制造工艺，探测器上的光敏元件(像元)往往伴随着杂质浓度不均，厚度不等，有效光敏面积做不到绝对平均等问题，像元之间的光电转换效率各不相同，对均匀辐射的景物的成像不均匀。另外，图像数据读出电路各通道之间的差，会导致图像出现呈列分布的固定条纹噪声。这就要求我们图像进行非均匀性校正，使图像得到更好的视觉效果。

常用的红外图像非均匀校正技术主要有定标法和场景法两种。定标法主要有两点校正法、多点校正法等；场景法主要有时域高通滤波法、恒定统计法、神经网络法等。但这两种方法在应用时都有很大的局限性。

近年来，国内外学者开始关注在在静态场景中，或者说在单帧图像中去除红外条纹非均匀性的方法。

Tendero和Gilles研究在单帧图像内去除图像上的非均匀性，提出了一种利用中值直方图均衡化的去条纹非均匀性校正算法。中值直方图最初是用来修正各相机中传感器增益之间的差异。若两幅图像的累计直方图分别是H₁和H₂，则其中间直方图的计算公式如下：

由于红外图像固定条纹噪声一般来说并不是孤立的，所以可以将上述方法用在红外图像上，根据相邻列的累计直方图计算某一列的中间直方图，用中间直方图来替代该列的累计直方图。通过这种方法使单幅图像上每一列的信息作用到其他列上，以进行单幅图像的非均匀性校正。该方法使用前提时图像列与列之间直方图相差不大。在图像中有复杂场景时，可能会出现校正效果不佳甚至出现带状效应，发生畸变等不理想的情况。

发明内容

本发明的目的在于提供一种基于中值直方图的单帧红外图像非均匀性校正方法，用于去除红外图像中的条纹噪声。

实现本发明目的的技术解决方案为：一种基于局部中值直方图的自适应红外图像去条纹算法，实现步骤如下：

步骤1、采集一幅像素数为M×N的原始红外图像o(i,j)，其中i∈{1,…,M}，j∈{1,…,N}， i表示图像上像素的列所在位置，j表示图像上像素的行所在位置，M表示采集图像的列数， N表示采集图像的行数；

步骤2、对于上述原始红外图像o(i,j)的第x列，以第x列为中心列构造一个大小为A×B 的滑动窗口，其中，A为窗口的列数，且为奇数，B为窗口的行数，计算窗口内场景复杂度μ：

其中，k表示像素灰度值，L表示最大灰度值，表示窗口中所有像素灰度值均值，p(k) 表示窗口中灰度值为k的像素的个数，μ越大，表示场景越复杂；

步骤3、逐像素地上下移动窗口，计算以第x列为中心列的所有窗口的场景复杂度，并进行比较，找出以第x列为中心列的场景复杂度最小的窗口；

步骤4、在以第x列为中心列的场景复杂度最小的窗口中，对第x列进行中值直方图均衡，得到校正后的像素灰度值d(x,j)；