```r
dem <- get_elev_raster(corfe_common_sssi, z = 14,clip="bbox")

contours<-rasterToContour(dem,nlevels=20)
sloper<-terrain(dem, opt='slope', unit='radians')
aspectr<-terrain(dem, opt='aspect',unit='radians')
hillshade<-hillShade(sloper,aspectr,angle=5)
```