rio-tiler 6.4.5__tar.gz → 6.4.7__tar.gz

{rio_tiler-6.4.5 → rio_tiler-6.4.7}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: rio-tiler
-Version: 6.4.5
+Version: 6.4.7
 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/

{rio_tiler-6.4.5 → rio_tiler-6.4.7}/rio_tiler/__init__.py

@@ -1,6 +1,6 @@
 """rio-tiler."""
 
-__version__ = "6.4.5"
+__version__ = "6.4.7"
 
 from . import (  # noqa
     colormap,

{rio_tiler-6.4.5 → rio_tiler-6.4.7}/rio_tiler/io/base.py

@@ -48,6 +48,18 @@ class SpatialMixin:
     def geographic_bounds(self) -> BBox:
         """Return dataset bounds in geographic_crs."""
         if self.crs == self.geographic_crs:
+            if self.bounds[1] > self.bounds[3]:
+                warnings.warn(
+                    "BoundingBox of the dataset is inverted (minLat > maxLat).",
+                    UserWarning,
+                )
+                return (
+                    self.bounds[0],
+                    self.bounds[3],
+                    self.bounds[2],
+                    self.bounds[1],
+                )
+
             return self.bounds
 
         try:

{rio_tiler-6.4.5 → rio_tiler-6.4.7}/rio_tiler/io/stac.py

@@ -2,6 +2,7 @@
 
 import json
 import os
+import warnings
 from typing import Any, Dict, Iterator, Optional, Set, Type, Union
 from urllib.parse import urlparse
 
@@ -306,7 +307,16 @@ class STACReader(MultiBaseReader):
                 for b in bands
                 if {"minimum", "maximum"}.issubset(b.get("statistics", {}))
             ]
-            if len(stats) == len(bands):
+            # check that stats data are all double and make warning if not
+            if (
+                stats
+                and all(isinstance(v, (int, float)) for stat in stats for v in stat)
+                and len(stats) == len(bands)
+            ):
                 info["dataset_statistics"] = stats
+            else:
+                warnings.warn(
+                    "Some statistics data in STAC are invalid, they will be ignored."
+                )
 
         return info

{rio_tiler-6.4.5 → rio_tiler-6.4.7}/rio_tiler/io/xarray.py

@@ -401,7 +401,7 @@ class XarrayReader(BaseReader):
 
         y, x = rowcol(ds.rio.transform(), ds_lon, ds_lat)
 
-        arr = ds[:, y[0], x[0]].to_masked_array()
+        arr = ds[:, int(y[0]), int(x[0])].to_masked_array()
         arr.mask |= arr.data == ds.rio.nodata
 
         return PointData(

{rio_tiler-6.4.5 → rio_tiler-6.4.7}/rio_tiler/reader.py

@@ -142,11 +142,17 @@ def read(
                 "crs": dst_crs,
                 "add_alpha": True,
                 "resampling": warp_resampling,
+                "dtype": src_dst.dtypes[0],
             }
 
             if nodata is not None:
                 vrt_params.update(
-                    {"nodata": nodata, "add_alpha": False, "src_nodata": nodata}
+                    {
+                        "nodata": nodata,
+                        "add_alpha": False,
+                        "src_nodata": nodata,
+                        "dtype": src_dst.dtypes[0],
+                    }
                 )
 
             if has_alpha_band(src_dst):
@@ -396,6 +402,7 @@ def part(
             "transform": vrt_transform,
             "width": vrt_width,
             "height": vrt_height,
+            "dtype": src_dst.dtypes[0],
         }
         if vrt_options:
             vrt_params.update(**vrt_options)

{rio_tiler-6.4.5 → rio_tiler-6.4.7}/rio_tiler/utils.py

@@ -334,6 +334,23 @@ def get_vrt_transform(
             src_height = round(w.height)
             src_width = round(w.width)
 
+        # Specific FIX when bounds and transform are inverted
+        # See: https://github.com/US-GHG-Center/veda-config-ghg/pull/333
+        elif (
+            src_dst.crs == WGS84_CRS
+            and dst_crs == WEB_MERCATOR_CRS
+            and (src_bounds[1] > 85.06 or src_bounds[3] < -85.06)
+        ):
+            warnings.warn(
+                "Adjusting dataset latitudes to avoid re-projection overflow",
+                UserWarning,
+            )
+            src_bounds[1] = min(src_bounds[1], 85.06)
+            src_bounds[3] = max(src_bounds[3], -85.06)
+            w = windows.from_bounds(*src_bounds, transform=src_dst.transform)
+            src_height = round(w.height)
+            src_width = round(w.width)
+
         dst_transform, _, _ = calculate_default_transform(
             src_dst.crs, dst_crs, src_width, src_height, *src_bounds
         )