【双评价学习笔记】城镇适宜性—气候评价

写在最前面：数读菌更新是真的很慢，请大家见谅啊。毕竟我也要工作，每天挣扎在贫困线上下，吃了上顿没下顿，也不能全心全意写公众号。写公众号也没有收益，全然靠兴趣坚持，所以大家给我一点点包容吧！不包容的话！那也没办法！哼！

城镇建设适宜性评价呢，目前更新了土地资源评价和水资源评价。今天写的是气候评价。主要用到的数据是气候数据集（废话），具体来说就是温度与湿度的数据。数据分享在文末。顺道说一下啊，前面几篇文章中的数据资料已经开放无条件分享了，获取方式可以看前面的文章，关键词没变，但是不需要条件了。

废话不多说了吧，毕竟打字也挺累的，直接开始吧。

----------------------------

第二篇  城镇建设适宜性评价

1

单项评价

1.1 土地资源评价

1.2 水资源评价

1.3 气候评价

说到气候评价，之前农业适宜性那边也有，具体的方法差不了多少，无非就是处理两个数据：温度和湿度。然后用指南里提到的公式算一下每个月的舒适度。

温度数据的处理。数读菌分享的数据是2010-2017年的日值温度数据。之前有小伙伴私信问我啥是日值。大概就是说每年每月每天的数据吧。数据量有点大，但是还好不算比较麻烦的处理。数据是按年月分割的txt文件，稍加处理导入arcgis，根据经纬度生成点就行。

2010-2017日值温度数据

气候评价的依据“舒适度”计算的时候要用到的是每个月的平均温度，因此，在把数据导入arcgis中之后，需要计算每个月的平均温度。计算方法很简单，使用【融合】工具把生成的点要素融合，在“融合字段”复选框中选中“区站号”和“月”（这里是我自己的字段命名，懂了就好，具体自己看着办），在“统计字段”复选框中选择“平均温度”（这里的平均温度是指日平均温度，也是我自己的字段名，分享的数据是没有字段名称的，需要在处理数据的时候自己添加，按照说明文档就可以）。简单说一下，这一步骤的操作是为了得到每一个气象站每个月的平均温度。“融合”工具不止是可以把拥有相同属性值的要素合并成一个，也可以用来统计（请划重点）。最后生成的新的要素属性表如下：

融合后的温度数据

融合后温度数据（要素）就可以用插值的方法生成每个月的平均温度图了。插值我不用多说了吧。但是这里有一个坑数读菌踩了一下，所以提醒下大家哦，指南里面的公式用的月均温度是华氏温度，所以最好在属性表里“计算字段”一下，把原有的摄氏度调整成华氏度之后再进行插值。数读菌一开始没有注意到，所以是在插值后利用“栅格计算器”调整的。

1月-12月月平均温度图（过程文件未调整分级）

哦对了，温度插值需要结合海拔校正，具体的方法在之前农业适宜性的气候评价部分写过了就不再赘述，具体可以自己去看之前的文章，放个传送门吧：

海拔校正后的1月-12月月平均温度图（简单分级）

湿度数据的处理。数读菌分享的数据是2010-2017年的日值湿度数据。湿度数据的处理没有温度那般复杂，少了几个步骤，最开始的操作却也是一样的。融合数据后进行插值：

1月-12月月平均湿度图（过程文件未调整分级）

舒适度计算。舒适度采用温湿指数表征，计算公式为：

式中：THI 为温湿指数，T 为月均温度(华氏温度)，f 是月均空气相对湿度(%)。根据之前得到的每月平均温度和每月平均空气相对湿度，利用“栅格计算器”可得到每月温湿指数：

1月-12月温湿指数图（简单分级）

按照舒适度分级标准对每个月温湿指数进行分级：

舒适度分级标准

1月-12月舒适度分级（过程数据色彩未统一）

区域舒适度评价。指南中说“取12个月舒适度等级的众数作为该区舒适度”。一开始数读菌看到众数，满心欢喜的以为用“像元统计数据”这个工具就可以得到最后的结果，可是最后生成的结果有很多NoData，查看帮助文档才发现问题出在哪里：

本以为可以直接计算众数的“像元统计”

像元统计中，众数的统计面临这么一个问题：“如果一个像元位置有多个众数值，则该位置的输出为 NoData。 ”这也就是一开始直接计算的时候有很多区域是空值的原因。

“像元统计”工作原理之众数计算

说实话，数读菌对此也没想到什么好办法。我采取的办法在我看来，肯定可以被别的更合理的办法取代。那先看一下我的方法。首先是用“等于频数”工具计算“基于单个像元计算一组栅格数据值等于另一个栅格数据值的次数”。大家有更好的方法还望不吝赐教。

“等于频数”工作原理

这一步的价值就是求得各区域舒适度分别为1、2、3、4、5、6、7的时候出现的频数。将生成的结果转为要素数据后添加字段，字段名为当前舒适度等级，取值为当前频数，然后进行“相交”操作。相交生成的属性表添加一个取名为“值”的字段。最后生成的属性表如下：

这一属性表就可以用来统计众数，由于一共也就7个字段比较，其实手操也不麻烦，比如图中选中的面代表着这个面舒适度为1的时候出现了8次，那这个面最后的值为1。而当不同舒适度频数相同时，我一般取的值以等级低的优先。最后统计完成后以“值”字段为value值转回栅格即可：

气候评价

本次分享的数据为温度日值和空气相对湿度日值。2010-2017年。下面是属性字段说明：

------------------------

嘛--今天就先到这里了~~~

温度和湿度数据后台回复“双评价5”

往期回顾（往期的数据的无条件获取均已开启哦，请点进去看关键词吧）：

土地资源评价（分享土壤质地数据、地形高程数据和研究边界数据）戳下面：

水资源评价、气候评价（分享长时间年值、近几年日值气候数据）戳下面：

环境评价、灾害评价中的干旱部分（分享1950年以来的日值降水、平均气温、最高气温、最低气温数据）戳下面：

灾害评价中的雨涝、高温、大风灾害部分（2010年-2017年的全国风速风向数据）戳下面：

生态评价戳下面：

农业集成评价戳下面：

T_T对上面文章的推倒重来戳下面（分享30m地形高程数据）：

城镇建设适宜性评价从这一篇开始，首先是城镇土地资源评价：

城镇适宜性评价中的水资源评价戳下面：

<END>

整理数据资料不易

喜欢这篇文章的话

记得帮我们分享点赞

特别特别喜欢的话

也可以赞赏啦

如需转载文章、投稿或者合作

请直接在公众号后台联系

真的特别欢迎投稿啊！！！！

使用关键词搜索历史文章请点击【阅读原文】