Delete app.py

app.py

@@ -1,456 +0,0 @@
-import numpy as np
-import pandas as pd
-import streamlit as st
-import matplotlib.pyplot as plt
-
-from zebris_extractor import extract_zebris_csv
-
-st.set_page_config(page_title="Zebris — Profil & Seuils", layout="wide")
-
-st.title("Zebris — Profil biomécanique & seuils individualisés")
-st.caption("Import d’un CSV Zebris brut → fiche profil + seuils personnalisés")
-
-with st.sidebar:
-    st.header("Import")
-    uploaded_file = st.file_uploader(
-        "Importer un fichier CSV Zebris",
-        type=["csv"],
-        accept_multiple_files=False,
-    )
-
-    if uploaded_file is not None and uploaded_file.size > 20 * 1024 * 1024:
-        st.error("Fichier trop volumineux (>20 MB)")
-        st.stop()
-
-    st.header("Contexte athlète")
-    volume_horaire = st.number_input(
-        "Volume horaire / semaine",
-        min_value=0.5,
-        max_value=40.0,
-        value=5.0,
-        step=0.5,
-    )
-
-if uploaded_file is None:
-    st.info("Importe un CSV Zebris brut pour afficher la fiche profil et les seuils.")
-    st.stop()
-
-try:
-    df_std, debug = extract_zebris_csv(uploaded_file)
-except Exception as e:
-    st.error(f"Erreur pendant l’extraction du CSV Zebris : {e}")
-    st.stop()
-
-if df_std.empty:
-    st.error("Aucune ligne exploitable n’a été trouvée dans le CSV.")
-    st.stop()
-
-athletes = sorted(df_std["Nom"].dropna().unique().tolist())
-selected_athlete = st.selectbox("Sportif", athletes)
-sub_df = df_std[df_std["Nom"] == selected_athlete].copy()
-
-if sub_df.empty:
-    st.error("Aucune donnée trouvée pour ce sportif.")
-    st.stop()
-
-sub_df = sub_df.sort_values("Vitesse (km/h)")
-speeds = sub_df["Vitesse (km/h)"].dropna().tolist()
-selected_speed = st.selectbox("Allure analysée (km/h)", speeds)
-row = sub_df[sub_df["Vitesse (km/h)"] == selected_speed].iloc[0]
-
-poids_csv = row["Poids (kg)"] if pd.notna(row["Poids (kg)"]) else np.nan
-poids_kg = st.number_input(
-    "Poids du sportif (kg)",
-    min_value=30.0,
-    max_value=150.0,
-    value=float(poids_csv) if pd.notna(poids_csv) else 70.0,
-    step=0.1,
-)
-poids_n = poids_kg * 9.81
-
-
-def avg(a, b):
-    if pd.isna(a) and pd.isna(b):
-        return np.nan
-    if pd.isna(a):
-        return float(b)
-    if pd.isna(b):
-        return float(a)
-    return (float(a) + float(b)) / 2
-
-
-def asym(a, b):
-    m = avg(a, b)
-    if pd.isna(m) or m == 0 or pd.isna(a) or pd.isna(b):
-        return np.nan
-    return abs(float(a) - float(b)) / m * 100
-
-
-def clamp_score(value, low, high, reverse=False):
-    if pd.isna(value):
-        return np.nan
-    score = (value - low) / (high - low) * 100
-    score = max(0, min(100, score))
-    return 100 - score if reverse else score
-
-
-def safe_mean(values):
-    vals = [v for v in values if pd.notna(v)]
-    if not vals:
-        return np.nan
-    return float(np.mean(vals))
-
-
-def compute_profile_metrics(row, poids_kg):
-    poids_n = poids_kg * 9.81
-
-    force_talon_moy = avg(row["Force talon G (N)"], row["Force talon D (N)"])
-    force_avant_moy = avg(row["Force avant-pied G (N)"], row["Force avant-pied D (N)"])
-    pression_talon_moy = avg(row["Pression talon G (N/cm²)"], row["Pression talon D (N/cm²)"])
-    cop_moy = avg(row["COP G (mm)"], row["COP D (mm)"])
-    transition_moy = avg(row["Transition G (s)"], row["Transition D (s)"])
-
-    asym_talon = asym(row["Force talon G (N)"], row["Force talon D (N)"])
-    asym_avant = asym(row["Force avant-pied G (N)"], row["Force avant-pied D (N)"])
-    asym_cop = asym(row["COP G (mm)"], row["COP D (mm)"])
-
-    diff_rotation = (
-        abs(float(row["Rotation G (°)"]) - float(row["Rotation D (°)"]))
-        if pd.notna(row["Rotation G (°)"]) and pd.notna(row["Rotation D (°)"])
-        else np.nan
-    )
-
-    force_talon_bw = force_talon_moy / poids_n if pd.notna(force_talon_moy) and poids_n else np.nan
-    ratio_talon_avant = (
-        force_talon_moy / force_avant_moy
-        if pd.notna(force_talon_moy) and pd.notna(force_avant_moy) and force_avant_moy != 0
-        else np.nan
-    )
-
-    # Impact corrigé : fallback sur force brute si BW indisponible
-    impact_force_score = (
-        clamp_score(force_talon_bw, 0.6, 1.1)
-        if pd.notna(force_talon_bw)
-        else clamp_score(force_talon_moy, 400, 750)
-    )
-    impact_pressure_score = clamp_score(pression_talon_moy, 15, 30)
-
-    impact = safe_mean([
-        0.6 * impact_force_score if pd.notna(impact_force_score) else np.nan,
-        0.4 * impact_pressure_score if pd.notna(impact_pressure_score) else np.nan,
-    ])
-    impact = round(impact) if pd.notna(impact) else np.nan
-
-    dynamique = safe_mean([
-        0.6 * clamp_score(row["Cadence (pas/min)"], 150, 185),
-        0.4 * clamp_score(row["Contact (%)"], 68, 76, reverse=True),
-    ])
-    dynamique = round(dynamique) if pd.notna(dynamique) else np.nan
-
-    sym_components = [x for x in [asym_talon, asym_avant, asym_cop, diff_rotation] if pd.notna(x)]
-    symetrie = round(100 - min(100, np.mean(sym_components) * 2.5)) if sym_components else np.nan
-
-    deroule = safe_mean([
-        0.5 * clamp_score(cop_moy, 210, 260),
-        0.5 * clamp_score(transition_moy, 0.05, 0.09, reverse=True),
-    ])
-    deroule = round(deroule) if pd.notna(deroule) else np.nan
-
-    attaque = "mixte"
-    if pd.notna(ratio_talon_avant):
-        if ratio_talon_avant > 1.05:
-            attaque = "talon"
-        elif ratio_talon_avant < 0.95:
-            attaque = "avant-pied"
-
-    return {
-        "force_talon_moy": force_talon_moy,
-        "force_avant_moy": force_avant_moy,
-        "pression_talon_moy": pression_talon_moy,
-        "cop_moy": cop_moy,
-        "transition_moy": transition_moy,
-        "asym_talon": asym_talon,
-        "asym_avant": asym_avant,
-        "asym_cop": asym_cop,
-        "diff_rotation": diff_rotation,
-        "force_talon_bw": force_talon_bw,
-        "ratio_talon_avant": ratio_talon_avant,
-        "impact": impact,
-        "dynamique": dynamique,
-        "symetrie": symetrie,
-        "deroule": deroule,
-        "attaque": attaque,
-    }
-
-
-def build_summary(row, metrics):
-    impact_txt = "marqué" if pd.notna(metrics["impact"]) and metrics["impact"] >= 70 else "modéré" if pd.notna(metrics["impact"]) and metrics["impact"] >= 45 else "faible"
-    dyn_txt = "bonne" if pd.notna(metrics["dynamique"]) and metrics["dynamique"] >= 70 else "moyenne" if pd.notna(metrics["dynamique"]) and metrics["dynamique"] >= 45 else "faible"
-    sym_txt = "satisfaisante" if pd.notna(metrics["symetrie"]) and metrics["symetrie"] >= 70 else "perfectible"
-    der_txt = "favorable" if pd.notna(metrics["deroule"]) and metrics["deroule"] >= 70 else "intermédiaire" if pd.notna(metrics["deroule"]) and metrics["deroule"] >= 45 else "à surveiller"
-
-    return (
-        f"À {row['Vitesse (km/h)']} km/h, {row['Nom']} présente une attaque {metrics['attaque']}, "
-        f"un impact {impact_txt}, une dynamique {dyn_txt}, une symétrie {sym_txt} "
-        f"et un déroulé {der_txt}."
-    )
-
-
-def draw_radar(metrics):
-    labels = ["Impact", "Dynamique", "Symétrie", "Déroulé"]
-    values = [
-        metrics["impact"] if pd.notna(metrics["impact"]) else 0,
-        metrics["dynamique"] if pd.notna(metrics["dynamique"]) else 0,
-        metrics["symetrie"] if pd.notna(metrics["symetrie"]) else 0,
-        metrics["deroule"] if pd.notna(metrics["deroule"]) else 0,
-    ]
-    values += values[:1]
-    angles = np.linspace(0, 2 * np.pi, len(labels), endpoint=False).tolist()
-    angles += angles[:1]
-
-    fig = plt.figure(figsize=(5, 5))
-    ax = plt.subplot(111, polar=True)
-    ax.plot(angles, values, linewidth=2)
-    ax.fill(angles, values, alpha=0.25)
-    ax.set_xticks(angles[:-1])
-    ax.set_xticklabels(labels)
-    ax.set_ylim(0, 100)
-    ax.set_yticks([25, 50, 75, 100])
-    ax.set_title("Radar biomécanique", pad=20)
-    return fig
-
-
-def draw_evolution(df, poids_kg):
-    data = []
-    for _, r in df.sort_values("Vitesse (km/h)").iterrows():
-        m = compute_profile_metrics(r, poids_kg)
-        data.append({
-            "Vitesse": r["Vitesse (km/h)"],
-            "Impact": m["impact"],
-            "Dynamique": m["dynamique"],
-            "Symétrie": m["symetrie"],
-            "Déroulé": m["deroule"],
-        })
-
-    evo = pd.DataFrame(data)
-    fig, ax = plt.subplots(figsize=(8, 4))
-    for col in ["Impact", "Dynamique", "Symétrie", "Déroulé"]:
-        ax.plot(evo["Vitesse"], evo[col], marker="o", label=col)
-    ax.set_ylim(0, 100)
-    ax.set_xlabel("Vitesse (km/h)")
-    ax.set_ylabel("Score /100")
-    ax.set_title("Évolution avec l’allure")
-    ax.legend()
-    ax.grid(True, alpha=0.3)
-    return fig
-
-
-def compute_thresholds(poids_kg, volume_horaire):
-    poids_n = poids_kg * 9.81
-
-    if volume_horaire <= 3:
-        charge = "faible"
-        force_bw_low, force_bw_high = 0.75, 1.00
-        pression_low, pression_high = 20, 26
-        cadence_low, cadence_high = 160, 172
-        contact_low, contact_high = 70, 74
-        flight_low, flight_high = 26, 30
-        asym_low, asym_high = 6, 10
-        rotation_low, rotation_high = 6, 10
-    elif volume_horaire <= 6:
-        charge = "modérée"
-        force_bw_low, force_bw_high = 0.70, 0.95
-        pression_low, pression_high = 20, 25
-        cadence_low, cadence_high = 165, 176
-        contact_low, contact_high = 69, 73
-        flight_low, flight_high = 27, 31
-        asym_low, asym_high = 5, 9
-        rotation_low, rotation_high = 5, 9
-    else:
-        charge = "élevée"
-        force_bw_low, force_bw_high = 0.65, 0.90
-        pression_low, pression_high = 20, 24
-        cadence_low, cadence_high = 168, 180
-        contact_low, contact_high = 68, 72
-        flight_low, flight_high = 28, 32
-        asym_low, asym_high = 4, 8
-        rotation_low, rotation_high = 4, 8
-
-    return {
-        "charge": charge,
-        "poids_n": poids_n,
-        "force_bw_low": force_bw_low,
-        "force_bw_high": force_bw_high,
-        "force_n_low": force_bw_low * poids_n,
-        "force_n_high": force_bw_high * poids_n,
-        "pression_low": pression_low,
-        "pression_high": pression_high,
-        "cadence_low": cadence_low,
-        "cadence_high": cadence_high,
-        "contact_low": contact_low,
-        "contact_high": contact_high,
-        "flight_low": flight_low,
-        "flight_high": flight_high,
-        "asym_low": asym_low,
-        "asym_high": asym_high,
-        "rotation_low": rotation_low,
-        "rotation_high": rotation_high,
-    }
-
-
-metrics = compute_profile_metrics(row, poids_kg)
-summary = build_summary(row, metrics)
-thresholds = compute_thresholds(poids_kg, volume_horaire)
-
-tab_profil, tab_seuils = st.tabs(["Profil biomécanique", "Seuils individualisés"])
-
-with tab_profil:
-    c1, c2, c3, c4 = st.columns(4)
-    with c1:
-        st.metric("Impact", f"{metrics['impact']}/100" if pd.notna(metrics["impact"]) else "N/A")
-    with c2:
-        st.metric("Dynamique", f"{metrics['dynamique']}/100" if pd.notna(metrics["dynamique"]) else "N/A")
-    with c3:
-        st.metric("Symétrie", f"{metrics['symetrie']}/100" if pd.notna(metrics["symetrie"]) else "N/A")
-    with c4:
-        st.metric("Déroulé", f"{metrics['deroule']}/100" if pd.notna(metrics["deroule"]) else "N/A")
-
-    left, right = st.columns([1.2, 1])
-
-    with left:
-        st.subheader("Carte d’identité biomécanique")
-        st.write(summary)
-
-        indicators = pd.DataFrame(
-            {
-                "Indicateur": [
-                    "Poids",
-                    "Cadence",
-                    "Contact",
-                    "Flight",
-                    "Force talon moyenne",
-                    "Force talon normalisée",
-                    "Pression talon moyenne",
-                    "Asymétrie talon",
-                    "COP moyen",
-                    "Différence rotation",
-                    "Attaque",
-                ],
-                "Valeur": [
-                    f"{poids_kg:.1f} kg",
-                    f"{row['Cadence (pas/min)']:.1f} pas/min" if pd.notna(row["Cadence (pas/min)"]) else "N/A",
-                    f"{row['Contact (%)']:.1f} %" if pd.notna(row["Contact (%)"]) else "N/A",
-                    f"{row['Flight (%)']:.1f} %" if pd.notna(row["Flight (%)"]) else "N/A",
-                    f"{metrics['force_talon_moy']:.1f} N" if pd.notna(metrics["force_talon_moy"]) else "N/A",
-                    f"{metrics['force_talon_bw']:.2f} BW" if pd.notna(metrics["force_talon_bw"]) else "N/A",
-                    f"{metrics['pression_talon_moy']:.1f} N/cm²" if pd.notna(metrics["pression_talon_moy"]) else "N/A",
-                    f"{metrics['asym_talon']:.1f} %" if pd.notna(metrics["asym_talon"]) else "N/A",
-                    f"{metrics['cop_moy']:.1f} mm" if pd.notna(metrics["cop_moy"]) else "N/A",
-                    f"{metrics['diff_rotation']:.1f}°" if pd.notna(metrics["diff_rotation"]) else "N/A",
-                    metrics["attaque"],
-                ],
-            }
-        )
-        st.dataframe(indicators, hide_index=True, use_container_width=True)
-
-    with right:
-        st.subheader("Radar biomécanique")
-        st.pyplot(draw_radar(metrics), use_container_width=True)
-
-    st.subheader("Évolution avec l’allure")
-    st.pyplot(draw_evolution(sub_df, poids_kg), use_container_width=True)
-
-with tab_seuils:
-    r1, r2, r3 = st.columns(3)
-    with r1:
-        st.metric("Poids", f"{poids_kg:.1f} kg")
-    with r2:
-        st.metric("Poids en Newton", f"{thresholds['poids_n']:.1f} N")
-    with r3:
-        st.metric("Charge", thresholds["charge"])
-
-    impact_df = pd.DataFrame({
-        "Variable": [
-            "Force talon normalisée",
-            "Force talon convertie",
-            "Pression talon",
-        ],
-        "Zone basse / faible": [
-            f"< {thresholds['force_bw_low']:.2f} BW",
-            f"< {thresholds['force_n_low']:.1f} N",
-            f"< {thresholds['pression_low']} N/cm²",
-        ],
-        "Zone attendue": [
-            f"{thresholds['force_bw_low']:.2f} à {thresholds['force_bw_high']:.2f} BW",
-            f"{thresholds['force_n_low']:.1f} à {thresholds['force_n_high']:.1f} N",
-            f"{thresholds['pression_low']} à {thresholds['pression_high']} N/cm²",
-        ],
-        "Zone haute / élevée": [
-            f"> {thresholds['force_bw_high']:.2f} BW",
-            f"> {thresholds['force_n_high']:.1f} N",
-            f"> {thresholds['pression_high']} N/cm²",
-        ],
-    })
-
-    dynamique_df = pd.DataFrame({
-        "Variable": [
-            "Cadence",
-            "Temps de contact",
-            "Temps de vol",
-        ],
-        "Zone basse / faible": [
-            f"< {thresholds['cadence_low']} pas/min",
-            f"< {thresholds['contact_low']} %",
-            f"< {thresholds['flight_low']} %",
-        ],
-        "Zone attendue": [
-            f"{thresholds['cadence_low']} à {thresholds['cadence_high']} pas/min",
-            f"{thresholds['contact_low']} à {thresholds['contact_high']} %",
-            f"{thresholds['flight_low']} à {thresholds['flight_high']} %",
-        ],
-        "Zone haute / élevée": [
-            f"> {thresholds['cadence_high']} pas/min",
-            f"> {thresholds['contact_high']} %",
-            f"> {thresholds['flight_high']} %",
-        ],
-    })
-
-    symetrie_df = pd.DataFrame({
-        "Variable": [
-            "Asymétrie force talon",
-            "Asymétrie force avant-pied",
-            "Asymétrie COP",
-            "Différence rotation G/D",
-        ],
-        "Zone faible": [
-            f"< {thresholds['asym_low']} %",
-            f"< {thresholds['asym_low']} %",
-            f"< {thresholds['asym_low']} %",
-            f"< {thresholds['rotation_low']}°",
-        ],
-        "Zone modérée": [
-            f"{thresholds['asym_low']} à {thresholds['asym_high']} %",
-            f"{thresholds['asym_low']} à {thresholds['asym_high']} %",
-            f"{thresholds['asym_low']} à {thresholds['asym_high']} %",
-            f"{thresholds['rotation_low']} à {thresholds['rotation_high']}°",
-        ],
-        "Zone marquée": [
-            f"> {thresholds['asym_high']} %",
-            f"> {thresholds['asym_high']} %",
-            f"> {thresholds['asym_high']} %",
-            f"> {thresholds['rotation_high']}°",
-        ],
-    })
-
-    s1, s2, s3 = st.tabs(["Impact", "Dynamique", "Symétrie"])
-    with s1:
-        st.dataframe(impact_df, hide_index=True, use_container_width=True)
-    with s2:
-        st.dataframe(dynamique_df, hide_index=True, use_container_width=True)
-    with s3:
-        st.dataframe(symetrie_df, hide_index=True, use_container_width=True)
-
-    st.write(
-        "Ces seuils sont individualisés à partir du poids et du volume horaire d’entraînement hebdomadaire. "
-        "Ils constituent une base de référence personnelle avant toute analyse ou recommandation."
-    )