GIS-LAB

Географические информационные системы и дистанционное зондирование


Форумы GIS-Lab.info

Геоинформационные системы (ГИС) и Дистанционное зондирование Земли


Номера тайлов по полигону

Вопросы общего характера по ГИС и дистанционному зондированию, не связанные с конкретным ПО.

Номера тайлов по полигону

Сообщение ANAT01 » 16 июл 2017, 17:59

Как получить список номеров тайлов на которые попадает полигон (не bbox, а по границе полигона)?

Тут в принципе написано про пересчет lat,long но как пройтись по границе?!
  • 0

Аватара пользователя
ANAT01
Активный участник
 
Зарегистрирован: 18 апр 2012
 
Сообщения: 153
Репутация: 39

Re: Номера тайлов по полигону

Сообщение Александр Мурый » 16 июл 2017, 20:40

Может быть, код из плагина Tile Cover сможет помочь?
  • 0

Редактор материалов, модератор форума
Александр Мурый
званий не имею
 
Зарегистрирован: 26 сен 2009
Откуда: Москва (на карте)
 
Сообщения: 4912
Репутация: 676
Статьи: 3
Проекты: 5/1

Re: Номера тайлов по полигону

Сообщение Zverik » 17 июл 2017, 11:08

pip install polytiles
polytiles -p polygon.poly -x tiles.lst -z minzoom maxzoom
  • 3

Zverik
Активный участник
 
Зарегистрирован: 11 сен 2011
Откуда: Москва (на карте)
 
Сообщения: 139
Репутация: 86
Статьи: 2

Re: Номера тайлов по полигону

Сообщение ANAT01 » 17 июл 2017, 18:36

Спасибо!!! А то я уж думал штудировать shapely/fiona и алгоритмы, снились кошмары и пробивал пот. :)

Для полноты картины:

- Для работы polytiles нужен файл в формате Osmosis polygon filter file format.
- poly умеют готовить несколько скриптов, я выбрал ogr2poly.py

Код: Выделить всё
# готовим окружение
virtualenv venv
. ./venv/bin/activate
pip install mapnik
pip install polytiles

# готовим poly
python ogr2poly.py mypolygon.geojson

# готовим список тайлов
polytiles -p mypolygon_0.poly -x tiles.lst -z 19 20


Результат еще не успел проверить...

PS:
Александр Мурый писал(а):Может быть, код из плагина Tile Cover сможет помочь?

Он делает по bbox. В принципе так можно делать, а потом находить пересечение с заданной фигурой, но это не совсем то что нужно.

UPD:

Проверил визуально

Здесь построена сетка тайлов с помощью Tile Cover и слинкована с таблицей polytiles. "Пересечение" данных выделено желтым.
Изображение
  • 1

Аватара пользователя
ANAT01
Активный участник
 
Зарегистрирован: 18 апр 2012
 
Сообщения: 153
Репутация: 39


Вернуться в Общие вопросы

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 4


(Геокруг)

© GIS-Lab и авторы, 2002-2013. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).