grdrotatergmt – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

grdrotater – Endliche Rotationsrekonstruktion des geografischen Gitters

ZUSAMMENFASSUNG

grdrotieren ingrdfile rot_file|lon/lat/Winkel outgrdfile [ rotoutline ] [ Polygondatei ] [ ]
[ Region ] [ ] [ Alter ] [ [Grad des ] ] [ -b] [ -d] [ -h] [
-i] [ -n] [ -:[i|o]]

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

BESCHREIBUNG

Grdrotater liest ein geografisches Raster und rekonstruiert es anhand einer vollständigen Rekonstruktion
Rotationen. Optional kann der Benutzer ein Clipping-Polygon im Mehrsegmentformat bereitstellen;
Dann wird nur der Teil des Gitters innerhalb des Polygons zur Bestimmung der Rekonstruktion verwendet
Region. Die Umrisse des rekonstruierten Bereichs werden ebenfalls zurückgegeben, sofern sie gedreht wurden
Region ist nicht der ganze Globus.

ERFORDERLICH ARGUMENTE

ingrdfile
Name einer Gitterdatei in geografischen Koordinaten (Längen- und Breitengrad).

-ERotfile
Geben Sie eine Datei mit Rotationsparametern an. Diese Datei muss jeweils einen Datensatz enthalten
Drehung; jeder Datensatz muss das folgende Format haben:

lon lat tbeginnen [stopp] Winkel [ Khat a b c d e f g df ]

woher tbeginnen und tstop sind in Myr und lon lat Winkel sind in Grad. tbeginnen und
tstop sind die Altersstufen des alten und des jungen Endes einer Etappe. Wenn tstop ist nicht vorhanden in
der Datensatz dann eine totale Rekonstruktionsrotation erwartet wird und tstop ist implizit
auf 0 gesetzt und sollte für keinen der Datensätze in der Datei angegeben werden. Wenn eine
Kovarianzmatrix C für die Rotation verfügbar ist, muss sie in einem Format angegeben werden
unter Verwendung der neun optionalen Begriffe, die in Klammern aufgeführt sind. Hier, C = (g/Khat)*[ a b d; b c
e; d e f ] welche Shows C bestehend aus drei Zeilenvektoren. Wenn die Freiheitsgrade
(df) bei der Anpassung ist die Drehung 0 oder nicht, wenn sie auf 10000 gesetzt ist. Leerzeilen und
Datensätze, deren erste Spalte # enthält, werden ignoriert. Sie können ein führendes + voranstellen
an den Dateinamen an, um anzugeben, dass Sie die Drehungen umkehren möchten. Alternative 1: Geben
der Dateiname besteht aus zwei Platten-IDs, die durch einen Bindestrich getrennt sind (z. B. PAC-MBL) und we
wird diese Rotation stattdessen aus der GPlates-Rotationsdatenbank extrahieren. Wir kehren zurück und
Fehler, wenn die Rotation nicht gefunden werden kann. Alternative 2: Angeben lon/lat/WinkelDh
der Längengrad, der Breitengrad und der Öffnungswinkel (alle in Grad und durch / getrennt) für
eine einzige vollständige Rekonstruktionsrotation.

-Goutgrdfile
Name des Ausgabegitters. Dies ist das Raster mit den entsprechend rekonstruierten Daten
angegebene Drehung. Wenn mehr als eine Rekonstruktionszeit impliziert ist, dann outgrdfile
muss einen C-Format-Bezeichner enthalten, um eine Gleitkommazahl zu formatieren (Rekonstruktion).
Zeit) zum Texten.

OPTIONAL ARGUMENTE

-Drotoutline
Name der Gitterpolygon-Umrissdatei. Dies stellt den Umriss des Rasters dar
auf den angegebenen Zeitpunkt rekonstruiert. Wenn mehr als eine Rekonstruktionszeit vorhanden ist
dann impliziert rotoutline muss einen C-Formatbezeichner enthalten, um ein Floating zu formatieren
Punktnummer (Rekonstruktionszeit) zum Text. Wenn nur einmal impliziert ist und -D is
nicht gesetzt, dann schreiben wir das Polygon nach stdout (aber siehe -N).

-FPolygondatei
Geben Sie eine geschlossene Polygondatei mit mehreren Segmenten an, die den Innenbereich des beschreibt
Raster, das projiziert werden soll [Standard projiziert das gesamte Raster].

-N Den gedrehten Polygonumriss nicht ausgeben [Standardmäßig wird er in stdout oder in a geschrieben
Datei über -D].

-R[Einheit]Westen/Osten/Süden/Norden[/zmin/zmax][R]
Westen, Osten, Süden und Norden Geben Sie die Region von Interesse an, und Sie können angeben
sie in Dezimalgrad oder im Format [+-]dd:mm[:ss.xxx][W|E|S|N]. Anhängen r wenn niedriger
Anstelle von w/e/s/n werden die linken und oberen rechten Kartenkoordinaten angegeben. Die Zwei
Abkürzungen -Rg und -Rd stehen für globale Domäne (0/360 und -180/+180 in Längengrad
jeweils mit -90/+90 Breitengrad). Geben Sie alternativ für die Rastererstellung
RCodelon/lat/nx/ny, Wobei Code ist eine 2-Zeichen-Kombination aus L, C, R (für links,
Mitte oder rechts) und T, M, B für oben, Mitte oder unten. zB BL für unten links.
Dies gibt an, an welchem ​​Punkt einer rechteckigen Region die lon/lat Koordinate bezieht sich
zu, und die Rastermaße nx und ny mit Rasterabständen über -I wird verwendet, um zu erstellen
die entsprechende Region. Geben Sie alternativ den Namen einer vorhandenen Rasterdatei an
und den -R Einstellungen (und ggf. Rasterabstand) werden aus dem Raster kopiert.
Die richtigen -REinheit erwartet projizierte (kartesische) Koordinaten, die mit ausgewählten kompatibel sind -J
und wir projizieren umgekehrt, um die tatsächliche rechteckige geographische Region zu bestimmen. Zum
perspektivische Ansicht (-p), optional anhängen /zmin/zmax. Bei perspektivischer Ansicht
(-p), ein Z-Bereich (zmin, zmax) kann angehängt werden, um die dritte Dimension anzugeben. Dies
muss nur bei Verwendung des -Jz Option, nicht wenn nur die -p .
Im letzteren Fall wird eine perspektivische Ansicht der Ebene ohne Drittel aufgetragen
Abmessungen.

-S Überspringen Sie die Drehung des Gitters, drehen Sie einfach den Polygonumriss (erfordert -F wenn nein
Raster ist vorhanden).

-TAlter Legt die gewünschten Rekonstruktionszeiten fest. Fügen Sie für ein einzelnes Mal die gewünschte Zeit hinzu.
Geben Sie für einen äquidistanten Bereich der Rekonstruktionszeiten an -TAnfang/halt/inc. or
-TAnfang/halt/npunkte+. Für einen nicht äquidistanten Satz von Rekonstruktionszeiten bitte
Übergeben Sie sie über die erste Spalte in einer Datei, z. B. -TAltersdatei. Wenn nein -T Option ist
gegeben und -E Wenn wir eine Rotationsdatei angegeben haben, setzen wir die Rotationsdateizeiten mit gleich
die Wiederaufbauzeiten.

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

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

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

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

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

-n[b|c|l|n][+a][+bBC][+c][+tSchwelle] (Mehr ...)
Wählen Sie den Interpolationsmodus für Gitter.

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

FOLGEN OF GRID RESAMPLING

Beim erneuten Abtasten oder Abtasten von Gittern werden verschiedene Algorithmen verwendet (siehe -n) das kann zu
mögliche Verzerrungen oder unerwartete Ergebnisse in den neu abgetasteten Werten. Ein erwarteter Effekt
beim Resampling mit Splines ist die Tendenz, dass die neuen Resampling-Werte leicht über
die globalen Min/Max-Grenzen des ursprünglichen Rasters. Wenn dies nicht akzeptabel ist, können Sie auferlegen
Clipping der neu abgetasteten Wertewerte, damit sie die eingegebenen Min/Max-Werte nicht um . überschreiten
Hinzufügen +c auf Ihre -n .

GEODETIK VERSUS GEOZENTRISCH KOORDINIEREN

Alle sphärischen Drehungen werden auf geozentrische Koordinaten angewendet. Dies bedeutet, dass eingehende
Datenpunkte und Gitter werden als geodätische Koordinaten angesehen und müssen zuerst
in geozentrische Koordinaten umgewandelt. Dann werden Drehungen angewendet, und die endgültige
rekonstruierte Punkte werden wieder in geodätische Koordinaten umgewandelt. Dieses Standardverhalten
kann umgangen werden, wenn die Ellipsoideinstellung PROJ_ELLIPSOID auf Kugel geändert wird.

Beispiele:

Um die durch das Raster topo.nc und den Polygonumriss clip_path.d definierten Daten zu drehen, verwenden Sie a
Gesamtrekonstruktionsrotation mit Pol bei (135.5, -33.0) und einem Rotationswinkel von 37.3
Grad und bikubische Interpolation, versuchen Sie es

gmt grdrotater topo.nc -E135.5/-33/37.3 -V -Fclip_path.d -Grot_topo.nc > rot_clip_path.d

Um das gesamte Gitter faa.nc mithilfe der Rotationsdatei wieder auf 32 Ma zu drehen rotations.txt und eine
Bilineare Interpolation, versuchen Sie es

gmt grdrotater faa.nc -Erotations.txt -T32 -V -Grot_faa.nc -nl > rot_faa_path.d

Versuchen Sie es einfach, um zu sehen, wie sich der Umriss des Gitters large.nc nach derselben Drehung darstellt

gmt grdrotater large.nc -Erotations.txt -T32 -V -S \| psxy -Rg -JH180/6i -B30 -W0.5p \| gv -

Um das Gitter topo.nc mithilfe der Rotationsdatei wieder auf 100 Ma zu drehen rotations.txt und
Fordern Sie alle 10 Millionen Jahre eine Rekonstruktion an, wobei sowohl Raster als auch Umrisse unter Dateinamen gespeichert werden
Von Vorlagen ableiten, versuchen Sie es

gmt grdrotater topo.nc -Erotations.txt -T10/100/10 -V -Grot_topo_%g.nc -Drot_topo_path_%gd

Angenommen, Sie haben „gridA.nc“ und „gridB.nc“ gedreht und jede Drehung auf die darin enthaltenen Knoten beschränkt
Polygone polyA.d bzw. polyB.d unter Verwendung der Rotation A = (123W,22S,16,4) und
Drehung B = (108W, 16S, -14.5), was die gedrehten Gitter rot_gridA.nc und rot_gridB.nc ergibt. Zu
Um den Überlappungsbereich zwischen den gedrehten Gittern zu bestimmen, verwenden wir grdmath:

gmt grdmath 1 rot_gridA.nc ISNAN SUB 1 rot_gridB.nc ISNAN SUB 2 EQ = Overlap.nc

Das Raster „overlap.nc“ hat jetzt 1 in den Überlappungsbereichen und 0 an anderen Stellen. Du kannst es benutzen
als Maske oder verwenden Sie grdcontour -D um ein Polygon (dh eine Kontur) zu extrahieren.

ANMERKUNG

GMT vertreibt das EarthByte-Rotationsmodell Global_EarthByte_230-0Ma_GK07_AREPS.rot. Zu
Verwenden Sie eine alternative Rotationsdatei, erstellen Sie einen Umgebungsparameter namens GPLATES_ROTATIONS
die auf eine alternative Rotationsdatei verweist.

Verwenden Sie grdrotatergmt online über die Dienste von onworks.net