daoas/fulddaekning/GeoPointHelper.java

package dk.daoas.fulddaekning;


public class GeoPointHelper {


        // Latitude (horizonal), longitude(vertical) so
        // 1 degree latitude is ~ 111320 meters, since the distance between the
        // horizonal lines is always the same
        // 1 degree longitude is ~111320 meters at equator but gets shorter as we
        // get closer to the poles.
        // so 1 degree longitude is 64.5 km at denmarks southern point
        // (gedser=54.55,11.95)
        // and is 59.4km at northern point (skagen = 57.75,10.65)

        public static double kmToLatitude(double km) {
                return km / 111.320;
        }

        public static double kmToLongitude(double km) {// denne er kun ca
                return km / 62.0;
        }
        
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

        /* ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: */
        /* :: This function converts decimal degrees to radians : */
        /* ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: */
        private static double deg2rad(double deg) {
                return (deg * Math.PI / 180.0);
        }

        /* ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: */
        /* :: This function converts radians to decimal degrees : */
        /* ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: */
        private static double rad2deg(double rad) {
                return (rad * 180 / Math.PI);
        }

        // http://www.geodatasource.com/developers/java
        private static double distanceHaversine(double lat1, double lon1,
                        double lat2, double lon2) {
                double theta = lon1 - lon2;
                double dist = Math.sin(deg2rad(lat1)) * Math.sin(deg2rad(lat2))
                                + Math.cos(deg2rad(lat1)) * Math.cos(deg2rad(lat2))
                                * Math.cos(deg2rad(theta));
                dist = Math.acos(dist);
                dist = rad2deg(dist);
                dist = dist * 60 * 1.1515;

                // indtil nu er dist i miles - så vi omregner lige til km
                dist = dist * 1.609344;
                return (dist);
        }

        public static double beregnAfstand(Adresse p1, Adresse p2) {
                if (p1.latitude == p2.latitude && p1.longitude == p2.longitude) {
                        return 0.00; //Haversine går i baglås hvis p1 og p2 er samme position
                }
                
                
                return distanceHaversine(p1.latitude, p1.longitude, p2.latitude, p2.longitude);
        }
        
}
1	torben	2585	package dk.daoas.fulddaekning;
2
3
4	torben	2748
5	torben	2585	public class GeoPointHelper {
6
7
8
9			// Latitude (horizonal), longitude(vertical) so
10			// 1 degree latitude is ~ 111320 meters, since the distance between the
11			// horizonal lines is always the same
12			// 1 degree longitude is ~111320 meters at equator but gets shorter as we
13			// get closer to the poles.
14			// so 1 degree longitude is 64.5 km at denmarks southern point
15			// (gedser=54.55,11.95)
16			// and is 59.4km at northern point (skagen = 57.75,10.65)
17
18			public static double kmToLatitude(double km) {
19			return km / 111.320;
20			}
21
22			public static double kmToLongitude(double km) {// denne er kun ca
23			return km / 62.0;
24			}
25
26			////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
27
28			/* ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: */
29			/* :: This function converts decimal degrees to radians : */
30			/* ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: */
31			private static double deg2rad(double deg) {
32			return (deg * Math.PI / 180.0);
33			}
34
35			/* ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: */
36			/* :: This function converts radians to decimal degrees : */
37			/* ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: */
38			private static double rad2deg(double rad) {
39			return (rad * 180 / Math.PI);
40			}
41
42			// http://www.geodatasource.com/developers/java
43			private static double distanceHaversine(double lat1, double lon1,
44			double lat2, double lon2) {
45			double theta = lon1 - lon2;
46			double dist = Math.sin(deg2rad(lat1)) * Math.sin(deg2rad(lat2))
47			+ Math.cos(deg2rad(lat1)) * Math.cos(deg2rad(lat2))
48			* Math.cos(deg2rad(theta));
49			dist = Math.acos(dist);
50			dist = rad2deg(dist);
51			dist = dist * 60 * 1.1515;
52
53			// indtil nu er dist i miles - så vi omregner lige til km
54			dist = dist * 1.609344;
55			return (dist);
56			}
57
58	torben	2707	public static double beregnAfstand(Adresse p1, Adresse p2) {
59	torben	2585	if (p1.latitude == p2.latitude && p1.longitude == p2.longitude) {
60			return 0.00; //Haversine går i baglås hvis p1 og p2 er samme position
61			}
62
63
64			return distanceHaversine(p1.latitude, p1.longitude, p2.latitude, p2.longitude);
65			}
66
67			}