计算球面上经纬度坐标方法比较

计算球面上的两点(坐标为经纬度)之间的距离可以直接通过公式计算得到，也可以先将经纬度坐标转化为墨卡托投影坐标来,然后用平面中两点之间的距离公式来计算。

在网上找了一些代码，然后简单进行了测试，发现前者精度更高：

 资料来源：http://0414.iteye.com/blog/2039199

　　　　　http://blog.sina.com.cn/s/blog_8ab3b2cf0100xd69.html

package com.suncreate.spatialquery.web.utils;

public class LatitudeLontitudeUtil2 {  

    // http://blog.charlee.li/location-search/  

    /** 地球半径 */  
    private static final double EARTH_RADIUS = 6371393;  

    static double M_PI = Math.PI;

    public LatitudeLontitudeUtil2() {  

    }  

    //经纬度转墨卡托
    public static double[] lonLat2Mercator(double lon,double lat)
    {
        double[] xy = new double[2];
        double x = lon *20037508.342789/180;
        double y = Math.log(Math.tan((90+lat)*M_PI/360))/(M_PI/180);
        y = y *20037508.34789/180;
        xy[0] = x;
        xy[1] = y;
        return xy;
    }
    
    //点到点距离算法一
    public static double getDistance2(double lat0, double lng0, double lat1, double lng1) {  
        double[] xy = new double[2];
        
        System.out.println("----原始经纬度坐标----");
        System.out.println("lat0:"+lat0+",lng0:"+lng0);
        System.out.println("lat1:"+lat1+",lng1:"+lng1);
        
        xy =LatitudeLontitudeUtil2.lonLat2Mercator(lng0, lat0);
        double lat0m= xy[1];
        double lng0m= xy[0];
        
        xy =LatitudeLontitudeUtil2.lonLat2Mercator(lng1, lat1);
        double lat1m= xy[1];
        double lng1m= xy[0];
        
        System.out.println("----转化后墨卡托坐标----");
        System.out.println("lat0m:"+lat0m+",lng0m:"+lng0m);
        System.out.println("lat1m:"+lat1m+",lng1m:"+lng1m);
         
        double distance = 0; 
        distance= Math.sqrt((lat0m - lat1m) * (lat0m - lat1m) + (lng0m - lng1m) * (lng0m - lng1m));
        return distance;  
    }  
    //点到点距离算法二
    public static double getDistance(double lat0, double lng0, double lat1, double lng1) {  

        lat0 = Math.toRadians(lat0);  
        lat1 = Math.toRadians(lat1);  
        lng0 = Math.toRadians(lng0);  
        lng1 = Math.toRadians(lng1);  

        double dlng = Math.abs(lng0 - lng1);  
        double dlat = Math.abs(lat0 - lat1);  
        double h = hav(dlat) + Math.cos(lat0) * Math.cos(lat1) * hav(dlng);  
        double distance = 2 * EARTH_RADIUS * Math.asin(Math.sqrt(h));  

        return distance;  
    }  
    
    public static double hav(double theta) {  
        double s = Math.sin(theta / 2);  
        return s * s;  
    }  
    
    public static void main(String[] args) {  
        //测试用例
        //117.222009,31.815101   祁门路天鹅湖左
        //117.238718,31.815132   祁门路天鹅湖右
        
        double lat = 31.815101;  
        double lng = 117.222009;  
        
        double lat1 = 31.815132 ;  
        double lng1 = 117.238718;  
        
        double d = LatitudeLontitudeUtil2.getDistance(lat, lng, lat1, lng1);  
        double d2 = LatitudeLontitudeUtil2.getDistance2(lat, lng, lat1, lng1);
        
        System.out.println(d); 
        System.out.println(d2); 
        
        //百度地图测距约为1.6km,可知直接计算球面距离较为准确，而直接先将经纬度转化为墨卡托坐标再求距离不准
    }  
}