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

Raw unadjusted radar rainfall<br />

targets : (INTERNAL) array of floats<br />

Coordinate pairs for locations on which the final adjustment product is interpolated Defaults<br />

to None. In this case, the output locations will be identical to the radar coordinates<br />

rawatobs : (INTERNAL) array of floats<br />

For internal use from AdjustBase.xvalidate only (defaults to None)<br />

ix : (INTERNAL) array of integers<br />

For internal use from AdjustBase.xvalidate only (defaults to None)<br />

<strong>wradlib</strong>.adjust.AdjustMFB.xvalidate<br />

AdjustMFB.xvalidate(obs, raw)<br />

Leave-One-Out Cross Validation, applicable to all gage adjustment classes.<br />

This method will be inherited to other Adjust classes. It should thus be applicable to all adjustment procedures<br />

without any modification. This way, the actual adjustment procedure has only to be defined once in the __call__<br />

method.<br />

The output of this method can be evaluated by using the verify.ErrorMetrics class.<br />

Parameters obs : array of floats<br />

raw : array of floats<br />

Returns obs : array of floats<br />

valid observations at those locations which have a valid radar observation<br />

estatobs : array of floats<br />

estimated values at the valid observation locations<br />

<strong>wradlib</strong>.adjust.AdjustMultiply<br />

class <strong>wradlib</strong>.adjust.AdjustMultiply(obs_coords, raw_coords, nnear_raws=9, stat=’median’,<br />

mingages=5, minval=0.0, Ipclass=, **ipargs)<br />

Gage adjustment using a multiplicative error model<br />

First, an instance of AdjustMultiply has to be created. Calling this instance then does the actual adjustment.<br />

The motivation behind this performance. In case the observation points are always the same for different time<br />

steps, the computation of neighbours and invserse distance weights only needs to be performed once during<br />

initialisation.<br />

AdjustMultiply automatically takes care of invalid gage or radar observations (e.g. NaN, Inf or other typical<br />

missing data flags such as -9999. However, in case e.g. the observation data contain missing values, the computation<br />

of the inverse distance weights needs to be repeated in __call__ which is at the expense of performance.<br />

Parameters obs_coords : array of float<br />

coordinate pairs of observations points<br />

raw_coords : array of float<br />

coordinate pairs of raw (unadjusted) field<br />

nnear_raws : integer<br />

3.13. Gage adjustment 69

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

Saved successfully!

Ooh no, something went wrong!