Расчет NDVI в GRASS 6.4.3
Аннотация:в статье представлена инструкция по расчету нормализованного относительного индекса растительности NDVI с помощью программы GRASS 6.4.3
Общие сведения
NDVI (Normalized Difference Vegetation Index) - нормализованный относительный индекс растительности - простой количественный показатель количества фотосинтетически активной биомассы (обычно называемый вегетационным индексом). Один из самых распространенных и используемых индексов для решения задач, использующих количественные оценки растительного покрова. Подробнее о самом индексе в статье "NDVI - [теория] и практика"http://gis-lab.info/qa/ndvi.html
Импорт данных в проект
Для дальнейшей работы необходимо импортировать растровые слои. Для расчёта NDVI нам потребуются файлы с 3 и 4 каналом.
Выбираем «Определить новую область» - «Привязанный файл»
нажимаем «Войти в GRASS», видим на экране три окна
Задаём имя выходной растровой карты
Нажимаем «Запустить». В окне «Вывод» видим сообщение
Повторяем операцию для создания слоя L8_4.
После импорта растров, но перед расчетом NDVI Точная подгонка региона под охват и разрешение выбранного растра:
g.region rast=растр -a
Флаг -a позволяет подогнать с максимальной точностью, без этого флага, подгонка разрешения выполняется с округлением. Этот подход предпочтительней.
Другой вар.: подогнать охват под выбранный растр и задать разрешение равное ровно 30 м. (в исходных растрах Landsat пространственное разрешение м.б. не 30 м, а, например, 29.85 КОД: ВЫДЕЛИТЬ ВСЁ g.region rast=растр res=30.0 -a
Посмотреть параметры текущего региона:
КОД: ВЫДЕЛИТЬ ВСЁ
g.region -p
Расчет NDVI в GRASS
Расчёт 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) производится по формуле:
Открываем калькулятор растров, в строки А и B добавляем слои L8_3 и L8_4. В формулу вставляем значения из метафайла, обозначаем имя выходного слоя, нажимаем Запустить.Если всё корректно выполнено,то в окне Вывод вы получите сообщение о том, что операция выполнена и создан новый слой.
Повторяем операцию для второго слоя.
Теперь раcсчитываем NDVI. Для этого в калькулятор добавляем слои, созданные по результатам радиометрической коррекции, записываем формулу и новый слой L8_ndvi
На этом расчет вегетационного индекса NDVI в программе GRASS окончен.
Ссылки по теме: 1. Знакомство с ГИС GRASS Колесов Д. А., Мурый А. А.http://gisconf.ru/talks/lect17-kolesov-grass.pdf 2. NDVI - [теория] и практика http://gis-lab.info/qa/ndvi.html 4. Структура базы геоданных ГИС GRASS и способы ее создания http://wiki.gis-lab.info/w/%D0%A1%D1%82%D1%80%D1%83%D0%BA%D1%82%D1%83%D1%80%D0%B0_%D0%B1%D0%B0%D0%B7%D1%8B_%D0%B3%D0%B5%D0%BE%D0%B4%D0%B0%D0%BD%D0%BD%D1%8B%D1%85_%D0%93%D0%98%D0%A1_GRASS_%D0%B8_%D1%81%D0%BF%D0%BE%D1%81%D0%BE%D0%B1%D1%8B_%D0%B5%D0%B5_%D1%81%D0%BE%D0%B7%D0%B4%D0%B0%D0%BD%D0%B8%D1%8F 5.http://ru.wikipedia.org/wiki/GRASS_(%D0%93%D0%98%D0%A1)