问题引出：

今天遇到经纬度坐标转换距离的工作，根据网站登录者的IP确定登录者目前的位置信息，将其经纬度信息与所有的营业厅的经纬度进行对比，网页上显示出距离登录者最近的营业厅地址，本打算就做一个二维坐标距离转换，上网一搜，经纬度距离的计算远远没有想象的那样简单，最终得到如下的公式计算。

经纬度坐标的计算公式为：

1. Lat1 Lung1 表示A点经纬度，Lat2 Lung2 表示B点经纬度；

2. a=Lat1 – Lat2 为两点纬度之差  b=Lung1 -Lung2 为两点经度之差；

3. 6378.137为地球半径，单位为千米；

计算出来的结果单位为千米。

从google maps的脚本里扒了段代码，是用来计算两点间经纬度距离

1 private const double EARTH_RADIUS = 6378.137;2 private static double rad(doubled)3 {4 return d * Math.PI / 180.0;5 }6 public static double GetDistance(double lat1, double lng1, double lat2, doublelng2)7 {8 double radLat1 =rad(lat1);9 double radLat2 =rad(lat2);10 double a = radLat1 -radLat2;11 double b = rad(lng1) -rad(lng2);12 double s = 2 * Math.asin(Math.sqrt(Math.row(Math.sin(a/2),2) +

13 Math.cos(radLat1)*Math.cos(radLat2)*Math.pow(Math.sin(b/2),2)));14 s = s *EARTH_RADIUS;15 s = Math.round(s * 10000) / 10000;16 returns;17 }

在Android中可采用如下代码获取距离

1 public double getDistance(double lat1, double lon1, double lat2, doublelon2)2 {3 float[] results=new float[1];4 Location.distanceBetween(lat1, lon1, lat2, lon2, results);5 return results[0];6 }

在其他设备若没有类似android的Location的distanceBetween方法开采用如下代码获取

这个计算得出的结果是英里,如果要转换成公里,需要乘以1.609344,若是海里需要乘以0.8684

1 double distance(double lat1, double lon1, double lat2, doublelon2)2 {3 double theta = lon1 -lon2;4 double dist = Math.sin(deg2rad(lat1)) *Math.sin(deg2rad(lat2))5 + Math.cos(deg2rad(lat1)) *Math.cos(deg2rad(lat2))6 *Math.cos(deg2rad(theta));7 dist =Math.acos(dist);8 dist =rad2deg(dist);9 double miles = dist * 60 * 1.1515;10 returnmiles;11 }12 //将角度转换为弧度

13 static double deg2rad(doubledegree)14 {15 return degree / 180 *Math.PI;16 }17 //将弧度转换为角度

18 static double rad2deg(doubleradian)19 {20 return radian * 180 /Math.PI;21 }

计算方位角pab，

其中lat_a, lng_a是A的纬度和经度； lat_b, lng_b是B的纬度和经度。代码如下：

1 private double gps2d(double lat_a, double lng_a, double lat_b, doublelng_b)2 {3 double d = 0;4 lat_a=lat_a*Math.PI/180;5 lng_a=lng_a*Math.PI/180;6 lat_b=lat_b*Math.PI/180;7 lng_b=lng_b*Math.PI/180;8

9 d=Math.sin(lat_a)*Math.sin(lat_b)+Math.cos(lat_a)*Math.cos(lat_b)* Math.cos(lng_b-lng_a);10 d=Math.sqrt(1-d*d);11 d=Math.cos(lat_b)*Math.sin(lng_b-lng_a)/d;12 d=Math.asin(d)*180/Math.PI;13 //d = Math.round(d*10000);

14 returnd;15 }

地图应用接口中使用的经纬度单位是NTU，您需要把其它单位的经纬度值转化过来，以下为常用的转换关系：

基本转换：

NTU = 度*100000

度 = NTU/100000

例如：

经度 = 116.21345°

纬度 = 39.445875°

NTU经度：116.21345 * 100000 = 11621345(NTU)

NTU纬度：39.445874*100000=3944587(NTU)

度分转换：

将度分单位数据转换为度单位数据

度=度+分/60

例如：

经度 = 116°20.12’

纬度 = 39°12.34’

经度 = 116 + 20.12 / 60 = 116.33533°

纬度 = 39 + 12.34 / 60 = 39.20567°

NTU经度 = 116.33533 * 100000 = 11633533(NTU)

NTU纬度 = 39.20567 * 100000 = 3920567(NTU)

度分秒转换：

将度分秒单位数据转换为度单位数据

度 = 度 + 分 / 60 + 秒 / 60 / 60

例如：

经度 = 116°20’43”

纬度 = 39°12’37”

经度 = 116 + 20 / 60 + 43 / 60 / 60 = 116.34528°

纬度 = 39 + 12 / 60 + 37 / 60 / 60 = 39.21028°

NTU经度 = 116.34528 * 100000 = 11634528(NTU)

NTU纬度 = 39.21028 * 100000 = 3921028(NTU)

其格式的经纬值先转换度，再采用上面的运算。

与经纬度有关的基础知识：

NTU 相当于十万分之一度。

如何计算经度方向距离

比如: 点A的纬度为3995400，点B的纬度为3995300，则这两个点在沿着纬线的方向相差100米。

点A的经度为11695400，点B的经度为11695300，则这两个点在沿着经线的方向相差大约为77米,

该值的计算方法可以用近似公式：

经度方向距离 = 经度差 * cos(纬度值) = 100 * cos(39) = 77米。

如何计算两点间距离

比如：点A的经度为11695400，纬度为3995400。点B的经度为11695300，纬度为3995300。

公式:两点间距离 = [ (A点经度 - B点经度)^2 + (A点纬度 - B点纬度)^2 ] ^ (1/2) = [ (11695400 - 11695300)^2

+ (3995400 - 3995300)^2 ] ^(1/2) =(10000+10000) ^ (1/2) =141米