Table des matières

Intégration sous ArcGIS des données SMF disponibles sur le site de la NOAA

Solution 1 : par zone

Pour récupérer les données par zone c'est ici

Rentrez une zone de 1° par 1° et dans les options de la grille mettre 120 m en taille de cellule. Rentrez votre mail et quelques heures plus tard vous recevrez un lien pour télécharger vos données.

(afin de vérifier la qualité des données vous pouvez réaliser une première extraction des données sur toute votre zone avec une résolution faible)

Téléchargez l'archive au format autochartxxxx.tar.gz, décompressez là et récupérez à l'intérieur le fichier “output.grd”

Sous ArcMAP ouvrez l'outils “Générer une couche raster NetCDF” (dans les outils multidimensionnels) pour convertir la grille au format grid d'ArcGIS.

Une fois les différentes dalles de votre zone chargées sous ArcMap utilisez l'outil “Mosaïque vers nouveau raster” pour les combiner.

Solution 2 par mission (indispensable si vous voulez également la réflectivité)

Pour récupérer les données par mission c'est ici

Calcul de la taille de la grille à générer

Avant de convertir les données téléchargées, vous devez calculer les dimensions TX et TY en coordonnées projetées de la grille XY à générer. Pour cela, utilisez l’emprise géographique de votre zone d’intérêt qui doit être strictement carrée.

TY = (LatN - LatS) * D * T TX = (LonE - LonW) * D * T

T est la taille de la maille élémentaire (en mètres). Elle est obligatoirement identique en X et en Y.

D est la valeur en mètres de 1° de latitude à l’équateur. Généralement D sera arrondi à 111000 mais vous pourrez moduler cette valeur.

Le NORD est le référentiel en Y, l’EST le référentiel en X. Pensez donc à convertir vos coordonnées (ex : 8°W = 352°E ; 10°S = -10°N).

Exemple pour une zone comprise entre 37,9°N/38°N et 9,3°W/9,4°W avec une résolution souhaitée (taille de la maille élémentaire) de 20m :

TY = (38-37,9)*111000*10=555 TX = (360-9,4-360-9,3)*111000*20=555

Fichier Excel pour faire les calculs

Décompression des données sous Linux (DGO5)

Création de la liste des fichiers

Pour les format .mbXX

Dans un terminal mettez vous dans le répertoire contenant les fichiers mbXX. Lancez les commandes suivantes :

Pour les autres formats (plus anciens ex. .r0x)
  1. ls | grep .*$ > tmplist
  2. Je n'ai pas trouvé comment faire marcher mbdatalist, il faut faire la liste à la main en rentrant le format (ex pour URI Sea Beam Format ID: 14). Liste des formats ici.

Conversion des données au format ascii d’ArcGIS

Pour générer le fichier ascii au format ArcGIS il faut utiliser la fonction mbgrid.

Pour la bathymétrie

mbgrid -A2 -R/LonW/LonE/LatS/LatN -CS -DTX/TY -G4 -Idatalist -Onomfichier

Pour la réflectivité

mbgrid -A3 -R/LonW/LonE/LatS/LatN -CS -DTX/TY -G4 -Idatalist -Onomfichier

Exemples de grilles qui fonctionnent :

Script pour lancer plusieurs traitement sous linux

Sous linux créer un fichier script.sh

#!/bin/bash

cd /Sedim2/seb/Noaa/KNOX04RR/multibeam/data/version1/MB/

ls -1

ls | grep .mb56$ > tmplist

mbdatalist -F-1 -I tmplist > datalist

mbgrid -A2 -R/64/79/5/20 -C5 -D16650/16650 -G4 -Idatalist -OmntKNOX04RRa

A lancer avec ./script.sh

Conversion des données ascii en raster

Une fois le fichier ascii récupéré sous Windows, utilisez l’outil de conversion ASCII vers raster de l’ArcToolbox d’ArcCatalog ou d’ArcMap. Veillez à renseignez FLOAT pour le type de données en sortie. Appliquez le système géodésique WGS84 au raster créé en allant dans les propriétés de référence spatiale de ce dernier.