pycontrails 0.57.0__tar.gz → 0.58.0__tar.gz

{pycontrails-0.57.0 → pycontrails-0.58.0}/.github/workflows/release.yaml

@@ -42,11 +42,11 @@ jobs:
       - name: Build wheels
         uses: pypa/cibuildwheel@v3.2
         env:
-          CIBW_BUILD: cp310-* cp311-* cp312-* cp313-*
+          CIBW_BUILD: cp311-* cp312-* cp313-* cp314-*
           CIBW_SKIP: '*-win32 *-manylinux_i686 *-musllinux*'
           CIBW_BUILD_VERBOSITY: 3
           CIBW_ARCHS_MACOS: x86_64 arm64
-          CIBW_TEST_SKIP: '*-macosx_arm64'
+          CIBW_TEST_SKIP: '*-macosx_arm64 cp314-*'
           # Completely isolate tests to prevent cibuildwheel from importing the
           # source instead of the wheel. This happens when tests/__init__.py is read.
           CIBW_TEST_EXTRAS: "complete,dev"

{pycontrails-0.57.0 → pycontrails-0.58.0}/.github/workflows/test.yaml

@@ -38,7 +38,7 @@ jobs:
       fail-fast: false
       matrix:
         os: [ubuntu-latest, windows-latest]
-        pyversion: ['3.10', '3.11', '3.12', '3.13']
+        pyversion: ['3.11', '3.12', '3.13']
     runs-on: ${{ matrix.os }}
 
     steps:

{pycontrails-0.57.0 → pycontrails-0.58.0}/CHANGELOG.md

@@ -1,5 +1,21 @@
 # Changelog
 
+## 0.58.0
+
+### Features
+
+- Build wheels for python 3.14. These are not yet tested in the CI (not all dependencies support python 3.14 yet).
+
+### Breaking changes
+
+- Drop support for python 3.10 per [NEP 29](https://numpy.org/neps/nep-0029-deprecation_policy.html).
+- Require scipy >= 1.12 to better support the `PycontrailsRegularGridInterpolator` interface.
+
+### Fixes
+
+- When instantiated with `copy=True`, `MetDataset` and `MetDataArray` now remove duplicate dimension values, preserving only the first occurrence of each. Previously, duplicates were silently retained (such duplicates are not expected with typical gridded weather data). With `copy=False`, a `ValueError` is raised if duplicate dimension values are present.
+- Update `PycontrailsRegularGridInterpolator` for improved compatibility with `scipy.interpolate.RegularGridInterpolator`. All methods now delegate to SciPy’s implementation, except for `method="linear"`, which retains a fast-path optimized in `PycontrailsRegularGridInterpolator`.
+
 ## 0.57.0
 
 ### Features

{pycontrails-0.57.0/pycontrails.egg-info → pycontrails-0.58.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pycontrails
-Version: 0.57.0
+Version: 0.58.0
 Summary: Python library for modeling aviation climate impacts
 Author-email: "Contrails.org" <py@contrails.org>
 License-Expression: Apache-2.0
@@ -13,23 +13,23 @@ Classifier: Development Status :: 4 - Beta
 Classifier: Intended Audience :: Science/Research
 Classifier: Operating System :: OS Independent
 Classifier: Programming Language :: Python :: 3
-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: Programming Language :: Python :: 3.14
 Classifier: Programming Language :: Python :: 3 :: Only
 Classifier: Topic :: Software Development :: Libraries :: Python Modules
 Classifier: Topic :: Scientific/Engineering :: Atmospheric Science
 Classifier: Topic :: Scientific/Engineering :: GIS
 Classifier: Typing :: Typed
-Requires-Python: >=3.10
+Requires-Python: >=3.11
 Description-Content-Type: text/markdown
 License-File: LICENSE
 License-File: NOTICE
 Requires-Dist: dask>=2022.3
 Requires-Dist: numpy>=1.22
 Requires-Dist: pandas>=2.0
-Requires-Dist: scipy>=1.10
+Requires-Dist: scipy>=1.12
 Requires-Dist: typing-extensions>=4.5; python_version < "3.12"
 Requires-Dist: xarray>=2022.3
 Provides-Extra: complete
@@ -140,7 +140,7 @@ Documentation and examples available at [py.contrails.org](https://py.contrails.
 
 ### Install with pip
 
-You can install pycontrails from PyPI with `pip` (Python 3.10 or later required):
+You can install pycontrails from PyPI with `pip` (Python 3.11 or later required):
 
 ```bash
 $ pip install pycontrails

{pycontrails-0.57.0 → pycontrails-0.58.0}/README.md

@@ -25,7 +25,7 @@ Documentation and examples available at [py.contrails.org](https://py.contrails.
 
 ### Install with pip
 
-You can install pycontrails from PyPI with `pip` (Python 3.10 or later required):
+You can install pycontrails from PyPI with `pip` (Python 3.11 or later required):
 
 ```bash
 $ pip install pycontrails

{pycontrails-0.57.0 → pycontrails-0.58.0}/docs/api.rst

@@ -83,13 +83,15 @@ ARCO ERA5
     datalib.ecmwf.arco_era5
 
 
-GOES
-""""
+Geostationary Satellites
+""""""""""""""""""""""""
 
 .. autosummary::
     :toctree: api/
 
     datalib.goes
+    datalib.himawari
+    datalib.geo_utils
 
 
 Low Earth Orbit Satellites

{pycontrails-0.57.0 → pycontrails-0.58.0}/docs/install.rst

@@ -17,7 +17,7 @@ The conda-forge package includes all optional runtime dependencies.
 pip install
 -----------
 
-With Python 3.10 or later, install the latest release from PyPI using ``pip``:
+With Python 3.11 or later, install the latest release from PyPI using ``pip``:
 
 .. code-block:: bash
 
@@ -27,7 +27,8 @@ With Python 3.10 or later, install the latest release from PyPI using ``pip``:
     # install with all optional dependencies
     $ pip install "pycontrails[complete]"
 
-Wheels are currently built and tested for python 3.10 - 3.13 on Linux, macOS, and Windows.
+Wheels are currently built for python 3.11 - 3.14 and tested for
+python 3.11 - 3.13 on Linux, macOS, and Windows.
 
 Install the latest development version directly from GitHub:
 

{pycontrails-0.57.0 → pycontrails-0.58.0}/docs/notebooks/GOES.ipynb

@@ -30,7 +30,7 @@
     "\n",
     "By default, any GOES data will be cached on disk. Set `cachestore=None` to disable caching when defining the `GOES` object.\n",
     "\n",
-    "We download data for channels 1, 2, and 3. These are needed for creating a true color image."
+    "We download data for bands 1, 2, and 3. These are needed for creating a true color image."
    ]
   },
   {
@@ -692,7 +692,7 @@
     }
    ],
    "source": [
-    "handler = goes.GOES(region=\"conus\", channels=(\"C01\", \"C02\", \"C03\"))\n",
+    "handler = goes.GOES(region=\"conus\", bands=(\"C01\", \"C02\", \"C03\"))\n",
     "\n",
     "# Download the data\n",
     "da = handler.get(\"2023-02-05T18:00:00\")\n",
@@ -3849,7 +3849,7 @@
     "\n",
     "The ash color scheme was originally developed to visualize volcanic ash. It is also useful for visualizing contrails.\n",
     "\n",
-    "We download data for channels 11, 14, and 15 to create an ash color image."
+    "We download data for bands 11, 14, and 15 to create an ash color image."
    ]
   },
   {
@@ -3858,7 +3858,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "handler = goes.GOES(region=\"conus\", channels=(\"C11\", \"C14\", \"C15\"))\n",
+    "handler = goes.GOES(region=\"conus\", bands=(\"C11\", \"C14\", \"C15\"))\n",
     "da = handler.get(\"2023-02-09T18:00:00\")\n",
     "\n",
     "rgb, src_crs, src_extent = goes.extract_visualization(da, color_scheme=\"ash\")"

{pycontrails-0.57.0 → pycontrails-0.58.0}/pycontrails/_version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '0.57.0'
-__version_tuple__ = version_tuple = (0, 57, 0)
+__version__ = version = '0.58.0'
+__version_tuple__ = version_tuple = (0, 58, 0)
 
-__commit_id__ = commit_id = 'g7b8b60b87'
+__commit_id__ = commit_id = 'g142ff7355'

{pycontrails-0.57.0 → pycontrails-0.58.0}/pycontrails/core/aircraft_performance.py

@@ -550,7 +550,7 @@ class AircraftPerformance(Model):
             if self.met is None:
                 cond = np.isnan(u) & np.isnan(v)
             else:
-                met_level_max = self.met.data["level"][-1].item()  # type: ignore[union-attr]
+                met_level_max = self.met.data["level"][-1].item()
                 cond = self.source.level > met_level_max
 
             # We DON'T overwrite the original u and v arrays already attached to the source

{pycontrails-0.57.0 → pycontrails-0.58.0}/pycontrails/core/cache.py

@@ -189,7 +189,7 @@ class DiskCacheStore(CacheStore):
         self,
         cache_dir: str | pathlib.Path | None = None,
         allow_clear: bool = False,
-    ):
+    ) -> None:
         if cache_dir is None:
             # Avoid unnecessary import of platformdirs (called in _get_user_cache_dir)
             cache_dir = os.getenv("PYCONTRAILS_CACHE_DIR") or _get_user_cache_dir()
@@ -461,7 +461,7 @@ class GCPCacheStore(CacheStore):
         timeout: int = 300,
         show_progress: bool = False,
         chunk_size: int = 64 * 262144,
-    ):
+    ) -> None:
         try:
             from google.cloud import storage
         except ModuleNotFoundError as e:

{pycontrails-0.57.0 → pycontrails-0.58.0}/pycontrails/core/fleet.py

@@ -5,12 +5,7 @@ from __future__ import annotations
 import sys
 import warnings
 from collections.abc import Iterable
-from typing import Any, NoReturn
-
-if sys.version_info >= (3, 11):
-    from typing import Self
-else:
-    from typing_extensions import Self
+from typing import Any, NoReturn, Self
 
 if sys.version_info >= (3, 12):
     from typing import override
@@ -122,7 +117,7 @@ class Fleet(Flight):
         # Set default fl_attrs if not provided
         fl_attrs = fl_attrs or {}
         for flight_id in groups.index:
-            fl_attrs.setdefault(flight_id, {})  # type: ignore[call-overload]
+            fl_attrs.setdefault(flight_id, {})
 
         extra = fl_attrs.keys() - groups.index
         if extra:

{pycontrails-0.57.0 → pycontrails-0.58.0}/pycontrails/core/flight.py

@@ -6,18 +6,13 @@ import enum
 import logging
 import sys
 import warnings
-from typing import TYPE_CHECKING, Any, NoReturn
+from typing import TYPE_CHECKING, Any, NoReturn, Self
 
 if sys.version_info >= (3, 12):
     from typing import override
 else:
     from typing_extensions import override
 
-if sys.version_info >= (3, 11):
-    from typing import Self
-else:
-    from typing_extensions import Self
-
 
 import numpy as np
 import numpy.typing as npt
@@ -2138,7 +2133,7 @@ def segment_rocd(
     T_correction[:-1] = (air_temperature[:-1] + air_temperature[1:]) / (T_isa[:-1] + T_isa[1:])
     T_correction[-1] = np.nan
 
-    return T_correction * out  # type: ignore[return-value]
+    return T_correction * out
 
 
 def _resample_to_freq_or_time(

{pycontrails-0.57.0 → pycontrails-0.58.0}/pycontrails/core/interpolation.py

@@ -26,57 +26,74 @@ class PycontrailsRegularGridInterpolator(scipy.interpolate.RegularGridInterpolat
 
     This class is a thin wrapper around the
     :class:`scipy.interpolate.RegularGridInterpolator` in order to make typical
-    ``pycontrails`` use-cases more efficient.
+    ``pycontrails`` linear interpolation use-cases more performant:
 
-    #. Avoid ``RegularGridInterpolator`` constructor validation. In :func:`interp`,
-    parameters are carefully crafted to fit into the intended form, thereby making
-    validation unnecessary.
+    #. Avoid ``RegularGridInterpolator`` constructor validation when `method="linear"`.
+       In :func:`interp`, parameters are carefully crafted to fit into the intended form,
+       thereby making validation unnecessary.
     #. Override the :meth:`_evaluate_linear` method with a faster implementation. See
-    the :meth:`_evaluate_linear` docstring for more information.
+       the :meth:`_evaluate_linear` docstring for more information.
 
-    This class should not be used directly. Instead, use the :func:`interp` function.
+    **This class should not be used directly. Instead, use the ``interp`` function.**
 
     .. versionchanged:: 0.40.0
 
         The :meth:`_evaluate_linear` method now uses a Cython implementation. The dtype
         of the output is now consistent with the dtype of the underlying :attr:`values`
 
+    .. versionchanged:: 0.58.0
+
+        Any ``method`` other than ``"linear"`` now uses the
+        :class:`scipy.interpolate.RegularGridInterpolator` implementation. This
+        allows for greater flexibility in the ``method`` parameter.
+
     Parameters
     ----------
     points : tuple[npt.NDArray[np.floating], ...]
         Coordinates of the grid points.
     values : npt.NDArray[np.floating]
         Grid values. The shape of this array must be compatible with the
-        coordinates. An error is raised if the dtype is not ``np.float32``
-        or ``np.float64``.
+        coordinates.
     method : str
         Passed into :class:`scipy.interpolate.RegularGridInterpolator`
     bounds_error : bool
         Passed into :class:`scipy.interpolate.RegularGridInterpolator`
     fill_value : float | np.float64 | None
         Passed into :class:`scipy.interpolate.RegularGridInterpolator`
+
+    See Also
+    --------
+    scipy.interpolate.RegularGridInterpolator
+    interp
     """
 
     def __init__(
         self,
         points: tuple[npt.NDArray[np.floating], ...],
         values: npt.NDArray[np.floating],
+        *,
         method: str,
         bounds_error: bool,
         fill_value: float | np.float64 | None,
-    ):
-        if values.dtype not in (np.float32, np.float64):
-            msg = f"values must be a float array, not {values.dtype}"
-            raise ValueError(msg)
-
+    ) -> None:
+        if method != "linear" or values.dtype not in (np.float32, np.float64):
+            # Slow path: use parent class
+            super().__init__(
+                points,
+                values,
+                method=method,
+                bounds_error=bounds_error,
+                fill_value=fill_value,
+            )
+            return
+
+        # Fast path: no validation
         self.grid = points
         self.values = values
-        # TODO: consider supporting updated tensor-product spline methods
-        # see https://github.com/scipy/scipy/releases/tag/v1.13.0
-        self.method = _pick_method(scipy.__version__, method)
+        self.method = method
         self.bounds_error = bounds_error
         self.fill_value = fill_value
-        self._spline = None
+        self._spline = None  # XXX: setting private attribute on RGI
 
     def _prepare_xi_simple(self, xi: npt.NDArray[np.floating]) -> npt.NDArray[np.bool_]:
         """Run looser version of :meth:`_prepare_xi`.
@@ -103,7 +120,7 @@ class PycontrailsRegularGridInterpolator(scipy.interpolate.RegularGridInterpolat
 
             return np.zeros(xi.shape[0], dtype=bool)
 
-        return self._find_out_of_bounds(xi.T)
+        return self._find_out_of_bounds(xi.T)  # XXX: calling private method on RGI
 
     def __call__(
         self, xi: npt.NDArray[np.floating], method: str | None = None
@@ -130,7 +147,7 @@ class PycontrailsRegularGridInterpolator(scipy.interpolate.RegularGridInterpolat
             return super().__call__(xi, method)
 
         out_of_bounds = self._prepare_xi_simple(xi)
-        xi_indices, norm_distances = rgi_cython.find_indices(self.grid, xi.T)
+        xi_indices, norm_distances = self._find_indices(xi.T)  # XXX: calling private method on RGI
 
         out = self._evaluate_linear(xi_indices, norm_distances)
         return self._set_out_of_bounds(out, out_of_bounds)
@@ -223,45 +240,6 @@ class PycontrailsRegularGridInterpolator(scipy.interpolate.RegularGridInterpolat
         raise ValueError(msg)
 
 
-def _pick_method(scipy_version: str, method: str) -> str:
-    """Select an interpolation method.
-
-    For scipy versions 1.13.0 and later, fall back on legacy implementations
-    of tensor-product spline methods. The default implementations in 1.13.0
-    and later are incompatible with this class.
-
-    Parameters
-    ----------
-    scipy_version : str
-        scipy version (major.minor.patch)
-
-    method : str
-        Interpolation method. Passed into :class:`scipy.interpolate.RegularGridInterpolator`
-        as-is unless ``scipy_version`` is 1.13.0 or later and ``method`` is ``"slinear"``,
-        ``"cubic"``, or ``"quintic"``. In this case, ``"_legacy"`` is appended to ``method``.
-
-    Returns
-    -------
-    str
-        Interpolation method adjusted for compatibility with this class.
-    """
-    if method == "linear":
-        return method
-
-    try:
-        version = scipy_version.split(".")
-        major = int(version[0])
-        minor = int(version[1])
-    except (IndexError, ValueError) as exc:
-        msg = f"Failed to parse major and minor version from {scipy_version}"
-        raise ValueError(msg) from exc
-
-    reimplemented_methods = ["slinear", "cubic", "quintic"]
-    if major > 1 or ((major == 1 and minor >= 13) and method in reimplemented_methods):
-        return method + "_legacy"
-    return method
-
-
 def _floatize_time(
     time: npt.NDArray[np.datetime64], offset: np.datetime64
 ) -> npt.NDArray[np.floating]:
@@ -431,7 +409,7 @@ def interp(
         Include ``indices`` and ``return_indices`` experimental parameters.
         Currently, nan values in ``longitude``, ``latitude``, ``level``, or ``time``
         are always propagated through to the output, regardless of ``bounds_error``.
-        In other words, a ValueError for an out of bounds coordinate is only raised
+        In other words, a ``ValueError`` for an out of bounds coordinate is only raised
         if a non-nan value is out of bounds.
 
     .. versionchanged:: 0.40.0
@@ -480,9 +458,9 @@ def interp(
 
     See Also
     --------
-    - :meth:`MetDataArray.interpolate`
-    - :func:`scipy.interpolate.interpn`
-    - :class:`scipy.interpolate.RegularGridInterpolator`
+    pycontrails.MetDataArray.interpolate
+    scipy.interpolate.interpn
+    scipy.interpolate.RegularGridInterpolator
     """
     if localize:
         coords = {"longitude": longitude, "latitude": latitude, "level": level, "time": time}
@@ -579,7 +557,7 @@ def _linear_interp_with_indices(
     if indices is None:
         assert xi is not None, "xi must be provided if indices is None"
         out_of_bounds = interp._prepare_xi_simple(xi)
-        xi_indices, norm_distances = rgi_cython.find_indices(interp.grid, xi.T)
+        xi_indices, norm_distances = interp._find_indices(xi.T)
         indices = RGIArtifacts(xi_indices, norm_distances, out_of_bounds)
 
     out = interp._evaluate_linear(indices.xi_indices, indices.norm_distances)
@@ -606,7 +584,7 @@ class EmissionsProfileInterpolator:
 
     This class simply wraps :func:`numpy.interp` with fixed values for the
     ``xp`` and ``fp`` arguments. Unlike :class:`xarray.DataArray` interpolation,
-    the `numpy.interp` automatically clips values outside the range of the
+    the :func:`numpy.interp` automatically clips values outside the range of the
     ``xp`` array.
 
     Parameters

{pycontrails-0.57.0 → pycontrails-0.58.0}/pycontrails/core/met.py

@@ -26,15 +26,11 @@ from typing import (
     Any,
     Generic,
     Literal,
+    Self,
     TypeVar,
     overload,
 )
 
-if sys.version_info >= (3, 11):
-    from typing import Self
-else:
-    from typing_extensions import Self
-
 if sys.version_info >= (3, 12):
     from typing import override
 else:
@@ -211,14 +207,22 @@ class MetBase(ABC, Generic[XArrayType]):
         Raises
         ------
         ValueError
-            If one of the coordinates is not sorted.
+            If one of the coordinates is not sorted or contains duplicate values.
         """
         indexes = self.indexes
-        if not np.all(np.diff(indexes["time"]) > np.timedelta64(0, "ns")):
-            raise ValueError("Coordinate 'time' not sorted. Instantiate with 'copy=True'.")
-        for coord in self.dim_order[:3]:  # exclude time, the 4th dimension
-            if not np.all(np.diff(indexes[coord]) > 0.0):
-                raise ValueError(f"Coordinate '{coord}' not sorted. Instantiate with 'copy=True'.")
+        for coord in self.dim_order:
+            arr = indexes[coord]
+            d = np.diff(arr)
+            zero = np.zeros((), dtype=d.dtype)  # ensure same dtype
+
+            if np.any(d <= zero):
+                if np.any(d == zero):
+                    msg = f"Coordinate '{coord}' contains duplicate values."
+                else:
+                    msg = f"Coordinate '{coord}' not sorted."
+
+                msg += " Instantiate with 'copy=True'."
+                raise ValueError(msg)
 
     def _validate_transpose(self) -> None:
         """Check that data is transposed according to :attr:`dim_order`."""
@@ -271,6 +275,10 @@ class MetBase(ABC, Generic[XArrayType]):
             Auxiliary coordinates (altitude and air_pressure) are now cast to the same
             dtype as the underlying grid data.
 
+        .. versionchanged:: 0.58.0
+
+            Duplicate dimension values are dropped, keeping the first occurrence.
+
 
         Parameters
         ----------
@@ -297,6 +305,18 @@ class MetBase(ABC, Generic[XArrayType]):
         # sortby to ensure each coordinate has ascending order
         self.data = self.data.sortby(list(self.dim_order), ascending=True)
 
+        # Drop any duplicated dimension values
+        indexes = self.indexes
+        for coord in self.dim_order:
+            arr = indexes[coord]
+            d = np.diff(arr)
+            zero = np.zeros((), dtype=d.dtype)  # ensure same dtype
+
+            if np.any(d == zero):
+                # Remove duplicates
+                filt = np.r_[True, d > zero]  # prepend True keeps the first occurrence
+                self.data = self.data.isel({coord: filt})
+
         if not self.is_wrapped:
             # Ensure longitude is contained in interval [-180, 180)
             # If longitude has value at 180, we might not want to shift it?
@@ -1087,7 +1107,7 @@ class MetDataset(MetBase):
             out[key] = da.values.ravel()  # type: ignore[index]
 
         if transfer_attrs:
-            out.attrs.update(self.attrs)  # type: ignore[arg-type]
+            out.attrs.update(self.attrs)
 
         return out
 
@@ -1300,7 +1320,7 @@ class MetDataset(MetBase):
         coords: dict[str, np.ndarray] = {}
         for key, val in input_data.items():
             dtype = "datetime64[ns]" if key == "time" else COORD_DTYPE
-            arr: np.ndarray = np.asarray(val, dtype=dtype)  # type: ignore[call-overload]
+            arr: np.ndarray = np.asarray(val, dtype=dtype)
 
             if arr.ndim == 0:
                 arr = arr.reshape(1)
@@ -1889,7 +1909,7 @@ class MetDataArray(MetBase):
         if not self.binary:
             raise NotImplementedError("proportion method is only implemented for binary fields")
 
-        return self.data.sum().values.item() / self.data.count().values.item()  # type: ignore[operator]
+        return self.data.sum().values.item() / self.data.count().values.item()
 
     def find_edges(self) -> Self:
         """Find edges of regions.
@@ -2598,9 +2618,9 @@ def _extract_2d_arr_and_altitude(
     except KeyError:
         altitude = None
     else:
-        altitude = round(altitude)  # type: ignore[call-overload]
+        altitude = round(altitude)
 
-    return arr, altitude  # type: ignore[return-value]
+    return arr, altitude
 
 
 def downselect(data: XArrayType, bbox: tuple[float, ...]) -> XArrayType:

{pycontrails-0.57.0 → pycontrails-0.58.0}/pycontrails/core/polygon.py

@@ -238,7 +238,7 @@ def _contours_to_polygons(
                 latitude=latitude,
                 precision=precision,
                 buffer=buffer,
-                i=child_i,  # type: ignore[arg-type]
+                i=child_i,
             )
 
             candidate = shapely.Polygon(polygon.exterior, [h.exterior for h in holes])
@@ -354,11 +354,11 @@ def find_multipolygon(
         return shapely.MultiPolygon()
 
     assert len(hierarchy) == 1
-    hierarchy = hierarchy[0]  # type: ignore[index]
+    hierarchy = hierarchy[0]
 
     polygons = _contours_to_polygons(
         contours,  # type: ignore[arg-type]
-        hierarchy,  # type: ignore[arg-type]
+        hierarchy,
         min_area,
         convex_hull,
         epsilon,