计算两个几个形状 Geometry 的面积占比_org.locationtech.jts.geom.polygon.getarea()方法得出的面积-优快云博客

本文链接：https://blog.youkuaiyun.com/qiaobing1226/article/details/126995246

本文介绍了一种使用LocationTech库计算两个多边形交集面积占总面积比例的方法，并提供了将坐标从WGS84转换到墨卡托投影进行精确面积计算的技术细节。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

1.使用locationtech自带方法（默认墨卡托坐标）

assert geometry1 != null;
assert geometry != null;Geometry intersection = geometry1.intersection(geometry);
double proportion = 0;
if (intersection != null) {
    //面积占比
    proportion = intersection.getArea() / geometry.getArea();
}
System.out.printf("面积占比为：",proportion);

2.坐标转换成墨卡托坐标

import org.geotools.geometry.jts.JTS;
import org.geotools.referencing.CRS;
import org.locationtech.jts.geom.Geometry;
import org.locationtech.jts.geom.Polygon;
import org.opengis.referencing.FactoryException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.operation.MathTransform;

    /**
     * 计算出多边形面积，如果为点或者线的话，结算结果为0（Ploygon不为0）
     *
     * @param geometry 多边形对象
     * @return 面积 单位（平方米）
     */
    private double getArea(Geometry geometry) throws FactoryException, TransformException {
        
         // Pseudo-Mercator(墨卡托投影)
        String decode57 = "EPSG:3857";
        // WGS84(一般项目中常用的是CSR:84和EPSG:4326)
        String decode84 = "EPSG:4326";
        CoordinateReferenceSystem sourceCRS = CRS.decode(decode84);
        CoordinateReferenceSystem targetCRS = CRS.decode(decode57);

        MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS, false);
        Geometry geometryMercator = JTS.transform(geometry, transform);
        return geometryMercator.getArea();
    }