# NOT RUN {
{
require(raster)
require(fieldRS)
# read raster data
r <- brick(system.file("extdata", "ndvi.tif", package="fieldRS"))
# read field data
data(fieldData)
data(fieldDataTS)
a.ts <- analyseTS(as.data.frame(fieldDataTS$weighted.mean), fieldData$crop)
# extract reference profiles
rp <- as.data.frame(do.call(rbind, lapply(a.ts$y.statistics, function(i) {i$median})))
# compare labels
cl <- compareLabel(as.data.frame(fieldDataTS$weighted.mean), rp, fieldData$crop, a.ts$labels)
}
# }
Run the code above in your browser using DataLab