Практика использования raster tools
- oldbay
- Участник
- Сообщения: 53
- Зарегистрирован: 14 июл 2015, 11:46
- Репутация: 19
- Откуда: Санкт-Петербург
- Контактная информация:
Практика использования raster tools
Добрый день.
Достаточно давно опубликовал на gis-lab статью Атмосферная коррекция по методу DOS где упомянул свой растровый калькулятор. За прошедшее время калькулятор превратился в "швейцарский нож", о использовании которого хочу рассказать в новой статье:
Статья опубликована: Практика использования raster tools
Исходный код raster_tools: https://github.com/oldbay/raster_tools
Исходный код тестовых примеров: https://gitlab.com/oldbay/raster_tools_examples
Достаточно давно опубликовал на gis-lab статью Атмосферная коррекция по методу DOS где упомянул свой растровый калькулятор. За прошедшее время калькулятор превратился в "швейцарский нож", о использовании которого хочу рассказать в новой статье:
Статья опубликована: Практика использования raster tools
Исходный код raster_tools: https://github.com/oldbay/raster_tools
Исходный код тестовых примеров: https://gitlab.com/oldbay/raster_tools_examples
Последний раз редактировалось oldbay 30 июл 2018, 23:54, всего редактировалось 1 раз.
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Re: Практика использования raster tools
Речь идёт о пакете osgeo? Если да, то в чём отличие-то? Резать растр можно и osgeo.gdal.Warp(), например.В течении продолжительного времени развиваю инструментарий для препарирования георастров. Который является набором рецептов упрощающих взаимодействие с гибкой но сложной структурой python gdal.
Если программа обновляется, то почему версия питона такая древняя?
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
- oldbay
- Участник
- Сообщения: 53
- Зарегистрирован: 14 июл 2015, 11:46
- Репутация: 19
- Откуда: Санкт-Петербург
- Контактная информация:
Re: Практика использования raster tools
Спасибо за отклик!
п.с:
А вот gdal.Warp() точно не использую - не помню сейчас чем он меня не устраивал. Завершающий этап вырезания работает через указание Band.ReadAsArray() индекса угла области и размера самой области.
Вы совершенно правы - всё что делает raster_tools это использует python инструментарий osgeo gdal в удобной мне синтаксической обёртке. Это обычный сахар, надстройка, прослойка Ничего нового не изобретено - просто "заточил" gdal так как нравится и использую, никому ничего не навязывая.
п.с:
А вот gdal.Warp() точно не использую - не помню сейчас чем он меня не устраивал. Завершающий этап вырезания работает через указание Band.ReadAsArray() индекса угла области и размера самой области.
Согласен: python 2.7 потихоньку отмирает(вот и новую джангу выпустили без его поддержки). Со временем портирую на 3-ю версию, там не должно возникнуть особых сложностей, но пока нет мотивации - вся разработка у меня пока идёт на 2.7.
- Эдуард Казаков
- Гуру
- Сообщения: 547
- Зарегистрирован: 23 апр 2014, 17:11
- Репутация: 532
- Откуда: Planet Earth
- Контактная информация:
Re: Практика использования raster tools
Спасибо за статью.
Сейчас многое из этого решается благодаря тому, что в python-gdal >= 2.0 утилиты Warp, Translate, Rasterize и другие доступны как функции с возможностью хранить результат в памяти, но всё равно есть что подсмотреть.
Сейчас многое из этого решается благодаря тому, что в python-gdal >= 2.0 утилиты Warp, Translate, Rasterize и другие доступны как функции с возможностью хранить результат в памяти, но всё равно есть что подсмотреть.
- oldbay
- Участник
- Сообщения: 53
- Зарегистрирован: 14 июл 2015, 11:46
- Репутация: 19
- Откуда: Санкт-Петербург
- Контактная информация:
Re: Практика использования raster tools
Спасибо за отзыв!
raster_tools начал писать ещё под python-gdal 1.х. Сейчас медленно переехал в 2.х . В принципе можно прошерстить код и вывести в 2 ветки - под разные версии gdal, но пока такая задача не стоит.Эдуард Казаков писал(а): ↑01 июн 2018, 17:58Сейчас многое из этого решается благодаря тому, что в python-gdal >= 2.0 утилиты Warp, Translate, Rasterize и другие доступны как функции с возможностью хранить результат в памяти, но всё равно есть что подсмотреть.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 4 гостя