Edzer Pebesma
As this is being discussed as a feature request for `sf`, can we think of an easy way of doing this without introducing a dependency on `spdep`?
> is good enough to use st_distance and choose neighbours as polygons closer than snap. IIRC, `st_is_within_distance()` is an optimised function for that, in any case fast in `s2`.
My guess: `buffer` computes a geometry on a sphere, `spin` rotates it in R2, in a Cartesian space, not on a sphere.
It got rotated in a Cartesian space.
```r library(terra) # terra 1.7.78 library(sf) # Linking to GEOS 3.12.1, GDAL 3.8.4, PROJ 9.3.1; sf_use_s2() is TRUE xlon sf::st_geometry() -> x x[[1]] # POLYGON ((8.660254 -95, 68.66025 8.923048, -51.33975...
This seems to now work: ```r library(stars) hu
> The thing is that read_mdim() doesn't read multiple files as time series. Can you give an example where that doesn't work?
Yes, but `read_mdim(file_list)` also works, even with `proxy=TRUE`. I'll look into you other issue.
```r hu hu st_crop(hu, st_bbox(c(xmin = 20, ymin = 25, xmax = 40, ymax = 50))) |> st_as_stars() stars object with 4 dimensions and 1 attribute attribute(s), summary of first...
> With read_stars(): Shouldn't be dropped dimension with only one level? Arguably no, because you loose information (the value of that single level). You can put an `adrop()` in the...