From 5ea0dc1bf7c78c9cbc112a45e39cbf20a6af60a6 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sat, 18 Feb 2023 12:51:27 -0600
Subject: [PATCH] snapshot...
---
modules/util/viirs_l1b_l2.py | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/modules/util/viirs_l1b_l2.py b/modules/util/viirs_l1b_l2.py
index c0e24df4..f62e5c5c 100644
--- a/modules/util/viirs_l1b_l2.py
+++ b/modules/util/viirs_l1b_l2.py
@@ -64,8 +64,10 @@ 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, np.logical_and(0.10 < grd_k, grd_k < 0.90), False)
- if np.sum(keep)/num_keep < 0.30:
+ # keep = np.where(keep, np.logical_and(0.10 < grd_k, grd_k < 0.90), False)
+ keep = np.where(keep, grd_k < 0.50, False)
+ frac_keep = np.sum(keep)/num_keep
+ if 0.30 < frac_keep < 0.70:
return None
grd_k = np.where(np.invert(keep), 0, grd_k)
return grd_k
@@ -89,7 +91,8 @@ def process_cld_opd_(grd_k, dum):
return None
grd_k = np.where(np.invert(keep), 0, grd_k)
keep = np.where(keep, np.logical_and(0.1 < grd_k, grd_k < 158.0), False)
- if np.sum(keep)/num_keep < 0.50:
+ frac_keep = np.sum(keep)/num_keep
+ if frac_keep < 0.50:
return None
return grd_k
--
GitLab