[发明专利]一种基于GIS的多准则决策分析绘制滑坡敏感性图的方法

权利要求书

1.一种基于GIS的多准则决策分析绘制滑坡敏感性图的方法，其特征在于：包括以下步骤：

S1：基于数字高程模型(DEM)绘制岩性、坡度、坡向、土地覆盖、排水密度、地形湿度指数、高程、坡长和到道路的距离等9个参数的专题地图，使用ArcGIS软件包注册到通用横向墨卡托投影系统；

S2：使用层次分析方法(AHP)方法对滑坡相关参数进行标准化，利用Idrisi Taiga软件的“模糊”模块，基于模糊隶属函数对滑坡因子图进行标准化；

S3：采用基于GIS的多准则分析法中的层次分析法(AHP)确定岩性，坡度，坡向等九个标准的权重，分析其相对重要性；

S4：基于WLC函数，使用Idrisi Taiga软件MCE模块中可用的WLC方法计算每个评估因子；

S5:将各个标准滑坡因子图通过多准则决策分析结果进行彼此堆叠，形成多层图像，生成滑坡敏感性图。

2.如权利要求1所述的基于GIS的多准则决策分析绘制滑坡敏感性图的方法，其特征在于：步骤S1基于数字高程模型(DEM)绘制岩性、坡度、坡向、土地覆盖、排水密度、地形湿度指数、高程、坡长和到道路的距离等9个参数的专题地图，使用ArcGIS软件包注册到通用横向墨卡托投影系统包括以下内容：

S101：绘制岩性专题地图，可根据目标区域当地地质勘察部门发布的比例尺地图进行绘制；

S102：绘制坡度专题地图，可使用各个不同比例尺的地形图通过数字高程模型(DEM)从1∶25生成；

S103：绘制坡向专题地图，与坡向相关的参数，如阳光照射、干燥风、降雨(饱和度)和不连续性，是触发滑坡的重要因素，通过DEM图像用于计算每个像素的纵横比值，以构建纵横比图，生成坡向专题地图；

S104：绘制土地覆盖专题地图，利用卫星数据制作了30米分辨率的土地覆盖图，确定了覆盖目标区域的土地覆盖类型；

S105：绘制排水密度的专题地图：排水密度是指流域内所有溪流和河流的总长度除以流域总面积，由以下公式计算，

D_y＝∑L/A

式中，D_y为排水密度，L为河流长度，A为集水区面积。为每个集水区估计DY，即，同一集水区内所有像素的排水密度值均为常数。使用ArcGIS软件中的适当算法，从DEM生成目标区域的排水密度图；

S106：绘制地形湿度指数专题地图：使用无凹陷DEM算法去除光栅图像中的汇，根据生成的DEM图像确定多个流向后，生成流量累积面积(A_S)和坡度切线(tanβ)图像，目标区域的地形湿度指数TWI是使用以下公式创建的，TWI＝in(A_S/tanβ)，β为坡度；

S107：绘制高程专题地图：基于目标区域的DEM高程数据，按500m等距进行分类，绘制高程专题地图；

S108：绘制坡长专题地图：根据目标区域的DEM绘制坡长图，然后以25m等间距将其划分为若干个等级；

S109：绘制距道路距离专题地图：于距离道路中心线25米的距离间隔，通过缓冲区分析从现有形状文件创建到路线图的距离；

S1010：将专题地图使用ArcGIS软件包注册到通用横向墨卡托投影系统，转化为滑坡因子图。

3.如权利要求2所述的基于GIS的多准则决策分析绘制滑坡敏感性图的方法，其特征在于：步骤S2使用层次分析方法(AHP)方法对滑坡相关参数进行标准化，利用Idrisi Taiga软件的“模糊”模块，基于模糊隶属函数对滑坡因子图进行标准化包括以下内容：

由于代表滑坡相关参数或因素的输入层处于不同的尺度或间隔，因此使用AHP方法对其进行标准化。将利用Idrisi Taiga软件的“模糊”模块，基于模糊隶属函数对滑坡因子图进行标准化。在因子图的标准化中，每个因子都使用了S形模糊隶属函数和用户定义函数。坡度、TWI和排水密度图像采用单调递增的S形模糊隶属函数；距道路距离，坡长和高程图像采用单调递减的S形模糊隶属函数；用户定义的隶属函数用于岩性、土地覆盖和坡向图像。