Apply histo-cut filter to retrieve boolean array that highlights clutter as well as beam blockage¶

fig=pl.figure(figsize=(14,8))ax=fig.add_subplot(121)ax,pm=plot_ppi(np.log(yearsum),ax=ax)pl.title("Logarithm of annual precipitation sum")pl.colorbar(pm,shrink=0.75)ax=fig.add_subplot(122)ax,pm=plot_ppi(mask.astype(np.uint8),ax=ax)pl.title("Map of execptionally low and high values\n(clutter and beam blockage)")pl.colorbar(pm,shrink=0.75)