Как узнать город из координаты lat long?

Question

[Oleg]

Доброго времени суток.
Имеется база openstreetmap одной страны. Как узнать город(не поселок, область, а именно город) из координаты lat long в postgres-postgis?

Answer 1

[Fixid]

Web api есть? Или только доступ к самой бд? Есть ли таблица placex?

Comments

[Oleg]

Api никакого нет(nominatim не пойдет), собираюсь писать свою прослойку. БД состоит из таблиц: line, nodes, point, polygon, rels, roads, ways с префиксом planet_osm.

[Oleg]

я могу узнать область по координате, а вот как сделать в одном запросе чтобы город выдавало - непонятно.
К примеру:
select name, admin_level,place from planet_osm_polygon
INNER JOIN (SELECT ST_TRANSFORM(ST_SetSRID(ST_POINT(51.521120,47.400039), 4326), 3857) AS point ) AS p
ON ST_Within(point, way)

Выдаст "Атырауская область", но города Атырау нет в результате. Если поставить координату 51.883669,47.108317, то можно сделать выборку по place = 'city', только вот для предыдущей координаты результат пустой.

[Fixid]

Oleg: тогда я вас очень очень огорчу (сам уже прошел этот путь)
Текущие данные в базе сырые и без индексации. Быстро и эффективно искать не получится.
Если хочется делать по нормальному, то придется запустить индексацию через nominatin/osm2pgsql, которая создаст новую структуру (старые таблицы можно будет удалить), где все объекты проиндексированы и легко найти что угодно и где угодно

[Fixid]

Oleg: не все города имеют place = 'city', большие города с полигоном обычно являются административными территориями с point, где place = 'city'.

т.е объект может быть комбинированным из административного полигона и связанного города точки

[Oleg]

Fixid: импорт в postgres делал именно через osm2pgsql. Теперь нужно делать индексацию? были ворнинги на индексацию при импорте, но их никак решить не удалось. На счет place - понял.

[Fixid]

Oleg: индексация обычно делается автоматически сразу после импорта, рекомендую взять nominatim с гита и импортировать базу с помощью него, апач и прочее ставить не обязательно. Там есть правильные бинарники osm2pgsql

Индексация может занимать очень очень много времени (до месяца). Но если это только одна страна и есть 4 нормальных ядра и ssd, то хватит и 3-4 часа

[Oleg]

Fixid: Благодарю за помощь. Буду пробовать nominatim, после отпишу свое решение в SQL для поиска название города.

Как раз есть ssd и 4 ядра :)

[Fixid]

Oleg: там уже искать будет очень просто, placex будет каталогом с полигонами и названиями городов. Если что, то там будут импортированы функции мультиязыкового поиска с приоритетом

[freeExec]

Да у вас координаты кривые, сравните где город, и где вы ищете.
Надо 47.11,51.92, город на месте.

[Oleg]

freeExec: Нет, специально поставил координату туда для примера. Плохой пример оказался :)