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

snapshot...

parent 5ea0dc1b
No related branches found
No related tags found
No related merge requests found
......@@ -64,10 +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)
keep = np.where(keep, grd_k < 0.50, False)
# 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)
frac_keep = np.sum(keep)/num_keep
if 0.30 < frac_keep < 0.70:
if 0.40 < frac_keep < 0.60:
return None
grd_k = np.where(np.invert(keep), 0, grd_k)
return grd_k
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment