写在最前面

数读菌作为失踪人口，偶尔又出现了

今天的这篇文章，其实在5月初就写完了

但是一直没发

要问什么原因的话

其实就是懒癌发作，好不容易逼着自己学习了下

临到最后关头

只是上传一下的功夫也拖延到现在

以至于文中很多内容已经过时了

比如双评价的征求意见稿已经到7月版了

文中还是4月那一版为基础的

但是本着，写了不能浪费的精神

我还是鼓起勇气把这篇文章发出来

没准还能赚些广告费

明天可以买一把青菜吃

（手动微笑脸）

哦对

这两个月我一直没登陆

所以很多消息没回（也可能没消息）

反正我也看不到

看到了我也回不了

所以有什么后台留言可以再发一次

按照我的规律，发文的一周时间左右

我是肯定在的

（继续手动微笑）

什么是“双评价”？

◆ ◆ ◆

作为一个最近很热的词（和7月的天气比起来，还是7月天气更热），“双评价”是指资源环境承载能力和国土空间开发适宜性评价。引用《双评价技术指南（汇报论证稿2019年4月）》的原文：

资源环境承载能力和国土空间开发适宜性评价是国土空间规划编制的前提和基础，也是国土空间规划编制过程中系列研究分析的重要组成部分。通过资源环境承载能力评价，按照生态保护、农业生产、城镇建设三个不同功能指向，确定资源环境承载能力等级，总结区域资源环境禀赋特点，分析其优势与短板。通过国土空间开发适宜性评价，划分生态保护重要性等级分区，农业生产、城镇建设适宜性等级分区，识别国土空间开发保护的问题和风险，综合分析未来适宜国土空间开发保护的潜力，为完善主体功能区战略，编制国土空间规划，科学划定生态保护红线、永久基本农田、城镇开发边界，统一行使国土空间用途管制和生态保护修复提供基础支撑。

为什么要写这系列文章？

◆ ◆ ◆

因为自身技术无法跟上时代的焦虑吧。这系列文章是自己学习双评价的学习笔记，属于边学边写，所以不是教程，错误估计也不少，希望能和大家一起交流。本来想建个微信群，专门用来讨论国土空间规划的相关技术，但是鉴于我不知道会更新几篇文章，所以这个计划暂时先等我有一定的产出后再提上日程吧。（话是这么说啦，这系列文章可能成为还没开始就已经结束的典型案例，哎，我尽量吧，尽量能继续更。。。吧。。。）

学习的依据或者参考是什么？

◆ ◆ ◆

主要参考《双评价技术指南（汇报论证稿2019年4月）》，如果有更新的或者正式指南当然就会参考那个了。至少我不知道有什么专门的教程，所以只能以这个指南为主，参考一些已有的成果和报告咯。慢慢摸索咯。文中斜体引用的无特别说明的皆为引用自《双评价技术指南（汇报论证稿2019年4月）》。（嗯，我知道过时了）

基础资料怎么获取？

◆ ◆ ◆

其实说到底，双评价这个东西，就要依靠大量的空间资料，在没有资料的情况下，几乎不可能做出什么，而且部分资料属于涉密。我也只能学习下技术方法，对于基础资料，尽量利用网络上能找到的一些资源，找不到相关资料的，找另一些可以某种程度上替代的资料，或者找一些成果图件空间化作为分析基础。当然这样精度会有问题，不过作为练习也不会有什么问题。每一期文章我也会把用到的资源分享给大家。必须要注意：分享的资料只针对于练习，本身没有任何价值，我也不承担任何责任。（下次吧，这次的源文件找不到了，==嘻嘻嘻）

最后强调，不是教程不是教程，只是学习笔记，欢迎指正交流。

第〇篇 前期准备
1城市选取

作为练习城市的选择，还是主要考虑数据的可获取，以及有没有相关基础。其中相关基础是指手头是否能找到已经做的双评价成果可以参考。最好这个城市还要有一些复杂的自然条件，这样也可以多一些练习的机会。最后，我要比较熟悉这个城市，这样至少可以从直观的角度来看待练习成果。

综上，最后我选择长三角某县级市作为练习城市。该城市地处江南水乡，毗邻长江，境内地势低平、水网交织，地处温带，属亚热带季风性海洋气候。该市经济社会发展水平较高，历史文化资源较为丰富，具备一定的练习价值。

2确定坐标

根据技术指南：

评价工作运行环境采用地理信息系统软件，统一采用2000国家大地坐标系统，高斯-克吕格投影，陆域部分采用1985 国家高程基准，海域部分采用理论深度基准面高程基准。

首先网上不可知来源的数据一般都是WGS84坐标，要转成2000国家大地坐标系统，Z坐标选择85基准。

==我简直就是太天真了，以为这样就能成功了，结果开始的第一步就陷入了痛苦。

首先网上不可知来源的数据一般都是WGS84坐标，要转成2000国家大地坐标系统，Z坐标选择85基准。讲真，这特么的要怎么转换？？？？还要参数的嘛！！！（这一步先pass吧）

不管了==先凑合着用吧，应该不影响后面的吧，等知道怎么转换了再说吧（真是一个敷衍的人）。所以以下的所有地理坐标暂时都统一成WGS84吧。

第一篇 资源环境承载能力评价

1土地资源评价

1.1 高程、坡度分析

土地资源评价别看说的很复杂，说白了就是传统的高程分析和坡度分析，再加上土壤质地分类（地形复杂地区考虑地表起伏度）进行土地资源评价。一个个来说。

高程分级：

高程数据我们可以用地理空间数据云的30M空间分辨率DEM来替代。按照当地高程最大值最小值的情况进行合理分级，并进行“重分类——栅格转面——汇总”，可以得到不同高程区间土地面积及所占比例，同时将栅格转面后的面要素与区划数据进行“相交”地理处理，可以得到不同区域各类高程区间面积与比例。

坡度分级：

坡度的生成，可以很简单的用ArcGIS中的坡度工具得到。这里要注意的是，数据框需要是投影坐标，不然用坡度工具得到的结果会很奇怪（至少我这边是这样）。技术指南中的表述：

基于数字地形图，计算栅格单元的坡度，按≤2°、2~6°、6~15°、15~25°、>25°生成坡度分级图。

我们就用这个标准来分级。

就分级结果而言，某市地势平坦，全市85%的区域坡度在2°以下，平地面积极大，十分适宜土地开发建设。坡度超过15°区域仅有0.28%，主要集中于Y林场的丘陵低山区域，且坡度分布格局与高程总体大致相符。

土地资源评价与分级：

按照技术指南，土地资源评价与分级以坡度分级结果为基础，结合高程和土壤质地，将土地资源的可利用程度划分为高、较高、中等、较低、低5 种类型。土壤数据我暂时没找到适合的替代数据，暂时先不考虑，先把高程作为参考量。技术指南中的描述比较拗口：

高程≥4900m区域，农业和城镇土地资源等级直接取最低等级；高程在3000~4900m之间的，将坡度分级降2级作为农业土地资源等级、将坡度分级降1级作为城镇土地资源等级；高程在2000~3000m之间的，将坡度分级降1级作为农业土地资源等级。

我列了个表不知道对不对：

土地资源评价与分级的高程对应关系表

由于研究的某市高程均在1000m以下，土地资源评价与分级直接以坡度分级为基础。

地形起伏度：

地形起伏度的计算同样基于数字高程，利用“焦点统计”工具，由于我使用的高程数据是30m*30m精度，因此邻域采用15*15，统计类型要选择Range（计算最大值与最小值之差）。

地形起伏度的结果用于修正基于坡度的土地资源分级，对应关系如下：

土地资源评价与分级的地形起伏度对应关系表

湿地评价：

湿地评价这一项，在技术指南的土地资源评价中未列出，但是对于具体的研究对象来说，位于江南水乡，河网密布，湿地的影响巨大，并且，某同一区域的城市双评价工作中，也把这一项单独放在这里，所以我这边也一起考虑。湿地结果的修正，提取河流、湖泊、水库水面，对土地资源分级初步结果进行修正。将河湖水面区域土地资源评价结果直接降为最低等级。

河湖水系的基础数据我利用GLC30的2010年地表覆盖数据进行提取。（之前有一篇公众号文章写过这个数据源）

土地资源综合评价：

由于高程都处于同一对应关系的大区间内，因此综合评价时可以不考虑，也就是说，以坡度分级为基础，综合考虑地形起伏与江河湖海进行修正。

具体操作上，考虑到存在多区间的修正，本来想的是用栅格镶嵌直接取最小值的方法不可行了（为什么取最小值，因为土地资源评价以较低结果为最终结果）。最后思考的方法是栅格转面后，使用“相交”工具进行切割，生成含多个属性字段的要素。

第一步：分别对坡度基础的土地评价、地形起伏评价和湿地评价进行栅格转面。

第二步：三个要素图层相交。

第三步：地形起伏修正，选择100~200m的要素后字段计算“坡度分级结果减1”；选择≥200m的要素后字段计算“坡度分级结果减1”；湿地修正，选中江河湖海水系要素，将选中要素的坡度分级结果赋值为1。（1.2.3.4.5对应坡度基础上的土地资源评价分级的低、较低、中等、较高和高）

第四步：结果修正，选中坡度分级结果<1的要素，赋值为1（由于第三步，因此可能存在小于1的情况）