一种区域浅表滑坡时空分异规律的ESDA分析法制造技术

本发明专利技术公开了一种区域浅表滑坡时空分异规律的ESDA分析方法。首先获取滑坡易发区多源时空数据，提取滑坡时空演化痕迹和变化环境下的滑坡孕灾因子信息，对因子时空分布进行ESDA分析；其次，计算时空综合易发性指数，选取最优的半变异函数模型进行指数空间估值与表面拟合；最后，基于ESDA法的易发性指数时空分异规律图谱分析。本发明专利技术解决了现有方法中易发性评估领域遇到的准确度不高、滑坡分异机制不明确对易发性影响顾及不够等瓶颈问题。方法充分考虑到了滑坡时空分异现象及规律的挖掘，在确定滑坡易发区与频发点更准确，方法的复杂程度更低，易于对多尺度空间尤其是大中尺度滑坡易发性问题展开分析。

A ESDA analysis method for temporal and spatial differentiation of surface shallow landslides

The invention discloses a ESDA analysis method for the spatio-temporal differentiation law of a shallow landslide in a region. First get the landslide prone area of multi-source spatiotemporal data extraction, landslide Landslide Disaster factor information evolution traces and changing environment, the ESDA factor analysis of temporal and spatial distribution; secondly, calculating the spatio-temporal susceptibility index, semi variogram model for selecting the optimal index space valuation and surface fitting; finally, ESDA method easy time of index variation based on spectrum analysis. The invention solves the bottleneck problems of the low accuracy, the unclear mechanism of the landslide differentiation and the insufficient influence on the influence of the occurrence in the field of the existing evaluation methods. Mining method considering the landslide spatial differentiation phenomenon and regularity, in determining landslide prone areas and frequent point method is more accurate, less complicated, easy to scale space especially large and medium scale landslide susceptibility analysis.

【技术实现步骤摘要】
一种区域浅表滑坡时空分异规律的ESDA分析法
本专利技术涉及滑坡易发性评估领域，特别涉及一种区域浅表滑坡时空分异规律的ESDA分析方法。
技术介绍
近年来，由于极端气候和复杂地形条件影响国内外滑坡等地质灾害频发，其中暴雨或短时特大暴雨等强对流天气极易诱发大规模的山体滑坡，危害性极大。由于信息熵、神经网络等空间数据驱动方法与数学模型往往忽视了滑坡空间结构差异的存在，未考虑计算图层内的空间数据自相关性与多图层数据的空间互相关性，因此空间数据叠置运算过程中易引入误差，给易发性评估精确性带来影响。针对我国西部地区滑坡分布广泛的国情，考虑到目前滑坡风险评估方法并不完善，准确率低。对区域滑坡的分布规律，尤其是时空分布规律并不明确的现状，本专利技术提出一种区域浅表滑坡时空分异规律的新分析手段。本专利技术通过一种区域浅表滑坡时空分异规律的ESDA分析方法(探索性空间数据分析法)，是在滑坡综合易发性指数和时空分布数据处理的基础上采用一种更精确的空间数据分析方法，该方法充分考虑到滑坡易发区的真实演化规律，基于时空发展趋势(历年、逐季)及易发性时空结构分析模型获得数据的时空集散规律。该方法的过程包括：首先获取滑坡易发区多源时空数据，提取滑坡时空演化痕迹和变化环境下的滑坡孕灾因子信息，基于探索性数据分析和变异函数对滑坡敏感性指数的空间分布特征进行检验，发现数据存在的空间变化趋势并找出数据内部的空间分异规律。其次，计算时空综合易发性指数，选取最优的半变异函数模型对易发性指数进行空间估值与表面拟合，对区域滑坡易发性进行更精细的描述和预测。最后利用区域滑坡编目数据对易发性区划进行改正和精度评价，基于ESDA法的易发性指数时空分异规律图谱分析实现区域浅表滑坡时空分异规律的发掘。
技术实现思路
专利技术的目的是提供一种区域浅表滑坡时空分异规律的ESDA分析方法，针对历史发生的滑坡分布数据特点进行自主性探索性空间分析，能够结合高时空分辨率光学影像提取滑坡演化信息与时空演化痕迹，提取关键变化环境下的孕灾因子信息，用于滑坡时空分异规律因素计算，为滑坡时空分异规律发掘和孕灾模式发现提供一种更准确的技术手段。本专利技术公开了一种区域浅表滑坡时空分异规律的ESDA分析方法。包括区域浅表滑坡分布指数的计算、空间变异性指数计算模型、空间估值拟合指数模型。方法的特征在于基于不断更新的历史滑坡数据，将滑坡面点状化，提取滑坡灾害点的易发性关键信息，构建采用ESDA法分析检验下的滑坡时空综合易发性指数模型，保障时空数据的有效利用与合理检验。具体步骤为：(1)获取滑坡易发区多源时空数据：多时相遥感数据、地形地质数据、人类工程活动数据，并进行几何精校正、影像增强、空间数据矢量化、影像裁剪处理。(2)时空数据处理与信息提取：采用GIS叠置分析工具，将解译的滑坡灾害信息与基础地理地质信息结合，提取灾害点的空间属性信息。(3)滑坡孕灾因子时空分布分析：应用探索性空间数据分析对历史滑坡点的各因子属性进行信息挖掘，通过距离变程分析孕灾因子的时间与空间变异性，探索滑坡孕灾因子的时空分布与不同因素对滑坡发育的内在影响规律。建立评价指标体系，利用面积频率比模型和层次分析法确定滑坡易发性评价指标权重。(4)计算时空综合易发性指数并制图，利用ESDA法分析滑坡时空综合易发性指数的分布规律、总体趋势、空间自相关及空间变异程度；通过距离变程分析、块金系数变化分析不同时期滑坡易发性指数的时空分异规律与变化趋势。(5)空间估值与表面拟合：采用最优表面拟合函数估值模型对滑坡综合易发性指数进行空间估值与表面拟合，空间插值得到空间分布状态图。(6)通过对比同期历史滑坡分布数据以及采用ROC曲线等方式评价精度，进行易发性区划改正。(7)时空分异规律图谱：针对不同时期滑坡易发性分布图谱，利用ESDA法分析其潜在的易发性及其时空分异演化趋势。本专利技术解决了现有方法中易发性评估领域遇到的准确度不高、滑坡分异机制不明确对易发性影响顾及不够等瓶颈问题。方法充分顾及到了滑坡时空分异现象及规律的挖掘，在确定滑坡易发区与频发点更准确，精度更高，方法的复杂程度更低，易于对多尺度空间尤其是大中尺度滑坡易发性问题进行分析。附图说明图1为本专利技术技术流程图。图2为本专利技术方法细节图。图3为本专利技术滑坡综合易发性指数QQplot分布图。图4为本专利技术滑坡综合易发指数空间半变异函数图。具体实施方式实施例:为了更清楚地说明本专利技术的技术方案，下面结合具体实施例对本专利技术做详细说明，以下实施例有助于本领域的技术人员更好的理解本专利技术。应当指出的是，本领域内的其他技术人员在不脱离于本专利技术和实例的基础上演变而来的其他实例，都属于本专利技术的保护范围。图1所示一种区域浅表滑坡时空分异规律的ESDA分析法方法的技术流程图，其详细与关键步骤如下：(1)时空数据处理与信息提取：典型浅表滑坡区近年滑坡识别与编目，面向高精度光学遥感数据的滑坡解译调查，滑坡灾害遥感监测数据来源，采用GoogleEarth、我国高分系列卫星可见光红外波段提取滑坡和植被信息，使用美国ASTERDEM作为滑坡调查与监测评估的背景环境评估数据；收集完善工区地形数据、地质数据、气象数据、人类工程活动数据，筛选有价值信息并从中挖掘可用信息。(2)滑坡发生后，坡体上的岩土体沿滑动面下滑，滑动过程中地表原有物质被掩埋和冲毁，通过对比灾前和灾后两个不同时期的遥感影像，识别因滑坡而改变的地面区域，划出滑坡破坏范围。(3)利用GIS软件对地形图、地质图以及遥感影像进行几何精校正，将滑坡面点状化。结合ASTER-GDEM2地形数据开展解译标志的建立，获取滑坡体的位置、面积、体积、类型、滑移的距离/方向等信息。结合滑坡点数据，从已知地灾系统电子地图采集和导出滑坡灾害点名称和坐标。(4)根据不同时序的典型滑坡遥感影像图，应用GIS软件对多时相遥感影像进行解译，获得土地覆被/利用数据；结合ASTER-GDEM2地形数据处理，获得工作区地形、地貌、植被、人类活动空间数据，完成孕灾环境遥感信息提取。应用全球ASTERDEM2数据集，从遥感和地形数据中提取分析水文因子和地貌因子数据；基于RS、GNSS、GIS技术集成的滑坡空间特征参数提取，反演获得地表浅表关键参数信息；遥感图像灾害环境参数解译、微地貌、与工程地质岩组含水量、NDVI、土地覆盖指数、河流侵蚀指数、治理工程指数、归一化水体指数(MNDWI)。(5)滑坡孕灾因子时空分布分析：对于不同尺度、不同时间分辨率的遥感影像数据、地图数据、文字资料和实地采集数据，通过最大值、最小值、平均值、中数、偏态和峰度等指标及参数直接反应样点数据的属性，基于ESDA直方图分析，统计分析各滑坡孕灾因子对滑坡发生的敏感性程度与频率分布规律。基于因子敏感性分析筛选相对重要的评价因子构建滑坡易发性评价指标体系。(6)基于历史因子分布图建立滑坡孕灾环境静态评价指标和动态评价指标，静态指标包括地层岩性、断层分布密度，动态指标包括坡度、高程、坡向、曲率、多年平均降雨量、归一化植被指数、土地利用、河流侵蚀、人类工程活动等主要致灾因子信息。(7)计算时空综合易发性指数并制图，对评估指标体系进行优化与权重确定，计算评估分值，结合面积频率比模型和层次分析法确定因子权重，具体步骤如下：1、构造判断矩阵，2、根据判本文档来自技高网...

【技术保护点】
一种区域浅表滑坡时空分异规律的ESDA分析法，其特征在于具体步骤为：（1）获取滑坡易发区多源时空数据：多时相遥感数据、地形地质数据、人类工程活动数据，并进行几何精校正、影像增强、空间数据矢量化、影像裁剪处理；（2）时空数据处理与信息提取：采用GIS叠置分析工具，将解译的滑坡灾害信息与基础地理地质信息结合，提取灾害点的空间属性信息；（3）滑坡孕灾因子时空分布分析：应用探索性空间数据分析对历史滑坡点的各因子属性进行信息挖掘，通过距离变程分析孕灾因子的时间与空间变异性，探索滑坡孕灾因子的时空分布与不同因素对滑坡发育的内在影响规律；建立评价指标体系，利用面积频率比模型和层次分析法确定滑坡易发性评价指标权重；（4）计算时空综合易发性指数并制图，利用ESDA法分析滑坡时空综合易发性指数的分布规律、总体趋势、空间自相关及空间变异程度；通过距离变程分析、块金系数变化分析不同时期滑坡易发性指数的时空分异规律与变化趋势；（5） 空间估值与表面拟合：采用最优表面拟合函数估值模型对滑坡综合易发性指数进行空间估值与表面拟合，空间插值得到空间分布状态图；（6）通过对比同期历史滑坡分布数据以及采用ROC曲线等方式评价精度，进行易发性区划改正；（7） 时空分异规律图谱：针对不同时期滑坡易发性分布图谱，利用ESDA法分析其潜在的易发性及其时空分异演化趋势。...

【技术特征摘要】
1.一种区域浅表滑坡时空分异规律的ESDA分析法，其特征在于具体步骤为：（1）获取滑坡易发区多源时空数据：多时相遥感数据、地形地质数据、人类工程活动数据，并进行几何精校正、影像增强、空间数据矢量化、影像裁剪处理；（2）时空数据处理与信息提取：采用GIS叠置分析工具，将解译的滑坡灾害信息与基础地理地质信息结合，提取灾害点的空间属性信息；（3）滑坡孕灾因子时空分布分析：应用探索性空间数据分析对历史滑坡点的各因子属性进行信息挖掘，通过距离变程分析孕灾因子的时间与空间变异性，探索滑坡孕灾因子的时空分布与不同因素对滑坡发育的内在影响规律；建立评价指标体系，利用面积频率比模型...

【专利技术属性】
技术研发人员：李百寿，高玉久，
申请(专利权)人：桂林理工大学，
类型：发明
国别省市：广西,45