pysteps.postprocessing.ensemblestats.excprob

pysteps.postprocessing.ensemblestats.excprob(X, X_thr, ignore_nan=False)

For a given forecast ensemble field, compute exceedance probabilities for the given intensity thresholds.

Parameters:

X : array_like

Array of shape (k,m,n,…) containing an k-member ensemble of forecasts with shape (m,n,…).

X_thr : float or a sequence of floats

Intensity threshold(s) for which the exceedance probabilities are computed.

ignore_nan : bool

If True, ignore nan values.

Returns:

out : ndarray

Array of shape (len(X_thr),m,n) containing the exceedance probabilities for the given intensity thresholds. If len(X_thr)=1, the first dimension is dropped.