260 likes | 417 Views
CASH. C ontribution de l’ A ltimétrie S patiale à l’ H ydrologie. Franck MERCIER Ouan-Zan ZANIFE Bruno PICARD CLS, Direction Océanographie Spatiale. Retracking Topex. 278 cycles traités 27 cycles récupérables (10 %). Changement d’orbite. Retracking Topex.
E N D
CASH Contribution de l’Altimétrie Spatiale à l’Hydrologie Franck MERCIER Ouan-Zan ZANIFE Bruno PICARD CLS, Direction Océanographie Spatiale
Retracking Topex 278 cycles traités 27 cycles récupérables (10 %) Changement d’orbite
Retracking Topex Accès aux sorties du retracking (FD-GDR: format ESA pour ENVISAT): ftp ftp.cls.fr User: ftp Password: email cd pub/oceano/cash/L2 cd BBB (BBB = AMZ ou PLT, CGO, GAG, MKG, YLW, DNB, LNA) get cXXX_BBB.tar.gz (XXX = n° cycle) Contient toutes les traces sur le bassin choisi pour le cycle considéré, à raison d’un fichier par trace usage interne au projet CASH !
Situation chaîne de traitement • Implémentation Jason: • Développement sous SUN/SOLARIS validé • Portage sous Linux validé • résultats préliminaires sur trace 063 (Manaus) peu encourageants • tests sur 4 traces supplémentaires (076, 165, 241 et 254) sur le Rio Negro • réunion avec le CNES (N. Picot) pour approfondir le problème • altimètre décroche pas de formes d’onde • modification logiciel bord ? • investigations complémentaires: visualisation de la télémesure pb dans les algos de compression/décompression de la forme d’onde
Analyse du retracking: JASON (2 cycles en août-sep 2005) Amazone: Negro
Analyse du retracking: JASON (4 cycles en été 2005:125, 126, 134, 135) Amazone: Negro , Trace 241 Retracking Ice2
Visualisation TM et SGDR (T063, C134 - 28 Aout 2005) CASH TM CCI CMA Proto EnviSat L1B L2 1.0 SGDR APID
Visualisation TM et SGDR (T063, C134 - 28 Aout 2005) Amazone: Manaus
Visualisation TM et SGDR (T063, C134 - 28 Aout 2005) Amazone: Manaus
Visualisation TM et SGDR (T063, C134 - 28 Aout 2005) Amazone: Manaus
Visualisation TM et SGDR (T063, C134 - 28 Aout 2005) Amazone: Manaus
Visualisation TM et SGDR (T063, C134 - 28 Aout 2005) Amazone: Manaus
Suite Jason • Retracking de 20 traces par cycle, de juin 2002 à décembre 2002 • Décision pour le traitement de l’archive Jason
…….. FD-GDR Trace t Cycle c FD-GDR Trace t’ Cycle c FD-GDR Trace t Cycle c’ FD-GDR Trace t’ Cycle c’ Situation chaîne de traitement • Données altimétriques Données hydrologiques • Schéma retenu pour implémentation à CLS: • Procédure développée et testée sur Trace 063 (Manaus) en novembre • Paramétrage à finaliser dès la livraison, par le LEGOS, de la liste des SV: toujours en attente … Mis à disposition sur Internet SV N Série tempo NetCDF Etape 2: Calcul Série Temporelle hauteur d’eau Développement d’outils spécifiques à partir d’outils existants à CLS, éprouvés en océano (éditing, lissage, filtrage, projection sur trace moyenne … ) Table de mesures CLS …. …. Etape 1: Remplissage Table de mesures SV M Série tempo NetCDF
H (m) c54 c352 c12 ST(t) t Données altimétriques données hydrologiques • Point de départ : • Nb données géoréférencées pour N passages au-dessus de la SV (Nb ≠ pour chaque cycle) • ST(t) = série temporelle de départ • st(t) = série temporelle pour 1 cycle donné • But : • 1 valeur par cycle • STR(t) : série temporelle réduite
H (m) c54 c352 c12 t Données altimétriques données hydrologiques • Traitement 0 • Suppression « manuelle » des valeurs extrêmes
Données altimétriques données hydrologiques • Traitement 0 • Suppression « manuelle » des valeurs extrêmes • Traitement 1 / Etape 1 • Suppression des valeurs extrêmes restantes: • filtre : H <= moy(ST(t)) + 15 m (paramétrable) H (m) c54 c352 c12 15 m t
H (m) c54 c352 c12 15 m t Données altimétriques données hydrologiques • Traitement 0 • Suppression « manuelle » des valeurs extrêmes • Traitement 1 / Etape 1 • Suppression des valeurs extrêmes restantes: • filtre : H <= moy(ST(t)) + 15 m • Traitement 1 / Etape 2 • ajustement linéaire staj(t) des st(t) de chaque cycle H (m) c54 st(t) staj(t) t
H (m) c54 st(t) > 2m > 2m t Données altimétriques données hydrologiques • Traitement 1 / Etape 3 • filtre : st(t) – staj(t) <= 2 m (paramétrable) : on obtient str(t) • Traitement 1 / Etape 4 • STR(t) = moy(str(t)) pour chaque cycle • Quality = stdev(str(t)) pour chaque cycle H (m) c54 str(t) t
Données altimétriques données hydrologiques • Traitement 1 / Etape 3 • filtre : st(t) – staj(t) <= 2 m : on obtient str(t) • Traitement 1 / Etape 4 • STR(t) = moy(str(t)) pour chaque cycle • Quality = stdev(str(t)) pour chaque cycle • Traitement 1 / Etape 5 • SI Nb<=2 après Filtres Etape 1 & 3 ALORS • STR(ti) = moyenne pondérée str(t) et (STR(ti-1) + STR(ti+1) )/2 • & Quality = 0 H (m) H (m) c54 c54 c352 c352 STR(t) STR(t) c12 c12 t
Données altimétriques données hydrologiques: données alti CASH: 2800 points GDR: 1600 points
Données altimétriques données hydrologiques: données hydro CASH et GDR: extraction entre -3.21 et -3.14 CASH_2: extraction entre -3.20 et -3.15
Données altimétriques données hydrologiques: données hydro
Données altimétriques données hydrologiques: données hydro