1.一种针对滑坡的InSAR形变监测结果解释方法,其特征在于,包括如下步骤:S1、获取待解释滑坡区域预设监测周期内的InSAR地表形变监测结果;所述InSAR地表形变监测结果包括:InSAR形变速率结果和各期次形变量结果;所述预设监测周期包括多个监测期次;S2、获取覆盖所述待解释滑坡区域的光学影像数据和数字高程数据;根据所述光学影像数据,分析所述待解释滑坡区域内的典型地物、植被界限、特征地形和平面形态,圈定所述待解释滑坡区域中滑坡体的初始边界;并根据所述InSAR形变速率结果对所述初始边界进行修正,确定所述滑坡体的边界范围;S3、对所述InSAR形变速率结果进行统计和区段划分,在所述边界范围内对所述滑坡体进行形变分级区域划分,生成形变分区结果;S4、根据所述各期次形变量结果,结合所述数字高程数据中的地形信息,获取所述滑坡体的空间演变规律;所述地形信息包括:坡度和坡向;S5、根据所述滑坡体的形变分区结果和空间演变规律,得到所述待解释滑坡区域内的重点监测点位信息。
2.如权利要求1所述的一种针对滑坡的InSAR形变监测结果解释方法,其特征在于,所述步骤S1还包括:采用地理编码的方式将所述InSAR地表形变监测结果转化为具有经纬度信息的结果并进行投影。
3.如权利要求2所述的一种针对滑坡的InSAR形变监测结果解释方法,其特征在于,所述步骤S2包括:S21、获取覆盖所述待解释滑坡区域的光学影像数据和数字高程数据;根据所述光学影像数据,分析所述待解释滑坡区域内的典型地物、植被界限、特征地形和平面形态,圈定所述待解释滑坡区域中滑坡体的初始边界;并对圈定的结果进行矢量化处理;S22、对矢量化后的初始边界图层进行投影转换,确保所述初始边界和所述InSAR形变速率结果的投影坐标系一致;S23、将相同坐标系下的所述InSAR形变速率结果和所述初始边界进行叠加,圈定所述初始边界外围预设临近区域的部分区域并形成矢量边界,确保所述部分区域和所述初始边界的投影坐标系一致;所述部分区域的InSAR形变速率绝对值大于预设数值;S24、将新圈定的矢量边界和所述初始边界进行合并,生成所述滑坡体的边界范围。
4.如权利要求1所述的一种针对滑坡的InSAR形变监测结果解释方法,其特征在于,所述步骤S3包括:S31、根据预设形变速率绝对值阈值,在所述边界范围内对所述滑坡体进行形变分级区域划分,得到不同形变区段对应的滑坡体区域;S32、对所述不同形变区段对应的滑坡体区域的边界进行矢量化处理,得到形变分级区域划分矢量图,生成形变分区结果。
5.如权利要求4所述的一种针对滑坡的InSAR形变监测结果解释方法,其特征在于,所述步骤S31包括:根据预设形变速率绝对值阈值K毫米/年和M毫米/年,分别以小于-M毫米/年、-M毫米/年至-K毫米/年、-K毫米/年至K毫米/年、K毫米/年至M毫米/年、大于M毫米/年为区间,在所述边界范围内对所述滑坡体进行形变分级区域划分,得到不同形变区段对应的滑坡体区域;其中,0<K<M。
6.如权利要求1所述的一种针对滑坡的InSAR形变监测结果解释方法,其特征在于,所述步骤S4包括:S41、根据所述数字高程数据,计算所述滑坡体的坡度值和坡向值;将所述坡向值根据预设规则转换为罗盘值;S42、根据所述形变分区结果,分别对所述坡度值和坡向值进行裁剪,得到不同形变分区内的坡度数据集和坡向数据集;S43、根据所述形变分区结果以及所述InSAR形变速率结果,分别计算每个所述形变分区内的InSAR形变速率平均值,并根据卫星成像角度和斜坡方向的几何关系,对所述InSAR形变速率平均值进行转换,得到所述滑坡体沿斜坡方向的形变速率;S44、根据每个所述形变分区的斜坡方向、斜坡方向形变速率量级、斜坡方向形变速率正负情况和各期次形变变化量,确定每个所述形变分区的滑动方向,以及各个所述形变分区之间的相对关系,得到所述滑坡体的空间演变规律。
7.如权利要求6所述的一种针对滑坡的InSAR形变监测结果解释方法,其特征在于,分别通过如下公式计算所述滑坡体的坡度值和坡向值:上式中,slope表示坡度值;aspect表示坡向值;f x 表示所述数字高程数据中的每个像元值在X方向的变化率;f y 表示所述数字高程数据中的每个像元值在Y方向的变化率;atan表示求反正切值;atan2表示求方位角。
8.如权利要求6所述的一种针对滑坡的InSAR形变监测结果解释方法,其特征在于,通过如下公式对所述InSAR形变速率平均值进行转换:上式中, V slope 表示所述滑坡体沿斜坡方向的形变速率; V los 表示每个所述形变分区内的InSAR形变速率平均值; c 表示斜坡方向和卫星视线方向夹角的余弦值。