diff --git a/modules/util/viirs_l1b_l2.py b/modules/util/viirs_l1b_l2.py index a1ff00554f22e5c2910b29a7d948571e33f01013..6e9f238be47f3909f838a5003163e2a70c4492ad 100644 --- a/modules/util/viirs_l1b_l2.py +++ b/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)