国土空间规划双评价操作教程-农业土地资源评价

本系列教程完全原创，笔名beatniu，转载请联系作者获取授权！！！

本系列基于自然资源部2019.3月和2020年1月印发的两版《资源环境承载能力和国土空间开发适宜性评价指南（试行）》（简称“双评价”），因为计算公式与方法在2020.1月版的少之又少，可以说不限制模型方法，只有规定大体的方向，但是都可以基于2019.6月的指标、模型、方法进行计算，因此一些定义、目标、流程、制图规范都基于新版的，但是方法都基于旧版的。

首先说明基于2020.1月版的“双评价”的定义、目标：

定义：（1）“资源环境承载力评价”：在一定发展阶段，经济技术水平和生产生活方式，一定地域范围内资源环境要素能够支撑的农业生产、城镇建设等人类活动的最大规模。

（2）“国土空间开发适宜性评价”： 在维系生态系统健康的前提下，综合考虑资源环境要素和区位条件以及特定国土空间，进行农业生产城镇建设等人类活动的适宜程度。

目标：分析区域资源禀赋与环境条件，研判国土空间开发利用问题和风险，识别生态保护极重要区（含生态系统服务功能极重要区和生态极脆弱区），明确农业生产、城镇建设的最大合理规模和适宜空间，为编制国土空间规划，优化国土空间开发保护格局，完善区域主题功能定位，划定三条控制线，实施国土空间生态修复和国土综合整治重大工程提供基础依据，促成形成以生态优先、绿色发展为导向的高质量发展新路子。

然后说明于2020.1月版的“双评价”的工作流程以及要求：

工作流程：

通过工作流程可看出，本系列将分为五大步：

第一步，首先要进行本底评价，也就是以往版本双评价中的1土地资源评价、2水资源评价、3气候评价、4生态评价、5环境评价、6灾害评价六个单项评价。需要提出的是，虽然新版的双评价似乎可以直接开始生态保护重要性评价，但是本系列依然先从6个本底评价开始，各位可以当作一个旧版的研究评价。

第二步：生态保护重要性评价。生态保护省级评价从区域生态安全底线出发，在陆海全域，评价水源涵养、水土保持、生物多样性维护、防风固沙、海岸防护等生态系统服务功能重要性，以及水土流失、石漠化、土地沙化、海岸侵蚀及沙源流失等生态敏感性，综合形成极重要区和重要区。

第三步：农业生产适宜性评价。省级在生态保护极重要区以外的区域，开展种植业、畜牧业、渔业等农业生产适宜性评价，识别农业生产适宜区和不适宜区。

城镇建设适宜性评价。省级在生态保护极重要区以外的区域，优先考虑环境安全、粮食安全和地质安全等底线要求，识别城镇建设不适宜区。沿海地区针对海洋开发利用活动开展评价。

第四步：承载规模评价。基于现有经济技术水平和生产生活方式，以水资源、空间约束等为主要约束，缺水地区重点考虑水平衡，分别评价各评价单元可承载农业生产、城镇建设的最大合理规模。各地可结合环境质量目标、污染物排放标准和总量控制等因素，评价环境容量对农业生产、城镇建设约束要求。按照短板原理，取各约束条件下的最小值作为可承载的最大合理规模。

第五步：综合分析，如资源禀赋分析、问题和风险识别、潜力分析和情景分析等。

最后说明于2020.1月版的“双评价”的制图规范：

要求：评价统一使用2000国家大地坐标系（CGCS2000），高斯-克吕格投影，陆域部分采用1985国家高程基准，海域部分采用理论深度基准面高程基准。制图规范、精度等参考省级国土空间规划要求。

制图规范：

那么再参考《省级国土空间规划编制指南》，明确制图规范与精度的要求：

字体：汉字宋体，西字Times New Roma。

比例尺：1:50万或1:100万。

指北针：

图例：图幅内左下角或右下角。

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

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

那么我们先从土地资源评价开始，这样第一步就开始了！

1-1-1 广东省农业土地资源评价

公式：农业指向的土地资源评价：[农业耕作条件]＝f ([坡度], [土壤质地])

农业指向的步骤要求（坐标系与精度与双评价要求不同，主要操作步骤相同）：

坡度：利用 DEM，计算地形坡度，按≤2°、2~6°、6~15°、15~25°、>25°划分为平地、平坡地、缓坡地、缓陡坡地、陡坡地 5 个等级，生成坡度分级图。

土壤质地：将土壤的粉砂含量≥80%区域，农业土地资源直接取最低等；60%≤粉砂土含量＜80%的区域，将坡度分级降 1 级作为农业土地资源等级。

首先我们需要的数据：DEM数据、土壤质地数据。

广东省DEM数据如下，空间分辨率为90m：

首先坡度分析，得到slope：

重分类，（至于为什么取50,40,30,20,10，在土壤质地分类后会讲到）得到nyslope，那么50,40,30,20,10依次代表平地、平坡地、缓坡地、缓陡坡地、陡坡地：

接下来需要根据土壤质地进行分类，那么在2020.1月版的双评价中，没有指明划分标准，因此参照2019.6月版的：

可以看出主要是从粉砂的含量进行划分，那么需要使用到HWSD世界土壤数据库，然后使用到其中的clay字段，clay指黏土，那么粗略的可以认为，土壤中若粉砂≥80%，则黏土≤20%，若粉砂≥60%，则黏土≤40%，因此，将土壤质地按照黏土含量进行分级，≤20%，取最低，20%

首先按掩膜提取广东省域的，然后需要进行投影以及重采样至90m分辨率，再次不重复类似操作：

然后进行连接，连接HWSD_DATA的MU_GLOBAL：

HWSD_DATA中的CLAY就是我们所需的黏土含量，然后我们新建字段，赋值后移除连接：

然后我们就可以根据得到的CLAY字段进行重分类了：

然后基于土壤质地，我们将坡度分级的农业资源等级进行降级，我们使用栅格计算器进行相加：

那么很明显的，我们会得到11,12,13,21,22,23…51,52,53这样的数值，那么个位数为1的，很明显就是黏土含量最少，也就是粉砂含量最多的，直接取最低等级；个位数为2的，降一级处理；个位数为3的，不降级，懂了吧！因此再进行一个重分类：

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

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