Расчет NDVI в GRASS 6.4.3
Аннотация:в статье представлена инструкция по расчету вегетационного индекса NDVI с помощью программы GRASS 6.4.3
Запускаем 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) производится по формуле: