Работа с геоданными?

Question

[DDwrt100]

Добрый день. Ответьте пожалуйста на вопрос по работе геолокацией. Я никогда не сталкивался с подобными задачами, и некоторые вещи вызывают затруднение.
Ситуация:
У меня есть сырые данные в базе данных, в этой базе есть два столбца(которые называются широта и долгота).
В виде запроса приходит маршрут в виде массива (широта долгота). Мне надо отобрать записи , которые относятся к этому маршруту.
Как мне правильно построить запрос ?
Так же у меня вызывает вопросы техника расчета принадлежности одной точки к другой. У кого есть ссылки на какие нибудь толковые статьи где разбирается с нуля все эта тематика ?

UPD: ok. у меня сейчас там clickhouse от яндекса, но и postgress тоже есть.
Собственно до чего я сейчас додумался:
Я вижу в документации базы данных, создание гео индекса
SELECT geoToH3(37.79506683, 55.71290588, 15) as h3Index
в текущих координатах это должно получиться :
644325524701193974 что то подобное.
далее я беру координаты из запроса и тоже конвертируя их в индекс могу сравнить принадлежность к чему ? Как сравнить единичный индекс с маршрутом ?
Как это делается в реальных проектах? Кто может рассказать?

Comments

[Владислав Лысков]

используй GeoJSON(TOPOJSON)

Answer 1

[Александр]

Лучше всего когда точка сохранена как объект геометрии - в postgresql есть возможность определения пересечения объектов геометрии.

А если два отдельных поля - то надо выполнять расчет расстояния и если оно меньше N метров то точка нужная.
для коротких дистанций подходит расчет прямоугольного треугольника, для больших - формула Хаверсина https://habr.com/ru/post/179157/

Понятно что расстояние небольшое и можно считать по простому - через r^2 = (x0-x)^2+(y0-y)^2

если есть геометрические функции - то точка - это окружность с радиусом в N метров которая пересекается с ломанной линией пути.

Если их нет, то придется весь расчет включить в запрос и он получится неприлично длинный. Приблизительно - можно найти сравнив все точки пути с точками записей найдя минимальное расстояние.

Answer 2

[freeExec]

Избавиться от колонок широта и долгота и воспользоваться пространственным расширением для вашей БД. А дальше строите буферы вокруг маршрута и получаете те точки которые в него попали.
На Гислабе полно статей по PostGIS, тема не простая, с наскока можно не осилить
gis-lab.info/docs/postgis/manual/ch06.html

Индекс строят не H3 - это больше для точечных объектов, а GIST
А запросы получаются примерно таких планов, достать дома вдоль дороги не далее 50 "псевдометров":

SELECT hw.highway, hw.name, bld."addr:housenumber", bld.building, hw.way, bld.way FROM planet_osm_line AS hw 
INNER JOIN planet_osm_polygon AS bld ON ST_Intersects(ST_Buffer(hw.way, 50), bld.way)
WHERE hw.highway IS NOT NULL AND bld."addr:housenumber" IS NOT NULL

LIMIT 30;

Comments

[DDwrt100]

Правильно ли я понял?
1) для работы с геометрией нам необходимо специальное расширение базы данных данных.
Это расширение содержит в себе методы для работы с геометрией?
2) К базе данных нам необходимо добавить еще одну базу данных, которая содержит координаты области(пусть будет планеты)
3) В случае если у меня просто числовые колонки широты и долготы , то мне необходимо их конвертировать в специальную яйчейку, которую назовем гео. А если у меня нет специального поля гео в базе данных, то мне нужно как то посчитать в нужный формат эти два столбца ?

[freeExec]

DDwrt100, 1) Нужно
2) Тут я не совсем понял что она содержит
3) Можно генерить на лету, но смысла в этом нет, т.к. данные не меняются и логично преобразовать один раз всё.

[DDwrt100]

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

[freeExec]

DDwrt100, Всё равно я вас не понимаю о каких географических областях вы говорите.
Вот вы можете сказать точка принадлежит фигуре или нет? Какие географические области вы использовали для ответа?

[DDwrt100]

freeExec, Я понял что Вы имеете ввиду.

[DDwrt100]

freeExec, Можно еще один вопрос?
с Postgress все понятно, расширяем плагином и радуемся жизни.
Проблема с Clickhouse.
Как я понял ClickHouse не имеет такого плагина, однако имеет несколько функций для обработки геоданных.
https://clickhouse.yandex/docs/en/query_language/f...
вот описанные методы.
За что цепляется глаз это перевод точек в различные геохеши
geohashEncode
geoToH3
Возможно ли как то сравнивать Hash точки с hash линии.
Что то типа запроса , возми хэш записи(геокординат), и сравни с хэшем моего маршрута, и скажи мне относится эта запись к этому маршруту ?

[freeExec]

DDwrt100, С геохешем тоже возможно получится, т.к. он позволяет сказать близко ли расположены объекты. Но для линий придётся генерить кучу точек вдоль неё.