fill value

worldbank · Dec 20, 2023 · 938a3fc · 938a3fc
1 parent 344f9c1
commit 938a3fc
Show file tree

Hide file tree

Showing 2 changed files with 1,417 additions and 543 deletions.
diff --git a/examples/quality-assessment.ipynb b/examples/quality-assessment.ipynb
diff --git a/src/blackmarble/raster.py b/src/blackmarble/raster.py
@@ -117,11 +117,13 @@ def h5_to_geotiff(
         data = dataset[:]
 
         data = np.where(data == fill_value, np.nan, data)
-        data = scale_factor * dataset[:] + offset
-        qf = qf[:]
+        data = scale_factor * data + offset
 
-        for val in quality_flag_rm:
-            data = np.where(qf == val, np.nan, data)
+        if len(quality_flag_rm) > 0:
+            qf = qf[:]
+
+            for val in quality_flag_rm:
+                data = np.where(qf == val, np.nan, data)
 
         # Get geospatial metadata (coordinates and attributes)
         height, width = data.shape