计算鼓楼区人口密度

2019-08-01 本文已影响0人机智的柠檬

1、读取含有POI点的CSV文件，将点在Arcmap中进行展示。

image.png

该数据集中，共有205567个POI，分为8类，裁剪出鼓楼区的POI
ArcToolBox--分析工具--提取分析--裁剪
输入数据：all_point
裁剪数据：鼓楼区矢量图
裁剪结果：

image.png

2、分别提取出鼓楼区各个类别的POI

image.png

将结果导出：

image.png

导出8类数据，分别加载已导出数据

3、找出每个格网中包含的POI 计算每个格网的中心POI点

加载格网数据，利用之前文章中算出的POI类别与人口密度的相关系数将格网人口按比例分配给每一类

image.png

计算公式为：a类人口=格网人口数*系数a/总系数
对每个格网进行迭代，找出每个格网中所包含的POI 并将其信息输出成JSON格式
具体模型为：

image.png
输出的JSON名称为格网编号，代表每个格网中的POI信息。
对每个POI经纬度取其平均值，计算出每个格网的中心POI

并输出为CSV文件
具体代码如下：

import json 
import csv
list1 = []
for i in range(1,5548):
    list2 = []
    filename = 'E:\\人口密度\\dian\\gezidian\\' + str(i)+'.json'
    with open(filename,'r') as f:
        js = json.load(f)
    features = js['features']
    print(i)
    lng = 0
    lat = 0
    count = len(features)
    for l in features:
        lng += l['attributes']['Field5']
        lat += l['attributes']['Field6']
    if count != 0:
        lng_ave = lng/count
        lat_ave = lat/count
        list2.append(i)
        list2.append(lng_ave)
        list2.append(lat_ave)
        list1.append(list2)
    else:
        lng_ave = -1
        lat_ave = -1
with open('E:\\人口密度\\0731\\zb.csv','w',encoding='utf-8',newline='') as f:
    csv_file = csv.writer(f)
    for l in list1:
        csv_file.writerow(l)