基于频率比模型采用arcgis的滑坡易发性评估

参考视频：

1. ArcGIS滑坡易发性评价_哔哩哔哩_bilibili
2. landslides susceptibility assessemnt using frequently ration approach /滑坡敏感性制图_哔哩哔哩_bilibili
3. Frequency Ratio Method part-2 How to calculate landslide susceptibility in Arcgis (youtube.com)

频率比模型

地形因子处理

确定研究区，建立面文件AOI，在文件目录，右键-新建面，开始编辑，点击绘制草图，绘制面。
按掩膜提取，使用研究区shp裁剪Dem.tif
使用投影栅格将dem从WGS84转换为墨卡托投影
使用要素-》投影将AOI.shp从WGS84转换为墨卡托投影
基于墨卡托投影的dem使用3D分析-栅格表面-》坡度，计算坡度；
基于墨卡托投影的dem使用3D分析-栅格表面-》坡向，计算坡向；
基于墨卡托投影的dem使用3D分析-栅格表面-》曲率，计算曲率或者剖面曲率、平面曲率；
计算地形起伏度：
ArcGIS地形起伏度提取及出图设置-CSDN博客
打开-spatial analyst工具-邻域分析-焦点统计，输入栅格为裁剪的DEM，统计类型为MAXINUM，即计算邻域内的最大值

再进行一次焦点统计操作，这次统计类型为MINNUM，即计算邻域内的最小值

打开栅格计算器，输入的代码为使刚才邻域分析得到的最大值图层减去最小值图层

以dem的像元大小为基准进行后续处理
坡度、坡向、曲率、地形起伏度的像元大小，范围和dem应该相同，
判断这几个tif的栅格总数应该相同。

shp文件与缓冲区

断层描绘shp线文件：
地理校准，使用至少三个点进行地理校准
Shp建立缓冲区
地图处理-》缓冲区
1000m，2000m，3000m，缓冲区……


擦除：相当于取反集。
用3000m擦除研究区得到大于3000m的区域，用2000擦除3000的到2000_{3000，用1000擦除2000的到1000}2000。



合并缓冲区：合并1000, 1000_2000，20003000，>3000的shp, 得到合并的四级缓冲区，

缓冲区会大于研究区，地图处理-》裁剪，用研究区裁剪缓冲区，将裁剪后的缓冲区使用

面转栅格转为tif文件，tif文件若大于研究的栅格数目，多次使用按掩膜提取进行处理。

地理校准

描绘边界线——要素转面，用地图处理-》裁剪，用研究区裁剪转为的面区，面转栅格
滑坡面积shp文件，为滑坡区，用滑坡面积shp擦除研究区得到非滑坡区，使用地图处理-》合并，合并滑坡区与非滑坡区，使用面转栅格，将合并的滑坡区与非滑坡区转化为tif，使用重采样分类数据，将滑坡区标记为1，非滑坡区标记为2。

填补缺失值

对所有因子统计栅格总数，若存在有栅格总数明显多或者明显少，进行多次掩膜提取或者填充缺失值。
栅格计算器：

Con(IsNull(“raster”), FocalStatistics(“raster”, NbrRectangle(5,5, “CELL”), “MEAN”), “raster”)
“raster”为目标栅格

【ArcGIS】栅格缺失值填补_栅格数据有空白怎么补上-CSDN博客

对所有因子进行统一重分类：
搜索重分类，点击栅格，重分类使用自然间断点分类或者手动间隔重分类；
不勾选缺失值填充nodata

打开所有重分类因子，将重分类因子栅格统计复制到excel：
全选—复制所选项，将全部因子如此操作

滑坡区域在每一子类所占的面积：
仅将滑坡区转化为tif文件，并添加一个字段object，将其值设置为1。
面积制表：
选中子类因子，选中仅滑坡tif，值选择object，导出表。
计算频率比-相对比或者信息量

PIEXL OF COUNT=OBIEC_`1/(cell_x*cell_y)

对已经重分类的各类栅格数据，重新赋值RF，这里RF是取整之后的RF，用PR乘以各个重分类后的因子再相加即可得到易发性值，对易发性栅格进行重分类。

利用模糊隶属度进行归一化

易发性结果

基于栅格计算器和模糊隶属度归一化后可以得到归一化的滑坡易发性结果

出图处理：

对连续dem计算山体阴影

ArcGIS晕渲图的制作（解决山体阴影问题——平坦地区凸点过于突出）_arcgis山体阴影不明显-CSDN博客
3D分析-栅格表面-》坡度，计算坡度；
叠加了山体阴影的地形与滑坡