Как найти ближайшую улицу к GPS точке?

Question

[beduin01]

Есть GPS координаты. Нужно найти ближайшую улицу к каждой GPS точке. Нашел пример подобного запроса:

SELECT p.id AS point_id, b.id AS road_id, ST_DISTANCE(b.geom, p.geom) AS distance
FROM points p
CROSS JOIN LATERAL (
    SELECT r.id AS id, r.geom AS geom
    FROM roads r
    ORDER BY r.geom <-> p.geom 
    LIMIT 1
) b;

Никак не могу понять как мне его использовать. К примеру у меня есть координаты:
lon, lat 37.72308, 55.47957

Мне их в `ST_DISTANCE` вписать? Если да то в какой форме?

Answer 1

[freeExec]

Да, прямо туда ST_POINT.
Но вообще странный запрос. Обычно строят буфер в 100, 200 и т.д. метров, и смотрят какая дорога в него попавшая ближе.
А можно и без буфера если нужно только название:

SELECT osm_id, highway, name, tags, ST_DISTANCE(way, point) AS Distance
  FROM public.planet_osm_line
	INNER JOIN (SELECT ST_TRANSFORM( ST_SetSRID(ST_POINT(48.323, 54.263), 4326), 900913) AS point ) AS p
		ON ST_DWithin(point, way, 200)
  WHERE highway IS NOT NULL AND name IS NOT NULL  
  ORDER BY Distance
  LIMIT 1;

Comments

[beduin01]

А не могли бы пример запроса показать? не получается что-то написать

[beduin01]

Если можно то лучше естественно с буфером.

[freeExec]

beduin01: Добавил пример.

[beduin01]

freeExec: спасибо, но у вас видимо структура БД какая-то другая т.к. у меня допустим `way` нет. Я выкинул разные переводы из проекции в проекцию и в итоге у меня вылетает ошибка: "ERROR: function st_distance(geometry) does not exist". Если не сложно можете дать какой-то другой пример.

Вот так моя БД выглядит: https://snag.gy/s8F9AC.jpg

[freeExec]

beduin01: Потому что есть только st_distance(geometry, geometry)

[beduin01]

freeExec: так, а не подскажите на что мне way заменить согласно моим данным? на geometry?

[freeExec]

beduin01: Если у вас геометрия лежит в geometry, то да.

[beduin01]

Не получается :( находят дорогу не там где она должна быть :(

[freeExec]

beduin01: У вас координаты lon, lat передаются в точку?
А так же нужно учитывать, что результат distance идёт в единицах проекции. Она у вас в метрах, скорее всего нет.

[beduin01]

В плане передаются в точку? У меня набор GPS координат.
Проекция же у OSM 4326.
Как это исправить можно?

[freeExec]

Вот 37.72 это у вас широта или долгота?
Да, сырые данные ОСМ это 4326(lon\lat). Перевод в другую проекцию ST_TRANSFORM.
Можете посчитать дистанцию через ST_Distance_Sphere она возвращает в метрах.

[beduin01]

Все разобрался. Мой косяк был. К сожалению оказалось, что ST_DISTANCE не позвоялет использовать пространственный индекс и нужно использовать ST_DWithin. Не могли бы подсказать как с помощью него запрос делать?