Blog SIG & Territoires

Dans l’article précédent nous avons vu comment Intégrer des prévisions météo marines et terrestres gratuites dans ArcGis.
Nous allons voir ici comment les intégrer dans un projet QGis. Il s’agit de données de prévisions marines à 8 jours. Ces données sont fournies par les sociétés NASCA et ACTIMAR, via les produits enav-forecast. Outre les logiciels qui utilisent ces prévisions et qui permettent leur consultation et analyse (enav-viewer, enav-forecast, enav-forecast pour iPhone, enav-Pilot), ces entreprises mettent aussi à votre disposition un accès aux données par le biais d’une connexion WFS.

Quelles données sont disponibles?

En ce qui concerne le type de données, vous accédez aux prévisions à 8 jours , toutes les 3 heures, pour les paramètres suivants:

• hauteur (en m),direction (en °) et période des vagues (en s)
• mer croisée (direction des vagues en °,direction de la houle en °)
• vitesse (en nds) et direction (en °) du vent moyen
• vitesse (en nds) des rafales de vent
• vitesse (en nds) et direction (en °) du courant de surface
• température (en °C) de la surface de la mer
• température (en °C) de l’air
• pluviométrie (en mm)

Quelle maille de données?

La maille des prévisions gratuites est de 0,5°. D’autres mailles sont disponibles mais sous condition d’abonnement (voir les abonnements).
Les prévisions couvrent le monde entier.

Fréquence de mise à jour

Les prévisions sont mises à jour quatre fois par jour.

Établir la connexion WFS

Pour établir une connexion WFS vous cliquez sur le bouton « Ajouter une couche WFS », puis cliquer sur le bouton « Nouveau » pour créer les paramètres de connexion.

Pour le nom tapez un nom suffisamment parlant : Météo Marine NASCA, par exemple.
Pour l’URL, rentrez la chaîne suivante :
Cliquez sur OK pour fermer la configuration.
La connexion apparaît maintenant dans la liste des connexions possibles.

Avant de vous connecter

Ouvrez un projet avec le SCR Géographique WGS84. Affichez un fond de carte qui vous permette de zoomer sur la zone qui vous intéresse. Quand vous allez vous connecter sur le serveur WFS, QGis va envoyer une requête avec les coordonnées des quatre coins de votre affichage actuel. Le serveur enverra les prévisions trouvées à l’intérieur de ce rectangle vers votre machine.

Si vous êtes en projection et non en géographique, ni QGis ni le serveur feront une transformation de coordonnées.Vous n’aurez aucun résultat en retour, souvent, pas même un message d’erreur.

Il faut donc vérifier que ce que vous allez transmettre sont bien des coordonnées en degrés de latitude et longitude. Pour cela, simplement, regardez la barre d’état en bas de la fenêtre QGis: vous avez un champ « coordonnées ». Les valeurs de ce champ doivent être comprises entre -180 et 180. Si c’est le cas vous êtes bien en coordonnées géographiques.

Si vous êtes bien en coordonnées géographiques mais que vous ne zoomez pas sur la zone d’intérêt, vous risquez d’avoir une zone trop grande pour la requête. Outre le fait que cela prendra beaucoup de temps, les serveurs WFS possèdent toujours une limite en nombre d’enregistrements. Si elle est atteinte, le serveur envoie les enregistrements trouvés, puis s’arrête. Vous n’aurez pas de message d’avertissement, simplement les données seront tronquées.

Connexion au serveur WFS

Une nouvelle connexion apparaît maintenant dans la liste des connexions proposées quand vous cliquez sur le bouton « Ajouter une couche WFS ». Cliquez sur le bouton « Connexion ».

Vous n’allez utiliser que la couche : view_meteo99

Cliquez sur cette ligne pour la sélectionner. Décochez la case « Cache feature ».
Cliquez sur « Ajouter ».

Affichage des données des prévisions

Après quelques instants, la nouvelle couche est ajoutée dans la fenêtre cartographique de QGis.

Chaque « point » de prévisions est en réalité une collection de 66 points. Pour pouvoir visualiser les prévisions pour une heure donnée, vous devrez effectuer une sélection sur l’ensemble des données de chaque point. Faire ce travail sur les requêtes vers le serveur est trop compliqué et cela ne vaut pas la peine. La méthode la plus efficace est d’afficher l’ensemble des données dans QGis et de créer une nouvelle couche sur votre poste de travail.

Vous pouvez le faire directement en cliquant sur la légende de la couche et en sélectionnant « Sauvegarder sous… ». Mais ne vous précipitez pas, il y a d’autres petits détails à régler, et nous allons le faire d’un seul coup.

Les connexions WFS avec QGis ont un petit défaut: selon la première valeur trouvée pour un champ attributaire de la table , il définit le type de donnée. Souvent on se retrouve avec des champs numériques dans la table d’origine du serveur, mais qui sont interprétés comme des chaînes de texte par QGis, en local. Avec les prévisions marines, le problème ne se pose pas pour les valeurs de force et vitesse, qui sont toujours de s valeurs avec le point décimal. Par contre pour les valeurs de direction, qui sont des entiers, le problème se pose.
Vous aurez des valeurs de direction sous forme de chaîne de texte.Si pour un simple affichage cela ne change pas grande chose, pour afficher des flèches orientées (courants, vagues, vent) il est impossible de demander une rotation du symbole avec un champ texte. Il faut absolument que le champ direction soit numérique.

Pour résoudre ce problème, il faut créer un nouveau champ numérique et utiliser la calculatrice de champ pour calculer les valeurs à partir des chaînes de texte, et ce, pour les trois champs de direction (vagues, vent, courant).

Travail fastidieux. Donc, il faut trouver une solution moins fastidieuse. Nous avons construit un modèle de traitement avec le « Modeleur graphique » de QGis. C’est un petit fichier, que vous téléchargez en cliquant ici.

Décompressez le fichier dans un répertoire de votre choix.

Ouvrez le Modeleur Graphique:

Dans la fenêtre du modeleur, cliquez sur le premier bouton pour ouvrir le modèle téléchargé. Pointez sur le fichier meteoNasca.model.
Le modèle se charge.

Pour l’exécuter, cliquer sur le bouton de droite « Lancer le modèle ».

Vérifiez que le fichier en entrée est bien la couche « view_meteo99 ».
Parcourez les répertoires pour définir l’emplacement et le nom du fichier shape qui sera créé sur votre poste.Sélectionnez « Save to file »
Cliquez sur Run
La nouvelle couche est chargée dans la fenêtre cartographique. Vous pouvez maintenant supprimer la couche WFS, ça évite de chargements chaque fois que vous changer l’emprise de la vue.

Symbologie des données des prévisions

Pour ce qui est des données simples (température, pluviométrie,…) la seule chose à ne pas oublier c’est que vous avez 66 points superposés à chaque endroit de prévision. Vous devez passer par une « requête de sous-ensemble » pour n’afficher qu’un pas de temps à la fois :
Propriétés de la couche -> Onglet « Général » -> Sous-ensemble (cliquez sur constructeur de requête)

Vous avez deux attributs, le premier pour la date du jour, le deuxième pour l’heure du jour. ATTENTION! les heures sont en GMT pas en heure locale.

Avant d’appliquer les symbologies, n’oubliez pas de définir la requête de sélection. Si non, vous aurez un empilement de symboles sur chaque point.

Pour les données complexes (vent, vagues, courant) vous avez au moins deux valeurs associées : direction et force. Le moyen classique d’afficher ce type de données est d’utiliser une flèche orientée selon la direction du phénomène et de lui attribuer soit une taille, soit une couleur différente selon la force.

Première chose à savoir: Comment est codée la direction du phénomène?
Toutes les directions présentes dans les attributs de la table indiquent la direction à partir de laquelle vient le phénomène (la direction d’où il vient).
Pour orienter les flèches, il faut appliquer au symbole une rotation. Si vous choisissez un symbole « flèche vers le haut » et vous rentrez comme champ de rotation seulement le nom du champ direction, vos flèches indiqueront la direction contraire! Il faut, dans ce cas, utiliser une expression, pour la rotation, de type [champ direction]- 180.

Nous avons préparé trois fichiers de couche avec une symbologie standard qui devraient faire l’affaire. Vous pouvez les télécharger en cliquant ici.
Vous avez, dans le fichier compressé, trois fichiers de styles de couches de Qgis (vent, courant et vagues) et un dossier avec les symboles utilisés pour le vent sur les cartes météo (barbules).
Les trois premiers vous pouvez les extraire dans le répertoire de votre choix. Par contre le dossier « barbules » doit être copié dans le répertoire d’installation de QGis/apps/qgis/svg . Si non, les symboles ne seront pas accessibles pour les fichiers de style.

La symbologie pour le vent est la suivante:

Les symboles utilisés sont les barbules standard, avec des incréments de 5 en 5 nœuds jusqu’à 65 nœuds.
Le champ rotation est déjà corrigé pour avoir la bonne orientation des barbules.

La symbologie pour le courant est la suivante:

Les symboles utilisés sont des flèches , avec des incréments de 1 nœud. La taille des flèches est croissante en fonction de la vitesse. De plus les vitesses jusqu’à 3 nœuds sont en bleu, entre 3 et 5 sont en orange et à partir de 6 nœuds en rouge.
Le champ rotation est corrigé de 180° pour avoir la bonne orientation des flèches.
Pensez à ajouter dans la requête de sélection, en plus de la date et heure, un filtre sur la direction et vitesse pour ne pas afficher les valeurs nulles à terre (CourantVal<>0 AND ICurDir<>0).

La symbologie pour les vagues est la suivante:

Les symboles utilisés sont des flèches , avec des incréments de 1 mètre. Les hauteurs jusqu’à 3 mètres sont en vert, entre 4 et 7 sont en jaune et à partir de 7 nœuds en orange-rouge.
Le champ rotation est corrigé de 180° pour avoir la bonne orientation des flèches.
Pensez à ajouter dans la requête de sélection, en plus de la date et heure, un filtre sur la direction et vitesse pour ne pas afficher les valeurs nulles à terre (HsVal<>0 AND IHsDir<>0).
Les étiquettes indiquent la hauteur significative des vagues en mètres, et la période des vagues en secondes.

Si vous avez des questions, n’hésitez pas à laisser un commentaire!

Si cet article vous a intéressé et que vous pensez qu'il pourrait bénéficier à d'autres personnes, n'hésitez pas à le partager sur vos réseaux sociaux en utilisant les boutons ci-dessous. Votre partage est apprécié !

Atilio Francois