发明名称 一种基于信号子空间技术的单通道星载SAR动目标检测方法
摘要 本发明提供了一种基于信号子空间技术的单通道星载SAR动目标检测方法,能够有效抑制残余杂波,降低虚警概率。第一步:将参考图像矩阵进行奇异值分解;第二步:求参考图像杂波子空间的标准正交基;第三步:求解待检测图像在参考图像杂波子空间的投影;根据参考图像杂波子空间的标准正交基求出f2(x,y)在f1(x,y)子空间的投影;第四步:求解待检测图像杂波抑制后的图像;第五步:用参考图像块的均值将杂波抑制后的图像归一化。
申请公布号 CN103336271A 申请公布日期 2013.10.02
申请号 CN201310098718.4 申请日期 2013.03.26
申请人 中国科学院电子学研究所 发明人 张丹丹;许睿鹏;胡东辉;仇晓兰;雷斌;付琨;丁赤飚
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 高燕燕
主权项 1.一种基于信号子空间技术的单通道星载SAR动目标检测方法,其特征在于,包括以下步骤:第一步:将参考图像矩阵进行奇异值分解;设滤波窗口大小为(2n<sub>x</sub>+1)·(2n<sub>y</sub>+1),图像大小为L<sub>x</sub>·L<sub>y</sub>,<maths num="0001"><![CDATA[<math><mrow><msub><mi>F</mi><mn>1</mn></msub><mo>&Element;</mo><msup><mi>C</mi><mrow><msub><mi>L</mi><mi>x</mi></msub><msub><mi>L</mi><mi>y</mi></msub><mo>&times;</mo><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>n</mi></mrow><mi>x</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>n</mi></mrow><mi>y</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msup><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>F</mi><mn>2</mn></msub><mo>&Element;</mo><msup><mi>C</mi><mrow><msub><mi>L</mi><mi>x</mi></msub><msub><mi>L</mi><mi>y</mi></msub><mo>&times;</mo><mn>1</mn></mrow></msup></mrow></math>]]></maths>分别是f<sub>1</sub>(x,y)和f<sub>2</sub>(x,y)得所有像素组成的数据矩阵,即:F<sub>1</sub>={[f<sub>1</sub>(x<sub>i</sub>-m,y<sub>j</sub>-n),m=-n<sub>x</sub>...n<sub>x</sub>,n=-n<sub>y</sub>...n<sub>y</sub>],i=1,...,L<sub>x</sub>,j=1,...,L<sub>y</sub>}F<sub>2</sub>={f<sub>2</sub>(x<sub>i</sub>,y<sub>j</sub>),i=1,...,L<sub>x</sub>,j=1,...,L<sub>y</sub>}然后对F<sub>1</sub>进行奇异值分解得到:F<sub>1</sub>(x,y)=USV<sup>H</sup>;第二步:求参考图像杂波子空间的标准正交基;上述奇异值分解后,正交矩阵U的前(2n<sub>x</sub>+1)·(2n<sub>y</sub>+1)列为杂波子空间的标准正交基,记为<img file="FDA00002964190100013.GIF" wi="546" he="109" />第三步:求解待检测图像在参考图像杂波子空间的投影;根据参考图像杂波子空间的标准正交基求出f<sub>2</sub>(x,y)在f<sub>1</sub>(x,y)子空间的投影,即:<img file="FDA00002964190100014.GIF" wi="478" he="134" />第四步:求解待检测图像杂波抑制后的图像;将<img file="FDA00002964190100015.GIF" wi="44" he="67" />重排得到<img file="FDA00002964190100016.GIF" wi="236" he="98" />则f<sub>2</sub>(x,y)相对于f<sub>1</sub>(x,y)变化的图像为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>f</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>f</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>f</mi><mo>^</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>第五步:用参考图像块的均值将杂波抑制后的图像归一化,归一化后的图像信号为:<maths num="0004"><![CDATA[<math><mrow><msub><mi>D</mi><mrow><mi>c</mi><mo>_</mo><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>f</mi><mrow><mi>c</mi><mo>_</mo><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mrow><mi>E</mi><mrow><mo>(</mo><msub><mi>f</mi><mrow><mn>1</mn><mo>_</mo><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>其中,f<sub>c_k</sub>(x,y)表示信号子空间处理后的第k块变化图像,E(·)表示均值。
地址 100190 北京市海淀区北四环西路19号