Пример использования утилит GDAL для вычисления яркостных характеристик снимков (на примере выявления рубок по зимним данным Landsat): различия между версиями
Gornak (обсуждение | вклад) Нет описания правки |
Gornak (обсуждение | вклад) Нет описания правки |
||
| Строка 1: | Строка 1: | ||
{{Статья|Черновик}} | {{Статья|Черновик}} | ||
{{Аннотация|Пример решения технических задач, возникающих при выявлении изменений отражательных свойств поверхности по разновременным материалам ДЗЗ}} | {{Аннотация|Пример решения технических задач, возникающих при выявлении изменений отражательных свойств поверхности по разновременным материалам ДЗЗ}} | ||
Одним из способов выявления изменений между сериями разновременных снимков является проведение различных математических операций между значениями яркости пикселей. Настоящая статья создана с целью помочь начинающим пользователям освоить основы техники данного приема с использованием свободного ПО и является опытом её автора. Для примера выбрана достаточно простая практическая задача и соответствующие ей исходные данные. | Одним из способов выявления изменений между сериями разновременных снимков является проведение различных математических операций между значениями яркости пикселей. Настоящая статья создана с целью помочь начинающим пользователям освоить основы техники данного приема с использованием свободного ПО и является опытом её автора. Для примера выбрана достаточно простая практическая задача и соответствующие ей исходные данные. | ||
== | ==Используемое ПО== | ||
Командная строка Windows 7, Python27, GDAL 1.11.1 | Командная строка Windows 7, Python27, GDAL 1.11.1 | ||
==Исходные данные== | ==Исходные данные== | ||
| Строка 10: | Строка 11: | ||
<pre>gdalinfo rastr</pre> | <pre>gdalinfo rastr</pre> | ||
*здесь и далее приведен пример синтаксиса командной строки | *здесь и далее приведен пример синтаксиса командной строки | ||
rastr - полный путь к набору растровых данных | *rastr - полный путь к набору растровых данных | ||
Результат работы утилиты, будет представлен в виде текста в командной строке: | Результат работы утилиты, будет представлен в виде текста в командной строке: | ||
[[Файл:Ginfextent.jpg| | [[Файл:Ginfextent.jpg|лево|300|Вывод значений координат углов изображения(экстент) в командной строке]] | ||
Версия от 09:25, 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 - полный путь к набору растровых данных
Результат работы утилиты, будет представлен в виде текста в командной строке:
