Guest User

Untitled

a guest
Jan 17th, 2019
76
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.18 KB | None | 0 0
  1. # terrain attributes computed with SAGA-GIS
  2. library(RSAGA)
  3. # set up rsaga environment
  4. work_env <- rsaga.env(path = "path/saga-6.3.0_x64/saga-6.3.0_x64", modules = "path/saga-6.3.0_x64/saga-6.3.0_x64/tools")
  5.  
  6. rsaga.get.libraries(path = work_env$modules)
  7. rsaga.get.modules(env = work_env,'grid_filter')
  8. rsaga.get.modules(libs = "io_gdal", env = work_env)
  9. rsaga.get.usage(lib = "io_gdal", module = "Export Raster", env = work_env)
  10. #rsaga.geoprocessor
  11.  
  12. rsaga.slope(in.dem = "path\dem.tif",
  13. out.slope = "path\slope",
  14. unit.slope = "degrees",
  15. out.aspect = "path\aspect",
  16. unit.aspect = "degrees",
  17. out.cgene = "path\general_curvature",
  18. out.clong = "path\longitudinal_curvature",
  19. out.ccros = "path\crosssectional_curvature",
  20. env = work_env)
  21.  
  22. Warning message:
  23. running command '"path/saga-6.3.0_x64/saga-6.3.0_x64/saga_cmd.exe"
  24. -f=q ta_morphometry "Slope, Aspect, Curvature" -ELEVATION
  25. "pathdem.tif" -SLOPE "path/slope.sgrd" -ASPECT
  26. "path/aspect.sgrd" -C_GENE "path/general_curvature.sgrd" -C_LONG
  27. "path/longitudinal_curvature.sgrd" -C_CROS
  28. "path/crosssectional_curvature.sgrd" -METHOD "6" -UNIT_SLOPE "1"
  29. -UNIT_ASPECT "1"' had status 5
Add Comment
Please, Sign In to add comment