原文转自:
http://blog.youkuaiyun.com/xjanker2/article/details/5404487
public class FCDTools {
public final static double PI = 3.14159265354;
private final static double D2R = 0.017453 ;
private final static double a2 = 6378137.0;
private final static double e2 = 0.006739496742337;
public static double Distance(FCD_Point pt1, FCD_Point pt2) {
if(pt1.lng == pt2.lng && pt1.lat == pt2.lat ) {
return 0.0;
} else {
double fdLambda =(pt1.lng - pt2.lng) * D2R;
double fdPhi = (pt1.lat - pt2.lat) * D2R;
double fPhimean = ((pt1.lat + pt2.lat) / 2.0) * D2R;
double fTemp = 1 - e2 * (Math.pow (Math.sin(fPhimean),2));
double fRho = (a2 * (1 - e2)) / Math.pow (fTemp, 1.5);
double fNu = a2 / (Math.sqrt(1 - e2 * (Math.sin(fPhimean) * Math.sin(fPhimean))));
double fz = Math.sqrt (Math.pow(Math.sin(fdPhi / 2.0), 2) +
Math.cos(pt2.lat * D2R) * Math.cos( pt1.lat*D2R ) * Math.pow( Math.sin(fdLambda / 2.0),2));
fz = 2 * Math.asin(fz);
double fAlpha = Math.cos(pt2.lat * D2R) * Math.sin(fdLambda) * 1 / Math.sin(fz);
fAlpha = Math.asin (fAlpha);
double fR = (fRho * fNu) / ((fRho * Math.pow( Math.sin(fAlpha),2)) + (fNu * Math.pow( Math.cos(fAlpha),2)));
return fz * fR;
}
}
}
测试:
System.out.println(FCDTools.Distance(new FCD_Point(117.107277,31.980298),
new FCD_Point(117.524757, 31.888227)));
// 输出40.7公里
本文介绍了一种基于经纬度的地理坐标距离计算方法,通过公式实现两点间的大地球面距离计算,输出示例展示了算法在实际应用中的使用。
657

被折叠的 条评论
为什么被折叠?



