Получить ближайший путь при задании (широта, долгота) с помощью OSM в Postgres с Postgis с помощью SQL-запроса - PullRequest
Новая
       35

Получить ближайший путь при задании (широта, долгота) с помощью OSM в Postgres с Postgis с помощью SQL-запроса

0 голосов
/ 05 сентября 2018

Учитывая набор (широта, длина), я пытаюсь найти максимальную скорость, используя "max_speed", а тип улицы - "шоссе".

Я загрузил свою базу данных (Postgres и Postgis) следующим образом: 

$ osm2pgsql -c -d gis --slim -C 50000  /var/lib/postgresql/data/germany-latest.osm.pbf


Самый близкий связанный вопрос, который я смог найти, был Как запросить все магазины вокруг определенной долготы / широты с помощью osm-postgis? . Я взял запрос и подключил (широта, долгота), который я нашел в Google Maps для центра города Мюнхен (так как пост был также связан с центром города Мюнхен, и у меня есть карта для Германии). Результат оказывается пустым. 

gis=# SELECT name, shop FROM planet_osm_point WHERE ST_DWithin(way ,ST_SetSrid(ST_Point(48.137969, 11.573829), 900913), 100);
 name | shop 
------+------
(0 rows)

Кроме того, когда я смотрю на planet_osm_nodes, который содержит пары (lat, long) напрямую, я не получаю результатов: 

gis=# SELECT * FROM planet_osm_nodes WHERE ((lat BETWEEN 470000000 AND 490000000) AND (lon BETWEEN 100000000 AND 120000000)) LIMIT 10;
 id | lat | lon | tags 
----+-----+-----+------
(0 rows)

Я проверил данные в моей базе данных: 

gis=# SELECT COUNT(*) FROM planet_osm_point;
  count  
---------
 9924531
(1 row)

и 

gis=# SELECT COUNT(*) FROM planet_osm_nodes;
   count   
-----------
 288597897
(1 row)

Так что в идеале мой вопрос будет
Q: Как мне найти "максимальную скорость" и "шоссе" по заданному набору (широта, долгота)
В качестве альтернативы, мои вопросы:
Q: Как мне заставить работать запрос из другого поста переполнения стека?

Мое лучшее предположение заключается в том, что мне нужно каким-то образом преобразовать мои (широта, долгота) или что по какой-то причине у меня просто неверные данные.

Редактировать: добавленные данные образца по запросу: 

gis=# SELECT * FROM planet_osm_point LIMIT 1;
  osm_id   | access | addr:housename | addr:housenumber | addr:interpolation | admin_level | aerialway | aeroway | amenity | area | barrier | bicycle | brand | bridge | boundary | building | capital | construction | covered | culvert |
 cutting | denomination | disused | ele | embankment | foot | generator:source | harbour | highway  | historic | horse | intermittent | junction | landuse | layer | leisure | lock | man_made | military | motorcar | name | natural | off
ice | oneway | operator | place | poi | population | power | power_source | public_transport | railway | ref | religion | route | service | shop | sport | surface | toll | tourism | tower:type | tunnel | water | waterway | wetland | wi
dth | wood | z_order |                        way                         
-----------+--------+----------------+------------------+--------------------+-------------+-----------+---------+---------+------+---------+---------+-------+--------+----------+----------+---------+--------------+---------+---------+
---------+--------------+---------+-----+------------+------+------------------+---------+----------+----------+-------+--------------+----------+---------+-------+---------+------+----------+----------+----------+------+---------+----
----+--------+----------+-------+-----+------------+-------+--------------+------------------+---------+-----+----------+-------+---------+------+-------+---------+------+---------+------------+--------+-------+----------+---------+---
----+------+---------+----------------------------------------------------
 304070863 |        |                |                  |                    |             |           |         |         |      |         |         |       |        |          |          |         |              |         |         |
         |              |         |     |            |      |                  |         | crossing |          |       |              |          |         |       |         |      |          |          |          |      |         |    
    |        |          |       |     |            |       |              |                  |         |     |          |       |         |      |       |         |      |         |            |        |       |          |         |   
|      |         | 010100002031BF0D0048E17A94F19F2941CDCCCCDCC60D5741
(1 row)

и 

gis=# SELECT * FROM planet_osm_nodes LIMIT 1;
   id   |    lat    |   lon    | tags 
--------+-----------+----------+------
 234100 | 666501948 | 80442755 | 
(1 row)

Редактировать 2: было упоминание о "SRID", поэтому я добавил пример данных из другой таблицы: 

gis=# SELECT * FROM spatial_ref_sys LIMIT 1;
 srid | auth_name | auth_srid |                                                                                                                                                                    srtext                                  
                                                                                                                                  |                                           proj4text                                            
------+-----------+-----------+------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
----------------------------------------------------------------------------------------------------------------------------------+------------------------------------------------------------------------------------------------
 3819 | EPSG      |      3819 |  GEOGCS["HD1909",DATUM["Hungarian_Datum_1909",SPHEROID["Bessel 1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],TOWGS84[595.48,121.69,515.35,4.115,-2.9383,0.853,-3.408],AUTHORITY["EPSG","1024"]],PR
IMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","3819"]] | +proj=longlat +ellps=bessel +towgs84=595.48,121.69,515.35,4.115,-2.9383,0.853,-3.408 +no_defs 
(1 row)

1 Ответ

0 голосов
/ 05 сентября 2018

Геометрия в PostGIS имеет другой порядок (широта долготы), сначала идет долгота, а не широта. Также, если вы хотите преобразовать точку из одного SRID в другой, используйте st_transfrom (), а не ST_SetSrid. ST_Transform обеспечивает точное преобразование ваших данных из одной системы координат в другую. 

select st_astext(st_transform(ST_SetSrid(ST_Point(11.573829,48.137969), 4326),900913))

ST_SetSrid - просто измените SRID для объекта. 

select st_astext((ST_SetSrid(ST_Point(11.573829,48.137969),900913)

Итак, вы должны изменить свой SQL таким образом 

SELECT name, shop 
  FROM planet_osm_point
 WHERE ST_DWithin(way,st_transform(ST_SetSrid(ST_Point(11.573829,48.137969), 4326),900913), 100);
...