```r
dorset_srtm<- raster::crop(srtm,dorset)

dslope<-terrain(dorset_srtm, opt='slope', unit='radians')
daspect<-terrain(dorset_srtm, opt='aspect',unit='radians')
dhill<-hillShade(dslope,daspect,angle=5)

plot(dhill, col = gray(0:100 / 100), legend = FALSE)
plot(dorset_srtm, col = terrain.colors(25), alpha = 0.4, legend = TRUE, add = TRUE) 
plot(dorset,add=TRUE)
```

<img src="OS_GIS_files/figure-html/unnamed-chunk-118-1.png" width="672" />