pysteps.verification.spatialscores.intensity_scale

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

[CRS04], [RL08], [EWW+13]