Sorry, I meant that, when you use Rtree, the "intersects" predicate is always True for each polygon obtained with the "intersection" method of Rtree.
So if two polygons don't intersect each other, if their "minimum bounding rectangle" intersects each other, the "intersects" predicate is True. 2009/3/20 Pascal Leroux <[email protected]> > Hi all, > > With the geos library version 3.1, it seems that the "intersects" predicate > (used with polygons) always returns True. For more than 135,000 couples of > polygons (got with the "intersection" method of Rtree), I always got True > but about 125,000 intersections were empty geometry collections. > > When I use the version 3.0.3, it doesn't happen : each time "intersects" > returns True, the intersection is not empty. > > I couldn't reproduce this issue in C with version 3.1 : when GEOSintersects > return 1 (True), the intersection returned by GEOSIntersection is not empty. > > All these tests have been done on Ubuntu with Shapely version 1.0.11 > > Here's a couple of polygons for which intersects return True with an empty > intersection > > POLYGON ((599801.1221290760440752 2546550.0634658345952630, > 599806.0258591863093898 2546538.0948937772773206, 599780.3041890829335898 > 2546526.9730210378766060, 599777.5322963981889188 2546499.6322209336794913, > 599807.2772084334865212 2546496.2815793519839644, 599804.4406587256817147 > 2546476.6445081932470202, 599760.1544660580111668 2546483.7799901803955436, > 599769.8198137159924954 2546536.2887050621211529, 599801.1221290760440752 > 2546550.0634658345952630)) > > POLYGON ((599758.8006132490700111 2546478.1623347853310406, > 599772.7430304150329903 2546474.0790138356387615, 599760.4723994010128081 > 2546410.1387482765130699, 599767.0911384929204360 2546408.3941536685451865, > 599762.7293910720618442 2546391.5474267113022506, 599740.1731945683714002 > 2546396.7629298437386751, 599758.8006132490700111 2546478.1623347853310406)) > > Pascal >
_______________________________________________ Community mailing list [email protected] http://lists.gispython.org/mailman/listinfo/community
