playing with elevatr public AWS data not authed

playing_with_elevatr.R

library(rgdal)
library(elevatr)

# this is the prjection system given as example in the documentation
ll_prj <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0"

# Im using this projection system based on reading this document:
# https://www.nceas.ucsb.edu/~frazier/RSpatialGuides/OverviewCoordinateReferenceSystems.pdf
us_prj = "+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0"

# SpatialPoints is from 'rgdal' package
# this is centered on minneapolis
mans_sp <- SpatialPoints(coordinates(data.frame(y=-93.289865, x=44.965937)), CRS(us_prj))

# get raster data, Z is the zoom level
elevation <- get_elev_raster(mans_sp, z = 12, src = "aws")

# plot the raster data
plot(elevation)