Edzer Pebesma
Not sure whether this is meant by "extension types", but on the vector side, an extension of 2-D simple features as they're used by GEOS and GDAL, is the `POLYGON...
Since you closed the issue, did you solve it?
Thanks for the clarification!
I can look into this if you provide a reproducible example.
``` singleSide: logical; if ‘TRUE’, single-sided buffers are returned for linear geometries, in which case negative ‘dist’ values give buffers on the right-hand side, positive on the left; see details...
It seems that the algorithm fails if the buffer self-intersects; this works: ```r library(sf) # Linking to GEOS 3.12.1, GDAL 3.8.4, PROJ 9.3.1; sf_use_s2() is TRUE nc = st_read(system.file("shape/nc.shp", package...
A possible cause outside the package code may be that `geos` uses another version of `GEOS` (3.11.1) than `sf` (3.12.1). The 3.12.0 [release notes](https://github.com/libgeos/geos/blob/main/NEWS.md) mention work on single sided buffers.
Why don't you install the binary package from CRAN?
Maybe ask at r-sig-mac?
Thanks for the clear problem description! ```r st_extract(thetao[,,,1], dd, time_column = "yrmon") ``` gives you the correct values for the first depth, but doing this for all depths (your example)...