Question

我想让geotools告诉我两点之间的距离。我提供度数的积分，我想以米为单位回到距离。我按照这个网页进行了说明：

http://docs.geotools.org/stable/userguide/library/api/jts.html

不幸的是，当我运行它时，会抛出异常，因为它表示我的纬度超出了范围。

这是我的代码：

public double measureDistanceInMeters(double[] coord1, double[] coord2) {
    System.out.println("latitude: "+coord1[0]+", longitude: "+coord1[1]);
    System.out.println("latitude: "+coord2[0]+", longitude: "+coord2[1]);
    CoordinateReferenceSystem crs = DefaultGeocentricCRS.CARTESIAN;
    // latitude first, longitude second.
    Coordinate start = new Coordinate(coord1[0],coord1[1]);
    Coordinate end = new Coordinate(coord2[0],coord2[1]);
    double distance = -1;
    try {
        distance = JTS.orthodromicDistance(start, end, crs);
    } catch (TransformException e1) {
        e1.printStackTrace();
    }
    return distance;

当我运行它时，控制台会打印出来：

latitude: 38.8951, longitude: -77.0367
latitude: 40.7127, longitude: -74.0059

JUnit失败，这是堆栈跟踪：

java.lang.IllegalArgumentException: Latitude �N is out of range (±90°).
    at org.geotools.referencing.GeodeticCalculator.checkLatitude(GeodeticCalculator.java:389)
    at org.geotools.referencing.GeodeticCalculator.setStartingGeographicPoint(GeodeticCalculator.java:550)
    at org.geotools.referencing.GeodeticCalculator.setStartingPosition(GeodeticCalculator.java:591)
    at org.geotools.geometry.jts.JTS.orthodromicDistance(JTS.java:635)
    at fungle.funfinder.data.geo.GeoToolsRuler.measureDistanceInMeters(GeoToolsRuler.java:49)
    at fungle.funfinder.data.geo.GeoRulerBase.getDistance(GeoRulerBase.java:15)
    at fungle.funfinder.data.geo.GeoToolsRulerTest.test(GeoToolsRulerTest.java:11)
    at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
    at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:62)
    at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43)
    at java.lang.reflect.Method.invoke(Method.java:483)
    at org.junit.runners.model.FrameworkMethod$1.runReflectiveCall(FrameworkMethod.java:47)
    at org.junit.internal.runners.model.ReflectiveCallable.run(ReflectiveCallable.java:12)
    at org.junit.runners.model.FrameworkMethod.invokeExplosively(FrameworkMethod.java:44)
    at org.junit.internal.runners.statements.InvokeMethod.evaluate(InvokeMethod.java:17)
    at org.junit.runners.ParentRunner.runLeaf(ParentRunner.java:271)
    at org.junit.runners.BlockJUnit4ClassRunner.runChild(BlockJUnit4ClassRunner.java:70)
    at org.junit.runners.BlockJUnit4ClassRunner.runChild(BlockJUnit4ClassRunner.java:50)
    at org.junit.runners.ParentRunner$3.run(ParentRunner.java:238)
    at org.junit.runners.ParentRunner$1.schedule(ParentRunner.java:63)
    at org.junit.runners.ParentRunner.runChildren(ParentRunner.java:236)
    at org.junit.runners.ParentRunner.access$000(ParentRunner.java:53)
    at org.junit.runners.ParentRunner$2.evaluate(ParentRunner.java:229)
    at org.junit.runners.ParentRunner.run(ParentRunner.java:309)
    at org.eclipse.jdt.internal.junit4.runner.JUnit4TestReference.run(JUnit4TestReference.java:50)
    at org.eclipse.jdt.internal.junit.runner.TestExecution.run(TestExecution.java:38)
    at org.eclipse.jdt.internal.junit.runner.RemoteTestRunner.runTests(RemoteTestRunner.java:467)
    at org.eclipse.jdt.internal.junit.runner.RemoteTestRunner.runTests(RemoteTestRunner.java:683)
    at org.eclipse.jdt.internal.junit.runner.RemoteTestRunner.run(RemoteTestRunner.java:390)
    at org.eclipse.jdt.internal.junit.runner.RemoteTestRunner.main(RemoteTestRunner.java:197)

所以我觉得我在这里误解了一些东西。 Geotools非常庞大，即使它有文档，也很难找到你正在寻找的文档，因为它非常庞大和复杂。

我猜我在这里使用了错误的坐标参考系统，或者作为一个arg传递到我不应该的地方。我不知道我做错了什么。请帮忙。

Answer 1

您使用的是错误的坐标系，基本上您告诉JTS您的坐标需要转换为纬度/经度然后使用，以便它们超出范围。

将代码更改为包括：

    CoordinateReferenceSystem crs = null;
    try {
        crs = CRS.decode("epsg:4326");
    } catch (NoSuchAuthorityCodeException e) {
        e.printStackTrace();
    } catch (FactoryException e) {
        e.printStackTrace();
    }

我得到了328739.66米或328Km的答案，这听起来是正确的（如果不是那么你可能需要将纬度和经度值交换成圆形）。

您还需要添加

    <dependency>
        <groupId>org.geotools</groupId>
        <artifactId>gt-epsg-hsql</artifactId>
        <version>${geotools.version}</version>
    </dependency>

到你的pom.xml文件或找不到坐标系的定义。

GeoTools JTS orthodromicDistance抛出纬度超出范围的异常

1 个答案: