# \donttest{
library(terra)
r <- rast(system.file("ex/elev.tif", package="terra"))
r <- r[[1]] / max(global(r, "max", na.rm=TRUE)[,1])
diss <- raster.mds(r)
diss.med <- raster.mds(r, window.median = TRUE)
opar <- par(no.readonly=TRUE)
par(mfrow=c(2,2))
plot(r)
title("Elevation")
plot( focal(r, w = matrix(1, nrow=5, ncol=5), fun = var) )
title("Variance")
plot(diss)
title("MDS")
plot(diss.med)
title("Median MDS")
par(opar)
# }
Run the code above in your browser using DataLab