pysteps.verification.spatialscores.binary_mse

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

Compute the MSE of the binary error as a function of spatial scale.

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 (m, n) containing the forecast field.

X_o : array_like

Array of shape (m, n) 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.

return_scales : bool, optional

Whether to return the spatial scales resulting from the wavelet decomposition.

Returns:

SS : array

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

scales : list, optional

If *return_scales*=True, return the spatial scales in pixels resulting from the wavelet decomposition.

References

[CRS2004]