发明名称 基于多层次迭代可视化优化的全聚焦合成孔径透视成像方法
摘要 发明公开了一种基于多层次迭代可视化优化的全聚焦合成孔径透视成像方法,用于解决现有基于分形的红外弱小目标检测方法检测效果差的技术问题。技术方案是首先对于一幅M×N大小的红外图像I,首先进行中值滤波,然后计算局部熵图HI。获取目标的粗略位置(r,c)之后,以该点为中心在原图I上下左右各扩散e个单位,提取目标在原始红外图像I的小区域Q。对提取的R×C大小的目标区域Q利用分形组合特征方法计算分维像D,精确提取目标。本发明将分形维数、分形几何空间变化率和分形几何误差有机的集合来表示图像的分形特征,得到的分形组合特征能更好的抑制自然背景对人造目标检测的影响,同时也不会显著的增加时间复杂度,在红外小目标检测方面能取得更好的效果。
申请公布号 CN103632370A 申请公布日期 2014.03.12
申请号 CN201310646293.6 申请日期 2013.12.02
申请人 西北工业大学 发明人 张艳宁;郗润平;张福俊
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于多层次迭代可视化优化的全聚焦合成孔径透视成像方法,其特征在于包括以下步骤:步骤一、对于一幅M×N大小的红外图像I,首先进行中值滤波,然后计算局部熵图H<sub>I</sub>;取m×m的滑动窗口W<sub>1</sub>,按从左到右,从上到下的方式遍历图像I,计算图像的局部熵图H<sub>I</sub>;记点(i,j)是窗口W<sub>1</sub>中心点坐标,则该点的熵值计算如下<![CDATA[<math><mrow><msub><mi>H</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>a</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>b</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>p</mi><mi>ab</mi></msub><mi>log</mi><msub><mi>p</mi><mi>ab</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA0000428451150000012.GIF" wi="637" he="156" />W<sub>1</sub>(a,b)、W<sub>1</sub>(q,w)分别是滑动窗口对应点的灰度值;得到局部熵图H<sub>I</sub>之后,计算得到最小熵值对应的坐标(r,c),获取目标的粗略位置;步骤二、获取目标的粗略位置(r,c)之后,以该点为中心在原图I上下左右各扩散e个单位,提取目标在原始红外图像I的小区域Q;记区域Q在原图I中左上角的坐标为(r<sub>1</sub>,c<sub>1</sub>),左下角的坐标为(r<sub>2</sub>,c<sub>2</sub>),则扩散后的坐标计算如下:<img file="FDA0000428451150000013.GIF" wi="1594" he="164" /><img file="FDA0000428451150000014.GIF" wi="1598" he="158" />式中,M,N分别是I的高度、宽度,由以上步骤提取得到R×C大小的目标区域Q;步骤三、对提取的R×C大小的目标区域Q利用分形组合特征方法计算分维像D,精确提取目标;取一个n×n的滑动窗口W<sub>2</sub>,按从左到右,从上到下的方式遍历区域Q;把W<sub>2</sub>划分成s×s的小块,s是介于1到n/2的之间的整数,得到t=n/s;记点(i,j)是窗口W<sub>2</sub>中心点坐标,求得窗口在不同尺度s下的盒子数B<sub>s</sub>(i,j);当s=4时,t=2,此时W<sub>2</sub>被划分四个4×4的小块,记为<img file="FDA0000428451150000016.GIF" wi="336" he="67" />分别获取四个小块的最大灰度值l<sup>λ</sup>和最小灰度值k<sup>λ</sup>,则<![CDATA[<math><mrow><msub><mi>B</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>&lambda;</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></munderover><mrow><mo>(</mo><msup><mi>l</mi><mi>&lambda;</mi></msup><mo>/</mo><mi>s</mi><mo>-</mo><msup><mi>k</mi><mi>&lambda;</mi></msup><mo>/</mo><mi>s</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>由此计算在尺度s=4的盒子数;对s=2时重复上面的操作,求取B<sub>2</sub>(i,j);然后利用最小二乘法线性拟合点对(log(t),log(B<sub>s</sub>(i,j))),其中t=4,2;s=2,4.即满足公式(6)log(B<sub>s</sub>)=Flog(t)+c    (6)线性拟合得到的F表示点(i,j)的分形维数,c的绝对值表示分形几何空间变化率,另有分形拟合误差表示如下:<![CDATA[<math><mrow><mi>&epsiv;</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>g</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><mo>[</mo><mi>log</mi><mrow><mo>(</mo><msubsup><mi>B</mi><mi>s</mi><mi>g</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mi>F</mi><mi>log</mi><mrow><mo>(</mo><msup><mi>t</mi><mi>g</mi></msup><mo>)</mo></mrow><mo>-</mo><mo>|</mo><mi>c</mi><mo>|</mo><msup><mo>]</mo><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(7)表示在不同尺度下的拟合误差和;由此计算分形组合特征维数d<img file="FDA0000428451150000022.GIF" wi="1223" he="166" />则窗口中心点(i,j)的分形组合特征维数为d,得到D(i,j)=d.当滑动窗口W<sub>2</sub>扫描整个区域图像Q后,得到分维图像D;最后对D进行二值分割,形态学滤波,得到目标的分割图,完成目标提取。
地址 710072 陕西省西安市友谊西路127号