Changeset 2570 in josm
 Timestamp:
 20091204T16:30:52+01:00 (12 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/src/org/openstreetmap/josm/data/coor/LatLon.java
r2512 r2570 6 6 import java.text.DecimalFormat; 7 7 import java.text.NumberFormat; 8 9 import static java.lang.Math.*; 8 10 9 11 import org.openstreetmap.josm.Main; … … 121 123 /** 122 124 * Computes the distance between this lat/lon and another point on the earth. 123 * Uses spherical law of cosines formula, not Haversine.125 * Uses Haversine formular. 124 126 * @param other the other point. 125 127 * @return distance in metres. 126 128 */ 127 129 public double greatCircleDistance(LatLon other) { 128 return (Math.acos( 129 Math.sin(Math.toRadians(lat())) * Math.sin(Math.toRadians(other.lat())) + 130 Math.cos(Math.toRadians(lat()))*Math.cos(Math.toRadians(other.lat())) * 131 Math.cos(Math.toRadians(other.lon()lon()))) * 6378135); 130 double R = 6378135; 131 double sinHalfLat = sin(toRadians(other.lat()  this.lat()) / 2); 132 double sinHalfLon = sin(toRadians(other.lon()  this.lon()) / 2); 133 double d = 2 * R * asin( 134 sqrt(sinHalfLat*sinHalfLat + 135 cos(toRadians(this.lat()))*cos(toRadians(other.lat()))*sinHalfLon*sinHalfLon)); 136 // For points opposite to each other on the sphere, 137 // rounding errors could make the argument of asin greater than 1 138 // (This should almost never happen.) 139 if (java.lang.Double.isNaN(d)) { 140 System.err.println("Error: NaN in greatCircleDistance"); 141 d = PI * R; 142 } 143 return d; 132 144 } 133 145
Note: See TracChangeset
for help on using the changeset viewer.