http://www.blogjava.net/sinoly/archive/2007/01/22/95406.html
研究了一下GoogleMap上的根据经纬度求地球表面两点间距离的实现， 用java实现了一把，对我国境内的Beijing54, Xian80,WGS84三种坐标系的空间距离计算感觉这个实现是比较准确的。当然，这里只是个人的感觉而已。。。可能我还没有遇到那种对精度要求非常严格的地方，或许，本身我这个方式就是错误的。。。呵呵。。。这里只是作个记号先。。。各位达人多多指教
package com.geotools.test;
/**
*
* CopyRight (C) All rights reserved.
*
* WuHan Inpoint Information Technology Development,Inc.
*
* Author sinoly
Project Name: PostGeo
*
* @version 1.0 2006-11-13
*
*
Base on : JDK1.5
*
*/
public class GeoUtils {
public enum GaussSphere{
Beijing54,
Xian80,
WGS84,
}
private static double Rad(double d){
return d * Math.PI / 180.0;
}
public double DistanceOfTwoPoints(double lng1,double lat1,double lng2,double lat2,
GaussSphere gs){
double radLat1 = Rad(lat1);
double radLat2 = Rad(lat2);
double a = radLat1 - radLat2;
double b = Rad(lng1) - Rad(lng2);
double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a/2),2) +
Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(Math.sin(b/2),2)));
s = s * (gs == GaussSphere.WGS84 ? 6378137.0 : (gs == GaussSphere.Xian80 ? 6378140.0 : 6378245.0));
s = Math.round(s * 10000) / 10000;
return s;
}
}
posted on 2007-01-22 20:38 sinoly
(#)
http://www.blogjava.net/sinoly/archive/2007/01/22/95406.html
研究了一下GoogleMap上的根据经纬度求地球表面两点间距离的实现， 用java实现了一把，对我国境内的Beijing54, Xian80,WGS84三种坐标系的空间距离计算感觉这个实现是比较准确的。当然，这里只是个人的感觉而已。。。可能我还没有遇到那种对精度要求非常严格的地方，或许，本身我这个方式就是错误的。。。呵呵。。。这里只是作个记号先。。。各位达人多多指教
package com.geotools.test;
/**
*
* CopyRight (C) All rights reserved.
*
* WuHan Inpoint Information Technology Development,Inc.
*
* Author sinoly
Project Name: PostGeo
*
* @version 1.0 2006-11-13
*
*
Base on : JDK1.5
*
*/
public class GeoUtils {
public enum GaussSphere{
Beijing54,
Xian80,
WGS84,
}
private static double Rad(double d){
return d * Math.PI / 180.0;
}
public double DistanceOfTwoPoints(double lng1,double lat1,double lng2,double lat2,
GaussSphere gs){
double radLat1 = Rad(lat1);
double radLat2 = Rad(lat2);
double a = radLat1 - radLat2;
double b = Rad(lng1) - Rad(lng2);
double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a/2),2) +
Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(Math.sin(b/2),2)));
s = s * (gs == GaussSphere.WGS84 ? 6378137.0 : (gs == GaussSphere.Xian80 ? 6378140.0 : 6378245.0));
s = Math.round(s * 10000) / 10000;
return s;
}
}
posted on 2007-01-22 20:38 sinoly
(#)