发明名称 基于高分辨率合成孔径雷达图像的舰船检测方法
摘要 本发明提出了一种基于高分辨率合成孔径雷达SAR图像的舰船检测方法,主要解决现有检测与鉴别技术应用于高分辨率图像时,整体检测性能差的问题。其实现步骤是:将原始图像分解为一组子孔径图像;利用多幅子孔径图像之间的相位相关性信息和各个子孔径图像的幅度信息获取二值目标检测结果;对该结果进行后处理以提取候选目标切片;对候选目标切片进行目标与背景的分割后,提取三个现有鉴别特征和一个基于协方差矩阵的新特征,生成特征向量;使用适用于小样本问题的K近邻鉴别器完成目标鉴别。本发明具有目标检测精度与鉴别性能高的优点,可用于高分辨率SAR图像中海面舰船目标的检测。
申请公布号 CN102208017A 申请公布日期 2011.10.05
申请号 CN201110140973.1 申请日期 2011.05.30
申请人 西安电子科技大学 发明人 王英华;刘宏伟
分类号 G06K9/00(2006.01)I;G06K9/62(2006.01)I;G01S13/90(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于高分辨率合成孔径雷达图像的舰船检测方法,包括:(1)预先筛选步骤:1a)将原始的单视复数型SAR图像分解为N幅子孔径图像,2≤N≤16;1b)对原始图像中的每一个像素点s,计算N幅子孔径图像中对应像素点处的N个随机变量s<sub>i</sub>,i=1,2,...,N之间的归一化相关矩阵R<sub>s</sub>,并计算检测统计量图像中对应像素点处的值ρ<sub>IC</sub>:ρ<sub>IC</sub>=log[&lt;|s<sub>1</sub>|<sup>2</sup>&gt;&lt;|s<sub>2</sub>|<sup>2</sup>&gt;...&lt;|s<sub>N</sub>|<sup>2</sup>&gt;(1-|R<sub>s</sub>|)]其中,|s<sub>i</sub>|表示s<sub>i</sub>的模值,|R<sub>s</sub>|表示R<sub>s</sub>的行列式,&lt;g&gt;表示空间平均,log(g)表示取自然对数;将检测统计量图像中每一个像素点的值ρ<sub>IC</sub>与根据检测虚警率设定的全局门限T<sub>ρ</sub>相比,如果ρ<sub>IC</sub>>T<sub>ρ</sub>,则二值检测结果图像I<sub>bw</sub>中对应像素点赋值为1且被视为目标像素,否则赋值为0且被视为背景杂波;1c)对二值检测结果图像I<sub>bw</sub>进行后处理,以提取出N<sub>c</sub>幅候选目标切片<img file="FDA0000064658510000011.GIF" wi="141" he="57" />其中G<sub>i</sub>表示第i个候选目标切片,它是从原始图像中切割出来的一幅L<sub>c</sub>×L<sub>c</sub>大小的单视复数图;(2)目标鉴别步骤:2a)对每一幅候选切片图像C<sub>i</sub>,进行目标与背景的分割,得到一幅分割标记图像,其中目标像素点都具有标记值1,背景像素点都具有标记值0,将标记值为1的像素点集合记为R<sub>tg</sub>,像素点数目记为N<sub>tg</sub>;2b)根据像素点集合R<sub>tg</sub>中的像素点数目和像素点的值,提取以下四个鉴别特征:2b1)提取标准差特征:<maths num="0001"><![CDATA[<math><mrow><msub><mi>f</mi><mn>1</mn></msub><mo>=</mo><msqrt><mfrac><mrow><msub><mi>S</mi><mn>2</mn></msub><mo>-</mo><msubsup><mi>S</mi><mn>1</mn><mn>2</mn></msubsup><mo>/</mo><msub><mi>N</mi><mi>tg</mi></msub></mrow><mrow><msub><mi>N</mi><mi>tg</mi></msub><mo>-</mo><mn>1</mn></mrow></mfrac></msqrt></mrow></math>]]></maths>其中,<maths num="0002"><![CDATA[<math><mrow><msub><mi>S</mi><mn>1</mn></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>s</mi><mo>&Element;</mo><msub><mi>R</mi><mi>tg</mi></msub></mrow></munder><msub><mrow><mn>10</mn><mi>log</mi></mrow><mn>10</mn></msub><msub><mi>I</mi><mi>s</mi></msub></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>S</mi><mn>2</mn></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>s</mi><mo>&Element;</mo><msub><mi>R</mi><mi>tg</mi></msub></mrow></munder><msup><mrow><mo>(</mo><msub><mrow><mn>10</mn><mi>log</mi></mrow><mn>10</mn></msub><msub><mi>I</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>I<sub>s</sub>表示像素点集合R<sub>tg</sub>中的像素点s的强度值;2b2)提取分形维特征:<img file="FDA0000064658510000023.GIF" wi="441" he="119" />其中M<sub>1</sub>是最亮点二值图像中标记值为1的像素数目,M<sub>2</sub>是使用2×2像素大小的盒子覆盖最亮点二值图像中标记值为1的像素点所需要的最少盒子数目;该最亮点二值图像是通过将切片图像C<sub>i</sub>中强度值最大的50个像素点标记为1,其余像素点标记为0而得到;2b3)提取排列能量比特征:<img file="FDA0000064658510000024.GIF" wi="416" he="210" />其中k个最亮的像素点取为R<sub>tg</sub>中前5%强度值最大的像素点;2b4)提取基于协方差矩阵描述符的新特征:<img file="FDA0000064658510000025.GIF" wi="504" he="62" />其中,||g||<sub>2</sub>表示欧氏范数,C<sub>R</sub>是使用R<sub>tg</sub>中的像素点值计算得到的协方差矩阵,<img file="FDA0000064658510000026.GIF" wi="36" he="44" />是训练集、验证集以及测试集中所有候选切片对应的协方差矩阵的均值,<img file="FDA0000064658510000027.GIF" wi="142" he="56" />表示均值<img file="FDA0000064658510000028.GIF" wi="36" he="44" />处的对数映射,<img file="FDA0000064658510000029.GIF" wi="157" he="56" />表示均值<img file="FDA00000646585100000210.GIF" wi="36" he="44" />处的向量算子;2c)根据已提取的四个特征,生成特征向量f=[f<sub>1</sub> f<sub>2</sub> f<sub>3</sub> f<sub>4</sub>]<sup>T</sup>;2d)基于特征向量f,使用K近邻鉴别器计算鉴别统计量Dt(f),并将Dt(f)与根据验证样本集设定的鉴别门限Tt相比,当Dt(f)≤Tt时,特征向量f对应的切片图像则认为是真实目标切片,否则为虚警。
地址 710071 陕西省西安市太白南路2号