这就是当你从 stackexchange 盲目复制代码时会发生的情况问题没有阅读question它基于哪个解释了原因。
我一直回答这个问题(并且发布这样的代码)提问者试图使用以度为单位的纬度/经度坐标来测量以米为单位的短距离。您问题中显示的技巧会创建一个自动 UTM 投影,以“AUTO:42001”位(在您的情况下为 52N 13E)之后指定的位置为中心 - 这需要是您感兴趣的区域的中心,因此在您的无论如何,这些值可能是错误的。
But您对孟买到德班的一个小区域不感兴趣,这是绕地球的一条重要路线,因此您需要考虑到地球表面的曲率。此外,您也不会尝试做一些困难的事情,而 JTS 是唯一的进程来源(例如缓冲)。在这种情况下,您应该使用GeodeticCalculator
它使用 C. F. F. Karney, Algorithms for geodesics, J. Geodesy 87, 43–55 (2013) 的库考虑地球的形状。
无论如何,足够的解释将来没有人会读,这里是代码:
public static void main(String[] args) {
DefaultGeographicCRS crs = DefaultGeographicCRS.WGS84;
if (args.length != 4) {
System.err.println("Need 4 numbers lat_1 lon_1 lat_2 lon_2");
return;
}
GeometryFactory geomFactory = new GeometryFactory();
Point[] points = new Point[2];
for (int i = 0, k = 0; i < 2; i++, k += 2) {
double x = Double.valueOf(args[k]);
double y = Double.valueOf(args[k + 1]);
if (CRS.getAxisOrder(crs).equals(AxisOrder.NORTH_EAST)) {
System.out.println("working with a lat/lon crs");
points[i] = geomFactory.createPoint(new Coordinate(x, y));
} else {
System.out.println("working with a lon/lat crs");
points[i] = geomFactory.createPoint(new Coordinate(y, x));
}
}
double distance = 0.0;
GeodeticCalculator calc = new GeodeticCalculator(crs);
calc.setStartingGeographicPoint(points[0].getX(), points[0].getY());
calc.setDestinationGeographicPoint(points[1].getX(), points[1].getY());
distance = calc.getOrthodromicDistance();
double bearing = calc.getAzimuth();
Quantity<Length> dist = Quantities.getQuantity(distance, SI.METRE);
System.out.println(dist.to(MetricPrefix.KILO(SI.METRE)).getValue() + " Km");
System.out.println(dist.to(USCustomary.MILE).getValue() + " miles");
System.out.println("Bearing " + bearing + " degrees");
}
Giving:
working with a lon/lat crs
POINT (72.8777 19.076)
POINT (31.049999 -29.883333)
7032.866960793305 Km
4370.020928274692 miles
Bearing -139.53428618565218 degrees