[发明专利]机载干涉合成孔径雷达的干涉相位仿真方法

权利要求书

1.一种机载干涉合成孔径雷达的干涉相位仿真方法，包括：

步骤A：在地距坐标系下，设置仿真场景的DEM以及仿真的机载InSAR系统参数；

步骤B：根据仿真场景的DEM和机载InSAR系统参数，判断该仿真场景中的叠掩区域、阴影区域和正常区域，用标记函数进行标记；

步骤C：在每一个方位向上，分别求出叠掩区域、阴影区域和正常区域的理想干涉相位；

步骤D：设置仿真场景的相干系数值γ；

步骤E：对每一个方位向上的非阴影区域的距离门，根据相干系数值γ和理想干涉相位值φ_ideal生成有噪声的仿真干涉相位φ_simu；步骤E包括：

步骤E1：根据相干系数值γ计算均值为0时的干涉相位概率密度函数为：

-π＜φ＜π

步骤E2：根据干涉相位概率密度函数计算其分布函数F(φ)：

步骤E3：计算F(φ)的反函数G(Φ)；

步骤E4：利用Matlab中的rand函数生成服从区间[0，1]上均匀分布的随机数Φ，将其代入反函数G(Φ)，由此生成干涉相位噪声为：

φ_noise＝G(Φ)

步骤E5：根据理想的干涉相位φ_ideal和干涉相位噪声φ_noise，生成仿真的干涉相位为：

φ_simu＝wrap{φ_ideal+φ_noise}

其中，wrap{}为干涉相位缠绕算子。

2.根据权利要求1所述的机载干涉合成孔径雷达的干涉相位仿真方法，所述步骤C包括：

步骤C1：对判断为阴影区域的距离门，设置其干涉相位为服从区间(-π，π]上均匀分布的随机噪声：

φ_shadow＝2π·rand-π

其中，rand为服从区间[0，1]上均匀分布的随机数，利用Matlab中的rand函数生成；

步骤C2：对判断为正常区域的距离门，根据机载InSAR成像几何关系，计算出该距离门的理想干涉相位φ_ideal；

步骤C3：对判断为叠掩区域的距离门，根据机载InSAR成像几何关系，计算出n个目标点对应的理想干涉相位为φ_{ideal_1}，...，φ_{ideal_n}，并通过叠加方式得到该距离门的理想干涉相位φ_ideal，其中n为投影到该距离门内的目标点个数；所述步骤C3中，通过叠加方式得到该距离门的理想干涉相位φ_ideal的方法为：

判断n个目标点中满足的目标点为k₁，...，k_l，l为n个目标点中位于迎坡面上的目标点个数，其中，为n个目标点处斜距R₁关于地距r_g的偏导数，则该距离门叠加后的理想干涉相位为：

其中，c_p为迎坡面上目标点k_p的干涉相位加权值，其介于2至5之间。

3.根据权利要求1所述的机载干涉合成孔径雷达的干涉相位仿真方法，所述步骤A中：

仿真场景的DEM根据干涉相位仿真的目的选择采用实际场景的DEM，或通过计算机仿真的手段生成；

所述仿真的机载InSAR系统参数包括：波长、斜距分辨率、主天线高度、基线长度和基线角。