Форумы GIS-Lab.info

Добрый день.
Достаточно давно опубликовал на gis-lab статью Атмосферная коррекция по методу DOS где упомянул свой растровый калькулятор. За прошедшее время калькулятор превратился в "швейцарский нож", о использовании которого хочу рассказать в новой статье:
Статья опубликована: Практика использования raster tools
Исходный код raster_tools: https://github.com/oldbay/raster_tools
Исходный код тестовых примеров: https://gitlab.com/oldbay/raster_tools_examples

В течении продолжительного времени развиваю инструментарий для препарирования георастров. Который является набором рецептов упрощающих взаимодействие с гибкой но сложной структурой python gdal.

Речь идёт о пакете osgeo? Если да, то в чём отличие-то? Резать растр можно и osgeo.gdal.Warp(), например.

Если программа обновляется, то почему версия питона такая древняя?

Спасибо за отклик!

rhot писал(а): ↑

01 июн 2018, 16:23

Речь идёт о пакете osgeo? Если да, то в чём отличие-то? Резать растр можно и osgeo.gdal.Warp(), например.

Вы совершенно правы - всё что делает raster_tools это использует python инструментарий osgeo gdal в удобной мне синтаксической обёртке. Это обычный сахар, надстройка, прослойка Ничего нового не изобретено - просто "заточил" gdal так как нравится и использую, никому ничего не навязывая.
п.с:
А вот gdal.Warp() точно не использую - не помню сейчас чем он меня не устраивал. Завершающий этап вырезания работает через указание Band.ReadAsArray() индекса угла области и размера самой области.

rhot писал(а): ↑

01 июн 2018, 16:23

Если программа обновляется, то почему версия питона такая древняя?

Согласен: python 2.7 потихоньку отмирает(вот и новую джангу выпустили без его поддержки). Со временем портирую на 3-ю версию, там не должно возникнуть особых сложностей, но пока нет мотивации - вся разработка у меня пока идёт на 2.7.

Спасибо за статью.
Сейчас многое из этого решается благодаря тому, что в python-gdal >= 2.0 утилиты Warp, Translate, Rasterize и другие доступны как функции с возможностью хранить результат в памяти, но всё равно есть что подсмотреть.

Эдуард Казаков писал(а): ↑

01 июн 2018, 17:58

Спасибо за статью.

Спасибо за отзыв!

Эдуард Казаков писал(а): ↑

01 июн 2018, 17:58

Сейчас многое из этого решается благодаря тому, что в python-gdal >= 2.0 утилиты Warp, Translate, Rasterize и другие доступны как функции с возможностью хранить результат в памяти, но всё равно есть что подсмотреть.

raster_tools начал писать ещё под python-gdal 1.х. Сейчас медленно переехал в 2.х . В принципе можно прошерстить код и вывести в 2 ветки - под разные версии gdal, но пока такая задача не стоит.