基于资源环境承载力和国土空间开发适宜性双评价的建设土地资源评价操作步骤

因为本系列只是起研究模拟作用，因此以广东省（省级）为例，地理坐标系使用WGS1984（许多网上得到的数据都是WGS84，所以不使用CGCS2000），投影使用WGS_1984_UTM_Zone_49N，空间分辨率设为90m。

首先，对广东省水资源、土地资源、气候、生态、环境、灾害进行评价。

今天的操作是建设指向的土地资源评价！

1-1-2 广东省城镇建设土地资源评价

然后我们进行城镇建设功能指向的土地资源评价：

首先看公式：

城镇指向的土地资源评价：[城镇建设条件]＝f ([坡度], [高程], [地形起伏度])

其次看步骤要求（坐标系与精度与双评价要求不同，主要操作步骤相同，参考双评价2019.6月版）：

第一步：坡度：利用DEM，计算地形坡度，一般按≤3°、3~8°、8~15°、15~25°、>25°生成坡度分级图。

第二步：土地资源评价与分级。以坡度分级结果为基础，结合高程，划分城镇建设条件高、较高、中等、较低、低 5 级。将高程≥5000m 区域，城镇土地资源等级直接取最低等级；高程在3500~5000m 之间的，将坡度分级降1 级作为城镇土地资源等级。

第三步：地形复杂地区评价结果修正。在地形起伏剧烈的地区，进一步通过地形起伏度指标对城镇土地资源等级进行修正。

利用DEM 邻域分析功能计算地形起伏度，邻域范围通常采用20公顷左右（如50m×50m 栅格建议采用9×9 邻域，30m×30m栅格建议采用15×15 邻域），对于地形起伏度>200m 的区域，将评价结果降 2 级作为城镇土地资源等级，地形起伏度在100~200m 之间的，将评价结果降1 级作为城镇土地资源等级。

首先，我们按照坡度进行重分类，得到jsslope:

接着本来我们是需要根据高程，对坡度分级的建设土地资源分级进行修正，但是广东省内最高点不足2000m，因此略过第二步，直接进行第三步起伏度修正。

那么首先运用DEM进行起伏度分析（使用焦点统计工具），需要明白领域范围采用20公顷的话，我们的DEM精度为90m，因此采用5*5的邻域，得到range：

然后进行重分类，>200m的降2级，100~200m之间的降1级，采用和农业土地资源相同的重分类方法，得到range2：

这个时候再使用栅格计算器，进行相加，可以预见，会产生10,11,12,20,21,22…50,51,52，那么，个位数为0的区域，即不变等级，个位数为1的，下降1级，个位数为2的，下降2级，得到jstdpj：

然后我们再次进行重分类，得到jstdpj1：

那么很明显的，5-1就是城镇建设指向的土地资源等级从高-低了，我们根据制图规范、配色以及比例尺、指北针，进行一个出图：

这样我们就得到了建设指向的土地资源分级图，当然由于坐标、数据精度等问题，与实际情况可能有很大出入，因此本文只是步骤流程的演示，并不对结果的准确性进行保证哦！