RStoolbox (version 0.1.10)

tasseledCap: Tasseled Cap Transformation

Description

Calculates brightness, greenness and wetness from multispectral imagery. Currently implemented Landsat 4 TM, Landsat 5 TM, Landsat 7ETM+, Landsat 8 OLI and MODIS.

Usage

tasseledCap(img, sat, ...)

Arguments

img

RasterBrick or RasterStack. Input image. Band order must correspond to sensor specifications (see Details and Examples)

sat

Character. Sensor; one of: c("Landsat4TM", "Landsat5TM", "Landsat7ETM", "Landsat8OLI", "MODIS"). Case is irrelevant.

...

Further arguments passed to writeRaster.

Value

Returns a RasterBrick with the thee bands: brigthness, greenness, and (soil) wetness.

Details

Currently implemented: Landsat 4 TM, Landsat 5 TM, Landsat 7ETM+, Landsat 8 OLI and MODIS. Input data must be in top of atmosphere reflectance. Moreover, bands must be provided in ascending order as listed in the table below. Irrelevant bands, such as Landsat Thermal Bands must be omitted. Required bands are:

sat	bands	coefficients	data unit
Landsat4TM	1,2,3,4,5,7	Crist 1985	reflectance
Landsat5TM	1,2,3,4,5,7	Crist 1985	reflectance
Landsat7ETM	1,2,3,4,5,7	Huang 2002	reflectance
Landsat8OLI	2,3,4,5,6,7	Baig 2014	reflectance
MODIS	1,2,3,4,5,6,7	Lobser 2007	reflectance

References

Crist (1985) "A TM Tasseled Cap Equivalent Transformation for Reflectance Factor Data." Remote Sensing of Environment 17 (3): 301-306

Huang et al. (2002) "Derivation of a Tasselled Cap Transformation Based on Landsat 7 At-Satellite Reflectance." International Journal of Remote Sensing 23 (8): 1741-1748

Baig et al. (2014) "Derivation of a Tasselled Cap Transformation Based on Landsat 8 At-Satellite Reflectance." Remote Sensing Letters 5 (5): 423-431.

Lobser et al. (2007) "MODIS Tasselled Cap: Land Cover Characteristics Expressed through Transformed MODIS Data." International Journal of Remote Sensing 28 (22): 5079-5101.

Examples

# NOT RUN {
library(raster)
data(lsat)

## Run tasseled cap (exclude thermal band 6)
lsat_tc <- tasseledCap(lsat[[c(1:5,7)]], sat = "Landsat5TM")
lsat_tc
plot(lsat_tc)
# }