grdtrackgmt – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

grdtrack – Beispielgitter an angegebenen (x,y)-Standorten

ZUSAMMENFASSUNG

grdspur [ xyfile ] grd1 grd2 ... [ f|p|m|r|R[+l] ] [ Länge[u]/ds[Abstand][+a] ] [dDatei
] [ Linie[,Linie,...][+aaz][+d][+iinc.[u]][+lLänge[u]][+nnp][+oaz][+rradius[u] ] [ ] [
Region ] [ Methode/Modifikatoren ] [ [radius[u]][+e|p]] [ [Grad des ] ] [ ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -n] [
-o] [ -s] [ -:[i|o]]

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

BESCHREIBUNG

grdspur liest eine oder mehrere Grid-Dateien (oder eine Sandwell/Smith-IMG-Datei) und eine Tabelle (aus
Datei oder Standardeingabe; aber siehe -E ausnahmsweise) mit (x,y) [oder (lon,lat)] Positionen in
die ersten beiden Spalten (es können weitere Spalten vorhanden sein). Es interpoliert das/die Gitter am
Positionen in der Tabelle und schreibt die Tabelle mit den hinzugefügten interpolierten Werten als (eins) aus
oder mehr) neue Spalten. Alternative (-C), wird die Eingabe als Liniensegmente betrachtet und
Wir erstellen orthogonale Kreuzprofile an jedem Datenpunkt oder mit äquidistantem Abstand
und probieren Sie das/die Gitter entlang dieser Profile aus. Ein bikubischer [Standard], bilinearer, B-Spline oder
nächster Nachbar (vgl -n) Interpolation verwendet wird, erfordert Randbedingungen an der
Grenzen der Region (vgl -n; Standardmäßig werden „natürliche“ Bedingungen verwendet (zweite partielle Ableitung).
(Normale zur Kante ist Null), es sei denn, das Gitter wird automatisch als periodisch erkannt.)

ERFORDERLICH ARGUMENTE

-GRasterdatei
grddatei ist eine 2D-Binärgitterdatei mit der Funktion f(x,y). Wenn das angegebene Raster
im Sandwell/Smith Mercator-Format vorliegen, müssen Sie eine durch Kommas getrennte Liste anhängen
Argumente, die eine Skala zum Multiplizieren der Daten (normalerweise 1 oder 0.1) und den Modus enthalten
die für Folgendes stehen: (0) IMG-Dateien ohne Einschränkungscode geben Daten zurück
an allen Punkten, (1) IMG-Datei mit codierten Einschränkungen, Rückgabedaten an allen Punkten, (2)
IMG-Datei mit codierten Einschränkungen, Rückgabe von Daten nur an eingeschränkten Punkten und NaN
an anderer Stelle, und (3) Img-Datei mit codierten Einschränkungen, 1 bei Einschränkungen und 0 zurückgeben
an anderer Stelle und optional den maximalen Breitengrad in der IMG-Datei [80.738]. Sie können es wiederholen
-G so oft, wie Sie Gitter haben, die Sie testen möchten. Alternativ verwenden Sie -G+lListe
um eine Liste von Dateinamen zu übergeben. Die Raster werden abgetastet und die Ergebnisse werden im ausgegeben
Auftrag gegeben. (Siehe GRID-DATEIFORMAT unten.)

OPTIONAL ARGUMENTE

xyfile Dies ist eine ASCII-Datei (oder Binärdatei, siehe -Bi) Datei, in der die ersten beiden Spalten die (x,y) enthalten
Positionen, an denen der Benutzer den 2D-Datensatz abtasten möchte.

-Af|pm|r|R
Für Spur-Resampling (wenn -C gesetzt ist) können wir auswählen, wie dies durchgeführt werden soll.
Anhängen f um die ursprünglichen Punkte beizubehalten, aber bei Bedarf Zwischenpunkte hinzuzufügen [Standard],
m as f, aber zuerst Meridian (entlang y) dann parallel (entlang x), p as f, Aber
Folgen Sie zuerst dem Parallelen (entlang y), dann dem Meridian (entlang x), r erneut abtasten
äquidistante Standorte; Eingabepunkte sind nicht unbedingt in der Ausgabe enthalten, und
R as r, aber passen Sie den angegebenen Abstand genau an die Schienenlänge an. Zum Schluss anhängen
+l wenn Entfernungen entlang der Loxodrome gemessen werden sollen. Wird ignoriert, es sei denn
-C wird eingesetzt.

-CLänge[u]/ds[Abstand][+a]
Verwenden Sie Eingabeliniensegmente, um eine äquidistante und (optional) gleichmäßig beabstandete Linie zu erstellen
Satz von Kreuzungsprofilen, entlang derer wir das/die Gitter abtasten. [Standardmäßig werden nur Stichproben durchgeführt
das/die Gitter an den Eingabeorten]. Geben Sie zwei Längenskalen an, die steuern, wie
Die Probenahme erfolgt: Länge legt die volle Länge jedes Querprofils fest, während ds
ist der Abtastabstand entlang jedes Querprofils. Optional anhängen /Abstand für
ein äquidistanter Abstand zwischen den Querprofilen [Standardmäßig werden Querprofile an der Stelle aufgestellt
Eingabekoordinaten]. Standardmäßig haben alle Querprofile die gleiche Richtung. Anhängen
+a um die Richtung der Querprofile zu ändern. Hängen Sie passende Einheiten an Länge;
Es legt die verwendete Einheit fest ds [und Abstand] (Siehe EINHEITEN unten). Die Standardeinheit für
Bei geographischen Gittern handelt es sich um Meter, bei kartesischen Gittern um die Benutzereinheit. Die Ausgabe
Spalten werden sein lon, lat, dist, Azimut, z1, z2... zn (Die zi sind die abgetasteten
Werte für jeden der n Gitter)

-DdDatei
Im Konzert mit -C Wir können die (möglicherweise neu abgetasteten) Originalzeilen in der Datei speichern
dDatei [Standardmäßig werden nur die Querprofile gespeichert]. Die Spalten werden sein lon, lat, dist,
Azimut, z1, z2, ... (abgetasteter Wert für jedes Gitter)

-ELinie[,Linie,...][+aaz][+d][+iinc.[u]][+lLänge[u]][+nnp][+oaz][+rradius[u]
Anstatt die eingegebenen Gleiskoordinaten zu lesen, geben Sie Profile über Koordinaten und an
Modifikatoren. Das Format jedes einzelnen Linie is Anfang/halt, Wobei Anfang or halt sind entweder
lon/lat (x/y für kartesische Daten) oder ein zweistelliger XY-Schlüssel, der das verwendet
Ausrichtungsformat im „pstext“-Stil, um einen Punkt auf der Karte anzugeben
[LCR][BMT]. Darüber hinaus können Sie Z-, Z+ verwenden, um das globale Minimum und Maximum zu bezeichnen
Positionen im Raster (nur verfügbar, wenn nur ein Raster angegeben ist). Statt zwei
Koordinaten können Sie einen Ursprung und einen davon angeben +a, +o, oder auch +r. Sie können anhängen
+iinc.[u], um das Abtastintervall festzulegen; Wenn nicht angegeben, verwenden wir standardmäßig die Hälfte
Mindestrasterabstand. Der +a legt den Azimut eines Profils gegebener Länge fest
beginnend am angegebenen Ursprung, while +o zentriert das Profil auf den Ursprung; beide
erfordern +l. Bei kreisförmiger Probenahme bitte angeben +r um einen Kreis mit gegebenem Radius zu definieren
zentriert auf den Ursprung; Diese Option erfordert beides +n or +idem „Vermischten Geschmack“. Seine +nnp setzt die
gewünschte Punktzahl, während +lLänge gibt die Gesamtlänge des Profils an.
Anhängen +d um die Streckenlängen nach den Koordinaten auszugeben. Hinweis: Nein
Die Trackdatei wird gelesen. Beachten Sie außerdem, dass nur eine Distanzeinheit ausgewählt werden kann.
Die Angabe unterschiedlicher Einheiten führt zu einem Fehler. Wenn keine Einheiten angegeben sind, wir
Standardmäßig werden Großkreisentfernungen in km verwendet (sofern geografisch). Wenn Sie mit arbeiten
geografische Daten, denen Sie - (Flacher Erdboden) oder + (Geodätisch) voranstellen können inc., Länge, oder auch
radius um den Modus der Entfernungsberechnung zu ändern [Großkreis].

-N Do nicht Überspringen Sie Punkte, die außerhalb des Bereichs des Gitters bzw. der Gitter liegen. [Nur Standardausgabe
Punkte innerhalb der Gitterdomäne].

-R[Einheit]xMin/xmax/ymin/ymax[r] (Mehr ...)
Geben Sie die Region von Interesse an.

-SMethode/Modifikatoren
In Verbindung mit -CBerechnen Sie ein einzelnes gestapeltes Profil aus allen Profilen
jedes Segment. Hängen Sie an, wie die Stapelung berechnet werden soll: a = Mittel (Durchschnitt), m =
Median, p = Modus (maximale Wahrscheinlichkeit), l = niedriger, L = niedriger, aber nur berücksichtigen
positive Werte, u = oben, U = oberer Wert, aber berücksichtigen Sie nur negative Werte [a]. Das
Modifikatoren Steuern Sie die Ausgabe. Wählen Sie eine oder mehrere dieser Möglichkeiten: +a : Anhängen
gestapelte Werte für alle Querprofile. +d : Stapelabweichungen an alle anhängen
Querprofile. +d : Datenreste (Daten – Stapel) an alle Querprofile anhängen.
+s[Datei]: Gestapeltes Profil speichern unter Datei [grdtrack_stacked_profile.txt]. +cTatsache :
Berechnen Sie die Hüllkurve des gestapelten Profils als +/- Tatsache *Abweichung [2]. Anmerkungen: (1)
Abweichungen hängen davon ab Methode und sind st.dev (a), L1-Skala (m und p) oder Halbbereich
(oben-unten)/2. (2) Die gestapelte Profildatei enthält 1 plus Gruppen von 4–6
Spalten, eine Gruppe für jedes Stichprobenraster. Die erste Spalte enthält den Querabstand,
während die ersten 4 in einer Gruppe gestapelte Werte, Abweichungen, Mindestwerte und Höchstwerte enthalten
Wert. Wenn Methode eines von a|m|p dann schreiben wir auch die untere und obere Konfidenz
Grenzen (vgl +c). Wenn einer oder mehrere von +a, +d und +r werden verwendet, dann hängen wir die an
Ergebnisse bis zum Ende jeder Zeile für alle Querprofile. Die Bestellung ist immer gestapelt
Wert (+a), gefolgt von Abweichungen (+d) und Residuen (+r). Wenn mehr als einer
Beim Abtasten des Rasters wird diese Folge von 1–3 Spalten für jedes Raster wiederholt.

-T[radius[u]][+e|p]
Zur Verwendung mit normalem Raster-Sampling und beschränkt auf ein einzelnes Nicht-IMG-Raster. Wenn
Der dem Eingabepunkt am nächsten gelegene Knoten ist NaN. Suchen Sie nach außen, bis wir ihn finden
Suchen Sie den nächstgelegenen Nicht-NaN-Knoten und melden Sie stattdessen diesen Wert. Geben Sie optional eine Suche an
Radius, der die Betrachtung auf Punkte innerhalb dieser Entfernung von der Eingabe beschränkt
Punkt. Um die Position des nächstgelegenen Knotens und seine Entfernung von der Eingabe zu melden
Punkt, anhängen +e. Um stattdessen den Eingabepunkt durch die Koordinaten des zu ersetzen
Nächster Knoten, anhängen +p.

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

-Z Schreiben Sie nur die abgetasteten Z-Werte aus [Standard schreibt alle Spalten].

-: Schaltet zwischen der Eingabe/Ausgabe (Längengrad, Breitengrad) und (Breitengrad, Längengrad) um.
[Standard ist (Längengrad, Breitengrad)].

-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 aus. [Standard ist eins mehr als die 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).

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

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

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

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

GRID FILE FORMATEN

Standardmäßig schreibt GMT das Raster als Floats mit einfacher Genauigkeit in einer COARDS-Beschwerde netCDF
Datei Format. GMT ist jedoch in der Lage, Rasterdateien in vielen anderen häufig verwendeten Raster zu erstellen
Dateiformate und erleichtert auch das sogenannte "Packen" von Rastern, das Ausschreiben von Gleitkomma
Daten als 1- oder 2-Byte-Ganzzahlen. Um Genauigkeit, Skalierung und Offset anzugeben, sollte der Benutzer
füge das Suffix hinzu =id[/Treppe/Offset[/nan]], wo id ist eine zweibuchstabige Kennung des Rasters
Art und Genauigkeit, und Treppe und Offset sind optional Skalierungsfaktor und Offset zu sein
auf alle Rasterwerte angewendet, und nan ist der Wert, der verwendet wird, um fehlende Daten anzuzeigen. Falls
die beiden Charaktere id ist nicht vorgesehen, wie in =/Treppe als ein id=nf wird angenommen. Wann
Leseraster wird das Format in der Regel automatisch erkannt. Wenn nicht, das gleiche Suffix
kann zu den Dateinamen des Eingaberasters hinzugefügt werden. Sehen grdconvert und Abschnitt Grid-Datei-Format der
GMT Technische Referenz und Kochbuch für weitere Informationen.

Beim Lesen einer netCDF-Datei, die mehrere Raster enthält, liest GMT standardmäßig die
das erste 2-dimensionale Raster, das in dieser Datei gefunden werden kann. Um GMT dazu zu bringen, eine andere zu lesen
mehrdimensionale Variable in der Rasterdatei, anhängen ?Varname zum Dateinamen, wobei
Varname ist der Name der Variablen. Beachten Sie, dass Sie möglicherweise der besonderen Bedeutung entkommen müssen
of ? in Ihrem Shell-Programm, indem Sie einen umgekehrten Schrägstrich davor setzen oder das
Dateiname und Suffix zwischen Anführungszeichen oder doppelten Anführungszeichen. Die ?Varname Suffix kann auch verwendet werden
für Ausgaberaster, um einen anderen Variablennamen als den Standard anzugeben: "z". Sehen
grdconvert und Abschnittsmodifikatoren-für-CF und Grid-Datei-Format des GMT Technical
Referenz- und Kochbuch für weitere Informationen, insbesondere zum Lesen von 3-,
4- oder 5-dimensionale Raster.

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 .

HINWEISE

Wenn sich ein Interpolationspunkt nicht auf einem Knoten des Eingabegitters befindet, dann ein NaN an einem beliebigen Knoten in
Die Umgebung um den Punkt ergibt ein interpoliertes NaN. Bikubisch
Interpolation [Standard] ergibt kontinuierliche erste Ableitungen, erfordert jedoch eine Umgebung von
4 Knoten mal 4 Knoten. Bilineare Interpolation [-n] verwendet nur eine 2 x 2-Nachbarschaft, aber
ergibt nur Kontinuität nullter Ordnung. Verwenden Sie bikubisch, wenn Glätte wichtig ist. Verwenden
bilinear, um die Ausbreitung von NaNs zu minimieren, oder weniger Schwelle.

Beispiele:

Zum Abtasten der Datei hawaii_topo.nc entlang der SEASAT-Strecke track_4.xyg (Eine ASCII-Tabelle).
Enthält Längengrad, Breitengrad und die von SEASAT abgeleitete Schwerkraft, gefolgt von einer Kopfzeile
aufzeichnen):

grdtrack track_4.xyg -Ghawaii_topo.nc -h > track_4.xygt

Zum Beispiel die Sandwell/Smith IMG-Formatdatei topo.8.2.img (2 Minuten vorhergesagte Bathymetrie).
auf einem Mercator-Gitter) und das Altersgitter age.3.2.nc von Muller et al. entlang der Längen- und Breitengradkoordinaten
in der Datei cruise_track.xy angegeben, versuchen Sie es

grdtrack cruise_track.xy -Gtopo.8.2.img,1,1 -Gage.3.2.nc > Depths-Age.d

Zum Beispiel die Datei grav.18.1.img im IMG-Format von Sandwell/Smith (1 Minute).
ein Mercator-Gitter) entlang 100 km langer Querprofile, die orthogonal zum Liniensegment sind
angegeben in der Datei track.xy, Errichtung von Querprofilen alle 25 km und Bemusterung des Rasters
Versuchen Sie es alle 3 km

grdtrack track.xy -Ggrav.18.1.img,0.1,1 -C100k/3/25 -Ar > xprofiles.txt

Um das Raster data.nc entlang einer Linie von der unteren linken zur oberen rechten Ecke abzutasten,
Versuchen Sie es mit einem Rasterabstand von 1 km und Ausgabeabständen

grdtrack -ELB/RT+i1k+d -Gdata.nc > Profiles.txt

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