Brent Pedersen wrote: > On Wed, Jun 25, 2008 at 8:52 AM, Jaakko Salli <[EMAIL PROTECTED]> wrote: >> Hi all, >> >> I have encountered strange problem with Polygon.contains(Point). It >> seems it stops functioning correctly after being called 21 times. At >> that point, instead of checking if point is inside the polygon, it just >> check if it is inside polygon's bounding box. >> >> Here is code to reproduce the issue (I have tried it with Python 2.5.2 >> and Shapely SVN trunk): >> >> from shapely.geometry import Point, Polygon >> >> # Form a 'L' shaped polygon >> points = [(0.0, 0.0), >> (1.0, 0.0), >> (1.0, 1.0), >> (2.0, 1.0), >> (2.0, 2.0), >> (0.0, 2.0), >> (0.0, 0.0)] >> >> poly = Polygon(points) >> >> # Create point that should be outside polygon, but >> # inside its bounding box >> pt_outside_but_inside_bbox = Point(1.5, 0.5) >> >> # Create point that should be outside polygon, >> # even outside its bounding box >> pt_outside_even_bbox = Point(2.5, 0.5) >> >> # At iteration number 22, poly.contains effectively becomes >> # poly_bounding_box.contains >> for i in range(1,100): >> res1 = poly.contains(pt_outside_but_inside_bbox) >> res2 = poly.contains(pt_outside_even_bbox) >> print '#%i: poly.contains(pt_outside_but_inside_bbox) = %s'%(i, res1) >> print ' poly.contains(pt_outside_even_bbox) = %s'%(res2) >> if res1 != False: >> print "poly.contains failed at try #%i!"%(i) >> break >> >> >> Thanks, >> Jaakko Salli >> >> >> _______________________________________________ >> Community mailing list >> [email protected] >> http://lists.gispython.org/mailman/listinfo/community >> > > thankfully, i dont see this with r1117. that scared me. > even upped the range to 100K and no problems. > fwiw, i'm using geos 3.0.0, not an rc... I'm using the geos DLL provided with the Windows installer (also using Windows XP SP2, if that matters).
Jaakko _______________________________________________ Community mailing list [email protected] http://lists.gispython.org/mailman/listinfo/community
