发明名称 双通道调频连续波SAR系统的动目标检测成像方法
摘要 双通道调频连续波SAR系统的动目标检测成像方法,雷达系统录取回波数据,对回波数据进行去除残余视频相位处理,进而对数据进行多普勒频移补偿,同时对二通道变换数据进行相位偏差补偿,并对两通道数据进行消减处理;对杂波相消后的图像进行徙动校正及方位去斜处理,实现目标粗成像,检测运动目标并逐个提取;将提取的目标反变换到原始数据域,估计运动参数来构造动目标精确的方位去斜函数,进行去除残余视频相位处理及方位精确去斜处理及距离徙动校正,构造模糊函数进行补偿,然后进行Keystone变换,完成动目标的精确成像。本发明方法有效解决了连续波体制下运动参数未知时的快速动目标成像问题,极大地提高信杂噪比和目标检测概率。
申请公布号 CN103744068A 申请公布日期 2014.04.23
申请号 CN201410026200.4 申请日期 2014.01.21
申请人 西安电子科技大学 发明人 周峰;王金伟;梁舒幸;王海兵
分类号 G01S7/41(2006.01)I;G01S7/36(2006.01)I;G01S13/89(2006.01)I;G01S13/50(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 广东秉德律师事务所 44291 代理人 杨焕军
主权项 1.双通道调频连续波SAR系统的动目标检测成像方法,其特征在于,包括以下步骤:雷达系统录取回波数据,包括一通道回波数据<img file="FDA0000459527780000011.GIF" wi="170" he="84" />和二通道回波数据<img file="FDA0000459527780000012.GIF" wi="208" he="87" />步骤1、对回波数据进行杂波抑制及动目标检测;步骤1a、对一通道回波数据<img file="FDA0000459527780000013.GIF" wi="182" he="84" />和二通道回波数据<img file="FDA0000459527780000014.GIF" wi="182" he="88" />分别进行距离向傅里叶变换及去除残余视频相位处理,然后变换回二维时域,得到一通道变换数据<img file="FDA0000459527780000015.GIF" wi="246" he="100" />和二通道变换数据<img file="FDA0000459527780000016.GIF" wi="261" he="95" />步骤1b、将一通道变换数据<img file="FDA0000459527780000017.GIF" wi="250" he="97" />和二通道变换数据<img file="FDA0000459527780000018.GIF" wi="238" he="95" />分别进行方位向傅里叶变换及多普勒频移补偿,同时对二通道变换数据进行相位偏差补偿,然后分别进行方位向逆傅里叶变换,得到一通道补偿数据S<sub>if1</sub>和二通道补偿数据S<sub>if2</sub>;步骤1c、将一通道补偿数据S<sub>if1</sub>和二通道补偿数据S<sub>if2</sub>相减,得到杂波相消后的图像S<sub>12</sub>;步骤1d、对杂波相消后的图像S<sub>12</sub>进行徙动校正,得到校正图像S<sub>12r</sub>;步骤1e、将校正图像S<sub>12r</sub>变换到距离频域并进行方位去斜处理,然后进行方位向傅里叶变换,得到目标图像S<sub>12l</sub>;步骤1f、对目标图像S<sub>12l</sub>进行运动目标检测,对动目标进行逐个提取,得到与检测到的目标对应的目标检测图像S<sub>i</sub>,i=1,…,M,M为检测到的动目标个数,记录目标检测图像S<sub>i</sub>中动目标的距离位置R<sub>i</sub>;步骤2、对动目标逐个进行精确成像;步骤2a、将第i个动目标的目标检测图像S<sub>i</sub>逆变换到原始数据域,得到原始数据域图像S<sub>mi</sub>,<maths num="0001"><![CDATA[<math><mrow><msub><mi>S</mi><mi>mi</mi></msub><mo>=</mo><msub><mi>IFFT</mi><msub><mi>f</mi><mi>r</mi></msub></msub><mrow><mo>(</mo><msub><mi>FFT</mi><msub><mi>f</mi><mi>r</mi></msub></msub><mrow><mo>(</mo><msub><mi>IFFT</mi><msub><mi>f</mi><mi>a</mi></msub></msub><mrow><mo>(</mo><msub><mi>FFT</mi><msub><mi>f</mi><mi>a</mi></msub></msub><mrow><mo>(</mo><msub><mi>IFFT</mi><msub><mi>f</mi><mi>r</mi></msub></msub><mrow><mo>(</mo><msub><mi>IFFT</mi><msub><mi>f</mi><mi>a</mi></msub></msub><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>*</mo><mi>conj</mi><mrow><mo>(</mo><msub><mi>H</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>*</mo><mi>conj</mi><mrow><mo>(</mo><msub><mi>H</mi><mi>RMC</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>*</mo><mi>conj</mi><mrow><mo>(</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>*</mo><mi>conj</mi><mrow><mo>(</mo><msub><mi>H</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>conj(·)表示取共轭运算;步骤2b、估计原始数据域图像S<sub>mi</sub>中动目标的多普勒中心偏移量f<sub>dci</sub>,计算动目标的基带内径向速度<img file="FDA0000459527780000019.GIF" wi="250" he="128" />f<sub>dci</sub>表示第i个动目标的多普勒中心偏移量,λ为信号波长;步骤2c、将原始数据域图像S<sub>mi</sub>变换到距离频域,得到距离频域图像S<sub>mri</sub>,估计距离频域图像S<sub>mri</sub>中动目标的斜线斜率K<sub>i</sub>,计算动目标的速度模糊数N<sub>i</sub>;步骤2d、利用动目标的基带径向速度v<sub>bi</sub>、动目标的速度模糊数N<sub>i</sub>以及动目标的距离位置R<sub>i</sub>构造该动目标的方位去斜函数H<sub>am</sub>,对距离频域图像S<sub>mri</sub>进行去除残余视频相位处理后,通过方位去斜函数H<sub>am</sub>对距离频域图像S<sub>mri</sub>进行方位精确去斜处理,然后变换至二维时域,得到二维时域图像S<sub>mrci</sub>;步骤2e、将二维时域图像S<sub>mrci</sub>变换至方位频域,进行多普勒频移补偿,然后变换回方位时域,得到方位时域图像S<sub>mrmi</sub>;步骤2f、对方位时域图像S<sub>mrmi</sub>进行距离徙动校正,并构造模糊函数H<sub>blur</sub>进行补偿,然后进行Keystone变换,最后变换到二维频域,得到第i个动目标的目标成像图S<sub>fi</sub>,完成动目标的精确成像;<maths num="0002"><![CDATA[<math><mrow><msub><mi>S</mi><mi>fi</mi></msub><mo>=</mo><msub><mi>FFT</mi><msub><mi>f</mi><mi>a</mi></msub></msub><mrow><mo>(</mo><msub><mi>FFT</mi><msub><mi>f</mi><mi>r</mi></msub></msub><mrow><mo>(</mo><mi>Keystone</mi><mrow><mo>(</mo><msub><mi>S</mi><mi>mrmi</mi></msub><mo>*</mo><msub><mi>H</mi><mi>RMC</mi></msub><mo>*</mo><msub><mi>H</mi><mi>blur</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>H</mi><mi>blur</mi></msub><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mfrac><mrow><mn>4</mn><mi>&pi;&gamma;</mi><msub><mi>N</mi><mi>i</mi></msub><msub><mi>&lambda;R</mi><mi>ref</mi></msub><mi>PRF</mi></mrow><msup><mi>c</mi><mn>2</mn></msup></mfrac><msub><mi>t</mi><mi>m</mi></msub><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><msub><mrow><mn>2</mn><mi>N</mi></mrow><mi>i</mi></msub><mi>&pi;PRF</mi></mrow><msub><mi>f</mi><mi>c</mi></msub></mfrac><mrow><mo>(</mo><mi>&gamma;</mi><mover><mi>t</mi><mo>^</mo></mover><mo>+</mo><msub><mi>f</mi><mi>c</mi></msub><mo>)</mo></mrow><msub><mi>t</mi><mi>m</mi></msub><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>N</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mi>PRF</mi><mover><mi>t</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中,R<sub>ref</sub>为参考距离,Keystone(·)表示进行Keystone变换;重复步骤2a至2f,直至完成对检测到的所有动目标的精确成像,以上步骤中的i=1,…,M。
地址 710068 陕西省西安市太白南路2号
您可能感兴趣的专利