目录

上一主题

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。