Kalkulator rastrów cz.2 – Indeksy wegetacyjne (spektralne)

   W zastosowaniach gisowych często wykorzystywane są efekty teledetekcyjnego przetwarzania zobrazowań. Do jednych z najbardziej popularnych „produktów” należą indeksy wegetacyjne, które są wskaźnikami oceny stanu roślinności, poza tym są także używane do korekty obrazów spektralnych. Obecnie istnieje całkiem spora liczba indeksów spektralnych, z których każdy jest pomocny przy określaniu konkretnych zjawisk. Warto dodać, że wartości indeksów zawierają się w przedziałach liczbowych od – 1,0 do 1,0 lub od 0 do 255.

    W QuantumGIS indeksy możemy obliczyć przy użyciu kalkulatora rastrów. W poniższym przykładzie wykorzystałem zobrazowania pochodzące z LANDSATA 7, które przedstawiają okolice zbiornika Żelazny Most. Wykorzystamy wskaźnik NDVI (Normalized Difference Vegetation Index) obliczany w następujący sposób:

NDVI=(NIR-RED)/(NIR+RED), gdzie

      NIR – kanał podczerwieni

      RED – kanał czerwieni

     Po otworzeniu obrazów dla kanału czerwonego (3) i bliskiej podczerwieni (4) przechodzimy do Raster > Kalkulator rastra...

     W oknie kalkulatora ustawiamy parametry warstwy wyjściowej: nazwa, zasięg, a także format wyjściowy. Jako wyrażenie kalkulatora rastrowego stosujemy powyższy wzór dodając obrazy podwójnym kliknięciem. Po napisaniu formuły zatwierdzamy i otrzymujemy raster wyświetlony w skali szarości, aby lepiej zobaczyć efekt zmieniamy paletę na pseudokolorową.

     Wartości poszczególnych rastrów możemy zobaczyć instalując wtyczkę ValueTools, która je pokazuje w zadokowanym oknie. W przykładzie wskazany piksel ma nastęujące wartości dla rastra NDVI – 0,13; dla RED - 55; dla NIR – 72.

     Oczywiście wykorzystując algorytmy dla innych indeksów spektralnych możemy je w ten sam sposób obliczyć, a tym samym określić gatunkowość drzewostanów (liściaste lub iglaste – MSI) czy lokalizację bagien i mokradeł – VI 3.