On Tue, Mar 15, 2011 at 2:19 AM, Alejandro Betancourt <[email protected]> wrote: > Hello, I'm Alejandro Betancourt, mathematical engineer and student of the > master in applied mathematics. Now i'm working on my tesis and i'm using > "shapely" in order to handle geometric operations between polygons. > I'm writing to you in order to report a possible bug in shapely. I'm getting > the bug when i tried to get the difference between two simple polygons > (attached). > The polygons are attached and the lines bellow are the way to get the error. >>>> import cPickle >>>> p1,p2 = cPickle.load(open("polygon.cp","rb")) >>>> p1.difference(p2) > Traceback (most recent call last): > File "<stdin>", line 1, in <module> > File "/usr/lib/pymodules/python2.6/shapely/geometry/base.py", line 314, in > difference > return geom_factory(self.impl['difference'](self, other)) > File "/usr/lib/pymodules/python2.6/shapely/topology.py", line 55, in > __call__ > "This operation produced a null geometry. Reason: unknown") > shapely.geos.TopologicalError: This operation produced a null geometry. > Reason: unknown > I hope this helps you in order to make better shapely.
Hi Alejandro, I apologize for being late with a reply. If you will you reattach the polygons in WKT format, I will be happy to take a look. Cheers, -- Sean _______________________________________________ Community mailing list [email protected] http://lists.gispython.org/mailman/listinfo/community
