Calculate distance with JTS

Question

I'm trying to get the distance between two points with JTS. The problem is that this function with JTS

new Coordinate(Lon1,Lat1).distance(new Coordinate(Lon2,Lat2))

Is giving a different result from this function used in postgis(wich provides the good distance)

ST_Distance_Sphere(point1,point2)

possible duplicate of java vividsolutions jts wgs-84 distance to meters — nicolas-f, Jul 30 '15 at 15:20

score 8 · Accepted Answer · edited Apr 29 '20 at 05:40

The short answer is you can't do that unless your points are very close together and you want the answer in degrees. JTS knows nothing about units or the curvature of the earth. So you need to pull in some GeoTools jars that do know about such things. Then you can create a method like:

private void calculateDistance(CoordinateReferenceSystem crs, Point[] points) {
    if (crs == null) {
        crs = default_crs;
    }

    double distance = 0.0;
    try {
        distance = JTS.orthodromicDistance(
            points[0].getCoordinate(),
            points[1].getCoordinate(),
            crs
        );
    } catch (TransformException e) {
        // TODO Auto-generated catch block
        e.printStackTrace();
    }

    Measure<Double, Length> dist = Measure.valueOf(distance, SI.METER);
    System.out.println(dist.doubleValue(SI.KILOMETER) + " Km");
    System.out.println(dist.doubleValue(NonSI.MILE) + " miles");
}

Which will take care of all the hard maths for you. The full example program can be seen here.

Mike T · Answer 2 · 2015-07-31T00:26:14.600

4

You could use the Java implementation of GeographicLib to solve the inverse geodesic problem on WGS84.

import net.sf.geographiclib.*;

...

    private double calculateDistance(
            double lon1, double lat1, double lon2, double lat2) {
        GeodesicData g = Geodesic.WGS84.Inverse(lat1, lon1, lat2, lon2);
        return g.s12;  // distance in metres
    }

The length of the geodesic, which represents the shortest distance in metres between the two points, is in g.s12.

edited Jul 31 '15 at 00:26

answered Jul 30 '15 at 23:25

Mike T

42,095
10
126
187

1

Note that GeoTools is using the same library under the hood – Ian Turton Feb 09 '18 at 15:00
Thanks for the note @IanTurton, I was debating if I should pull in all of GeoTools just for this or if I could just use this smaller lib instead. – carusot42 Nov 22 '19 at 17:07

score 1 · Answer 3 · answered Nov 22 '19 at 17:11

In case anyone wants a Scala implementation using GeographicLib and JTS points:

Add to build.sbt

libraryDependencies += "net.sf.geographiclib" % "GeographicLib-Java" % "1.50"

Distance.scala

import net.sf.geographiclib.Geodesic

object Distance {
  def calculateDistance(pt1: Point, pt2: Point): Double = {
     Geodesic.WGS84.Inverse(pt1.getY, pt1.getX, pt2.getY, pt2.getX).s12
  }
}

score 0 · Answer 4 · answered Jul 31 '15 at 08:42

Like the others already said, JTS does know nothing about geodesy or projections or coordinate systems. JTS does "linear geometry on the 2-dimensional Cartesian plane" (http://tsusiatsoftware.net/jts/main.html).

If you do not want to pull in another library and if the strict mathematical solution on a sphere is enough, you can simply implement the haversine formula in your code. It is very straight-forward (even if Wikipedia makes it look complicated if you are not used to mathematical formulas like me): https://en.wikipedia.org/wiki/Haversine_formula & http://rosettacode.org/wiki/Haversine_formula#Java

The earth is not a sphere though, so if you need high accuracy you probably don't want to do this.

Calculate distance with JTS

4 Answers4