Edzer Pebesma

Follow-up email to gdal-dev: https://lists.osgeo.org/pipermail/gdal-dev/2022-August/056145.html

We point out [here](https://keen-swartz-3146c4.netlify.app/spherical.html#validity-on-the-sphere) that there are cases that can't be valid on both R2 and S2. Another one is Fiji, which is cut in two. The admin boundaries distributed...

This now terminates, but requires the user to understand the problem. ```r options(s2_oriented=TRUE) library(stars) # Loading required package: abind # Loading required package: sf # Linking to GEOS 3.10.2, GDAL...

Works without problems here, also in rstudio (ubuntu 20.04, RStudio 2022.07.1+554), also no issues when run under valgrind. What is your `sessionInfo()` after loading `sf`?

`terra` copied some C++ code from `sf` for reading/writing OGR data, so it might be in there, or it is in GDAL (but I'm running GDAL 3.4.3 too) or in...

Great; I'm looking at @dblodgett-usgs stack trace for sf, and am completely puzzled why `CPL_get_z_range` seems to get called from within `sf_from_ogrlayer`. There's no reason for it, but there isn't...

@rhijmans maybe you could take a look what is going on here, being familiar with windows debugging?

What was returned before, and before what or when? If anything this would be a change in GEOS, which maybe returned a `GEOMETRYCOLLECTION EMPTY` before?

the [AAIGrid](https://gdal.org/drivers/raster/aaigrid.html) GDAL driver doesn't support creation; I also don't think it is a useful data format for creation of anything serious; I would consider a PR but will not...

That is a well-known property of the Gaussian variogram - adding a small nugget often helps.