shapely 1.8a1__tar.gz → 1.8a2__tar.gz

{Shapely-1.8a1 → Shapely-1.8a2}/CHANGES.txt

@@ -1,6 +1,37 @@
 Changes
 =======
 
+1.8a2 (TBD)
+-----------
+
+Python version support:
+
+Shapely 1.8 will support only Python versions >= 3.6.
+
+New features:
+
+- The STRtree nearest*() methods now take an optional argument that
+  specifies exclusion of the input geometry from results (#1115).
+- A GeometryTypeError has been added to shapely.errors and is consistently
+  raised instead of TypeError or ValueError as in version 1.7. For backwards
+  compatibility, the new exception will derive from TypeError and Value error
+  until version 2.0 (#1099).
+- The STRtree class constructor now takes an optional second argument, a
+  sequence of objects to be stored in the tree. If not provided, the sequence
+  indices of the geometries will be stored, as before (#1112).
+- The STRtree class has new query_geoms(), query_items(), nearest_geom(), and
+  nearest_item() methods (#1112). The query() and nearest() methods remain as
+  aliases for query_geoms() and nearest_geom().
+
+Bug fixes:
+
+- We no longer attempt to load libc to get the free function on Linux, but get
+  it from the global symbol table.
+- GEOS error messages printed when GEOS_getCoordSeq() is passed an empty
+  geometry are avoided by never passing an empty geometry (#1134).
+- Python's builtin super() is now used only as described in PEP 3135 (#1109).
+- Only load conda GEOS dll if it exists (on Windows) (#1108).
+
 1.8a1 (2021-03-03)
 ------------------
 
@@ -8,6 +39,10 @@ Shapely 1.8.0 will be a transitional version. There are a few bug fixes and new
 features, but it is mainly about warning of the upcoming changes in 2.0.0.
 Several more pre-releases before 1.8.0 are expected.
 
+Python version support:
+
+Shapely 1.8 will support only Python versions >= 3.5 (#884).
+
 Deprecations:
 
 The following functions and geometry attributes and methods will be removed in

{Shapely-1.8a1 → Shapely-1.8a2}/CREDITS.txt

@@ -24,6 +24,7 @@ Shapely is written by:
 * David Swinkels <davidswinkelss@gmail.com>
 * Denis Rykov <rykovd@gmail.com>
 * Erwin Sterrenburg <e.w.sterrenburg@gmail.com>
+* Felix Divo <4403130+felixdivo@users.noreply.github.com>
 * Felix Yan <felixonmars@archlinux.org>
 * Filipe Fernandes <ocefpaf@gmail.com>
 * Frédéric Junod <frederic.junod@camptocamp.com>

{Shapely-1.8a1/Shapely.egg-info → Shapely-1.8a2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: Shapely
-Version: 1.8a1
+Version: 1.8a2
 Summary: Geometric objects, predicates, and operations
 Home-page: https://github.com/Toblerity/Shapely
 Author: Sean Gillies
@@ -60,9 +60,9 @@ Description: =======
         Requirements
         ============
         
-        Shapely 1.7 requires
+        Shapely 1.8 requires
         
-        * Python 2.7, >=3.5
+        * Python >=3.6
         * GEOS >=3.3
         
         Installing Shapely
@@ -197,6 +197,7 @@ Description: =======
         * David Swinkels <davidswinkelss@gmail.com>
         * Denis Rykov <rykovd@gmail.com>
         * Erwin Sterrenburg <e.w.sterrenburg@gmail.com>
+        * Felix Divo <4403130+felixdivo@users.noreply.github.com>
         * Felix Yan <felixonmars@archlinux.org>
         * Filipe Fernandes <ocefpaf@gmail.com>
         * Frédéric Junod <frederic.junod@camptocamp.com>
@@ -295,6 +296,37 @@ Description: =======
         Changes
         =======
         
+        1.8a2 (TBD)
+        -----------
+        
+        Python version support:
+        
+        Shapely 1.8 will support only Python versions >= 3.6.
+        
+        New features:
+        
+        - The STRtree nearest*() methods now take an optional argument that
+          specifies exclusion of the input geometry from results (#1115).
+        - A GeometryTypeError has been added to shapely.errors and is consistently
+          raised instead of TypeError or ValueError as in version 1.7. For backwards
+          compatibility, the new exception will derive from TypeError and Value error
+          until version 2.0 (#1099).
+        - The STRtree class constructor now takes an optional second argument, a
+          sequence of objects to be stored in the tree. If not provided, the sequence
+          indices of the geometries will be stored, as before (#1112).
+        - The STRtree class has new query_geoms(), query_items(), nearest_geom(), and
+          nearest_item() methods (#1112). The query() and nearest() methods remain as
+          aliases for query_geoms() and nearest_geom().
+        
+        Bug fixes:
+        
+        - We no longer attempt to load libc to get the free function on Linux, but get
+          it from the global symbol table.
+        - GEOS error messages printed when GEOS_getCoordSeq() is passed an empty
+          geometry are avoided by never passing an empty geometry (#1134).
+        - Python's builtin super() is now used only as described in PEP 3135 (#1109).
+        - Only load conda GEOS dll if it exists (on Windows) (#1108).
+        
         1.8a1 (2021-03-03)
         ------------------
         
@@ -302,6 +334,10 @@ Description: =======
         features, but it is mainly about warning of the upcoming changes in 2.0.0.
         Several more pre-releases before 1.8.0 are expected.
         
+        Python version support:
+        
+        Shapely 1.8 will support only Python versions >= 3.5 (#884).
+        
         Deprecations:
         
         The following functions and geometry attributes and methods will be removed in
@@ -1068,13 +1104,12 @@ Classifier: Intended Audience :: Science/Research
 Classifier: License :: OSI Approved :: BSD License
 Classifier: Operating System :: OS Independent
 Classifier: Programming Language :: Python :: 3
-Classifier: Programming Language :: Python :: 3.5
 Classifier: Programming Language :: Python :: 3.6
 Classifier: Programming Language :: Python :: 3.7
 Classifier: Programming Language :: Python :: 3.8
 Classifier: Programming Language :: Python :: 3.9
 Classifier: Topic :: Scientific/Engineering :: GIS
-Requires-Python: >=3.5
+Requires-Python: >=3.6
+Provides-Extra: all
 Provides-Extra: vectorized
 Provides-Extra: test
-Provides-Extra: all

{Shapely-1.8a1 → Shapely-1.8a2}/README.rst

@@ -50,9 +50,9 @@ See the manual for more examples and guidance.
 Requirements
 ============
 
-Shapely 1.7 requires
+Shapely 1.8 requires
 
-* Python 2.7, >=3.5
+* Python >=3.6
 * GEOS >=3.3
 
 Installing Shapely

{Shapely-1.8a1 → Shapely-1.8a2/Shapely.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: Shapely
-Version: 1.8a1
+Version: 1.8a2
 Summary: Geometric objects, predicates, and operations
 Home-page: https://github.com/Toblerity/Shapely
 Author: Sean Gillies
@@ -60,9 +60,9 @@ Description: =======
         Requirements
         ============
         
-        Shapely 1.7 requires
+        Shapely 1.8 requires
         
-        * Python 2.7, >=3.5
+        * Python >=3.6
         * GEOS >=3.3
         
         Installing Shapely
@@ -197,6 +197,7 @@ Description: =======
         * David Swinkels <davidswinkelss@gmail.com>
         * Denis Rykov <rykovd@gmail.com>
         * Erwin Sterrenburg <e.w.sterrenburg@gmail.com>
+        * Felix Divo <4403130+felixdivo@users.noreply.github.com>
         * Felix Yan <felixonmars@archlinux.org>
         * Filipe Fernandes <ocefpaf@gmail.com>
         * Frédéric Junod <frederic.junod@camptocamp.com>
@@ -295,6 +296,37 @@ Description: =======
         Changes
         =======
         
+        1.8a2 (TBD)
+        -----------
+        
+        Python version support:
+        
+        Shapely 1.8 will support only Python versions >= 3.6.
+        
+        New features:
+        
+        - The STRtree nearest*() methods now take an optional argument that
+          specifies exclusion of the input geometry from results (#1115).
+        - A GeometryTypeError has been added to shapely.errors and is consistently
+          raised instead of TypeError or ValueError as in version 1.7. For backwards
+          compatibility, the new exception will derive from TypeError and Value error
+          until version 2.0 (#1099).
+        - The STRtree class constructor now takes an optional second argument, a
+          sequence of objects to be stored in the tree. If not provided, the sequence
+          indices of the geometries will be stored, as before (#1112).
+        - The STRtree class has new query_geoms(), query_items(), nearest_geom(), and
+          nearest_item() methods (#1112). The query() and nearest() methods remain as
+          aliases for query_geoms() and nearest_geom().
+        
+        Bug fixes:
+        
+        - We no longer attempt to load libc to get the free function on Linux, but get
+          it from the global symbol table.
+        - GEOS error messages printed when GEOS_getCoordSeq() is passed an empty
+          geometry are avoided by never passing an empty geometry (#1134).
+        - Python's builtin super() is now used only as described in PEP 3135 (#1109).
+        - Only load conda GEOS dll if it exists (on Windows) (#1108).
+        
         1.8a1 (2021-03-03)
         ------------------
         
@@ -302,6 +334,10 @@ Description: =======
         features, but it is mainly about warning of the upcoming changes in 2.0.0.
         Several more pre-releases before 1.8.0 are expected.
         
+        Python version support:
+        
+        Shapely 1.8 will support only Python versions >= 3.5 (#884).
+        
         Deprecations:
         
         The following functions and geometry attributes and methods will be removed in
@@ -1068,13 +1104,12 @@ Classifier: Intended Audience :: Science/Research
 Classifier: License :: OSI Approved :: BSD License
 Classifier: Operating System :: OS Independent
 Classifier: Programming Language :: Python :: 3
-Classifier: Programming Language :: Python :: 3.5
 Classifier: Programming Language :: Python :: 3.6
 Classifier: Programming Language :: Python :: 3.7
 Classifier: Programming Language :: Python :: 3.8
 Classifier: Programming Language :: Python :: 3.9
 Classifier: Topic :: Scientific/Engineering :: GIS
-Requires-Python: >=3.5
+Requires-Python: >=3.6
+Provides-Extra: all
 Provides-Extra: vectorized
 Provides-Extra: test
-Provides-Extra: all

{Shapely-1.8a1 → Shapely-1.8a2}/Shapely.egg-info/SOURCES.txt

@@ -127,7 +127,6 @@ tests/test_singularity.py
 tests/test_snap.py
 tests/test_split.py
 tests/test_strtree.py
-tests/test_strtree_nearest.py
 tests/test_styles.py
 tests/test_substring.py
 tests/test_svg.py

Shapely-1.8a2/VERSION.txt

@@ -0,0 +1 @@
+1.8a2

{Shapely-1.8a1 → Shapely-1.8a2}/_vendor/packaging/_structures.py

@@ -4,7 +4,7 @@
 from __future__ import absolute_import, division, print_function
 
 
-class Infinity(object):
+class Infinity:
 
     def __repr__(self):
         return "Infinity"
@@ -36,7 +36,7 @@ class Infinity(object):
 Infinity = Infinity()
 
 
-class NegativeInfinity(object):
+class NegativeInfinity:
 
     def __repr__(self):
         return "-Infinity"

{Shapely-1.8a1 → Shapely-1.8a2}/_vendor/packaging/markers.py

@@ -40,7 +40,7 @@ class UndefinedEnvironmentName(ValueError):
     """
 
 
-class Node(object):
+class Node:
 
     def __init__(self, value):
         self.value = value
@@ -254,7 +254,7 @@ def default_environment():
     }
 
 
-class Marker(object):
+class Marker:
 
     def __init__(self, marker):
         try:

{Shapely-1.8a1 → Shapely-1.8a2}/_vendor/packaging/requirements.py

@@ -72,7 +72,7 @@ NAMED_REQUIREMENT = \
 REQUIREMENT = stringStart + NAMED_REQUIREMENT + stringEnd
 
 
-class Requirement(object):
+class Requirement:
     """Parse a requirement.
 
     Parse a given requirement string into its parts, such as name, specifier,

{Shapely-1.8a1 → Shapely-1.8a2}/_vendor/packaging/specifiers.py

@@ -503,7 +503,7 @@ class Specifier(_IndividualSpecifier):
                 return False
 
         # Ensure that we do not allow a local version of the version mentioned
-        # in the specifier, which is techincally greater than, to match.
+        # in the specifier, which is technically greater than, to match.
         if prospective.local is not None:
             if Version(prospective.base_version) == Version(spec.base_version):
                 return False
@@ -589,7 +589,7 @@ def _pad_version(left, right):
 class SpecifierSet(BaseSpecifier):
 
     def __init__(self, specifiers="", prereleases=None):
-        # Split on , to break each indidivual specifier into it's own item, and
+        # Split on , to break each individual specifier into it's own item, and
         # strip each item to remove leading/trailing whitespace.
         specifiers = [s.strip() for s in specifiers.split(",") if s.strip()]
 

{Shapely-1.8a1 → Shapely-1.8a2}/_vendor/packaging/version.py

@@ -39,7 +39,7 @@ class InvalidVersion(ValueError):
     """
 
 
-class _BaseVersion(object):
+class _BaseVersion:
 
     def __hash__(self):
         return hash(self._key)

{Shapely-1.8a1 → Shapely-1.8a2}/docs/design.rst

@@ -3,7 +3,7 @@ Design Notes
 ============
 
 Shapely provides classes that implement, more or less, the interfaces in the
-OGC's simple feature acess specification [1]_. The classes are defined in
+OGC's simple feature access specification [1]_. The classes are defined in
 similarly named modules under ``shapely.geometry``: ``Point`` is in
 ``shapely.geometry.point``, ``MultiPolygon`` is in
 ``shapely.geometry.multipolygon``. These classes derive from

{Shapely-1.8a1 → Shapely-1.8a2}/docs/manual.rst

@@ -1290,7 +1290,7 @@ Shapely's binary :meth:`within` predicate.
 
 .. code-block:: python
 
-  class Within(object):
+  class Within:
       def __init__(self, o):
           self.o = o
       def __lt__(self, other):
@@ -1570,7 +1570,7 @@ Constructive Methods
 Shapely geometric object have several methods that yield new objects not
 derived from set-theoretic analysis.
 
-.. method:: object.buffer(distance, resolution=16, cap_style=1, join_style=1, mitre_limit=5.0)
+.. method:: object.buffer(distance, resolution=16, cap_style=1, join_style=1, mitre_limit=5.0, single_sided=False)
 
   Returns an approximate representation of all points within a given `distance`
   of the this geometric object.
@@ -1647,6 +1647,28 @@ With a `resolution` of 1, the buffer is a square patch.
   >>> q.area
   200.0
 
+You may want a buffer only on one side. You can achieve this effect with
+`single_sided` option.
+
+The side used is determined by the sign of the buffer distance:
+
+- a positive distance indicates the left-hand side
+- a negative distance indicates the right-hand side
+
+.. code-block:: pycon
+
+  >>> line = LineString([(0, 0), (1, 1), (0, 2), (2, 2), (3, 1), (1, 0)])
+  >>> left_hand_side = line.buffer(0.5, single_sided=True)
+  >>> right_hand_side = line.buffer(-0.3, single_sided=True)
+
+.. plot:: code/buffer_single_side.py
+
+Figure 10. Single sided buffer of 0.5 left hand (left) and of 0.3 right hand (right).
+
+The single-sided buffer of point geometries is the same as the regular buffer.
+The End Cap Style for single-sided buffers is always ignored, and forced to
+the equivalent of `CAP_STYLE.flat`.
+
 Passed a `distance` of 0, :meth:`buffer` can sometimes be used to "clean" self-touching
 or self-crossing polygons such as the classic "bowtie". Users have reported
 that very small distance values sometimes produce cleaner results than 0. Your
@@ -1690,7 +1712,7 @@ Buffering splits the polygon in two at the point where they touch.
 
 .. plot:: code/convex_hull.py
 
-Figure 10. Convex hull (blue) of 2 points (left) and of 6 points (right).
+Figure 11. Convex hull (blue) of 2 points (left) and of 6 points (right).
 
 .. attribute:: object.envelope
 
@@ -1722,7 +1744,7 @@ Figure 10. Convex hull (blue) of 2 points (left) and of 6 points (right).
 
 .. plot:: code/minimum_rotated_rectangle.py
 
-Figure 11. Minimum rotated rectangle for a multipoint feature (left) and a
+Figure 12. Minimum rotated rectangle for a multipoint feature (left) and a
 linestring feature (right).
 
 .. method:: object.parallel_offset(distance, side, resolution=16, join_style=1, mitre_limit=5.0)
@@ -1762,7 +1784,7 @@ linestring feature (right).
 
 .. plot:: code/parallel_offset.py
 
-Figure 12. Three styles of parallel offset lines on the left side of a simple
+Figure 13. Three styles of parallel offset lines on the left side of a simple
 line string (its starting point shown as a circle) and one offset on the right
 side, a multipart.
 
@@ -1770,7 +1792,7 @@ The effect of the `mitre_limit` parameter is shown below.
 
 .. plot:: code/parallel_offset_mitre.py
 
-Figure 13. Large and small mitre_limit values for left and right offsets.
+Figure 14. Large and small mitre_limit values for left and right offsets.
 
 .. method:: object.simplify(tolerance, preserve_topology=True)
 
@@ -1797,7 +1819,7 @@ Douglas-Peucker algorithm [6]_ is used.
 
 .. plot:: code/simplify.py
 
-Figure 14. Simplification of a nearly circular polygon using a tolerance of 0.2
+Figure 15. Simplification of a nearly circular polygon using a tolerance of 0.2
 (left) and 0.5 (right).
 
 .. note::
@@ -1925,7 +1947,7 @@ preserved or supported by 3D affine transformations.
 
   .. plot:: code/rotate.py
 
-  Figure 15. Rotation of a `LineString` (gray) by an angle of 90°
+  Figure 16. Rotation of a `LineString` (gray) by an angle of 90°
   counter-clockwise (blue) using different origins.
 
 .. function:: shapely.affinity.scale(geom, xfact=1.0, yfact=1.0, zfact=1.0, origin='center')
@@ -1967,7 +1989,7 @@ preserved or supported by 3D affine transformations.
 
   .. plot:: code/scale.py
 
-  Figure 16. Scaling of a gray triangle to blue result: a) by a factor of 1.5
+  Figure 17. Scaling of a gray triangle to blue result: a) by a factor of 1.5
   along x-direction, with reflection across y-axis; b) by a factor of 2 along
   x-direction with custom origin at (1, 1).
 
@@ -1999,7 +2021,7 @@ preserved or supported by 3D affine transformations.
 
   .. plot:: code/skew.py
 
-  Figure 17. Skewing of a gray "R" to blue result: a) by a shear angle of 20°
+  Figure 18. Skewing of a gray "R" to blue result: a) by a shear angle of 20°
   along the x-direction and an origin at (1, 1); b) by a shear angle of 30°
   along the y-direction, using default origin.
 
@@ -2186,6 +2208,33 @@ using functions in the :mod:`shapely.ops` module.
      <shapely.geometry.linestring.LineString object at 0x...>,
      <shapely.geometry.linestring.LineString object at 0x...>]
 
+Efficient Rectangle Clipping
+----------------------------
+The :func:`~shapely.ops.clip_by_rect` function in `shapely.ops` returns the
+portion of a geometry within a rectangle.
+
+.. function:: shapely.ops.clip_by_rect(geom, xmin, ymin, xmax, ymax)
+
+    The geometry is clipped in a fast but possibly dirty way. The output is
+    not guaranteed to be valid. No exceptions will be raised for topological
+    errors.
+
+    `New in version 1.7.`
+
+    Requires GEOS 3.5.0 or higher
+
+.. code-block:: python
+
+  >>> from shapely.geometry import Polygon
+  >>> from shapely.ops import clip_by_rect
+  >>> polygon = Polygon(
+         shell=[(0, 0), (0, 30), (30, 30), (30, 0), (0, 0)],
+         holes=[[(10, 10), (20, 10), (20, 20), (10, 20), (10, 10)]],
+      )
+  >>> clipped_polygon = clip_by_rect(polygon, 5, 5, 15, 15)
+  >>> print(clipped_polygon.wkt)
+  POLYGON ((5 5, 5 15, 10 15, 10 10, 15 10, 15 5, 5 5))
+
 Efficient Unions
 ----------------
 
@@ -2835,7 +2884,7 @@ Or a simple placemark-type object:
 
 .. code-block:: pycon
 
-  >>> class GeoThing(object):
+  >>> class GeoThing:
   ...     def __init__(self, d):
   ...         self.__geo_interface__ = d
   >>> thing = GeoThing({"type": "Point", "coordinates": (0.0, 0.0)})

{Shapely-1.8a1 → Shapely-1.8a2}/setup.cfg

@@ -2,9 +2,6 @@
 testpaths = tests
 addopts = --ignore _vendor --ignore shapely/examples
 
-[bdist_wheel]
-universal = 1
-
 [egg_info]
 tag_build = 
 tag_date = 0

{Shapely-1.8a1 → Shapely-1.8a2}/setup.py

@@ -82,7 +82,7 @@ if 'all' in sys.warnoptions:
     log.level = logging.DEBUG
 
 
-class GEOSConfig(object):
+class GEOSConfig:
     """Interface to config options from the `geos-config` utility
     """
 
@@ -197,7 +197,6 @@ setup_args = dict(
         'License :: OSI Approved :: BSD License',
         'Operating System :: OS Independent',
         'Programming Language :: Python :: 3',
-        'Programming Language :: Python :: 3.5',
         'Programming Language :: Python :: 3.6',
         'Programming Language :: Python :: 3.7',
         'Programming Language :: Python :: 3.8',
@@ -205,7 +204,7 @@ setup_args = dict(
         'Topic :: Scientific/Engineering :: GIS',
     ],
     cmdclass           = {},
-    python_requires    = '>=3.5',
+    python_requires    = '>=3.6',
     extras_require     = extra_reqs,
     package_data={
         'shapely': ['shapely/_geos.pxi']},

Shapely-1.8a2/shapely/__init__.py

@@ -0,0 +1 @@
+__version__ = "1.8a2"

{Shapely-1.8a1 → Shapely-1.8a2}/shapely/affinity.py

@@ -1,6 +1,7 @@
 """Affine transforms, both in general and specific, named transforms."""
 
 from math import sin, cos, tan, pi
+from shapely.errors import GeometryTypeError
 
 __all__ = ['affine_transform', 'rotate', 'scale', 'skew', 'translate']
 
@@ -92,7 +93,7 @@ def affine_transform(geom, matrix):
         return type(geom)([affine_transform(part, matrix)
                            for part in geom.geoms])
     else:
-        raise ValueError('Type %r not recognized' % geom.type)
+        raise GeometryTypeError('Type %r not recognized' % geom.type)
 
 
 def interpret_origin(geom, origin, ndim):

{Shapely-1.8a1 → Shapely-1.8a2}/shapely/algorithms/polylabel.py

@@ -3,7 +3,7 @@ from ..geos import TopologicalError
 from heapq import heappush, heappop
 
 
-class Cell(object):
+class Cell:
     """A `Cell`'s centroid property is a potential solution to finding the pole
     of inaccessibility for a given polygon. Rich comparison operators are used
     for sorting `Cell` objects in a priority queue based on the potential

{Shapely-1.8a1 → Shapely-1.8a2}/shapely/coords.py

@@ -11,7 +11,7 @@ from shapely.geos import lgeos
 from shapely.topology import Validating
 
 
-class CoordinateSequence(object):
+class CoordinateSequence:
     """
     Iterative access to coordinate tuples from the parent geometry's coordinate
     sequence.
@@ -41,15 +41,22 @@ class CoordinateSequence(object):
         self.__p__ = parent
 
     def _update(self):
-        self._ndim = self.__p__._ndim
-        self._cseq = lgeos.GEOSGeom_getCoordSeq(self.__p__._geom)
+        if self.__p__.is_empty:
+            self._ndim = None
+            self._cseq = None
+        else:
+            self._ndim = self.__p__._ndim
+            self._cseq = lgeos.GEOSGeom_getCoordSeq(self.__p__._geom)
 
     def __len__(self):
         self._update()
-        cs_len = c_uint(0)
-        if self._cseq:
-            lgeos.GEOSCoordSeq_getSize(self._cseq, byref(cs_len))
-        return cs_len.value
+        if self._cseq is None:
+            return 0
+        else:
+            cs_len = c_uint(0)
+            if self._cseq:
+                lgeos.GEOSCoordSeq_getSize(self._cseq, byref(cs_len))
+            return cs_len.value
 
     def __iter__(self):
         self._update()
@@ -74,7 +81,7 @@ class CoordinateSequence(object):
         m = self.__len__()
         has_z = self._ndim == 3
         if isinstance(key, int):
-            if key + m < 0 or key >= m:
+            if self._cseq is None or key + m < 0 or key >= m:
                 raise IndexError("index out of range")
             if key < 0:
                 i = m + key
@@ -88,17 +95,20 @@ class CoordinateSequence(object):
             else:
                 return (dx.value, dy.value)
         elif isinstance(key, slice):
-            res = []
-            start, stop, stride = key.indices(m)
-            for i in range(start, stop, stride):
-                lgeos.GEOSCoordSeq_getX(self._cseq, i, byref(dx))
-                lgeos.GEOSCoordSeq_getY(self._cseq, i, byref(dy))
-                if has_z:
-                    lgeos.GEOSCoordSeq_getZ(self._cseq, i, byref(dz))
-                    res.append((dx.value, dy.value, dz.value))
-                else:
-                    res.append((dx.value, dy.value))
-            return res
+            if self._cseq is None:
+                return []
+            else:
+                res = []
+                start, stop, stride = key.indices(m)
+                for i in range(start, stop, stride):
+                    lgeos.GEOSCoordSeq_getX(self._cseq, i, byref(dx))
+                    lgeos.GEOSCoordSeq_getY(self._cseq, i, byref(dy))
+                    if has_z:
+                        lgeos.GEOSCoordSeq_getZ(self._cseq, i, byref(dz))
+                        res.append((dx.value, dy.value, dz.value))
+                    else:
+                        res.append((dx.value, dy.value))
+                return res
         else:
             raise TypeError("key must be an index or slice")