Пример использования утилит GDAL для вычисления яркостных характеристик снимков (на примере выявления рубок по зимним данным Landsat)
Пример решения технических задач, возникающих при выявлении изменений отражательных свойств поверхности по разновременным материалам ДЗЗ
Одним из способов выявления изменений между сериями разновременных снимков является проведение различных математических операций между значениями яркости пикселей. Настоящая статья создана с целью помочь начинающим пользователям освоить основы техники данного приема с использованием свободного ПО и является опытом её автора. Для примера выбрана достаточно простая практическая задача и соответствующие ей исходные данные.
Используемое ПО
Командная строка Windows 7, Python27, GDAL 1.11.1
Исходные данные
Зимние безоблачные снимки Landsat 8 OLI панхроматический(_B8) канал одной сцены (отмечу что для моего региона в средней полосе Европейской части, удалось найти порядка по 1 такому снимку за весь период работы сенсора(с апреля 2013 г.) в периоде февраль-март за 2015 и 2016 г.)
Создание общей области покрытия
Для корректной работы утилиты gdal_calc.py необходимо, чтобы входящие растры имели одинаковый экстент. Получить сведения об экстенте можно используя утилиту gdalinfo.exe:
gdalinfo rastr
- здесь и далее приведен пример синтаксиса командной строки
- rastr - полный путь к набору растровых данных
Результат работы утилиты, будет представлен в виде текста в командной строке:
