Merge pull request pytroll#2824 from ludwigVonKoopa/rioxarray_imports

add import rioxarray where readers actually need them

pyproject.toml

@@ -41,6 +41,7 @@ avhrr_l1b_gaclac = ["pygac >= 1.3.0"]
 modis_l1b = ["pyhdf", "python-geotiepoints >= 1.1.7"]
 geocat = ["pyhdf"]
 goci2 = ["netCDF4 >= 1.1.8"]
+generic_image = ["rasterio", "rioxarray"]
 acspo = ["netCDF4 >= 1.1.8"]
 clavrx = ["netCDF4 >= 1.1.8"]
 viirs_l1b = ["netCDF4 >= 1.1.8"]

satpy/readers/generic_image.py

@@ -31,6 +31,7 @@
 import dask.array as da
 import numpy as np
 import rasterio
+import rioxarray  # noqa: F401  # need by xarray with the engine rasterio
 import xarray as xr
 from pyresample import utils
 
@@ -76,6 +77,9 @@ def read(self):
         if hasattr(dataset, "crs") and dataset.crs is not None:
             self.area = utils.get_area_def_from_raster(dataset)
 
+        # xarray use the engine 'rasterio' to open the file, but
+        #   its actually rioxarray used in the backend.
+        #   however, error is not explicit enough (see https://github.com/pydata/xarray/issues/7831)
         data = xr.open_dataset(self.finfo["filename"], engine="rasterio",
                                chunks={"band": 1, "y": CHUNK_SIZE, "x": CHUNK_SIZE}, mask_and_scale=False)["band_data"]
         if hasattr(dataset, "nodatavals"):

satpy/readers/sar_c_safe.py

@@ -46,6 +46,7 @@
 import defusedxml.ElementTree as ET
 import numpy as np
 import rasterio
+import rioxarray  # noqa F401  # xarray open_dataset use engine rasterio, which use rioxarray
 import xarray as xr
 from dask import array as da
 from geotiepoints.geointerpolator import lonlat2xyz, xyz2lonlat