Return data for points and multi points

[jaseeverett]

@clairedavies and I regularly have need to extract data for points in the ocean. For example at mooring locations, CPR data points etc. It seems that oceandatr doesn't support the use of points. I understand this is slightly out of scope for you original intention, but it would certainly open up the package to new users.

For example we would like to get bathymetry data for multiple locations.

library(tidyverse)

dat = data.frame(lon = c(156, 158, 160, 162), lat = -42) %>% 
  sf::st_as_sf(coords = c("lon", "lat"), crs = "EPSG:4326")

out <- oceandatr::get_bathymetry(spatial_grid = dat, raw = TRUE)
#> This may take seconds to minutes, depending on grid size
#> Error in h(simpleError(msg, call)): error in evaluating the argument 'x' in selecting a method for function 'mask': [crop] extents do not overlap

^{Created on 2024-07-03 with reprex v2.1.0}

Is it possible to return the nearest neighbour for point data.

[jflowernet]

I like this idea. As you said, could increase the usefulness of the package. It would be good to have this work for any data request, i.e. not just bathymetry, so I will need to make changes to spatialgridr for this to work consistently.

I'm going to be working though issues over the coming weeks, so will try to get to this.

[clairedavies]

Thanks Jason, we would also potentially use it for other datasets too if it worked with some of the other functions you have, i.e. distance from shore and ecoregions.