CN113642849A - 考虑空间分布特征的地质灾害危险性综合评价方法及装置 - Google Patents
考虑空间分布特征的地质灾害危险性综合评价方法及装置 Download PDFInfo
- Publication number
- CN113642849A CN113642849A CN202110814485.8A CN202110814485A CN113642849A CN 113642849 A CN113642849 A CN 113642849A CN 202110814485 A CN202110814485 A CN 202110814485A CN 113642849 A CN113642849 A CN 113642849A
- Authority
- CN
- China
- Prior art keywords
- geological disaster
- evaluation
- risk
- disaster
- historical
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 136
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000010801 machine learning Methods 0.000 claims abstract description 34
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 31
- 238000007637 random forest analysis Methods 0.000 claims abstract description 24
- 230000000694 effects Effects 0.000 claims abstract description 23
- 238000007781 pre-processing Methods 0.000 claims abstract description 19
- 238000012706 support-vector machine Methods 0.000 claims abstract description 15
- 238000007477 logistic regression Methods 0.000 claims abstract description 12
- 239000006185 dispersion Substances 0.000 claims abstract description 6
- 238000012549 training Methods 0.000 claims description 46
- 238000011160 research Methods 0.000 claims description 35
- 238000012360 testing method Methods 0.000 claims description 33
- 238000000605 extraction Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 239000013598 vector Substances 0.000 claims description 6
- 230000001788 irregular Effects 0.000 claims description 5
- 230000002776 aggregation Effects 0.000 claims description 4
- 238000004220 aggregation Methods 0.000 claims description 4
- 230000001131 transforming effect Effects 0.000 claims description 2
- 230000002265 prevention Effects 0.000 abstract description 9
- 238000013439 planning Methods 0.000 abstract description 3
- 238000010276 construction Methods 0.000 description 7
- 238000003066 decision tree Methods 0.000 description 4
- 230000008030 elimination Effects 0.000 description 4
- 238000003379 elimination reaction Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000013179 statistical model Methods 0.000 description 4
- 238000007689 inspection Methods 0.000 description 3
- 238000011835 investigation Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000012502 risk assessment Methods 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/211—Selection of the most significant subset of features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Bioinformatics & Computational Biology (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Economics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Software Systems (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑空间分布特征的地质灾害危险性综合评价方法。针对历史地质灾害点的空间集聚与分散特征,分别提出利用聚类算法提取区域聚类属性作为评价指标以及基于渔网网格的数据预处理方法,基于逻辑回归(LR)、支持向量机(SVM)、梯度提升树(GBDT)与随机森林(RF)的多机器学习算法构建模型,通过对比模型预测精度与地图评价效果,确定最优模型并绘制危险性地图,同时给出检验技术可靠性的实验方案。本发明方法能显著提高模型评估精度、增强模型评价性能,生成预测更准、质量更优的地质灾害危险性地图,从而为灾害风险防治规划工作提供契合实情的决策依据。
Description
技术领域
本发明属于地质灾害工程技术领域,具体涉及一种考虑空间分布特征的地质灾害危险性综合评价方法 及装置。
背景技术
近年来,全球气候变化加剧与极端天气事件频发背景下,我国的地质灾害呈现多发趋势,局部地区可 能加重。地质灾害突发性、群发性、复杂性、不确定性的特点给灾前勘查防治与灾后应急处置工作带来诸 多困难,对人民生命财产、区域经济发展、生态环境安全均产生严重威胁。因此,详实精细的地质灾害危 险性评价结果尤为重要,为揭示地质灾害危险性区划特征、开展区域防治规划建设提供了可靠依据。同时 使地质灾害防治对策的制定更具针对性,对合理调配防灾减灾资源、切实减轻地质灾害危害具有重要意义。
评估主体方面,当前地质灾害危险性评价多局限在崩塌或滑坡等单一灾种,缺乏对崩塌滑坡泥石流为 代表的多灾种地质灾害危险性的综合评价,考虑到实际多种地质灾害普遍共存的情况,单灾种成果的全面 性、适用性有所欠缺。然而,由于多灾种地质灾害的时空构成与成因机制更为复杂,现有综合地质灾害危 险性评估的效果与精度较之单灾种偏低,仍需提出新方法予以改进。
评价技术层面,可分为基于物理或统计的方法。物理模型通常基于力学分析评估斜坡稳定性,蕴含灾 害发生机理,能提供灾害强度信息且预测精度高,但需要精细化的数据支持,仅适用于局部小尺度的区域 研究。统计模型多基于数据驱动,预测发生概率、提供灾害空间信息,常用的定量评价方法包括信息量模 型、层次分析法、频率比法、证据权法、确定性系数法及多种方法的耦合分析等,可信度高,应用范围广。 但这类传统统计模型对地质灾害与影响要素间存在的复杂非线性关系的模拟能力较弱,机器学习模型能有 效解决该问题,同时克服了指标选取、权重确定过程中主观性强的缺点,已逐步成为传统统计模型的替代 方法。常用的机器学习模型包括逻辑回归(LR)、支持向量机(SVM)、随机森林(RF)、人工神经网 络(ANN)等。不同研究区自然环境与成灾条件有别,导致各算法在不同研究区的适用性存在差异。因此, 有必要依据研究区实际灾情因地制宜地选择最为合适的机器学习模型。
已有技术的评价因子一般包含地形地质、气象水文、土地覆盖、人类活动等方面,但忽视了地质灾害 的空间分布特征。不少研究均表明,不论是活动断裂抑或强降雨诱发的地质灾害均具有明显的群发性特征, 在空间分布上既有分散性又有聚集性。这种聚集效应表现为地质灾害的多发地段或区域存在灾害群,反映 了区域灾害的易发性程度本质上是由区域性地质构造以及降水、地震、地形地貌等多条件共同制约形成, 因此灾害点的集聚性隐含了潜在的地质灾害危险性信息。然而现有技术并没有将灾害的集聚性特征作为评价指标构建模型,一定角度而言损失了可用信息。而地质灾害的分散性表现为发生地点分散、覆盖面广, 呈现“星罗棋布”的空间格局,体现了地质灾害的随机性与偶发性,这种随机性降低了模型训练精度,对 危险性评估带来了“噪声”干扰,进而使评估结果与实际灾情存在显著偏差,但是鲜有评估方案对此进行 分析、处理。
发明内容
针对传统多灾种地质灾害危险性评价精度不高,模型构建过程中忽视空间集聚性与分散性特征对评估 有效性的潜在影响,以及如何挖掘地质灾害空间分布所内含的同区块间均质性与不同区块间异质性信息等 问题,本发明提供一种考虑空间分布特征的地质灾害危险性综合评价方法及装置。
本发明的目的在于提供一种考虑空间分布特征的地质灾害危险性综合评价方法及装置。
一种考虑空间分布特征的地质灾害危险性评价方法,包括:
步骤A:利用聚类算法依据历史地质灾害情况数据提取地质灾害聚类属性作为区域地质灾害危险性评 价的评价因子;
步骤B:对历史地质灾害情况数据进行预处理,基于渔网格网过滤历史地质灾害噪声点,生成区域地 质灾害危险性评价训练集及测试集;
步骤C:选取地质灾害危险性评价的其他评价因子,采用不规则研究区的空间属性叠加手段提取所述 其他评价因子;
步骤D:构建多机器学习模型,并依据步骤A中确定的所述地质灾害聚类属性和步骤C中确定的所述 其他评价因子、以及步骤B中获得的所述训练集训练所述模型;输出地质灾害危险性地图;
以及,确定最优地质灾害危险性地图,并依据所述最优地质灾害危险性地图对地质灾害危险性进行评价。
进一步地,步骤A中,所述地质灾害聚类属性的提取流程为:
①针对历史地质灾害情况数据输入搜索距离、最小聚类要素个数,经OPTICS算法计算得到灾害点 的聚类属性;
②根据每个历史地质灾害点的位置构建泰森多边形,并将各所述灾害点的聚类属性赋值给相应的泰 森多边形;
③将泰森多边形矢量图层转栅格,得到聚类属性评价因子,与步骤C中获得的所述其他评价因子一 同输入步骤D中的所述模型训练;
④整理不同搜索距离与最小聚类要素个数参数组合下各所述模型在所述训练集上的训练精度,对比 得出最佳方案,即最优搜索距离与最优最小聚类要素个数,并以此得到所述地质灾害聚类属性。
进一步地,步骤C中,还包括基于因子共线性分析方法对所述其他评价因子的选取进行合理性评价。
进一步地,步骤D中,基于逻辑回归(LR)、支持向量机(SVM)、梯度提升树(GBDT)与随机森 林(RF)的多机器学习算法构建所述模型。
进一步地,步骤D中,还包括基于步骤B中获得的所述测试集对各所述模型的预测精度进行评价。其 中,预测精度评价指标包括:精确率(Precision)、召回率(Recall)、正确率(Accuracy)和F1分值(F1 score),ROC曲线与AUC面积。
进一步地,还包括步骤E:对所述地质灾害危险性地图的评价效果识别与评估质量检验。从全局尺度 上不同危险区的空间分布格局与面积占比,以及局部尺度上空间细节刻画水平、渐变性与斑块化特征角度, 对比评价不同所述模型的评价效果,结合历史灾害点落入各层级危险区的比例即危险区对历史灾害的解释 率进一步评估所述地图的质量。依据对所述地质灾害危险性地图的评价效果识别与评估质量检验的结果, 确定所述最优地质灾害危险性地图。
进一步地,还包括步骤F:地质灾害空间分布特征考量的可靠性与实用性检验,包括厘定空间聚集特 征对所述模型性能的影响、厘定空间分散特征对所述模型性能的影响。
进一步地,对地质灾害危险性进行评价时包括依据发生概率的数值分布与不同危险区面积分布分析各 级水平的地质灾害危险性。
一种考虑空间分布特征的地质灾害危险性评价装置,包括:
历史地质灾害数据模块110,用于存储区域历史地质灾害数据;
地质灾害聚类属性提取模块120,用于利用聚类算法依据所述历史地质灾害数据模块110中存储的历 史地质灾害情况提取地质灾害聚类属性因子;
数据预处理模块130,用于对所述历史地质灾害数据模块110中存储的历史灾害情况数据进行预处理, 基于渔网格网过滤历史地质灾害噪声点,并生成区域地质灾害危险性评价训练集及测试集;
评价因子提取模块140,用于提取选取的地质灾害危险性评价的其他评价因子;
机器学习模块150,用于多机器学习模型的构建、训练,依据所述地质灾害聚类属性因子和所述其他 评价因子、所述训练集训练所述模型,并利用所述模型输出地质灾害危险性地图;
地质灾害危险性评价模块160,用于依据所述地质灾害危险性地图进行地质灾害危险性水平分析。
相对于现有技术,本发明的有益效果:
考虑了空间分布特征在地质灾害危险性评价问题中的潜在影响,分别揭示了空间集聚性特征与分散性 特征对提高预测精度、增强评价性能具有显著的正向作用与负面影响,对此提出了利用聚类算法提取区域 地质灾害聚类属性作为评价指标以及基于渔网网格的数据预处理方法以进一步改进现有技术。该方法能有 效提升地质灾害危险性地图的预测精度与评价质量,使高低危险区识别更精准、对历史灾害的可解释性更 强,体现出地质灾害危险性的空间群发性、区间异质性及区内渐变性特征,利于防灾资源的管控与调配,能为政府部门的灾害风险普查及防灾减灾工作提供更契合灾害实情的参考。
附图说明
为了更完整地理解本发明,现在将参考结合附图的以下描述,其中:
图1为根据本发明实施例的地质灾害危险性评价方法的流程图;
图2中:a)为本发明实施例研究区的OPTICS算法的聚类结果;b)为本发明实施例研究区的基于点簇 聚类结果构建的泰森多边形;
图3中:a)为本发明实施例研究区不同最小聚类要素个数下各算法的训练精度(搜索距离固定为5km); b)为本发明实施例研究区不同搜索距离下各算法的训练精度(最小聚类个数固定为6个);
图4中:a)为本发明实施例研究区的渔网格网过滤历史灾害点源后的结果;b)为本发明实施例研究区 的基于渔网格网生成的训练集;c)为本发明实施例研究区的基于渔网格网生成的测试集;
图5中:a)为本发明实施例研究区的各模型的预测精度指标结果;b)为本发明实施例研究区的各模型 的ROC曲线及AUC面积;
图6为本发明实施例研究区四种机器学习模型绘制的地质灾害危险性评价图;
图7中:a)为本发明实施例研究区的各模型不同危险区的面积占比;b)为本发明实施例研究区的各模 型不同危险区的灾害点数占比;
图8中:a)为本发明实施例研究区的随机森林模型各指标的重要度;b)为本发明实施例研究区的剔除 聚类属性指标后各模型的ROC曲线及AUC面积;c)为本发明实施例研究区的剔除聚类属性指标后的随机 森林地质灾害危险性地图;d)为本发明实施例研究区的剔除聚类属性指标前后模型预测效果对比;
图9中:a)为本发明实施例研究区的不采用渔网进行数据预处理时的训练集与测试集分布;b)为本发 明实施例研究区的不采用渔网进行数据预处理时各模型的ROC曲线及AUC面积;c)为本发明实施例研究 区的不采用渔网进行数据预处理时的随机森林地质灾害危险性地图;d)为本发明实施例研究区的采用渔网 进行数据预处理时前后模型预测效果对比;
图10中:a)为本发明实施例研究区的各县所有栅格单元地质灾害发生概率的数据分布情况;b)为本发 明实施例研究区的各县不同危险性区域的栅格数量分布情况。
具体实施方式
为说明清楚本发明的目的、技术细节及有效应用,使之便于本领域普通技术人员理解与实施,下面将 结合本发明实施例及附图作进一步的详细阐述。显然,此处描述的实施示例仅用于说明和解释本发明,并 不用于限定本发明。基于本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有 其他实施例,都属于本发明保护的范围。
本发明提供一种考虑空间分布特征的地质灾害危险性评价方法,下面以贵州省毕节市的地质灾害危险 性评价为例,说明本发明的技术内容。图1示意性地示出本发明地质灾害危险性评价方法的流程图,本发 明所述的地质灾害危险性评价方法主要包括以下步骤:
步骤A、利用聚类算法依据历史地质灾害情况提取地质灾害聚类属性作为区域地质灾害危险性评价的 重要评价指标因子。
具体而言,采用OPTICS聚类算法对毕节市历史灾害点进行聚类,该聚类算法原理及实现可参见 K.P.Agrawal等发表的文章“基于OPTICS的时空聚类技术的开发与验证”(Information Sciences,第369 卷,2016年11月10日,第388-401页;https://doi.org/10.1016/j.ins.2016.06.048)。
该算法有两个重要的调节参数:搜索距离与最小聚类要素个数,不同参数组合决定了不同尺度的聚类 结果。因此采用控制变量法对搜索距离与最小聚类要素个数进行测试,通过比较不同参数组合下各机器学 习算法的模型训练精度来确定最优参数组合。主要流程为:
①针对历史地质灾害情况输入搜索距离、最小聚类要素个数,经OPTICS算法计算得到灾害点的聚 类属性(参见图2.a);
②根据每个历史地质灾害点的位置构建泰森多边形,并将各所述灾害点的聚类属性赋值给相应的泰 森多边形(参见图2.b);
③将泰森多边形矢量图层转栅格,得到聚类属性评价因子,与步骤C中获得的其他14个评价因子一 同输入步骤D中的机器学习模型训练;
④整理不同搜索距离与最小聚类要素个数参数组合下各机器学习模型在训练集上的训练精度,对比 得出最佳方案(参见图3),即最优搜索距离与最优最小聚类要素个数。
图3说明本实施例的最佳参数组合为:搜索距离为5km,最小聚类要素个数为6。以该参数组合下的 结果即图2.b作为毕节市地质灾害的聚类属性评价因子。
步骤B:对历史地质灾害情况数据进行预处理,基于渔网格网过滤历史地质灾害噪声点,生成区域地 质灾害危险性评价训练集及测试集,其中,所述训练集及所述测试集均包括多个灾害样本点及多个非灾害 样本点。
首先,针对历史地质灾害情况以步骤A中确定的所述最优搜索距离长度为邻域边界(本实施例即为 5km),生成渔网切割研究区,当该渔网网格内的历史灾害点数不少于2个时保留,否则剔除该网格内的 灾害点(参见图4.a)。过滤后,本实施例中有288个渔网格网存在多个灾害点,共计灾害点1003个,随 机划分70%(702个)的灾害点作为训练集的灾害样本点,30%(301个)作为测试集的灾害样本点。
其次,为了保证样本的平衡性,需生成等量非灾害样本点用于训练与测试。基于上述渔网过滤方法, 以包含多灾害点的288个渔网格网为缓冲区,非灾害样本点随机生成在该缓冲区之外(参见图4.b及图4.c)。
步骤C:选取地质灾害危险性评价的其他评价因子,采用不规则研究区的空间属性叠加手段提取评价 因子,基于因子共线性分析方法对评价因子的选取进行合理性评价。
评价因子选取包括高程、坡度、坡向、平面曲率、剖面曲率的地形要素,距断层距离、岩性、土壤类 型的地质要素,年均降雨量、汇流累积量、距河流距离的水文要素,土地利用、归一化植被指数的土地覆 盖要素,距道路距离的人类活动要素以及聚类属性的历史灾害要素。
评价因子提取过程:以本实施例研究区的最小外接矩形提取各指标栅格的ASCII值(空值为-9999), 转化为列向量,再提取所有空值栅格的行列位置,若各评价因子的空值栅格位置完全一致,则等价于不规 则研究区的边界及内部栅格空间位置完全吻合,可以提取灾害点对应各评价因子的属性值。
评价因子选取合理性评价:基于Pearson相关系数计算相关系数矩阵以判断评价因子间是否存在成对 共线性,基于方差膨胀因子(VIF)与容忍度(TOL)值判断评价因子间是否存在多重共线性。Pearson相 关系数要求两变量为连续变量,故剔除类别变量分析共线性。同时根据容忍度(TOL)计算VIF。本实施 例各评价因子间的相关系数均小于0.7,各评价因子的TOL在0.356~0.994间,VIF在1.006~2.807间,表 明评价因子间不存在共线性,评价因子选取合理、可行。
步骤D:基于不同机器学习算法的原理构建多机器学习模型,并依据步骤A中确定的所述聚类属性评 价因子、步骤C中确定的所述其他评价因子数据集、以及步骤B中获得的所述训练集训练模型;基于步骤 B中获得的测试集对各所述多机器学习模型的预测精度进行评价;输出地质灾害危险性地图。
其中,包括:
步骤D1:基于逻辑回归(LR)、支持向量机(SVM)、梯度提升树(GBDT)与随机森林(RF)等 多机器学习算法构建多机器学习模型,输入步骤B中获得的所述训练集与依据步骤A中确定的所述聚类属 性评价因子、步骤C中确定的其他评价因子数据集,训练基于不同多机器学习算法的所述多机器学习模型。
其中,所述基于不同多机器学习算法的多机器学习模型的构建具体包括:
①逻辑回归(LR)的函数表达式为:
式中:P为地质灾害发生概率,P∈[0,1];α为常数项;β为反映该因子对P影响程度大小的偏回归系数;i 为评价因子种类数量。
P/(1-P)为地质灾害发生与不发生的概率比值,对其取自然对数得到:
ln[P/(1-P)]=α+β1x1+β2x2+…+βixi (2)
②支持向量机(SVM)的目标是寻求一个n维特征空间且法向量为w的超平面,即wTx+b=0,将2 类数据点正确分隔且使分类间隔最大,其构建过程为:给定训练样本集D={(xi,yi)},yi∈{-1,1},间隔 2/||w||最大化等价于最小化||w||2/2,由此得到SVM的基本型:
min(||w||2/2),s.t.yi(wTxi+b)≥1(i=1,2,…,m) (3)
式(3)是一个凸二次规化问题,约束条件为线性约束,可通过引入拉格朗日乘子得到其对偶问题, 这里不再赘述,最终可求得最优分类超平面:
式中:αi为拉格朗日乘子分量,xi为各评价因子构成的特征向量,K(xi,x)为核函数,b为偏置值。
SVM输出的是一个无标定量,其取值范围不定。为了便于对比不同算法的预测效果,需要将输出值 映射到[0,1]间,方法为:
p(x)=1/(1+eAf(x)+B) (5)
式中:p(x)为地质灾害危险性值,p(x)∈[0,1];A,B为待定系数,可依据最大似然法估定。
③梯度提升树(GBDT)是以CART回归树为基分类器的一种迭代决策树算法,其构建过程为:
式中:f(x)为最终的强分类器;T为最大迭代次数;ctj为在叶节点j上的损失函数负梯度方向的最佳拟合值; Rtj(j=1,2,…,J)为第t个回归树对应的叶节点区域;J为回归树t的叶节点数。输出值同样按照式(5)映射 到[0,1]。
④随机森林(RF)是以决策树为基本单元的一种集成学习算法。其构建过程为:首先利用bootstrap 抽样从训练集D中抽取k个样本(含m个特征),每个样本中随机选取n(n≤m)个特征构成特征空间,对 k个样本分别建立k个决策树模型,选择最优特征尽可能生长,得到k种分类结果{h1(X),h2(X),…,hk(X)}, 对每种结果进行投票表决或取均值作为最终分类结果。
步骤D2:基于步骤B中获得的所述测试集评价各所述多机器学习模型的预测精度;预测精度评价指 标包括精确率(Precision)、召回率(Recall)、正确率(Accuracy)和F1分值(F1 score),ROC曲线与 AUC面积。它们的定义如下:
Precision=TP/(TP+FP) (7)
Recall=TP/(TP+FN) (8)
Accuracy=(TP+TN)/(TP+FP+TN+FN) (9)
F1 socre=2*Precision*Recall/(Precision+Recall) (10)
式中:TP(True Positive,真阳性)是正确分类的地质灾害样本的数量;TN(TrueNegative,真阴性)是 正确分类的非地质灾害样本的数量;FP(False Positive,假阳性)是误分类的非地质灾害样本的数量;FN (False Negative,假阴性)是误分类的地质灾害样本的数量。
依据实施例毕节市14个地质灾害评价因子及聚类属性因子训练模型,基于测试集评价模型性能,模 型均在Matlab2019b中编写运行。表1列出各模型测试结果的混淆矩阵。
表1各模型测试集上的混淆矩阵
各项统计指标(参见图5.a)表明RF模型在训练集及测试集上表现均优于其他模型,其正确率达76%, 精确率为78%。RF与GBDT的F1分值明显更高,说明以决策树为基分类器的集成式算法性能明显优于 LR、SVM等传统机器学习算法。四种模型的线下面积AUC值均达到0.78以上,整体预测精度较高,其 中RF最高达0.825(参见图5.b)。结合表1,GBDT模型在灾害点识别效果方面与随机森林相近,但其 在非灾害点上预测性能欠。
步骤D3:基于上述模型预测研究区内所有栅格单元的地质灾害危险性概率值(p),按照不同危险性 层级水平划分评价阈值:极低危险区(0≤p≤0.1)、低危险区(0.1<p≤0.3)、中危险区(0.3<p≤0.6)、 高危险区(0.6<p≤0.9)、极高危险区(0.9<p≤1),输出地质灾害危险性地图。
具体而言,按照上述方式生成毕节市地质灾害危险性评价结果,如图6所示。
步骤E:对所述地质灾害危险性地图的评价效果识别与评估质量检验。从全局尺度上不同危险区的空 间分布格局与面积占比,以及局部尺度上空间细节刻画水平、渐变性与斑块化特征角度,对比评价不同模 型的评价效果,结合历史灾害点落入各层级危险区的比例即危险区对历史灾害的解释率(也即评估精度), 进一步评估地图质量。
从全局尺度、局部尺度、评估精度三个角度论证所述地质灾害危险性地图的评价效果与评估质量。结 合图6的地质灾害危险性地图与图7各模型不同危险区的统计情况,具体方法为:
①全局尺度上,RF地图中不同层级危险区的空间分布与面积占比更合理。例如极高危险区在高值区 的各县中均有适量分布,而LR和SVM地图中,极高危险区在毕节市最北侧、大方县西南侧、黔西县中 部东北侧集中成块分布,呈现局部夸大的特点,与真实灾情有偏,尤其是大方县西南侧过于高估,实际灾 害点仅有6个。GBDT地图的中危险区面积偏少,仅占17%(参见图6.a);
②局部尺度上,RF空间细节刻画更为丰富,在高值区域保有危险性从极高到高再到中的渐变性特点, 贴切实际,相反其他三个模型斑块化特征更明显,高、低危险区内的描述较为粗糙;
③评估精度上,RF模型中面积占比25%的高危区(包含高危险区和极高危险区)囊括了75%的历史 灾害点,同时面积占比47%的低危区(包含低危险区和极低危险区)仅包含了12%的灾害点(参见图7), 由此说明该危险性地图能客观有效地反映真实灾情。其他模型的表现精度欠佳,如GBDT高危区面积占比 更大(31%),但仅解释了67%的历史灾害,而SVM的低危区(63%)过多地解释了30%的历史灾害。
步骤F:地质灾害空间分布特征考量的可靠性与实用性检验。
其中包括:
步骤F1:厘定空间聚集特征对模型性能的影响。从上述机器学习模型的评价指标因子数据集中剔除聚 类属性因子并重新建模,通过对比剔除前后的精度评价指标、危险性地图的空间分布与评价效果来衡量空 间聚集特征的实际效用。
为了定量分析空间聚集特征对模型预测性能的影响,从原来的四种模型评价因子中剔除聚类属性因子 并重新建模。以RF模型为例,剔除聚类属性因子后,其训练精度由78.8%降至72.0%,测试集上AUC面 积由0.825降至0.793(参见图8.b),模型正确率、精确率及F1分值均显著降低(表2),由此表明空间 集聚特征能有效提高模型评估精度、增强模型评价性能。
表2剔除聚类属性前后模型在测试集上的混淆矩阵(以RF为例)
同样统计剔除后的危险性地图中不同危险区的面积占比及不同危险区内的灾害点占比。发现剔除后高 危区与低危区的面积占比显著下降,分别由原来的25%、47%下降至18%、34%(参见图8.d),部分灾害 多发地(如纳雍县南侧、织金县东南侧)未能识别、局部灾害少发处(如大方县南侧)过高识别(参见图 8.c),存在估计有偏、重点危险区把控不准的问题。再者,剔除后高危区对实际历史灾害点的解释率下降 了15%(原来75%,剔除后60%),低危区尽管对灾害点的解释率降低了4%(原来12%,剔除后8%), 但面积减少了13%(参见图8.d),可见剔除聚类属性后模型的实际评测性能有较大削减。故空间集聚特征对提升模型的性能具有显著的正向效用。
步骤F2:厘定空间分散特征对模型性能的影响。不对历史灾害点数据做任何处理,直接划分训练集与 测试集,输入模型重新评估,通过对比数据预处理前后的混淆矩阵与精度评价指标、危险性地图的空间分 布与评价效果来评估空间分散特征的负面作用。
为了验证基于渔网格网数据预处理方法的可行性,不对历史灾害点数据做任何处理,将灾害点仍按照 70%、30%的比例划分训练集与测试集,同时生成等量的非灾害点(参见图9.a)。结果表明,四种算法的 精度均显著下降,RF模型的训练精度仅达66.3%,测试集上AUC面积降至0.716(参见图9.b),正确率、 召回率等精度统计指标均大幅度下降近10%(参见表3),可见基于渔网的数据处理方案对于提升模型精 度非常有效。
表3数据预处理前后模型在测试集上的预测精度(以RF为例)
同样统计不采用渔网数据时的危险性地图中不同危险区的面积占比及不同危险区内的灾害点占比。发 现高危区与低危区面积显著下降,分别减少了6%、13%,极高易发区对历史灾害的解释率大幅度降低(降 低11%)(参见图9.d),中易发区占比陡增至47%,空间分布上存在低危区高估、高危区低估的问题(参 见图9.c),且未体现出渐变性特征。故本发明提出的渔网数据预处理方案能够显著提高模型预测精度、 优化空间评价结果。
步骤G:地质灾害危险性地图的确定及县级地质灾害危险性水平分析。依据上述评价及检验过程,以 评价质量最佳的模型生成结果作为最终的地质灾害危险性地图,并依据发生概率的数值分布与不同危险区 面积分布进一步地分析县级水平的地质灾害危险性。
结合各模型在测试集上的预测精度和最终的危险性地图评价效果,本实施例以RF的预测结果作为最 终的毕节市地质灾害危险性地图。拟从发生概率数值统计与危险区面积占比角度进一步分析县级水平的地 质灾害危险性。
基于各县地质灾害发生概率的数值分布情况(参见图10.a),可依据由高到低的危险性水平将各县划 分为三个梯队:第一梯队(纳雍县、黔西县、金沙县、七星关区)、第二梯队(织金县、大方县、赫章县)、 第三梯队(威宁县)。从小提琴宽度(参见图10.a)及扇形面积占比(参见图10.b)可以看出,第一梯队 的区县具有高危区面积多、占比大的特点,应予以重点监测防范。第二梯队的区县整体危险性有所降低, 但高危区在局部区域的聚集值得关注,以赫章县为代表,其极高危险区面积依然广阔(参见图10.b),因 此在灾害防御与治理环节中,这些局部区块应成为相关部门的重点考察对象。第三梯队的威宁县领域广阔, 中部地带海拔高、灾害较少,具有鲜明的危险性两极分化的特点,类似地,其东南侧与西北侧的部分高危 区应是未来防灾减灾规划管控的主要工作重心。
根据本发明实施例的考虑空间分布特征的地质灾害危险性评价装置100,包括历史地质灾害数据模块 110、地质灾害聚类属性提取模块120、数据预处理模块130、评价因子提取模块140、机器学习模块150、 以及地质灾害危险性评价模块160。该评价装置100可以执行上述参考图1描述的评价方法。
具体的,历史地质灾害数据模块110,用于针对特定区域的历史地质灾害数据的输入与存储。
地质灾害聚类属性提取模块120,用于利用聚类算法依据所述历史地质灾害数据模块110中存储的历 史地质灾害情况提取地质灾害聚类属性因子。
更进一步地,提取地质灾害聚类属性因子时,将不同的搜索距离与最小聚类要素个数参数组合下得到 的地质灾害聚类属性因子输入机器学习模块150中,通过对比不同参数组合下各机器学习模型在训练集上 的训练精度,对比得出最佳方案,即最优搜索距离与最优最小聚类要素个数,从而最终确定地质灾害聚类 属性因子。
数据预处理模块130,用于对所述历史地质灾害数据模块110中存储的历史灾害情况数据进行预处理, 并生成区域地质灾害危险性评价训练集及测试集,其中,所述训练集及所述测试集均包括多个灾害样本点 及多个非灾害样本点。
更进一步地,对历史灾害情况数据进行预处理时,基于渔网格网过滤历史地质灾害噪声点。
评价因子提取模块140,用于提取其他评价因子,及对所述其他评价因子的选取进行合理性评价。针 对选定的所述其他评价因子,采用不规则研究区的空间属性叠加手段提取,基于因子共线性分析方法对所 述其他评价因子的选取进行合理性评价。
机器学习模块150,用于多机器学习模型的构建、训练、预测精度评价,并采用所述多机器学习模型 输出地质灾害危险性地图。
更进一步地,基于不同机器学习算法的原理构建多机器学习模型,构建模型时,输入数据预处理模块 130中获得的所述训练集与地质灾害聚类属性提取模块120中确定的所述聚类属性评价因子、评价因子提 取模块140中提取的所述其他评价因子数据集,训练基于不同多机器学习算法的所述多机器学习模型。
更进一步地,基于数据预处理模块130中获得的所述测试集评价各所述多机器学习模型的预测精度。
更进一步地,基于所述多机器学习模型预测研究区内所有栅格单元的地质灾害危险性概率值(p), 按照不同危险性层级水平划分评价阈值:极低危险区(0≤p≤0.1)、低危险区(0.1<p≤0.3)、中危 险区(0.3<p≤0.6)、高危险区(0.6<p≤0.9)、极高危险区(0.9<p≤1),输出地质灾害危险性 地图。
更进一步地,机器学习模块150中还可以包括地质灾害危险性地图评估效果识别与检验模块,用于对 所述地质灾害危险性地图的评价效果识别与评估质量检验,得到最优地质灾害危险性地图。
地质灾害危险性评价模块160,用于县级地质灾害危险性水平分析;在所述最优地质灾害危险性地图 的基础上,依据发生概率的数值分布与不同危险区面积分布进一步地分析县级水平的地质灾害危险性。
最后应说明的是:以上实施例的描述尽管较为详细,但仅用以说明本发明的技术方案,而非对其限制, 本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,依然可以对前 述实施例中的技术细则进行修改,或者对其中部分技术特征进行等同替换,而这些修改或者替换,均落入 本发明的保护范围之内。
Claims (12)
1.一种考虑空间分布特征的地质灾害危险性评价方法,包括:
步骤A:利用聚类算法依据历史地质灾害情况数据提取地质灾害聚类属性作为区域地质灾害危险性评价的评价因子;
步骤B:对历史地质灾害情况数据进行预处理,基于渔网格网过滤历史地质灾害噪声点,生成区域地质灾害危险性评价训练集及测试集;
步骤C:选取地质灾害危险性评价的其他评价因子,采用不规则研究区的空间属性叠加手段提取所述其他评价因子;
步骤D:构建多机器学习模型,并依据步骤A中确定的所述地质灾害聚类属性和步骤C中确定的所述其他评价因子、以及步骤B中获得的所述训练集训练所述模型;输出地质灾害危险性地图;
以及,确定最优地质灾害危险性地图,并依据所述最优地质灾害危险性地图对地质灾害危险性进行评价。
2.根据权利要求1所述的方法,其特征在于,步骤A中,所述地质灾害聚类属性的提取流程为:
①针对历史地质灾害情况数据输入搜索距离、最小聚类要素个数,经OPTICS算法计算得到灾害点的聚类属性;
②根据每个历史地质灾害点的位置构建泰森多边形,并将各所述灾害点的聚类属性赋值给相应的泰森多边形;
③将泰森多边形矢量图层转栅格,得到聚类属性评价因子,与步骤C中获得的所述其他评价因子一同输入步骤D中的所述模型训练;
④整理不同搜索距离与最小聚类要素个数参数组合下各所述模型在所述训练集上的训练精度,对比得出最佳方案,即最优搜索距离与最优最小聚类要素个数,并以此得到所述地质灾害聚类属性。
3.根据权利要求1所述的方法,其特征在于,步骤C中,还包括基于因子共线性分析方法对所述其他评价因子的选取进行合理性评价。
4.根据权利要求1所述的方法,其特征在于,步骤D中,基于逻辑回归(LR)、支持向量机(SVM)、梯度提升树(GBDT)与随机森林(RF)的多机器学习算法构建所述模型。
5.根据权利要求4所述的方法,其特征在于,步骤D中,还包括基于步骤B中获得的所述测试集对各所述模型的预测精度进行评价。
6.根据权利要求5所述的方法,其特征在于,预测精度评价指标包括:精确率(Precision)、召回率(Recall)、正确率(Accuracy)和F1分值(F1 score),ROC曲线与AUC面积。
7.根据权利要求1、4-6中任一项所述的方法,其特征在于,还包括步骤E:对所述地质灾害危险性地图的评价效果识别与评估质量检验。
8.根据权利要求7所述的方法,其特征在于,从全局尺度上不同危险区的空间分布格局与面积占比,以及局部尺度上空间细节刻画水平、渐变性与斑块化特征角度,对比评价不同所述模型的评价效果,结合历史灾害点落入各层级危险区的比例即危险区对历史灾害的解释率进一步评估所述地图的质量。
9.根据权利要求7所述的方法,其特征在于,依据对所述地质灾害危险性地图的评价效果识别与评估质量检验的结果,确定所述最优地质灾害危险性地图。
10.根据权利要求1所述的方法,其特征在于,还包括步骤F:地质灾害空间分布特征考量的可靠性与实用性检验,包括厘定空间聚集特征对所述模型性能的影响、厘定空间分散特征对所述模型性能的影响。
11.根据权利要求1所述的方法,其特征在于,对地质灾害危险性进行评价时包括依据发生概率的数值分布与不同危险区面积分布分析各级水平的地质灾害危险性。
12.一种考虑空间分布特征的地质灾害危险性评价装置,包括:
历史地质灾害数据模块(110),用于存储区域历史地质灾害数据;
地质灾害聚类属性提取模块(120),用于利用聚类算法依据所述历史地质灾害数据模块(110)中存储的历史地质灾害情况提取地质灾害聚类属性因子;
数据预处理模块(130),用于对所述历史地质灾害数据模块(110)中存储的历史灾害情况数据进行预处理,基于渔网格网过滤历史地质灾害噪声点,并生成区域地质灾害危险性评价训练集及测试集;
评价因子提取模块(140),用于提取选取的地质灾害危险性评价的其他评价因子;
机器学习模块(150),用于多机器学习模型的构建、训练,依据所述地质灾害聚类属性因子和所述其他评价因子、所述训练集训练所述模型,并利用所述模型输出地质灾害危险性地图;
地质灾害危险性评价模块(160),用于依据所述地质灾害危险性地图进行地质灾害危险性水平分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110814485.8A CN113642849B (zh) | 2021-07-19 | 2021-07-19 | 考虑空间分布特征的地质灾害危险性综合评价方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110814485.8A CN113642849B (zh) | 2021-07-19 | 2021-07-19 | 考虑空间分布特征的地质灾害危险性综合评价方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113642849A true CN113642849A (zh) | 2021-11-12 |
CN113642849B CN113642849B (zh) | 2023-10-20 |
Family
ID=78417706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110814485.8A Active CN113642849B (zh) | 2021-07-19 | 2021-07-19 | 考虑空间分布特征的地质灾害危险性综合评价方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113642849B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113989859A (zh) * | 2021-12-28 | 2022-01-28 | 江苏苏宁银行股份有限公司 | 一种防刷机设备指纹相似度识别方法和装置 |
CN114091756A (zh) * | 2021-11-23 | 2022-02-25 | 国家海洋环境预报中心 | 一种基于泰森多边形的乡镇级海啸风险评估方法 |
CN114117827A (zh) * | 2022-01-24 | 2022-03-01 | 四川省公路规划勘察设计研究院有限公司 | 一种雁列式断层模型的随机模拟及参数敏感度分析方法 |
CN114595629A (zh) * | 2022-03-01 | 2022-06-07 | 合肥工业大学 | 基于gbdt-lr和信息量模型的滑坡易发性评价方法 |
CN115049053A (zh) * | 2022-06-20 | 2022-09-13 | 航天宏图信息技术股份有限公司 | 基于TabNet网络的黄土地区滑坡易发性评估方法 |
CN115271564A (zh) * | 2022-09-29 | 2022-11-01 | 淄博市交通建设发展中心 | 一种公路边坡灾害空间危险性区划方法及终端机 |
CN115328373A (zh) * | 2022-08-02 | 2022-11-11 | 北京市地质灾害防治研究所 | 一种确定性系数计算方法及系统 |
CN115688404A (zh) * | 2022-10-24 | 2023-02-03 | 中国科学院空天信息创新研究院 | 基于svm-rf模型的降雨型滑坡预警方法 |
CN117540830A (zh) * | 2024-01-05 | 2024-02-09 | 中国地质科学院探矿工艺研究所 | 基于断层分布指数的泥石流易发性预测方法、装置及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111582386A (zh) * | 2020-05-11 | 2020-08-25 | 四川师范大学 | 一种基于随机森林的地质灾害多灾种综合风险评价方法 |
CN112819207A (zh) * | 2021-01-19 | 2021-05-18 | 武汉中地云申科技有限公司 | 基于相似性度量的地质灾害空间预测方法、系统及存储介质 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006266858A (ja) * | 2005-03-24 | 2006-10-05 | Railway Technical Res Inst | 硬質堆積岩地域の災害発生危険度評価方法 |
CN102645650A (zh) * | 2012-03-06 | 2012-08-22 | 北京北科安地科技发展有限公司 | 一种基于D-InSAR干涉差分的滑坡动态识别及监测技术 |
KR101663679B1 (ko) * | 2014-11-17 | 2016-10-10 | 서울대학교산학협력단 | 기하학적 증폭계수를 고려한 지진시 산사태 재해도 도출방법 및 도출시스템 |
CN107463991A (zh) * | 2017-06-28 | 2017-12-12 | 西南石油大学 | 一种基于坡体单元与机器学习的区域滑坡危险性评价方法 |
CN111858803B (zh) * | 2020-07-06 | 2021-09-21 | 东华理工大学 | 一种滑坡地灾风险区划图生成方法 |
CN112561297B (zh) * | 2020-12-09 | 2024-07-05 | 中南大学 | 不同地形下邻域影响的地质灾害风险评价方法 |
-
2021
- 2021-07-19 CN CN202110814485.8A patent/CN113642849B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111582386A (zh) * | 2020-05-11 | 2020-08-25 | 四川师范大学 | 一种基于随机森林的地质灾害多灾种综合风险评价方法 |
CN112819207A (zh) * | 2021-01-19 | 2021-05-18 | 武汉中地云申科技有限公司 | 基于相似性度量的地质灾害空间预测方法、系统及存储介质 |
Non-Patent Citations (3)
Title |
---|
YUMIAO WANG ETAL: ""A hybrid model considering spatial heterogeneity for landslide susceptibility mapping in Zhejiang Province, China"" * |
YUMIAO WANG ETAL: ""A hybrid model considering spatial heterogeneity for landslide susceptibility mapping in Zhejiang Province, China"", 《CATENA》 * |
杨斌: ""一种不确定NNSB-OPTICS聚类算法在滑坡危险性预测研究"" * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114091756A (zh) * | 2021-11-23 | 2022-02-25 | 国家海洋环境预报中心 | 一种基于泰森多边形的乡镇级海啸风险评估方法 |
CN114091756B (zh) * | 2021-11-23 | 2022-07-05 | 国家海洋环境预报中心 | 一种基于泰森多边形的乡镇级海啸风险评估方法 |
CN113989859A (zh) * | 2021-12-28 | 2022-01-28 | 江苏苏宁银行股份有限公司 | 一种防刷机设备指纹相似度识别方法和装置 |
CN113989859B (zh) * | 2021-12-28 | 2022-05-06 | 江苏苏宁银行股份有限公司 | 一种防刷机设备指纹相似度识别方法和装置 |
CN114117827A (zh) * | 2022-01-24 | 2022-03-01 | 四川省公路规划勘察设计研究院有限公司 | 一种雁列式断层模型的随机模拟及参数敏感度分析方法 |
CN114117827B (zh) * | 2022-01-24 | 2022-04-15 | 四川省公路规划勘察设计研究院有限公司 | 一种雁列式断层模型的随机模拟及参数敏感度分析方法 |
CN114595629A (zh) * | 2022-03-01 | 2022-06-07 | 合肥工业大学 | 基于gbdt-lr和信息量模型的滑坡易发性评价方法 |
CN115049053B (zh) * | 2022-06-20 | 2023-03-24 | 航天宏图信息技术股份有限公司 | 基于TabNet网络的黄土地区滑坡易发性评估方法 |
CN115049053A (zh) * | 2022-06-20 | 2022-09-13 | 航天宏图信息技术股份有限公司 | 基于TabNet网络的黄土地区滑坡易发性评估方法 |
CN115328373A (zh) * | 2022-08-02 | 2022-11-11 | 北京市地质灾害防治研究所 | 一种确定性系数计算方法及系统 |
CN115328373B (zh) * | 2022-08-02 | 2024-01-12 | 北京市地质灾害防治研究所 | 一种确定性系数计算方法及系统 |
CN115271564A (zh) * | 2022-09-29 | 2022-11-01 | 淄博市交通建设发展中心 | 一种公路边坡灾害空间危险性区划方法及终端机 |
CN115271564B (zh) * | 2022-09-29 | 2022-12-27 | 淄博市交通建设发展中心 | 一种公路边坡灾害空间危险性区划方法及终端机 |
CN115688404A (zh) * | 2022-10-24 | 2023-02-03 | 中国科学院空天信息创新研究院 | 基于svm-rf模型的降雨型滑坡预警方法 |
CN117540830A (zh) * | 2024-01-05 | 2024-02-09 | 中国地质科学院探矿工艺研究所 | 基于断层分布指数的泥石流易发性预测方法、装置及介质 |
CN117540830B (zh) * | 2024-01-05 | 2024-04-12 | 中国地质科学院探矿工艺研究所 | 基于断层分布指数的泥石流易发性预测方法、装置及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113642849B (zh) | 2023-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113642849A (zh) | 考虑空间分布特征的地质灾害危险性综合评价方法及装置 | |
Polykretis et al. | Comparison and evaluation of landslide susceptibility maps obtained from weight of evidence, logistic regression, and artificial neural network models | |
Cao et al. | Multi-geohazards susceptibility mapping based on machine learning—a case study in Jiuzhaigou, China | |
Can et al. | Landslide susceptibility mapping at Ovacık-Karabük (Turkey) using different artificial neural network models: comparison of training algorithms | |
Jochem et al. | Identifying residential neighbourhood types from settlement points in a machine learning approach | |
Guzzetti et al. | Estimating the quality of landslide susceptibility models | |
Wan | A spatial decision support system for extracting the core factors and thresholds for landslide susceptibility map | |
CN115688404B (zh) | 基于svm-rf模型的降雨型滑坡预警方法 | |
CN110929939B (zh) | 一种基于聚类-信息量耦合模型下的滑坡灾害易发性空间预测方法 | |
Wan | Entropy-based particle swarm optimization with clustering analysis on landslide susceptibility mapping | |
CN113487105B (zh) | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 | |
Ge et al. | A comparison of five methods in landslide susceptibility assessment: a case study from the 330-kV transmission line in Gansu Region, China | |
CN114330812A (zh) | 基于机器学习的滑坡灾害风险评估方法 | |
CN109118004B (zh) | 一种工程构筑选址适宜区预测方法 | |
CN112907113A (zh) | 一种考虑空间相关性的植被变化成因识别方法 | |
Okalp et al. | Landslide susceptibility assessment in medium-scale: case studies from the major drainage basins of Turkey | |
Huang et al. | Modelling landslide susceptibility prediction: a review and construction of semi-supervised imbalanced theory | |
Zou et al. | Quantitative assessment of geological hazard risk with different hazard indexes in mountainous areas | |
Ilmy et al. | Application of machine learning on google earth engine to produce landslide susceptibility mapping (case study: Pacitan) | |
CN116258279B (zh) | 基于综合赋权的滑坡易发性评价方法和装置 | |
CN116110210B (zh) | 复杂环境下数据驱动的滑坡灾害辅助决策方法 | |
CN110674471A (zh) | 基于GIS与Logistic回归模型的泥石流易发性预测方法 | |
Gouda et al. | Data mining for weather and climate studies | |
Bai et al. | Deep semantic segmentation for rapid extraction and spatial-temporal expansion variation analysis of China’s urban built-up areas | |
Asmael | A GIs Based Weight of Evidence for Prediction Urban Growth of Baghdad City by Using Remote Sensing Data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |