[发明专利]InSAR相位时序处理方法、电子设备及存储介质

权利要求书

1.一种InSAR相位时序处理方法，其特征在于，所述方法包括：

获取配准好的N+1幅SAR图像，其中，所述N+1幅SAR图像中的每幅图像存在至少一幅可干涉的SAR图像，N为大于1的整数；

从所述N+1幅SAR图像中选取目标干涉对，每个所述目标干涉对包括两幅可干涉且基线值小于等于预设基线阈值的SAR图像；

基于所述目标干涉对，形成基线子集，其中，每个所述基线子集包括可接连干涉的多幅SAR图像；

对每个所述基线子集中的每个目标干涉对进行差分干涉处理，得到差分干涉图；

根据预设的解缠算子对从所述差分干涉图选取的高相干像元进行相位解缠处理，得到解缠结果，所述解缠结果包括所述高相干像元解缠后的干涉相位；

对所述解缠结果进行滤波，得到滤波后的解缠结果；

基于所述滤波后的解缠结果，建立所述目标干涉对中两幅SAR图像的形变模型；

通过岭回归算法，对所述形变模型进行数据反演，得到所述目标干涉对中两幅SAR图像的形变结果。

2.根据权利要求1所述的方法，其特征在于，在获取配准好的N+1幅SAR图像之前，所述方法还包括：

根据时间序列依次获得N+1幅待配准SAR图像；

从所述N+1幅待配准SAR图像中选取任一图像作为主图像，其他图像为辅图像，其中，所述其他图像为所述N+1幅待配准SAR图像中除去所述主图像之外的图像；

对所述主图像和所述辅图像进行图像配准，得到N+1幅配准好的SAR图像。

3.根据权利要求1所述的方法，其特征在于，对每个所述基线子集中的每个目标干涉对进行差分干涉处理，得到差分干涉图，包括：

从所述基线子集中选择任一个目标干涉对作为当前干涉对；

对所述当前干涉对中的图像进行复共轭相乘，得到原始干涉图；

对所述原始干涉图进行平滑滤波，得到滤波后的原始干涉图；

基于DEM模拟的相位或从另一目标干涉对获取的地形相位，对所述滤波后的原始干涉图进行差分处理，得到所述当前干涉对的所述差分干涉图，其中，所述另一目标干涉对为从所述基线子集中选择的除去所述当前干涉对之外的干涉对；

遍历所有所述基线子集中的每个目标干涉以作为新的当前干涉对，并重复步骤对所述当前干涉对中的图像进行复共轭相乘，得到原始干涉图，至步骤基于DEM模拟的相位或从另一目标干涉对获取的地形相位，对所述滤波后的原始干涉图进行差分处理，以得到每个目标干涉对的差分干涉图。

4.根据权利要求1所述的方法，其特征在于，所述形变模型为：

BV＝δφ

其中，B指大小为M*N的矩阵，B矩阵的每行对应目标干涉对中的两幅图像之间的时间跨度，B矩阵的每一列对应时间相邻的两个SAR图像的同一像元点的时间间隔；V指待求解的形变相位速率；δφ指所述差分干涉图中的所述高相干像元解缠后的干涉相位。

5.根据权利要求4所述的方法，其特征在于，通过岭回归算法，对所述形变模型进行数据反演，得到所述目标干涉对中两幅SAR图像的形变结果，包括：

通过岭回归算法，对所述形变模型进行数据反演，得到所述形变相位速率的岭估计，表示为：

其中，k指岭参数；I指单位矩阵；

基于所述目标干涉对的相干点，通过L曲线法确定所述岭参数，并根据所述岭估计的计算公式确定所述目标干涉对的所述形变相位速率；

根据所述目标干涉对中两幅SAR图像的时间间隔和所述形变相位速率，确定所述目标干涉对中两幅SAR图像的形变量，并标记发生形变的图区的坐标，以作为所述形变结果。

6.根据权利要求1-5中任一项所述的方法，其特征在于，所述方法还包括：

当所述形变结果中的形变量超过预设阈值时，发出告警提示。

7.根据权利要求1-5中任一项所述的方法，其特征在于，所述基线值包括时间基线值和/或空间垂直基线值。