1         /**
 2     　　  * 将角度转化为弧度
 3     　　*/
 4         public static double radians(double d)
 5         {
 6             return d * Math.PI / 180.0;
 7         }
 8         /**
 9           * 根据两点经纬度坐标计算直线距离
10           * <p>
11           * S = 2arcsin√sin²(a/2)+cos(lat1)*cos(lat2)*sin²(b/2)￣*6378.137
12           * <p>
13           * 1. lng1 lat1 表示A点经纬度，lng2 lat2 表示B点经纬度；<br>
14           * 2. a=lat1 – lat2 为两点纬度之差  b=lng1 -lng2 为两点经度之差；<br>
15           * 3. 6378.137为地球赤道半径，单位为千米；
16           *
17           * @param lng1 点1经度
18           * @param lat1 点1纬度
19           * @param lng2 点2经度
20           * @param lat2 点2纬度
21           * @return 距离，单位千米(KM)
22           * @see <a href="https://zh.wikipedia.org/wiki/%E5%8D%8A%E6%AD%A3%E7%9F%A2%E5%85%AC%E5%BC%8F">半正矢(Haversine)公式</a>
23           */
24         public static double getDistanceFrom2LngLat(double lng1, double lat1, double lng2, double lat2)
25         {
26             //将角度转化为弧度
27             double radLng1 = radians(lng1);
28             double radLat1 = radians(lat1);
29             double radLng2 = radians(lng2);
30             double radLat2 = radians(lat2);
31 
32             double a = radLat1 - radLat2;
33             double b = radLng1 - radLng2;
34 
35             return 2 * Math.Asin(Math.Sqrt(Math.Sin(a / 2) * Math.Sin(a / 2) + Math.Cos(radLat1) * Math.Cos(radLat2) * Math.Sin (b / 2) * Math.Sin (b / 2))) * 6378.137;
36         }

 

 

 1 SELECT
 2     *,
 3     6378.138 * 2 * ASIN(
 4       SQRT(
 5         POW(
 6           SIN(
 7             (
 8               '.$lat.' * PI() / 180 - lat * PI() / 180
 9             ) / 2
10           ), 2
11         ) + COS('.$lat.' * PI() / 180) * COS(lat * PI() / 180) * POW(
12           SIN(
13             (
14               '.$lng.' * PI() / 180 - lng * PI() / 180
15             ) / 2
16           ), 2
17         )
18       )
19     ) *1000 AS distance
20 FROM
21     distance
22 ORDER BY
23     distance ASC