From bd96427aecad7b4f9085e0edb31e7c0b0d2af321 Mon Sep 17 00:00:00 2001 From: tomrink <rink@ssec.wisc.edu> Date: Sat, 18 Feb 2023 14:07:00 -0600 Subject: [PATCH] snapshot... --- modules/util/viirs_l1b_l2.py | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) diff --git a/modules/util/viirs_l1b_l2.py b/modules/util/viirs_l1b_l2.py index a1ff0055..6e9f238b 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) -- GitLab