¿Cómo encontrar la ciudad más cercana para un PDI en una base de datos OpenStreetMap / PostGIS?

9

Estoy interesado en encontrar datos adicionales para los PDI en la base de datos OSM.

Actualmente puedo usar el siguiente SQL para obtener el nombre + long + lat para todos los POI

SELECT name, 
         x(transform(way, 4326)) AS lon, 
         y(transform(way, 4326)) AS lat
  FROM planet_osm_point 
  WHERE tourism='museum'
UNION
  SELECT name, 
         x(centroid(transform(way, 4326))) AS lon, 
         y(centroid(transform(way, 4326))) AS lat
  FROM planet_osm_polygon 
  WHERE tourism='museum'

También me gustaría tener el nombre de la ciudad / ciudad más cercana o la ciudad / ciudad en la que se encuentra el PDI.

postgis openstreetmap point-of-interest Catalin
fuente
¿Ya tiene configurado un sistema de geocodificación inversa (ubicación para nombrar)?
BradHards

Respuestas:

0

Respuesta extraída de una pregunta relacionada :

Bastante fácil suponiendo que tienes los lat / lngs.

  1. Convierta lat / lng a la cuadrícula de datos osm
  2. Calcule la distancia a todos los PDI que le interesan en la base de datos.
  3. Elija la fila con la menor distancia.
  4. Además, es posible que desee limitar su búsqueda a un radio de interés (digamos 5 KM ..)

Como ejemplo, busquemos las 5 estaciones de servicio más cercanas a nuestro PDI.

SELECT osm_id,
       name,
       brand,
OPERATOR,
       ST_Distance(barabara_point.way, st_transform( st_setsrid(st_makepoint(36.768676, -1.289927), 4326), 900913)) AS distance
FROM barabara_point
WHERE ST_DWithin(barabara_point.way, st_transform( st_setsrid(st_makepoint(36.768676, -1.289927), 4326), 900913), 5000)
  AND amenity = 'fuel'
ORDER BY distance;

Conjunto resultante

-[ RECORD 1 ]--------------
osm_id   | 1334386299
name     |
brand    | Kobil
operator |
distance | 244.544617525396
-[ RECORD 2 ]--------------
osm_id   | 915386940
name     | BP
brand    |
operator |
distance | 1012.76398108684
-[ RECORD 3 ]--------------
osm_id   | 1916188670
name     | Shell
brand    |
operator |
distance | 1160.48349020106
-[ RECORD 4 ]--------------
osm_id   | 558923017
name     | Kenol
brand    |
operator |
distance | 1310.22002620899
-[ RECORD 5 ]--------------
osm_id   | 30092081
name     | Caltex
brand    |
operator |
distance | 1419.61102322244
Lmwangi
fuente