DOC: Enhance description of STRtree to indicate that geometries are indexed and queried based on extents

Expected behavior and actual behavior.

STRTree returns a collision where there is none. A minimal case:

Steps to reproduce the problem.

from shapely import wkt
from shapely.strtree import STRtree

poly_query = wkt.loads('POLYGON ((58.27 44.72, 56.96 47.18, '
                       '58.47 47.98, 59.78 45.52, '
                       '58.27 44.72))')
poly_in_tree = wkt.loads('POLYGON ((57.94 44.73, 57.94 44.42, '
                         '56.02 44.42, 56.02 44.73, '
                         '57.94 44.73))')

query_tree = STRtree([poly_in_tree])
collide_objects = query_tree.query_geoms(poly_query)

assert not poly_in_tree.intersects(poly_query)
assert 0 == len(collide_objects)

Operating system

Windows

Shapely version and provenance

Shapely 1.8.2 via conda.

[brendan-ward]

axis-aligned-bounding-box

correct. We can try to make the description of STRtree more clear in this respect.

[adishavit]

I’m using Shapely 1.8, so I have to do the filtering manually.