La mareggiata vista dagli ormeggi

Funzione Python

import numpy as np
import pandas as pd

def ampiezze(x, soglia):
    """Ampiezze dei semicicli tra due inversioni successive della tensione."""
    d = np.diff(x)
    inversioni = x[np.r_[True, d[1:] * d[:-1] < 0, True]]
    a = np.abs(np.diff(inversioni))
    return a[a >= soglia]

def transform(inputs, params, ctx):
    d = inputs["rows"].sort_values("ts")
    d["inizio"] = pd.to_datetime(d["ts"]).dt.floor(f"{params['finestra_minuti']}min")
    righe = []
    for inizio, w in d.groupby("inizio"):
        a = ampiezze(w["tensione_l3_kn"].to_numpy(), params["soglia_ciclo_kn"])
        righe.append({
            "inizio": inizio,
            "hs_m": 4 * w["onda_m"].std(ddof=0),
            "rollio_rms_deg": float(np.sqrt(np.mean(w["rollio_deg"] ** 2))),
            "tensione_max_l3_kn": w["tensione_l3_kn"].max(),
            "danno_l3": 0.5 * float(np.sum(a ** params["esponente_sn"])),
        })
    out = pd.DataFrame(righe)
    out["danno_l3_cumulato"] = out["danno_l3"].cumsum()
    return {"rows": out}