From f97351f5513cbaf84c952b4560607679ce361866 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 16 Mar 2023 11:29:01 -0500
Subject: [PATCH] snapshot...
---
modules/util/viirs_l1b_l2.py | 13 +++++++++++++
1 file changed, 13 insertions(+)
diff --git a/modules/util/viirs_l1b_l2.py b/modules/util/viirs_l1b_l2.py
index 90bea111..203fd843 100644
--- a/modules/util/viirs_l1b_l2.py
+++ b/modules/util/viirs_l1b_l2.py
@@ -85,6 +85,19 @@ def process_cld_prob(grd_k):
return grd_k
+def process_cld_prob_2(grd_k):
+ keep = np.invert(np.isnan(grd_k))
+ num_keep = np.sum(keep)
+ keep_clr = np.where(keep, grd_k < 0.20, False)
+ keep_cld = np.where(keep, grd_k > 0.80, False)
+ frac_clr = np.sum(keep_clr)/num_keep
+ frac_cld = np.sum(keep_cld)/num_keep
+ if not (frac_clr >= 0.35 and frac_cld >= 0.35):
+ return None
+ grd_k = np.where(np.invert(keep), 0, grd_k) # Convert NaN to 0
+ return grd_k
+
+
def process_cld_opd(grd_k):
keep = np.invert(np.isnan(grd_k))
num_keep = np.sum(keep)
--
GitLab