grdgradientgmt – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

grdgradient – ​​Richtungsableitung oder -gradient aus einem Gitter berechnen

ZUSAMMENFASSUNG

gradient in_grdfile out_grdfile [ azim[/Azim2] ] [ [a][c][o][n] ] [
[s|p]Azim/Elev[/Umgebungs-/diffundieren/spiegelnd/leuchten] ] [ Flagge ] [ [e][t][amp][/Sigma[/Offset]]
] [ Region ] [ Slopefile ] [ [Grad des ] ] [ -fg ] [ -n]

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

BESCHREIBUNG

gradient kann verwendet werden, um die Richtungsableitung in einer bestimmten Richtung zu berechnen (-A),
oder die Richtung (-S) [und die Größe (-D)] des Vektorgradienten der Daten.

Die geschätzten Werte in der ersten/letzten Zeile/Spalte der Ausgabe hängen von den Randbedingungen ab (siehe
-L).

ERFORDERLICH ARGUMENTE

in_grdfile
2D-Gitterdatei, aus der die Richtungsableitung berechnet wird. (Siehe GRID-DATEIFORMATE
unten).

-Gout_grdfile
Name der Ausgabegitterdatei für die Richtungsableitung. (Siehe GRID-DATEIFORMATE
unten).

OPTIONAL ARGUMENTE

-Aazim[/Azim2]
Azimutale Richtung für eine Richtungsableitung; azim ist der Winkel in x,y
Ebene gemessen in Grad positiv im Uhrzeigersinn von Norden (der +y-Richtung) in Richtung
Osten (die +x-Richtung). Das Negativ der Richtungsableitung,
-[dz/dx*sin(azim) + dz/dy*cos(azim)], gefunden; Negation ergibt positive Werte
wenn die Steigung von z(x,y) bergab verläuft azim Richtung, der richtige Sinn für
die Beleuchtung eines Bildes schattieren (siehe grdimage und grdview) durch eine Lichtquelle
über der x,y-Ebene, die vom azim Richtung. Geben Sie optional zwei an
Azimute, -Aazim/Azim2, in diesem Fall sind die Gradienten in jeder dieser Richtungen
berechnet und der größere Betrag bleibt erhalten; das ist nützlich für
Beleuchtung von Daten mit zwei Richtungen linienförmiger Strukturen, z. B. -A0/270
Leuchtet von Norden (oben) und Westen (links).

-D[a][c][o][n]
Ermitteln Sie die Richtung des positiven (steigenden) Gradienten der Daten. Stattdessen
Finden Sie den Aspekt (die Gefällerichtung) und verwenden Sie ihn -da. Standardmäßig sind Wegbeschreibungen
gemessen im Uhrzeigersinn von Norden, as azim in -A über. Anhängen c konventionell zu verwenden
Kartesische Winkel, gemessen gegen den Uhrzeigersinn aus der positiven x-Richtung (Osten).
Anhängen o um Ausrichtungen (0-180) statt Richtungen (0-360) zu melden. Anhängen n
um allen Winkeln 90 Grad hinzuzufügen (z. B. um lokale Schläge auf die Oberfläche zu erzielen).

-E[s|p]Azim/Elev[/Umgebungs-/diffundieren/spiegelnd/leuchten]
Berechnen Sie die für die Verwendung geeignete Lambertsche Strahldichte grdimage und grdviewdem „Vermischten Geschmack“. Seine
Die Lambertsche Reflexion geht von einer idealen Oberfläche aus, die das gesamte Licht reflektiert
trifft darauf und die Oberfläche erscheint aus allen Blickrichtungen gleich hell. azim
und Schüler sind der Azimut und die Höhe des Lichtvektors. Optional liefern Umgebungs-
diffundieren spiegelnd leuchten Dabei handelt es sich um Parameter, die die Reflexionseigenschaften steuern
der Oberfläche. Standardwerte sind: 0.55/0.6/0.4/10 Um einige der Werte zu belassen
unangetastet, geben Sie = als neuen Wert an. Zum Beispiel -E60/30/=/0.5 setzt die azim Schüler
und diffundieren auf 60, 30 und 0.5 und belässt die anderen Reflexionsparameter
unberührt. Anhängen s einen einfacheren Lambertschen Algorithmus zu verwenden. Beachten Sie das bei diesem Formular
Sie müssen lediglich die Azimut- und Höhenparameter angeben. Anhängen p zu verwenden, die
Stückweise lineare Peucker-Approximation (einfacher, aber schnellerer Algorithmus; in diesem Fall
die azim und Schüler sind fest auf 315 und 45 Grad verdrahtet. Dies bedeutet, dass auch wenn Sie
Geben Sie andere Werte an, diese werden ignoriert.)

-LFlagge Randbedingung Flagge könnte sein x or y or xy Zeigt an, dass die Daten periodisch im Bereich von liegen
x oder y oder beides, oder Flagge könnte sein g Angabe der geografischen Bedingungen (x und y sind).
Lon und Lat). [Standardmäßig werden „natürliche“ Bedingungen verwendet (zweite partielle Ableitungsnormale).
zur Kante ist Null).]

-Netz][amp][/Sigma[/Offset]]
Normalisierung. [Standard: keine Normalisierung.] Die tatsächlichen Farbverläufe g sind versetzt und
skaliert, um normalisierte Farbverläufe zu erzeugen gn mit einer maximalen Ausgangsgröße von amp.
If amp ist nicht angegeben, Standard amp = 1. Wenn Offset nicht gegeben ist, wird es auf die gesetzt
Durchschnitt von g. -N Erträge gn = amp * (g - Offset)/max(abs(g - Offset)). -Ne
normalisiert sich unter Verwendung einer kumulativen Laplace-Verteilungsausbeute gn = amp * (1.0 -
exp(sqrt(2) * (g - Offset)/ Sigma)) wo Sigma wird unter Verwendung der L1-Norm von geschätzt
(g - Offset), wenn es nicht angegeben ist. -Nt normalisiert sich unter Verwendung eines kumulativen Cauchy
Verteilung ertragreich gn = (2 * amp / PI) * atan( (g - Offset)/ Sigma) woher Sigma
wird unter Verwendung der L2-Norm von (g - Offset), wenn es nicht angegeben ist.

-R[Einheit]xMin/xmax/ymin/ymax[r] (Mehr ...)
Geben Sie die Region von Interesse an. Verwendung der -R Mit dieser Option wird ein Unterabschnitt von ausgewählt
in_grdfile Netz. Wenn dieser Unterabschnitt über die Grenzen des Rasters hinausgeht, wird nur der
gemeinsame Region extrahiert werden.

-SSlopefile
Name der Ausgabegitterdatei mit skalaren Größen von Gradientenvektoren. Erfordert -D
aber macht -G optional.

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

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

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.

HINWEISE

Wenn du nicht weißt was -N Optionen zum Erstellen einer Intensitätsdatei für grdimage or
grdview, ein guter erster Versuch ist -Ne0.6

Normalerweise sind 255 Farbtöne für Visualisierungszwecke mehr als ausreichend. Sie können 75 % der Festplatte einsparen
Speicherplatz durch Anhängen von =nb/a an den Namen der Ausgabedatei out_grdfile.

Wenn Sie mehrere beleuchtete Karten von Teilregionen eines großen Datensatzes erstellen möchten, und Sie
Wenn Sie möchten, dass die Beleuchtungseffekte auf allen Karten konsistent sind, verwenden Sie die -N Option und
den gleichen Wert liefern Sigma und Offset zu gradient für jede Karte. Eine gute Vermutung ist
Offset = 0 und Sigma gefunden von grdinfo -L2 or -L1 angewendet auf einen nicht normalisierten Gradienten grd.

Wenn Sie einfach das brauchen x- oder y-Ableitungen des Gitters, Verwendung grdmath.

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.

Beispiele:

So erstellen Sie eine Datei zur Beleuchtung der Daten in geoid.nc mit exp-normalisierten Farbverläufen in
der Bereich [-0.6,0.6] imitiert Lichtquellen in Nord- und Westrichtung:

gmt grdgradient geoid.nc -A0/270 -Ggradients.nc=nb/a -Ne0.6 -V

So finden Sie die Azimutausrichtungen des Meeresbodengewebes in der Datei topo.nc:

gmt grdgradient topo.nc -Dno -Gazimuths.nc -V

REFERENZEN

Horn, BKP, Hill-Shading and the Reflectance Map, Proceedings of the IEEE, Bd. 69, Nr.
1, Januar 1981, S. 14–47. (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)

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