24.12.2013 Views

wradlib Documentation - Bitbucket

wradlib Documentation - Bitbucket

wradlib Documentation - Bitbucket

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>wradlib</strong> <strong>Documentation</strong>, Release 0.1.1<br />

All raw data is provided by DWD<br />

2.4 Converting reflectivities<br />

Reflectivity (Z) and precipitation rate (R) can be related in form of a power law R=a*Z**b. The parameters a and b<br />

depend on the type of precipitation (i.e. drop size distribution and water temperature).<br />

Because the ZR-relationship is based on reflectivities in Z and the encoded radar data is given in dBZ we have to take<br />

the antilogarithm of the dBZ-values after reading the DX-data (Read and plot raw DWD radar data):<br />

import <strong>wradlib</strong> as wrl<br />

# set the location of your data<br />

datadir = ’D:/THIS/IS/MY/DATA/DIRECTORY/’<br />

data_dBZ = wrl.io.readDX(datadir + ’raa00-dx_10908-200608180225-fbg---bin’)[0]<br />

data_Z = wrl.trafo.idecibel(data_dBZ)<br />

After taking the antilogarithm the scale of low precipitation rates appears distinctly more squeezed:<br />

wrl.vis.polar_plot(data_dBZ, title = ’Reflectivity: Radarscan Feldberg 18-08-2006 02:25’,<br />

unit = ’dBZ’, colormap = ’spectral’)<br />

wrl.vis.polar_plot(data_Z, title = ’Reflectivity: Radarscan Feldberg 18-08-2006 02:25’,<br />

unit = ’Z’, colormap = ’spectral’)<br />

22 Chapter 2. Tutorials

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!