shapely 1.0.3__tar.gz → 1.0.5__tar.gz

Shapely-1.0.5/CHANGES.txt

@@ -0,0 +1,24 @@
+All tickets are children of http://trac.gispython.org/projects/PCL/ticket/
+
+1.0.5 (2008-05-20)
+------------------
+- Added access to GEOS polygonizer function.
+- Raise exception when insufficient coordinate tuples are passed to LinearRing
+  constructor (#164).
+
+1.0.4 (2008-05-01)
+------------------
+- Disentangle Python and topological equality (#163).
+- Add shape(), a factory that copies coordinates from a geo interface provider.
+  To be used instead of asShape() unless you really need to store coordinates
+  outside shapely for efficient use in other code.
+- Cache GEOS geometries in adapters (#163).
+
+1.0.3 (2008-04-09)
+------------------
+- Do not release GIL when calling GEOS functions (#158).
+- Prevent faults when chaining multiple GEOS operators (#159).
+
+1.0.2 (2008-02-26)
+------------------
+- Fix loss of dimensionality in polygon rings (#155).

{Shapely-1.0.3 → Shapely-1.0.5}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 1.0
 Name: Shapely
-Version: 1.0.3
+Version: 1.0.5
 Summary: Geospatial geometries, predicates, and operations
 Home-page: http://trac.gispython.org/projects/PCL/wiki/Shapely
 Author: Sean Gillies

{Shapely-1.0.3 → Shapely-1.0.5}/Shapely.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 1.0
 Name: Shapely
-Version: 1.0.3
+Version: 1.0.5
 Summary: Geospatial geometries, predicates, and operations
 Home-page: http://trac.gispython.org/projects/PCL/wiki/Shapely
 Author: Sean Gillies

{Shapely-1.0.3 → Shapely-1.0.5}/Shapely.egg-info/SOURCES.txt

@@ -18,6 +18,7 @@ shapely/__init__.py
 shapely/geos.py
 shapely/iterops.py
 shapely/kml.py
+shapely/ops.py
 shapely/predicates.py
 shapely/topology.py
 shapely/wkb.py
@@ -32,6 +33,7 @@ shapely/geometry/multipoint.py
 shapely/geometry/multipolygon.py
 shapely/geometry/point.py
 shapely/geometry/polygon.py
+shapely/geometry/proxy.py
 tests/Array.txt
 tests/GeoInterface.txt
 tests/IterOps.txt
@@ -48,6 +50,7 @@ tests/__init__.py
 tests/attribute-chains.txt
 tests/binascii_hex.txt
 tests/invalid_intersection.txt
+tests/polygonize.txt
 tests/rungrind.dist
 tests/test_doctests.py
 tests/test_threads.py

{Shapely-1.0.3 → Shapely-1.0.5}/manual/manual.txt

@@ -4,16 +4,15 @@ The Shapely Manual
 
 :Author: Sean Gillies
 :address: sgillies@frii.com
-    http://zcologia.com/sgillies
-
-:revision: 1.0
-:date: 2 January 2008
-
-:copyright: This work is licensed under a `Creative Commons Attribution 3.0 United States License`__.
+:revision: 1.0.5
+:date: 20 May 2008
+:copyright: This work is licensed under a `Creative Commons Attribution 3.0
+  United States License`__.
 
 .. __: http://creativecommons.org/licenses/by/3.0/us/
 
-:abstract: This document describes the Shapely Python package for programming with geospatial geometries.
+:abstract: This document describes the Shapely Python package for programming
+  with geospatial geometries.
 
 .. sectnum::
 
@@ -342,6 +341,32 @@ not:
   <shapely.geometry.polygon.Polygon object at ...>
 
 
+Other Operations
+----------------
+
+Polygonization
+++++++++++++++
+
+shapely.ops.polygonize(lines) : iterator
+  Returns an iterator over polygons constructed from the *lines* iterator. The
+  elements of *lines* may be Shapely geometries, objects that provide the geo
+  interface, or Numpy arrays or Python sequences shaped like LineStrings.
+
+.. code-block:: python
+
+  >>> from shapely.ops import polygonize
+  >>> lines = [
+  ...     ((0, 0), (1, 1)),
+  ...     ((0, 0), (0, 1)),
+  ...     ((0, 1), (1, 1)),
+  ...     ((1, 1), (1, 0)),
+  ...     ((1, 0), (0, 0))
+  ...     ]
+  >>> result = polygonize(lines)
+  >>> list(result.geoms)
+  [<shapely.geometry.polygon.Polygon object at ...>, <shapely.geometry.polygon.Polygon object at ...>]
+
+
 Unary Spatial Predicates
 ------------------------
 
@@ -798,6 +823,9 @@ Or a simple placemark-type object:
   >>> list(shape.coords)
   [(0.0, 0.0)]
 
+If you want to copy coordinate data to a new geometry, use the
+*shapely.geometry.shape* function instead.
+
 .. _Python geo interface: http://trac.gispython.org/projects/PCL/wiki/PythonGeoInterface
 
 

{Shapely-1.0.3 → Shapely-1.0.5}/setup.py

@@ -11,7 +11,7 @@ if version_info[:2] < (2,5):
 readme_text = file('README.txt', 'rb').read()
 
 setup(name          = 'Shapely',
-      version       = '1.0.3',
+      version       = '1.0.5',
       description   = 'Geospatial geometries, predicates, and operations',
       license       = 'BSD',
       keywords      = 'geometry topology',

{Shapely-1.0.3 → Shapely-1.0.5}/setup_windows.py

@@ -11,7 +11,7 @@ if version_info[:2] < (2,5):
 readme_text = file('README.txt', 'rb').read()
 
 setup(name          = 'Shapely',
-      version       = '1.0.3',
+      version       = '1.0.5',
       description   = 'Geospatial geometries, predicates, and operations',
       license       = 'BSD',
       keywords      = 'geometry topology',

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/geometry/__init__.py

@@ -1,4 +1,4 @@
-from geo import asShape
+from geo import shape, asShape
 from point import Point, asPoint
 from linestring import LineString, asLineString
 from polygon import Polygon, asPolygon
@@ -6,4 +6,3 @@ from multipoint import MultiPoint, asMultiPoint
 from multilinestring import MultiLineString, asMultiLineString
 from multipolygon import MultiPolygon, asMultiPolygon
 from collection import GeometryCollection
-

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/geometry/base.py

@@ -207,11 +207,29 @@ class GeometrySequence(object):
                 max = l
 
 
+class HeterogeneousGeometrySequence(GeometrySequence):
+
+    def __init__(self, parent):
+        self.__p__ = parent
+        self._geom = parent._geom
+        self._ndim = parent._ndim
+
+    def next(self):
+        if self._index < self.__len__():
+            sub = lgeos.GEOSGetGeometryN(self._geom, self._index)
+            g = geom_factory(sub)
+            g._owned = True
+            self._index += 1
+            return g
+        else:
+            raise StopIteration 
+    
+
 def exceptNull(func):
     """Decorator which helps avoid GEOS operations on null pointers."""
     def wrapper(*args, **kwargs):
         # self is the first arg
-        if args[0]._geom is None:
+        if not args[0]._geom:
             raise ValueError, "Null geometry supports no operations"
         return func(*args, **kwargs)
     return wrapper
@@ -220,7 +238,7 @@ def exceptEitherNull(func):
     """Decorator which avoids GEOS operations on one or more null pointers."""
     def wrapper(*args, **kwargs):
         # self is the first arg
-        if args[0]._geom is None or args[1]._geom is None:
+        if not args[0]._geom or not args[1]._geom:
             raise ValueError, "Null geometry supports no operations"
         return func(*args, **kwargs)
     return wrapper
@@ -248,12 +266,6 @@ class BaseGeometry(object):
     def __str__(self):
         return self.to_wkt()
 
-    def __eq__(self, other):
-        return self.equals(other)
-    
-    def __ne__(self, other):
-        return not self.equals(other)
-
     # To support pickling
 
     def __reduce__(self):
@@ -465,3 +477,29 @@ class BaseGeometry(object):
         
         return (minx, miny, maxx, maxy)
 
+
+class BaseMultiPartGeometry(BaseGeometry):
+
+    @property
+    def ctypes(self):
+        raise NotImplementedError, \
+        "Multi-part geometries have no ctypes representations"
+
+    @property
+    def __array_interface__(self):
+        """Provide the Numpy array protocol."""
+        raise NotImplementedError, \
+        "Multi-part geometries do not themselves provide the array interface"
+
+    def _get_coords(self):
+        raise NotImplementedError, \
+        "Sub-geometries may have coordinate sequences, but collections do not"
+
+    def _set_coords(self, ob):
+        raise NotImplementedError, \
+        "Sub-geometries may have coordinate sequences, but collections do not"
+
+    @property
+    def coords(self):
+        raise NotImplementedError, \
+        "Multi-part geometries do not provide a coordinate sequence"

Shapely-1.0.5/shapely/geometry/collection.py

@@ -0,0 +1,38 @@
+"""
+Geometry collections.
+"""
+
+from shapely.geometry.base import BaseMultiPartGeometry
+from shapely.geometry.base import HeterogeneousGeometrySequence, exceptNull
+
+
+class GeometryCollection(BaseMultiPartGeometry):
+
+    """A geometry collection.
+    """
+
+    def __init__(self):
+        BaseMultiPartGeometry.__init__(self)
+
+    @property
+    def __geo_interface__(self):
+        geometries = []
+        for geom in self.geoms:
+            geometries.append(geom.__geo_interface__)
+        return dict(type='GeometryCollection', geometries=geometries)
+
+    @property
+    @exceptNull
+    def geoms(self):
+        return HeterogeneousGeometrySequence(self)
+
+
+# Test runner
+def _test():
+    import doctest
+    doctest.testmod()
+
+
+if __name__ == "__main__":
+    _test()
+

Shapely-1.0.5/shapely/geometry/geo.py

@@ -0,0 +1,64 @@
+"""
+Geometry factories based on the geo interface.
+"""
+
+from point import Point, asPoint
+from linestring import LineString, asLineString
+from polygon import Polygon, asPolygon
+from multipoint import MultiPoint, asMultiPoint
+from multilinestring import MultiLineString, asMultiLineString
+from multipolygon import MultiPolygon, asMultiPolygon
+
+
+def shape(context):
+    """Return a new, independent geometry with coordinates *copied* from the
+    context.
+    """
+    if hasattr(context, "__geo_interface__"):
+        ob = context.__geo_interface__
+    else:
+        ob = context
+    geom_type = ob.get("type").lower()
+    if geom_type == "point":
+        return Point(ob["coordinates"])
+    elif geom_type == "linestring":
+        return LineString(ob["coordinates"])
+    elif geom_type == "polygon":
+        return Polygon(ob["coordinates"][0], ob["coordinates"][1:])
+    elif geom_type == "multipoint":
+        return MultiPoint(ob["coordinates"])
+    elif geom_type == "multilinestring":
+        return MultiLineString(ob["coordinates"])
+    elif geom_type == "multipolygon":
+        return MultiPolygon(ob["coordinates"])
+    else:
+        raise ValueError, "Unknown geometry type: %s" % geom_type
+
+def asShape(context):
+    """Adapts the context to a geometry interface. The coordinates remain
+    stored in the context.
+    """
+    if hasattr(context, "__geo_interface__"):
+        ob = context.__geo_interface__
+    else:
+        ob = context
+
+    try:
+        geom_type = ob.get("type").lower()
+    except AttributeError:
+        raise ValueError, "Context does not provide geo interface"
+
+    if geom_type == "point":
+        return asPoint(ob["coordinates"])
+    elif geom_type == "linestring":
+        return asLineString(ob["coordinates"])
+    elif geom_type == "polygon":
+        return asPolygon(ob["coordinates"][0], ob["coordinates"][1:])
+    elif geom_type == "multipoint":
+        return asMultiPoint(ob["coordinates"])
+    elif geom_type == "multilinestring":
+        return asMultiLineString(ob["coordinates"])
+    elif geom_type == "multipolygon":
+        return asMultiPolygon(ob["coordinates"])
+    else:
+        raise ValueError, "Unknown geometry type: %s" % geom_type

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/geometry/linestring.py

@@ -7,6 +7,7 @@ from ctypes import ArgumentError
 
 from shapely.geos import lgeos
 from shapely.geometry.base import BaseGeometry, exceptNull
+from shapely.geometry.proxy import CachingGeometryProxy
 
 
 def geos_linestring_from_py(ob, update_geom=None, update_ndim=0):
@@ -15,8 +16,12 @@ def geos_linestring_from_py(ob, update_geom=None, update_ndim=0):
         array = ob.__array_interface__
         assert len(array['shape']) == 2
         m = array['shape'][0]
-        n = array['shape'][1]
-        assert m >= 2
+        if m < 2:
+            raise ValueError, "LineStrings must have at least 2 coordinate tuples"
+        try:
+            n = array['shape'][1]
+        except IndexError:
+            raise ValueError, "Input %s is the wrong shape for a LineString" % str(ob)
         assert n == 2 or n == 3
 
         # Make pointer to the coordinate array
@@ -41,8 +46,11 @@ def geos_linestring_from_py(ob, update_geom=None, update_ndim=0):
             dy = c_double(cp[n*i+1])
             dz = None
             if n == 3:
-                dz = c_double(cp[n*i+2])
-        
+                try:
+                    dz = c_double(cp[n*i+2])
+                except IndexError:
+                    raise ValueError, "Inconsistent coordinate dimensionality"
+
             # Because of a bug in the GEOS C API, 
             # always set X before Y
             lgeos.GEOSCoordSeq_setX(cs, i, dx)
@@ -53,8 +61,12 @@ def geos_linestring_from_py(ob, update_geom=None, update_ndim=0):
     except AttributeError:
         # Fall back on list
         m = len(ob)
-        n = len(ob[0])
-        assert m >= 2
+        if m < 2:
+            raise ValueError, "LineStrings must have at least 2 coordinate tuples"
+        try:
+            n = len(ob[0])
+        except TypeError:
+            raise ValueError, "Input %s is the wrong shape for a LineString" % str(ob)
         assert n == 2 or n == 3
 
         # Create a coordinate sequence
@@ -74,7 +86,10 @@ def geos_linestring_from_py(ob, update_geom=None, update_ndim=0):
             dy = c_double(coords[1])
             dz = None
             if n == 3:
-                dz = c_double(coords[2])
+                try:
+                    dz = c_double(coords[2])
+                except IndexError:
+                    raise ValueError, "Inconsistent coordinate dimensionality"
         
             # Because of a bug in the GEOS C API, 
             # always set X before Y
@@ -157,18 +172,18 @@ class LineString(BaseGeometry):
     coords = property(BaseGeometry._get_coords, _set_coords)
 
 
-class LineStringAdapter(LineString):
+class LineStringAdapter(CachingGeometryProxy, LineString):
 
     """Adapts a Python coordinate pair or a numpy array to the line string
     interface.
     """
     
     context = None
-    #__geom = None
     _owned = False
 
     def __init__(self, context):
         self.context = context
+        self.factory = geos_linestring_from_py
 
     @property
     def _ndim(self):
@@ -182,14 +197,6 @@ class LineStringAdapter(LineString):
             # Fall back on list
             return len(self.context[0])
 
-    @property
-    def _geom(self):
-        """Keeps the GEOS geometry in synch with the context."""
-        if self.__geom__ is not None:
-            lgeos.GEOSGeom_destroy(self.__geom__)
-        self.__geom__, n = geos_linestring_from_py(self.context)
-        return self.__geom__
-
     @property
     def __array_interface__(self):
         """Provide the Numpy array protocol."""
@@ -220,4 +227,3 @@ def _test():
 
 if __name__ == "__main__":
     _test()
-

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/geometry/multilinestring.py

@@ -7,6 +7,7 @@ from ctypes import byref, c_double, c_int, c_void_p, cast, POINTER, pointer
 from shapely.geos import lgeos
 from shapely.geometry.base import BaseGeometry, GeometrySequence, exceptNull
 from shapely.geometry.linestring import LineString, geos_linestring_from_py
+from shapely.geometry.proxy import CachingGeometryProxy
 
 
 def geos_multilinestring_from_py(ob):
@@ -115,18 +116,18 @@ class MultiLineString(BaseGeometry):
         return GeometrySequence(self, LineString)
 
 
-class MultiLineStringAdapter(MultiLineString):
+class MultiLineStringAdapter(CachingGeometryProxy, MultiLineString):
 
     """Adapts sequences of sequences or numpy arrays to the multilinestring
     interface.
     """
     
     context = None
-    #__geom = None
     _owned = False
 
     def __init__(self, context):
         self.context = context
+        self.factory = geos_multilinestring_from_py
 
     @property
     def _ndim(self):
@@ -140,14 +141,6 @@ class MultiLineStringAdapter(MultiLineString):
             # Fall back on list
             return len(self.context[0][0])
 
-    @property
-    def _geom(self):
-        """Keeps the GEOS geometry in synch with the context."""
-        if self.__geom__ is not None:
-            lgeos.GEOSGeom_destroy(self.__geom__)
-        self.__geom__ = geos_multilinestring_from_py(self.context)[0]
-        return self.__geom__
-
 
 def asMultiLineString(context):
     """Factory for MultiLineStringAdapter instances."""
@@ -162,4 +155,3 @@ def _test():
 
 if __name__ == "__main__":
     _test()
-

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/geometry/multipoint.py

@@ -7,6 +7,7 @@ from ctypes import byref, c_double, c_int, c_void_p, cast, POINTER, pointer
 from shapely.geos import lgeos
 from shapely.geometry.base import BaseGeometry, GeometrySequence, exceptNull
 from shapely.geometry.point import Point, geos_point_from_py
+from shapely.geometry.proxy import CachingGeometryProxy
 
 
 def geos_multipoint_from_py(ob):
@@ -137,18 +138,18 @@ class MultiPoint(BaseGeometry):
         return GeometrySequence(self, Point)
         
 
-class MultiPointAdapter(MultiPoint):
+class MultiPointAdapter(CachingGeometryProxy, MultiPoint):
 
     """Adapts a Python coordinate pair or a numpy array to the multipoint
     interface.
     """
     
     context = None
-    #__geom__ = None
     _owned = False
 
     def __init__(self, context):
         self.context = context
+        self.factory = geos_multipoint_from_py
 
     @property
     def _ndim(self):
@@ -162,14 +163,6 @@ class MultiPointAdapter(MultiPoint):
             # Fall back on list
             return len(self.context[0])
 
-    @property
-    def _geom(self):
-        """Keeps the GEOS geometry in synch with the context."""
-        if self.__geom__ is not None:
-            lgeos.GEOSGeom_destroy(self.__geom__)
-        self.__geom__, n = geos_multipoint_from_py(self.context)
-        return self.__geom__
-
     @property
     def __array_interface__(self):
         """Provide the Numpy array protocol."""
@@ -192,4 +185,3 @@ def _test():
 
 if __name__ == "__main__":
     _test()
-

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/geometry/multipolygon.py

@@ -7,6 +7,7 @@ from ctypes import byref, c_double, c_int, c_void_p, cast, POINTER, pointer
 from shapely.geos import lgeos
 from shapely.geometry.base import BaseGeometry, GeometrySequence, exceptNull
 from shapely.geometry.polygon import Polygon, geos_polygon_from_py
+from shapely.geometry.proxy import CachingGeometryProxy
 
 
 def geos_multipolygon_from_py(ob):
@@ -101,18 +102,18 @@ class MultiPolygon(BaseGeometry):
         return GeometrySequence(self, Polygon)
 
 
-class MultiPolygonAdapter(MultiPolygon):
+class MultiPolygonAdapter(CachingGeometryProxy, MultiPolygon):
 
     """Adapts sequences of sequences or numpy arrays to the multipolygon
     interface.
     """
     
     context = None
-    #__geom = None
     _owned = False
 
     def __init__(self, context):
         self.context = context
+        self.factory = geos_multipolygon_from_py
 
     @property
     def _ndim(self):
@@ -126,14 +127,6 @@ class MultiPolygonAdapter(MultiPolygon):
             # Fall back on list
             return len(self.context[0][0][0])
 
-    @property
-    def _geom(self):
-        """Keeps the GEOS geometry in synch with the context."""
-        if self.__geom__ is not None:
-            lgeos.GEOSGeom_destroy(self.__geom__)
-        self.__geom__ = geos_multipolygon_from_py(self.context)[0]
-        return self.__geom__
-
 
 def asMultiPolygon(context):
     """Factory for MultiLineStringAdapter instances."""
@@ -148,4 +141,3 @@ def _test():
 
 if __name__ == "__main__":
     _test()
-

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/geometry/point.py

@@ -10,6 +10,7 @@ from ctypes import cast, POINTER
 from shapely.geos import lgeos, DimensionError
 from shapely.geometry.base import BaseGeometry, CoordinateSequence
 from shapely.geometry.base import exceptNull
+from shapely.geometry.proxy import CachingGeometryProxy
 
 
 def geos_point_from_py(ob, update_geom=None, update_ndim=0):
@@ -199,17 +200,16 @@ class Point(BaseGeometry):
     coords = property(BaseGeometry._get_coords, _set_coords)
 
 
-class PointAdapter(Point):
+class PointAdapter(CachingGeometryProxy, Point):
 
     """Adapts a Python coordinate pair or a numpy array to the point interface.
     """
 
-    context = None
-    __geom = None
     _owned = False
 
     def __init__(self, context):
         self.context = context
+        self.factory = geos_point_from_py
 
     @property
     def _ndim(self):
@@ -223,14 +223,6 @@ class PointAdapter(Point):
             # Fall back on list
             return len(self.context)
 
-    @property
-    def _geom(self):
-        """Keeps the GEOS geometry in synch with the context."""
-        if self.__geom is not None:
-            lgeos.GEOSGeom_destroy(self.__geom)
-        self.__geom, n = geos_point_from_py(self.context)
-        return self.__geom
-
     # TODO: reimplement x, y, z properties without calling invoking _geom
 
     @property
@@ -262,4 +254,3 @@ def _test():
 
 if __name__ == "__main__":
     _test()
-

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/geometry/polygon.py

@@ -7,6 +7,7 @@ from ctypes import byref, c_double, c_int, c_void_p, cast, POINTER, pointer
 from shapely.geos import lgeos
 from shapely.geometry.base import BaseGeometry, exceptNull
 from shapely.geometry.linestring import LineString, LineStringAdapter
+from shapely.geometry.proxy import PolygonProxy
 
 
 def geos_linearring_from_py(ob, update_geom=None, update_ndim=0):
@@ -16,7 +17,8 @@ def geos_linearring_from_py(ob, update_geom=None, update_ndim=0):
         assert len(array['shape']) == 2
         m = array['shape'][0]
         n = array['shape'][1]
-        assert m >= 2
+        if m < 3:
+            raise ValueError, "A LinearRing must have at least 3 coordinate tuples"
         assert n == 2 or n == 3
 
         # Make pointer to the coordinate array
@@ -75,8 +77,9 @@ def geos_linearring_from_py(ob, update_geom=None, update_ndim=0):
         # Fall back on list
         m = len(ob)
         n = len(ob[0])
-        assert m >= 2
-        assert n == 2 or n == 3
+        if m < 3:
+            raise ValueError, "A LinearRing must have at least 3 coordinate tuples"
+        assert (n == 2 or n == 3)
 
         # Add closing coordinates if not provided
         if ob[0][0] != ob[-1][0] or ob[0][1] != ob[-1][1]:
@@ -282,10 +285,7 @@ def geos_polygon_from_py(shell, holes=None):
         if holes:
             ob = holes
             L = len(ob)
-            try:
-                N = len(ob[0][0])
-            except:
-                import pdb; pdb.set_trace()
+            N = len(ob[0][0])
             assert L >= 1
             assert N == 2 or N == 3
 
@@ -395,7 +395,7 @@ class Polygon(BaseGeometry):
             }
 
 
-class PolygonAdapter(Polygon):
+class PolygonAdapter(PolygonProxy, Polygon):
 
     """Adapts sequences of sequences or numpy arrays to the polygon
     interface.
@@ -408,6 +408,8 @@ class PolygonAdapter(Polygon):
     def __init__(self, shell, holes=None):
         self.shell = shell
         self.holes = holes
+        self.context = (shell, holes)
+        self.factory = geos_polygon_from_py
 
     @property
     def _ndim(self):
@@ -421,14 +423,6 @@ class PolygonAdapter(Polygon):
             # Fall back on list
             return len(self.shell[0])
 
-    @property
-    def _geom(self):
-        """Keeps the GEOS geometry in synch with the context."""
-        if self.__geom__ is not None:
-            lgeos.GEOSGeom_destroy(self.__geom__)
-        self.__geom__ = geos_polygon_from_py(self.shell, self.holes)[0]  
-        return self.__geom__
-
 
 def asPolygon(shell, holes=None):
     """Factory for PolygonAdapter instances."""
@@ -442,4 +436,3 @@ def _test():
 
 if __name__ == "__main__":
     _test()
-

Shapely-1.0.5/shapely/geometry/proxy.py

@@ -0,0 +1,48 @@
+"""
+Proxy for coordinates stored outside Shapely geometries.
+"""
+
+from shapely.geos import lgeos
+
+
+class CachingGeometryProxy(object):
+
+    context = None
+    factory = None
+    __geom__ = None
+    _gtag = None
+
+    def __init__(self, context):
+        self.context = context
+
+    @property
+    def _geom(self):
+        """Keeps the GEOS geometry in synch with the context."""
+        gtag = self.gtag()
+        if gtag != self._gtag:
+            if self.__geom__ is not None:
+                lgeos.GEOSGeom_destroy(self.__geom__)
+            self.__geom__, n = self.factory(self.context)
+        elif self.__geom__ is None:
+            self.__geom__, n = self.factory(self.context)
+        self._gtag = gtag
+        return self.__geom__
+        
+    def gtag(self):
+        return hash(repr(self.context))
+
+
+class PolygonProxy(CachingGeometryProxy):
+
+    @property
+    def _geom(self):
+        """Keeps the GEOS geometry in synch with the context."""
+        gtag = self.gtag()
+        if gtag != self._gtag:
+            if self.__geom__ is not None:
+                lgeos.GEOSGeom_destroy(self.__geom__)
+            self.__geom__, n = self.factory(self.context[0], self.context[1])
+        elif self.__geom__ is None:
+            self.__geom__, n = self.factory(self.context[0], self.context[1])
+        self._gtag = gtag
+        return self.__geom__

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/iterops.py

@@ -1,4 +1,5 @@
 """
+Iterative forms of operations.
 """
 
 from ctypes import c_char_p, c_size_t
@@ -36,7 +37,7 @@ class BinaryPredicateIterator(object):
             except TypeError:
                 this_geom = item
                 ob = this_geom
-            if this_geom._geom is None:
+            if not this_geom._geom:
                 raise ValueError, "Null geometry supports no operations"
             retval = self.fn(geom._geom, this_geom._geom)
             if retval == 2:
@@ -53,4 +54,3 @@ within = BinaryPredicateIterator(lgeos.GEOSWithin)
 contains = BinaryPredicateIterator(lgeos.GEOSContains)
 overlaps = BinaryPredicateIterator(lgeos.GEOSOverlaps)
 equals = BinaryPredicateIterator(lgeos.GEOSEquals)
-

Shapely-1.0.5/shapely/ops.py

@@ -0,0 +1,29 @@
+from shapely.geos import lgeos
+from shapely.geometry.base import geom_factory, BaseGeometry
+from shapely.geometry import asShape, asLineString
+from ctypes import byref, c_void_p
+
+def shapeup(ob):
+    if isinstance(ob, BaseGeometry):
+        return ob
+    else:
+        try:
+            return asShape(ob)
+        except ValueError:
+            return asLineString(ob)
+
+def polygonize(iterator):
+    """Creates polygons from a list of LineString objects.
+    """
+    lines = [shapeup(ob) for ob in iterator]
+    geom_array_type = c_void_p * len(lines)
+    geom_array = geom_array_type()
+    for i, line in enumerate(lines):
+        geom_array[i] = line._geom
+    product = lgeos.GEOSPolygonize(byref(geom_array), len(lines))
+    collection = geom_factory(product)
+    for g in collection.geoms:
+        clone = lgeos.GEOSGeom_clone(g._geom)
+        g = geom_factory(clone)
+        g._owned = False
+        yield g

{Shapely-1.0.3 → Shapely-1.0.5}/shapely/predicates.py

@@ -27,7 +27,7 @@ class BinaryPredicate(object):
         return self
 
     def __call__(self, other):
-        if self.context._geom is None or other._geom is None:
+        if not self.context._geom or not other._geom:
             raise ValueError, "Null geometry supports no operations"
         return bool(self.fn(self.context._geom, other._geom))
 
@@ -49,7 +49,7 @@ class UnaryPredicate(object):
         self.fn.errcheck = errcheck
 
     def __get__(self, obj, objtype=None):
-        if obj._geom is None:
+        if not obj._geom:
             raise ValueError, "Null geometry supports no operations"
         return bool(self.fn(obj._geom))
     

{Shapely-1.0.3 → Shapely-1.0.5}/tests/LineString.txt

@@ -97,8 +97,8 @@ Now, the inverse
 
     >>> la.__array_interface__ == la.context.__array_interface__
     True
-    >>> as = asarray(la)
-    >>> as
+    >>> pas = asarray(la)
+    >>> pas
     array([[ 1.,  2.],
            [ 3.,  4.]])
 

{Shapely-1.0.3 → Shapely-1.0.5}/tests/MultiPoint.txt

@@ -71,8 +71,8 @@ Adapter
   >>> geoma.__array_interface__ == geoma.context.__array_interface__
   True
   >>> from numpy import asarray
-  >>> as = asarray(geoma)
-  >>> as
+  >>> pas = asarray(geoma)
+  >>> pas
   array([[ 1.,  2.],
          [ 3.,  4.]])
 

{Shapely-1.0.3 → Shapely-1.0.5}/tests/Point.txt

@@ -99,8 +99,8 @@ Adapter
   >>> pa.__array_interface__ == pa.context.__array_interface__
   True
   >>> from numpy import asarray
-  >>> as = asarray(pa)
-  >>> as
+  >>> pas = asarray(pa)
+  >>> pas
   array([ 1.,  2.])
 
   Adapt a coordinate list to a point
@@ -123,8 +123,8 @@ Adapter
   Now, the inverse
   
   >>> ai = pa.__array_interface__
-  >>> as = asarray(pa)
-  >>> as
+  >>> pas = asarray(pa)
+  >>> pas
   array([ 1.,  4.])
 
 Test Non-operability of Null geometry

{Shapely-1.0.3 → Shapely-1.0.5}/tests/Polygon.txt

@@ -139,4 +139,10 @@ Check that we can set coordinates of a null geometry
   >>> print r_null.length # doctest: +ELLIPSIS
   3.414...
 
+Error handling
+--------------
 
+  >>> p = Polygon([[1,2], [2, 3]])
+  Traceback (most recent call last):
+  ...
+  ValueError: A LinearRing must have at least 3 coordinate tuples

{Shapely-1.0.3 → Shapely-1.0.5}/tests/Predicates.txt

@@ -28,13 +28,7 @@ Binary Predicates
   >>> point.equals(Point(-1.0, -1.0))
   False
 
-  >>> point == Point(0.0, 0.0)
-  True
-
-  >>> point == Point(-1.0, -1.0)
-  False
-
-  >>> point != Point(-1.0, -1.0)
+  >>> point.equals(Point(0.0, 0.0))
   True
 
 Unary Predicates

Shapely-1.0.5/tests/polygonize.txt

@@ -0,0 +1,29 @@
+Polygonizing
+============
+
+    >>> from shapely.geometry import LineString, Point
+    >>> from shapely.ops import polygonize
+    >>> lines = [
+    ...     LineString(((0, 0), (1, 1))),
+    ...     LineString(((0, 0), (0, 1))),
+    ...     LineString(((0, 1), (1, 1))),
+    ...     LineString(((1, 1), (1, 0))),
+    ...     LineString(((1, 0), (0, 0))),
+    ...     LineString(((5, 5), (6, 6))),
+    ...     Point(0, 0),
+    ...     ]
+    >>> result = polygonize(lines)
+    >>> list(result)
+    [<shapely.geometry.polygon.Polygon object at ...>, <shapely.geometry.polygon.Polygon object at ...>]
+
+    >>> lines2 = [
+    ...     ((0, 0), (1, 1)),
+    ...     ((0, 0), (0, 1)),
+    ...     ((0, 1), (1, 1)),
+    ...     ((1, 1), (1, 0)),
+    ...     ((1, 0), (0, 0)),
+    ...     ((5, 5), (6, 6)),
+    ...     ]
+    >>> result2 = polygonize(lines2)
+    >>> list(result2)
+    [<shapely.geometry.polygon.Polygon object at ...>, <shapely.geometry.polygon.Polygon object at ...>]

Shapely-1.0.3/CHANGES.txt

@@ -1,10 +0,0 @@
-All tickets are children of http://trac.gispython.org/projects/PCL/ticket/
-
-1.0.3 (2008-04-09)
-------------------
-- Do not release GIL when calling GEOS functions (#158).
-- Prevent faults when chaining multiple GEOS operators (#159).
-
-1.0.2 (2008-02-26)
-------------------
-- Fix loss of dimensionality in polygon rings (#155).

Shapely-1.0.3/shapely/geometry/collection.py

@@ -1,24 +0,0 @@
-"""
-Geometry collections.
-"""
-
-from shapely.geometry.base import BaseGeometry
-
-class GeometryCollection(BaseGeometry):
-
-    """A geometry collection.
-    """
-
-    def __init__(self):
-        BaseGeometry.__init__(self)
-
-
-# Test runner
-def _test():
-    import doctest
-    doctest.testmod()
-
-
-if __name__ == "__main__":
-    _test()
-

Shapely-1.0.3/shapely/geometry/geo.py

@@ -1,35 +0,0 @@
-
-from point import PointAdapter
-from linestring import LineStringAdapter
-from polygon import PolygonAdapter
-from multipoint import MultiPointAdapter
-from multilinestring import MultiLineStringAdapter
-from multipolygon import MultiPolygonAdapter
-
-
-def asShape(context):
-    """Adapts the context to a geometry interface. The coordinates remain
-    stored in the context.
-    """
-    if hasattr(context, "__geo_interface__"):
-        ob = context.__geo_interface__
-    else:
-        ob = context
-
-    geom_type = ob.get("type").lower()
-
-    if geom_type == "point":
-        return PointAdapter(ob["coordinates"])
-    elif geom_type == "linestring":
-        return LineStringAdapter(ob["coordinates"])
-    elif geom_type == "polygon":
-        return PolygonAdapter(ob["coordinates"][0], ob["coordinates"][1:])
-    elif geom_type == "multipoint":
-        return MultiPointAdapter(ob["coordinates"])
-    elif geom_type == "multilinestring":
-        return MultiLineStringAdapter(ob["coordinates"])
-    elif geom_type == "multipolygon":
-        return MultiPolygonAdapter(ob["coordinates"])
-    else:
-        raise ValueError, "Unknown geometry type: %s" % geom_type
-