Edzer Pebesma
There are still several tests failing for me on this PR: ``` R CMD build sf --no-build-vignettes --resave-data * checking for file ‘sf/DESCRIPTION’ ... OK * preparing ‘sf’: * checking...
The shapefile open options are found [here](https://gdal.org/drivers/vector/shapefile.html#open-options).
The `st_read` docs say: `option: driver dependent dataset open options`. I pointed you to these options; they do not include `-dialect`.
... then please correct the issue title.
> without spatialite we will can't read gpkg files without spatialite GDAL can read gpkg files; the only dependency is sqlite3. The [ogr sqlite dialect](https://gdal.org/user/sql_sqlite_dialect.html#sql-sqlite-dialect) comes "preferrably with spatialite support",...
Installing sf from compiled GDAL, PROJ and GEOS in Red Hat Enterprise Linux RHEL 8.9 (no admin user)
Congrats on getting this done, and I'd agree that this is a very hard and painful procedure.
Thanks for the interesting benchmarks; most of this is well known though. I do consider (but do not promise) to work on the points-only case (see e.g. the pointx branch)...
Yes, but the number of polygons returned is the same. Maybe some grid cells touching the pole become triangles in geographical coordinates, and loose one node, hence the smaller number...
Please modify your reproducible example such that it exposes the error; `geneMetaUsed` is not known.