Расчет NDVI в GRASS 6.4.3

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
Эта страница является черновиком статьи.


Аннотация:в статье представлена инструкция по расчету вегетационного индекса NDVI с помощью программы GRASS 6.4.3

Запускаем Grass в его графической оболочке, через терминал Tc|Tk GUI, Grass ndvi.jpg получаем окно Grass ndvi2.jpg

Выбираем «Определить новую область» - «Привязанный файл» 

NDV oblast.jpg нажимаем «Войти в GRASS», видим на экране три окна Ndvi grass 3.jpg Для дальнейшей работы необходимо импортировать растровые слои. Для расчёта NDVI нам потребуются файлы с 3 и 4 каналом. Import sloya.jpg Задаём имя выходной растровой карты Import rastra.jpg Нажимаем «Запустить». В окне «Вывод» видим сообщение Vwvod1.jpg Повторяем операцию для создания слоя L8_4. Расчёт NDVI производится по формуле (http://gis-lab.info/qa/ndvi.html) Ndvi формула.gif где, 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) производится по формуле: