Расчет NDVI в GRASS 6.4.3: различия между версиями
мНет описания правки |
Нет описания правки |
||
| Строка 7: | Строка 7: | ||
получаем окно | получаем окно | ||
[[Файл:Grass_ndvi2.jpg]] | [[Файл:Grass_ndvi2.jpg]] | ||
Выбираем «Определить новую область» - «Привязанный файл» | Выбираем «Определить новую область» - «Привязанный файл» | ||
[[Файл:NDV_oblast.jpg]] | [[Файл:NDV_oblast.jpg]] | ||
нажимаем «Войти в GRASS», видим на экране три окна | нажимаем «Войти в GRASS», видим на экране три окна | ||
[[Файл:Ndvi_grass_3.jpg]] | [[Файл:Ndvi_grass_3.jpg]] | ||
Для дальнейшей работы необходимо импортировать растровые слои. Для расчёта NDVI нам потребуются файлы с 3 и 4 каналом. | Для дальнейшей работы необходимо импортировать растровые слои. Для расчёта NDVI нам потребуются файлы с 3 и 4 каналом. | ||
[[Файл:Import sloya.jpg]] | [[Файл:Import sloya.jpg]] | ||
Задаём имя выходной растровой карты | Задаём имя выходной растровой карты | ||
[[Файл:Import rastra.jpg]] | [[Файл:Import rastra.jpg]] | ||
Нажимаем «Запустить». В окне «Вывод» видим сообщение | Нажимаем «Запустить». В окне «Вывод» видим сообщение | ||
[[Файл:Vwvod1.jpg]] | [[Файл:Vwvod1.jpg]] | ||
Версия от 09:08, 1 апреля 2014
Аннотация:в статье представлена инструкция по расчету вегетационного индекса 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) производится по формуле: