# NOT RUN {
library(raster)
data(elev)
# calculate tpi and plot
tpi7 <- tpi(elev, scale=7)
tpi025 <- tpi(elev, win = "circle", scale=0.025)
tpi025.zc <- tpi(elev, win = "circle", scale=0.025,
zero.correct = TRUE)
par(mfrow=c(2,2))
plot(elev, main="original raster")
plot(tpi7, main="tpi 7x7")
plot(tpi025, main="tpi Circular window d=0.025")
plot(tpi025, main="tpi Circular window d=0.025, zero correct")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab