Hello Sean, > Well, I was unable to reproduce because I was getting bus errors ;)
Aah... I see :) > In case you missed some recent project news: what I've been calling [...] Yes I received the emails, and noticed that probably in 1.2 there will be prepared geometries :) I was trying 1.1 just to give them a try! Cheers, Mario On Fri, Sep 18, 2009 at 1:44 PM, Sean Gillies <[email protected]> wrote: > Well, I was unable to reproduce because I was getting bus errors ;) > > In case you missed some recent project news: what I've been calling > 1.1 is going to be 2.0. Aron Bierbaum is leading the work on a 1.2 > version that will make the most of GEOS 3.1 features while remaining > backwards compatible. 2.0 is going to be late because I'm tied up with > other work at ISAW for a few more weeks. > > Cheers, > > On Sep 18, 2009, at 1:28 PM, Mario Ceresa wrote: > >> Hello Sean, >> thanks for your mail. Here I can reproduce every time. Sometimes when >> I try with the attached code, instead of throwing an exception, it >> segfaults in: >> >> #0 GEOSGeomTypeId_r (extHandle=0x21c2530, g1=0x1) at geos_ts_c.cpp: >> 1191 >> 1191 return g1->getGeometryTypeId(); >> >> If you need more information I'll be more than happy to help! >> >> A presto, >> >> Mario >> >> On Fri, Sep 18, 2009 at 12:59 PM, Sean Gillies >> <[email protected]> wrote: >>> >>> On Sep 17, 2009, at 6:48 PM, Mario Ceresa wrote: >>> >>>> Hi everybody! >>>> I was curious to see new Shapely 1.1 in action and installed it in a >>>> virtualenv as explained in >>>> >>>> http://sgillies.net/blog/900/diving-into-shapely-1-1/ >>>> >>>> but when I run the following code: >>>> >>>> from shapely.geometry import Point, Polygon >>>> from random import random >>>> spots = [Point(random()*2.0-0.5, random()*2.0-1.0).buffer(0.1) for i >>>> in xrange(200)] >>>> triangle = Polygon(((0.0, 0.0), (1.0, 1.0), (1.0, -1.0))) >>>> x = [s for s in spots if triangle.intersects(s)] >>>> >>>> I get: >>>> >>>> Traceback (most recent call last): >>>> File "prepared.py", line 5, in <module> >>>> x = [s for s in spots if triangle.intersects(s)] >>>> File "Download/try-shapely/lib/python2.5/site-packages/ >>>> Shapely-1.1a4-py2.5.egg/shapely/geometry/base.py", >>>> line 191, in intersects >>>> return self._binaryGeometryPredicates.intersects(self, other) >>>> File "Download/try-shapely/lib/python2.5/site-packages/ >>>> shapely.geos-1.0a4-py2.5.egg/shapely/geos/predicates.py", >>>> line 63, in __call__ >>>> raise PredicateError, "Failed to evaluate %s in context %s and >>>> %s" >>>> % (self.func.__name__, repr(context), repr(other)) >>>> shapely.geos.PredicateError: Failed to evaluate GEOSIntersects_r in >>>> context <shapely.geometry.polygon.Polygon object at 0x19e2310> and >>>> <shapely.geometry.polygon.Polygon object at 0x19e25d0> >>>> Exception exceptions.AttributeError: "'NoneType' object has no >>>> attribute 'GEOSWKBWriter_destroy'" in <bound method >>>> GEOSWKBWriter.__del__ of <shapely.geos.wkb.GEOSWKBWriter object at >>>> 0x1a63750>> ignored >>>> Exception exceptions.AttributeError: "'NoneType' object has no >>>> attribute 'GEOSWKBReader_destroy'" in <bound method >>>> GEOSWKBReader.__del__ of <shapely.geos.wkb.GEOSWKBReader object at >>>> 0x1a636d0>> ignored >>>> >>>> Any idea on what I'm doing wrong? I use Geos 3.1.0 on a Fedora 10 >>>> x64 system. >>>> >>>> Thanks and regards, >>>> >>>> Mario >>> >>> Hi Mario, >>> >>> GEOS is complaining that it can't intersect two geometries. It's >>> likely that there's a bug in shapely.geos, either in the code that >>> computes the GEOS geometry from the Python coordinate sequence >>> arrays, >>> or in the caching of the geometry. >>> >>> I'm unable to reproduce the bug right now, but will look into it. >>> >>> Cheers, >>> >>> -- >>> Sean >>> >>> _______________________________________________ >>> Community mailing list >>> [email protected] >>> http://lists.gispython.org/mailman/listinfo/community >>> >> <prepared.py>_______________________________________________ >> Community mailing list >> [email protected] >> http://lists.gispython.org/mailman/listinfo/community > > -- > Sean > > _______________________________________________ > Community mailing list > [email protected] > http://lists.gispython.org/mailman/listinfo/community > _______________________________________________ Community mailing list [email protected] http://lists.gispython.org/mailman/listinfo/community
