9

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

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

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

我正在建立这个答案的第一部分。

上述答案的核心是:

直线10公里是:

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

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

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

我的班级在那里挂着,似乎工作了 10 公里。但是,我无法使其与任何其他距离一起使用。

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

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

这没有用。

你能解释一下为什么这种方法不适用于 5 公里(以及我尝试过的其他距离)吗?

下面是我的课程代码。显然开关是垃圾。任何帮助表示赞赏。

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;
    } 
}
4

2 回答 2

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);
}

公式取自这个问题。

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

$lngDeg = $distanceInKm/6371 * 180 / M_PI
于 2012-07-20T18:51:04.677 回答
0

正如另一个答案指出的那样,您真的想使用考虑到用户纬度的 Haversine 公式。

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

这个另一个答案有大量的资源可以在 PHP 或 SQL 中完成:

MySQL 大圆距离(Haversine 公式)

于 2012-07-20T19:06:35.507 回答