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($calculate);
}

附一个在线计算两点距离的地址：

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