在公差范围内溶解QGIS（GEOS）中的多边形

我有一个多边形矢量层，代表洪水淹没的区域。该层的制作方式使得某些区域由多个较小的接触多边形组成（请参见下面的屏幕快照，标记为“原始数据”）。

我想溶解它们接触的多边形。如果我在使用GEOS库的QGIS fTools中使用了溶解地理处理工具，则会得到一些意外的结果（请参阅“使用QGIS fTools溶解”）。似乎在这些问题中提到，GEOS并没有使用阈值来确定多边形是否在接触：

http://trac.osgeo.org/qgis/ticket/3126

http://trac.osgeo.org/qgis/ticket/3184

http://hub.qgis.org/issues/2806

在ArcGIS中，相同的方法可以提供预期的结果（请参阅“使用ArcGIS 10.1分解”）。

显然，这是GEOS的“功能”，而不是错误。有什么办法可以修复我必须产生预期结果的数据？对我来说，轻微的精度损失（几厘米？）是可以接受的。理想情况下，解决方案是可以将其写入QGIS插件中，作为预处理步骤。

相关问题（无答案）：如何在QGIS中设置地理处理工具的公差值？

以shapefile链接到示例数据：https : //www.dropbox.com/s/c8ns7wwfv289ad0/sample.zip

在这种情况下，您似乎需要通过首先修复拓扑来帮助溶解工具。这是您在GRASS中执行此操作的方法。

加载shapefile时，GRASS会看到它：

每个十字标记一个区域。您会看到边界线上有一些十字：这些是由边界线之间的微小偏移引起的条形多边形。

1. 删除由真实多边形之间的间隙形成的条形多边形：

   v.clean input=sample output=tempA --o tool=rmarea thresh=0.001

2. 捕捉边界线：

   v.clean input=tempA output=tempB --o tool=snap thresh=0.001

3. 在新的相交处断开边界线：

   v.clean input=tempB output=tempA --o tool=break thresh=0.001

4. 删除现在重复的边界线：

   v.clean input=tempA output=tempB --o tool=rmdupl thresh=0.001

5. 溶解多边形

   v.dissolve input=tempB output=cleaned --o column=ID

等等：

使用GRASS GIS命令v.dissolve（在Sextante插件-> GRASS命令->矢量中）

如果我快速重现您的示例：

v.dissolve的结果：

说明

这个问题很容易理解。shapefile没有拓扑，如果两个区域共享一个公共边界，则该边界将被数字化两次，并且也将重复存储。

当您尝试编辑shapefile时，您会清楚地看到以下内容：

GRASS GIS（作为ArcGIS）是一种拓扑gis。在GRASS中导入shapefile时，结果层具有拓扑结构：两个区域之间的公共边界存在一次并被共享。

因此，当您要分解shapefile（分解两个边界）时，可能会遇到问题（在QGIS，GvSIG，OpenJUMP或旧的ArcView 3.x中），请参阅ArcUserI：Understanding Topology和Shapefiles。

在GRASS中，v.dissolve旨在消除具有相同类别编号的区域之间的公共边界