海上阳光亮带干扰下的红外图像补偿方法技术

一种海上阳光亮带干扰下的红外图像补偿方法。其步骤为：(1)读入一幅阳光亮带干扰下的红外图像；(2)利用中值滤波，对输入的阳光亮带干扰下的红外图像进行去噪；(3)通过最大类间方差Otsu算法，划分受阳光亮带干扰区域与未受亮带干扰区域；(4)计算补偿模板参数；(5)补偿图像中受阳光亮带干扰区域；(6)对阳光亮带干扰下的红外图像进行增益变换处理。本发明专利技术对单幅图像进行处理，对阳光亮带干扰区域面积无特殊要求，具有用范围广、计算复杂度低，实时性高的特点。将受阳光亮带干扰的区域进行填充替换，使得本发明专利技术具有了图像信息利用率高的优点。本发明专利技术可用于图像处理、目标识别等技术领域。

【技术实现步骤摘要】
海上阳光亮带干扰下的红外图像补偿方法
本专利技术属于图像处理
，更进一步涉及红外图像处理
中的海上阳光亮带干扰下的红外图像补偿方法。本专利技术可对红外图像重新量化和增强，提高海上阳光亮带干扰下的目标识别能力。
技术介绍
复杂海面背景情况下，由于成像器件、大气、拍摄角度、时间、气象等诸多因素影响，以及不同的海面波浪状况对光照的反射能力不同，使得红外图像的亮度、对比度等信息存在很大变化，对目标识别和跟踪造成很大影响。当水面波动时，太阳光照射在海浪的顶端出现闪光点，形成太阳耀光。尤其是在清晨和傍晚，太阳天顶角小，对海面照射范围大，成像系统受干扰影响观察视角范围广，大片的太阳耀光形成阳光亮带，在这两个时间段观测识别海上目标极易受到大面积阳光亮带的影响，红外探测系统CCD图像传感器临近饱和状态，此时获取的图像整体亮度较高，对比度极低且受探测系统的强噪声干扰，目标与背景难以区分，使得传统的目标侦查和识别难以有效工作。受阳光亮带干扰严重的红外图像，要实现目标的有效识别，关键问题在于减弱海面阳光反射对红外目标成像的影响。目前解决该类问题，通常有两种方法：一种是将海上的阳光反射当作杂波噪声进行滤波处理，这种方法实时性强，但当图像中存在大片亮带时该算法失效。另一种是通过采集序列图像，比对干扰及非干扰条件的图像特点，获取关注区域的目标特征，该方法需要预处理大量的序列图像，影响侦查识别的实时性。华中科技大学在其申请的专利“一种红外图像海面目标检测识别定位方法”(专利申请号200910063265.5，公开号CN101609504A)中提出了一种红外图像海面目标检测识别定位方法。该方法通过基于同质变换的方式将背景和噪声预处理，消除复杂背景下其他热源的干扰；为了克服目标出现块状效应，利用单峰分割方法提取目标疑似区域。再在疑似区域上提取出目标。该方法存在的不足是，在海面有阳光亮带时，探测器视场变亮，目标不是最亮点，且与背景灰度差别很小。利用单峰分割时，灰度直方图中背景与目标灰度相似，单拐点不明显，导致分割阈值不准确，分离背景和目标较为困难。北京航空航天大学在其申请的专利“一种基于偏振遥感探测的波浪水面太阳耀光剥离方法”(专利申请号201210509290.3，公开号CN102998667A)中提出了一种基于偏振遥感探测的波浪水面太阳耀光剥离方法。该方法通过特定条件采集序列图像，对获取的不同相位延迟量的偏振强度图像通过对定标仪器矩阵求逆，解算出四个Stokes参量图像、偏振度图像与偏振角图像，基于时间序列分析、频率域降噪与图像融合等相关数据处理技术，对四幅波浪水面太阳耀光区偏振度图像进行处理，从而达到波浪水面太阳耀光剥离的目的。该方法存在的不足是，对特定条件下获取的待处理图像要求高，要求快速切换偏振片方位角，以至于研究对象不发生太大空间位置移动，且每张偏振度图像的饱和区域不能覆盖目标物区域一半以上的面积比例，导致图像数据利用率低，并且需要输入序列图像进行处理，工作量大，影响目标识别效率。
技术实现思路
本专利技术的目的在于克服上述现有技术的不足，提供一种海上阳光亮带干扰下的红外图像补偿方法。本方法可以明显减弱海上阳光亮带干扰对红外目标成像的影响，增强图像对比度，提高目标识别效率。实现本方法的基本思路是：利用中值滤波对输入的阳光亮带干扰下的红外图像进行去噪，通过最大类间方差Otsu算法，划分受阳光亮带干扰区域与未受亮带干扰区域，计算补偿模板参数，利用补偿模板对受阳光亮带干扰区域的像素进行补偿，最后对图像进行增益变换，从而达到对海上阳光亮带干扰下的红外图像的补偿。为实现上述目的，本方法的具体步骤如下：(1)读入一幅阳光亮带干扰下的红外图像；(2)利用中值滤波，对输入的阳光亮带干扰下的红外图像进行去噪；(3)利用最大类间方差Otsu算法，得到经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像的阈值，将经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像的像素灰度值大于或等于图像阈值的点作为受阳光亮带干扰区域，将像素灰度值小于图像阈值的点作为未受阳光亮带干扰区域；(4)计算补偿模板参数：(4a)统计经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像未受阳光亮带干扰区域中的像素个数；(4b)按照下式，计算未受阳光亮带干扰区域中所有像素灰度的平均值：其中，表示未受阳光亮带干扰区域中所有像素灰度的平均值，S表示未受阳光亮带干扰区域，xij表示未受阳光亮带干扰区域S中横坐标为i、纵坐标为j处像素的灰度值，n表示未受阳光亮带干扰区域S中的像素个数；(4c)按照下式，计算未受阳光亮带干扰区域中所有像素灰度的方差：其中，σ2表示未受阳光亮带干扰区域中所有像素灰度的方差，n表示未受阳光亮带干扰区域S中的像素个数，xij表示未受阳光亮带干扰区域S中横坐标为i、纵坐标为j处像素的灰度值，表示未受阳光亮带干扰区域所有像素灰度的平均值；(5)补偿图像中受阳光亮带干扰区域：(5a)创建一个M行N列的矩阵模板，行数、列数分别与受阳光亮带干扰的红外图像纵向像素个数、横向像素个数相同，矩阵元素服从高斯分布；(5b)用模板矩阵中的矩阵元素值替代受阳光亮带干扰区域中像素的灰度值；(6)按照下式，将经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像进行增益变换处理：其中，yij表示经过增益变换后阳光亮带干扰下的红外图像中横坐标为i、纵坐标为j处像素的灰度值，G表示未受阳光亮带干扰时红外探测器的系统增益，H表示受阳光亮带干扰时红外探测器的系统增益，xij表示阳光亮带干扰下的红外图像中横坐标为i，纵坐标为j处像素的灰度值。本专利技术与现有技术相比具有以下优点：第一，由于本专利技术在对海上阳光亮带干扰下的红外图像进行补偿时，使用单幅图像进行处理，且对阳光亮带干扰区域面积无特殊要求，克服了现有技术中需要处理序列图像和对图像特征要求高的不足，使得本专利技术具有了使用范围广、计算复杂度低，实时性高的特点。第二，由于本专利技术在对海上阳光亮带干扰下的红外图像进行补偿时，将受阳光亮带干扰的区域进行填充替换，克服了现有技术中将海上阳光反射作为噪声去除丢失图像信息的不足，使得本专利技术具有了图像信息利用率高的优点。附图说明图1为本专利技术的流程图；图2为本专利技术实施例1的输入图像和效果图像；图3为本专利技术实施例2的输入图像和效果图像。具体实施方式下面结合附图对本专利技术做进一步的详细描述。参照附图1，本专利技术的具体步骤如下：步骤1，读入图像。本专利技术读入的一幅阳光亮带干扰下的红外图像是选用的中波热像仪在2013年7月20日在山东青岛拍摄的红外图像。本专利技术的实施例1和实施例2中，输入的阳光亮带干扰下的红外图像图2(a)、图3(a)，输出的图像图2(b)、图3(b)的图像分辨率为240×320像素。步骤2，图像去噪。利用中值滤波，对输入的阳光亮带干扰下的红外图像进行去噪。中值滤波方法的步骤如下：第一步，将中值滤波器的滑动窗口设置为3×3个像素的正方形；第二步，用正方形窗口沿着阳光亮带干扰下的红外图像行的方向逐像素滑动，在每一次滑动期间，将正方形窗口内的所有像素按照像素值从小到大进行排序，选取排序结果的中间值，替代正方形窗口中心位置的像素灰度值；第三步，判断是否处理完阳光亮带干扰下的红外图像中的所有像素，若是，中值滤波完成，得到去本文档来自技高网...

【技术保护点】
一种海上阳光亮带干扰下的红外图像补偿方法，包括以下步骤：(1)读入一幅阳光亮带干扰下的红外图像；(2)利用中值滤波，对输入的阳光亮带干扰下的红外图像进行去噪；(3)利用最大类间方差Otsu算法，得到阳光亮带干扰下的红外图像的阈值，将阳光亮带干扰下的红外图像的像素灰度值大于或等于图像阈值的点，作为受阳光亮带干扰区域，将像素灰度值小于图像阈值的点作为未受阳光亮带干扰区域；(4)计算补偿模板参数：(4a)统计阳光亮带干扰下的红外图像未受阳光亮带干扰区域中的像素个数；(4b)按照下式，计算未受阳光亮带干扰区域中所有像素灰度的平均值：x‾=Σ(i,j)∈Snxijn]]>其中，表示未受阳光亮带干扰区域中所有像素灰度的平均值，S表示未受阳光亮带干扰区域，xij表示未受阳光亮带干扰区域S中横坐标为i、纵坐标为j处像素的灰度值，n表示未受阳光亮带干扰区域S中的像素个数；(4c)按照下式，计算未受阳光亮带干扰区域中所有像素灰度的方差：σ2=1nΣ(i,j)∈Sn(xij-x‾)2]]>其中，σ2表示未受阳光亮带干扰区域中所有像素灰度的方差，n表示未受阳光亮带干扰区域S中的像素个数，xij表示未受阳光亮带干扰区域S中横坐标为i、纵坐标为j处像素的灰度值，表示未受阳光亮带干扰区域所有像素灰度的平均值；(5)补偿图像中受阳光亮带干扰区域：(5a)创建一个M行N列的矩阵模板，行数、列数分别与受阳光亮带干扰的红外图像纵向像素个数、横向像素个数相同，矩阵元素服从高斯分布；(5b)用模板矩阵中的矩阵元素值替代受阳光亮带干扰区域中像素的灰度值；(6)按照下式，将阳光亮带干扰下的红外图像进行增益变换处理：yij=GHxij]]>其中，yij表示经过增益变换后阳光亮带干扰下的红外图像中横坐标为i、纵坐标为j处像素的灰度值，G表示未受阳光亮带干扰时红外探测器的系统增益，H表示受阳光亮带干扰时红外探测器的系统增益，xij表示阳光亮带干扰下的红外图像中横坐标为i，纵坐标为j处像素的灰度值。...

【技术特征摘要】
1.一种海上阳光亮带干扰下的红外图像补偿方法，包括以下步骤：(1)读入一幅阳光亮带干扰下的红外图像；(2)利用中值滤波，对输入的阳光亮带干扰下的红外图像进行去噪；(3)利用最大类间方差Otsu算法，得到经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像的阈值，将经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像的像素灰度值大于或等于图像阈值的点，作为受阳光亮带干扰区域，将像素灰度值小于图像阈值的点作为未受阳光亮带干扰区域；(4)计算补偿模板参数：(4a)统计经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像未受阳光亮带干扰区域中的像素个数；(4b)按照下式，计算未受阳光亮带干扰区域中所有像素灰度的平均值：其中，表示未受阳光亮带干扰区域中所有像素灰度的平均值，S表示未受阳光亮带干扰区域，xij表示未受阳光亮带干扰区域S中横坐标为i、纵坐标为j处像素的灰度值，n表示未受阳光亮带干扰区域S中的像素个数；(4c)按照下式，计算未受阳光亮带干扰区域中所有像素灰度的方差：其中，σ2表示未受阳光亮带干扰区域中所有像素灰度的方差，n表示未受阳光亮带干扰区域S中的像素个数，xij表示未受阳光亮带干扰区域S中横坐标为i、纵坐标为j处像素的灰度值，表示未受阳光亮带干扰区域所有像素灰度的平均值；(5)补偿图像中受阳光亮带干扰区域：(5a)创建一个M行N列的矩阵模板，行数、列数分别与受阳光亮带干扰的红外图像纵向像素个数、横向像素个数相同，矩阵元素服从高斯分布；(5b)用模板矩阵中的矩阵元素值替代受阳光亮带干扰区域中像素的灰度值；(6)按照下式，将经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像进行增益变换处理：其中，yij表示经过增益变换后阳光亮带干扰下的红外图像中横坐标为i、纵坐标为j处像素的灰度值，G表示未受阳光亮带干扰时红外探测器的系统增益，H表示受阳光亮带干扰时红外探测器的系统增益，xij表示阳光亮带干扰下的红外图像中横坐标为i，纵坐标为j处像素的灰度值。2.根据权利要求1所述的海上阳光亮带干扰下的红外图像补偿方法，其特征在于...

【专利技术属性】
技术研发人员：黄曦，柴国贝，陆瑶，张建奇，王晓蕊，刘德连，董维科，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61