处理人口栅格数据学习笔记

我这次做的人口栅格数据是用的Landscan的数据（ORNL LandScan Viewer - Oak Ridge National Laboratory）可以直接根据需求下载相应的年份，但是还是跟统计年鉴当中的人口数据有些出入，所以想用统计年鉴当中的人口数据修正下Landscan的数据。Landscan的数据我不知道是不是下载下来可以直接使用，因为我国的人口统计数据通常以行政区为单元逐级统计和汇总而来，所以本着严谨的科学态度修正下，看一下是否实际情况。

Landscan的数据介绍可以去其他人博客里搜一下。我就简单介绍下Landscan的数据：这是个全球的人口数据，分辨率是1km*1km，可以去上面连接直接下载。时间是2000年到现在，已经更新到2023年。

因为我的研究年份是从1992年开始，所以用了中科院的人口数据（资源科学数据中心中国人口空间分布https://www.resdc.cn/data.aspx?DATAID=251）注意：这个网站的人口数据只有这样的年份数据：1995、2000、2005、2010、2015和2019年6期，所以还是不大好用，如果研究年份刚好是这样的年份可以使用。

接下来我叙述下处理的步骤：

第一步，当然是下载数据，下载下来解压后是这样的数据，有tif文件可以直接在arcgis里提取自己想要的区域，提要是最好先建立个缓冲区，这样边界处就没有太大的缝隙

第二步，是找统计年鉴里的人口数据，最好是用个excle表格。因为几十年来有行政区的变化，人口数据有些不太统一，最好是根据你自己的感兴趣区域（矢量边界）里进行统一。这里1992年没有人口栅格数据，后来用的1995年，而青岛市城阳区的设立是在1994年，刚好绕过了这个问题，所以又换成了1995年的数据（中科院）

 

第三步，利用argis工具箱里的以表格显示分区统计，统计各个市、区的栅格值的和，统计类型选sum

统计好后可以跟统计年鉴中的数据做个比较，可以发现还是有些出入，但是数量级是在同一个数量级，个别区域会差的数量比较大，计算他的修正系数：统计年鉴数据/Landscan数据=修正系数

在矢量边界中添加个字段，将修正系数填进去，注意对应的区域，有一次没注意对应区域，怎么算怎么不对，最后检查了一遍才找到问题所在。添加字段的时候注意是浮点型，整型的话会变成0，后面转为栅格后的栅格数据也是浮点型栅格

第四步，就是用面转栅格把修正系数转为栅格，以便于后面用栅格计算器计算，如图所示

第五步，利用栅格计算器把修正系数和landscan数据相乘，就得到结果啦。在属性里可以看到得到的结果也是浮点型，可以后续在转为整型，转为整型的目的是可以打开属性表，打开属性表可以看一下最大的那个栅格是否复合实际情况。我这里是有一个非常大非常大的值，我把它当作异常值处理了，如果没有异常值的话到这一步应该就结束了，有异常值的话可以用栅格计算器处理一下把它赋值成其他的值。

第六步，利用栅格计算器里的条件函数处理异常值

con（“你的栅格”>123456789，**，“你的栅格”），**为你想把这个异常值赋成哪个值，这个异常值的确定我是用envi里的统计值确定的。最大的这个29491我感觉这个值有点异常，而且只有这一个像元是这个值，就赋为了其他的值。

至此，我认为我处理的这个人口栅格感觉还是比原来直接下载下来的数据正确一点的，起码从统计上的数据来讲比较符合，没有相差太大。

另，免责声明：我也是在学习阶段，如有同学按照我的方法出错，本人概不负责，如果有同学感觉这个方法不太靠谱，有更好的方法欢迎留言讨论。