RGISTools (version 1.0.2)

varEVI: Calculate the enhanced vegetation index (EVI)

Description

varEVI computes the enhanced vegetation index (EVI) from the blue, near-infrared (NIR) and red bands.

Usage

varEVI(blue, red, nir, scfun = function(r) {
    r
})

Arguments

blue

a raster with the blue band of the capture.

red

a raster with the red band of the capture.

nir

a raster with the NIR band of the capture.

scfun

a function to re-scale the original pixel values into reflectance (0-1).

Value

An EVI image in raster format.

Details

The enhanced vegetation index (EVI) is a vegetation indicator that improves the sensitivity towards high biomass densities compared to NDVI huete2002overviewRGISTools (See varNDVI). This function is used within ls7FolderToVar, ls8FolderToVar, modFolderToVar and senFolderToVar.

References

huete2002overviewRGISTools

Examples

# NOT RUN {
# path to the cropped and cutted MODIS images for the region of Navarre
wdir <- system.file("ExNavarreVar", package = "RGISTools")
# list all the tif files
files.mod <- list.files(wdir, pattern = "\\.tif$", recursive = TRUE, full.names = TRUE)
# print the MOD09 bands
getRGISToolsOpt("MOD09BANDS")
scale.factor <- 0.0001

# select the red, blue and nir bands
img.mod.red <- raster(files.mod[1]) * scale.factor
img.mod.blue <- raster(files.mod[3]) * scale.factor
img.mod.nir <- raster(files.mod[2]) * scale.factor
# calculate the EVI without scale
img.mod.evi <- varEVI(img.mod.blue, img.mod.red, img.mod.nir)
# calculate the EVI scaling 0-1
img.mod.evi.2 <- varEVI(img.mod.blue, img.mod.red, img.mod.nir,scfun=getRGISToolsOpt("MOD09SCL"))
img.mod.evi.12 <- stack(img.mod.evi,img.mod.evi.2)
# plot the image
spplot(img.mod.evi.12,col.regions=rev(terrain.colors(20)),at = c(seq(0,1,0.05)))
# }