指南中土地资源单项评价要求:
利用全域DEM计算地形坡度,按≤3°、38°、815°、15~25°、>25°生成坡度分级图,将城镇建设土地资源划分为高、较高、中等、较低、低5级。
对于高程≥5000m的区域,土地资源等级直接取最低等;高程在3500~5000m之间的,将坡度分级降1级作为土地资源等级。
计算地形起伏度。邻域范围通常采用20公顷左右(如50m×50m栅格建议采用9×9邻域,30m×30m栅格建议采用15×15邻域),对于地形起伏度>200m的区域,将初步评价结果降2级,地形起伏度在100~200m之间的,将初步评价结果降1级作为城镇土地资源等级。
过程:首先计算坡度,使用【重分类】工具对坡度数据重分类为5级得到土地资源初评结果,然后用【栅格计算器】工具结合高程进行修正,接着用【焦点统计】工具计算地形起伏度,再次用【栅格计算器】结合地形起伏度进行第二次修正,最后得到土地资源评价结果。
数据:成都市30mdem
地形坡度计算评价
过程就不记录啦,【表面分析】-【坡度】,扔进dem计算坡度得到左图,【重分类】按阈值分为5级得到右图。接下来是根据高程进行修正。
地形高程进行修正
按照分级阈值用【重分类】对dem进行分级,根据高程分为3级,小于3500为3,3500-5000为2,5000以上为1。
接着用【栅格计算器】将两组数据相加,用高程去修正坡度
计算得到下图左图结果,然后根据指南和表格,碰到小于3500m的初评结果等级不变,碰到3500-5000m的降1级,碰到5000m以上的直接取最低等级1。用【重分类】据此再分类,得到高程修正的结果,右图。
地形起伏度计算并修正
【邻域分析】-【焦点统计】将dem扔进去计算地形起伏度,领域设置根据指南要求设置为15×15(下面的图我写错了是15才对,懒得重新截图了,注意),统计类型选择RANGE。
起伏度计算出来了,接着还是进行重分类,按照指南中的要求根据计算出来的结果分级,并重分类为3类。
(套娃)接着用【栅格计算器】将高程修正结果与地形起伏度分类结果两组数据相加,用起伏度去修正高程的修正。碰到起伏度大于200的,降2级作为最终评价结果,碰到起伏度在100-200的,降1级作为最终评价结果,碰到100以下的就是等级不变。
步骤截图不再重复,最终结果如下