Mercurial > dive4elements > gnv-client
comparison gnv-artifacts/src/main/java/de/intevation/gnv/utils/DistanceCalculator.java @ 185:5fc8f41669a6
Added Calculation of the Distance in the unit km for Geodetic-Coordinates according to the Specification GT_0030.004
gnv-artifacts/trunk@229 c6561f87-3c4e-4783-a992-168aeb5c3f6f
| author | Tim Englich <tim.englich@intevation.de> |
|---|---|
| date | Fri, 16 Oct 2009 12:37:00 +0000 |
| parents | |
| children | da1499a464b9 |
comparison
equal
deleted
inserted
replaced
| 184:5d050f6de41d | 185:5fc8f41669a6 |
|---|---|
| 1 /** | |
| 2 * | |
| 3 */ | |
| 4 package de.intevation.gnv.utils; | |
| 5 | |
| 6 import com.vividsolutions.jts.geom.Point; | |
| 7 | |
| 8 /** | |
| 9 * @author Tim Englich <tim.englich@intevation.de> | |
| 10 * | |
| 11 */ | |
| 12 public class DistanceCalculator { | |
| 13 | |
| 14 private final static double flattening = 1.0 / 298.257233563; | |
| 15 | |
| 16 private final static double earthRadius = 6378137.0 / 1000.0 ; | |
| 17 | |
| 18 /** | |
| 19 * Constructor | |
| 20 */ | |
| 21 public DistanceCalculator() { | |
| 22 } | |
| 23 | |
| 24 public double calculateDistance(Point p1, Point p2){ | |
| 25 double resultValue = 0.0; | |
| 26 | |
| 27 double b1 = p1.getY(); | |
| 28 double b2 = p2.getY(); | |
| 29 | |
| 30 double l1 = p1.getX(); | |
| 31 double l2 = p2.getX(); | |
| 32 | |
| 33 | |
| 34 double F = (b1 + b2) / 2.0; | |
| 35 double G = (b1 - b2) / 2.0; | |
| 36 double l = (l1 - l2) / 2.0; | |
| 37 | |
| 38 F = (Math.PI / 180.0) * F; | |
| 39 G = (Math.PI / 180.0) * G; | |
| 40 l = (Math.PI / 180.0) * l; | |
| 41 | |
| 42 double S = ((Math.sin(G) * Math.sin(G)) * ((Math.cos(l) * Math.cos(l))))+ | |
| 43 ((Math.cos(F) * Math.cos(F)) * ((Math.sin(l) * Math.sin(l)))); | |
| 44 | |
| 45 double C = ((Math.cos(G) * Math.cos(G)) * ((Math.cos(l) * Math.cos(l))))+ | |
| 46 ((Math.sin(F) * Math.sin(F)) * ((Math.sin(l) * Math.sin(l)))); | |
| 47 | |
| 48 double w = Math.atan(Math.sqrt((S/C))); | |
| 49 | |
| 50 double D = 2.0 * w * earthRadius; | |
| 51 | |
| 52 double R = Math.sqrt((S*C)) / w; | |
| 53 | |
| 54 double H1 = (3.0 * R - 1.0 ) / (2.0 * C); | |
| 55 double H2 = (3.0 * R + 1.0 ) / (2.0 * S); | |
| 56 | |
| 57 resultValue = D * (1 + (flattening * H1 * (Math.sin(F) * Math.sin(F)) * | |
| 58 (Math.cos(G) * Math.cos(G))) - | |
| 59 (flattening * H2 * (Math.cos(F) * Math.cos(F)) * | |
| 60 (Math.sin(G) * Math.sin(G)))); | |
| 61 | |
| 62 return resultValue; | |
| 63 } | |
| 64 | |
| 65 } |