[发明专利]一种幂律模型InSAR对流层延迟改正方法

权利要求书

1.一种基于ERA-5再分析资料的幂律模型InSAR对流层延迟改正方法，其特征在于，通过以下步骤实现：

S1：收集N幅SAR单视复数影像的数据集，选取其中一幅作为主影像，其余的N-1幅影像作为辅影像；

S2：将N-1幅辅影像与主影像进行影像配准与重采样，预处理后，生成n幅干涉图；利用外部的DEM进行处理后，去除地形相位，最终生成差分干涉图；

S3：进行相干点的选取，筛选出高相干目标点；

S4：将所选的相干点目标进行相位解缠，得到每一个相干点解缠后的差分干涉相位；

S5：利用ERA-5气象再分析资料数据集进行对流层延迟改正，计算研究区对流层的干延迟和湿延迟；

S6：分析ERA-5再分析资料数据集改正的对流层延迟相位与地形高程之间的关系，用得到的对流层延迟相位和高程的关系拟合出所需要的幂律参数α和参考高程h_c，最后利用幂律模型进行对流层延迟改正；

S7：经过改正垂直分层延迟后，生成InSAR时间序列结果。

2.根据权利要求1所述一种基于ERA-5再分析资料的幂律模型InSAR对流层延迟改正方法，其特征在于，所述S5中的具体技术方法为：

采用表达式：计算研究区内干延迟导致的大气折射率，计算研究区内湿延迟导致大气的折射率，计算对流层大气折射率N；其中，K₁＝0.776Kpa^-1，K₂＝0.716Kpa^-1，K₃＝3.75*10⁵K²pa^-1，P为总气压，T为温度，e为水蒸气分压；

通过对上述折射率积分得到对流层干延迟和湿延迟：

采用表达式：计算研究区内雷达视线向对流层干延迟，

计算研究区内雷达视线向对流层湿延迟，

计算研究区内雷达视线向对流层总延迟，计算干涉图的对流层延迟；其中，为伪距增加到相位延迟的转换因子，θ为入射角，h_top为对流层天顶，h为地表高程。

3.根据权利要求1所述一种基于ERA-5再分析资料的幂律模型InSAR对流层延迟改正方法，其特征在于，所述S6中采用表达式：K为估计对流层延迟相位与参考高程之间的尺度因子；为常数偏差项，是高度h_c时的延迟相位。