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 />

References<br />

[Jacobi2011]<br />

3.12.4 <strong>wradlib</strong>.atten.correctAttenuationConstrained<br />

<strong>wradlib</strong>.atten.correctAttenuationConstrained(gateset, a_max=0.000167, a_min=2.33e-05,<br />

b_max=0.7, b_min=0.2, na=30, nb=5,<br />

l=1.0, mode=’error’, constraints=None,<br />

constr_args=None, diagnostics={})<br />

Gate-by-Gate attenuation correction based on the iterative approach of Stefan Kraemer [Kraemer2008] with a<br />

generalized and arbitrary number of constraints.<br />

Parameters gateset : array<br />

Multidimensional array, where the range gates (over which iteration has to be performed)<br />

are supposed to vary along the last dimension so, e.g., for a set of l radar<br />

images stored in polar form with m azimuths and n range-bins the input array’s shape<br />

can be either (l,m,n) or (m,l,n).<br />

Data has to be provided in decibel representation of reflectivity [dBZ].<br />

a_max : float<br />

initial value for linear coefficient of the k-Z relation ( k = a ∗ Z b ). Per default set to<br />

1.67e-4.<br />

a_min : float<br />

minimal allowed linear coefficient of the k-Z relation ( k = a ∗ Z b ) in the downwards<br />

iteration of a in case of signal overflow (sum of signal and attenuation exceeds the<br />

threshold thrs). Per default set to 2.33e-5.<br />

b : float<br />

exponent of the k-Z relation ( k = a ∗ Z b ). Per default set to 0.7.<br />

n : integer<br />

l : float<br />

number of iterations from a_max to a_min. Per default set to 30.<br />

length of a range gate [km]. Per default set to 1.0.<br />

mode : string<br />

Controls how the function reacts in case of signal overflow (sum of signal and attenuation<br />

exceeds the threshold thrs). Possible values:<br />

‘warn’ : emit a warning through the module’s logger but continue execution<br />

‘zero’ : set offending gates to 0.0<br />

‘nan’ : set offending gates to nan<br />

Per default set to ‘zero’. Any other mode will raise an Exception.<br />

constraints : list<br />

list of constraint functions. The signature of these functions has to be constraint_function(gateset,<br />

k, *‘constr_args‘). Their return value must be a boolean array<br />

of shape gateset.shape[:-1] set to True for beams, which do not fulfill the constraint.<br />

62 Chapter 3. Library Reference

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

Saved successfully!

Ooh no, something went wrong!