# NOT RUN {
{
require(raster)
# read raster data
r <- raster(system.file("extdata", "ndvi.tif", package="fieldRS")) > 2000
r[r == 0] <- NA
# filter image (3x3 erosion)
or <- pixFilter(r, 1, "erode")
plot(r)
plot(or)
#' # filter image (3x3 dilation)
or <- pixFilter(r, 1, "dilate")
plot(r)
plot(or)
}
# }
Run the code above in your browser using DataLab