Hey Sean, I just rebuilt geos from Paul's 3.1.0rc1.tar.bz2 and was excited to play around with the cascaded union with a big watershed dataset that's always been a bear.
I am on mac 10.5.6 and used your script from: http://sgillies.net/blog/870/a-more-perfect-union-continued It seems to hit the last feature and then crashes python: [....] cxx type MultiPolygon Assertion failed: (!"should never be reached"), function itemsTree, file AbstractSTRtree.cpp, line 358. Abort trap Then I tried on Bjorn's border shapefile (v3) (TM_WORLD_BORDERS-0.3) and also on the last feature it fails, but with this exception: Traceback (most recent call last): File "<stdin>", line 1, in <module> File "/Library/Python/2.5/site-packages/Shapely-1.1.0-py2.5.egg/ shapely/ops.py", line 48, in cascaded_union return geom_factory(lgeos.GEOSUnionCascaded(collection)) File "/Library/Python/2.5/site-packages/Shapely-1.1.0-py2.5.egg/ shapely/geometry/base.py", line 35, in geom_factory raise ValueError, "No Shapely geometry can be created from this null value" ValueError: No Shapely geometry can be created from this null value Seems odd huh? Let me know if you have any ideas for further testing. Thanks! Dane On Jan 26, 2009, at 11:11 AM, Sean Gillies wrote: > So far so good. > > http://sgillies.net/blog/869/efficient-batch-operations-for-shapely > > Anyone want to join me as a downstream user and tester of GEOS 3.1? > > -- > Sean Gillies > [email protected] > http://sgillies.net > > _______________________________________________ > Community mailing list > [email protected] > http://lists.gispython.org/mailman/listinfo/community _______________________________________________ Community mailing list [email protected] http://lists.gispython.org/mailman/listinfo/community
