24.12.2013 Views

wradlib Documentation - Bitbucket

wradlib Documentation - Bitbucket

wradlib Documentation - Bitbucket

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

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

constr_args : list<br />

list of lists, which are to be passed to the individual constraint functions using the *args<br />

mechanism (len(constr_args) == len(constraints))<br />

diagnostics : dictionary<br />

dictionary of variables, which are usually not returned by the function but may be of<br />

interest for research or during testing. Defaults to {}, in which case no diagnostics are<br />

generated. If a dictionary with certain keys is passed to the function, the respective<br />

diagnostics are generated. Currently implemented diagnostics:<br />

• ‘a’ - returns the values of the a coefficient of the k-Z relation, which was used to<br />

calculate the attenuation for the respective beam as a np.array. The shape of the<br />

returned array will be gateset.shape[:-1].<br />

Returns k : array<br />

Array with the same shape as gateset containing the calculated attenuation [dB] for<br />

each range gate.<br />

Raises AttenuationOverflowError :<br />

Exception, if not all constraints are satisfied even with the smallest possible linear coefficient<br />

(a_min) and no handling mode is set.<br />

Examples<br />

>>> # Implementing the original Hitschfeld & Bordan (1954) algorithm with<br />

>>> # otherwise default parameters<br />

>>> k = correctAttenuationConstrained(gateset, n=1, mode=’nan’)<br />

>>> # Implementing the basic Kraemer algorithm<br />

>>> k = correctAttenuationConstrained(gateset,<br />

... mode=’nan’,<br />

... constraints=[constraint_dBZ],<br />

... constr_args=[[59.0]])<br />

>>> # Implementing the PIA algorithm by Jacobi et al.<br />

>>> k = correctAttenuationConstrained(gateset,<br />

... mode=’nan’,<br />

... constraints=[constraint_dBZ,<br />

... constraint_PIA],<br />

... constr_args=[[59.0],<br />

... [20.0]])<br />

3.12.5 <strong>wradlib</strong>.atten.constraint_dBZ<br />

<strong>wradlib</strong>.atten.constraint_dBZ(gateset, k, thrs_dBZ)<br />

Constraint callback function for correctAttenuationConstrained. Selects beams, in which at least one pixel<br />

exceeds thrs_dBZ [dBZ].<br />

3.12.6 <strong>wradlib</strong>.atten.constraint_PIA<br />

<strong>wradlib</strong>.atten.constraint_PIA(gateset, k, thrs_PIA)<br />

Constraint callback function for correctAttenuationConstrained. Selects beams, in which the path integrated<br />

attenuation exceeds thrs_PIA.<br />

3.12. Attenuation Correction 63

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

Saved successfully!

Ooh no, something went wrong!