发明名称 一种基于模糊距离的红外小目标检测方法
摘要 本发明为有效检测复杂背景下的红外小目标,公开了一种基于模糊距离的红外小目标检测方法,涉及数字图像处理技术领域。本发明首先针对小目标的出现会引起局部纹理发生较大变化这一特征,提出一种模糊距离概念,从而将局部纹理的变化转化为模糊距离的度量;其次针对小目标的尺寸会随成像距离的改变而发生相应变化这一特点,提出一种多尺度模糊距离及多尺度模糊距离图,能剔除大量背景杂波和噪声干扰;然后通过迭代运算,有效抑制残留背景和噪声,增强目标;最后利用自适应阈值检测目标,该检测方法简单且有效。
申请公布号 CN105869156A 申请公布日期 2016.08.17
申请号 CN201610177589.1 申请日期 2016.03.25
申请人 中国科学院武汉物理与数学研究所 发明人 周欣;邓鹤;孙献平;刘买利;叶朝辉
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 武汉宇晨专利事务所 42001 代理人 李鹏;王敏锋
主权项 一种基于模糊距离的红外小目标检测方法,其特征在于,包括以下步骤:步骤1、初始化相关参数:设置最大迭代次数L,其中L为正整数;初始化迭代次数索引k=0;设置最大局部窗口大小m×n,其中m和n均为大于1的正奇数;步骤2、求解红外图像I每个像素点的模糊距离,包括以下步骤:步骤2.1、获得单帧红外图像I每一个像素点(x,y)的邻域空间集{Ω<sub>l</sub>|l=1,2,...,s},其中s=min{0.5·(m‑1),0.5·(n‑1)},Ω<sub>l</sub>的大小为(2l+1)×(2l+1),像素点(x,y)的邻域空间Ω<sub>l</sub>的定义为Ω<sub>l</sub>={(p,q)|max(|p‑x|,|q‑y|)≤l},(p,q)是邻域空间Ω<sub>l</sub>内的像素点;步骤2.2、计算每一个像素点(x,y)的各个邻域空间Ω<sub>l</sub>内像素的灰度均值D<sub>l</sub>(x,y):<maths num="0001"><math><![CDATA[<mrow><msub><mi>D</mi><mi>l</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mo>#</mo><msub><mi>&Omega;</mi><mi>l</mi></msub></mrow></mfrac><munder><mo>&Sigma;</mo><mrow><mo>(</mo><mi>a</mi><mo>,</mo><mi>b</mi><mo>)</mo><mo>&Element;</mo><msub><mi>&Omega;</mi><mi>l</mi></msub></mrow></munder><mi>I</mi><mrow><mo>(</mo><mi>a</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mo>,</mo><mi>l</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>s</mi></mrow>]]></math><img file="FDA0000950102000000011.GIF" wi="670" he="111" /></maths>其中,#Ω<sub>l</sub>表示邻域空间Ω<sub>l</sub>内像素点的数目,I(a,b)表示邻域空间Ω<sub>l</sub>内像素点(a,b)处的灰度值。步骤2.3、计算每一个像素点(x,y)所对应的最大邻域空间Ω<sub>s</sub>与其它各个邻域空间Ω<sub>i</sub>,i=1,2,...,s‑1之间的模糊距离E<sub>i</sub>:<maths num="0002"><math><![CDATA[<mrow><msub><mi>E</mi><mi>i</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><mrow><mi>e</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mrow><mo>(</mo><msup><mi>e</mi><mrow><mn>1</mn><mo>-</mo><mrow><mo>|</mo><mrow><msub><mi>D</mi><mi>i</mi></msub><mo>+</mo><msub><mi>D</mi><mi>s</mi></msub></mrow><mo>|</mo></mrow><mo>&CenterDot;</mo><mrow><mo>|</mo><mrow><msub><mi>D</mi><mi>i</mi></msub><mo>-</mo><msub><mi>D</mi><mi>s</mi></msub></mrow><mo>|</mo></mrow></mrow></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000950102000000012.GIF" wi="558" he="102" /></maths>其中e为自然常数,D<sub>s</sub>表示最大邻域空间Ω<sub>s</sub>内像素的灰度均值,D<sub>i</sub>表示第i个邻域空间Ω<sub>i</sub>内像素的灰度均值;步骤3、求解多尺度模糊距离图:遍历红外图像I中每一个像素点,得到每一个像素点的多尺度模糊距离E(x,y),然后根据每一个像素点的多尺度模糊距离E(x,y)并通过归一化方法获得红外图像I的多尺度模糊距离图E;步骤4、迭代停止准则判断:迭代次数索引k加1,判断迭代次数索引k与最大迭代次数L之间的关系,若k&lt;L,把步骤3所获得的多尺度模糊距离图E作为新的红外图像I,返回步骤2;若k≥L,停止迭代,把步骤3所获得的多尺度模糊距离图E作为最终的滤波结果,进行步骤5;步骤5、求解自适应阈值T:对经过步骤4所获得的最终滤波结果,即多尺度模糊距离图E,求解自适应阈值T,并通过自适应阈值T对多尺度模糊距离图E进行二值化,检测出红外小目标。
地址 430071 湖北省武汉市武昌区小洪山西30号