pysteps.verification.spatialscores.binary_mse

pysteps.verification.spatialscores.binary_mse(X_f, X_o, thr, wavelet='haar')

Compute an intensity-scale verification as the MSE of the binary error.

This method uses PyWavelets for decomposing the error field between the forecasts and observations into multiple spatial scales.

Parameters:

X_f : array_like

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

X_o : array_like

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

thr : sequence

The intensity threshold for which to compute the verification.

wavelet : str, optional

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

Returns:

SS : array

One-dimensional array containing the binary MSE for each spatial scale.

spatial_scale : list

References

[CRS04]