| 主权项 |
一种一种干涉相位仿真方法，其特征在于，所述方法包括：步骤1：分别利用设定的n个假定高程h<sub>0，1</sub>…h<sub>0，n</sub>，生成主卫星雷达图像坐标系中栅格点的地理坐标系中的前向地理编码多项式系数{c<sub>k</sub>}和{d<sub>k</sub>}，其中，<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>ζ</mi><mo>=</mo><msub><mi>c</mi><mn>0</mn></msub><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><msub><mi>h</mi><mrow><mn>0</mn><mo>,</mo><mi>l</mi></mrow></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>c</mi><mi>n</mi></msub><msubsup><mi>h</mi><mrow><mn>0</mn><mo>,</mo><mi>l</mi></mrow><mi>n</mi></msubsup></mtd></mtr><mtr><mtd><mi>ξ</mi><mo>=</mo><msub><mi>d</mi><mn>0</mn></msub><mo>+</mo><msub><mi>d</mi><mn>1</mn></msub><msub><mi>h</mi><mrow><mn>0</mn><mo>,</mo><mi>l</mi></mrow></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>d</mi><mi>n</mi></msub><msubsup><mi>h</mi><mrow><mn>0</mn><mo>,</mo><mi>l</mi></mrow><mi>n</mi></msubsup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FSA0000119652280000011.GIF" wi="764" he="242" /></maths>(ζ，ξ)为栅格点对应的经纬度值，1≤l≤n；步骤2：根据目标点的主卫星雷达图像坐标计算成像时刻的主卫星位置矢量S<sub>1</sub>；步骤3：设定成像时刻目标点高程值h<sub>0</sub>的值；步骤4：通过步骤1得到的前向地理编码多项式系数{c<sub>k</sub>}和{d<sub>k</sub>}，得到高程值h<sub>0</sub>对应的目标点的地理坐标系中的经纬度值(ζ<sub>0</sub>，ξ<sub>0</sub>)；步骤：5：对DEM模型h＝z(ζ，ξ)进行内插，得到目标点对应的经纬度值(ζ<sub>0</sub>，ξ<sub>0</sub>)在DEM模型中的实际高程值h，并判断实际高程值h和设定高程值h<sub>0</sub>之差是否小于门限，若小于门限，则DEM模型匹配成功，进入步骤6，否则重新设定成像时刻目标点高程值h<sub>0</sub>的值，并返回步骤4；步骤6：根据DEM匹配成功时的实际高程值h，根据主卫星轨道参数、主卫星位置矢量S<sub>1</sub>，通过前向地理编码方法计算目标点的笛卡尔坐标；步骤7：根据步骤6得到的目标点的笛卡尔坐标，以及从卫星轨道参数，采用后向地理编码计算目标点在从图像中的坐标，并由此获得目标点在成像时刻所对应的从卫星的位置矢量S<sub>2</sub>；步骤8：根据主卫星的位置S<sub>1</sub>、从卫星的位置S<sub>2</sub>以及目标点的笛卡尔坐标，计算目标点对主从天线的波程差，得到干涉相位值。 |
