gdallocationinfo

The goal of gdallocationinfo is to look up values on a raster source.

Input is a character string to the GDAL readable source, and matrix of lon,lat.

Installation

You can install the development version of gdallocationinfo like so:

devtools::install_github("hypertidy/gdallocationinfo")

Example

This is a basic example, lookup GEBCO 2023 values on world cities:

library(gdallocationinfo)
#options(parallelly.fork.enable = TRUE, future.rng.onMisuse = "ignore")
#library(furrr); plan(multicore)

loc <- locationinfo("/vsicurl/https://gebco2023.s3.valeria.science/gebco_2023_land_cog.tif", 
                    as.matrix(maps::world.cities[1:1000, c("long", "lat")]))
#> 
#> Attaching package: 'purrr'
#> The following object is masked from 'package:base':
#> 
#>     %||%

#plan(sequential)

str(loc)
#>  num [1:1000] 79 86 140 17 331 ...
range(loc)
#> [1] -365 3848
maps::world.cities$name[which.max(loc)]
#> [1] "Achacachi"
maps::world.cities$name[which.min(loc)]
#> [1] "Ahe"

Code of Conduct

Please note that the gdallocationinfo project is released with a Contributor Code of Conduct. By contributing to this project, you agree to abide by its terms.