elevatr | An R package for accessing elevation data | Development Tools library

I hava another geodata question. I am trying to access the elevation data for specific points with the get_elev_point ()

I have about 160 points that I would like to get the elevation data for, so I would like to not do it one by one but I do not know how to retrieve them all at once.

I have made a data.frame of all my points where the first column x is longitude and the second column y is latitude and these are the only data in the data.frame.

I am using the elevatr package

Reproducible example:

location_1 <- data.frame (x=-7.37,y=5.775)

location_1_elev <-get_elev_point(locations = location_1, units="meters", prj = ll_prj, src = "aws")

When I do this, everything is fine and I get an elevation point back, but when I try to access several points at the same time I run into errors.

I took the earthquake data from R and transformed it into a data.frame that has only the longitude and latitiude. Then tried to access the elevation points via get_elev_points and got the error message:

data(quakes)

head(quakes)

locations <- data.frame(x = c(quakes$long, 1000), y = c(quakes$lat, 1000))

quakes_elev <-get_elev_point(locations = locations, units="meters", prj = ll_prj, src = "aws")

Error: API did not return tif

Do you have any tips how to make this happen, to be able to access multiple elevation points?

Thank you! ps.: Sorry for my clumsy asking, I am only learning now.

This is my first question so doing my best to make my question simple and reproducible.

I'm running into a problem with the R command geom_contour_fill from the metR library. The code below shows that geom_contour_fill fails to plot in some cases when geom_contour works just fine with the same input data. geom_contour_fill just returns a gray box in the plot pane.

However, by changing the zoom level to z = 9 in the get_elev_raster command in the code, geom_contour_fill works fine again.

I want to be able to use the higher zoom level and I've run into this problem randomly with other map extents, so can anyone see how to resolve the problem?

Thanks!

I am quite new to using R - and i am currently trying to create an elevation diagram of a study site using the below code: