如何使用NetTopologySuite合并多边形？解决孔洞引发的合并错误

嘿，这个问题我之前在处理NTS多边形合并时踩过好几次坑！那个side location conflict错误本质上是拓扑结构不合法导致的——要么是顶点精度误差引发的边重叠/交叉异常，要么是小孔洞和主多边形的拓扑关系不符合OGC规则，还有可能是多边形自身存在自相交、环方向错误这类问题。下面我给你拆解一下检测和解决的具体步骤：

一、先定位问题根源：检测拓扑错误 #

1. 用NTS自带工具做有效性检查 #

直接调用几何对象的IsValid方法快速判断是否合法，如果返回false，用GeometryValidator获取具体错误细节，比如自相交、无效环方向、重复顶点等：

var validator = new GeometryValidator();
var errors = validator.Validate(yourPolygon);
if (errors.Any())
{
    // 打印错误信息定位问题
    Console.WriteLine($"Invalid geometry issues: {string.Join(", ", errors)}");
}

2. 可视化排查关键坐标点 #

错误日志里给出的坐标(12.235768128349596, 44.387249570104892)是重点！把多边形导入QGIS这类GIS工具，放大到这个位置，看看是不是有：

• 极小的重叠边（因为浮点精度差导致看似重合实际不重合）
• 小孔洞的边和主多边形边几乎贴合
• 自相交的线段

3. 检查环的方向是否合规 #

NTS要求多边形的外环是逆时针方向，内环（孔洞）是顺时针方向，如果方向反了，合并时很容易触发拓扑冲突。可以用以下代码检查并修正：

// 检查外环方向，逆时针为合法
if (!polygon.Shell.IsCCW)
{
    polygon.Shell = (LinearRing)polygon.Shell.Reverse();
}

// 检查所有孔洞方向，顺时针为合法
foreach (var hole in polygon.Holes)
{
    if (hole.IsCCW)
    {
        hole = (LinearRing)hole.Reverse();
    }
}

二、针对性解决问题 #

1. 处理浮点精度问题（最常见原因） #

大部分拓扑冲突都是因为坐标的小数位数太多，导致计算时出现精度误差。你可以：

• 减少坐标精度：用PrecisionReduce方法把坐标保留到合适的小数位数（比如6位，对应米级精度的厘米级别）：

  var cleanedGeometry = yourPolygon.PrecisionReduce(6);
  

• 使用Snap机制合并：不要直接用Geometry.Union()，改用SnapIfNeededOverlayOp并设置合适的snap距离，让接近的顶点自动对齐：

  var snapTolerance = 0.0001; // 根据你的坐标单位调整，比如米的话设0.1毫米
  var unionResult = SnapIfNeededOverlayOp.Union(polygonA, polygonB, snapTolerance);
  

2. 修复无效几何 #

NTS有个“黑魔法”操作——Buffer(0)，它能自动修复绝大多数拓扑错误（自相交、环方向错误、重复顶点等），还能顺便移除极小的孔洞：

var fixedGeometry = yourPolygon.Buffer(0);

注意：这个操作可能会改变极小的几何细节，所以要根据你的业务场景调整，比如如果是高精度地图数据，先测试再用。

3. 手动清理小孔洞 #

如果明确是小孔洞导致的问题，可以先遍历多边形的内环，过滤掉面积小于阈值的孔洞，再进行合并：

double minHoleArea = 0.001; // 阈值根据你的坐标单位设置，比如平方米
var validHoles = yourPolygon.Holes
    .Where(hole => hole.Area > minHoleArea)
    .ToArray();
var cleanedPolygon = new Polygon(yourPolygon.Shell, validHoles);

4. 完整的安全合并示例 #

把上面的步骤整合起来，写一个健壮的合并方法：

using NetTopologySuite.Geometries;
using NetTopologySuite.Operation.Overlay.Snap;
using NetTopologySuite.Operation.Valid;

public Geometry MergePolygonsSafely(IEnumerable<Geometry> polygons, double snapTolerance = 0.0001, double minHoleArea = 0.001)
{
    // 先逐个清理每个多边形
    var cleanedGeometries = polygons.Select(g => CleanGeometry(g, snapTolerance, minHoleArea)).ToList();
    
    // 用SnapIfNeededOverlayOp逐步合并
    var unionResult = cleanedGeometries.Aggregate((g1, g2) => SnapIfNeededOverlayOp.Union(g1, g2, snapTolerance));
    
    return unionResult;
}

private Geometry CleanGeometry(Geometry geometry, double snapTolerance, double minHoleArea)
{
    // 修复无效几何
    if (!geometry.IsValid)
    {
        geometry = geometry.Buffer(0);
        // 再次检查，如果还是无效就抛出异常（或者做更细致处理）
        var validator = new GeometryValidator();
        var errors = validator.Validate(geometry);
        if (errors.Any())
        {
            throw new InvalidOperationException($"无法修复的无效几何: {string.Join(", ", errors)}");
        }
    }
    
    // 处理多边形的孔洞
    if (geometry is Polygon polygon)
    {
        var validHoles = polygon.Holes
            .Where(hole => hole.Area > minHoleArea)
            .Select(hole => (LinearRing)hole.PrecisionReduce(6))
            .ToArray();
        return new Polygon((LinearRing)polygon.Shell.PrecisionReduce(6), validHoles);
    }
    
    // 减少其他几何类型的精度
    return geometry.PrecisionReduce(6);
}

内容的提问来源于stack exchange，提问作者Luca Morelli