pysteps.postprocessing.ensemblestats.excprob

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.