无距离地理搜索

时间:2012-07-19 15:54:43

标签: php geolocation gis

这是我第一次看到地理位置。请原谅我的长话题。我希望很清楚。

我有一个数据库,其中包含lat&amp ;; LNG。我正在尝试编写一个PHP类,它定义了四个十进制坐标,形成一个边界框(然后检索框内的所有记录)。

我只需知道哪些位置在用户位置的 n km范围内。我不需要从用户到位置的距离。

我正在构建answer的第一部分。

以上答案的核心是:

  

直线10公里是:

on the latitude is equal to ~1'(minute)
on the longitude is equal to ~6'(minutes)
     

以此为基础,进行一些快速数学运算并在查询中添加   WHERE子句删除“盒子”之外的任何位置   通过添加缓冲区创建,假设为1'lat& 6'   长。

嗯,“做一些快速数学”帕特里克说。几天后我没能正确地完成数学计算。

我的班级在那里停留,似乎工作了10公里。但是,我没有让它与任何其他距离一起工作。

首先我移动到秒,所以10km是:

on the latitude is equal to 60 seconds
on the longitude is equal to 360 seconds

这很好用。

然而,我尝试了5km是:

on the latitude is equal to 30 seconds
on the longitude is equal to 180 seconds

这不起作用。

你能解释一下为什么这种方法在5km(以及我尝试的其他距离)下都不起作用吗?

以下是我的班级代码。显然这个开关是垃圾。任何帮助表示赞赏。

class coords 
{
/*
    Calculate a bounding box from a decimal coordinate and distance.
    The public vars are populated with minimum and maximum longitudes and latitudes which define the boundary.
*/
    public $dec_lat_min;
    public $dec_lat_max;
    public $dec_lng_min;
    public $dec_lng_max;

    function init($dec_lat, $dec_lng, $dist) 
    { 
        // This switch is a terrible way to allow multiple distances.

        // 10km = 1 min lat = 60 sec lat
        // 10km = 6 min lng = 360 sec lng

        // 5km = 30 sec lat
        // 5km = 180 sec lng

        // 1km = 6 sec lat
        // 1km = 36 sec lat

        // 500m = 3 sec lat
        // 500m = 18 sec lat

        switch($dist)
        {
            case 10: // 10km
                $sec_diff_lat = 60;
                $sec_diff_lng = 360;
                break;
            case 5: // 5km
                $sec_diff_lat = 30;
                $sec_diff_lng = 180;
                break;
            case 1: // 1km
                $sec_diff_lat = 6;
                $sec_diff_lng = 36;
                break;
            default: // 500m
                $sec_diff_lat = 3;
                $sec_diff_lng = 18;
                break;  
        }

        // Convert lat to DMS
        $dms_lat = $this->dec2dms($dec_lat);

        // Allow for western hemisphere (ie negative)
        $dms_lat['hem'] == '-' ? $h = -1 : $h = 1;

        // Populate min and max latitudes
        $this->dec_lat_min = $this->dms2dec($dms_lat['deg'],$dms_lat['min'],$dms_lat['sec']+(-1 * $sec_diff_lat * $h),$dms_lat['hem']);
        $this->dec_lat_max = $this->dms2dec($dms_lat['deg'],$dms_lat['min'],$dms_lat['sec']+($sec_diff_lat * $h),$dms_lat['hem']);

        $dms_lng = $this->dec2dms($dec_lng);

        $dms_lng['hem'] == '-' ? $h = -1 : $h = 1;

        $this->dec_lng_min = $this->dms2dec($dms_lng['deg'],$dms_lng['min'],$dms_lng['sec']+(-1 * $sec_diff_lng * $h),$dms_lng['hem']);
        $this->dec_lng_max = $this->dms2dec($dms_lng['deg'],$dms_lng['min'],$dms_lng['sec']+($sec_diff_lng * $h),$dms_lng['hem']);

    }

    function dec2dms($d) 
    {
        $d = (string)$d;

        // got dashes?
        if ($d[0] == "-") {
            $hem = '-';
            $dVal = substr($d,1);
        } else {
            $hem = '';
            $dVal = $d;
        }

        // degrees = degrees
        $dVals = explode('.', $dVal);
        $dmsDeg = $dVals[0];

        // * 60 = mins
        $dRemainder = ('0.'.$dVals[1]) * 60;
        $dmsMinVals = explode('.', $dRemainder);
        $dmsMin = $dmsMinVals[0];

        // * 60 again = secs
        $dMinRemainder = ('0.'.$dmsMinVals[1]) * 60;
        $dmsSec = round($dMinRemainder);

        return array("deg"=>$dmsDeg,"min"=>$dmsMin,"sec"=>$dmsSec, "hem"=>$hem);
    }

    function dms2dec($deg,$min,$sec,$hem) 
    {
        // find decimal latitude
        $d = $deg+((($min*60)+($sec))/3600);
        $d = $hem.$d;

        return $this->round100000($d);
    }

    function round100000($v) 
    {
        return round($v * 100000) / 100000;
    } 
}

2 个答案:

答案 0 :(得分:1)

  

直线10公里是:经度等于~6分钟

距离与经度之比不是常数 - 它取决于用户的纬度。您可以将此功能用于经度差异:

function distanceToLng($distance,$latDeg) {
    $latRad = $latDeg / 180 * M_PI;
    $R = 6371;
    return atan2(1*sin($distance/$R)*cos($latRad), cos($distance/$R) - sin($latRad)*sin($latRad)) * (180 / M_PI);
}

公式取自this问题。

您可以使用以下内容计算给定距离的纬度差:

$lngDeg = $distanceInKm/6371 * 180 / M_PI

答案 1 :(得分:0)

正如另一个答案所指出的,你真的想要使用考虑用户纬度的Haversine公式。

Haversine公式:http://en.wikipedia.org/wiki/Haversine_formula

这个其他答案有很多很好的资源可以在PHP或SQL中实现这一点:

MySQL Great Circle Distance (Haversine formula)