===== 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 [[http://www.ngdc.noaa.gov/autochart/|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. {{:noaa:2013-12-13_09h39_59.png|}} 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 [[http://maps.ngdc.noaa.gov/viewers/multibeam/|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 {{:noaa:mbgrid.xlsx|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 === == Pour les format .mbXX == 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 == Pour les autres formats (plus anciens ex. .r0x) == - ls | grep .*$ > tmplist - 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 [[http://www.ldeo.columbia.edu/res/pi/MB-System/html/mbsystem_formats.html|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 [[http://www.mbari.org/data/mbsystem/html/mbgrid.html|mbgrid]]. == Pour la bathymétrie == mbgrid -A2 -R/LonW/LonE/LatS/LatN -CS -DTX/TY -G4 -Idatalist -Onomfichier * 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"). == Pour la réflectivité == mbgrid -A3 -R/LonW/LonE/LatS/LatN -CS -DTX/TY -G4 -Idatalist -Onomfichier Exemples de grilles qui fonctionnent : * Mer d'Arabie : mbgrid -A2 -R/41/67/0/26 -C5 -D14430/14430 -G4 -Idatalist -Onomfichier * Golfe de Gascogne : mbgrid -A2 -R/345/347/45/47 -C5 -D1000/1000 -G4 -Idatalist -Onomfichier * Monde : mbgrid -A2 -R/90/270/-90/90 -C5 -D4000/4000 -G4 -Idatalist -Onomfichier * Golfe Bengal : mbgrid -A2 -R/79/94/3/18 -C5 -D1665/1665 -G4 -Idatalist -Onomfichier * Golfe Bengal : mbgrid -A2 -R/79/94/3/18 -C5 -D16650/16650 -G4 -Idatalist -Onomfichier * Réflectivité Golfe Bengal : mbgrid -A3 -R/79/94/3/18 -C5 -D16650/16650 -G4 -Idatalist -Onomfichier === 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.