# NOT RUN {
# Generate points to download raster tiles for
set.seed(123)
simulated_data <- data.frame(
id = seq(1, 100, 1),
lat = runif(100, 44.1114, 44.1123),
lng = runif(100, -73.92273, -73.92147)
)
# Create an sf object from our original simulated data
simulated_data_sf <- sf::st_as_sf(simulated_data, coords = c("lng", "lat"))
sf::st_crs(simulated_data_sf) <- sf::st_crs(4326)
# Download data!
downloaded_tiles <- get_tiles(simulated_data_sf, tempfile())
merged_file <- merge_rasters(
downloaded_tiles[[1]],
tempfile(fileext = ".tif")
)
# Create an overlay image
vector_to_overlay(simulated_data_sf, merged_file[[1]], na.rm = TRUE)
# }
# NOT RUN {
# }
Run the code above in your browser using DataCamp Workspace