Расчет NDVI в GRASS 6.4.3
Аннотация:в статье представлена инструкция по созданию проекта базы геоданных, области и набора данных в GRASS и расчету вегетационного индекса NDVI с помощью программы GRASS 6.4.3
GRASS (англ. Geographic Resources Analysis Support System — система для обработки географической информации. Аббревиатура складывается в английское слово grass — трава.) Геоинформационная система с открытым исходным кодом. Поддерживает большое количество форматов. Выпущены версии, работающие под многими операционными системами — Microsoft Windows, Apple Mac OS X, POSIX-совместимыми вообще и Linux в частности http://ru.wikipedia.org/wiki/GRASS_(%D0%93%D0%98%D0%A1).
Для того, чтоб начать любую работу в GRASS необходимо создать базу геоданных (DATABASE)
Запускаем Grass в его графической оболочке, через терминал Tc|Tk GUI,
получаем окно
Выбираем «Определить новую область» - «Привязанный файл»
нажимаем «Войти в GRASS», видим на экране три окна
Для дальнейшей работы необходимо импортировать растровые слои. Для расчёта NDVI нам потребуются файлы с 3 и 4 каналом.
Задаём имя выходной растровой карты
Нажимаем «Запустить». В окне «Вывод» видим сообщение
Повторяем операцию для создания слоя L8_4. Расчёт NDVI производится по формуле (http://gis-lab.info/qa/ndvi.html)
где NIR - отражение в ближней инфракрасной области спектра RED - отражение в красной области спектра. В связи с тем, что при регистрации данных Landsat физические величины яркости каналов сжимаются в диапазон [1,255] (преобразование L→DN),необходимо провести радиометрическую коррекцию/ Исходные максимальные и минимальные значения яркостей приводятся в файле метаданных (MTL) к снимку: LMAX_BAND3 (RADIANCE_MAXIMUM_BAND_3) = 705.61206 LMIN_BAND3(RADIANCE_MINIMUM_BAND_3) = -58.26967 LMAX_BAND4(RADIANCE_MAXIMUM_BAND_4) = 597.61023 LMIN_BAND3(RADIANCE_MINIMUM_BAND_4) = -49.35085 Обратное преобразование (DN→L) производится по формуле:
При написании статьи были использованы следующие материалы: 1. Знакомство с ГИС GRASS Колесов Д. А., Мурый А. А.http://gisconf.ru/talks/lect17-kolesov-grass.pdf 2. NDVI - [теория] и практика http://gis-lab.info/qa/ndvi.html