如何合并并简化Java中表示街道几何的2D多边形列表？

我在Java中有一个表示街道几何的2D多边形列表，类型为List<double[]>，每个数组以坐标对（x1,y1,x2,y2,...）定义多边形，示例为{x1,y1,x2,y2,x3,y3,x1,y1}。希望合并可组合的多边形，得到简化后的多边形列表。所有多边形均为简单多边形，无方向环或自相交。尝试使用JTS（Java Topology Suite）进行合并，但结果不符合预期，已多次确认多边形定义正确，但合并后仍出现异常结果。以下是使用的代码：

public class PolygonMerger {
    
    public static List<double[]> mergePolygons(List<double[]> polygons) {
        List<Geometry> jtsPolygons = new ArrayList<>();
        GeometryFactory factory = new GeometryFactory();

        for (double[] coords : polygons) {
            Coordinate[] coordinates = new Coordinate[coords.length / 2];
            for (int i = 0; i < coords.length; i += 2) {
                coordinates[i / 2] = new Coordinate(coords[i], coords[i + 1]);
            }
            
            // Ensure the polygon is closed
            if (!coordinates[0].equals(coordinates[coordinates.length - 1])) {
                Coordinate[] closedCoordinates = new Coordinate[coordinates.length + 1];
                System.arraycopy(coordinates, 0, closedCoordinates, 0, coordinates.length);
                closedCoordinates[coordinates.length] = coordinates[0]; // Close the polygon
                coordinates = closedCoordinates;
            }

            // Create the polygon
            Polygon polygon = factory.createPolygon(coordinates);
            
            // Validate and repair the polygon if needed
            if (!new IsValidOp(polygon).isValid()) {
                polygon = (Polygon) polygon.buffer(0); // Repair polygon
            }
            
            jtsPolygons.add(polygon);
        }
        
        // Union all polygons
        Geometry mergedGeometry = CascadedPolygonUnion.union(jtsPolygons);
        List<Polygon> mergedPolygons = new ArrayList<>();
        
        if (mergedGeometry instanceof Polygon) {
            mergedPolygons.add((Polygon) mergedGeometry);
        } else if (mergedGeometry instanceof MultiPolygon) {
            for (int i = 0; i < mergedGeometry.getNumGeometries(); i++) {
                mergedPolygons.add((Polygon) mergedGeometry.getGeometryN(i));
            }
        }
        
        return convertToDoubleArrays(mergedPolygons);
    }
    
    private static List<double[]> convertToDoubleArrays(List<Polygon> polygons) {
        List<double[]> polygonCoordsList = new ArrayList<>();
        
        for (Polygon polygon : polygons) {
            Coordinate[] coordinates = polygon.getExteriorRing().getCoordinates();
            double[] coordsArray = new double[coordinates.length * 2];
            for (int i = 0; i < coordinates.length; i++) {
                coordsArray[i * 2] = coordinates[i].x;
                coordsArray[i * 2 + 1] = coordinates[i].y;
            }
            
            polygonCoordsList.add(coordsArray);
        }
        
        return polygonCoordsList;
    }
}

可能的异常原因及修复方法 #

1. 坐标精度问题 #

JTS对浮点数精度极度敏感，原始坐标的微小误差（比如计算时的浮点精度丢失）可能导致多边形拓扑关系判断错误，比如本该相邻的多边形被判定为不相交，合并后出现碎片或空洞。

• 修复：创建GeometryFactory时指定统一精度模型，强制坐标归一化：

  // 保留4位小数的固定精度模型，可根据需求调整精度量级
  PrecisionModel pm = new PrecisionModel(PrecisionModel.FIXED, 10000);
  GeometryFactory factory = new GeometryFactory(pm);
  

2. 多边形环方向不规范 #

JTS要求多边形外环必须为逆时针方向，内环（若存在）为顺时针方向。即使多边形无自相交，方向错误也会导致拓扑运算逻辑紊乱，比如合并后生成错误的多边形范围。

• 修复：创建多边形后统一环的方向：

  // 确保外环为逆时针
  if (!polygon.getExteriorRing().isCCW()) {
      polygon = (Polygon) polygon.reverse();
  }
  

3. CascadedPolygonUnion算法局限性 #

CascadedPolygonUnion在处理大量或空间关系复杂的多边形时，容易出现拓扑错误。可以改用更稳定的UnaryUnionOp实现合并。

• 修复：替换合并逻辑：

  // 替换原CascadedPolygonUnion.union调用
  Geometry mergedGeometry = UnaryUnionOp.union(jtsPolygons);
  

4. 多边形修复逻辑不严谨 #

仅用buffer(0)修复无效多边形虽然常用，但可能过度修改原始形状，尤其是当多边形存在微小缝隙或重叠时。建议先明确无效原因，再针对性处理。

• 修复：增加无效原因检查，确保修复有效性：

  IsValidOp isValidOp = new IsValidOp(polygon);
  if (!isValidOp.isValid()) {
      ValidationError error = isValidOp.getValidationError();
      polygon = (Polygon) polygon.buffer(0);
      // 修复后再次验证，避免无效多边形进入合并流程
      if (!new IsValidOp(polygon).isValid()) {
          throw new IllegalArgumentException("无法修复无效多边形: " + error.getMessage());
      }
  }
  

5. 坐标转换冗余问题 #

JTS返回的多边形环会强制闭合（首尾坐标相同），转换回double[]时会保留该冗余坐标，若输入格式要求不闭合或与原始格式一致，可能引发后续处理异常。

• 修复：在转换方法中移除重复的首尾坐标：

  private static List<double[]> convertToDoubleArrays(List<Polygon> polygons) {
      List<double[]> polygonCoordsList = new ArrayList<>();
      
      for (Polygon polygon : polygons) {
          Coordinate[] coordinates = polygon.getExteriorRing().getCoordinates();
          int coordCount = coordinates.length;
          // 移除重复的首尾坐标
          if (coordinates[0].equals(coordinates[coordCount - 1])) {
              coordCount--;
          }
          double[] coordsArray = new double[coordCount * 2];
          for (int i = 0; i < coordCount; i++) {
              coordsArray[i * 2] = coordinates[i].x;
              coordsArray[i * 2 + 1] = coordinates[i].y;
          }
          polygonCoordsList.add(coordsArray);
      }
      return polygonCoordsList;
  }
  

验证步骤 #

1. 在转换为JTS多边形后，输出每个多边形的isValid()状态和环方向，确认输入转换后的多边形有效性。
2. 合并前输出所有多边形的边界框，验证它们的空间关系是否符合预期。
3. 合并后输出mergedGeometry的类型、坐标信息，定位异常是出在合并阶段还是转换阶段。

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