dem= rast(paste(system.file("extdata", package = "SurfRough"), "/trento1.tif",sep=""))
# iqr of slope in degrees
slope=terrain(dem, v="slope")
w=KernelCircular(3)
w
iqrSlope=iqrST(slope,w=w)
plot(iqrSlope)
Run the code above in your browser using DataLab