GPS两坐标点距离计算

工作中测试百度LBS定位接口的准确性，需要计算本公司的GPS点和调用百度接口得到的GPS坐标点的误差，从网上找了很多资料，发现下面的算法比较准确（测试坐标系wgs84）。

转自：
http://www.open-open.com/lib/view/open1388553235141.html

/**
*计算gps两坐标点的距离（单位：M）
**/
public static double getDistanceFromXtoY(double lat_a, double lng_a, double lat_b, double lng_b){
            double pk = (double)(180 / 3.14169);

            double a1 = lat_a / pk;
            double a2 = lng_a / pk;
            double b1 = lat_b / pk;
            double b2 = lng_b / pk;

            double t1 = Math.Cos(a1) * Math.Cos(a2) * Math.Cos(b1) * Math.Cos(b2);
            double t2 = Math.Cos(a1) * Math.Sin(a2) * Math.Cos(b1) * Math.Sin(b2);
            double t3 = Math.Sin(a1) * Math.Sin(b1);
            double tt = Math.Acos(t1 + t2 + t3);

            return 6366000 * tt;
        }

测试数据：
公司gps坐标点：38.9631423611111,121.585119166667
百度接口坐标点：38.963727756709325,121.58595020536029
坐标点距离:0.09687790472161027KM

公司gps坐标点：31.4088551666667,116.330832333333
百度接口坐标点：31.408822041273368,116.33178209779567
坐标点距离:0.09014046635371215

顺便提供一个gps好的工具网址：
http://www.gpsspg.com/maps.htm