Edzer Pebesma

See also https://twitter.com/paleolimbot/status/1315088110234087424

There's an unexported function that could be used ```r > stars:::colrow_from_xy(st_coordinates(x), s) [,1] [,2] 1 2 3 2 2 3 3 5 1 4 1 4 ``` and then combined,...

Great reprex, btw @btupper !

gives us ```r library(stars) # Loading required package: abind # Loading required package: sf # Linking to GEOS 3.10.2, GDAL 3.4.3, PROJ 8.2.0; sf_use_s2() is TRUE library(tibble) nr = 5...

Interestingly I get  with ``` > sessionInfo() R version 4.2.1 (2022-06-23) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu 22.04.1 LTS Matrix products: default BLAS: /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3 LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.20.so locale: [1] LC_CTYPE=en_US.UTF-8...

@paleolimbot could you take a look at this? Is something maybe recycled once too often?

Thanks! Should we make `edge_type = "undirected"` the default, for all transformers? New & old values: ```r library(sf) # Linking to GEOS 3.10.2, GDAL 3.4.3, PROJ 8.2.0; sf_use_s2() is TRUE...

Yes, great idea. Maybe it is easier to let the NA in; the current implementation also reports the time resolution.

Will spatial also be the vertical dimension? I like to have the option for only x/y raster; getting the x/y raster dimensions is done by ```r > attr(st_dimensions(r), "raster")$dimensions [1]...

Thanks for trying! This indicates that you are trying to load a zarr file that has BLOSC compression, but your GDAL version used to compile your sf package was not...