Qgs开发-拓扑分析
sf2gis@163.com
2014年10月22日
1 概述
拓扑分析一般是指地理元素之间的相关关系。一般是以基本的空间关系为基础,进行综合运算,计算元素之间的复杂关系,如空洞检查、重叠检查等。
参考:http://en.wikipedia.org/wiki/Geospatial_topology
http://en.wikipedia.org/wiki/DE-9IM
http://help.arcgis.com/zh-cn/arcgisdesktop/10.0/help/index.html#//006200000003000000
http://docs.qgis.org/2.0/en/docs/user_manual/plugins/plugins_topology_checker.html
2 原理
针对特定的几何关系,综合运行基本空间关系,检测出特定的目标元素。
3 方法
3.1 空洞检测(gaps):检测一个多边形是否存在空洞
1) 首先获取多边形geo的外接多边形,并做一个缓冲区buff,目的是最外层的多边形是一个完整的多边形。
2) 然后将这个buff与geo做difference运算,获取在buff中,而不在geo中的多边形,此时会生成最外层的一个多边形(第一部中如果不做buff,这里会有很多的外部多边形,无法与内容空洞区别)和内部的空洞。
3) 删除最外层的多边层,剩余的多边形就是内部的空洞。
3.2 边界融合:将一个多边形的边界合并到另一个多边的公共边上,删除相交部分,填补空洞
1) 首先删除相交部分:geo2与geo1做difference,从geo2中删除与geo1相交的部分。
2) 然后将geo2与geo1合并(combine),检测合并后多边形的空洞(参见空洞检测(gaps):检测一个多边形是否存在空洞)。
3) 将空洞与第一步中的diff合并:将所有空洞添加到geo2中。
4 示例
4.1 空洞检测detectGaps
/**
*@briefQgsVectorAnalysis::detectGaps
*
*getgeometrygaps.
*@notefromtopoTest.cpp
*@parampGeo
*@returngapsgeometry
*@authorsf2gis@163.com
*@date2014-10-2118:31:50
*/
QgsGeometry*QgsVectorAnalysis::detectGaps(QgsGeometry