mapprojectgmt - En ligne dans le Cloud

PROGRAMME:

Nom

mapproject - Effectuez des transformations de carte en avant et en arrière, des conversions de datum et de la géodésie

SYNOPSIS

projet de carte [ les tables ] paramètres région [ b|B|f|F|o|O[can0/lat0] ] [ [dx/dy] ] [ c|i|p ] [
[donnée] ] [ [unité] ] [ [x0/y0/][[+|-]unité][+|-] ] [ ] [ ligne.xy[/[+|-]unité][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [ [h]de[/à] ] [ [niveau] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o] ]

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

DESCRIPTION

projet de carte lit les positions (longitude, latitude) à partir de les tables [ou entrée standard] et
calcule les coordonnées (x,y) en utilisant la projection cartographique et les échelles spécifiées. En option, il
peut lire les positions (x,y) et calculer les valeurs (longitude, latitude) en faisant l'inverse
transformation. Ceci peut être utilisé pour transformer des points linéaires (x,y) obtenus en numérisant
une carte de projection connue aux coordonnées géographiques. Peut également calculer des distances le long
voie, à un point fixe ou à l'approche la plus proche d'une ligne. Enfin, peut être utilisé pour effectuer
diverses conversions de références. Des champs de données supplémentaires sont autorisés après les 2 premières colonnes
qui doit avoir (longitude,latitude) ou (x,y). Voir l'option -: sur la façon de lire
(latitude, longitude).

EST REQUIS ARGUMENTS

-Jparamètres (plus ...)
Sélectionnez la projection cartographique.

-R[unité]x min/xmax/ymin/ymax[r] (plus ...)
Spécifiez la région d'intérêt. Cas particulier de la projection UTM : Si -C est utilisé
ainsi que -R n'est pas donné, la région est définie pour coïncider avec la zone UTM donnée, donc
afin de préserver la solution ellipsoïdale complète (Voir RESTRICTIONS pour plus
information).

EN OPTION ARGUMENTS

table Un ou plusieurs ASCII (ou binaire, voir -bi[ncols][type]) fichier(s) de table de données contenant un
nombre de colonnes de données. Si aucun tableau n'est fourni, nous lisons à partir de l'entrée standard.

-Ab|B|f|F|o|O[can0/lat0]
-Un F calcule l'azimut (vers l'avant) à partir d'un point fixe long/lat à chaque point de données.
Utilisez -Un B pour obtenir l'azimut de retour des points de données au point fixe. Utilisation -Ao obtenir
orientations (-90/90) plutôt que des azimuts (0/360). Haut de casse F, B or O sera
convertir des latitudes géodésiques en latitudes géocentriques et estimer l'azimut des géodésiques
(en supposant que l'ellipsoïde actuel n'est pas une sphère). Si aucun point fixe n'est donné alors
nous calculons l'azimut (ou l'azimut arrière) à partir du point précédent.

-C[dx/dy]
Définissez le centre des coordonnées projetées sur le centre de projection de la carte [La valeur par défaut est
coin inférieur gauche]. En option, ajoutez des décalages dans les unités projetées à ajouter (ou
soustrait quand -I est défini) à (à partir de) les coordonnées projetées, telles que false
abscisse et ordonnée pour des zones de projection particulières [0/0]. L'unité utilisée pour le
offsets est l'unité de distance du tracé en vigueur (voir PROJ_LENGTH_UNIT) à moins que -F is
utilisé, auquel cas les décalages sont en mètres.

-Dc|je|p
Remplacer temporairement PROJ_LENGTH_UNIT et utiliser c (cm), i (pouce), ou p (points)
au lieu. Ne peut pas être utilisé avec -F.

-E[donnée]
Convertir de géodésique (lon, lat, hauteur) en Terre centrée sur la Terre fixe (ECEF)
(x,y,z) coordonnées (ajouter -I pour la conversion inverse). Ajouter un identifiant de référence (voir -Qd)
ou donner ellipsoïde:dx,dy,dz De ellipsoïde peut être un identifiant ellipsoïde (voir -Qe) ou
donné comme a[,*inv_f*], où a est le demi-grand axe et inv_f est l'inverse
aplatissement (0 si omis). Si donnée est - ou non donné, nous supposons WGS-84.

-F[unité]
Forcer la mise à l'échelle 1:1, c'est-à-dire la sortie (ou l'entrée, voir -I) les données sont en réel projeté
mètres. Pour spécifier d'autres unités, ajoutez l'unité souhaitée (voir UNITÉS). Sans -F,
la sortie (ou l'entrée, voir -I) sont dans les unités spécifiées par PROJ_LENGTH_UNIT (mais
sur le lien -D).

-G[x0/y0/][[+|-]unité][+|-]
Calculer les distances le long de la voie or au point facultatif réglé avec -Gx0/y0. Ajouter
l'unité de distance (voir UNITÉS), y compris c (distance cartésienne en utilisant l'entrée
coordonnées) ou C (distance cartésienne utilisant les coordonnées projetées). Les C unité
a besoin -R ainsi que -J à régler. Sans point fixe, nous calculons le cumul
distances le long de la voie. Ajouter - pour obtenir une distance incrémentielle entre les
points. Ajouter + pour spécifier le 2ème point via deux colonnes supplémentaires dans le fichier d'entrée.

-I Faire la transformation inverse, c'est-à-dire obtenir (longitude,latitude) à partir des données (x,y).

-Lligne.xy[/[+|-]unité][+]
Déterminez la distance la plus courte entre les points de données d'entrée et la ou les lignes indiquées dans
le fichier multisegment ASCII ligne.xy. La distance et les coordonnées du
le point le plus proche sera ajouté à la sortie sous forme de trois nouvelles colonnes. Ajoutez le
unité de distance (voir UNITÉS), y compris c (distance cartésienne utilisant les coordonnées d'entrée)
or C (distance cartésienne utilisant les coordonnées projetées). Les C l'unité nécessite -R ainsi que
-J à régler. Enfin, ajoutez + pour signaler l'identifiant du segment de ligne et la fraction
numéro de point au lieu de lon/lat du point le plus proche.

-N[a|c|g|m]
Convertir à partir des latitudes géodésiques (en utilisant l'ellipsoïde actuel ; voir PROJ_ELLIPSOID)
à l'une des quatre latitudes auxiliaires différentes (les longitudes ne sont pas affectées). Choisir
de authalique, conformel, géocentrique, et mlatitudes éridionales [géocentriques]. Utilisation -I
pour convertir des latitudes auxiliaires en latitudes géodésiques.

-Q[d|e Répertoriez tous les paramètres de projection. Pour lister uniquement les références, utilisez -Qd. Pour lister uniquement
ellipsoïdes, utilisez -Qe.

-S Supprimez les points qui tombent en dehors de la région.

-E]de[/à]
Coordonner les conversions entre les références de ainsi que à en utilisant le Molodensky standard
transformation. Utilisation -E si la 3ème colonne d'entrée a une hauteur au-dessus de l'ellipsoïde [Par défaut
suppose hauteur = 0, c'est-à-dire sur l'ellipsoïde]. Spécifiez les références à l'aide de l'ID de référence (voir
-Qd) ou donner ellipsoïde:dx,dy,dz De ellipsoïde peut être un identifiant ellipsoïde (voir -Qe) ou
donné comme a[,*inv_f*], où a est le demi-grand axe et inv_f est l'inverse
aplatissement (0 si omis). Si donnée est - ou non donné, nous supposons WGS-84. -T peut être
utilisé en conjonction avec -R -J pour changer le datum avant la projection des coordonnées
(ajouter -I pour appliquer la conversion de datum après la projection inverse). Sois sûr que
le paramètre PROJ_ELLIPSOID est correct pour votre cas.

-V[niveau] (plus ...)
Sélectionnez le niveau de verbosité [c].

-W[w|h]
Imprime la largeur et la hauteur de la carte sur la sortie standard. Aucun fichier d'entrée n'est lu. À seulement
afficher la largeur ou la hauteur, ajouter w or h, respectivement. Les unités de la
les dimensions peuvent être modifiées via -D.

-bi[ncols][t] (plus ...)
Sélectionnez l'entrée binaire native. [La valeur par défaut est de 2 colonnes d'entrée].

-bo[ncols][type] (plus ...)
Sélectionnez la sortie binaire native. [La valeur par défaut est identique à l'entrée].

-d[i|o]pas de données (plus ...)
Remplacer les colonnes d'entrée égales pas de données avec NaN et faire l'inverse en sortie.

-f[je|o]colinfo (plus ...)
Spécifiez les types de données des colonnes d'entrée et/ou de sortie.

-g[a]x|y|d|X|Y|D|[avec]z[+|-]écart[u] (plus ...)
Déterminez les lacunes dans les données et les sauts de ligne.

-h[je|o][n][+c][+d][+rremarque][+rtitre] (plus ...)
Ignorer ou produire des enregistrements d'en-tête.

-icols[l][sescaliers][ocompenser][, ] (plus ...)
Sélectionnez les colonnes d'entrée (0 est la première colonne).

-ocols[,...] (plus ...)
Sélectionnez les colonnes de sortie (0 est la première colonne).

-p[x|y|z]azim/élever[/niveau z][+wcan0/lat0[/z0]][+vx0/y0] (plus ...)
Sélectionnez la vue en perspective.

-s[cols][a|r] (plus ...)
Définir la gestion des enregistrements NaN.

- :[i|o] (plus ...)
Permutez la 1ère et la 2ème colonne en entrée et/ou en sortie.

-^ 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.

UNITÉS

Pour l'unité de distance de la carte, ajoutez unité d pour le degré d'arc, m pour la minute d'arc, et s pour arc
deuxième, ou e pour le compteur [Par défaut], f pour pied, k pour km, M pour mille terrestre, n pour nautique
mille, et u pour pied d'arpentage américain. Par défaut, nous calculons de telles distances en utilisant un
approximation avec les grands cercles. Préfixer - à une distance (ou l'unité est aucune distance n'est
donné) pour effectuer des calculs « Terre plate » (plus rapide mais moins précis) ou ajouter + à
effectuer des calculs géodésiques exacts (plus lents mais plus précis).

ASCII Format LA PRÉCISION

Les formats de sortie ASCII des données numériques sont contrôlés par les paramètres de votre gmt.conf
déposer. La longitude et la latitude sont formatées selon FORMAT_GEO_OUT, tandis que d'autres
les valeurs sont formatées selon FORMAT_FLOAT_OUT. Sachez que le format en vigueur peut
conduire à une perte de précision dans la sortie, ce qui peut conduire à divers problèmes en aval. Si
vous trouvez que la sortie n'est pas écrite avec suffisamment de précision, envisagez de passer en binaire
sortir (-bo si disponible) ou spécifiez plus de décimales à l'aide du paramètre FORMAT_FLOAT_OUT.

EXEMPLES

Pour convertir les coordonnées UTM en mètres en emplacements géographiques, étant donné un fichier utm.txt et
connaissant la zone UTM (et la zone ou l'hémisphère), essayez

gmt mapproject utm.txt -Ju+11/1:1 -C -I -F

Transformer un fichier avec (longitude,latitude) en (x,y) positions en cm sur un Mercator
grille pour une échelle donnée de 0.5 cm par degré, exécuter

gmt mapproject lonlatfile -R20/50/12/25 -Jm0.5c > xyfile

Pour transformer plusieurs fichiers 2 colonnes, binaires, double précision avec (latitude,longitude)
en positions (x,y) en pouces sur une grille de Mercator transverse (longitude centrale 75W) pour
échelle = 1:500000 et supprimez les points qui tomberaient en dehors de la zone de la carte, exécutez

gmt mapproject tracks.* -R-80/-70/20/40 -Jt-75/1:500000 -: -S -Di -bo -bi2 > tmfile.b

Pour convertir les coordonnées géodésiques (lon, lat, hauteur) dans le fichier old.dat du NAD27
Référence CONUS (ID de référence 131 qui utilise l'ellipsoïde Clarke-1866) à WGS 84, exécuter

gmt mapproject old.dat -Th131 > new.dat

Pour calculer la distance la plus proche (en km) entre chaque point dans le fichier d'entrée quakes.dat
et les segments de ligne donnés dans le fichier ASCII multisegment littoral.xy, exécutez

gmt mapproject quakes.dat -Lcoastline.xy/k > quake_dist.dat

RESTRICTIONS

La région d'entrée rectangulaire définie avec -R sera en général mappé dans un non-rectangulaire
la grille. Sauf si -C est défini, le point le plus à gauche sur cette grille a xvalue = 0.0, et le
le point le plus bas aura yvalue = 0.0. Ainsi, avant de numériser une carte, exécutez l'extrême
coordonnées cartographiques à travers projet de carte en utilisant l'échelle appropriée et voir quelles valeurs (x,y)
ils sont mappés dessus. Utilisez ces valeurs lors de la configuration de la numérisation afin d'avoir le
la transformation inverse fonctionne correctement, ou bien, utilisez awk pour mettre à l'échelle et déplacer le
(x,y) avant la transformation.

Pour certaines projections, une solution sphérique peut être utilisée bien que l'utilisateur ait sélectionné un
ellipsoïde. Cela se produit lorsque les utilisateurs -R paramètre implique une région qui dépasse le domaine
dans laquelle les développements en série ellipsoïdale sont valides. Ce sont les conditions : (1)
Conique conforme de Lambert (-JL) et Albers à aire égale (-JB) utilisera la solution sphérique
lorsque l'échelle de la carte dépasse 1.0E7. (2) Mercator transverse (-JT) et UTM (-JU) sera sera
utiliser la solution sphérique lorsque la limite ouest ou est donnée dans -R est plus que
10 degrés du méridien central, et (3) idem pour Cassini (-JC) mais avec une limite de
seulement 4 degrés.

ELLIPSODSDES ET SPHÉRDES

GMT utilisera des formules ellipsoïdales si elles sont implémentées et que l'utilisateur a sélectionné un
ellipsoïde comme forme de référence (voir PROJ_ELLIPSOID). L'utilisateur doit être conscient de quelques
pièges potentiels : (1) Pour certaines projections, telles que Transverse Mercator, Albers et
Conique conforme de Lambert, nous utilisons les expressions ellipsoïdales lorsque les zones cartographiées sont
petit, et passez aux expressions sphériques (et en remplaçant l'auxiliaire approprié
latitudes) pour les cartes plus grandes. Les formules ellipsoïdales sont utilisées comme suit : (a) Transverse
Mercator : lorsque tous les points sont à moins de 10 degrés du méridien central, (b) Projections coniques
lorsque la portée longitudinale est inférieure à 90 degrés, (c) projection de Cassini lorsque tous les points
sont à moins de 4 degrés du méridien central. (2) Lorsque vous essayez de faire correspondre des données historiques
données (par exemple, coordonnées obtenues avec une certaine projection et une certaine référence
ellipsoïde), vous pouvez constater que GMT donne des résultats légèrement différents. un probable
La source de cette inadéquation est que les calculs plus anciens utilisaient souvent des chiffres moins significatifs. Pour
exemple, les exemples de Snyder utilisent souvent l'ellipsoïde de Clarke 1866 (défini par lui comme ayant
un aplatissement f = 1/294.98). À partir de f, nous obtenons l'excentricité au carré égale à 0.00676862818
(c'est ce que GMT utilise), tandis que Snyder arrondit et utilise 0.00676866. Cette différence peut
donner des écarts de plusieurs dizaines de cm. Si vous avez besoin de reproduire des coordonnées projetées
avec cette excentricité légèrement différente, vous devez spécifier votre propre ellipsoïde avec le
mêmes paramètres que Clarke 1866, mais avec f = 1/294.97861076. Sachez également que les plus âgés
les données peuvent être référencées à des références différentes, et à moins que vous ne sachiez quelle donnée a été utilisée et
convertir toutes les données en une donnée commune, vous pouvez rencontrer des décalages de dizaines à des centaines de
mètres. (3) Enfin, sachez que PROJ_SCALE_FACTOR a certaines valeurs par défaut pour certains
projections, vous devrez donc peut-être remplacer le paramètre afin de faire correspondre les résultats produits
avec d'autres paramètres.

Utiliser mapprojectgmt en ligne à l'aide des services onworks.net