signing urls with planetary computers module because anonymous acess …

…to blob storage is not accepted anymore

docs/notebooks/gis.ipynb

@@ -849,6 +849,13 @@
     ")"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "We can also extract the surface properties for the same `gpd.GeoDataFrame` : "
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": 10,
@@ -972,6 +979,13 @@
     "xhgis.surface_properties(gdf)"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Again, for convenience, we can output the results in `xarray.Dataset` format : "
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": 11,

src/xhydro/gis.py

@@ -15,6 +15,7 @@
 import matplotlib.pyplot as plt
 import numpy as np
 import pandas as pd
+import planetary_computer
 import pystac_client
 import rasterio
 import rasterio.features
@@ -270,11 +271,9 @@ def _recursive_upstream_lookup(
 
 
 def _flatten(x, dim="time"):
-    # FIXME: assert statements should only be found in test code
-    assert isinstance(x, xr.DataArray)  # noqa: S101
-    if len(x[dim].values) > len(set(x[dim].values)):
-        x = x.groupby(dim).map(stackstac.mosaic)
-
+    if isinstance(x, xr.DataArray):
+        if len(x[dim].values) > len(set(x[dim].values)):
+            x = x.groupby(dim).map(stackstac.mosaic)
     return x
 
 
@@ -324,6 +323,7 @@ def surface_properties(
 
     catalog = pystac_client.Client.open(
         "https://planetarycomputer.microsoft.com/api/stac/v1",
+        modifier=planetary_computer.sign_inplace,
     )
 
     search = catalog.search(
@@ -489,6 +489,7 @@ def land_use_classification(
     """
     catalog = pystac_client.Client.open(
         "https://planetarycomputer.microsoft.com/api/stac/v1",
+        modifier=planetary_computer.sign_inplace,
     )
     collection = catalog.get_collection(collection)
     ia = ItemAssetsExtension.ext(collection)
@@ -552,6 +553,7 @@ def land_use_plot(
     """
     catalog = pystac_client.Client.open(
         "https://planetarycomputer.microsoft.com/api/stac/v1",
+        modifier=planetary_computer.sign_inplace,
     )
 
     collection = catalog.get_collection(collection)