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

说明书

技术领域

本发明涉及干涉相位图仿真领域，特别涉及一种干涉相位仿真方法。

背景技术

干涉相位的模拟是基于干涉相位与两次成像的路径差的比例关系来进行的。在已知DEM数据以及卫星轨道参数的前提下，干涉相位仿真的关键在于寻找目标点在雷达成像时刻的主从卫星的位置、以及地理坐标位置之间三者之间的对应点，从而可根据路径差及雷达波长求出干涉相位。

发明内容

有鉴于此，本发明提供一种新的干涉相位仿真方法，通过将星载SAR测绘带内成像几何规律等效为n阶多项式描述的目标点的地理坐标系中的经纬度值随高程值的变化规律，从而简化地理编码的运算过程，提高了运算效率，节约了运算时间。

本发明提供了一种干涉相位仿真方法，包括：

步骤1：分别利用设定的n个假定高程h_0，1…h_0，n，生成主卫星雷达图像坐标系中栅格点的地理坐标系中的前向地理编码多项式系数{c_k}和{d_k}，其中，

ζ=c0+c1h0,l+...+cnh0,lnξ=d0+d1h0,l+...+dnh0,ln,]]>(ζ，ξ)为栅格点对应的经纬度值，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₂；