Merge pull request #2923 from pnuu/bugfix-generic_image-float64

Fix data promotion in `generic_image` reader and `satpy.composites.add_bands`
pytroll · Oct 14, 2024 · b02fd62 · b02fd62
2 parents 2d7f5ed + b8991c0
commit b02fd62
Show file tree

Hide file tree

Showing 4 changed files with 304 additions and 264 deletions.
diff --git a/satpy/composites/__init__.py b/satpy/composites/__init__.py
@@ -985,6 +985,7 @@ def add_bands(data, bands):
         alpha = new_data[0].copy()
         alpha.data = da.ones((data.sizes["y"],
                               data.sizes["x"]),
+                             dtype=new_data[0].dtype,
                              chunks=new_data[0].chunks)
         # Rename band to indicate it's alpha
         alpha["bands"] = "A"

diff --git a/satpy/readers/generic_image.py b/satpy/readers/generic_image.py
@@ -147,7 +147,7 @@ def _mask_image_data(data, info):
         if not np.issubdtype(data.dtype, np.integer):
             raise ValueError("Only integer datatypes can be used as a mask.")
         mask = data.data[-1, :, :] == np.iinfo(data.dtype).min
-        data = data.astype(np.float64)
+        data = data.astype(np.float32)
         masked_data = da.stack([da.where(mask, np.nan, data.data[i, :, :])
                                 for i in range(data.shape[0])])
         data.data = masked_data