php 根據兩點的經緯度計算距離

最近用到了geohash（參考《geohash：用字符串實現附近地點搜索》）來實現附近搜索的功能，在《微信、陌陌架構方案分析》文中提供了一個php的geohash類直接能對經緯地址編碼成字符串,  而其中包含一個getDistance方法獲取兩點間的距離，但在sphinx中使用過濾附近地標的方法SetGeoAnchor發現返回的距離值有誤，遂google出如下方法，經測試可用。。

/**
*  @desc 根據兩點間的經緯度計算距離
*  @param float $lat 緯度值
*  @param float $lng 經度值
*/
 function getDistance($lat1, $lng1, $lat2, $lng2)
 {
     $earthRadius = 6367000; //approximate radius of earth in meters

     /*
       Convert these degrees to radians
       to work with the formula
     */

     $lat1 = ($lat1 * pi() ) / 180;
     $lng1 = ($lng1 * pi() ) / 180;

     $lat2 = ($lat2 * pi() ) / 180;
     $lng2 = ($lng2 * pi() ) / 180;

     /*
       Using the
       Haversine formula

       http://en.wikipedia.org/wiki/Haversine_formula

       calculate the distance
     */

     $calcLongitude = $lng2 - $lng1;
     $calcLatitude = $lat2 - $lat1;
     $stepOne = pow(sin($calcLatitude / 2), 2) + cos($lat1) * cos($lat2) * pow(sin($calcLongitude / 2), 2);  $stepTwo = 2 * asin(min(1, sqrt($stepOne)));
     $calculatedDistance = $earthRadius * $stepTwo;

     return round($calculatedDistance);
 }

附一個在線計算兩點距離的地址：

http://www.nhc.noaa.gov/gccalc.shtml