Greenness Vegetation Index (GEVI)

Source: R/gevi.R

Greenness Vegetation Index is obtained from the Tasseled Cap Transformation.

gevi(img, sat = "Landsat8OLI")

Arguments

img

It could be RasterStack or RasterBrick.

sat

Specify satellite and sensor type (Landsat5TM, Landsat7ETM or Landsat8OLI).

Details

Currently implemented for satellites such as Landsat-4 TM, Landsat-5 TM, Landsat-7 ETM+, Landsat-8 OLI and Sentinel2. The input data must be in top of atmosphere reflectance (toa). Bands required as input must be ordered as:

Type of SensorName of bands
Landsat4TMblue,green,red,nir,swir1,swir2
Landsat5TMblue,green,red,nir,swir1,swir2
Landsat7ETM+blue,green,red,nir,swir1,swir2
Landsat8OLIblue,green,red,nir,swir1,swir2
Landsat8OLI-Li2016coastal,blue,green,red,nir,swir1,swir2
Sentinel2MSIcoastal,blue,green,red,nir-1,mir-1,mir-2

Examples

library(ForesToolboxRS)
library(raster)

# Load an example dataset
data(img_l8)

# Tasseled-cap using Landsat8OLI
gevi_index <- gevi(img = img_l8 / 10000, sat = "Landsat8OLI")

# Improving the plot histogram
gevi_index[gevi_index > 1] <- NA
gevi_index[gevi_index < -1] <- NA
plot(gevi_index)