shapely 2.0.6__tar.gz → 2.0.7__tar.gz

{shapely-2.0.6 → shapely-2.0.7}/CHANGES.txt

@@ -1,6 +1,18 @@
 Changes
 =======
 
+2.0.7 (2025-01-30)
+------------------
+
+Bug fixes:
+
+- Fix compilation error on certain Linux platforms, such as Alpine (#1945).
+- Fixes crash when reading nonlinear geometry types (CircularString,
+  CompoundCurve, MultiCurve, CurvePolygon, MultiSurface) from WKB/WKT with
+  GEOS >= 3.13; these types are not yet supported in Shapely and now raise a
+  ``NotImplementedError`` (#2160)
+- Fix the ``project`` method to return a Python float (#2093).
+
 2.0.6 (2024-08-19)
 ------------------
 

{shapely-2.0.6 → shapely-2.0.7}/CITATION.cff

@@ -2,8 +2,8 @@ cff-version: 1.2.0
 message: "Please cite this software using these metadata."
 type: software
 title: Shapely
-version: "2.0.6"
-date-released: "2024-08-19"
+version: "2.0.7"
+date-released: "2025-01-30"
 doi: 10.5281/zenodo.5597138
 abstract: "Manipulation and analysis of geometric objects in the Cartesian plane."
 repository-artifact: https://pypi.org/project/Shapely

{shapely-2.0.6 → shapely-2.0.7}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.2
 Name: shapely
-Version: 2.0.6
+Version: 2.0.7
 Summary: Manipulation and analysis of geometric objects
 Author: Sean Gillies
 Maintainer: Shapely contributors
@@ -22,6 +22,7 @@ Classifier: Programming Language :: Python :: 3.9
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
 Classifier: Topic :: Scientific/Engineering :: GIS
 Requires-Python: >=3.7
 Description-Content-Type: text/x-rst
@@ -53,12 +54,6 @@ Shapely
    :alt: Github Actions status
    :target: https://github.com/shapely/shapely/actions?query=branch%3Amain
 
-.. Travis CI status -- https://travis-ci.com
-
-.. image:: https://travis-ci.com/shapely/shapely.svg?branch=main
-   :alt: Travis CI status
-   :target: https://travis-ci.com/github/shapely/shapely
-
 .. PyPI
 
 .. image:: https://img.shields.io/pypi/v/shapely.svg

{shapely-2.0.6 → shapely-2.0.7}/README.rst

@@ -14,12 +14,6 @@ Shapely
    :alt: Github Actions status
    :target: https://github.com/shapely/shapely/actions?query=branch%3Amain
 
-.. Travis CI status -- https://travis-ci.com
-
-.. image:: https://travis-ci.com/shapely/shapely.svg?branch=main
-   :alt: Travis CI status
-   :target: https://travis-ci.com/github/shapely/shapely
-
 .. PyPI
 
 .. image:: https://img.shields.io/pypi/v/shapely.svg

{shapely-2.0.6 → shapely-2.0.7}/pyproject.toml

@@ -39,6 +39,7 @@ classifiers = [
     "Programming Language :: Python :: 3.10",
     "Programming Language :: Python :: 3.11",
     "Programming Language :: Python :: 3.12",
+    "Programming Language :: Python :: 3.13",
     "Topic :: Scientific/Engineering :: GIS",
 ]
 requires-python = ">=3.7"

{shapely-2.0.6 → shapely-2.0.7}/shapely/_version.py

@@ -8,11 +8,11 @@ import json
 
 version_json = '''
 {
- "date": "2024-08-19T23:29:00+0200",
+ "date": "2025-01-30T18:00:49-0700",
  "dirty": false,
  "error": null,
- "full-revisionid": "5a4207d2fb74b25d7bb2fb5d04813f68a3a612f4",
- "version": "2.0.6"
+ "full-revisionid": "341209eebf905ad10d3050b67ab495129963dae8",
+ "version": "2.0.7"
 }
 '''  # END VERSION_JSON
 

{shapely-2.0.6 → shapely-2.0.7}/shapely/geometry/base.py

@@ -830,7 +830,9 @@ class BaseGeometry(shapely.Geometry):
 
         Alias of `project`.
         """
-        return shapely.line_locate_point(self, other, normalized=normalized)
+        return _maybe_unpack(
+            shapely.line_locate_point(self, other, normalized=normalized)
+        )
 
     def project(self, other, normalized=False):
         """Returns the distance along this geometry to a point nearest the
@@ -841,7 +843,9 @@ class BaseGeometry(shapely.Geometry):
 
         Alias of `line_locate_point`.
         """
-        return shapely.line_locate_point(self, other, normalized=normalized)
+        return _maybe_unpack(
+            shapely.line_locate_point(self, other, normalized=normalized)
+        )
 
     def line_interpolate_point(self, distance, normalized=False):
         """Return a point at the specified distance along a linear geometry

{shapely-2.0.6 → shapely-2.0.7}/shapely/tests/geometry/test_geometry_base.py

@@ -164,6 +164,10 @@ def test_array_argument_binary_predicates(op):
     expected = np.array([getattr(polygon, op)(p) for p in points], dtype=bool)
     np.testing.assert_array_equal(result, expected)
 
+    # check scalar
+    result = getattr(polygon, op)(points[0])
+    assert type(result) is bool
+
 
 @pytest.mark.parametrize(
     "op, kwargs",
@@ -188,6 +192,10 @@ def test_array_argument_binary_predicates2(op, kwargs):
     expected = np.array([getattr(polygon, op)(p, **kwargs) for p in points], dtype=bool)
     np.testing.assert_array_equal(result, expected)
 
+    # check scalar
+    result = getattr(polygon, op)(points[0], **kwargs)
+    assert type(result) is bool
+
 
 @pytest.mark.parametrize(
     "op",
@@ -207,6 +215,10 @@ def test_array_argument_binary_geo(op):
     expected = np.array([getattr(box, op)(g) for g in polygons], dtype=object)
     assert_geometries_equal(result, expected)
 
+    # check scalar
+    result = getattr(box, op)(polygons[0])
+    assert isinstance(result, (Polygon, MultiPolygon))
+
 
 @pytest.mark.parametrize("op", ["distance", "hausdorff_distance"])
 def test_array_argument_float(op):
@@ -218,23 +230,42 @@ def test_array_argument_float(op):
     expected = np.array([getattr(polygon, op)(p) for p in points], dtype="float64")
     np.testing.assert_array_equal(result, expected)
 
+    # check scalar
+    result = getattr(polygon, op)(points[0])
+    assert type(result) is float
 
-def test_array_argument_linear():
+
+@pytest.mark.parametrize("op", ["line_interpolate_point", "interpolate"])
+def test_array_argument_linear_point(op):
     line = LineString([(0, 0), (0, 1), (1, 1)])
     distances = np.array([0, 0.5, 1])
-    result = line.line_interpolate_point(distances)
+
+    result = getattr(line, op)(distances)
     assert isinstance(result, np.ndarray)
     expected = np.array(
         [line.line_interpolate_point(d) for d in distances], dtype=object
     )
     assert_geometries_equal(result, expected)
 
+    # check scalar
+    result = getattr(line, op)(distances[0])
+    assert isinstance(result, Point)
+
+
+@pytest.mark.parametrize("op", ["line_locate_point", "project"])
+def test_array_argument_linear_float(op):
+    line = LineString([(0, 0), (0, 1), (1, 1)])
     points = shapely.points([(0, 0), (0.5, 0.5), (1, 1)])
-    result = line.line_locate_point(points)
+
+    result = getattr(line, op)(points)
     assert isinstance(result, np.ndarray)
     expected = np.array([line.line_locate_point(p) for p in points], dtype="float64")
     np.testing.assert_array_equal(result, expected)
 
+    # check scalar
+    result = getattr(line, op)(points[0])
+    assert type(result) is float
+
 
 def test_array_argument_buffer():
     point = Point(1, 1)
@@ -244,3 +275,7 @@ def test_array_argument_buffer():
     assert isinstance(result, np.ndarray)
     expected = np.array([point.buffer(d) for d in distances], dtype=object)
     assert_geometries_equal(result, expected)
+
+    # check scalar
+    result = point.buffer(distances[0])
+    assert isinstance(result, Polygon)

{shapely-2.0.6 → shapely-2.0.7}/shapely/tests/test_io.py

@@ -157,6 +157,31 @@ def test_from_wkt_empty(wkt):
     assert shapely.to_wkt(geom) == wkt
 
 
+# WKT from https://github.com/libgeos/geos/blob/main/tests/unit/io/WKBReaderTest.cpp
+@pytest.mark.parametrize(
+    "wkt",
+    (
+        "CIRCULARSTRING(1 3,2 4,3 1)",
+        "COMPOUNDCURVE(CIRCULARSTRING(1 3,2 4,3 1),(3 1,0 0))",
+        "CURVEPOLYGON(COMPOUNDCURVE(CIRCULARSTRING(0 0,2 0,2 1,2 3,4 3),(4 3,4 5,1 4,0 0)),CIRCULARSTRING(1.7 1,1.4 0.4,1.6 0.4,1.6 0.5,1.7 1))",  # noqa: E501
+        "MULTICURVE((0 0,5 5),COMPOUNDCURVE((-1 -1,0 0),CIRCULARSTRING(0 0,1 1,2 0)),CIRCULARSTRING(4 0,4 4,8 4))",  # noqa: E501
+        "MULTISURFACE(CURVEPOLYGON(CIRCULARSTRING(0 0,4 0,4 4,0 4,0 0),(1 1,3 3,3 1,1 1)),((10 10,14 12,11 10,10 10),(11 11,11.5 11,11 11.5,11 11)))",  # noqa: E501
+    ),
+)
+def test_from_wkt_nonlinear_unsupported(wkt):
+    if shapely.geos_version >= (3, 13, 0):
+        with pytest.raises(
+            NotImplementedError,
+            match="Nonlinear geometry types are not currently supported",
+        ):
+            shapely.from_wkt(wkt)
+
+    else:
+        # prior to GEOS 3.13 nonlinear types were rejected by GEOS on read from WKT
+        with pytest.raises(shapely.errors.GEOSException, match="Unknown type"):
+            shapely.from_wkt(wkt)
+
+
 def test_from_wkb():
     expected = shapely.points(1, 1)
     actual = shapely.from_wkb(POINT11_WKB)
@@ -255,6 +280,36 @@ def test_from_wkb_empty(geom):
     assert shapely.to_wkb(geom) == wkb
 
 
+# WKB from https://github.com/libgeos/geos/blob/main/tests/unit/io/WKBReaderTest.cpp
+@pytest.mark.parametrize(
+    "wkb",
+    (
+        # "CIRCULARSTRING(1 3,2 4,3 1)",
+        "010800000003000000000000000000F03F0000000000000840000000000000004000000000000010400000000000000840000000000000F03F",
+        # "COMPOUNDCURVE(CIRCULARSTRING(1 3,2 4,3 1),(3 1,0 0))",
+        "01090000200E16000002000000010800000003000000000000000000F03F0000000000000840000000000000004000000000000010400000000000000840000000000000F03F0102000000020000000000000000000840000000000000F03F00000000000000000000000000000000",
+        # "CURVEPOLYGON(COMPOUNDCURVE(CIRCULARSTRING(0 0,2 0,2 1,2 3,4 3),(4 3,4 5,1 4,0 0)),CIRCULARSTRING(1.7 1,1.4 0.4,1.6 0.4,1.6 0.5,1.7 1))",  # noqa: E501
+        "010A0000200E1600000200000001090000000200000001080000000500000000000000000000000000000000000000000000000000004000000000000000000000000000000040000000000000F03F00000000000000400000000000000840000000000000104000000000000008400102000000040000000000000000001040000000000000084000000000000010400000000000001440000000000000F03F000000000000104000000000000000000000000000000000010800000005000000333333333333FB3F000000000000F03F666666666666F63F9A9999999999D93F9A9999999999F93F9A9999999999D93F9A9999999999F93F000000000000E03F333333333333FB3F000000000000F03F",
+        # "MULTICURVE((0 0,5 5),COMPOUNDCURVE((-1 -1,0 0),CIRCULARSTRING(0 0,1 1,2 0)),CIRCULARSTRING(4 0,4 4,8 4))",  # noqa: E501
+        "010B000000030000000102000000020000000000000000000000000000000000000000000000000014400000000000001440010900000002000000010200000002000000000000000000F0BF000000000000F0BF0000000000000000000000000000000001080000000300000000000000000000000000000000000000000000000000F03F000000000000F03F00000000000000400000000000000000010800000003000000000000000000104000000000000000000000000000001040000000000000104000000000000020400000000000001040",
+        # "MULTISURFACE(CURVEPOLYGON(CIRCULARSTRING(0 0,4 0,4 4,0 4,0 0),(1 1,3 3,3 1,1 1)),((10 10,14 12,11 10,10 10),(11 11,11.5 11,11 11.5,11 11)))",  # noqa: E501
+        "010C00000002000000010A000000020000000108000000050000000000000000000000000000000000000000000000000010400000000000000000000000000000104000000000000010400000000000000000000000000000104000000000000000000000000000000000010200000004000000000000000000F03F000000000000F03F000000000000084000000000000008400000000000000840000000000000F03F000000000000F03F000000000000F03F01030000000200000004000000000000000000244000000000000024400000000000002C40000000000000284000000000000026400000000000002440000000000000244000000000000024400400000000000000000026400000000000002640000000000000274000000000000026400000000000002640000000000000274000000000000026400000000000002640",
+    ),
+)
+def test_from_wkb_nonlinear_unsupported(wkb):
+    if shapely.geos_version >= (3, 13, 0):
+        with pytest.raises(
+            NotImplementedError,
+            match="Nonlinear geometry types are not currently supported",
+        ):
+            shapely.from_wkb(wkb)
+
+    else:
+        # prior to GEOS 3.13 nonlinear types were rejected by GEOS on read from WKB
+        with pytest.raises(shapely.errors.GEOSException, match="Unknown WKB type"):
+            shapely.from_wkb(wkb)
+
+
 def test_to_wkt():
     point = shapely.points(1, 1)
     actual = shapely.to_wkt(point)

{shapely-2.0.6 → shapely-2.0.7}/shapely.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.2
 Name: shapely
-Version: 2.0.6
+Version: 2.0.7
 Summary: Manipulation and analysis of geometric objects
 Author: Sean Gillies
 Maintainer: Shapely contributors
@@ -22,6 +22,7 @@ Classifier: Programming Language :: Python :: 3.9
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
 Classifier: Topic :: Scientific/Engineering :: GIS
 Requires-Python: >=3.7
 Description-Content-Type: text/x-rst
@@ -53,12 +54,6 @@ Shapely
    :alt: Github Actions status
    :target: https://github.com/shapely/shapely/actions?query=branch%3Amain
 
-.. Travis CI status -- https://travis-ci.com
-
-.. image:: https://travis-ci.com/shapely/shapely.svg?branch=main
-   :alt: Travis CI status
-   :target: https://travis-ci.com/github/shapely/shapely
-
 .. PyPI
 
 .. image:: https://img.shields.io/pypi/v/shapely.svg

{shapely-2.0.6 → shapely-2.0.7}/src/pygeom.c

@@ -23,6 +23,16 @@ PyObject* GeometryObject_FromGEOS(GEOSGeometry* ptr, GEOSContextHandle_t ctx) {
   if (type_id == -1) {
     return NULL;
   }
+
+  // Nonlinear types (CircularString, CompoundCurve, MultiCurve, CurvePolygon,
+  // MultiSurface are not currently supported
+  // TODO: this can be removed once these types are added to the type registry
+  if (type_id >= 8) {
+    PyErr_Format(PyExc_NotImplementedError,
+                 "Nonlinear geometry types are not currently supported");
+    return NULL;
+  }
+
   PyObject* type_obj = PyList_GET_ITEM(geom_registry[0], type_id);
   if (type_obj == NULL) {
     return NULL;

{shapely-2.0.6 → shapely-2.0.7}/src/ufuncs.c

@@ -2178,7 +2178,7 @@ static void polygonize_func(char** args, const npy_intp* dimensions, const npy_i
 
   GEOS_INIT;
 
-  GEOSGeometry** geoms = malloc(sizeof(void*) * dimensions[1]);
+  const GEOSGeometry** geoms = malloc(sizeof(void*) * dimensions[1]);
   if (geoms == NULL) {
     errstate = PGERR_NO_MALLOC;
     goto finish;