如何由距离计算经纬度?经纬度计算距离极坐标?根据两点经纬度计算距离(2)

方法1：由于地球是椭球体，这个太难算了，如果假设地球是球体，可以使用以下公式：设地球上某点的经度为A,纬度为B， 则这点的空间坐标是 x=cos(B)*cos(A) y=cos(B)*sin(A) z=sin(B) 设地球上两点的空间坐标分别为(x1,y1,z1),(x2,y2,z2) 则它们的夹角为 C=acos(x1*x2y1*y2z1*z2)，C是角度 则两地距离为 C/180*pi*R，其中R为地球平均半径6371 误差不超过1%

地球是一个近乎标准的椭球体，它的赤道半径为6378.140千米，极半径为 6356.755千米，平均半径6371.004千米。如果我们假设地球是一个完美的球体，那么它的半径就是地球的平均半径，记为R。如果以0度经线为基 准，那么根据地球表面任意两点的经纬度就可以计算出这两点间的地表距离（这里忽略地球表面地形对计算带来的误差，仅仅是理论上的估算）。设第一点A的经 纬度为(LonA, LatA)，第二点B的经纬度为(LonB, LatB)，按照0度经线的基准，东经取经度的正(Longitude)，西经取经度负(-Longitude)，北纬取90-纬度(90- Latitude)，南纬取90纬度(90Latitude)，则经过上述处理过后的两点被计为(MLonA, MLatA)和(MLonB, MLatB)。那么根据三角推导，可以得到计算两点距离的如下公式：

C = sin(MLatA)*sin(MLatB)*cos(MLonA-MLonB) cos(MLatA)*cos(MLatB) Distance = R*Arccos(C)*Pi/180

这里，R和Distance单位是相同，如果是采用6371.004千米作为半径，那么Distance就是千米为单位，如果要使用其他单位，比如mile，还需要做单位换算，1千米=0.621371192mile

如果仅对经度作正负的处理，而不对纬度作90-Latitude(假设都是北半球，南半球只有澳洲具有应用意义)的处理，那么公式将是：

C = sin(LatA)*sin(LatB) cos(LatA)*cos(LatB)*cos(MLonA-MLonB)

Distance = R*Arccos(C)*Pi/180

以上通过简单的三角变换就可以推出。

如果三角函数的输入和输出都采用弧度，那么公式还可以写作：

C = sin(LatA*Pi/180)*sin(LatB*Pi/180) cos(LatA*Pi/180)*cos(LatB*Pi/180)*cos((MLonA-MLonB)*Pi/180)

Distance = R*Arccos(C)*Pi/180

也就是：

C = sin(LatA/57.2958)*sin(LatB/57.2958) cos(LatA/57.2958)*cos(LatB/57.2958)*cos((MLonA-MLonB)/57.2958)

Distance = R*Arccos(C) = 6371.004*Arccos(C) kilometer = 0.621371192*6371.004*Arccos(C) mile = 3958.758349716768*Arccos(C) mile

private static final double EARTH_RADIUS = 6378137;//赤道半径(单位m) /** * 转化为弧度(rad) * */ private static double rad(double d) { return d * Math.PI / 180.0; } /** * 基于余弦定理求两经纬度距离 * @param lon1 第一点的精度 * @param lat1 第一点的纬度 * @param lon2 第二点的精度 * @param lat3 第二点的纬度 * @return 返回的距离，单位km * */ public static double LantitudeLongitudeDist(double lon1, double lat1,double lon2, double lat2) { double radLat1 = rad(lat1); double radLat2 = rad(lat2); double radLon1 = rad(lon1); double radLon2 = rad(lon2); if (radLat1 < 0) radLat1 = Math.PI / 2 + Math.abs(radLat1);// south if (radLat1 > 0) radLat1 = Math.PI / 2 - Math.abs(radLat1);// north if (radLon1 < 0) radLon1 = Math.PI * 2 - Math.abs(radLon1);// west if (radLat2 < 0) radLat2 = Math.PI / 2 + Math.abs(radLat2);// south if (radLat2 > 0) radLat2 = Math.PI / 2 - Math.abs(radLat2);// north if (radLon2 < 0) radLon2 = Math.PI * 2 - Math.abs(radLon2);// west double x1 = EARTH_RADIUS * Math.cos(radLon1) * Math.sin(radLat1); double y1 = EARTH_RADIUS * Math.sin(radLon1) * Math.sin(radLat1); double z1 = EARTH_RADIUS * Math.cos(radLat1); double x2 = EARTH_RADIUS * Math.cos(radLon2) * Math.sin(radLat2); double y2 = EARTH_RADIUS * Math.sin(radLon2) * Math.sin(radLat2); double z2 = EARTH_RADIUS * Math.cos(radLat2); double d = Math.sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2)+ (z1 - z2) * (z1 - z2)); //余弦定理求夹角 double theta = Math.acos((EARTH_RADIUS * EARTH_RADIUS + EARTH_RADIUS * EARTH_RADIUS - d * d) / (2 * EARTH_RADIUS * EARTH_RADIUS)); double dist = theta * EARTH_RADIUS; return dist; }

方法2：google地图提供的方法：

对上面的公式解释如下：
1.Lat1 Lung1 表示A点经纬度，Lat2 Lung2 表示B点经纬度；
2.a=Lat1 – Lat2 为两点纬度之差 b=Lung1 -Lung2 为两点经度之差；
3.6378.137为地球半径，单位为千米；
计算出来的结果单位为千米，若将半径改为米为单位则计算的结果单位为米。
计算精度与谷歌地图的距离精度差不多，相差范围在0.2米以下。
private static final double EARTH_RADIUS = 6378137;//赤道半径(单位m) /** * 转化为弧度(rad) * */ private static double rad(double d) { return d * Math.PI / 180.0; } /** * 基于googleMap中的算法得到两经纬度之间的距离,计算精度与谷歌地图的距离精度差不多，相差范围在0.2米以下 * @param lon1 第一点的精度 * @param lat1 第一点的纬度 * @param lon2 第二点的精度 * @param lat3 第二点的纬度 * @return 返回的距离，单位km * */ public static double GetDistance(double lon1,double lat1,double lon2, double lat2) { double radLat1 = rad(lat1); double radLat2 = rad(lat2); double a = radLat1 - radLat2; double b = rad(lon1) - rad(lon2); 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 * EARTH_RADIUS; //s = Math.round(s * 10000) / 10000; return s; }

以上就是关于经纬度计算距离的全部内容，相信你一定会非常满意。

本文来自电脑杂谈，转载请注明本文网址：
http://www.pc-fly.com/a/jisuanjixue/article-19181-2.html