Java中如何按指定距离放大带经纬度海拔坐标的多边形

核心思路 #

原始坐标为WGS84（EPSG:4326）球面经纬度坐标，1度经度对应的实际距离随纬度变化，无法直接在地理坐标系下做平面缓冲计算。需要先将多边形投影到适配当前区域的等距平面坐标系（如UTM坐标系，适配小范围区域），完成米级缓冲计算后再转换回WGS84经纬度坐标即可。

依赖引入 #

推荐使用业界成熟的空间运算库JTS Topology Suite做几何计算，配合GeoTools完成坐标系转换，Maven依赖如下：

<dependencies>
    <!-- JTS 空间几何运算库 -->
    <dependency>
        <groupId>org.locationtech.jts</groupId>
        <artifactId>jts-core</artifactId>
        <version>1.19.0</version>
    </dependency>
    <!-- GeoTools 坐标系转换工具 -->
    <dependency>
        <groupId>org.geotools</groupId>
        <artifactId>gt-main</artifactId>
        <version>28.2</version>
    </dependency>
    <dependency>
        <groupId>org.geotools</groupId>
        <artifactId>gt-referencing</artifactId>
        <version>28.2</version>
    </dependency>
</dependencies>

实现步骤 #

1. 解析原始坐标构造多边形 #

先将你提供的经纬度坐标解析为JTS的Polygon对象：

import org.locationtech.jts.geom.*;
import java.util.Arrays;

// 构造JTS几何工厂
GeometryFactory geometryFactory = new GeometryFactory(new PrecisionModel(), 4326);

// 解析示例坐标点
Coordinate[] coords = new Coordinate[]{
    new Coordinate(-0.4152899,39.4792891,0),
    new Coordinate(-0.4156305,39.4785997,0),
    new Coordinate(-0.4154294,39.478521,0),
    new Coordinate(-0.4150941,39.478935,0),
    new Coordinate(-0.4150914,39.4791048,0),
    new Coordinate(-0.4152899,39.4792891,0)
};
Polygon originalPolygon = geometryFactory.createPolygon(coords);

2. 坐标系转换+米级缓冲 #

将WGS84多边形投影到对应UTM坐标系，完成缓冲后再转换回WGS84：

import org.geotools.referencing.CRS;
import org.geotools.geometry.jts.JTS;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.operation.MathTransform;

// 计算当前多边形所在的UTM带号
int utmZone = (int) Math.floor((originalPolygon.getCentroid().getX() + 180) / 6) + 1;
// 判断南北半球，示例坐标纬度为北纬，所以用N对应的EPSG前缀326，南纬替换为327
String utmCrsCode = "EPSG:326" + utmZone; 

// 获取坐标系实例
CoordinateReferenceSystem wgs84 = CRS.decode("EPSG:4326");
CoordinateReferenceSystem utmCrs = CRS.decode(utmCrsCode);

// WGS84转UTM投影
MathTransform wgs84ToUtm = CRS.findMathTransform(wgs84, utmCrs);
Polygon utmPolygon = (Polygon) JTS.transform(originalPolygon, wgs84ToUtm);

// 执行缓冲，参数为要扩充的米数，比如要扩100米就填100
double bufferMeters = 100;
Polygon bufferedUtmPolygon = (Polygon) utmPolygon.buffer(bufferMeters, 8); // 第二个参数为圆滑度，值越大边越平滑

// UTM转WGS84得到最终缓冲后的经纬度多边形
MathTransform utmToWgs84 = CRS.findMathTransform(utmCrs, wgs84);
Polygon bufferedWgs84Polygon = (Polygon) JTS.transform(bufferedUtmPolygon, utmToWgs84);

// 输出结果坐标
Coordinate[] bufferedCoords = bufferedWgs84Polygon.getCoordinates();
Arrays.stream(bufferedCoords).forEach(coord -> 
    System.out.printf("%.7f,%.7f,%.0f%n", coord.x, coord.y, coord.getZ())
);

注意事项 #

• 如果多边形范围跨UTM分带，可自定义以多边形中心经线为中央经线的等距投影，避免跨带误差
• 缓冲运算默认会自动处理凹多边形尖角、自相交等异常情况，无需额外处理
• 海拔值可统一继承原始多边形的海拔参数，无需参与运算

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