We can add additional predicates to STRtree query to enable better parity of spatial indexes in geopandas (based on this) and sf in R.

Crossref: geopandas/geopandas#1271 (comment)

These should be easy to add:

• covers
• covered by
• contains properly

May require a bit more care:

• equals / equals exact (in part due to extra tolerance parameter for equals exact)

Would require more logic to invert the results:

• disjoint (internally would run intersects, then would return these subtracted from the tree indexes)