update for rio-tiler 8.0 (#1272)

* update for rio-tiler 8.0

* update changelog

* update extensions

* update xarray

* update mosaic

* update mosaic and app

Author: vincentsarago

.gitignore

@@ -108,4 +108,6 @@ ENV/
 cdk.out/
 deployment/k8s/titiler/values-test.yaml
 
-.vscode/
+.vscode/
+
+dev_notebooks/*

CHANGES.md

@@ -1,5 +1,41 @@
 # Release Notes
 
+## Unreleased
+
+### Misc 
+
+* update rio-tiler requirement to `>=8.0,<9.0`
+
+### titiler.core
+
+* add `band_description` attribute to `Point` output model (returned by /point endpoints)  **breaking change**
+
+### titiler.extensions
+
+* update rio-cogeo requirement to `7.0,<8.0`
+
+### titiler.mosaic
+
+* change Response model for `/point` endpoint  **breaking change**
+
+    ```python
+    # before
+    class Point(BaseModel):
+        coordinates: List[float]
+        values: List[Tuple[str, List[Optional[float]], List[str]]]
+
+    # now
+    class AssetPoint(BaseModel):
+        name: str
+        values: list[float | None]
+        band_names: list[str]
+        band_descriptions: list[str] | None = None
+
+    class Point(BaseModel):
+        coordinates: list[float]
+        assets: list[AssetPoint]
+    ```
+
 ## 0.26.0 (2025-11-25)
 
 ### titiler.xarray

pyproject.toml

@@ -49,7 +49,6 @@ dev = [
     "s3fs>=2025.2.0",
     "aiohttp",
     "requests",
-    "cogeo-mosaic>=8.2,<9.0",
     "pystac[validation]>=1.0.0,<2.0.0",
     "brotlipy",
     "boto3",
@@ -74,6 +73,12 @@ docs = [
 server = [
     "uvicorn",
 ]
+notebook = [
+    "IPython",
+    "jupyter",
+    "matplotlib",
+    "folium",
+]
 
 [project.urls]
 Homepage = 'https://developmentseed.org/titiler/'

src/titiler/application/tests/routes/test_mosaic.py

@@ -47,7 +47,6 @@ def test_info(app):
     response = app.get("/mosaicjson/info", params={"url": MOSAICJSON_FILE})
     assert response.status_code == 200
     body = response.json()
-    assert body["name"] == "mosaic"  # mosaic.name is not set
     assert body["quadkeys"] == []
     assert body["mosaic_minzoom"] == 7
     assert body["mosaic_maxzoom"] == 9
@@ -58,7 +57,6 @@ def test_info(app):
     assert response.headers["content-type"] == "application/geo+json"
     body = response.json()
     assert body["geometry"]
-    assert body["properties"]["name"] == "mosaic"  # mosaic.name is not set
     assert body["properties"]["quadkeys"] == []
     assert body["properties"]["mosaic_minzoom"] == 7
     assert body["properties"]["mosaic_maxzoom"] == 9
@@ -95,9 +93,9 @@ def test_point(app):
         )
     assert response.status_code == 200
     body = response.json()
-    assert len(body["values"]) == 1
-    assert body["values"][0][0].endswith(".tif")
-    assert body["values"][0][1] == [9943, 9127, 9603]
+    assert len(body["assets"]) == 1
+    assert body["assets"][0]["name"].endswith(".tif")
+    assert body["assets"][0]["values"] == [9943, 9127, 9603]
 
 
 def test_tile(app):

src/titiler/core/pyproject.toml

@@ -35,7 +35,7 @@ dependencies = [
     "numpy",
     "pydantic~=2.0",
     "rasterio",
-    "rio-tiler>=7.7,<8.0",
+    "rio-tiler>=8.0,<9.0",
     "morecantile",
     "simplejson",
     "typing_extensions>=4.6.1",

src/titiler/core/tests/test_rendering.py

@@ -86,20 +86,32 @@ def test_rendering():
         500: (100, 100, 100, 50),
         1000: (255, 255, 255, 255),
     }
-    d = numpy.ma.zeros((1, 256, 256), dtype="float32") + 1
-    d[0, 0:10, 0:10] = 500
-    d[0, 10:20, 10:20] = 1000
+    data = numpy.ma.zeros((1, 256, 256), dtype="float32") + 1
+    data.mask = False
+
+    data[0, 0, 0] = 0
+    data.mask[0, 0, 0] = True
+    data[0, 1:, 1:] = 1
+    data[0, 2:, 2:] = 500
+    data[0, 3:, 3:] = 1000
+
     content, media = render_image(
-        ImageData(d),
+        ImageData(data),
         output_format=ImageType.png,
         colormap=cm,
     )
     assert media == "image/png"
 
     with MemoryFile(content) as mem:
         with mem.open() as dst:
+            data_converted = dst.read()
             assert dst.count == 4
             assert dst.dtypes == ("uint8", "uint8", "uint8", "uint8")
-            assert dst.read()[:, 0, 0].tolist() == [100, 100, 100, 50]
-            assert dst.read()[:, 11, 11].tolist() == [255, 255, 255, 255]
-            assert dst.read()[:, 30, 30].tolist() == [0, 0, 0, 0]
+            # Masked from Original Mask | set to UINT8 (0)
+            assert data_converted[:, 0, 0].tolist() == [0, 0, 0, 0]
+            # masked from CMAP
+            assert data_converted[:, 1, 1].tolist() == [0, 0, 0, 0]
+            # Partially masked from CMAP
+            assert data_converted[:, 2, 2].tolist() == [100, 100, 100, 50]
+            # Non-masked from CMAP
+            assert data_converted[:, 3, 3].tolist() == [255, 255, 255, 255]

src/titiler/core/titiler/core/factory.py

@@ -1301,6 +1301,7 @@ def point(
                 "coordinates": [lon, lat],
                 "values": pts.array.tolist(),
                 "band_names": pts.band_names,
+                "band_descriptions": pts.band_descriptions,
             }
 
     ############################################################################

src/titiler/core/titiler/core/models/responses.py

@@ -18,9 +18,10 @@ class Point(BaseModel):
 
     """
 
-    coordinates: List[float]
-    values: List[Optional[float]]
-    band_names: List[str]
+    coordinates: list[float]
+    values: list[float | None]
+    band_names: list[str]
+    band_descriptions: list[str] | None = None
 
 
 InfoGeoJSON = Feature[Union[Polygon, MultiPolygon], Info]

src/titiler/core/titiler/core/utils.py

@@ -76,13 +76,16 @@ def render_image(  # noqa: C901
         image.apply_color_formula(color_formula)
 
     data, mask = image.data.copy(), image.mask.copy()
-    datatype_range = image.dataset_statistics or (dtype_ranges[str(data.dtype)],)
+    input_range = dtype_ranges[str(data.dtype)]
+    output_range = image.dataset_statistics or (input_range,)
 
     if colormap:
         data, alpha_from_cmap = apply_cmap(data, colormap)
+        output_range = (dtype_ranges[str(data.dtype)],)
         # Combine both Mask from dataset and Alpha band from Colormap
-        mask = numpy.bitwise_and(alpha_from_cmap, mask)
-        datatype_range = (dtype_ranges[str(data.dtype)],)
+        mask = numpy.where(
+            mask != input_range[0], alpha_from_cmap, output_range[0][0]
+        ).astype(data.dtype)
 
     # If output_format is not set, we choose between JPEG and PNG
     if not output_format:
@@ -105,7 +108,7 @@ def render_image(  # noqa: C901
             InvalidDatatypeWarning,
             stacklevel=1,
         )
-        data = rescale_array(data, mask, in_range=datatype_range)
+        data = rescale_array(data, mask, in_range=output_range)
 
     creation_options = {**kwargs, **output_format.profile}
     if output_format.driver == "GTiff":

src/titiler/extensions/pyproject.toml

@@ -34,7 +34,7 @@ dependencies = [
 
 [project.optional-dependencies]
 cogeo = [
-    "rio-cogeo>=5.0,<6.0",
+    "rio-cogeo>=7.0,<8.0",
 ]
 stac = [
     "rio-stac>=0.12,<0.13",