Hello Ben, I think you should take a look on the raster package, there you have the TRI as a function Best, Jorge -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/topographic-roughness-tp7127002p7127110.html Sent from the R-sig-geo mailing list archive at Nabble.com.