Take coverage weights in account when generating statistics from GeoJSON features

CHANGES.md

@@ -28,6 +28,9 @@
 
   - `/crop/{width}x{height}.{format}  [POST]` -> `/feature/{width}x{height}.{format} [POST]`
 
+- update `rio-tiler` requirement to `>=6.2.1`
+
+- Take coverage weights in account when generating statistics from GeoJSON features
 
 ## 0.14.1 (2023-09-14)
 

src/titiler/core/pyproject.toml

@@ -35,7 +35,7 @@ dependencies = [
     "numpy",
     "pydantic~=2.0",
     "rasterio",
-    "rio-tiler>=6.0,<7.0",
+    "rio-tiler>=6.2.1,<7.0",
     "morecantile>=5.0,<6.0",
     "simplejson",
     "typing_extensions>=4.6.1",

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

@@ -456,16 +456,25 @@ def geojson_statistics(
             with rasterio.Env(**env):
                 with self.reader(src_path, **reader_params) as src_dst:
                     for feature in fc:
+                        shape = feature.model_dump(exclude_none=True)
                         data = src_dst.feature(
-                            feature.model_dump(exclude_none=True),
+                            shape,
                             shape_crs=coord_crs or WGS84_CRS,
                             **layer_params,
                             **image_params,
                             **dataset_params,
                         )
 
+                        # Get the coverage % array
+                        coverage_array = data.get_coverage_array(
+                            shape,
+                            shape_crs=coord_crs or WGS84_CRS,
+                        )
+
                         stats = data.statistics(
-                            **stats_params, hist_options={**histogram_params}
+                            **stats_params,
+                            hist_options={**histogram_params},
+                            coverage=coverage_array,
                         )
 
                         feature.properties = feature.properties or {}