Как найти все точки в заданной координате с определенным радиусом?

Question

[Владимир Шишминцев]

Добрый день.

Подскажите как найти все точки по базе, зная координаты и радиус.

CREATE TABLE IF NOT EXISTS `points` (
  `point` int(11) DEFAULT NULL,
  `n` double DEFAULT NULL,
  `e` double DEFAULT NULL,
  KEY `n` (`n`),
  KEY `e` (`e`),
) ENGINE=MyISAM DEFAULT CHARSET=utf8;


INSERT INTO `points` (`point`, `n`, `e`) VALUES
(1, 55.75272901, 37.60548887),
(2, 55.75473895, 37.61351058),
(3, 55.75271691, 37.60612187),
(4, 55.75280693, 37.61898127);

То есть приходит ко мне координата 55.75272901:37.60548887 с радиусом 100 метров, и согласно её мне надо получить запись 1 и 3.

Answer 1

[egor_nullptr]

Для начала создать хранимую процедуру вычисления расстояния по координатам:

delimiter $$

CREATE FUNCTION GetDistanceBetweenGeoLocations
(
lat1 float,
long1 float,
lat2 float,
long2 float
)
RETURNS float
DETERMINISTIC

BEGIN

DECLARE distance float;
DECLARE earth_radius float;
DECLARE axis float;

SET earth_radius = 6371.009;

SET axis = (SIN(RADIANS(lat2-lat1)/2) * SIN(RADIANS(lat2-lat1)/2) +
COS(RADIANS(lat1)) * COS(RADIANS(lat2)) *
SIN(RADIANS(long2-long1)/2) * SIN(RADIANS(long2-long1)/2));

SET distance = earth_radius * (2 * ATN2(SQRT(axis), SQRT(1-axis)));

RETURN ROUND(distance, 3);

END;
$$

После этого можно искать:

select * from points
where GetDistanceBetweenGeoLocations(n, 55.75272901, e, 37.60548887) <= 0.1

Comments

[Владимир Шишминцев]

@egor_nullptr что за значение "0.1" ?
что-то не работает(

[egor_nullptr]

0.1 - это 100 метров в километрах

[Владимир Шишминцев]

@egor_nullptr mysql ругался на ATN2, я заменил на ATAN2 (по мануалу это аналог). Но результата нету - возвращает 0 записей

[egor_nullptr]

в теле функции замените вызов ATN2 на ATAN2,

$ select *, GetDistanceBetweenGeoLocations(n, 55.75272901, e, 37.60548887) as d from points;
+-------+-------------+-------------+-------------------+
| point | n           | e           | d                 |
+-------+-------------+-------------+-------------------+
|     1 | 55.75272901 | 37.60548887 |     2430.95703125 |
|     2 | 55.75473895 | 37.61351058 |    2430.337890625 |
|     3 | 55.75271691 | 37.60612187 |  2430.89404296875 |
|     4 | 55.75280693 | 37.61898127 | 2429.639892578125 |
+-------+-------------+-------------+-------------------+
4 rows in set (0.00 sec)

Расстояния до точек в последнем столбце, ~2.4км до каждой

[egor_nullptr]

ага вот оно, параметры в запросе местами надо немного поменять

select *, GetDistanceBetweenGeoLocations(n, e, 55.75272901, 37.60548887) as d from points;
+-------+-------------+-------------+---------------------+
| point | n           | e           | d                   |
+-------+-------------+-------------+---------------------+
|     1 | 55.75272901 | 37.60548887 |                   0 |
|     2 | 55.75473895 | 37.61351058 |  0.5495193600654602 |
|     3 | 55.75271691 | 37.60612187 | 0.03964652866125107 |
|     4 | 55.75280693 | 37.61898127 |  0.8443683385848999 |
+-------+-------------+-------------+---------------------+
4 rows in set (0.01 sec)

[Владимир Шишминцев]

@egor_nullptr благодарю за помощь

[mind3]

Норм решение через функцию в sql запросе. Но что делать если вдруг записей точек гораздо больше? Мне кажется такой запрос будет медленным для обработки большого количества точек.

Answer 2

[Spetros]

В школьном курсе геометрии есть формула описывающая окружность с центром в x0 и y0:

Если из этого равенства сделать неравенство, то мы можем получить формулу описывающую все точки находящиеся внутри нее, для которых это неравенство справедливо.

Comments

[Rsa97]

Это из планиметрии, для географии формулы другие.

Answer 3

[Sergei Borisov]

Тригонометрия в хранимых процедурах, конечно решение. Но вообще-то поддержка geo это большая тема, которая решается в СУБД отдельнвыми инструментами. Например, для mysql вам сюда: dev.mysql.com/doc/refman/5.0/en/gis-introduction.html