Como encontrar a cidade mais próxima de um PI em um banco de dados OpenStreetMap / PostGIS?

9

Estou interessado em encontrar dados adicionais para os POIs no banco de dados OSM.

Atualmente, posso usar o seguinte SQL para obter o nome + long + lat para todos os POIs

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'

Também gostaria de ter o nome da cidade / cidade mais próxima ou da cidade / local em que o PI está situado.

postgis openstreetmap point-of-interest Catalin
fonte
Você já tem um sistema de geocodificação reversa (local para nomear) configurado?
BradHards

Respostas:

0

Resposta retirada de uma pergunta relacionada :

Muito fácil, assumindo que você tenha o lat / lngs.

  1. Converter lat / lng no osm data srid
  2. Calcule a distância de todos os PIs de seu interesse no banco de dados.
  3. Escolha a linha com a menor distância.
  4. Além disso, você pode limitar sua pesquisa a um raio de interesse (por exemplo, 5 KM ..)

Como exemplo, vamos encontrar os 5 postos de gasolina mais próximos perto de nosso PI.

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 de resultados

-[ 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
fonte