Запрос MySQL для координат широты и долготы, которые находятся в пределах данного радиуса мили

Question

В настоящее время у меня есть таблица MySQL, которая имеет следующую структуру:

id         name     lon            lat  
-----      -----    -----------    -----------
1          Mark     -76.316528     40.036027
2          John     -95.995102     41.25716
3          Paul     -82.337036     29.645095
4          Dave     -82.337036     29.645095
5          Chris    -76.316528     40.036027

Делая это, я в настоящее время запрашиваю базу данных, чтобы увидеть, находится ли местоположение пользователя в радиусе определенной мили от заданного местоположения.

function Haversine($lat_from, $lon_from, $lat_to, $lon_to) {
    $radius = 6371000;
    $delta_lat = deg2rad($lat_to-$lat_from);
    $delta_lon = deg2rad($lon_to-$lon_from);

    $a = sin($delta_lat/2) * sin($delta_lat/2) +
        cos(deg2rad($lat_from)) * cos(deg2rad($lat_to)) *
        sin($delta_lon/2) * sin($delta_lon/2);
    $c = 2*atan2(sqrt($a), sqrt(1-$a));

    // Convert the distance from meters to miles
    return ceil(($radius*$c)*0.000621371);
}

// Set the given location to NYC
$my_lon = -73.9844;
$my_lat = 40.7590;

// Query the DB for all of the users
$sql = "SELECT * FROM users";
$result = mysqli_query($con, $sql)or die(mysqli_error($con));
$count = mysqli_num_rows($result);
$i = 0;

while($row = mysqli_fetch_assoc($result)) {
    $lon[$i] = $row['lon'];
    $lat[$i] = $row['lat'];

    $i++;
}

for($i=0;$i<$count;$i++) {
    // Calculate the distance between each user's location and my location
    $distance = Haversine($my_lat, $my_lon, $lat[$i], $lon[$i);

    if($distance < 50) {
        echo "Close enough";
    }
}

Это хорошо работает только с несколькими сотнями строк в таблице. Но теперь, когда у меня есть десятки тысяч строк, проверка такого количества строк оказалась очень трудоемкой. Мне интересно, есть ли способ использовать формулу Haversine только для запроса строк, которые находятся в радиусе 50 миль.

Answer 1

Сферический закон формулы косинусов
(37 и -122 - широта и долгота вашего центра радиуса)

SELECT id, ( 3959 * acos( cos( radians(37) ) * cos( radians( lat ) ) 
    * cos( radians( long ) - radians(-122) ) + sin( radians(37) ) * sin(radians(lat)) ) ) AS distance 
FROM myTable
HAVING distance < 50
ORDER BY distance 

Характеристики

• Самый быстрый
• Точность, аналогичная формуле Harvesine

Формула Haversine

SELECT id, 3956 * 2 * ASIN(SQRT(POWER(SIN((37 - abs(lat)) * pi()/180 / 2), 2)
       + COS(37 * pi()/180 ) * COS(abs(lat) * pi()/180)
       * POWER(SIN((-122 – long) * pi()/180 / 2), 2) )) as  distance
FROM myTable
HAVING distance < 50
ORDER BY distance

Характеристики

• Быстро
• Более устойчивы к ошибкам с плавающей запятой

Обратите внимание, что 3959 - это радиус Земли в милях. Радиус Земли в км: 6371

Вы можете найти больше информации здесь