Outils pour utilisateurs

Outils du site


noaa:start

Ceci est une ancienne révision du document !


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

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)

  • gunzip -Nv *.tgz
  • fichier.tgz: 57.4% – replaced with fichier.tar
  • tar xvf *.tar
  • déplacez vous dans le répertoire contenant tout les .mb, s'il sont compressés gunzip -Nv *.gz

Création de la liste des fichiers

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

  • ls -1
  • ls | grep .mb56$ > tmplist (Attention de mettre le bon n° de MB)
  • mbdatalist -F-1 -I tmplist > datalist

Conversion des données au format ascii d’ArcGIS

Générez le fichier ascii au format ArcGIS à partir de la commande :

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

Vous devez renseigner les termes mentionnés en gras :

LonW, LonE, LatS et LatN : bords de la grille S : taille maximale d’interpolation de Spline (nombre de lignes/colonnes interpolées) TX et TY : dimensions de la grille nomfichier : au nom du fichier ascii créé.

TX et TY sont calculés au préalable (cf. section “Calcul de la taille de la grille à générer”).

Exemple de grille qui fonctionnent :

  • Mer d'Arabie : mbgrid -R/41/67/0/26 -C5 -D14430/14430 -G4 -Idatalist -Omnt
  • Golfe de Gascogne : mbgrid -R/345/347/45/47 -C5 -D1000/1000 -G4 -Idatalist -Otest10
  • Monde : mbgrid -R/90/270/-90/90 -C5 -D4000/4000 -G4 -Idatalist -Otest
  • Golfe Bengal : mbgrid -R/79/94/3/18 -C5 -D1665/1665 -G4 -Idatalist -Grid1000m
  • Golfe Bengal : mbgrid -R/79/94/3/18 -C5 -D16650/16650 -G4 -Idatalist -Grid100m

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.

noaa/start.1387360199.txt.gz · Dernière modification : de zaragosi