grdfftgmt – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

grdfft – Führen Sie mathematische Operationen an Gittern im Wellenzahl- (oder Frequenz-)Bereich durch

ZUSAMMENFASSUNG

grdfft ingrid [ ingrid2 ] Outfile [ Azimut ] [ zebene ] [ [Treppe|g] ] [ [r|x|y][w[k]] ] [
[r|x|y]params ] [ [Treppe|g] ] [
[f|q|s|nx/ny][+a|d|h|l][+e|n|m][+tBreite][+w[Suffix]][+z[p]] ] [ Treppe ] [ [Grad des ] ] [ -fg
]

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

BESCHREIBUNG

grdfft nimmt die 2-D-Vorwärts-Fast-Fourier-Transformation und führt eine oder mehrere durch
mathematische Operationen im Frequenzbereich vor der Rücktransformation in den Raum
Domain. Es besteht die Möglichkeit, die Daten zu skalieren, bevor die neuen Werte in einen Ausgang geschrieben werden
Datei. Die horizontalen Abmessungen des Gitters werden in Metern angenommen. Geografisch
Raster können durch Angabe von verwendet werden -fg Option, die Grad in Meter skaliert. Wenn Sie haben
Gitter mit Abmessungen in km, Sie können dies mit in Meter ändern gredit oder skalieren Sie die
Ausgabe mit grdmath.

ERFORDERLICH ARGUMENTE

ingrid 2D-Binärgitterdatei, die bearbeitet werden soll. (Siehe GRID-DATEIFORMATE unten). Für
Kreuzspektraloperationen geben auch die zweite Gitterdatei an ingrd2.

-GOutfile
Geben Sie den Namen der Ausgabegitterdatei oder der 1D-Spektrumtabelle an (siehe -E). (Sehen
Rasterdateiformate unten).

OPTIONAL ARGUMENTE

-AAzimut
Nehmen Sie die Richtungsableitung im Azimut Richtung gemessen in Grad CW
von Norden.

-Czebene
Aufwärts (für zebene > 0) oder nach unten (z zebene < 0) Feld fortsetzen zebene
Meter.

-D[Treppe|g]
Differenzieren Sie das Feld, dh nehmen Sie d(Feld)/dz. Dies entspricht einer Multiplikation
durch kr im Frequenzbereich (kr ist die Radialwellenzahl). Hängen Sie eine Skala an
multiplizieren mit (kr * Treppe) stattdessen. Alternativ anhängen g um anzuzeigen, dass Ihr
Die Daten sind Geoidhöhen in Metern und die Ausgabe sollte Schwerkraftanomalien in mGal sein.
[Standard ist keine Skala].

-E[r|x|y][w[k]]
Schätzen Sie das Leistungsspektrum in radialer Richtung [r]. Ort x or y gleich nach
-E stattdessen das Spektrum in x- oder y-Richtung berechnen. Es gibt keine Rasterdatei
erstellt. Wenn ein Gitter gegeben ist, dann sind f (d. h. Frequenz oder Wellenzahl), Leistung[f],
und 1 Standardabweichung in Potenz[f] werden in die von eingestellte Datei geschrieben -G [Standardausgabe]. Wenn
Es sind zwei Gitter gegeben, wir schreiben f und 8 Größen: Xpower[f], Ypower[f], kohärent
Leistung[f], Rauschleistung[f], Phase[f], Admittanz[f], Verstärkung[f], Kohärenz[f]. Jede
Auf die Menge folgt eine eigene 1-Standard-Entwicklungsfehlerschätzung, daher ist die Ausgabe 17
Spalten breit. Anhängen w statt Frequenz Wellenlänge zu schreiben. Wenn Ihr Raster ist
Geografisch können Sie weitere anhängen k um Wellenlängen von Meter [Standard] auf zu skalieren
km.

-F[r|x|y]params
Filtern Sie die Daten. Ort x or y gleich nach -F filtern x or y nur Richtung;
Der Standardwert ist isotrop [r]. Wählen Sie zwischen einem kosinusförmig verjüngten Bandpass und einem Gaußschen
Bandpassfilter oder ein Butterworth-Bandpassfilter.

Kosinuskegel:
Geben Sie vier Wellenlängen an lc/lp/hp/hc in den richtigen Einheiten (siehe -fg) zu entwerfen a
Bandpassfilter: Wellenlängen größer als lc oder weniger als hc wird geschnitten,
Wellenlängen größer als lp und weniger als hp durchgelassen werden, und Wellenlängen
dazwischen ist kosinusförmig. Z.B, -F1000000 / 250000 / 50000 / 10000 -fg
wird Bandpass und schneidet Wellenlängen > 1000 km und < 10 km durch
Wellenlängen zwischen 250 km und 50 km. Um einen Hochpass- oder Tiefpassfilter zu erstellen,
Geben Sie Bindestriche (-) für ein hp/hc or lc/lp. Z.B, -Fx-/-/50/10 wird Tiefpass sein x,
Durchlassen von Wellenlängen > 50 und Zurückweisen von Wellenlängen < 10. -Fy1000/250 / - / -
wird Hochpass y, lässt Wellenlängen < 250 durch und unterdrückt Wellenlängen > 1000.

Gauß Bandpass:
Anhängen lo/hi, die beiden Wellenlängen in korrekten Einheiten (siehe -fg) zu entwerfen a
Bandpassfilter. Bei den gegebenen Wellenlängen gelten die Gewichte des Gaußschen Filters
0.5 sein. Um einen Hochpass- oder Tiefpassfilter zu erstellen, geben Sie einen Bindestrich (-) für ein hi
or lo Wellenlänge bzw. Z.B, -F-/30 wird die Daten mit a tiefpassen
Gaußfilter mit Halbgewicht bei 30, während -F400/- wird das übertreffen
Daten.

Butterworth Bandpass:
Anhängen lo/hi/Auftrag, die beiden Wellenlängen in korrekten Einheiten (siehe -fg) und das
Filterreihenfolge (eine Ganzzahl), um einen Bandpassfilter zu entwerfen. Zum gegebenen Zeitpunkt
Wellenlängen betragen die Butterworth-Filtergewichte 0.5. Einen Hochpass machen
oder Tiefpassfilter, geben Sie einen Bindestrich (-) für das ein hi or lo Wellenlänge,
bzw. Z.B, -F-/30/2 führt einen Tiefpass für die Daten mit 2. Ordnung durch
Butterworth-Filter, mit halbem Gewicht bei 30, während -F400/-/2 wird Hochpass sein
die Daten.

-ICH[Treppe|g]
Integrieren Sie das Feld, dh berechnen Sie integral_over_z (Feld * dz). Das ist
entspricht der Division durch kr im Frequenzbereich (kr ist die Radialwellenzahl).
Fügen Sie eine Skala hinzu, um durch (kr * zu dividieren) Treppe) stattdessen. Alternativ anhängen g zu
Geben Sie an, dass es sich bei Ihrem Datensatz um Schwerkraftanomalien in mGal handelt und dass die Ausgabe geoid sein sollte
Höhen in Metern. [Standard ist keine Skala].

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tBreite][+w[Suffix]][+z[p]]
Wählen oder erkundigen Sie sich nach geeigneten Rastermaßen für FFT und stellen Sie sie optional ein
Parameter. Steuern Sie die FFT-Dimension:
-Nf zwingt die FFT, die tatsächlichen Dimensionen der Daten zu verwenden.

-Nq Ich werde mich nach geeigneteren Dimensionen erkundigen, diese melden und dann fortfahren.

-Ns zeigt eine Liste optionaler Dimensionen an und wird dann beendet.

-Nnx/ny führt eine FFT für die Array-Größe durch nx/ny (muss >= Rasterdateigröße sein). Standard
wählt Dimensionen >= Daten, die die Geschwindigkeit und Genauigkeit der FFT optimieren. Wenn FFT
Abmessungen > Rasterdateiabmessungen, Daten werden erweitert und auf Null reduziert.

Steuern Sie die Trendbeseitigung von Daten: Fügen Sie Modifikatoren zum Entfernen eines linearen Trends hinzu:
+d: Trenddaten entfernen, dh den am besten passenden linearen Trend entfernen [Standard].

+a: Nur Mittelwert entfernen.

+h: Nur Mittelwert entfernen, also 0.5 * (Max + Min).

+l: Daten in Ruhe lassen.

Steuern Sie die Erweiterung und Verjüngung der Daten: Verwenden Sie Modifikatoren, um die Erweiterung zu steuern
und Tapering sind durchzuführen:
+e erweitert das Raster durch Auferlegen einer Kanten-Punkt-Symmetrie [Standard],

+m Erweitert das Raster durch Auferlegen einer Kantenspiegelsymmetrie

+n deaktiviert die Datenerweiterung.

Die Verjüngung erfolgt vom Datenrand zum FFT-Gitterrand [100 %]. Ändern
dieser Prozentsatz über +tBreite. Wenn +n wirksam ist, wird die Verjüngung angewendet
stattdessen an den Datenrändern, da keine Erweiterung verfügbar ist [0 %].

Kontrollieren Sie das Schreiben vorläufiger Ergebnisse: Für eine detaillierte Untersuchung können Sie das schreiben
Zwischengitter wird an die Vorwärts-FFT weitergeleitet; das dürfte so gewesen sein
trendbefreit, durch Punktsymmetrie entlang aller Kanten verlängert und verjüngt. Anhängen
+w[Suffix], aus dem/denen Ausgabedateinamen erstellt werden (d. h. ingrid_prefix.ext)
[konisch], wo ext ist Ihre Dateierweiterung. Abschließend können Sie das komplexe Raster speichern
erzeugt durch die Vorwärts-FFT durch Anhängen +z. Standardmäßig schreiben wir das echte und
imaginäre Komponenten zu ingrid_real.ext und ingrid_image.ext. Anhängen p , um zu speichern
stattdessen die polare Form von Betrag und Phase zu Dateien ingrid_mag.ext und
ingrid_Phase.ext.

-STreppe
Multiplizieren Sie jedes Element mit Treppe im Raumbereich (nach dem Frequenzbereich).
Operationen). [Standard ist 1.0].

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

-fg Geografische Raster (Längen- und Breitengrade) werden in Meter umgerechnet
über eine "Flat Earth"-Näherung unter Verwendung der aktuellen Ellipsoidparameter.

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

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.

GRID DISTANCE EINHEITEN

Wenn das Raster kein Meter als horizontale Einheit hat, fügen Sie +uEinheit zur Eingabedatei
Name, der von der angegebenen Einheit in Meter umgerechnet werden soll. Wenn Ihr Raster geografisch ist, konvertieren Sie
Entfernungen zu Metern durch Anlieferung -fg stattdessen.

ÜBERLEGUNGEN

netCDF COARDS-Gitter werden automatisch als geografisch erkannt. Für andere Gitter
Wählen Sie die geografischen Gitter aus, in denen Sie Grad in Meter umrechnen möchten -fg. Wenn die Daten
sich in der Nähe eines Pols befinden, sollten Sie in Betracht ziehen, die Rasterdatei auf ein Rechteck zu projizieren
Koordinatensystem mit grdprojekt

Beispiele:

Um nach oben zu gelangen, werden die magnetischen Anomalien auf Meereshöhe in der Datei mag_0.nc bis zu einer Höhe von 800 m fortgesetzt
Über dem Meeresspiegel:

gmt grdfft mag_0.nc -C800 -V -Gmag_800.nc

Um Geoidhöhen in m (geoid.nc) auf einem geografischen Gitter in die Schwerkraft in freier Luft umzuwandeln
Anomalien in mGal:

gmt grdfft geoid.nc -Dg -V -Ggrav.nc

Um Schwerkraftanomalien in mGal (faa.nc) in Abweichungen der Vertikalen (in) umzuwandeln
Mikroradiant) in der 038-Richtung müssen wir dann zuerst die Schwerkraft integrieren, um das Geoid zu erhalten
Nehmen Sie die Richtungsableitung und skalieren Sie schließlich das Bogenmaß in Mikro-Bogenmaß:

gmt grdfft faa.nc -Ig -A38 -S1e6 -V -Gdefl_38.nc

Zweite vertikale Ableitungen von Schwerkraftanomalien hängen mit der Krümmung der zusammen
Feld. Wir können diese als mGal/m^2 berechnen, indem wir zweimal differenzieren:

gmt grdfft Gravity.nc -D -D -V -Ggrav_2nd_derivative.nc

Um kreuzspektrale Schätzungen für gemeinsam registrierte Bathymetrie- und Schwerkraftgitter zu berechnen, und
Geben Sie das Ergebnis als Funktion der Wellenlängen in km an. Versuchen Sie es

gmt grdfft Bathymetry.nc Gravity.grd -Ewk -fg -V > cross_spectra.txt

Untersuchung des Pre-FFT-Gitters nach Trendentfernung, Punktsymmetriereflexion und Verjüngung
angewendet wurden, sowie die Speicherung der realen und imaginären Komponenten des Rohspektrums von
Versuchen Sie es mit den Daten in topo.nc

gmt grdfft topo.nc -N+w+z -fg -V

Sie können jetzt Diagramme der Daten in topo_taper.nc, topo_real.nc und topo_imag.nc erstellen.

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