Principe de la passerelle ArcGrid/Arcinfo / FTS63 DANS VISHYR

La passerelle est accessible dans le module Vishyr, par la commande :

Fichier/ Import-Export/ Calamar-Panther à fts63

 

 

 

 

 

Objectif

Constituer un fichier d'épisodes de pluie FTS63 à partir de fichiers radar au format ascii GeoTiff (Panthere/Antilope) ou au format Calamar®.

 

 

Principes de la passerelle

Un fichier radar est une matrice de valeurs de pluies cumulées sur un pas de temps donné, aux nœuds d'une grille régulière. Chaque nœud représente un pixel radar. Un fichier radar est donc caractérisé par sa résolution temporelle (généralement 5 minutes) et sa résolution spatiale (pixels généralement espacés d'1 km). Chaque pixel radar va être converti en un pluviomètre virtuel, dont les données seront écrites dans un fichier 63.

 

Le principe de la passerelle consiste à définir des listes de fichiers radar correspondant à des pas de temps successifs, chaque liste représentant un épisode distinct, associé à une date de départ (calculée à partir du 1^er élément de la liste) et une date de fin (calculée à partir du dernier élément de la liste). Si les pas de temps associés aux fichiers radar ne se succèdent pas exactement dans la liste, les pluies correspondant aux pas de temps manquants seront signalées par des lacunes.

 

Sauf mention contraire, tous les pixels contenus dans le fichier radar seront intégrés au fichier 63. Il est possible de sélectionner une partie seulement des pixels, soit par délimitation d'un cadre rectangulaire définis par les coordonnées de ses sommets, soit par utilisation d'un masque (correspondant par exemple au contour du bassin).

 

Formats des fichiers

 

Les fichiers radar lus sont soit des fichiers GeoTiff (Panthere, Antilope de Meteo-France), soit des fichiers Calamar® de la société Rhea. Pour être reconnus, les noms de ces fichiers doivent impérativement être de la forme :

- aaaaMMjjhhmm.grd ou autre extension (attention le . est obligatoire) pour les fichiers GeoTiff

- xxZZ-aaaaMMjjhhmm.dat ou autre extension (où ZZ, par exemple 05, indique le no de zone calamar, et xx représentent 2 caractères quelconques) pour les fichiers Calamar®

 

De plus, les en-têtes des fichiers GeoTiff doivent être de la forme (les valeurs numériques sont des exemples) :

NCOLS 530

NROWS 528

XLLCORNER 529425.307217329741

YLLCORNER 1663253.926465130411

CELLSIZE 1000.000000000000

 

NB : Les majuscules doivent être respectées, ainsi que l’espace unique entre le texte et la valeur numérique. Les en-têtes comportant XLLCENTER ou YLLCENTER sont également gérées.

 

Le type de ces fichiers est automatiquement déterminé. Il n'est pas possible de traiter simultanément :

-          des fichiers de types différents GeoTiff et Calamar©

-          des fichiers Calamar® associés à différentes zones

 

Limitations de traitement

 

Dans la version ATHYS 2.0, les limitations de traitement sont les suivantes :

-          Nombre maximum d’événements à constituer : 100

-          Nombre maximum de pas de temps d'un épisode : 50000

-          Nombre maximum de pixels du fichier radar : 80000

 

Les fichiers 63 en sortie du traitement ne doivent pas contenir plus de 9999 pixels ( = postes pluviométriques).

 

En pratique

La 1^ère opération consiste à sélectionner les fichiers radar associés à 1 épisode de pluie à

Cette icône active un navigateur qui permet de sélectionner les fichiers radar. Les options de sélection multiple sont utilisables dans ce navigateur.

 

Après avoir validé le choix des fichiers radar, la liste des fichiers radar apparaît dans la fenêtre de gauche. Il est possible d'additionner plusieurs groupes de fichiers radar (ré-activer + pour chaque nouveau groupe) pour constituer la liste. Il est également possible de supprimer des fichiers radar de la liste : on sélectionne les fichiers à supprimer et on active l'icône .

 

 

 

Une fois la liste des fichiers radar constituée, la 2^ème opération consiste à valider la définition de l'épisode

 

 

La validation de l'épisode réinitialise la liste, pour constituer l'épisode suivant. Simultanément, la liste des épisodes est mise à jour, et résume les caractéristiques des épisodes déjà constitués. De plus, pour les fichiers de type GeoTiff, les coordonnées des images sont automatiquement affichées lorsqu'on valide le 1^er épisode.

 

 

 

 

Dans le cas des fichiers Calamar®, le processus est identique, mais il faut préalablement déclarer le nom du fichier des zones Calamar®, définissant les coordonnées de chaque zone. Ces coordonnées doivent correspondre à la projetction géographique qui sera ultérieurement utilisée (une conversion des coordonnées est éventuellement nécessaire)

 

 

Exemple fichier zone Calamar©

 

<?xml version="1.0" ?>

<defzones>

<zone num="01" xlambert="738.7" ylambert="3185.1" nbcol="50" nblig="50" pas="1"/>

<zone num="02" xlambert="705.7" ylambert="3195.1" nbcol="50" nblig="50" pas="1"/>

<zone num="03" xlambert="697.7" ylambert="3241.1" nbcol="64" nblig="64" pas="1"/>

<zone num="04" xlambert="742.7" ylambert="3232.1" nbcol="64" nblig="64" pas="1"/>

<zone num="05" xlambert="768.7" ylambert="3268.1" nbcol="64" nblig="64" pas="1"/>

<zone num="06" xlambert="797.7" ylambert="3257.1" nbcol="64" nblig="64" pas="1"/>

<zone num="07" xlambert="806.7" ylambert="3214.1" nbcol="64" nblig="64" pas="1"/>

<zone num="08" xlambert="713.7" ylambert="3283.1" nbcol="64" nblig="64" pas="1"/>

<zone num="09" xlambert="744.9" ylambert="3326.1" nbcol="64" nblig="64" pas="1"/>

</defzones>

 

On réitère les opérations 1 et 2 pour chaque nouvel épisode à constituer.

 

Une fois l'ensemble des épisodes constitués, la 3^ème opération consiste à limite éventuellement le nombre de pixels radar à traiter. Cette réduction peut être réalisée :

-          en définissant une sous-région de la région qui correspond à l'image radar

-          en utilisant un fichier masque : seuls les pixels associés à des mailles dont le géocode est strictement positif doit être traité, les autres ne figureront pas dans le fichier 63. Le fichier masque doit être un fichier géographique. Le format Lamont (non géoréférencé) n'est pas accepté.

-          en utilisant un facteur d'échantillonnage N, qui ne retiendra qu'1 pixel sur NxN (1 pixel sur 4 dans le cas d'un facteur échantillonnage 2). La valeur du pixel échantillonné est égale à la moyenne des valeurs des pixels du bloc correspondant.

 

Ensuite, il s'agit de nommer les fichiers de sortie du traitement :

-          le fichier 63 de destination, contenant les données postes pluviométriques virtuels correspondant aux pixels sélectionnés, pour les épisodes définis

-          le fichier des coordonnées des pixels, contenant la liste des pixels sélectionnés, leur code numérique, leur nom, leurs coordonnées en X et en Y dans la projection géographique utilisée dans les fichiers radar, leur type (ici P pour Pluviomètre). Ce fichier est très utile ultérieurement, pour remplir automatiquement la liste des coordonnées des pixels radar, dans le menu Pluie-débit de Mercedes

 

NB : les codes numériques identifiant les pixels dans le fichier 63 sont automatiquement formés à partir de leurs coordonnées. Ces codes comprennent 10 chiffres, correspondant à xxxxxyyyyy, où xxxxx désignent les 5 1^ers chiffres de la cordonnées X, et yyyyy les 5 1^ers chiffres de la coordonnée Y

 

Enfin, le traitement est activé par le bouton . Avant exécution, le nombre de pixel à traiter est contrôlé, et doit être inférieur à 9999. Sinon, le traitement n'est pas effectué, et il est nécessaire de réduire le nombre de pixels (fenêtrage, masque ou échantillonnage) avant de relancer le traitement. Une confirmation est demandée, et si le fichier 63 existe déjà.

 

 

 

 

Après exécution, les options de traitement ne sont pas effacées, pour relancer un nouveau traitement après modifications. Le bouton permet de réinitialiser une nouvelle session.