```r
lon<--99.756
lat <-19.1

pnt<-data.frame(id=1,x=lon,y=lat)
pnt<-st_as_sf(pnt,coords=c("x","y"),crs=4326)
pnt %>% st_transform(3857)%>% st_buffer(5000) -> circle
elev <- get_elev_raster(circle, z = 12,clip="bbox")
sloper<-terrain(elev, opt='slope', unit='radians')
aspectr<-terrain(elev, opt='aspect',unit='radians')
hillshade<-hillShade(sloper,aspectr,angle=5)
plot(hillshade, col = gray(0:100 / 100), legend = FALSE)
plot(elev, col = terrain.colors(25), alpha = 0.3, legend = TRUE, add = TRUE)
```

<img src="006_terrain_analysis_files/figure-html/unnamed-chunk-9-1.png" width="672" />

```r
writeRaster(elev,file="rasters/volcanoes/nevado_toluca.tif",overwrite=TRUE)
```