php通过经纬度计算距离及其逆运算

有朋友遇到这样一个问题，自己数据库里面有一些商家，需要根据用户提交的地理位置找出一定范围内的商家来返回个用户。解决这个问题有两种思路。

第一种思路，从数据库提取出商家的经纬度，然后计算和用户经纬度间的距离，留下在距离范围内地。

这种方法貌似很科学，不过数据量稍微大一点就完蛋了。这里涉及到的是计算两个经纬度间的距离，方法的代码如下：

/**
 * [getdistance3 description]
 * @param $lng1 [经度1]
 * @param $lat1 [纬度1]
 * @param $lng2 [经度2]
 * @param $lat2 [纬度2]
 * @return $distance [距离（米）]
 */
function getdistance3($lng1, $lat1, $lng2, $lat2){
    $radLat1 = deg2rad($lat1);
    $radLat2 = deg2rad($lat2);
    $radLng1 = deg2rad($lng1);
    $radLng2 = deg2rad($lng2);
    $a = $radLat1 - $radLat2;
    $b = $radLng1 - $radLng2;
    $distance = 2*asin(sqrt(pow(sin($a/2),2)+cos($radLat1)*cos($radLat2)*pow(sin($b/2),2)))*6378137;
    return $distance;
}

这个算法的公式是Google map用的，6378137(m)是地球半径，如果觉得其他值更合适可以改。还有另外一个公式，写成代码替换上面的$distance是这样：

$distance = acos(sin($radLat1)*sin($radLat2)+cos($radLat1)*cos($radLat2)*cos($b))*637813;

通过三角恒等变换可以知道二者是恒等的。

这个方法用来算两点距离是没问题，不过要用来解决上面的给出特定范围内商家的问题不现实。

第二种思路，通过给出的点A的经纬度坐标和一个确定的单位距离，先算出距A单位距离的经纬度范围，再通过这个范围为条件从数据库获取满足条件的商家数据。

程序代码如下：

/**
 * [getRange description]
 * @param  $lng    [经度]
 * @param  $lat    [纬度]
 * @param  $raidus [半径(米)]
 * @return [array] $range  [经纬度范围]
 */
function getRange($lng, $lat, $raidus){
    $PI = 3.14159265;
    $latitude = $lat;
    $longitude = $lng;
    $degree = 40075016/360;
    $raidusMetre = $raidus;
    $radiusLat = $raidusMetre / $degree;
    $minLat = $latitude - $radiusLat;
    $maxLat = $latitude + $radiusLat;
    $mpdLng = $degree * cos($latitude * ($PI/180));
    $radiusLng = $raidusMetre / $mpdLng;
    $minLng = $longitude - $radiusLng;
    $maxLng = $longitude + $radiusLng;
    $range = array('minLng' => $minLng, 'maxLng' => $maxLng, 'minLat' => $minLat, 'maxLat' => $maxLat);
    return $range;
}

这里40075016(m)是地球周长，如果觉得其他值更合适也可以改。通过这个方法，得出了距给出的点（$lng, $lat）5000米以内的坐标范围：经度[$minLng, $maxLng]，纬度[$minLat, $maxLat]，本来想要的是一个半径为5000米的圆形范围，现在得出的是一个正方形范围，作为近似使用应该可以了。
接下来就只用到数据库筛选出经纬度在这个范围的商家出来，按照需要的数据格式返回给用户了。如果给定的半径范围内没有找到商家，可以扩大半径来搜出商家。
到此算是给出这个问题一个近似解决了，先暂时给他用着。