# This dataset is essentially constructed via:
nGrayLevels <- 15
require(raster)
ndvi <- raster(StrucDiv::ndvi)
ndvi15 <- cut(ndvi, breaks=seq(minValue(ndvi), maxValue(ndvi), len=nGrayLevels + 1),
include.lowest=TRUE, right=FALSE)
Run the code above in your browser using DataLab