But we need to into account that the radar signal travels along a beam.
Cumulative beam blockage (CBB) in one bin along a beam will always be at
least as high as the maximum PBB of the preceeding bins (see
wradlib.qual.cum_beam_block_frac)

Now we visualize - the average terrain altitude per radar bin - a beam
blockage map - interaction with terrain along a single beam

In [10]:

# just a little helper function to style x and y axes of our mapsdefannotate_map(ax,cm=None,title=""):ticks=(ax.get_xticks()/1000).astype(np.int)ax.set_xticklabels(ticks)ticks=(ax.get_yticks()/1000).astype(np.int)ax.set_yticklabels(ticks)ax.set_xlabel("Kilometers")ax.set_ylabel("Kilometers")ifnotcmisNone:pl.colorbar(cm,ax=ax)ifnottitle=="":ax.set_title(title)ax.grid()