发明名称 一种基于误差等效RD模型修正的SAR几何校正方法
摘要 本发明提供了一种基于误差等效RD模型修正的SAR几何校正方法,仅需很少的控制点,且对控制点的均匀分布特性要求不高;而且采用常规的数值求解方法即可得到很好的结果,无需迭代。第一步:建立基于误差等效RD模型修正的距离-多普勒方程,第二步:对第一步中得到的方程进行泰勒级数一级展开得到误差方程;第三步:采用最小二乘的方法求解第二步得到的误差方程;第四步:求得误差方程的解后,采用修正后的距离多普勒方程对每个像素进行几何定位,即实现了SAR的几何校正。
申请公布号 CN103235304A 申请公布日期 2013.08.07
申请号 CN201310098477.3 申请日期 2013.03.26
申请人 中国科学院电子学研究所 发明人 仇晓兰;刘佳音;尤红建;胡玉新;雷斌;付琨;丁赤飚
分类号 G01S13/90(2006.01)I;G01S7/40(2006.01)I 主分类号 G01S13/90(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 高燕燕
主权项 基于误差等效RD模型修正的SAR几何校正方法,其特征在于,包括以下步骤:第一步:建立基于误差等效RD模型修正的距离‑多普勒方程,如下:距离方程: <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>多普勒方程: <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> <mo>+</mo> <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> <mrow> <mi>p</mi> <msub> <mi>N</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>N</mi> <mi>i</mi> </msub> <mi>&delta;</mi> </mrow> <mi>r</mi> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>;</mo> </mrow>其中Rnear表示图像距离向起始像素对应的斜距,即近距;δr为图像的距离向像素间隔;λ为波长;i=1,...,K为控制点号,K为控制点总数;Ni为第i个控制点的距离向像素号;Mi是第i个控制点的方位向像素号;[xi,yi,zi]T为第i个控制点的三维位置坐标;fd为成像处理所采用的多普勒中心频率;[Xsi,Ysi,Zsi]T为SAR卫星与第i个控制点相对位移引起的多普勒频率为fd时,卫星的位置矢量;[Vsxi,Vsyi,Vszi]T为SAR卫星与第i个控制点相对位移引起的多普勒频率为fd时,卫星与目标的相对速度矢量;ΔR,p,Δfd,q为四个待定的误差等效模型参数;此处称ΔR+pNi为距离修正量,ΔR和p为距离模型修正参数;称Δfd+qMi为多普勒修正量,Δfd和q为多普勒模型修正参数;第二步:对第一步中得到的方程进行泰勒级数一级展开得到误差方程;第三步:采用最小二乘的方法求解第二步得到的误差方程;第四步:求得误差方程的解后,采用修正后的距离多普勒方程对每个像素进行几何定位,即实现了SAR的几何校正。
地址 100190 北京市海淀区北四环西路19号