目录

上一个主题

9.1. 缓冲区分析

下一个主题

9.3. QGIS交集分析(Intersect)


本教程由OSGeo中国中心维护

9.2. 分割分析

9.2.1. 什么是分割

区域创建是GIS的主要功能之一,除了缓冲区创建外,还包括分割。分割是一种设定影响范围的方法,考虑到多个物体的影响, 也称为泰森分裂。泰森分裂是将区域划分为由垂直平分线形成的多边形的方法,假设连接存在多个设施的区域中的任何两个 相邻设施的线段。

在本节中,将就利用点数据进行泰森分裂的解说。这里使用的数据是“全球地震数据.shp”,使用在缓冲区制作的操作解说中 使用的文件夹的“空间分析.qgs”中嵌入的层。此处创建的所有文件都存储在“空间分析”中。

9.2.2. 通过分割破碎来生成区域

首先,从QGIS主菜单中选择 Vector ‣ Geometry Tools ‣ Voronoi Polygons , 调出 "Voronoi Polygons" 窗口。在该窗口的选项卡菜单中选择 Parameters , 然后在 "Input Layer" 中选择 "world earthquake(EPSG:4326)" 。 Buffer region 保持 "0.000000" 。 单击右侧的 "..." 按钮并选择 Save file 。显示 "Save file" 窗口时,在“保存类型”中选择 "SHP files(.shp)" , 另外,指定文件夹应保存为“地震点_分割.shp”,然后单击 Save 按钮。 选中 "Open output file after running algorithm" 并单击 Run 按钮,如 图 9.5

_images/img_1_x4b.png

图 9.5 Voronoi Polygons设置窗口

处理结束后,"Layers"上就会出现"Voronoi Polygons"层,描绘出破碎的分割后的多边形,如 图 9.6 。 “地震点_分割.shp”在“空间分析”文件夹中创建,并在 "Layers" 中显示为 "Voronoi Polygons" 。 要更改图层名称,请右键单击 Layers 中的 "Voronoi Polygons" ,然后从弹出菜单中选择 Rename Layer 。 在这种情况下,可以在图层面板中更改图层名称,因此将 "Voronoi Polygons" 更改为“地震点_分割”。 在这里,在 Layers 中点击“地震点_分割”层时,将其移动到低于 "world earthquake(EPSG:4326)" 的位置。 此外,取消选中其他图层,以便它们不会显示在地图视图中。

_images/img_2_xdy.png

图 9.6 Voronoi Polygons结果生成图

9.2.3. 计算面积和周长

如果制作了Voronoi Polygons,则要求各个区域的面积和外周线的长度(外周长)。 在QGIS主菜单中选择 Vector ‣ Geometry Tools ‣ Add Geometry Attributes , 调出 "Add Geometry Attributes" 窗口,然后在窗口的选项卡菜单中选择 Prameters 。 在 "Input Layer" 中选择 “地震点_分割(EPSG:4326)” , 并选择 "Layer CRS" 作为 Calculate using ,如 图 9.7

_images/img_3_xat.png

图 9.7 Add Geometry Attributes参数设置

在 "Added geom info" 中,按右侧的 ... 按钮并选择。选择 "Save to file" 出现 "Save file" 窗口时, 在“保存类型”中选择 "SHP files(.shp)" ,此外,指定文件夹应保存为“地震点_分割2.shp”, 然后单击 Save 按钮。选中 "Open output file after running algorithm" 并单击 Run 按钮。 处理完成后,在 "Layers" 中创建名为 "Added geom info" 的图层。

文件夹创建了“地震点_分割2.shp”,这在 "Layers" 中以 "Added geom info" 的名义被显示。 右键单击 Layers 的 "Added geom info" 以更改该层名称,然后在出现的菜单中 选择 Rename Layer 。这样,就可以在 "Layers" 中更改层名称, 因此将 "Added geom info" 改为“地震点_分割2”。

在此,右键单击 Layers 的层“地震点_分割”,在出现的菜单中 选择 Open Attribut Table 看属性表,只记录我们前面输入的项目。 接着,用同样的操作查看层“地震点_分割2”的属性表的话, 新追加了 "area_2" 和 "perimeter" 的项目。该 "area_2" 是每个区域的面积(单位:m2), "perimeter" 是周长(单位:m),如 图 9.8

_images/img_4_xq3.png

图 9.8 Geometry Attributes计算结果

完成目前的工作后,选择 Project ‣ Save ,然后保存“空间分析2.qgs”。之后, 从主菜单 Project ‣ Exit QGIS 退出QGIS。