发明名称 一种基于图像质量最优的SAR多普勒模糊数估计方法
摘要 本发明涉及一种基于图像质量最优的合成孔径雷达(SAR)多普勒模糊数估计方法,属于SAR成像技术领域。本发明利用图像质量与多普勒模糊数估计偏差之间的函数关系,将使图像质量最好时对应的多普勒模糊数作为SAR多普勒模糊数的估计值。本发明采用了基于图像质量最优的SAR多普勒模糊数估计方法,由于SAR图像质量对多普勒模糊数估计误差十分敏感,并且该方法充分利用了图像质量与多普勒模糊数之间的定量关系,因此,对于低对比度和高对比度场景均具有较高的估计精度。
申请公布号 CN102520404A 申请公布日期 2012.06.27
申请号 CN201110390786.9 申请日期 2011.11.30
申请人 北京理工大学 发明人 丁泽刚;曾涛;吕争;刘荦锶;陈亮
分类号 G01S13/90(2006.01)I;G01S7/41(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 代理人
主权项 1.一种图像质量最优的SAR多普勒模糊数估计方法,其特征在于:其具体步骤如下:步骤一、按目标合成孔径时间长度获取距离压缩后的SAR回波数据,采用时域相关法得到多普勒中心小数部分估计值,用f<sub>dc_base</sub>表示;<maths num="0001"><![CDATA[<math><mrow><msub><mi>f</mi><mrow><mi>dc</mi><mo>_</mo><mi>base</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>N</mi><mi>r</mi></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>r</mi></msub></munderover><mo>{</mo><mfrac><msub><mi>F</mi><mi>a</mi></msub><mrow><mn>2</mn><mi>&pi;k</mi></mrow></mfrac><mo>&CenterDot;</mo><mi>arg</mi><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>a</mi></msub></munderover><mi>S</mi><mrow><mo>(</mo><msub><mi>&tau;</mi><mi>m</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><msup><mi>S</mi><mo>*</mo></msup><mrow><mo>(</mo><msub><mi>&tau;</mi><mi>m</mi></msub><mo>,</mo><mi>t</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,N<sub>r</sub>为距离压缩后的SAR回波数据的距离向采样点数,F<sub>a</sub>为脉冲重复频率,k为相关间隔,arg(·)为取复角操作,N<sub>a</sub>为距离压缩后的SAR回波数据的方位向采样点数,S(τ<sub>m</sub>,t)为距离压缩后的SAR回波数据中第m个距离单元的回波信号,1≤m≤Nr,τ<sub>m</sub>为距离压缩后的SAR回波数据中第m个距离单元的回波信号对应的距离向快时间,t为方位向慢时间;步骤二、根据SAR几何关系,通过公式(2)得到多普勒模糊数的粗估计值,用符号M<sub>amb0</sub>表示;然后,根据多普勒模糊数的粗估计值M<sub>amb0</sub>和预先设定的多普勒模糊数的估计偏差,用符号M<sub>ε</sub>表示,得到多普勒模糊数的搜索区间,多普勒模糊数用符号M<sub>amb</sub>表示,则多普勒模糊数的搜索区间为[M<sub>amb0</sub>-M<sub>ε</sub>,M<sub>amb0</sub>+M<sub>ε</sub>],即M<sub>amb</sub>∈[M<sub>amb0</sub>-M<sub>ε</sub>,M<sub>amb0</sub>+M<sub>ε</sub>];<maths num="0002"><![CDATA[<math><mrow><msub><mi>M</mi><mrow><mi>amb</mi><mn>0</mn></mrow></msub><mo>=</mo><mi>round</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>V</mi><mi>cos</mi><mi>&theta;</mi></mrow><mrow><mi>&lambda;</mi><msub><mi>F</mi><mi>a</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,round(·)为取整操作,V为载机速度,θ为波束中心对应的前斜角,λ为SAR的波长;步骤三、利用步骤一和步骤二得到的结果,根据公式(3),依次得到步骤二所述搜索区间[M<sub>amb0</sub>-M<sub>ε</sub>,M<sub>amb0</sub>+M<sub>ε</sub>]内的所有多普勒模糊数对应的多普勒中心,用符号f<sub>dc</sub>(M<sub>amb</sub>)表示;f<sub>dc</sub>(M<sub>amb</sub>)=f<sub>dc_base</sub>+M<sub>amb</sub>F<sub>a</sub>                            (3)步骤四、利用步骤三得到的多普勒中心f<sub>dc</sub>(M<sub>amb</sub>),对距离压缩后的SAR回波数据进行距离走动校正和方位压缩处理得到SAR图像,并计算其图像熵;步骤五、取步骤四得到的所有SAR图像熵的最小值,并将其对应的多普勒模糊数M<sub>amb</sub>作为距离压缩后的SAR回波数据的多普勒模糊数估计值,用符号<img file="FDA0000114311930000021.GIF" wi="105" he="55" />表示;<img file="FDA0000114311930000022.GIF" wi="1852" he="113" />经过上述步骤即可得到距离压缩后的SAR回波数据的基于图像质量最优的SAR多普勒模糊数估计值。
地址 100081 北京市海淀区中关村南大街5号