Vector geometry invalidity/frailty exposed as unrelated bugs fixed in GEOS and GDAL

[rsbivand]

Ongoing changes in vector geometry validity handling in GEOS and GDAL drivers used in R packages sf, rgdal and rgeos have led to unexpected changes or failures in R workflows and packages using these packages. Operations which previously appeared to succeed now fail, either noisily (good) or silently (bad).

• Invite users/developers to contribute cases

• Add other examples and links.

Shapefile affected: ESRI Shapefile driver not roundtripping in edge case

Intersection affected: GEOS 3.7.1 to 3.7.2 tightening of validity needed for operations

Plotting in tmap affected: tm_graticules() does not work for some cases

lidR problem with class inheritance

MODIS getTile() sp subset or raster::crop issues with PROJ 6.1 and GDAL 3.0

inlmisc Possible GEOS version vulnerability (3.7.2)

• Provide recommendations for checking/testing driver/GEOS versions in use and expected geometry handling outcomes

This issue is based on https://github.com/rsbivand/geostat19_talk; https://rsbivand.github.io/geostat19_talk/Discuss_issue_drafts.html