ocetrac.measures.calculate_intensity_metrics

ocetrac.measures.calculate_intensity_metrics(anomalies: DataArray, quantile_threshold: float = 0.9) → dict

Calculate intensity metrics from anomaly data.

Parameters

anomaliesxr.DataArray: Input anomaly data with dimensions (time, lat, lon)
quantile_thresholdfloat, optional: Quantile threshold for calculating extreme intensities (0-1), default 0.9

dict: Dictionary containing: - ‘cumulative_intensity’: Spatial sum time series - ‘mean_intensity_timeseries’: Spatial mean time series - ‘mean_intensity’: Global mean value - ‘max_intensity_timeseries’: Spatial max time series - ‘max_intensity’: Global maximum value - ‘std_intensity_timeseries’: Spatial std time series - ‘std_intensity’: Global std value - f’percentile_{int(quantile_threshold*100)}_intensity_timeseries’: Spatial quantile time series - f’percentile_{int(quantile_threshold*100)}_intensity’: Global quantile value - ‘quantile_threshold_used’: The quantile threshold applied