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

snapshot...

parent 5563614c
No related branches found
No related tags found
No related merge requests found
...@@ -48,12 +48,11 @@ def keep_tile(p_idx, tile): ...@@ -48,12 +48,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.30, False) # keep_clr = np.where(keep, grd_k < 0.30, False)
keep_cld = np.where(keep, grd_k > 0.70, 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.20 and frac_cld >= 0.20): if not (0.10 <= frac_cld <= 0.90):
if not frac_cld >= 0.15:
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