Пример использования утилит GDAL для вычисления яркостных характеристик снимков (на примере выявления рубок по зимним данным Landsat): различия между версиями

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

Версия от 08:56, 31 мая 2016

Эта страница является черновиком статьи.


Пример решения технических задач, возникающих при выявлении изменений отражательных свойств поверхности по разновременным материалам ДЗЗ

Одним из способов выявления изменений между сериями разновременных снимков является проведение различных математических операций между значениями яркости пикселей. Настоящая статья создана с целью помочь начинающим пользователям освоить основы техники данного приема с использованием свободного ПО и является опытом её автора. Для примера выбрана достаточно простая практическая задача и соответствующие ей исходные данные.

Использованное ПО

Командная строка Windows 7, Python27, GDAL 1.11.1

Исходные данные

Зимние безоблачные снимки Landsat 8 OLI панхроматический(_B8) канал одной сцены (отмечу что для моего региона в средней полосе Европейской части, удалось найти порядка по 1 такому снимку за весь период работы сенсора(с апреля 2013 г.) в периоде февраль-март за 2015 и 2016 г.)

Создание общей области покрытия

Для корректной работы утилиты gdal_calc.py необходимо, чтобы входящие растры имели одинаковый экстент. Получить сведения об экстенте можно используя утилиту gdalinfo.exe:

gdalinfo rastr
*здесь и далее приведен пример синтаксиса командной строки
rastr - полный путь к набору растровых данных