Как в PostGIS получить точки внутри заданной площади?

Question

[Павел]

У меня есть таблица с местами. В таблице есть колонка coordinates типа geography(POINT), координаты с гугл карт.
Задача выбрать места, которые входят в определенную область на карте (видимый на фронте участок гугл карт).

С фронта мне присылают координаты нижнего левого и верхнего правого угла, теоретически можно попробовать получить все четыре.

Не могу разобраться в доках к постгису.. смотрел на функции ST_Contains и ST_CoveredBy, но как-то не очень понял.. в доках примеры не очень ясные. А может не то смотрел.

Answer 1

[freeExec]

Из этих координат углов строите линию. Затем ищете пересечение охватывающих прямоугольников для данной линии и ваших точек.

SELECT * FROM table WHERE geom && 'LINESTRING(0 0, 3 3)'::geometry

П.С. не забудь правильно указать систему координат.

Comments

[Павел]

Супер. Выглядит как что-то, что работает) А где подробнее можно почитать, как ты получил именно такой запрос. И какая здесь будет разница в указании ::geometry или ::geography. Точки все-таки хранятся как география, а не геометрия.

П.С. не забудь правильно указать систему координат.

это ты про srid? вроде бы же по умолчанию используется 4326?

[freeExec]

Павел, Известно где - в мануале postgis.net/docs/manual-2.0/reference.html
Лучше явно указывать SRID, во избежания недопонимания. Особенно, что касается расстояний, а то указывают числа в метрах, когда система от них ждёт в градусах. А потом удивляются, где это на глобусе 500 градусов.

[freeExec]

Павел, Странно, ну обверните линию ещё в Box2D('LINESTRING(0 0, 3 3)'::geometry)

[Павел]

freeExec, Ну, я в итоге, из двух координат получил другие 2 вершины, и в LINESTRING просо вставил 4 пары координат. Теперь все как надо.

[Павел]

freeExec, а Box2D, видимо не работает с колонкой типа geography. Postgres выбрасывает ошибку

operator does not exist: geography && box2d

[freeExec]

Павел, Оператор && оперирует 2D-box'сами, поэтому нет нужны их создавать намеренно. Я сходил проверил, всё работает с одной линей. Единственное потребовалось, чтобы обе геометрии были в одной системе координат.

SELECT ST_AsText(ST_Transform(way, 4326)) AS wkt FROM public.planet_osm_point WHERE way && ST_Transform(ST_GeomFromText('LINESTRING(50 53, 51 54)', 4326), 3857)