# \donttest{
library(raster)
r <- raster(nrows=500, ncols=500, xmn=571823, xmx=616763,
ymn=4423540, ymx=4453690)
r[] <- runif(ncell(r), 1, 100)
r <- focal(r, focalWeight(r, 150, "Gauss") )
r.inv <- raster.invert(r)
opar <- par(no.readonly=TRUE)
par(mfrow=c(1,2))
plot(r, main="original raster")
plot(r.inv, main="inverted raster")
par(opar)
# }
Run the code above in your browser using DataLab