面积权重法就是根据每个地区每种土地利用分类对应的一个产业GDP的权重,来计算每个1km×1km的格网内各土地利用类型所对应的GDP值,然后将各产业GDP进行相加,便得到一个格网内的总GDP数据。
数据来源
主要包括分辨率为30m的土地利用分类数据(中科院科学数据库)、国内生产总值数据来源武汉统计局的统计年鉴以及知网中国经济社会大数据研究平台(湖北省经济社会发展统计数据库),以及武汉市各区行政矢量图(公众号上找到的)。
方法
本次实验按照面积权重法来进行GDP空间可视化。分产业进行,第一产业模型同文献[1],只是将其模型中的gij用GDP平均值来代替(统计各区GDP/统计各区各土地分类面积),与其相关的土地利用类型一级分类包括耕地、林地、草地以及水域。二三产业合起来一并进行可视化,模型同第一产业,与其相关的土地类型包括城乡工矿居民用地。 选择这样的原因主要是因为原则上只要统计的各分类面积值不出错的话,模型计算得到的结果和统计值几乎不会出现太大的偏差,能满足空间可视化的要求。
具体步骤
矢量图处理
由于网上找不到武汉市的带有各区的区名属性的完整市行政区矢量图,因此只能利用在公众号上找到的各区矢量图,自己在ArcMap中把十三个矢量图拼成一幅图,并添加字段“区名称”。(具体步骤可以百度)
空间化步骤
**1.**土地利用分类TIFF图的重分类:将一级分类的耕地、林地、水域和草地合并成一类A,城乡工矿居民用地合并成一类B。
**2.**利用各区的行政矢量图裁剪土地分类TIFF图,然后分别统计各区各土地分类类别的面积(属性表中,像元个数×900),单位是平方米。
**3.**结合统计好的各区的分产业GDP,分别统计农林牧渔产值和剩余产值。分别对应A类和B类,并在EXCEL中计算各区各类别的GDP权重。
**4.**将TIFF格式的土地利用图转化为矢量图、创建渔网(1km×1Km)、将土地利用矢量图+渔网+武汉市行政区图三者进行“相交”处理(!一定要注意统一投影坐标系,并将结果保存至文件夹中,否则无法“计算几何”求面积)
**5.**在相交成果图的属性表中添加权重字段,“按要素选择”各区输入之前计算的各权重值,然后将A类数据导出为一个新的dbase表,B类同理。然后分别对两个表的对“Shape_area”字段进行“计算几何”,选择“面积”,单位是“平方米”。然后添加字段“分产业GDP”,利用“字段计算器”,权重×面积即为分产业的GDP值。
**6.**分别计算完两类各自对应的GDP值后,将两个表都连接到相交成果的属性表上(!这时会出现“空”的值,可以将表导出,用WPS打开再保存,再导入到ArcMap中即可),然后添加字段“GDP”,利用字段计算器将两类计算的各产业生产总值进行相加即可。
**7.**最后利用符号系统显示即可得到成果图。
总结
①方法模型还是太简单了,只是操作繁琐了些。
②很多文献的方法难以直接通过阅读文献就复现,而且更难以达到文献中所模拟出的值的数据精度。
③统计值与模拟值在千元内相关性为1。精度很高。(第一产业生产总值差值相比较而言较大,因为使用的武汉市行政图是最新的,与各年份的土地利用图边缘并不完全重合,而边缘的几个区农牧渔业较发达,因此对第一产业生产总值影响较大)。
参考资料
[1] 刘红辉, 江东, 杨小唤,等. 遥感支持下的全国1KM格网GDP的空间化表达[C]// 全国地图学与GIS学术会议. 0.
[2] 黄莹,包安明,陈曦,刘海隆,杨光华.新疆天山北坡干旱区GDP时空模拟[J].地理科学进展,2009,28(04):494-502.