Greenness Vegetation Index is obtained from the Tasseled Cap Transformation.
gevi(img, sat = "Landsat8OLI")
| img | It could be RasterStack or RasterBrick. |
|---|---|
| sat | Specify satellite and sensor type (Landsat5TM, Landsat7ETM or Landsat8OLI). |
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 Sensor | Name of bands |
| Landsat4TM | blue,green,red,nir,swir1,swir2 |
| Landsat5TM | blue,green,red,nir,swir1,swir2 |
| Landsat7ETM+ | blue,green,red,nir,swir1,swir2 |
| Landsat8OLI | blue,green,red,nir,swir1,swir2 |
| Landsat8OLI-Li2016 | coastal,blue,green,red,nir,swir1,swir2 |
| Sentinel2MSI | coastal,blue,green,red,nir-1,mir-1,mir-2 |
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)