img2grdgmt - En ligne dans le Cloud

PROGRAMME:

Nom

img2grd - Extraire le sous-ensemble du fichier img au format Mercator ou Géographique

SYNOPSIS

img2grd fichier img fichier grd région type [ ] [ [minlat/maxlat] ] [ ] [ minutes ] [ ] [ navigation
] [ [escaliers] ] [ [niveau] ] [ maxlon ] [ -n]

Remarque: Aucun espace n'est autorisé entre l'indicateur d'option et les arguments associés.

DESCRIPTION

img2grd lit un fichier au format img, en extrait un sous-ensemble et l'écrit dans un fichier de grille. Les -M
L'option dicte si la projection sphérique de Mercator du fichier img est
conservée ou si une grille géographique doit être écrite en annulant la projection de Mercator. Si
la grille géographique est sélectionnée, vous pouvez également demander un rééchantillonnage sur la -R donné.

EST REQUIS ARGUMENTS

fichier img
Un fichier au format Mercator img tel que les champs de gravité marine ou de topographie des fonds marins
estimée à partir des données altimétriques satellitaires de Sandwell et Smith. Si l'utilisateur a défini
une variable d'environnement $GMT_DATADIR, puis img2grd va essayer de trouver fichier img in
$GMT_DATADIR; sinon il essaiera de s'ouvrir fichier img directement.

-Gfichier grd
fichier grd est le nom du fichier de grille de sortie.

-R[unité]ouest/est/sud/nord[/zmin/zmax][r]
ouest, est, sudet nord spécifiez la région d'intérêt, et vous pouvez spécifier
en degrés décimaux ou au format [+-]dd:mm[:ss.xxx][W|E|S|N]. Ajouter r si inférieur
les coordonnées de gauche et en haut à droite de la carte sont données au lieu de w/e/s/n. Les deux
raccourcis -Rg ainsi que -Rd signifie domaine global (0/360 et -180/+180 en longitude
respectivement, avec -90/+90 en latitude). Alternativement pour la création de grille, donnez
Rcodelon/lat/nx/ny, Où code est une combinaison de 2 caractères de L, C, R (pour gauche,
au centre ou à droite) et T, M, B pour le haut, le milieu ou le bas. par exemple, BL pour le coin inférieur gauche.
Ceci indique quel point sur une région rectangulaire le lon/lat la coordonnée fait référence
à, et les dimensions de la grille nx ainsi que ny avec des espacements de grille via -I est utilisé pour créer
la région correspondante. Vous pouvez également spécifier le nom d'un fichier de grille existant
et par -R les paramètres (et l'espacement de la grille, le cas échéant) sont copiés à partir de la grille.
En utilisant -Runité attend des coordonnées projetées (cartésiennes) compatibles avec celles choisies -J
et nous projetons inversement de déterminer la région géographique rectangulaire réelle. Pour
vue de perspective (-p), éventuellement ajouter /zmin/zmax. En cas de vue en perspective
(-p), une plage z (zmin, zmax) peut être ajouté pour indiquer la troisième dimension. Cette
doit être fait uniquement lors de l'utilisation du -Jz option, pas lorsque vous utilisez uniquement l'option -p option.
Dans ce dernier cas, une vue en perspective de l'avion est tracée, sans troisième
dimension.

EN OPTION ARGUMENTS

-C Définissez les coordonnées Mercator x et y par rapport au centre de projection [La valeur par défaut est
par rapport au coin inférieur gauche de la grille]. A besoin -M.

-RÉ[minlat/maxlat]
Utilisez la plage de latitude étendue -80.738/+80.738. Sinon, ajoutez
minlat/maxlat comme étendue de latitude du fichier img d'entrée. [La valeur par défaut est
-72.006/72.006]. Ce n'est généralement pas nécessaire puisque nous pouvons déterminer l'étendue à partir de
contrôle de la taille du fichier.

-E Peut être utilisé lorsque -M n'est pas défini pour forcer la grille finale à avoir exactement la même
région comme demandé avec -R. Par défaut, la région finale est une projection directe de
la région de Mercator d'origine et s'étendra généralement légèrement au-delà de la
plage de latitude demandée, et de plus l'incrément de grille en latitude ne
correspondre à l'incrément de longitude. Cependant, le rééchantillonnage supplémentaire introduit de petits
erreurs d'interpolation et ne doit être utilisé que si la grille de sortie doit correspondre à la
région demandée et avoir x_inc = y_inc. Dans ce cas, la région définie par -R doit être
donné en multiples de l'incrément (.par exemple, -R0/45/45/72).

-I Indiquer minutes comme la largeur d'un pixel img d'entrée en minutes de longitude.
[La valeur par défaut est 2.0]. Ce n'est généralement pas nécessaire car nous pouvons déterminer la taille des pixels à partir de
contrôle de la taille.

-M Sortez une grille de Mercator sphérique [La valeur par défaut est une grille géographique lon/lat]. Les
La projection sphérique de Mercator du fichier img est conservée, de sorte que la région -R
défini par l'utilisateur est légèrement modifié ; la région modifiée correspond aux arêtes
de pixels [ou groupes de navigation pixels]. L'en-tête du fichier de grille est défini de telle sorte que le x et
les longueurs de l'axe y représentent la distance des bords ouest et sud de l'image,
mesuré dans les unités par défaut de l'utilisateur, avec -Jm1 et l'ajusté -R. En réglant le
défaut PROJ_ ELLIPSOIDDE = Sphère, l'utilisateur peut faire des superpositions avec le -R
pour qu'ils correspondent. Voir EXEMPLES au dessous de. Le ajusté -R est aussi écrit dans la grille
remarque d'en-tête, afin qu'elle puisse être trouvée plus tard. Voir -C pour définir les coordonnées relatives à
centre de projection.

-Nnavigation Faites la moyenne des valeurs dans les pixels img d'entrée dans navigation by navigation carrés et créer
un pixel de sortie pour chacun de ces carrés. Si utilisé avec -T3 il rapportera une moyenne
contrainte entre 0 et 1. Si utilisé avec -T2 la sortie sera la valeur moyenne des données
ou NaN selon que la contrainte moyenne est > 0.5. navigation doit diviser également
dans les dimensions du fichier img en pixels. [Défaut 1 ne fait pas de moyenne].

-S[escaliers]
Multipliez les valeurs du fichier img par escaliers avant de stocker dans le fichier de grille. [La valeur par défaut est
1.0]. Pour les fichiers img récents : les fichiers img topo sont stockés dans des compteurs (corrigés) [-S1] ;
fichiers de gravité à l'air libre en mGal*10 [-S0.1 pour obtenir mGal] ; fichiers de déviation verticale dans
microradians*10 [-S0.1 pour obtenir des micro-radians], fichiers de gradient de gravité vertical dans
Éotvos*50 [-S0.02 pour obtenir des Eotvos, ou -S0.002 pour obtenir mGal/km]). Sinon escaliers est donné
nous essayons de déterminer l'échelle en examinant le nom du fichier pour trouver des indices.

-Ttype type gère l'encodage des informations de contrainte. type = 0 indique que non
ces informations sont encodées dans le fichier img (utilisé pour les versions antérieures à 1995 du
données de gravité) et obtient toutes les données. type > 0 indique que les informations de contrainte sont
encodé (versions 1995 et ultérieures (actuelles) des fichiers img) afin que l'on puisse
produire un fichier de grille comme suit : -T1 obtient des valeurs de données à tous les points, -T2 obtient des données
valeurs aux points contraints et NaN aux points interpolés ; -T3 obtient 1 à
points contraints et 0 aux points interpolés [La valeur par défaut est 1].

-V[niveau] (plus ...)
Sélectionnez le niveau de verbosité [c]. Particulièrement recommandé ici, car il est utile de voir
comment les coordonnées sont ajustées.

-Wmaxlon
Indiquer maxlon comme étendue de longitude maximale du fichier img d'entrée. Versions
depuis 1995 ont eu maxlon = 360.0, alors que certains fichiers antérieurs avaient maxlon = 390.0.
[La valeur par défaut est 360.0].

-^ or juste -
Imprime un court message sur la syntaxe de la commande, puis quitte (REMARQUE : sous Windows
utiliser juste -).

-+ or juste +
Imprimez un message d'utilisation détaillé (aide), y compris l'explication de tout
option spécifique au module (mais pas les options communes GMT), puis se ferme.

-? or aucune arguments
Imprimez un message d'utilisation (aide) complet, y compris l'explication des options, puis
sorties.

--version
Imprimer la version GMT et quitter.

--show-datadir
Affichez le chemin complet vers le répertoire de partage GMT et quittez.

GÉOGRAPHIQUE EXEMPLES

La -M L'option doit être exclue si vous avez besoin que la grille de sortie soit en zone géographique
coordonnées. Pour extraire des données dans la région -R-40/40/-70/-30 de monde_grav.img.7.2 ainsi que
reprojeter pour obtenir des coordonnées géographiques, vous pouvez essayer

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V

Étant donné que l'espacement des latitudes dans le fichier img est équidistant en unités Mercator, le
la grille résultante ne correspondra pas à celle spécifiée -R exactement, et l'espacement des latitudes ne sera pas
égal à l'espacement des longitudes. Si vous avez besoin d'une correspondance exacte avec votre -R et le même espacement
en longitude et latitude, utilisez le -E option:

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V

MERCATEUR EXEMPLES

Étant donné que les fichiers img sont dans une projection Mercator, vous ne devez PAS extraire une grille géographique
si votre plan est de faire une carte Mercator. Si vous avez fait cela, vous avez fini de projeter et
reprojection la grille, perdant les détails de courte longueur d'onde. Mieux vaut utiliser -M et tracer la grille
en utilisant une projection linéaire à la même échelle que la projection de Mercator souhaitée (voir GMT
Exemple 29). Pour extraire des données dans la région -R-40/40/-70/-30 de monde_grav.img.7.2, courir

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V

Notez que le -V l'option nous indique que la plage a été ajustée à
-R-40/40/-70.0004681551/-29.9945810754. Nous pouvons également utiliser grdinfo pour trouver que le fichier de grille
l'en-tête indique que sa région est -R0/80/0/67.9666667 C'est la plage de x,y que nous obtiendrons
à partir d'une projection sphérique de Mercator utilisant -R-40/40/-70.0004681551/-29.9945810754 et
-Jm1. Ainsi, pour prendre ship.lonlatgrav et l'utiliser pour échantillonner le merc_grav.nc, nous pouvons faire ceci :

gmt set PROJ_ELLIPSOID Sphère

gmt mapproject -R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | projet de carte gmt \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > ship.lonlatgravsat

Il est recommandé d'utiliser la méthode ci-dessus pour projeter et déprojeter les données dans un tel
une application, car il n'y a alors qu'une seule étape d'interpolation (en piste grd). Si un
essaie d'abord de convertir le fichier de grille en lon, lat, puis de l'échantillonner, il y a deux
étapes d'interpolation (en conversion et en échantillonnage).

Pour faire une grille lon, lat à partir de la grille ci-dessus, nous pouvons utiliser

gmt grdproject merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc

Dans certains cas, cela ne sera pas facile car le -R dans les deux systèmes de coordonnées peuvent ne pas s'aligner
bien. Lorsque cela se produit, nous pouvons également utiliser (en fait, il peut toujours être préférable d'utiliser)

gmt grd2xyz merc_grav.nc | projet de carte gmt \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
surface gmt -R-40/40/-70/70 -I2m -Ggrav.nc

Pour créer une carte Mercator de la région ci-dessus, supposons que notre valeur gmt.conf pour
PROJ_LENGTH_UNIT est en pouces. Ensuite, puisque le fichier merc_grav.nc ci-dessus est projeté avec -Jm1i
il fait 80 pouces de large. Nous pouvons faire une carte de 8 pouces de large en utilisant -Jx0.1i sur tous les programmes cartographiques
appliqué à cette grille (par exemple, grdcontour, grdimage, vue grd), puis pour les superpositions qui
travailler en lon,lat (par exemple, psxy, pcoast) nous pouvons utiliser le ci-dessus ajusté -R ainsi que -Jm0.1 pour obtenir
les deux systèmes pour correspondre.

Cependant, nous pouvons être plus intelligents que cela. Réalisant que le fichier img d'entrée avait des pixels 2.0
minutes de large (ou en vérifiant le nx et ny avec grdinfo merc_grav.nc) nous nous rendons compte que
merc_grav.nc a utilisé la pleine résolution du fichier img et il a 2400 par 2039 pixels, et
à 8 pouces de large, c'est 300 pixels par pouce. Nous décidons que nous n'en avons pas besoin d'autant et nous
sera satisfait de 100 pixels par pouce, nous voulons donc moyenner les données en 3 par 3
carrés. (Si nous voulons un tracé de contour, nous choisirons probablement de faire la moyenne des données beaucoup plus
(par exemple, 6 par 6) pour obtenir des contours lisses.) Puisque 2039 n'est pas divisible par 3, nous obtiendrons un
différent ajusté -R cette fois:

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V

Cette fois, nous trouvons que la région ajustée est -R-40/40/-70.023256525/-29.9368261101 et le
la sortie est de 800 par 601 pixels, une meilleure taille pour nous. Maintenant, nous pouvons créer un artificiel
fichier d'éclairage pour cela en utilisant dégradé:

gmt grdgradient merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

et si nous avons également un fichier CPT appelé "grav.cpt", nous pouvons créer une carte en relief ombrée de couleurs
comme ça:

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K > map.ps
gmt psbasemap -R-40/40/-70.023256525/-29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

Supposons que vous souhaitiez obtenir uniquement les valeurs de données contraintes d'un fichier img, en lat/lon
coordonnées. Puis cours img2grd les -T2 options, utiliser grd2xyz pour vider les valeurs, pipe
via grep -v NaN pour éliminer les NaN, et passer à travers projet de carte avec l'inverse
projection comme ci-dessus.

Utilisez img2grdgmt en ligne en utilisant les services onworks.net