rLandsat8 - an R interface to Landsat8

Normalized Difference Vegetation Index (NDVI)

The NDVI formula is:

$NDVI\quad =\quad \frac { \rho NIR\quad -\quad \rho Red }{ \rho NIR\quad +\quad \rho Red }$

Example

library(rLandsat8)

setwd("~/Downloads")
product  <- "LC82040322013219LGN00"
l <- ReadLandsat8(product)

ndvi <- ToNDVI(l)

Modified Normalized Difference Water Index (MNDWI)

The MNDWI formula is:

$MNDWI\quad =\quad \frac { \rho Green\quad -\quad \rho SWIR1 }{ \rho Green\quad +\quad \rho SWIR1 }$

Example

library(rLandsat8)

setwd("~/Downloads")
product  <- "LC82040322013219LGN00"
l <- ReadLandsat8(product)

ndvi <- ToMNDWI(l)

Land Surface Water Index (LSWI)

The LSWI formula is:

$LSWI\quad =\quad \frac { \rho NIR\quad -\quad \rho SWIR1 }{ \rho NIR\quad +\quad \rho SWIR1 }$

Example

library(rLandsat8)

setwd("~/Downloads")
product  <- "LC82040322013219LGN00"
l <- ReadLandsat8(product)

ndvi <- ToLSWI(l)

The Normalized Difference Vegetation Index (NDVI), Modified Normalized Difference Water Index (MNDWI) and Land Surface Water Index (LSWI) need to compute the TOA Reflectance, so you can choose to use the sun angle correction, if sun angle correction is required, specifying is.suncorrected = TRUE (default value is is.suncorrected = FALSE)

Example

library(rLandsat8)

setwd("~/Downloads")
product  <- "LC82040322013219LGN00"
l <- ReadLandsat8(product, is.suncorrected = TRUE)

ndvi <- ToMNDWI(l)