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... -brentp _______________________________________________ Community mailing list [email protected] http://lists.gispython.org/mailman/listinfo/community
