如何到达WKT区域？

Question

我使用jts.jar，package是com.vividsolutions.jts.geom，得到一个几何区域。但是getArea()方法没有给出正确的结果。我的代码是

public static void main(String[] args) throws ParseException {
        GeometryFactory geometryFactory = new GeometryFactory();
        WKTReader reader = new WKTReader(geometryFactory);
        String s = "MULTIPOLYGON (((114.273193 40.480272, 114.274645 
        MultiPolygon mpolygon = null;
            mpolygon = (MultiPolygon) reader.read(s);
            System.out.println(mpolygon.getArea());
    }

Answer 1

在寻找类似问题的答案时，我发现了这个问题，于是我想分享一下我的发现：

在解释您试图解决的问题时，从博客文章中获取的描述：

GIS上反复出现的问题之一是“JTS getArea()计算的面积是多少”。似乎不重要的是，多少次人们说，这些都不是答案，因为这取决于你的数据投影。这个问题不会消失，所以我编写了一些代码，应该给大多数多边形一个贴切的答案。 与大多数这些问题一样，诀窍是将多边形转换为平坦的笛卡儿平面(这是JTS工作的最佳位置)。我们可以使用GeoTools自动投影(假设多边形足够小)，然后简单地调用.getArea()方法。

解决方案：

public static void main(String[] args) throws Exception {
  GeometryFactory geometryFactory = new GeometryFactory();
  WKTReader reader = new WKTReader(geometryFactory);
  String s = "MULTIPOLYGON (((40 40, 20 45, 45 30, 40 40)), ((20 35, 10 30, 10 10, 30 5, 45 20, 20 35), (30 20, 20 15, 20 25, 30 20)))";
  MultiPolygon mpolygon = (MultiPolygon) reader.read(s);
  Point centroid = mpolygon.getCentroid();
  String autoCode = "AUTO:42001," + centroid.getX() + "," + centroid.getY();
  CoordinateReferenceSystem auto = CRS.decode(autoCode);
  MathTransform transform = CRS.findMathTransform(DefaultGeographicCRS.WGS84, auto);
  MultiPolygon projed = (MultiPolygon) JTS.transform(mpolygon, transform);
  Measure<Double, Area> measure = Measure.valueOf(projed.getArea(), SI.SQUARE_METRE);
  System.out.println(measure);
}

输出：

8.146055550674082E12 m²

参考资料：https://blog.ianturton.com/geotools,/projections/2017/08/01/area-of-a-polygon.html