详细讲解 “双评价”——农业土地资源评价
本系列基于自然资源部2019年6月和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%<clay<40%降一级，≥40%不降级。HWSD世界土壤数据库如下：
首先按掩膜提取广东省域的，然后需要进行投影以及重采样至90m分辨率，再次不重复类似操作：
然后进行连接，连接HWSD_DATA的MU_GLOBAL：
HWSD_DATA中的CLAY就是我们所需的黏土含量，然后我们新建字段，赋值后移除连接：
然后我们就可以根据得到的CLAY字段进行重分类了：
然后基于土壤质地，我们将坡度分级的农业资源等级进行降级，我们使用栅格计算器进行相加：
那么很明显的，我们会得到11,12,13,21,22,23…51,52,53这样的数值，那么个位数为1的，很明显就是黏土含量最少，也就是粉砂含量最多的，直接取最低等级；个位数为2的，降一级处理；个位数为3的，不降级，懂了吧！因此再进行一个重分类：
那么很明显的，5-1就是农业指向的土地资源等级从高-低了，我们根据制图规范、配色以及比例尺、指北针，进行一个出图：
这样我们就得到了农业指向的土地资源分级图，当然由于坐标、数据精度等问题，与实际情况可能有很大出入，因此本文只是步骤流程的演示，并不对结果的准确性进行保证哦！
- END -
收藏 | 国土空间规划“三线”划定详解
国土空间规划中的GIS应用——通识篇
ArcGIS Python从入门到精通，你需要这个文档
收藏 | 这可能是最全的国土空间规划合集了
ArcGIS在国土空间规划编制中的应用
7.72GB ARCGIS ERDAS ENVI干货教程网盘下载