主权项 |
基于误差等效RD模型修正的SAR几何校正方法,其特征在于,包括以下步骤:第一步:建立基于误差等效RD模型修正的距离‑多普勒方程,如下:距离方程: <mrow> <msub> <mi>F</mi> <mrow> <mn>1</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>ΔR</mi> <mo>,</mo> <mi>p</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>R</mi> <mi>near</mi> </msub> <mo>+</mo> <mi>Δ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>δ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>ΔR</mi> <mo>,</mo> <mi>p</mi> <mo>,</mo> <mi>Δ</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>Δ</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>λ</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>Δ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>δ</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的几何校正。 |