Shapely-1.1a4-py2.5 - Failed to evaluate GEOSIntersects_r

[Mario Ceresa]

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
_______________________________________________
Community mailing list
Community@...
http://lists.gispython.org/mailman/listinfo/community

[Sean Gillies-2]

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
Community@...
http://lists.gispython.org/mailman/listinfo/community

[Mario Ceresa]

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 <sean.gillies@...> 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
> Community@...
> http://lists.gispython.org/mailman/listinfo/community
>

_______________________________________________
Community mailing list
Community@...
http://lists.gispython.org/mailman/listinfo/community

[Sean Gillies-2]

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  
> <sean.gillies@...> 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
>> Community@...
>> http://lists.gispython.org/mailman/listinfo/community
>>
> <prepared.py>_______________________________________________
> Community mailing list
> Community@...
> http://lists.gispython.org/mailman/listinfo/community

--
Sean

_______________________________________________
Community mailing list
Community@...
http://lists.gispython.org/mailman/listinfo/community

[Mario Ceresa]

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 <sean.gillies@...> 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
>> <sean.gillies@...> 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
>>> Community@...
>>> http://lists.gispython.org/mailman/listinfo/community
>>>
>> <prepared.py>_______________________________________________
>> Community mailing list
>> Community@...
>> http://lists.gispython.org/mailman/listinfo/community
>
> --
> Sean
>
> _______________________________________________
> Community mailing list
> Community@...
> http://lists.gispython.org/mailman/listinfo/community
>

_______________________________________________
Community mailing list
Community@...
http://lists.gispython.org/mailman/listinfo/community