mapprojectgmt - Online in der Cloud

PROGRAMM:

NAME/FUNKTION

mapproject - Führen Sie vorwärts und inverse Kartentransformationen, Datumskonvertierungen und Geodäsie durch

ZUSAMMENFASSUNG

Kartenprojekt [ Tabellen ] Parameter Region [ b|B|f|F|o|O[can0/lat0] ] [ [dx/dy] ] [ c|i|p ] [
[Datum] ] [ [Einheit] ] [ [x0/y0/][[+|-]Einheit][+|-] ] [ ] [ Zeile.xy[/[+|-]Einheit][+] ] [
[a|c|g|m] ] [ [d|e ] [ ] [ [h]für[/zu] ] [ [Grad des ] ] [ [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

Hinweis: Zwischen dem Optionsflag und den zugehörigen Argumenten ist kein Leerzeichen zulässig.

BESCHREIBUNG

Kartenprojekt liest (Längengrad, Breitengrad) Positionen von Tabellen [oder Standardeingabe] und
berechnet (x,y)-Koordinaten mit der angegebenen Kartenprojektion und den angegebenen Maßstäben. Optional ist es
kann (x,y) Positionen lesen und (Längengrad, Breitengrad) Werte berechnen, indem man die Umkehrung macht
Transformation. Dies kann verwendet werden, um durch Digitalisierung erhaltene lineare (x,y) Punkte zu transformieren
eine Karte mit bekannter Projektion auf geographische Koordinaten. Kann auch Entfernungen entlang berechnen
Spur, zu einem festen Punkt oder der nächsten Annäherung an eine Linie. Schließlich kann verwendet werden, um durchzuführen
verschiedene Datumskonvertierungen. Nach den ersten 2 Spalten sind zusätzliche Datenfelder erlaubt
die haben muss (Längengrad, Breitengrad) oder (x,y). Siehe Option -: zum Lesen
(Breitengrad, Längengrad) Dateien.

ERFORDERLICH ARGUMENTE

-JParameter (Mehr ...)
Kartenprojektion auswählen.

-R[Einheit]xMin/xmax/ymin/ymax[r] (Mehr ...)
Geben Sie die Region von Interesse an. Sonderfall für die UTM-Projektion: Wenn -C verwendet wird
und -R nicht angegeben ist, wird die Region so eingestellt, dass sie mit der angegebenen UTM-Zone übereinstimmt
um die vollständige ellipsoide Lösung zu erhalten (Siehe EINSCHRÄNKUNGEN für mehr
Information).

OPTIONAL ARGUMENTE

Tabelle Ein oder mehrere ASCII (oder binär, siehe -Bi[ncols][tippe]) Datentabellendatei(en), die a . enthalten
Anzahl der Datenspalten. Wenn keine Tabellen angegeben sind, lesen wir von der Standardeingabe.

-Ab|B|f|F|o|O[can0/lat0]
-Af berechnet den (vorwärts gerichteten) Azimut vom Fixpunkt lang/lat zu jedem Datenpunkt.
Verwenden Sie die -Ab um den Azimut von Datenpunkten zum Fixpunkt zurückzubekommen. Verwenden -Ao bekommen
Ausrichtungen (-90/90) statt Azimute (0/360). Großbuchstaben F, B or O werden wir
Konvertieren von geodätischen in geozentrische Breiten und Schätzen des Azimuts der Geodäten
(vorausgesetzt, das aktuelle Ellipsoid ist keine Kugel). Wenn kein Fixpunkt angegeben ist, dann
Wir berechnen den Azimut (oder Back-Azimut) vom vorherigen Punkt.

-C[dx/dy]
Stellen Sie den Mittelpunkt der projizierten Koordinaten auf den Kartenprojektionspunkt ein [Standard ist
untere linke Ecke]. Fügen Sie optional Versätze in den projizierten Einheiten hinzu, die hinzugefügt werden sollen (oder
abgezogen, wenn -I gesetzt ist) zu (von) den projizierten Koordinaten, z. B. false
Rechts- und Hochwerte für bestimmte Projektionszonen [0/0]. Die Einheit für die
offsets ist die gültige Einheit der Plotentfernung (siehe PROJ_LENGTH_UNIT), es sei denn -F is
verwendet, wobei die Offsets in Metern angegeben sind.

-Dc|i|p
PROJ_LENGTH_UNIT vorübergehend überschreiben und verwenden c (cm), i (Zoll), oder p (Punkte)
stattdessen. Nicht verwendbar mit -F.

-E[Datum]
Konvertieren von geodätisch (lon, lat, height) in Earth Centered Earth Fixed (ECEF)
(x,y,z) Koordinaten (add -I für die Rückwandlung). Datums-ID anhängen (siehe -Qd)
oder gib Ellipsoid:dx,dy,dz woher Ellipsoid kann eine Ellipsoid-ID sein (siehe -Qeoder
gegeben als a[,*inv_f*], wobei a ist die große Halbachse und inv_f ist das Gegenteil
Abflachung (0 wenn weggelassen). Wenn Datum ist - oder nicht gegeben gehen wir von WGS-84 aus.

-F[Einheit]
1:1 Skalierung erzwingen, dh Ausgabe (oder Eingabe, siehe -I) Daten sind in der tatsächlichen Projektion
Meter. Um andere Einheiten anzugeben, hängen Sie die gewünschte Einheit an (siehe UNITS). Ohne -F,
die Ausgabe (oder Eingabe, siehe -I) sind in den von PROJ_LENGTH_UNIT angegebenen Einheiten (aber
sehen -D).

-G[x0/y0/][[+|-]Einheit][+|-]
Berechnen Sie Entfernungen entlang der Strecke or zum optionalen Punktsatz mit -Gx0/y0. Anhängen
die Entfernungseinheit (siehe UNITS), einschließlich c (Kartesischer Abstand mit Eingabe
Koordinaten) oder C (Kartesischer Abstand unter Verwendung projizierter Koordinaten). Die C Einheit
erfordert -R und -J eingestellt werden. Wenn kein Fixpunkt angegeben ist, berechnen wir kumulieren
Entfernungen entlang der Strecke. Anhängen - um einen inkrementellen Abstand zwischen aufeinanderfolgenden . zu erhalten
Punkte. Anhängen + um den 2. Punkt über zwei zusätzliche Spalten in der Eingabedatei anzugeben.

-I Führen Sie die inverse Transformation durch, dh erhalten Sie (Längengrad, Breitengrad) von (x,y)-Daten.

-LZeile.xy[/[+|-]Einheit][+]
Bestimmen Sie den kürzesten Abstand von den Eingabedatenpunkten zu den in angegebenen Linie(n).
die ASCII-Multisegment-Datei Zeile.xy. Die Entfernung und die Koordinaten der
Der nächste Punkt wird als drei neue Spalten an die Ausgabe angehängt. Anhängen der
Entfernungseinheit (siehe EINHEITEN), einschließlich c (Kartesischer Abstand mit Eingabekoordinaten)
or C (Kartesischer Abstand unter Verwendung projizierter Koordinaten). Die C Einheit erfordert -R und
-J eingestellt werden. Zum Schluss anhängen + um die Liniensegment-ID und den Bruch anzugeben
Punktnummer statt lon/lat des nächsten Punktes.

-N[a|c|g|m]
Von geodätischen Breitengraden umrechnen (unter Verwendung des aktuellen Ellipsoids; siehe PROJ_ELLIPSOID)
auf einen von vier verschiedenen zusätzlichen Breitengraden (Längengrade sind nicht betroffen). Wählen
für authalisch, cformell, geozentrisch und meridionale Breitengrade [geozentrisch]. Verwenden -I
von Hilfsbreiten in geodätische Breiten umzuwandeln.

-Q[d|e Listen Sie alle Projektionsparameter auf. Um nur Bezüge aufzulisten, verwenden Sie -Qd. Nur auflisten
Ellipsoide, verwenden -Qe.

-S Punkte unterdrücken, die außerhalb der Region liegen.

-NS]für[/zu]
Koordinatenumrechnungen zwischen Bezügen für und zu mit dem Standard Molodensky
Transformation. Verwenden -NS wenn die 3. Eingabespalte eine Höhe über dem Ellipsoid hat [Standard]
nimmt Höhe = 0 an, dh auf dem Ellipsoid]. Festlegen von Bezügen anhand der Bezugs-ID (siehe
-Qd) oder gib Ellipsoid:dx,dy,dz woher Ellipsoid kann eine Ellipsoid-ID sein (siehe -Qeoder
gegeben als a[,*inv_f*], wobei a ist die große Halbachse und inv_f ist das Gegenteil
Abflachung (0 wenn weggelassen). Wenn Datum ist - oder nicht gegeben gehen wir von WGS-84 aus. -T könnte sein
verwendet in Verbindung mit -R -J um das Datum vor der Koordinatenprojektion zu ändern
(hinzufügen -I um die Datumsumwandlung nach der inversen Projektion anzuwenden). Stelle sicher das
die PROJ_ELLIPSOID-Einstellung ist für Ihren Fall richtig.

-V[Grad des ] (Mehr ...)
Wählen Sie die Ausführlichkeitsstufe [c].

-W[b|h]
Druckt Kartenbreite und -höhe auf der Standardausgabe. Es werden keine Eingabedateien gelesen. Nur an
Breite oder Höhe ausgeben, anhängen w or h, bzw. Die Einheiten der
Abmessungen können geändert werden über -D.

-Bi[ncols][T] (Mehr ...)
Wählen Sie den nativen Binäreingang. [Standard sind 2 Eingabespalten].

-bo[ncols][tippe] (Mehr ...)
Wählen Sie die native Binärausgabe. [Standard ist gleich Eingabe].

-d[i|o]keine Daten (Mehr ...)
Ersetzen Sie Eingabespalten, die gleich sind keine Daten mit NaN und umgekehrt bei der Ausgabe.

-f[i|o]Colinfo (Mehr ...)
Geben Sie Datentypen von Eingabe- und/oder Ausgabespalten an.

-g[a]x|y|d|X|Y|D|[col]z[+|-]Lücke[u] (Mehr ...)
Bestimmen Sie Datenlücken und Zeilenumbrüche.

-h[i|o][n][+c][+d][+rAnmerkung][+rTitel] (Mehr ...)
Überspringen oder erzeugen Sie Kopfdatensätze.

-iSpalten[l][sTreppe][ÖOffset][,...] (Mehr ...)
Eingabespalten auswählen (0 ist die erste Spalte).

-oSpalten[,...] (Mehr ...)
Ausgabespalten auswählen (0 ist die erste Spalte).

-p[x|y|z]azim/Schüler[/zebene][+wcan0/lat0[/z0]][+vx0/y0] (Mehr ...)
Perspektivische Ansicht auswählen.

-S[Spalten][a|r] (Mehr ...)
Legen Sie die Behandlung von NaN-Datensätzen fest.

-:[i|o] (Mehr ...)
Vertauschen Sie die 1. und 2. Spalte bei der Eingabe und/oder Ausgabe.

-^ or nur -
Drucken Sie eine kurze Nachricht über die Syntax des Befehls und beenden Sie ihn (HINWEIS: unter Windows
benutze nur -).

-+ or nur +
Drucken Sie eine ausführliche Nutzungs-(Hilfe-)Nachricht, einschließlich der Erläuterungen zu allen
modulspezifische Option (aber nicht die allgemeinen GMT-Optionen), wird dann beendet.

-? or nicht Argumente
Drucken Sie eine vollständige Nutzungs-(Hilfe-)Nachricht, einschließlich der Erklärung der Optionen, dann
Ausgänge.

--Version
GMT-Version drucken und beenden.

--show-datadir
Vollständigen Pfad zum GMT-Freigabeverzeichnis drucken und beenden.

EINHEITEN

Für Kartenentfernungseinheit, anhängen Einheit d für Bogengrad, m für Bogenminute und s für Bogen
zweitens, oder e für Zähler [Standard], f für Fuß, k für km, M für Satzmeile, n für nautik
Meile, und u für US-Vermessungsfuß. Standardmäßig berechnen wir solche Abstände mit einem sphärischen
Näherung mit Großkreisen. Voranstellen - zu einer Entfernung (oder die Einheit ist keine Entfernung ist
gegeben), um "Flat Earth"-Berechnungen durchzuführen (schneller, aber weniger genau) oder voranzustellen + zu
genaue geodätische Berechnungen durchführen (langsamer, aber genauer).

ASCII FORMAT PRÄZISION

Die ASCII-Ausgabeformate numerischer Daten werden durch Parameter in Ihrem gmt.conf
Datei. Längen- und Breitengrad werden gemäß FORMAT_GEO_OUT formatiert, während andere
Werte werden nach FORMAT_FLOAT_OUT formatiert. Beachten Sie, dass das gültige Format
führen zu Präzisionsverlusten in der Ausgabe, die nachgelagert zu verschiedenen Problemen führen können. Wenn
Wenn Sie feststellen, dass die Ausgabe nicht mit ausreichender Genauigkeit geschrieben ist, sollten Sie auf Binär umschalten
Ausgang (-bo falls verfügbar) oder geben Sie weitere Dezimalstellen mit der Einstellung FORMAT_FLOAT_OUT an.

Beispiele:

Um UTM-Koordinaten in Metern in geografische Standorte umzuwandeln, geben Sie eine Datei utm.txt und
Wenn Sie die UTM-Zone (und Zone oder Hemisphäre) kennen, versuchen Sie es

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

Um eine Datei mit (Längengrad, Breitengrad) in (x,y) Positionen in cm auf einem Mercator umzuwandeln
Raster für einen gegebenen Maßstab von 0.5 cm pro Grad, Lauf

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

Um mehrere zweispaltige, binäre Dateien mit doppelter Genauigkeit zu transformieren mit (Breitengrad, Längengrad)
in (x,y) Positionen in Zoll auf einem Quer-Mercator-Gitter (mittlere Länge 75W) für
Maßstab = 1:500000 und unterdrücken Sie die Punkte, die außerhalb des Kartenbereichs liegen würden, ausführen

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

Um die geodätischen Koordinaten (lon, lat, height) in die Datei old.dat vom NAD27 umzuwandeln
CONUS-Datum (Datum ID 131, das das Clarke-1866-Ellipsoid verwendet) nach WGS 84, Lauf

gmt mapproject old.dat -Th131 > new.dat

Um die kürzeste Entfernung (in km) zwischen jedem Punkt in der Eingabedatei quakes.dat zu berechnen
und die in der Multisegment-ASCII-Datei Küstenlinie.xy angegebenen Liniensegmente, run

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

RESTRICTIONS

Der rechteckige Eingabebereich mit -R wird im Allgemeinen in ein nicht-rechteckiges
Netz. Wenn nicht -C gesetzt ist, hat der Punkt ganz links auf diesem Gitter xvalue = 0.0, und die
unterster Punkt hat yvalue = 0.0. Also, bevor Sie eine Karte digitalisieren, laufen Sie das Extrem
Kartenkoordinaten durch Kartenprojekt Verwenden Sie die entsprechende Skala und sehen Sie, welche (x,y)-Werte
sie sind abgebildet. Verwenden Sie diese Werte beim Einrichten der Digitalisierung, um die
inverse Transformation korrekt funktionieren oder alternativ verwenden awk skalieren und verschieben
(x,y)-Werte vor der Transformation.

Für einige Projektionen kann eine sphärische Lösung verwendet werden, obwohl der Benutzer eine gewählt hat
Ellipsoid. Dies geschieht, wenn die Benutzer -R Einstellung impliziert eine Region, die die Domäne überschreitet
in denen die ellipsoiden Reihenentwicklungen gültig sind. Dies sind die Bedingungen: (1)
Lambert-konformer Kegel (-J L) und Albers flächentreu (JB) verwendet die sphärische Lösung
wenn der Kartenmaßstab 1.0E7 überschreitet. (2) Quermercator (-JT) und UTM (-JU) wird
Verwenden Sie die sphärische Lösung, wenn entweder die West- oder Ostgrenze in angegeben ist -R ist mehr als
10 Grad vom Mittelmeridian und (3) dasselbe für Cassini (-JC) aber mit einer Grenze von
nur 4 Grad.

ELLIPSOIDE UND SPHEROIDE

GMT verwendet ellipsoide Formeln, wenn sie implementiert sind und der Benutzer eine . ausgewählt hat
Ellipsoid als Referenzform (siehe PROJ_ELLIPSOID). Der Benutzer muss einiges beachten
potenzielle Fallstricke: (1) Für einige Projektionen wie Transverse Mercator, Albers und
Lamberts konformer Kegelschnitt verwenden wir die Ellipsoidausdrücke, wenn die abgebildeten Flächen
klein, und wechseln Sie zu den sphärischen Ausdrücken (und ersetzen Sie die entsprechende Hilfsfunktion
Breitengrade) für größere Karten. Die Ellipsoidformeln werden wie folgt verwendet: (a) Quer
Mercator: Wenn alle Punkte innerhalb von 10 Grad vom Mittelmeridian liegen, (b) Kegelprojektionen
wenn der Längsbereich weniger als 90 Grad beträgt, (c) Cassini-Projektion, wenn alle Punkte
innerhalb von 4 Grad vom Mittelmeridian liegen. (2) Wenn Sie versuchen, einige historische
Daten (z. B. Koordinaten, die mit einer bestimmten Projektion und einer bestimmten Referenz erhalten wurden)
Ellipsoid) stellen Sie möglicherweise fest, dass GMT leicht abweichende Ergebnisse liefert. Eine wahrscheinlich
Ursache für diese Diskrepanz ist, dass bei älteren Berechnungen oft weniger signifikante Ziffern verwendet wurden. Zum
Zum Beispiel verwenden Snyders Beispiele oft das Ellipsoid von Clarke 1866 (von ihm definiert als mit
eine Abflachung f = 1/294.98). Aus f erhalten wir das Quadrat der Exzentrizität zu 0.00676862818
(dies wird von GMT verwendet), während Snyder abrundet und 0.00676866 verwendet. Dieser Unterschied kann
Abweichungen von mehreren zehn cm ergeben. Wenn Sie projizierte Koordinaten reproduzieren müssen
bei dieser etwas anderen Exzentrizität sollten Sie Ihr eigenes Ellipsoid mit dem
gleiche Parameter wie Clarke 1866, jedoch mit f = 1/294.97861076. Beachten Sie auch, dass ältere
Daten können auf unterschiedliche Datumsangaben bezogen sein und es sei denn, Sie wissen, welches Datum verwendet wurde und
Konvertieren aller Daten in ein gemeinsames Datum kann es zu Abweichungen von einigen Dutzend bis Hunderten kommen
Meter. (3) Beachten Sie schließlich, dass PROJ_SCALE_FACTOR für einige bestimmte Standardwerte hat
Projektionen, sodass Sie die Einstellung möglicherweise überschreiben müssen, um die erzeugten Ergebnisse abzugleichen
mit anderen Einstellungen.

Verwenden Sie mapprojectgmt online mit den onworks.net-Diensten