snapshot...

modules/util/viirs_l1b_l2.py

@@ -64,12 +64,11 @@ def process_cld_prob_(grd_k, dum):
     if num_keep / grd_k.size < 0.98:
         return None
     # hist_10 += np.histogram(grd_k.flatten(), range=[0.0, 1.0], bins=10)[0]
-    # keep = np.where(keep, grd_k < 0.50, False)
-    keep = np.where(keep, np.logical_or(0.15 > grd_k, grd_k < 0.85), False)
+    keep = np.where(keep, grd_k < 0.20, False)
     frac_keep = np.sum(keep)/num_keep
-    if 0.40 < frac_keep < 0.60:
+    if not (0.40 < frac_keep < 0.60):
         return None
-    grd_k = np.where(np.invert(keep), 0, grd_k)
+    grd_k = np.where(np.invert(keep), 0, grd_k)  # Convert NaN to 0
     return grd_k
 
 
@@ -141,7 +140,7 @@ def run_all(directory, out_directory, day_night='ANY', start=10):
             if len(data_train_tiles) == 0 and len(data_valid_tiles) == 0:
                 continue
 
-            if (f_cnt % 5) == 0:
+            if (f_cnt % 10) == 0:
                 num_valid_samples = 0
                 if len(data_valid_tiles) > 0:
                     data_valid = np.stack(data_valid_tiles)