Skip to content
Snippets Groups Projects
Commit f1170407 authored by tomrink's avatar tomrink
Browse files

snapshot...

parent dd11d98b
Branches
No related tags found
No related merge requests found
...@@ -88,11 +88,11 @@ def keep_tile(p_idx, tile): ...@@ -88,11 +88,11 @@ def keep_tile(p_idx, tile):
def process_cld_prob(grd_k): def process_cld_prob(grd_k):
keep = np.invert(np.isnan(grd_k)) keep = np.invert(np.isnan(grd_k))
num_keep = np.sum(keep) num_keep = np.sum(keep)
keep_clr = np.where(keep, grd_k < 0.20, False) keep_clr = np.where(keep, grd_k < 0.30, False)
keep_cld = np.where(keep, grd_k > 0.80, False) keep_cld = np.where(keep, grd_k > 0.70, False)
frac_clr = np.sum(keep_clr)/num_keep frac_clr = np.sum(keep_clr)/num_keep
frac_cld = np.sum(keep_cld)/num_keep frac_cld = np.sum(keep_cld)/num_keep
if not (frac_clr >= 0.23 and frac_cld >= 0.23): if not (frac_clr >= 0.22 and frac_cld >= 0.22):
return None return None
grd_k = np.where(np.invert(keep), 0, grd_k) # Convert NaN to 0 grd_k = np.where(np.invert(keep), 0, grd_k) # Convert NaN to 0
return grd_k return grd_k
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment