pysteps.verification.spatialscores.intensity_scale

pysteps.verification.spatialscores.intensity_scale(X_f, X_o, name, thrs, scales=None, wavelet='Haar')

Compute an intensity-scale verification score.

Parameters:
X_f : array_like

Array of shape (m, n) containing the forecast field.

X_o : array_like

Array of shape (m, n) containing the verification observation field.

name : string

A string indicating the name of the spatial verification score to be used:

Name Description
FSS Fractions skill score
BMSE Binary mean squared error
thrs : float or array_like

Scalar or 1-D array of intensity thresholds for which to compute the verification.

scales : float or array_like, optional

Scalar or 1-D array of spatial scales in pixels, required if *name*=”FSS”.

wavelet : str, optional

The name of the wavelet function to use in the BMSE. Defaults to the Haar wavelet, as described in Casati et al. 2004. See the documentation of PyWavelets for a list of available options.

Returns:
out : array_like

The two-dimensional array containing the intensity-scale skill scores for each spatial scale and intensity threshold.

References

[CRS2004], [RL2008], [EWWM2013]