[发明专利]基于InSAR多维形变信息的黄土滑坡类型及滑动模式分析方法

权利要求书

1.一种基于InSAR多维形变信息的黄土滑坡类型及滑动模式分析方法，其特征在于，所述的方法包括以下步骤：

步骤1，通过覆盖研究区域的T个SAR传感器采集的数据，获取滤波InSAR差分干涉图；

步骤2，对滤波InSAR差分干涉图进行相位解缠，将解缠后的滤波InSAR差分干涉图编码至地理坐标系，并采样至相同的坐标格网和空间分辨率，得到地理坐标系下的解缠干涉图；

步骤3，对地理坐标系下的解缠干涉图，计算研究区域的多维地表形变速率，并计算研究区域的多维累积地表形变；

针对于研究区域，按时间顺序选取多维累积地表形变构成多维地表形变时间序列；

步骤4，对获得的研究区域的多维地表形变速率及多维形变时间序列，结合研究区域的遥感影像、地形图进行滑坡形变机理分析，确定黄土滑坡的类型及滑动模式。

2.如权利要求1所述的基于InSAR多维形变信息的黄土滑坡类型及滑动模式分析方法，其特征在于，所述的获取滤波InSAR差分干涉图，包括：

通过覆盖研究区域的T个SAR传感器，获得覆盖研究区域的T个升降轨SAR影像数据集并利用无人机摄影测量方式获得外部DEM数据，对获取的T个升降轨SAR影像数据集中每一个SAR影像数据集的影像分别进行两两做差形成差分干涉图，然后将每张差分干涉图减去外部DEM数据，得到InSAR差分干涉图，对InSAR差分干涉图进行滤波，得到滤波InSAR差分干涉图。

3.如权利要求1所述的基于InSAR多维形变信息的黄土滑坡类型及滑动模式分析方法，其特征在于，所述的计算研究区域的多维地表形变速率，包括：

采用下式3计算研究区域的多维地表形变速率：

上式中，A¹到A^T分别表示1到T个升降轨影像数据集的系数矩阵，所述的多维地表变形速率包括V_N、V_E、V_U，其中V_N表示南北向的地表形变速率，V_E表示东西向的地表形变速率，V_U表示垂直方向的地表形变速率；到分别表示1到T个数据集的观测的相位值。

4.如权利要求1所述的基于InSAR多维形变信息的黄土滑坡类型及滑动模式分析方法，其特征在于，所述的计算研究区域的多维累积地表形变，包括：

采用下式4计算研究区域的多维累积地表形变：

上式中，表示T个SAR影像数据集总共的SAR影像数量，分别为南北、东西以及垂直向的累积地表形变，表示第在第t_i时刻地表的南北向形变速率；表示在第t_i时刻地表的东西向形变速率；表示在第t_i时刻地表的垂直向形变速率；Δt_i为t_i-1时刻与t_i时刻的时间间隔。