tura: NDVI time series raster brick

Description Usage Source References Examples

Description

Time series raster brick for an area in Tura kebele, Kafa Zone, SW Ethiopia. Values are NDVI rescaled by a factor of 10000. Data originate from the Landsat 5 TM and Landsat 7 ETM+ sensors. Original Landsat scene names can be found by typing names(tura). Dates are also contained in the z-dimension (getZ(tura); see also timeStack).

Usage

1

Source

http://earthexplorer.usgs.gov

References

DeVries B., Verbesselt J., Kooistra L. and Herold M. (2014). Robust Monitoring of Small-Scale Forest Disturbances in an Afromontane Forest Using Landsat Time Series. Remote Sensing of Environment, in review.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
data(tura)

# scene names
names(tura)
s <- getSceneinfo(names(tura))
s$year <- as.numeric(format(s$date, format = "%Y"))
hist(s$year, breaks = c(1980:2015))

# plot first 9 scenes
plot(tura, 1:9)
plot(tura, 1:9, main = getZ(tura))

# plot a time series of the 50th cell
plot(getZ(tura), tura[50], xlab="time", ylab="NDVI (x 10000)")

# make a new rasterBrick from only ETM+ scenes
x <- subsetRasterTS(tura, sensor = "ETM+")

# plot time series from the revised brick
plot(getZ(x), x[50], xlab="time", ylab="NDVI (x 10000)")

loicdtx/bfastSpatial documentation built on May 21, 2019, 7:37 a.m.