# NOT RUN {
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)")
# }
Run the code above in your browser using DataLab