本发明涉及一种滑坡灾害隐患识别方法，具体是指一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法。
背景技术：
合成孔径雷达干涉测量技术作为当前地表变形监测和区域尺度滑坡灾害监测与识别的有效技术，已经是滑坡灾害隐患识别必不可少的方法和手段，但目前国内对于该技术应用前、应用中和应用后的数据选择、方法流程及数据分析方法不甚清楚，大量相关工作人员利用该技术对于滑坡灾害隐患识别和分析效率太低。所以，一种提高滑坡灾害隐患识别的工作效率和经济效益的基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法体系成为人们亟待解决的问题。
技术实现要素：
本发明要解决的技术问题是现有技术对于合成孔径雷达干涉测量技术的应用不清楚、对于滑坡灾害隐患识别和分析效率太低等问题。
为解决上述技术问题，本发明提供的技术方案为：一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法，所述方法如下所示：
步骤1：滑坡灾害调查：主要调查内容包括区域地质构造、地层岩性、地形地貌、地表覆盖类型以及地质灾害分布，所述地质灾害分布与地质环境背景服务于数据处理分析中变形体的初步识别和验证，地形地貌服务于合成孔径雷达数据应用过程中可行性分析，地表覆盖类型用于最终可能变形信息分布的预判；
步骤2：合成孔径雷达数据查询：用于确定覆盖研究区、可获取、满足数据量需求的合成孔径雷达数据类型和相关几何参数；
步骤3:合成孔径雷达干涉测量可行性分析：是对区域可获取合成孔径雷达数据类型进行地形可视性分析，获取合成孔径雷达数据的地形可视区域和几何畸变区域，验证合成孔径雷达数据在山区变形监测的可行性；所述地形可行性分析是通过gis平台，利用地形坡度、坡向和卫星视线方向的几何特征结合山体阴影模型和栅格计算器，模拟合成孔径雷达数据地形可视性系数和叠掩及阴影所处位置，进一步结合滑坡数据库对研究区合成孔径雷达数据的滑坡可视性进行分析，评价合成孔径雷达数据的滑坡灾害监测能力，筛选适用于研究区、经济效益最高的合成孔径雷达数据类型；
步骤4：合成孔径雷达干涉测量技术优选：利用差分干涉测量、永久散射体干涉测量和小基线集干涉测量技术通过影像裁剪、基线估计、影像配准、干涉、去除地形相位、相位滤波、相位解缠、误差去除和地理编码获取wgs84坐标系下地表变形数据，通过区域分析、地面监测数据验证和野外考察验证合成孔径雷达干涉测量技术的优劣和变形监测结果的有效性，筛选适用于研究区的最优合成孔径雷达干涉测量技术；
步骤:5：滑坡灾害隐患早期识别：基于最优合成孔径雷达干涉测量技术获取的地表变形速率，结合历史滑坡遥感影像和地形地貌特征建立滑坡灾害隐患解译标志和变形阈值，对滑坡灾害隐患进行遥感解译和圈定，并制定野外考察路线，进行详细野外考察，验证和精炼滑坡灾害隐患识别结果；
步骤6：区域滑坡活动性和强度评价：基于长时间序列历史变形数据，对滑坡历史形变和当前形变进行比较分析，评价滑坡的活动性，利用gis平台空间分析功能提取滑坡最大变形速率，确定滑坡活动强度；
步骤7：滑坡变形模式分析：将滑坡变形剖面结合地质剖面确定滑坡空间变形模式，确定滑坡搬运区和堆积区；基于滑坡空间变形分布特征划分滑坡的二维空间变形分区，将滑坡不同位置的时间序列变形结合降雨、地震等数据分析滑坡时间变形模式；最后结合滑坡时间变形模式、空间变形模式和地形地貌特征分析滑坡运动机制。
进一步地，所述滑坡灾害调查的目的在于对研究区的地质环境背景和灾害现状进行深入的了解，认识地质灾害发育特征、孕灾背景、灾害分布规律等，尤其对区域典型灾害变形特征及其遥感影像光谱特征进行总结，服务于后期遥感数据处理过程中的人为参与和分析。
进一步地，所述合成孔径雷达数据查询的方式如下所示：通过欧空局软件、阿拉斯加卫星数据查询网站、哨兵数据查询网站和相关合成孔径雷达商业卫星代理商对研究区所有合成孔径雷达数据进行查询，获取研究区存档合成孔径雷达数据类型和覆盖范围。
进一步地，所述合成孔径雷达干涉测量可行性分析有针对性的选择下载、订购研究所需数据，减少数据处理量，提高数据处理和分析效率，提高合成孔径雷达数据利用率，降低购买合成孔径雷达数据的经济消耗。
进一步地，所述滑坡变形模式分析针对性的区分滑坡搬运和堆积区，确定滑坡运动部位，揭示滑坡加速过程，能够提高灾害体地质调查效率，提升灾害防治工程设计的准确性。
与现有技术相比，本发明的优点：本发明利用合成孔径雷达干涉测量技术在区域尺度开展滑坡灾害隐患早期识别的方法体系，该方法体系规范了在利用合成孔径雷达干涉测量技术识别区域潜在滑坡灾害隐患中的各个步骤和所采取的方法，能够提升滑坡灾害隐患识别的工作效率和经济效益，滑坡变形模式分析成果可有效提升滑坡灾害调查效率和防治工程设计准确性；使利用合成孔径雷达干涉测量技术开展滑坡灾害隐患识别的方法流程更加完善；本发明设计合理，值得大力推广。
附图说明
图1是一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法的流程图；
如图所示：sar数据为合成孔径雷达数据，insar为合成孔径雷达干涉测量，dinsar为差分干涉测量，psinsar为永久散射体干涉测量，sbas-insar为小基线集干涉测量。
具体实施方式
下面结合附图对本发明做进一步的详细说明。
结合附图1，对本发明进行详细介绍。
本发明在具体实施时提供了一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法，
所述方法如下所示：
步骤1：滑坡灾害调查：主要调查内容包括区域地质构造、地层岩性、地形地貌、地表覆盖类型以及地质灾害分布，所述地质灾害分布与地质环境背景服务于数据处理分析中变形体的初步识别和验证，地形地貌服务于合成孔径雷达数据应用过程中可行性分析，地表覆盖类型用于最终可能变形信息分布的预判；
步骤2：合成孔径雷达数据查询：用于确定覆盖研究区、可获取、满足数据量需求的合成孔径雷达数据类型和相关几何参数；
步骤3:合成孔径雷达干涉测量可行性分析：是对区域可获取合成孔径雷达数据类型进行地形可视性分析，获取合成孔径雷达数据的地形可视区域和几何畸变区域，验证合成孔径雷达数据在山区变形监测的可行性；所述地形可行性分析是通过gis平台，利用地形坡度、坡向和卫星视线方向的几何特征结合山体阴影模型和栅格计算器，模拟合成孔径雷达数据地形可视性系数和叠掩及阴影所处位置，进一步结合滑坡数据库对研究区合成孔径雷达数据的滑坡可视性进行分析，评价合成孔径雷达数据的滑坡灾害监测能力，筛选适用于研究区、经济效益最高的合成孔径雷达数据类型；
步骤4：合成孔径雷达干涉测量技术优选：利用差分干涉测量、永久散射体干涉测量和小基线集干涉测量技术通过影像裁剪、基线估计、影像配准、干涉、去除地形相位、相位滤波、相位解缠、误差去除和地理编码获取wgs84坐标系下地表变形数据，通过区域分析、地面监测数据验证和野外考察验证合成孔径雷达干涉测量技术的优劣和变形监测结果的有效性，筛选适用于研究区的最优合成孔径雷达干涉测量技术；
步骤:5：滑坡灾害隐患早期识别：基于最优合成孔径雷达干涉测量技术获取的地表变形速率，结合历史滑坡遥感影像和地形地貌特征建立滑坡灾害隐患解译标志和变形阈值，对滑坡灾害隐患进行遥感解译和圈定，并制定野外考察路线，进行详细野外考察，验证和精炼滑坡灾害隐患识别结果；
步骤6：区域滑坡活动性和强度评价：基于长时间序列历史变形数据，对滑坡历史形变和当前形变进行比较分析，评价滑坡的活动性，利用gis平台空间分析功能提取滑坡最大变形速率，确定滑坡活动强度；
步骤7：滑坡变形模式分析：将滑坡变形剖面结合地质剖面确定滑坡空间变形模式，确定滑坡搬运区和堆积区；基于滑坡空间变形分布特征划分滑坡的二维空间变形分区，将滑坡不同位置的时间序列变形结合降雨、地震等数据分析滑坡时间变形模式；最后结合滑坡时间变形模式、空间变形模式和地形地貌特征分析滑坡运动机制。
所述滑坡灾害调查的目的在于对研究区的地质环境背景和灾害现状进行深入的了解，认识地质灾害发育特征、孕灾背景、灾害分布规律等，尤其对区域典型灾害变形特征及其遥感影像光谱特征进行总结，服务于后期遥感数据处理过程中的人为参与和分析。
所述合成孔径雷达数据查询的方式如下所示：通过欧空局软件、阿拉斯加卫星数据查询网站、哨兵数据查询网站和相关合成孔径雷达商业卫星代理商对研究区所有合成孔径雷达数据进行查询，获取研究区存档合成孔径雷达数据类型和覆盖范围。
所述合成孔径雷达干涉测量可行性分析有针对性的选择下载、订购研究所需数据，减少数据处理量，提高数据处理和分析效率，提高合成孔径雷达数据利用率，降低购买合成孔径雷达数据的经济消耗。
所述滑坡变形模式分析针对性的区分滑坡搬运和堆积区，确定滑坡运动部位，揭示滑坡加速过程，能够提高灾害体地质调查效率，提升灾害防治工程设计的准确性。
本发明利用合成孔径雷达干涉测量技术在区域尺度开展滑坡灾害隐患早期识别的方法体系，该方法体系规范了在利用合成孔径雷达干涉测量技术识别区域潜在滑坡灾害隐患中的各个步骤和所采取的方法，能够提升滑坡灾害隐患识别的工作效率和经济效益，滑坡变形模式分析成果可有效提升滑坡灾害调查效率和防治工程设计准确性；使利用合成孔径雷达干涉测量技术开展滑坡灾害隐患识别的方法流程更加完善；本发明设计合理，值得大力推广。
以上对本发明及其实施方式进行了描述，这种描述没有限制性，实际的结构并不局限于此。总而言之如果本领域的普通技术人员受其启示，在不脱离本发明创造宗旨的情况下，不经创造性的设计出与该技术方案相似的结构方式及实施例，均应属于本发明的保护范围。
1.一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法，其特征在于：所述方法如下所示：
步骤1：滑坡灾害调查：主要调查内容包括区域地质构造、地层岩性、地形地貌、地表覆盖类型以及地质灾害分布，所述地质灾害分布与地质环境背景服务于数据处理分析中变形体的初步识别和验证，地形地貌服务于合成孔径雷达数据应用过程中可行性分析，地表覆盖类型用于最终可能变形信息分布的预判；
步骤2：合成孔径雷达数据查询：用于确定覆盖研究区、可获取、满足数据量需求的合成孔径雷达数据类型和相关几何参数；
步骤3:合成孔径雷达干涉测量可行性分析：是对区域可获取合成孔径雷达数据类型进行地形可视性分析，获取合成孔径雷达数据的地形可视区域和几何畸变区域，验证合成孔径雷达数据在山区变形监测的可行性；所述地形可行性分析是通过gis平台，利用地形坡度、坡向和卫星视线方向的几何特征结合山体阴影模型和栅格计算器，模拟合成孔径雷达数据地形可视性系数和叠掩及阴影所处位置，进一步结合滑坡数据库对研究区合成孔径雷达数据的滑坡可视性进行分析，评价合成孔径雷达数据的滑坡灾害监测能力，筛选适用于研究区、经济效益最高的合成孔径雷达数据类型；
步骤4：合成孔径雷达干涉测量技术优选：利用差分干涉测量、永久散射体干涉测量和小基线集干涉测量技术通过影像裁剪、基线估计、影像配准、干涉、去除地形相位、相位滤波、相位解缠、误差去除和地理编码获取wgs84坐标系下地表变形数据，通过区域分析、地面监测数据验证和野外考察验证合成孔径雷达干涉测量技术的优劣和变形监测结果的有效性，筛选适用于研究区的最优合成孔径雷达干涉测量技术；
步骤:5：滑坡灾害隐患早期识别：基于最优合成孔径雷达干涉测量技术获取的地表变形速率，结合历史滑坡遥感影像和地形地貌特征建立滑坡灾害隐患解译标志和变形阈值，对滑坡灾害隐患进行遥感解译和圈定，并制定野外考察路线，进行详细野外考察，验证和精炼滑坡灾害隐患识别结果；
步骤6：区域滑坡活动性和强度评价：基于长时间序列历史变形数据，对滑坡历史形变和当前形变进行比较分析，评价滑坡的活动性，利用gis平台空间分析功能提取滑坡最大变形速率，确定滑坡活动强度；
步骤7：滑坡变形模式分析：将滑坡变形剖面结合地质剖面确定滑坡空间变形模式，确定滑坡搬运区和堆积区；基于滑坡空间变形分布特征划分滑坡的二维空间变形分区，将滑坡不同位置的时间序列变形结合降雨、地震等数据分析滑坡时间变形模式；最后结合滑坡时间变形模式、空间变形模式和地形地貌特征分析滑坡运动机制。
2.根据权利要求1所述的一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法，其特征在于：所述滑坡灾害调查的目的在于对研究区的地质环境背景和灾害现状进行深入的了解，认识地质灾害发育特征、孕灾背景、灾害分布规律等，尤其对区域典型灾害变形特征及其遥感影像光谱特征进行总结，服务于后期遥感数据处理过程中的人为参与和分析。
3.根据权利要求1所述的一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法，其特征在于：所述合成孔径雷达数据查询的方式如下所示：通过欧空局软件、阿拉斯加卫星数据查询网站、哨兵数据查询网站和相关合成孔径雷达商业卫星代理商对研究区所有合成孔径雷达数据进行查询，获取研究区存档合成孔径雷达数据类型和覆盖范围。
4.根据权利要求1所述的一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法，其特征在于：所述合成孔径雷达干涉测量可行性分析有针对性的选择下载、订购研究所需数据，减少数据处理量，提高数据处理和分析效率，提高合成孔径雷达数据利用率，降低购买合成孔径雷达数据的经济消耗。
5.根据权利要求1所述的一种基于合成孔径雷达干涉测量的滑坡灾害隐患识别方法，其特征在于：所述滑坡变形模式分析针对性的区分滑坡搬运和堆积区，确定滑坡运动部位，揭示滑坡加速过程，能够提高灾害体地质调查效率，提升灾害防治工程设计的准确性。
技术总结