[发明专利]基于GIS的区域网格化人口密度空间化计算方法

摘要

本发明公开了一种基于GIS的区域网格化人口密度空间化计算方法，充分考虑影响人口密度的自然因素和社会经济因素的区域分布特征，以区域的最小行政单元村为单位建立统计模型，对研究区域进行100m*100m的网格划分，将村级统计模型反演到每一个网格中，以格网单元来描述人口分布；本发明所设计的基于GIS的区域网格化人口密度空间化计算方法能够更好的提高人口空间分析精度。

主权项

一种基于GIS的区域网格化人口密度空间化计算方法，其特征在于，包括如下步骤：步骤（1）：任意选取一个研究区域，统计该研究区域的行政村级人口数据，得出该研究区域每个行政村级人口数据Mn，n=1,2,3...m，其中，m定义为该研究区域中行政村的总数；步骤（2）：通过相关分析法得出与研究区域人口密度的影响因子Ai，i=1,2,3...q，其中，q定义为影响因子的个数；步骤（3）：建立人口数密度Ps和影响因子Ai的回归方程： Ps = Σ i = 1 q C i A i + μ 其中，Ci定义为回归系数，μ定义为误差；步骤（4）：选出类型为面数据集的影响因子，将面数据集影响因子的图层逐个与行政村图层B进行求交，得出每个行政村中图层中面数据集影响因子的面积areaα，α=1,2,3...s，其中，s定义为类型为面数据集的影响因子个数；步骤（5）：选择类型为点或线数据集的影响因子，将点或线数据集影响因子以半径ka作多态缓冲区，a=1,2,3...g，g定义为类型为点或线数据集的影响因子的个数，对多态缓冲区赋予权重Wa，将多态缓冲区图层逐个与行政村图层B进行求交，利用面积加权法计算点或线数据集影响因子在每个行政村中的面积areaβ=areaβ1×W1+areaβ2×W2...+areaβg×Wg其中，β=1,2,3...m，areaβ1,areaβ2,areaβ3,...areaβg表示在第β个行政村图层中分别以k1,k2,k3...kg为半径形成的多态缓冲区的面积；步骤（6）：经过步骤（4）和步骤（5）的处理后，共形成q个影响因子图层属性表，每个属性表中都有一个area字段属性值，将q个图层属性表中的area字段属性值都放到行政村图层B的属性表中；步骤（7）：将研究区域划分为S个区；步骤（8）：对行政村图层B的属性表中的area字段进行标准化，将标准化后的值代入步骤（3）中的回归方程，列出S个回归方程，解出系数Ci；步骤（9）：对研究区域进行网格划分，每个网格单元的划分标准为100m×100m，共划分为t个网格单元；步骤（10）：将影响因子的图层与网格图层进行求交分析，当影响因子图层为面要素时，直接将影响因子图层与网格面图层求交得出每个网格单元中影响因子的面积areaω，ω=1,2,3...s；当影响因子图层为点或线要素时，对点或线以半径rb作多态缓冲区，b=1,2,3...g，对多态缓冲区赋予权重Vb，将多态缓冲区图层逐个与网格图层进行求交，利用面积加权法计算影响因子在每个网格图层中的面积：areaθ=areaθ1×V1+areaθ2×V2...+areaθg×Vg，并将上述所有的area字段全部放到网格图层属性表中，其中，θ=1,2,3...t，areaθ1,areaθ2,areaθ3,...areaθg表示在第θ个网格图层中分别以r1,r2,r3...rg为半径形成的多态缓冲区的面积；步骤（11）：将网格图层属性表中的area字段进行标准化，将每个网格中的面积值代人步骤（8）中相应的回归方程中，得到每一网格单元中的人口数，对研究区域中的所有网格单元进行统计，获得该研究区域的人口密度。