rio-tiler 7.9.2__py3-none-any.whl → 8.0.0__py3-none-any.whl

rio_tiler/mosaic/backend.py

@@ -0,0 +1,242 @@
+"""rio_tiler.mosaic.backend: base Backend class."""
+
+import abc
+import logging
+from typing import Any, Type
+
+import attr
+from morecantile import TileMatrixSet
+from pydantic import BaseModel, ConfigDict, Field
+from rasterio.crs import CRS
+from rasterio.features import bounds as featureBounds
+from rasterio.features import geometry_mask
+
+from rio_tiler.constants import WEB_MERCATOR_TMS, WGS84_CRS
+from rio_tiler.errors import NoAssetFoundError, PointOutsideBounds
+from rio_tiler.io import BaseReader, MultiBandReader, MultiBaseReader, Reader
+from rio_tiler.models import ImageData, PointData
+from rio_tiler.mosaic import mosaic_reader
+from rio_tiler.tasks import multi_values_list
+from rio_tiler.types import BBox
+from rio_tiler.utils import CRS_to_uri, Timer, _validate_shape_input
+
+logger = logging.getLogger(__name__)
+
+
+class MosaicInfo(BaseModel):
+    """Mosaic info responses."""
+
+    bounds: BBox = Field(default=(-180, -90, 180, 90))
+    crs: str
+
+    model_config = ConfigDict(extra="allow")
+
+
+@attr.s
+class BaseBackend(BaseReader):
+    """Base Class for mosaic backend.
+
+    Attributes:
+        input (str): mosaic path.
+        tms (morecantile.TileMatrixSet, optional): TileMatrixSet grid definition. Defaults to `WebMercatorQuad`.
+        reader (rio_tiler.io.BaseReader): Dataset reader. Defaults to `rio_tiler.io.Reader`.
+        reader_options (dict): Options to forward to the reader config.
+        bounds (tuple): mosaic bounds (left, bottom, right, top). **READ ONLY attribute**.
+        crs (rasterio.crs.CRS): mosaic crs in which its bounds is defined. **READ ONLY attribute**.
+        minzoom (int): mosaic minimum zoom level. **READ ONLY attribute**.
+        maxzoom (int): mosaic maximum zoom level. **READ ONLY attribute**.
+
+    """
+
+    input: str = attr.ib()
+    tms: TileMatrixSet = attr.ib(default=WEB_MERCATOR_TMS)
+
+    reader: Type[BaseReader] | Type[MultiBaseReader] | Type[MultiBandReader] = attr.ib(
+        default=Reader
+    )
+    reader_options: dict = attr.ib(factory=dict)
+
+    bounds: BBox = attr.ib(init=False)
+    crs: CRS = attr.ib(init=False)
+
+    minzoom: int = attr.ib(init=False)
+    maxzoom: int = attr.ib(init=False)
+
+    @abc.abstractmethod
+    def assets_for_tile(self, x: int, y: int, z: int, **kwargs: Any) -> list[Any]:
+        """Retrieve assets for tile."""
+
+    @abc.abstractmethod
+    def assets_for_point(
+        self,
+        lng: float,
+        lat: float,
+        coord_crs: CRS | None = None,
+        **kwargs: Any,
+    ) -> list[Any]:
+        """Retrieve assets for point."""
+
+    @abc.abstractmethod
+    def assets_for_bbox(
+        self,
+        xmin: float,
+        ymin: float,
+        xmax: float,
+        ymax: float,
+        coord_crs: CRS | None = None,
+        **kwargs,
+    ) -> list[Any]:
+        """Retrieve assets for bbox."""
+
+    def asset_name(self, asset: Any) -> str:
+        """Get asset name."""
+        return str(asset)
+
+    def info(self) -> MosaicInfo:  # type: ignore
+        """Mosaic info."""
+        return MosaicInfo(
+            bounds=self.bounds,
+            crs=CRS_to_uri(self.crs) or self.crs.to_wkt(),
+        )
+
+    def point(  # type: ignore
+        self,
+        lon: float,
+        lat: float,
+        coord_crs: CRS = WGS84_CRS,
+        search_options: dict | None = None,
+        **kwargs: Any,
+    ) -> list[tuple[str, PointData]]:
+        """Get Point value from multiple assets."""
+        search_options = search_options or {}
+        mosaic_assets = self.assets_for_point(
+            lon, lat, coord_crs=coord_crs, **search_options
+        )
+        if not mosaic_assets:
+            raise NoAssetFoundError(f"No assets found for point ({lon},{lat})")
+
+        def _reader(
+            asset: Any, lon: float, lat: float, coord_crs: CRS, **kwargs
+        ) -> PointData:
+            with self.reader(asset, **self.reader_options) as src_dst:
+                return src_dst.point(lon, lat, coord_crs=coord_crs, **kwargs)
+
+        if "allowed_exceptions" not in kwargs:
+            kwargs.update({"allowed_exceptions": (PointOutsideBounds,)})
+
+        logger.info(
+            f"reading Point for {len(mosaic_assets)} assets with reader: {self.reader}"
+        )
+        return [
+            (self.asset_name(asset), pt)
+            for asset, pt in multi_values_list(
+                mosaic_assets, _reader, lon, lat, coord_crs, **kwargs
+            )
+        ]
+
+    def tile(  # type: ignore
+        self,
+        x: int,
+        y: int,
+        z: int,
+        search_options: dict | None = None,
+        **kwargs: Any,
+    ) -> tuple[ImageData, list[str]]:
+        """Get Tile from multiple assets."""
+        timings = []
+        with Timer() as t:
+            search_options = search_options or {}
+            mosaic_assets = self.assets_for_tile(x, y, z, **search_options)
+        timings.append(("search", round(t.elapsed * 1000, 2)))
+
+        if not mosaic_assets:
+            raise NoAssetFoundError(f"No assets found for tile {z}-{x}-{y}")
+
+        def _reader(asset: Any, x: int, y: int, z: int, **kwargs: Any) -> ImageData:
+            with self.reader(asset, tms=self.tms, **self.reader_options) as src_dst:
+                return src_dst.tile(x, y, z, **kwargs)
+
+        logger.info(
+            f"reading Tile for {len(mosaic_assets)} assets with reader: {self.reader}"
+        )
+        with Timer() as t:
+            img, asset_used = mosaic_reader(mosaic_assets, _reader, x, y, z, **kwargs)
+        timings.append(("mosaicking", round(t.elapsed * 1000, 2)))
+        img.metadata = {**img.metadata, "timings": timings}
+
+        asset_used = [self.asset_name(asset) for asset in asset_used]
+        return img, asset_used
+
+    def part(  # type: ignore
+        self,
+        bbox: BBox,
+        bounds_crs: CRS = WGS84_CRS,
+        search_options: dict | None = None,
+        **kwargs: Any,
+    ) -> tuple[ImageData, list[str]]:
+        """Create an Image from multiple assets for a bbox."""
+        xmin, ymin, xmax, ymax = bbox
+        timings = []
+
+        with Timer() as t:
+            search_options = search_options or {}
+            mosaic_assets = self.assets_for_bbox(
+                xmin,
+                ymin,
+                xmax,
+                ymax,
+                coord_crs=bounds_crs,
+                **search_options,
+            )
+        timings.append(("search", round(t.elapsed * 1000, 2)))
+
+        if not mosaic_assets:
+            raise NoAssetFoundError("No assets found for bbox input")
+
+        def _reader(asset: Any, bbox: BBox, bounds_crs: CRS, **kwargs: Any) -> ImageData:
+            with self.reader(asset, **self.reader_options) as src_dst:
+                return src_dst.part(bbox, bounds_crs=bounds_crs, **kwargs)
+
+        logger.info(
+            f"reading Part for {len(mosaic_assets)} assets with reader: {self.reader}"
+        )
+        with Timer() as t:
+            img, asset_used = mosaic_reader(
+                mosaic_assets, _reader, bbox, bounds_crs, **kwargs
+            )
+        timings.append(("mosaicking", round(t.elapsed * 1000, 2)))
+        img.metadata = {**img.metadata, "timings": timings}
+
+        asset_used = [self.asset_name(asset) for asset in asset_used]
+        return img, asset_used
+
+    def feature(  # type: ignore
+        self,
+        shape: dict,
+        shape_crs: CRS = WGS84_CRS,
+        search_options: dict | None = None,
+        **kwargs: Any,
+    ) -> tuple[ImageData, list[str]]:
+        """Create an Image from multiple assets for a GeoJSON feature."""
+        shape = _validate_shape_input(shape)
+        bbox = featureBounds(shape)
+
+        img, asset_used = self.part(
+            bbox,
+            bounds_crs=shape_crs,
+            search_options=search_options,
+            **kwargs,
+        )
+        img.array.mask = geometry_mask([shape], (img.height, img.width), img.transform)
+        return img, asset_used
+
+    ############################################################################
+    # Not Implemented methods
+    # BaseReader required those method to be implemented
+    def statistics(self):
+        """PlaceHolder for statistics."""
+        raise NotImplementedError
+
+    def preview(self):
+        """PlaceHolder for preview."""
+        raise NotImplementedError

rio_tiler/mosaic/reader.py

@@ -95,6 +95,7 @@ def mosaic_reader(  # noqa: C901
                 crs = img.crs
                 bounds = img.bounds
                 band_names = img.band_names
+                band_descriptions = img.band_descriptions
                 pixel_selection.cutline_mask = img.cutline_mask
                 pixel_selection.width = img.width
                 pixel_selection.height = img.height
@@ -135,6 +136,7 @@ def mosaic_reader(  # noqa: C901
                         crs=crs,
                         bounds=bounds,
                         band_names=band_names,
+                        band_descriptions=band_descriptions,
                         metadata={
                             "mosaic_method": pixel_selection.__class__.__name__,
                             "mosaic_assets_count": len(mosaic_assets),
@@ -154,6 +156,7 @@ def mosaic_reader(  # noqa: C901
             crs=crs,
             bounds=bounds,
             band_names=band_names,
+            band_descriptions=band_descriptions,
             metadata={
                 "mosaic_method": pixel_selection.__class__.__name__,
                 "mosaic_assets_count": len(mosaic_assets),
@@ -229,6 +232,7 @@ def mosaic_point_reader(
                 crs = pt.crs
                 coordinates = pt.coordinates
                 band_names = pt.band_names
+                band_descriptions = pt.band_descriptions
                 pixel_selection.width = 1
                 pixel_selection.height = 1
                 pixel_selection.count = pt.count
@@ -248,6 +252,7 @@ def mosaic_point_reader(
                         crs=crs,
                         coordinates=coordinates,
                         band_names=band_names,
+                        band_descriptions=band_descriptions,
                         metadata={
                             "mosaic_method": pixel_selection.__class__.__name__,
                             "mosaic_assets_count": len(mosaic_assets),
@@ -267,6 +272,7 @@ def mosaic_point_reader(
             crs=crs,
             coordinates=coordinates,
             band_names=band_names,
+            band_descriptions=band_descriptions,
             metadata={
                 "mosaic_method": pixel_selection.__class__.__name__,
                 "mosaic_assets_count": len(mosaic_assets),

rio_tiler/reader.py

@@ -34,7 +34,6 @@ from rio_tiler.utils import (
 class Options(TypedDict, total=False):
     """Reader Options."""
 
-    force_binary_mask: Optional[bool]
     nodata: Optional[NoData]
     vrt_options: Optional[Dict]
     resampling_method: Optional[RIOResampling]
@@ -76,7 +75,6 @@ def read(
     max_size: Optional[int] = None,
     indexes: Optional[Indexes] = None,
     window: Optional[windows.Window] = None,
-    force_binary_mask: bool = True,
     nodata: Optional[NoData] = None,
     vrt_options: Optional[Dict] = None,
     out_dtype: Optional[Union[str, numpy.dtype]] = None,
@@ -99,7 +97,6 @@ def read(
         vrt_options (dict, optional): Options to be passed to the rasterio.warp.WarpedVRT class.
         resampling_method (RIOResampling, optional): RasterIO resampling algorithm. Defaults to `nearest`.
         reproject_method (WarpResampling, optional): WarpKernel resampling algorithm. Defaults to `nearest`.
-        force_binary_mask (bool, optional): Cast returned mask to binary values (0 or 255). Defaults to `True`.
         unscale (bool, optional): Apply 'scales' and 'offsets' on output data value. Defaults to `False`.
         post_process (callable, optional): Function to apply on output data and mask values.
 
@@ -258,37 +255,33 @@ def read(
         # We only add dataset statistics if we have them for all the indexes
         dataset_statistics = stats if len(stats) == len(indexes) else None
 
-        # TODO: DEPRECATED, masked array are already using bool
-        if force_binary_mask:
-            pass
-
+        scales = numpy.array(dataset.scales)[numpy.array(indexes) - 1]
+        offsets = numpy.array(dataset.offsets)[numpy.array(indexes) - 1]
         if unscale:
             data = data.astype("float32", casting="unsafe")
-
-            # reshaped to match data
-            scales = numpy.array(dataset.scales)[numpy.array(indexes) - 1].reshape(
-                (-1, 1, 1)
-            )
-            offsets = numpy.array(dataset.offsets)[numpy.array(indexes) - 1].reshape(
-                (-1, 1, 1)
-            )
-
-            numpy.multiply(data, scales, out=data, casting="unsafe")
-            numpy.add(data, offsets, out=data, casting="unsafe")
+            numpy.multiply(data, scales.reshape((-1, 1, 1)), out=data, casting="unsafe")
+            numpy.add(data, offsets.reshape((-1, 1, 1)), out=data, casting="unsafe")
 
             # apply scale/offsets to stats
             if dataset_statistics:
-                scales = numpy.array(dataset.scales)[numpy.array(indexes) - 1].reshape(
-                    (-1, 1)
+                stats_array = numpy.array(dataset_statistics)
+                numpy.multiply(
+                    stats_array,
+                    scales.reshape((-1, 1)),
+                    out=stats_array,
+                    casting="unsafe",
                 )
-                offsets = numpy.array(dataset.offsets)[numpy.array(indexes) - 1].reshape(
-                    (-1, 1)
+                numpy.add(
+                    stats_array,
+                    offsets.reshape((-1, 1)),
+                    out=stats_array,
+                    casting="unsafe",
                 )
-                stats_array = numpy.array(dataset_statistics)
-                numpy.multiply(stats_array, scales, out=stats_array, casting="unsafe")
-                numpy.add(stats_array, offsets, out=stats_array, casting="unsafe")
                 dataset_statistics = [tuple(s) for s in stats_array.tolist()]
 
+            scales = numpy.zeros(len(indexes)) + 1.0
+            offsets = numpy.zeros(len(indexes))
+
         if post_process:
             data = post_process(data)
 
@@ -301,8 +294,12 @@ def read(
             bounds=out_bounds,
             crs=dataset.crs,
             band_names=[f"b{idx}" for idx in indexes],
+            band_descriptions=[dataset.descriptions[ix - 1] or "" for idx in indexes],
             dataset_statistics=dataset_statistics,
             metadata=dataset.tags(),
+            nodata=nodata,
+            scales=scales.tolist(),
+            offsets=offsets.tolist(),
         )
 
 
@@ -319,7 +316,6 @@ def part(
     minimum_overlap: Optional[float] = None,
     padding: Optional[int] = None,
     buffer: Optional[float] = None,
-    force_binary_mask: bool = True,
     nodata: Optional[NoData] = None,
     vrt_options: Optional[Dict] = None,
     out_dtype: Optional[Union[str, numpy.dtype]] = None,
@@ -456,7 +452,6 @@ def part(
             out_dtype=out_dtype,
             resampling_method=resampling_method,
             reproject_method=reproject_method,
-            force_binary_mask=force_binary_mask,
             unscale=unscale,
             post_process=post_process,
         )
@@ -501,7 +496,6 @@ def part(
             out_dtype=out_dtype,
             resampling_method=resampling_method,
             reproject_method=reproject_method,
-            force_binary_mask=force_binary_mask,
             unscale=unscale,
             post_process=post_process,
         )
@@ -511,6 +505,10 @@ def part(
             bounds=bounds,
             crs=img.crs,
             band_names=img.band_names,
+            band_descriptions=img.band_descriptions,
+            nodata=img.nodata,
+            scales=img.scales,
+            offsets=img.offsets,
             dataset_statistics=img.dataset_statistics,
             metadata=img.metadata,
         )
@@ -525,7 +523,6 @@ def part(
         out_dtype=out_dtype,
         resampling_method=resampling_method,
         reproject_method=reproject_method,
-        force_binary_mask=force_binary_mask,
         unscale=unscale,
         post_process=post_process,
     )
@@ -536,7 +533,6 @@ def point(
     coordinates: Tuple[float, float],
     indexes: Optional[Indexes] = None,
     coord_crs: CRS = WGS84_CRS,
-    force_binary_mask: bool = True,
     nodata: Optional[NoData] = None,
     vrt_options: Optional[Dict] = None,
     out_dtype: Optional[Union[str, numpy.dtype]] = None,
@@ -642,7 +638,6 @@ def point(
             window=window,
             out_dtype=out_dtype,
             resampling_method=resampling_method,
-            force_binary_mask=force_binary_mask,
             unscale=unscale,
             post_process=post_process,
         )
@@ -650,8 +645,12 @@ def point(
         return PointData(
             img.array[:, 0, 0],
             band_names=img.band_names,
+            band_descriptions=img.band_descriptions,
             coordinates=coordinates,
             crs=coord_crs,
-            metadata=dataset.tags(),
             pixel_location=(col, row),
+            nodata=img.nodata,
+            scales=img.scales,
+            offsets=img.offsets,
+            metadata=img.metadata,
         )

rio_tiler/tasks.py

@@ -99,3 +99,19 @@ def multi_values(
         asset: val
         for val, asset in filter_tasks(tasks, allowed_exceptions=allowed_exceptions)
     }
+
+
+def multi_values_list(
+    asset_list: Sequence,
+    reader: Callable,
+    *args: Any,
+    threads: int = MAX_THREADS,
+    allowed_exceptions: Optional[Tuple] = None,
+    **kwargs: Any,
+) -> list[tuple[Any, Any]]:
+    """Merge values returned from tasks."""
+    tasks = create_tasks(reader, asset_list, threads, *args, **kwargs)
+    return [
+        (asset, val)
+        for val, asset in filter_tasks(tasks, allowed_exceptions=allowed_exceptions)
+    ]

rio_tiler/utils.py

@@ -1,6 +1,7 @@
 """rio_tiler.utils: utility functions."""
 
 import math
+import time
 import warnings
 from io import BytesIO
 from typing import (
@@ -517,7 +518,7 @@ def linear_rescale(
     """
     imin, imax = in_range
     omin, omax = out_range
-    im = numpy.clip(image, imin, imax) - imin
+    im = numpy.clip(image, imin, imax, dtype=numpy.float64) - imin
     im = im / numpy.float64(imax - imin)
     return im * (omax - omin) + omin
 
@@ -893,3 +894,30 @@ def CRS_to_urn(crs: CRS) -> Optional[str]:
         return f"urn:ogc:def:crs:{authority}:{version}:{code}"
 
     return None
+
+
+# This code is copied from marblecutter
+#  https://github.com/mojodna/marblecutter/blob/master/marblecutter/stats.py
+# License:
+# Original work Copyright 2016 Stamen Design
+# Modified work Copyright 2016-2017 Seth Fitzsimmons
+# Modified work Copyright 2016 American Red Cross
+# Modified work Copyright 2016-2017 Humanitarian OpenStreetMap Team
+# Modified work Copyright 2017 Mapzen
+class Timer(object):
+    """Time a code block."""
+
+    def __enter__(self):
+        """Starts timer."""
+        self.start = time.time()
+        return self
+
+    def __exit__(self, ty, val, tb):
+        """Stops timer."""
+        self.end = time.time()
+        self.elapsed = self.end - self.start
+
+    @property
+    def from_start(self):
+        """Return time elapsed from start."""
+        return time.time() - self.start

{rio_tiler-7.9.2.dist-info → rio_tiler-8.0.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: rio-tiler
-Version: 7.9.2
+Version: 8.0.0
 Summary: User friendly Rasterio plugin to read raster datasets.
 Project-URL: Homepage, https://cogeotiff.github.io/rio-tiler/
 Project-URL: Documentation, https://cogeotiff.github.io/rio-tiler/
@@ -43,54 +43,30 @@ Keywords: COGEO,Cloud Optimized Geotiff,STAC,rasterio,slippy-map
 Classifier: Intended Audience :: Information Technology
 Classifier: Intended Audience :: Science/Research
 Classifier: License :: OSI Approved :: BSD License
-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.9
+Requires-Python: >=3.11
 Requires-Dist: attrs
 Requires-Dist: cachetools
 Requires-Dist: color-operations
 Requires-Dist: httpx
-Requires-Dist: morecantile<7.0,>=5.0
+Requires-Dist: morecantile<8.0,>=5.0
 Requires-Dist: numexpr
 Requires-Dist: numpy
 Requires-Dist: pydantic~=2.0
 Requires-Dist: pystac<2.0,>=1.9
 Requires-Dist: rasterio>=1.4.0
 Requires-Dist: typing-extensions
-Provides-Extra: benchmark
-Requires-Dist: pytest; extra == 'benchmark'
-Requires-Dist: pytest-benchmark; extra == 'benchmark'
-Provides-Extra: dev
-Requires-Dist: bump-my-version; extra == 'dev'
-Requires-Dist: pre-commit; extra == 'dev'
-Provides-Extra: docs
-Requires-Dist: griffe-inherited-docstrings>=1.0.0; extra == 'docs'
-Requires-Dist: mkdocs-jupyter>=0.24.5; extra == 'docs'
-Requires-Dist: mkdocs-material[imaging]>=9.5; extra == 'docs'
-Requires-Dist: mkdocs>=1.4.3; extra == 'docs'
-Requires-Dist: mkdocstrings[python]>=0.25.1; extra == 'docs'
-Requires-Dist: pygments; extra == 'docs'
 Provides-Extra: s3
 Requires-Dist: boto3; extra == 's3'
-Provides-Extra: test
-Requires-Dist: boto3; extra == 'test'
-Requires-Dist: h5netcdf; extra == 'test'
-Requires-Dist: morecantile<7.0,>=6.0; extra == 'test'
-Requires-Dist: pytest; extra == 'test'
-Requires-Dist: pytest-cov; extra == 'test'
-Requires-Dist: rioxarray; extra == 'test'
-Requires-Dist: vsifile>=0.2; extra == 'test'
-Requires-Dist: xarray; extra == 'test'
-Provides-Extra: tilebench
-Requires-Dist: pytest; extra == 'tilebench'
-Requires-Dist: tilebench; extra == 'tilebench'
 Provides-Extra: xarray
 Requires-Dist: rioxarray; extra == 'xarray'
 Requires-Dist: xarray; extra == 'xarray'
+Provides-Extra: zarr
+Requires-Dist: obstore; extra == 'zarr'
+Requires-Dist: zarr~=3.0; extra == 'zarr'
 Description-Content-Type: text/markdown
 
 # rio-tiler
@@ -264,17 +240,17 @@ At the low level, `rio-tiler` is *just* a wrapper around the [rasterio](https://
 You can install `rio-tiler` using pip
 
 ```bash
-$ python -m pip install -U pip
-$ python -m pip install -U rio-tiler
+python -m pip install -U pip
+python -m pip install -U rio-tiler
 ```
 
 or install from source:
 
 ```bash
-$ git clone https://github.com/cogeotiff/rio-tiler.git
-$ cd rio-tiler
-$ python -m pip install -U pip
-$ python -m pip install -e .
+git clone https://github.com/cogeotiff/rio-tiler.git
+cd rio-tiler
+python -m pip install -U pip
+python -m pip install -e .
 ```
 
 ## Plugins

{rio_tiler-7.9.2.dist-info → rio_tiler-8.0.0.dist-info}/RECORD

@@ -1,16 +1,16 @@
-rio_tiler/__init__.py,sha256=cZthX4uSnrRjzSzlllIPZcsJZZDBznuvk3N6NP7saaU,192
-rio_tiler/colormap.py,sha256=E2FYp6AysmnSAhj2901XUrgB2q-8OvO8zoe_hmM_UZU,11186
+rio_tiler/__init__.py,sha256=lKZDEOYy6ewX4dX9hskgOUWtXPSxhBHJ4fzw3Yvf6Zc,192
+rio_tiler/colormap.py,sha256=4sG9VIfDAHI-ckQgmci0-U-wL9ZU9exZJszgsrt8wBA,11310
 rio_tiler/constants.py,sha256=55i-7JZDupTXZdLgxL03KsgM4lAzuGuIVP1zZKktzp0,426
-rio_tiler/errors.py,sha256=dnaAHPc5VL7SNKXm00gzlrU_4XIuw5rwzD-9Q2kFJsc,2018
-rio_tiler/expression.py,sha256=EVHWgjMqIn4TTPfqT_HV1RWCQlxqOj3OtEr1iUDiFHE,2282
+rio_tiler/errors.py,sha256=GFAuE1AaSvx6dd0z5um9n3y1wVzUyQ5S8qY7_tXdrR8,2178
+rio_tiler/expression.py,sha256=8StB3kqW0EqsW48mu8tlnFNgq48Cyuaj0I-ahmg27Uc,2305
 rio_tiler/logger.py,sha256=RR8lnW3uVXkFkPa3nNJS_tTndmdiNNDVXpCDGDxGf0A,81
-rio_tiler/models.py,sha256=_4JaeOrGInIKJqLSWhM_EZPoB8X8RRPSrR1FIel7Lmk,29293
+rio_tiler/models.py,sha256=ABXuGNTtVI1bBZ90AdyXigwug2P7kg9lrDbLss0h61w,34921
 rio_tiler/profiles.py,sha256=EAx2JdcaOcMw5PZjxbCqQBXXWMac9mjtpHoVFPJEDNQ,1562
 rio_tiler/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-rio_tiler/reader.py,sha256=ELTCGRwl53yKAk1k-KvBD25GXnb1xEsSHh1ajSoVSH0,25170
-rio_tiler/tasks.py,sha256=tiQo24rZHBVxids40QkNInuMzDFEei08zYW3I_pkSVE,3157
+rio_tiler/reader.py,sha256=5JVlhJS8Ra8jPPVkjkRw2UuXOnG7mPtKxYlDC9vF6qI,25047
+rio_tiler/tasks.py,sha256=7MwIicDMSCScsbds6crWRa7lRwvvt2cbGQAxybkrRvk,3613
 rio_tiler/types.py,sha256=R9XOL-EEb_cwbD4SFNVmViXC84D2p4BEVNCXg0nlh7I,1684
-rio_tiler/utils.py,sha256=Xcx7Ycdy52A3t0OrMIw5lMC2t-rWqZke3iySqsuCXjk,31314
+rio_tiler/utils.py,sha256=BwTHM_qmwRl0QEHwf7qQa22zdQPWeuhUHzCDqlUDkU8,32156
 rio_tiler/cmap_data/__init__.py,sha256=8FtVmfpTjXlvhxQ5QesN0UC1m_B3MuF3LbGbhMC5Rw4,1039
 rio_tiler/cmap_data/accent.npy,sha256=Qde1ldOoXghe4L05v1QbVvnMA1ldwNjKWPf5xCBbmI4,1152
 rio_tiler/cmap_data/accent_r.npy,sha256=ba-GWSMSPRAcm9CGzlXJeNG4ABbBHDCV602hAWV2idc,1152
@@ -225,18 +225,20 @@ rio_tiler/cmap_data/ylorrd.npy,sha256=9ImXljw40oe60w8uV4EMDPY4aFFVkGbyCBi6SlTX83
 rio_tiler/cmap_data/ylorrd_r.npy,sha256=K5uiHNHbLxV5SizyT09cSVAxldE-BW5GpOXxUp7UsTE,1152
 rio_tiler/experimental/__init__.py,sha256=dgW817h5U9OQt4V68L912g3DivPfhoWKvS-mI83sqio,328
 rio_tiler/experimental/vsifile.py,sha256=t1tJviPmFkjWIbr86BEnsK7pU4iiqLh1pKYZyNNb1AY,960
+rio_tiler/experimental/zarr.py,sha256=S9PJcv1JwoMI3sBRg0z2uEj5rk8KvwlR0-k7kA27rJ4,11775
 rio_tiler/io/__init__.py,sha256=_L4iILm6vSiJ14GEDDOvkuUHRtbWC9oqx6Bu8PxHhvA,270
-rio_tiler/io/base.py,sha256=nyaAtfodM-5fXd1WsGG5Qh544PNv9Hy4FLQ_nPj8hw8,52154
-rio_tiler/io/rasterio.py,sha256=4j4ffxsBHD7wg2R1p6x6JE_7UZjMZQGx29qtO6JEs8U,29636
+rio_tiler/io/base.py,sha256=3aKWc_1eGe_6Gzn5v8XOmZqUodfYsgytVLvNISqi5-g,53199
+rio_tiler/io/rasterio.py,sha256=JcjHJYmntBx4yyLOnTOJsZ3631TNsXFHbVGwK2MPnS8,29190
 rio_tiler/io/stac.py,sha256=kIEW4F71PouXF-Ubpz-VVXujnp8LqftptPKDo_J3BTw,12831
-rio_tiler/io/xarray.py,sha256=geOGIyhwlKdG4jguuW5AQl657sPmNGRrSJFSk5UEs-Q,26401
+rio_tiler/io/xarray.py,sha256=jISNXCV2lBjs5qOYSLvSfrD7v6fO69sdsb7KC3fatUQ,26915
 rio_tiler/mosaic/__init__.py,sha256=Yj6CKpnFl8PJhLSp-a55wo33hKZ8-6OOBJtWA1HZVy8,118
-rio_tiler/mosaic/reader.py,sha256=KSk6DaUJRN4Ye8lJ-V12sl5i2FN6y9o4wS1GzVSmtgU,10694
+rio_tiler/mosaic/backend.py,sha256=OI4hOM2No30gQVSIyAAwgx6ntEvw0Rk9neTeOmwstrw,8200
+rio_tiler/mosaic/reader.py,sha256=4sgl8IlV7zOtcfCt6XWw5UZLMCP0EgNbR2vRXObEGDA,11029
 rio_tiler/mosaic/methods/__init__.py,sha256=tgkXM9skaTLXIm5QFoheOEznQXM97KGflcAWHfkrt1g,612
 rio_tiler/mosaic/methods/base.py,sha256=9YZJWVRwH5Fk9KO9q5CW52Q8Mf60tAJ21oM4ixEDXBo,1424
 rio_tiler/mosaic/methods/defaults.py,sha256=P-zSVVgvQ2oRiMchC3bMqIYLTX7H3743jsD_VmqStxM,7642
-rio_tiler-7.9.2.dist-info/METADATA,sha256=h03DPeFZMtPvtd3etzxHFLhm2XeRShbk8eTCK9MWD9k,12233
-rio_tiler-7.9.2.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
-rio_tiler-7.9.2.dist-info/licenses/AUTHORS.txt,sha256=FCVd4Tjg-8syl0ZugCunpXER8X2-XonW2ZfllyTnRvE,158
-rio_tiler-7.9.2.dist-info/licenses/LICENSE,sha256=vq8Tt4KoYQT9JxAjQ4yXMmmhFYRTsBRgrOj-ao-bC5o,1517
-rio_tiler-7.9.2.dist-info/RECORD,,
+rio_tiler-8.0.0.dist-info/METADATA,sha256=wIZwjY9mXGOW3o3KvcLq-UreK9LX433GBkdvbhFv6yM,11159
+rio_tiler-8.0.0.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+rio_tiler-8.0.0.dist-info/licenses/AUTHORS.txt,sha256=FCVd4Tjg-8syl0ZugCunpXER8X2-XonW2ZfllyTnRvE,158
+rio_tiler-8.0.0.dist-info/licenses/LICENSE,sha256=vq8Tt4KoYQT9JxAjQ4yXMmmhFYRTsBRgrOj-ao-bC5o,1517
+rio_tiler-8.0.0.dist-info/RECORD,,