trig.trans(slope = 48.146, aspect = 360.000)
library(raster)
elev <- getData('alt', country='CHE')
sa <- terrain(elev, opt=c("slope", "aspect"), unit="degrees")
scosa <- overlay(sa[[1]], sa[[2]], fun = trig.trans)
Run the code above in your browser using DataLab