img2grdgmt – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

img2grd – Teilmenge der IMG-Datei im Mercator- oder Geographic-Format extrahieren

ZUSAMMENFASSUNG

img2grd Bilddatei grddatei Region tippe [ ] [ [Minlat/Maxlat] ] [ ] [ Minuten ] [ ] [ Navigation
] [ [Treppe] ] [ [Grad des ] ] [ maxlon ] [ -n]

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

BESCHREIBUNG

img2grd liest eine Datei im IMG-Format, extrahiert eine Teilmenge und schreibt sie in eine Grid-Datei. Der -M
Die Option bestimmt, ob die sphärische Mercator-Projektion der IMG-Datei vorhanden ist oder nicht
erhalten bleibt oder ob ein geografisches Gitter durch Rückgängigmachen der Mercator-Projektion geschrieben werden soll. Wenn
Wenn das geografische Raster ausgewählt ist, können Sie auch eine Neuberechnung mit dem genauen Raster anfordern -R gegeben.

ERFORDERLICH ARGUMENTE

Bilddatei
Eine Datei im Mercator-IMG-Format, z. B. die Felder Meeresschwere oder Meeresbodentopographie
geschätzt aus Satellitenhöhenmesserdaten von Sandwell und Smith. Wenn der Benutzer dies festgelegt hat
eine Umgebungsvariable $GMT_DATADIR und dann img2grd werde versuchen zu finden Bilddatei in
$GMT_DATADIR; Andernfalls wird versucht, sich zu öffnen Bilddatei direkt.

-Ggrddatei
grddatei ist der Name der Ausgabe-Grid-Datei.

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

OPTIONAL ARGUMENTE

-C Legen Sie die x- und y-Mercator-Koordinaten relativ zum Projektionszentrum fest [Standard ist
relativ zur unteren linken Ecke des Rasters]. Erfordert -M.

-D[Minlat/Maxlat]
Nutzen Sie den erweiterten Breitengradbereich -80.738/+80.738. Alternativ anhängen
Minlat/Maxlat als Breitengrad der Eingabe-IMG-Datei. [Standard ist
-72.006/72.006]. Normalerweise nicht erforderlich, da wir das Ausmaß anhand dessen ermitteln können
Überprüfung der Dateigröße.

-E Kann verwendet werden, wenn -M ist nicht so eingestellt, dass das endgültige Raster genau dasselbe hat
Region wie gewünscht mit -R. Standardmäßig ist die endgültige Region eine direkte Projektion von
die ursprüngliche Mercator-Region und erstreckt sich typischerweise leicht darüber hinaus
angeforderten Breitengradbereich, und darüber hinaus ist die Rasterinkrementierung im Breitengrad nicht der Fall
mit dem Längengradinkrement übereinstimmen. Das zusätzliche Resampling führt jedoch zu geringen Auswirkungen
Interpolationsfehler und sollten nur verwendet werden, wenn das Ausgaberaster mit dem übereinstimmen muss
angeforderte Region und haben x_inc = y_inc. In diesem Fall die von festgelegte Region -R muß
angegeben in Vielfachen des Inkrements (z. B. -R0/45/45/72).

-I Zeigen Minuten als Breite eines Eingabebildpixels in Längengradminuten.
[Standard ist 2.0]. Normalerweise nicht erforderlich, da wir die Pixelgröße daraus bestimmen können
Kontrolle der Größe.

-M Gibt ein sphärisches Mercator-Gitter aus [Standard ist ein geografisches Längen-/Breitengradgitter]. Der
Die sphärische Mercator-Projektion der IMG-Datei bleibt erhalten, sodass die Region -R
Der vom Benutzer festgelegte Wert wird geringfügig geändert. Der geänderte Bereich entspricht den Kanten
von Pixeln [oder Gruppen von Navigation Pixel]. Der Header der Rasterdatei ist so eingestellt, dass x und
Die Längen der y-Achse stellen den Abstand vom West- und Südrand des Bildes dar.
gemessen in benutzerdefinierten Standardeinheiten, mit -Jm1 und die angepasst -R. Durch Einstellen der
Standard PROJ_ ELLIPSOID = Kugel, der Benutzer kann Überlagerungen mit den angepassten erstellen -R
damit sie zusammenpassen. Sehen Beispiele: unter. Die angepasst -R steht auch im Raster
Kopfzeilenbemerkung, damit sie später gefunden werden kann. Sehen -C Koordinaten relativ zu festlegen
Projektionszentrum.

-NNavigation Mitteln Sie die Werte in den eingegebenen Bildpixeln Navigation by Navigation Quadrate und erstellen
ein Ausgabepixel für jedes dieser Quadrate. Bei Verwendung mit -T3 Es wird ein Durchschnitt gemeldet
Einschränkung zwischen 0 und 1. Bei Verwendung mit -T2 Die Ausgabe erfolgt als durchschnittlicher Datenwert
oder NaN, je nachdem, ob die durchschnittliche Einschränkung > 0.5 ist. Navigation muss sich gleichmäßig aufteilen
in die Abmessungen der Bilddatei in Pixel umwandeln. [Standard 1 führt keine Mittelung durch].

-S[Treppe]
Multiplizieren Sie die Werte der IMG-Datei mit Treppe vor dem Speichern in der Rasterdatei. [Standard ist
1.0]. Für aktuelle IMG-Dateien: IMG-Topodateien werden in (korrigierten) Metern gespeichert [-S1];
Freiluft-Schwerkraftdateien in mGal*10 [-S0.1, um mGal zu erhalten]; vertikale Ablenkungsdateien in
Mikroradiant*10 [-S0.1, um Mikroradianten zu erhalten], vertikale Schwerkraftgradientendateien in
Eotvos*50 [-S0.02, um Eotvos zu erhalten, oder -S0.002, um mGal/km zu erhalten]). Wenn nein Treppe ist gegeben
Wir versuchen den Maßstab zu bestimmen, indem wir den Dateinamen auf Hinweise untersuchen.

-Ttippe tippe kümmert sich um die Kodierung von Einschränkungsinformationen. tippe = 0 zeigt an, dass nein
Solche Informationen sind in der IMG-Datei kodiert (wird für Versionen vor 1995 verwendet).
Schwerkraftdaten) und ruft alle Daten ab. tippe > 0 gibt an, dass die Einschränkungsinformationen vorhanden sind
verschlüsselt (1995 und spätere (aktuelle) Versionen der IMG-Dateien), so dass man
Erstellen Sie eine Rasterdatei wie folgt: -T1 erhält Datenwerte an allen Punkten, -T2 erhält Daten
Werte an eingeschränkten Punkten und NaN an interpolierten Punkten; -T3 bekommt 1 bei
eingeschränkte Punkte und 0 an interpolierten Punkten [Standard ist 1].

-V[Grad des ] (Mehr ...)
Wählen Sie die Ausführlichkeitsstufe [c] aus. Hier besonders empfehlenswert, da es hilfreich ist, es zu sehen
wie die Koordinaten angepasst werden.

-Wmaxlon
Zeigen maxlon als maximale Längenausdehnung der Eingabe-IMG-Datei. Versionen
seit 1995 maxlon = 360.0, während einige frühere Dateien dies hatten maxlon = 390.0.
[Standard ist 360.0].

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

GEOGRAFIE Beispiele:

Das -M Die Option sollte ausgeschlossen werden, wenn das Ausgaberaster geografisch sein soll
Koordinaten. Um Daten in der Region zu extrahieren -R-40/40/-70/-30 ab world_grav.img.7.2 und
Neuprojektieren, um geografische Koordinaten zu erhalten, können Sie versuchen

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V

Da der Breitengradabstand in der IMG-Datei in Mercator-Einheiten äquidistant ist, ist der
Das resultierende Raster stimmt nicht mit dem angegebenen überein -R genau und der Breitengradabstand nicht
gleich dem Längengradabstand. Wenn Sie eine genaue Übereinstimmung mit Ihrem benötigen -R und der gleiche Abstand
Verwenden Sie für Längen- und Breitengrad die -E Option:

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V

Merkator Beispiele:

Da es sich bei den IMG-Dateien um eine Mercator-Projektion handelt, sollten Sie KEIN geografisches Gitter extrahieren
wenn Sie vorhaben, eine Mercator-Karte zu erstellen. Wenn Sie das getan haben, ist die Projektion beendet und
Das Raster wird neu projiziert, wodurch kurzwellige Details verloren gehen. Besser zu verwenden -M und zeichnen Sie das Raster
Verwenden einer linearen Projektion mit demselben Maßstab wie die gewünschte Mercator-Projektion (siehe GMT
Beispiel 29). Um Daten in der Region zu extrahieren -R-40/40/-70/-30 ab world_grav.img.7.2, Lauf

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V

Beachten Sie, dass die -V Option sagt uns, dass der Bereich angepasst wurde
-R-40/40/-70.0004681551/-29.9945810754. Wir können auch verwenden grdinfo um herauszufinden, dass es sich um die Rasterdatei handelt
Der Header zeigt die Region an -R0/80/0/67.9666667 Dies ist der Bereich von x,y, den wir erhalten werden
aus einer sphärischen Mercator-Projektion mit -R-40/40/-70.0004681551/-29.9945810754 und
-Jm1. Um also „ship.lonlatgrav“ zu nehmen und damit merc_grav.nc abzutasten, können wir Folgendes tun:

gmt set PROJ_ELLIPSOID Sphere

gmt mapproject -R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | GMT-Kartenprojekt \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > Ship.lonlatgravsat

Es wird empfohlen, die oben beschriebene Methode zum Projizieren und Deprojizieren der Daten zu verwenden
eine Anwendung, denn dann gibt es nur einen Interpolationsschritt (in grdspur). Wenn man
Versuchen Sie zunächst, die Rasterdatei in lon,lat zu konvertieren und sie dann abzutasten. Es gibt zwei
Interpolationsschritte (bei der Konvertierung und beim Sampling).

Um ein Längen- und Breitengradgitter aus dem obigen Gitter zu erstellen, können wir es verwenden

gmt grdproject merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc

In einigen Fällen wird dies nicht einfach sein -R in den beiden Koordinatensystemen möglicherweise nicht übereinstimmen
Also. In diesem Fall können wir auch Folgendes verwenden (tatsächlich ist es möglicherweise immer besser):

gmt grd2xyz merc_grav.nc | GMT-Kartenprojekt \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
gmt-Oberfläche -R-40/40/-70/70 -I2m -Ggrav.nc

Um eine Mercator-Karte der oben genannten Region zu erstellen, nehmen wir unseren gmt.conf-Wert für an
PROJ_LENGTH_UNIT ist Zoll. Dann wird da oben die Datei merc_grav.nc mit projiziert -Jm1i
es ist 80 Zoll breit. Wir können eine Karte mit einer Breite von 8 Zoll erstellen, indem wir verwenden -Jx0.1i auf allen Kartenprogrammen
auf dieses Raster angewendet (z. B. Grundkontur, grdimage, grdview), und dann für Overlays which
Arbeit in Lon, Lat (z. B. Pseudo, pscoast) können wir die oben angepassten verwenden -R und -Jm0.1 zu bekommen
dass die beiden Systeme zusammenpassen.

Wir können jedoch schlauer sein. Erkennen, dass die Eingabe-IMG-Datei Pixel 2.0 hatte
Minuten breit (oder Überprüfung von nx und ny mit grdinfo merc_grav.nc), das ist uns klar
merc_grav.nc hat die volle Auflösung der IMG-Datei verwendet und sie hat 2400 x 2039 Pixel
Bei einer Breite von 8 Zoll sind das 300 Pixel pro Zoll. Wir entscheiden, dass wir nicht so viele brauchen und wir
wird mit 100 Pixeln pro Zoll zufrieden sein, also wollen wir die Daten in 3 mal 3 mitteln
Quadrate. (Wenn wir ein Konturdiagramm wünschen, werden wir die Daten wahrscheinlich stärker mitteln
(z. B. 6 mal 6, um glatte Konturen zu erhalten.) Da 2039 nicht durch 3 teilbar ist, erhalten wir a
unterschiedlich angepasst -R diesmal:

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V

Dieses Mal finden wir die angepasste Region -R-40/40/-70.023256525/-29.9368261101 und die
Die Ausgabe beträgt 800 x 601 Pixel, eine bessere Größe für uns. Jetzt können wir ein Künstliches erstellen
Beleuchtungsdatei für diese Verwendung gradient:

gmt grdgradient merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

Und wenn wir auch eine CPT-Datei namens „grav.cpt“ haben, können wir eine farbschattierte Reliefkarte erstellen
so was:

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K > map.ps
gmt psbasemap -R-40/40/-70.023256525/-29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

Angenommen, Sie möchten nur die eingeschränkten Datenwerte aus einer IMG-Datei in Breiten- und Längengraden abrufen
Koordinaten. Dann renne img2grd an. Nach der Installation können Sie HEIC-Dateien mit der -T2 Option, verwenden grd2xyz Um die Werte zu sichern, Pipe
durch grep -v NaN, um NaNs zu eliminieren, und weiterleiten Kartenprojekt mit der Umkehrung
Projektion wie oben.

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