点文件从矢量文件中抽取数据[spatial join]

时间 2021-06-11 标签 02地理信息系统导论学习

1.首先在矢量文件的属性表中计算人口密度。高亮矢量文件，右键点击【Open Attribute Table】打开属性表。新建字段’area1‘,数据类型设置为’double‘,点击确定。

创建一个字段“density1”（字段类型为double）。选择density1列，用filed calculator计算人口密度。"density1 = dp001001/area1"；

2.图层叠置操作。把生成的渔网图“grid20_1105”和有“idcensus”图层合并为一个图层“union”。形成了若干个细碎的小多边形。

3.计算"union"图层的每个细碎的小多边形的人口数目。(因为渔网图“grid20_1105”是规则的矩形，“idcensus”图层是不规则的多边形，两个图层叠置之后生成的面是细碎的多边形，只有“idcensus”图层是有人口密度的属性的，而叠置之后的多边形需要重新计算人口密度，在union.shp中创建一个字段”area2“（double类型），计算多边形的面积（平方公里）在union.shp中创建一个字段”pop“（double类型），计算多边形的人口数（pop = area2 * density1）)

4.计算渔网中每个方形网格的人口数（spatial join）。target为grdi20.shap
join feature为union.shp
选择 Join one-to-one方式
输出图层的属性：一定把union.shp的pop加入，其它可以删除（保留也可以）。注意：一定要高亮pop，右击鼠标键弹出快捷菜单，选择merge rules中的sum。
sum表示一个方形格的人口数是由其对应的几个细小多边形的人口数的汇总而成。
match选项：一定要选contains，表示target的方形多边形(grid20_1105.shp)和join-feature的不规则多边形（union.shp）的空间关系。这里将一个方形多边形同它内含的几个细小多边形建立了对应关系。

5.转为raster文件。【polygon to raster】，属性字段选pop作为栅格的像元值,像元大小为20000米，同时注意设置environment环境变量，如：output coorodinate, processing extent.

6.采样。通过search命令，可以搜索到”Extract Multi Values to Points“命令。位于toolbox空间分析的extract命令中。将加载矢量的样本点samples.shp文件（请自行用create random points生成）和人口栅格文件，确定。检查samples.shp，会发现增加了一个新列popgrd，内含样本点在对应popgrd像元上的数值。