国土空间规划双评价（一）：土地资源单项评价

理论：
（1）评价时需要扣除河流、湖泊等水域，利用DEM计算地形坡度，按≤3°、3-8°、8-15°、15~25°和＞25°生成坡度分级图，将城镇建设用地资源划分为高、较高、中等、较低、低五级。
（2）高程≥5000m的区域，土地等级取最低等（比如我在第一步中将等级划分为高，但高程≥5000m，直接归为低）；3500-5000m的，将坡度分级降1级作为土地资源等级。
（3）计算地形起伏度。邻域范围通常采用20公顷 （20w平方米哦，所以50m50m栅格建议采用99也就是81个格子），对于地形起伏度＞200m的区域，将初步评价结果（前两步评定完的结果）降2级，地形起伏度在100~200m之间的，初步评价结果降1级作为城镇土地资源等级。

操作过程：
一、坡度分级
利用DEM数据生成坡度（spatial analyst—表面分析—坡度）；
重分类：spatial analyst—表面分析—坡度
高、较高、中等、较低、低分别对应5, 4，3，2，1

二、利用高程进行修正
（1）按照理论（2）对高程进行重分类
1,2分别对应0-3500m，3500-4250m
（2）利用栅格计算器进行修正
第一次进行降级
Con(“DEM_reclass.tif” == 2,“slope_reclass.tif” -1,“slope_reclass.tif”)
第二次把从1级降为0的拉回1
Con(“land_modify11.tif” == 0,1,“land_modify11.tif”)

三、计算地形起伏度
(1)计算地形起伏度（spatial analyst—邻域分析—焦点统计）
备注：3030m—>20公顷（1515个格子）

(2)对地形起伏度重分类（1，2，3分别表示＞200,100-200，＜100m）
(3) 同样采用栅格计算器进行降级
将重分类后的地形起伏度*30+高程校正后的数据（“起伏度_reclass.tif” * 10+“land_modify.tif”）

31-35表示地形起伏度＜100m，不需要降级
21-25表示地形起伏度介于[100,200]之间，降1级
11-15表示地形起伏度大于200，降2级
重新分类