Расчет 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)

$NDVI={\frac {NIR-RED}{NIR+RED}}$

где 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) производится по формуле: