如何使用Java实现2D Alpha Shapes算法提取二维点集轮廓

直接用JTS Topology Suite内置的Alpha Shapes实现就行，不用从零手写算法，自己写很容易在边界判断、共线点处理上出bug，这个库是Java生态里几何计算最成熟的，稳定性经过大量生产场景验证。

依赖引入 #

Maven项目在pom.xml里加如下依赖：

<dependency>
    <groupId>org.locationtech.jts</groupId>
    <artifactId>jts-core</artifactId>
    <version>1.19.0</version>
</dependency>

Gradle项目引入方式：

implementation 'org.locationtech.jts:jts-core:1.19.0'

核心实现代码 #

JTS 1.15及以上版本已经内置了Alpha Shape计算能力，封装好的工具类直接用就行，返回的坐标序列可以直接传给OpenCV做绘制：

import org.locationtech.jts.algorithm.hull.AlphaShape;
import org.locationtech.jts.geom.Coordinate;
import org.locationtech.jts.geom.Geometry;
import org.locationtech.jts.geom.GeometryFactory;
import org.locationtech.jts.geom.MultiPoint;
import java.util.List;

public class PointSetContourUtil {
    private static final GeometryFactory GEOM_FACTORY = new GeometryFactory();

    /**
     * 提取二维点集的Alpha轮廓
     * @param pointList 输入点集，每个元素为[x,y]格式的坐标数组，坐标系和你采集点的坐标系一致
     * @param alpha 轮廓紧致度参数：值越小轮廓越贴合点的凹凸细节，值越大轮廓越趋近于凸包
     * @return 按顺序排列的轮廓点坐标，可直接用于绘制
     */
    public static Coordinate[] getContour(List<double[]> pointList, double alpha) {
        // 转换输入点为JTS要求的几何格式
        Coordinate[] coordinates = pointList.stream()
                .map(point -> new Coordinate(point[0], point[1]))
                .toArray(Coordinate[]::new);
        MultiPoint multiPoint = GEOM_FACTORY.createMultiPointFromCoords(coordinates);

        // 计算Alpha Shape
        AlphaShape alphaShape = new AlphaShape(multiPoint);
        alphaShape.setAlpha(alpha);
        Geometry resultShape = alphaShape.getResult();

        // 返回轮廓边界的点序列
        return resultShape.getBoundary().getCoordinates();
    }
}

针对512*512尺寸图像点集的调优说明 #

• alpha参数初始值可以先设为15，根据出图效果调整：如果轮廓毛刺太多、把离群噪声点也包进去了，就把alpha往大调（比如调到20_{30）；如果轮廓太圆滑，丢失了凹陷细节，就把alpha往小调（比如调到8}12）
• 如果采集的点集有明显离群噪声，先做一步离群点过滤再跑算法，不然轮廓会被离群点拉得变形
• 拿到返回的Coordinate[]后，直接遍历取x、y值转成整数，传入OpenCV的Imgproc.polylines方法就能画出闭合轮廓，不需要额外做坐标转换

踩坑提醒：不要用网上零散的手写Alpha Shapes代码，大部分没处理共线点、轮廓自交、孔洞判断的逻辑，跑实际采集的点集很容易出现轮廓断裂、形状异常的问题。

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