[发明专利]一种干涉相位仿真方法

摘要

本发明提供一种干涉相位仿真方法，根据星载SAR成像几何在局部范围内变化缓慢的现象，通过将星载SAR测绘带内成像几何规律等效为n阶多项式描述的目标点的地理坐标系中的经纬度值随高程值的变化规律，生成多项式系数，能够以极小的相位精度损失换来计算复杂度的进一步降低，简化地理编码的运算过程，提高了运算效率，节约了运算时间。

主权项

一种一种干涉相位仿真方法，其特征在于，所述方法包括：步骤1：分别利用设定的n个假定高程h_0，1…h_0，n，生成主卫星雷达图像坐标系中栅格点的地理坐标系中的前向地理编码多项式系数{c_k}和{d_k}，其中，$\left\{\begin{array}{c}\mathrm{\ζ}=c_0+c_1{h}_{0,l}+\·\·\·+c_n{h}_{0,l}^n\\ \mathrm{\ξ}=d_0+d_1{h}_{0,l}+\·\·\·+d_n{h}_{0,l}^n\end{array}\right.,$(ζ，ξ)为栅格点对应的经纬度值，1≤l≤n；步骤2：根据目标点的主卫星雷达图像坐标计算成像时刻的主卫星位置矢量S₁；步骤3：设定成像时刻目标点高程值h₀的值；步骤4：通过步骤1得到的前向地理编码多项式系数{c_k}和{d_k}，得到高程值h₀对应的目标点的地理坐标系中的经纬度值(ζ₀，ξ₀)；步骤：5：对DEM模型h＝z(ζ，ξ)进行内插，得到目标点对应的经纬度值(ζ₀，ξ₀)在DEM模型中的实际高程值h，并判断实际高程值h和设定高程值h₀之差是否小于门限，若小于门限，则DEM模型匹配成功，进入步骤6，否则重新设定成像时刻目标点高程值h₀的值，并返回步骤4；步骤6：根据DEM匹配成功时的实际高程值h，根据主卫星轨道参数、主卫星位置矢量S₁，通过前向地理编码方法计算目标点的笛卡尔坐标；步骤7：根据步骤6得到的目标点的笛卡尔坐标，以及从卫星轨道参数，采用后向地理编码计算目标点在从图像中的坐标，并由此获得目标点在成像时刻所对应的从卫星的位置矢量S₂；步骤8：根据主卫星的位置S₁、从卫星的位置S₂以及目标点的笛卡尔坐标，计算目标点对主从天线的波程差，得到干涉相位值。