节选自百度地图API:
若开发者使用的是百度地图或定位API,且版本在1.3.5以后的,
路线规划提供了获取路线距离的方法,见MKRoutePlan 类的 getDistance 方法。
如果是计算任意两点的距离,自2.0.0版本开始,Android SDK为开发者提供了计算距离的接口(DistanceUtil),具体使用方法如下:
GeoPoint p1LL = new GeoPoint(39971802, 116347927);
GeoPoint p2LL = new GeoPoint(39892131, 116498555);
double distance = DistanceUtil.getDistance(p1LL, p2LL);
若未使用百度地图或定位SDK,或是地图、定位SDK版本在1.3.5以下的,可以选用以下两种计算方式:
一、利用勾股定理计算,适用于两点距离很近的情况;
二、按标准的球面大圆劣弧长度计算,适用于距离较远的情况;
static double DEF_PI = 3.14159265359; // PI
static double DEF_2PI = 6.28318530712; // 2*PI
static double DEF_PI180 = 0.01745329252; // PI/180.0
static double DEF_R = 6370693.5; // radius of earth
/**
* 第一种方式,按勾股定律求结果
* @param lon1
* @param lat1
* @param lon2
* @param lat2
* @return
* @author kazeik.chen QQ:77132995 2014-4-1下午4:30:09
* TODO kazeik@163.com
*/
public double GetShortDistance(double lon1, double lat1, double lon2,
double lat2) {
double ew1, ns1, ew2, ns2;
double dx, dy, dew;
double distance;
// 角度转换为弧度
ew1 = lon1 * DEF_PI180;
ns1 = lat1 * DEF_PI180;
ew2 = lon2 * DEF_PI180;
ns2 = lat2 * DEF_PI180;
// 经度差
dew = ew1 - ew2;
// 若跨东经和西经180 度,进行调整
if (dew > DEF_PI)
dew = DEF_2PI - dew;
else if (dew < -DEF_PI)
dew = DEF_2PI + dew;
dx = DEF_R * Math.cos(ns1) * dew; // 东西方向长度(在纬度圈上的投影长度)
dy = DEF_R * (ns1 - ns2); // 南北方向长度(在经度圈上的投影长度)
// 勾股定理求斜边长
distance = Math.sqrt(dx * dx + dy * dy);
return distance;
}
/**
* 第二种方式,按大圆劣弧求距离
* @param lon1
* @param lat1
* @param lon2
* @param lat2
* @return
* @author kazeik.chen QQ:77132995 2014-4-1下午4:30:30
* TODO kazeik@163.com
*/
public double GetLongDistance(double lon1, double lat1, double lon2,
double lat2) {
double ew1, ns1, ew2, ns2;
double distance;
// 角度转换为弧度
ew1 = lon1 * DEF_PI180;
ns1 = lat1 * DEF_PI180;
ew2 = lon2 * DEF_PI180;
ns2 = lat2 * DEF_PI180;
// 求大圆劣弧与球心所夹的角(弧度)
distance = Math.sin(ns1) * Math.sin(ns2) + Math.cos(ns1)
* Math.cos(ns2) * Math.cos(ew1 - ew2);
// 调整到[-1..1]范围内,避免溢出
if (distance > 1.0)
distance = 1.0;
else if (distance < -1.0)
distance = -1.0;
// 求大圆劣弧长度
distance = DEF_R * Math.acos(distance);
return distance;
}
//以下是测试
double mLat1 = 39.90923; // point1纬度
double mLon1 = 116.357428; // point1经度
double mLat2 = 39.90923;// point2纬度
double mLon2 = 116.397428;// point2经度
double distance = GetShortDistance(mLon1, mLat1, mLon2, mLat2);