发明名称 海上阳光亮带干扰下的红外图像补偿方法
摘要 一种海上阳光亮带干扰下的红外图像补偿方法。其步骤为:(1)读入一幅阳光亮带干扰下的红外图像;(2)利用中值滤波,对输入的阳光亮带干扰下的红外图像进行去噪;(3)通过最大类间方差Otsu算法,划分受阳光亮带干扰区域与未受亮带干扰区域;(4)计算补偿模板参数;(5)补偿图像中受阳光亮带干扰区域;(6)对阳光亮带干扰下的红外图像进行增益变换处理。本发明对单幅图像进行处理,对阳光亮带干扰区域面积无特殊要求,具有用范围广、计算复杂度低,实时性高的特点。将受阳光亮带干扰的区域进行填充替换,使得本发明具有了图像信息利用率高的优点。本发明可用于图像处理、目标识别等技术领域。
申请公布号 CN104268831B 申请公布日期 2017.03.29
申请号 CN201410465310.0 申请日期 2014.09.12
申请人 西安电子科技大学 发明人 黄曦;柴国贝;陆瑶;张建奇;王晓蕊;刘德连;董维科
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 田文英;王品华
主权项 一种海上阳光亮带干扰下的红外图像补偿方法,包括以下步骤:(1)读入一幅阳光亮带干扰下的红外图像;(2)利用中值滤波,对输入的阳光亮带干扰下的红外图像进行去噪;(3)利用最大类间方差Otsu算法,得到经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像的阈值,将经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像的像素灰度值大于或等于图像阈值的点,作为受阳光亮带干扰区域,将像素灰度值小于图像阈值的点作为未受阳光亮带干扰区域;(4)计算补偿模板参数:(4a)统计经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像未受阳光亮带干扰区域中的像素个数;(4b)按照下式,计算未受阳光亮带干扰区域中所有像素灰度的平均值:<maths num="0001"><math><![CDATA[<mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><mi>S</mi></mrow><mi>n</mi></munderover><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow><mi>n</mi></mfrac></mrow>]]></math><img file="FDA0001113014170000011.GIF" wi="254" he="207" /></maths>其中,<img file="FDA0001113014170000012.GIF" wi="31" he="55" />表示未受阳光亮带干扰区域中所有像素灰度的平均值,S表示未受阳光亮带干扰区域,x<sub>ij</sub>表示未受阳光亮带干扰区域S中横坐标为i、纵坐标为j处像素的灰度值,n表示未受阳光亮带干扰区域S中的像素个数;(4c)按照下式,计算未受阳光亮带干扰区域中所有像素灰度的方差:<maths num="0002"><math><![CDATA[<mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo><mo>&Element;</mo><mi>S</mi></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0001113014170000013.GIF" wi="463" he="149" /></maths>其中,σ<sup>2</sup>表示未受阳光亮带干扰区域中所有像素灰度的方差,n表示未受阳光亮带干扰区域S中的像素个数,x<sub>ij</sub>表示未受阳光亮带干扰区域S中横坐标为i、纵坐标为j处像素的灰度值,<img file="FDA0001113014170000014.GIF" wi="36" he="55" />表示未受阳光亮带干扰区域所有像素灰度的平均值;(5)补偿图像中受阳光亮带干扰区域:(5a)创建一个M行N列的矩阵模板,行数、列数分别与受阳光亮带干扰的红外图像纵向像素个数、横向像素个数相同,矩阵元素服从高斯分布;(5b)用模板矩阵中的矩阵元素值替代受阳光亮带干扰区域中像素的灰度值;(6)按照下式,将经过中值滤波处理后的去噪的阳光亮带干扰下的红外图像进行增益变换处理:<maths num="0003"><math><![CDATA[<mrow><msub><mi>y</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mi>G</mi><mi>H</mi></mfrac><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow>]]></math><img file="FDA0001113014170000021.GIF" wi="206" he="119" /></maths>其中,y<sub>ij</sub>表示经过增益变换后阳光亮带干扰下的红外图像中横坐标为i、纵坐标为j处像素的灰度值,G表示未受阳光亮带干扰时红外探测器的系统增益,H表示受阳光亮带干扰时红外探测器的系统增益,x<sub>ij</sub>表示阳光亮带干扰下的红外图像中横坐标为i,纵坐标为j处像素的灰度值。
地址 710071 陕西省西安市太白南路2号
您可能感兴趣的专利