人口空间化（一） 面积权重法实现人口空间化

大约两年前我从我老师那拿到这个方向的时候，全程一脸懵逼，然后我在CSDN搜到一版在ArcGIS中利用“行政单元面积权重法”实现人口数据格网化，我用这个方法做了我第一版人口空间化结果。

所谓的“面积权重法就是根据你选的要素占格网的面积来分配每个格网上的人口数量，比如你选择土地利用数据中的建设用地作为人口分布指标（也就是说你认为人口都分布在建设用地上）。具体过程如下（以北京作为研究区，土地利用数据作为建模因素）：

1、相关性分析，确定土地利用数据中哪几类地表与人口分布相关性较大（SPSS可以做）；

2、确定好建模因素后，比如选择的是居住区、工矿区、草地，按属性提取将这几类地表从土地利用数据中提取出来；

3、统计各区总人口和这三类地表的面积，在SPSS中做多元回归分析，确定每类地表覆盖类型的建模系数，其中统计总人口：spaital analyst->区域分析——按表格显示分区统计——统计类型选sum；统计三类地表面积：spaital analyst->区域分析—>面积制表

                                               

类似这种形式，但你拟合的不一定是个多元线性的，这个看具体情况；

 

4、生成渔网，如果你做的分辨率比较高，这里生成渔网会很麻烦，而且渔网太大的话无法生成，经常说你跑了3个小时（一点也不夸张，我曾经一晚上都没弄出来过）结果告诉你生成不了，新建一个字段pop_density

那应该怎么做的：1、尽量避开这玩意，你明确你的目的是啥，其实不用生成渔网就做的方法很多；

                             2、如果一定要做，那给你点忠告，把他保存在数据库（.gdb）中，这样虽不会加快速度，但最后一定会生成；

                             3、你要生成的是面状渔网，记得生成样本点一定要勾上，不然做出来也没法用

5、统计每个格网上三类地表的面积（方法同3）

6、将你生成的渔网和北京行政区划联合（主菜单栏—地理处理—union（联合））

7、用字段计算器计算pop_density，将这个字段转为栅格

8、最后由于模型误差，你模拟的总人口和实际可能有差，这个时候用下面这个公式校正一下即可

       校正后人口=人口普查总人口/模拟总人口*7生成的栅格人口；

备注：这个方法容易出问题的地方：

1、格网生成太麻烦

2、第三步拟合出来的结果很有可能有些系数为负数，这时候你人口拟合出来的就是负值，不符合实际物理意义。