PostGis вопросы по запросам к БД

Suleyma · 2014-03-11T10:55:34.000Z

defond:

При этом, названий аэродромов нет. Т.е. у меня их нет в таблице, или я неправильно запрос сформировал? Или неправильно сформировались таблицы при импорте данных?

Может просто наименований нет? У первых попавшихся трех. Выведите еще колонку osm_id, посмотрю у себя, у меня local развернут.

defond:

А зачем ключик -s?

Чтобы потом обновления накатывать.

Surly · 2014-03-11T11:04:43.000Z

defond:

Пишет, что таблица planet_osm_polygon не существует и пропускается. С остальными так же. Это нормально?

Это нормально, просто сообщение об ошибке сформулировано неудачно. На самом деле таблица, которая не существует – создаётся.
Причём, это поведение описано в документе http://switch2osm.org/serving-tiles/manually-building-a-tile-server-12-04/ , раздел “Importing data into the database”:

The first part of the osm2pgsql output looks scary, but is normal:
…
NOTICE: table “planet_osm_point” does not exist, skipping
…
Don’t be concerned by the NOTICE: entries above. All normal.

defond · 2014-03-11T11:14:30.000Z

Понял, спасибо огромное!

defond · 2014-03-11T11:16:23.000Z

Вопрос - есть вариант sql-запросом рассчитать:

Расстояние между двумя точками
Азимут из точки на какую-то другую точку (извините за каламбур, подвисаю на работе)

Suleyma · 2014-03-11T11:26:35.000Z

defond:

Вопрос - есть вариант sql-запросом рассчитать:

Расстояние между двумя точками

Азимут из точки на какую-то другую точку (извините за каламбур, подвисаю на работе)

Конечно.

ST_Azimuth — Returns the angle in radians from the horizontal of the vector defined by pointA and pointB

ST_Distance — For geometry type Returns the 2-dimensional cartesian minimum distance (based on spatial ref) between two geometries in projected units. For geography type defaults to return spheroidal minimum distance between two geographies in meters.

defond · 2014-03-11T13:08:08.000Z

Спасибо!

defond · 2014-03-11T13:37:18.000Z

При расчете дистанции, нужно указывать (как я понял) район, в котором проводятся расчеты, чтобы учесть искривление поверхности Земли. Т.е. вот это:


units in planar degrees 4326 is WGS 84 long lat unit=degrees

где их можно посмотреть?

defond · 2014-03-11T14:09:02.000Z

Чего то я плыву совсем…

Делаю запрос:


SELECT 'point' AS type,ST_X(way) AS lon,ST_Y(way) AS lat,aeroway,name FROM planet_osm_point 
WHERE aeroway IN ('airport','aerodrome')    
UNION
SELECT 'polygon' AS type,ST_X(ST_Centroid(way)) AS lon,ST_Y(ST_Centroid(way)) AS lat,aeroway,name FROM planet_osm_polygon
WHERE aeroway IN ('airport','aerodrome') LIMIT 50;

Получил я координаты:


"polygon";4359829.4534671;5627019.51677918;"aerodrome";"Международный аэропорт Краснодар"

А что мне показывает этот запрос? О_о

Координаты “Международный аэропорт Краснодар” - 45°02′04″ с. ш. 39°09′00″ в. д.

Вот страничка Вики - Аэропорт Краснодар

Мягко говоря, я не понимаю…

defond · 2014-03-11T14:50:29.000Z

Добрый человек подсказал:

Это EPSG:3857 тебе нужна EPSG:4326.

Таким образом нужно получить данные в EPSG:4326. Ищу пример запроса.

Alexandr_Zeinalov · 2014-03-11T16:51:47.000Z

Собственно вместо way надо использовать ST_Transform(way,4326) и получим геометрию в 4326.

defond · 2014-03-12T06:48:57.000Z

Все работает! Спасибо огромное!

Вот такой запрос нужен для получения координат в “правильном” виде. (вдруг кому пригодится)


SELECT 'point' AS type,ST_X(ST_Transform(way,4326)) AS lon,ST_Y(ST_Transform(way,4326)) AS lat,aeroway,name FROM planet_osm_point 
WHERE aeroway IN ('airport','aerodrome')    
UNION
SELECT 'polygon' AS type,ST_X(ST_Centroid(ST_Transform(way,4326))) AS lon,ST_Y(ST_Centroid(ST_Transform(way,4326))) AS lat,aeroway,name FROM planet_osm_polygon
WHERE aeroway IN ('airport','aerodrome') LIMIT 3;

Suleyma · 2014-03-12T07:50:24.000Z

Мне кажется, более естественно сразу загружать в 4326 (--proj 4326). Оно и запросы проще будут.
select name from planet_osm_polygon where st_dwithin(way, st_transform(setsrid(makepoint(39.0900,45.0204),4326),900913)::geometry,1);
или
select name from planet_osm_polygon where st_dwithin(way, setsrid(makepoint(39.09,45.02),4326)::geometry,1);

Suleyma · 2014-03-12T08:02:03.000Z

defond:

При расчете дистанции, нужно указывать (как я понял) район, в котором проводятся расчеты, чтобы учесть искривление поверхности Земли. Т.е. вот это:
units in planar degrees 4326 is WGS 84 long lat unit=degrees
где их можно посмотреть?

Вообще удобнее работать с типом geography, для него всякие там ST_Distance, ST_DWithin, пользуют метры.

ST_Distance(ST_Transform(way,4326)::geography, SetSRID(MakePoint(lon,lat),4326)::geography)

defond · 2014-03-12T08:25:48.000Z

Suleyma:

Мне кажется, более естественно сразу загружать в 4326 (--proj 4326). Оно и запросы проще будут.
select name from planet_osm_polygon where st_dwithin(way, st_transform(setsrid(makepoint(39.0900,45.0204),4326),900913)::geometry,1);
или
select name from planet_osm_polygon where st_dwithin(way, setsrid(makepoint(39.09,45.02),4326)::geometry,1);

Спасибо. Попробую.

Только если можно - объясняйте чуть проще, я не специалист вообще, тяжеловато разобраться. Через пару недель будет проще, а сейчас тяжело еще.

akks · 2014-03-17T14:35:50.000Z

Я попробовал скомпилировать свежий osm2pgsql на Visual Studio.

Пришлось поправлять исходники, pbf поддержать не получилось (кто-нибудь видел живой protobuf-c для Windows?).

Кому интересно, можете потестировать: https://dl.dropboxusercontent.com/u/63393258/osm2pgsql_win_unstable.zip
Исходный код (переписал под cmake): https://github.com/alex85k/osm2pgsql. (нужен VS2013, он более совместим c C99).
(может затребовать runtime: http://www.microsoft.com/ru-ru/download/details.aspx?id=40784)

Suleyma · 2014-03-20T08:26:19.000Z

Может подскажет кто. Почему osm2pgsql может не загрузить точку из данных? Смотрю загрузку osmosis-ом точка есть, osm2pgsql - нет. Соответственно и дорога поличилась обрезана. например node - 2246232221, way - 215216040

akks · 2014-03-20T09:05:10.000Z

Suleyma:

Может подскажет кто. Почему osm2pgsql может не загрузить точку из данных? Смотрю загрузку osmosis-ом точка есть, osm2pgsql - нет. Соответственно и дорога получилась обрезана. например node - 2246232221, way - 215216040

64-битные ID в некоторых сборках по умолчанию выключены, особенно если под Windows

Кто-нибудь знает рабочую сборку, кроме моей экспериментальной из предыдущего поста?

olamedia · 2014-03-26T15:40:28.000Z

Suleyma:

Может подскажет кто. Почему osm2pgsql может не загрузить точку из данных? Смотрю загрузку osmosis-ом точка есть, osm2pgsql - нет. Соответственно и дорога поличилась обрезана. например node - 2246232221, way - 215216040

А там случаем не составная дорога? один и тот же osm_id могут иметь несколько рядов в таблицах

Suleyma · 2014-03-28T04:43:39.000Z

olamedia:

А там случаем не составная дорога? один и тот же osm_id могут иметь несколько рядов в таблицах

Да нет, обычная дорога. В JOSM немного её обновил, теперь нормально. а было так:

Suleyma · 2014-03-28T04:53:05.000Z

поторопился с нормально

vs