根据经纬度算出距离

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

java主要代码如下：

public static double getDistance()
    {  
        for(int i=0;i<list.size()-1;i++){
            BDLocation p1 = list.get(i);
            BDLocation p2 = list.get(i+1);
            double lat1 = (Math.PI/180)*(p1.getLatitude());
            double lon1 = (Math.PI/180)*(p1.getLongitude());
            double lat2 = (Math.PI/180)*(p2.getLatitude());
            double lon2 = (Math.PI/180)*(p2.getLongitude());
            double R = 6371.004;//地球半径，单位KM
            distance += Math.acos(Math.sin(lat1)*Math.sin(lat2)+Math.cos(lat1)*Math.cos(lat2)*Math.cos(lon2-lon1))*R*1000;
        } 
        return distance;            //单位：m
    }
   
    public static void main(String[] args) {
        list=new ArrayList<BDLocation>();    
        BDLocation bd1=new BDLocation(38.891819,115.521052);
        BDLocation bd2=new BDLocation(38.891901,115.521109);
        BDLocation bd3=new BDLocation(38.891994,115.521194);
        list.add(bd1);
        list.add(bd2);
        list.add(bd3);
        System.out.print(getDistance());
    }

完整代码如下：请点击这里

参考：http://www.cnblogs.com/ycsfwhh/archive/2010/12/20/1911232.html