GIS i okolice: QGIS, OrfeoToolbox i indeksy spektralne

27 wrz 2013

QGIS, OrfeoToolbox i indeksy spektralne

Niedawno pisałem o możliwości obliczenia indeksów wegetacyjnych (spektralnych) w Quantum GIS przy użyciu Kalkulatora Rastrów. W programie istnieje również inny sposób określenia indeksów. W aplikacji, mamy dostęp do potężnego narzędzia – Sextante, które jest zintegrowane z QGIS-em pod nazwą Geoprocesing (dawniej była to wtyczka). Wśród jego modułów znajduje się Orfeo Toolbox, dzięki któremu możemy korzystać z dużej liczby algorytmów w celu przetwarzania obrazów teledetekcyjnych (opis jak skonfigurować moduł Orfeo Toolbox znajduje się w poście Teledetekcja w QGIS). W OT obliczanie indeksów wegetacyjnych dla wczytanego zobrazowania zaimplementowano w Feature Extraction > Radiometric Indices.

Widoczne indeksy zostały podzielone na 3 grupy:

Wegetacja:

NDVI - Normalized Difference Vegetation Index
TNDVI – Transformed Normalized Difference Vegetation Index
RVI - Ratio Vegetation Index
SAVI - Soil Adjusted Vegetation Index
TSAVI - Transformed Soil Adjusted Vegetation Index
MSAVI - Modified Soil Adjusted Vegetation Index
MSAVI2 - Modified Soil Adjusted Vegetation Index 2
GEMI - Global Environment Monitoring Index
IPVI - Infrared Percentage Vegetation Index
LAIFromNDVILog - Leaf Area Index from NDVI
LAIFromRefLinear - Leaf Area Index from reﬂectance
LAIFromNDVIFormo - Leaf Area Index from reﬂectance. Dotyczy tylko Formosata 2.

Woda:

NDWI - Normalized Difference Water Index
NDWI2 - Normalized Difference Water Index
MNDWI - Modified Normalized Difference Water Index
NDPI - Normalized Difference Pond Index
NDTI - Normalized Difference Turbidity Index

Gleba:

RI - Redness Index
CI - Color Index
BI - Brilliance Index
BI2 - Brilliance Index

W celu obliczenia interesującego nas indeksu w oknie wczytujemy obraz wielospektralny (jak wykonać taki obraz w QuantumGIS z poszczególnych kanałów spektralnych opisałem w poście Kompozycje barwne w QGIS) i zaznaczamy odpowiednią kolejność zakresów naszego zobrazowania np. dla Landsata 7 kanałem czerwonym będzie – 3, a kanałem podczerwonym – 4. Następnie z rozwijalnej listy wybieramy indeks wegetacyjny. Po określeniu miejsca zapisu pliku wynikowego i wygenerowaniu finalnego obrazu mamy do dyspozycji scenę z obliczonym indeksem.