发明名称 一种基于误差等效RD模型修正的SAR几何校正方法
摘要 本发明提供了一种基于误差等效RD模型修正的SAR几何校正方法,仅需很少的控制点,且对控制点的均匀分布特性要求不高;而且采用常规的数值求解方法即可得到很好的结果,无需迭代。第一步:建立基于误差等效RD模型修正的距离-多普勒方程,第二步:对第一步中得到的方程进行泰勒级数一级展开得到误差方程;第三步:采用最小二乘的方法求解第二步得到的误差方程;第四步:求得误差方程的解后,采用修正后的距离多普勒方程对每个像素进行几何定位,即实现了SAR的几何校正。
申请公布号 CN103235304B 申请公布日期 2014.11.05
申请号 CN201310098477.3 申请日期 2013.03.26
申请人 中国科学院电子学研究所 发明人 仇晓兰;刘佳音;尤红建;胡玉新;雷斌;付琨;丁赤飚
分类号 G01S13/90(2006.01)I;G01S7/40(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 高燕燕
主权项 基于误差等效RD模型修正的SAR几何校正方法,其特征在于,包括以下步骤:第一步:建立基于误差等效RD模型修正的距离‑多普勒方程,如下:距离方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>F</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>&Delta;R</mi><mo>,</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>R</mi><mi>near</mi></msub><mo>+</mo><mi>&Delta;R</mi><mo>+</mo><mi>p</mi><msub><mi>N</mi><mi>i</mi></msub><mo>+</mo><msub><mi>N</mi><mi>i</mi></msub><mi>&delta;r</mi><mo>-</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>si</mi></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Y</mi><mi>si</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Z</mi><mi>si</mi></msub><mo>-</mo><msub><mi>z</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA0000555440760000011.GIF" wi="1647" he="99" /></maths>多普勒方程:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>F</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>&Delta;R</mi><mo>,</mo><mi>p</mi><mo>,</mo><mi>&Delta;</mi><msub><mi>f</mi><mi>d</mi></msub><mo>,</mo><mi>q</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>f</mi><mi>d</mi></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>f</mi><mi>d</mi></msub><mo>+</mo><mi>q</mi><msub><mi>M</mi><mi>i</mi></msub><mo>+</mo><mfrac><mn>2</mn><mi>&lambda;</mi></mfrac><mfrac><mrow><msub><mi>V</mi><mi>sxi</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>si</mi></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>V</mi><mi>syi</mi></msub><mrow><mo>(</mo><msub><mi>Y</mi><mi>si</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>V</mi><mi>szi</mi></msub><mrow><mo>(</mo><msub><mi>Z</mi><mi>si</mi></msub><mo>-</mo><msub><mi>z</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>R</mi><mi>near</mi></msub><mo>+</mo><mi>&Delta;R</mi><mo>+</mo><mi>p</mi><msub><mi>N</mi><mi>i</mi></msub><mo>+</mo><msub><mi>N</mi><mi>i</mi></msub><mi>&delta;r</mi></mrow></mfrac><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA0000555440760000012.GIF" wi="1708" he="148" /></maths>其中R<sub>near</sub>表示图像距离向起始像素对应的斜距,即近距;δr为图像的距离向像素间隔;λ为波长;i=1,...,K为控制点号,K为控制点总数;N<sub>i</sub>为第i个控制点的距离向像素号;M<sub>i</sub>是第i个控制点的方位向像素号;[x<sub>i</sub>,y<sub>i</sub>,z<sub>i</sub>]<sup>T</sup>为第i个控制点的三维位置坐标;f<sub>d</sub>为成像处理所采用的多普勒中心频率;[X<sub>si</sub>,Y<sub>si</sub>,Z<sub>si</sub>]<sup>T</sup>为SAR卫星与第i个控制点相对位移引起的多普勒频率为f<sub>d</sub>时,卫星的位置矢量;[V<sub>sxi</sub>,V<sub>syi</sub>,V<sub>szi</sub>]<sup>T</sup>为SAR卫星与第i个控制点相对位移引起的多普勒频率为f<sub>d</sub>时,卫星与目标的相对速度矢量;ΔR,p,Δf<sub>d</sub>,q为四个待定的误差等效模型参数;此处称ΔR+pN<sub>i</sub>为距离修正量,ΔR和p为距离模型修正参数;称Δf<sub>d</sub>+qM<sub>i</sub>为多普勒中心频率修正量,Δf<sub>d</sub>和q为多普勒模型修正参数;第二步:对第一步中得到的方程进行泰勒级数一级展开得到误差方程;第三步:采用最小二乘的方法求解第二步得到的误差方程;第四步:求得误差方程的解后,采用修正后的距离多普勒方程对每个像素进行几何定位,即实现了SAR的几何校正。
地址 100190 北京市海淀区北四环西路19号