if (FALSE) {
library("RSAGA")
data(dem)
# Print the DEM header:
dem$header
# Write the DEM to a SAGA grid:
write.sgrd(data = dem, file = "dem", header = dem$header, env = env)
# Calculate slope of DEM:
rsaga.slope(in.dem = "dem", out.slope = "slope", method = "poly2zevenbergen")
}
Run the code above in your browser using DataLab