img2grdgmt - Online in de cloud

PROGRAMMA:

NAAM

img2grd - Extraheer een subset van het img-bestand in Mercator- of geografisch formaat

KORTE INHOUD

img2grd imgbestand grdbestand regio type dan: [ ] [ [minlat/maxlat] ] [ ] [ minuten ] [ ] [ navigeren
] [ [schaal] ] [ [niveau] ] [ maxlon ] [ -n]

Opmerking: Er is geen spatie toegestaan ​​tussen de optievlag en de bijbehorende argumenten.

PRODUCTBESCHRIJVING

img2grd leest een bestand in img-formaat, extraheert een subset en schrijft deze naar een rasterbestand. De -M
optie bepaalt of de sferische Mercator-projectie van het img-bestand al dan niet is
behouden blijft of dat er een geografisch raster moet worden geschreven door de Mercator-projectie ongedaan te maken. Als
geografisch raster is geselecteerd, kunt u ook een herbemonstering naar het exacte raster aanvragen -R gegeven.

VERPLICHT ARGUMENTEN

imgbestand
Een Mercator img-bestand, zoals de mariene zwaartekracht of de topografievelden van de zeebodem
geschat op basis van satelliethoogtemetergegevens door Sandwell en Smith. Als de gebruiker dit heeft ingesteld
een omgevingsvariabele $GMT_DATADIRdan img2grd zal proberen te vinden imgbestand in
$GMT_DATADIR; anders zal het proberen te openen imgbestand direct.

-Ggrdbestand
grdbestand is de naam van het uitvoerrasterbestand.

-R[eenheid]west/oosten/zuiden/noorden[/zmin/zmax][R]
west, oosten, zuiden en noorden specificeer de regio van belang, en u kunt specificeren
ze in decimale graden of in [+-]dd:mm[:ss.xxx][W|E|S|N] formaat. toevoegen r indien lager
links en rechtsboven kaartcoördinaten worden gegeven in plaats van w/e/s/n. De twee
steno -Rg en -Rd staat voor globaal domein (0/360 en -180/+180 in lengtegraad
respectievelijk, met -90/+90 in breedte). Als alternatief voor het maken van een raster, geef
Rcodelon/lat/nx/ny, Waar code is een combinatie van 2 tekens van L, C, R (voor links,
midden of rechts) en T, M, B voor boven, midden of onder. bijv. BL voor linksonder.
Dit geeft aan welk punt op een rechthoekig gebied de lon/lat coördinaat verwijst
to, en de afmetingen van het raster nx en ny met rasterafstanden via -I wordt gebruikt om te creëren
de bijbehorende regio. U kunt ook de naam van een bestaand rasterbestand specificeren
en -R instellingen (en rasterafstand, indien van toepassing) worden gekopieerd uit het raster.
gebruik -Reenheid verwacht geprojecteerde (Cartesiaanse) coördinaten die compatibel zijn met gekozen -J
en we projecteren omgekeerd om het werkelijke rechthoekige geografische gebied te bepalen. Voor
perspectief (-p), eventueel toevoegen /zmin/zmax. In het geval van perspectiefweergave
(-p), een z-bereik (zmin, zmax) kan worden toegevoegd om de derde dimensie aan te geven. Dit
hoeft alleen te worden gedaan bij gebruik van de -Jz optie, niet bij gebruik van alleen de -p optie.
In het laatste geval wordt een perspectivisch aanzicht van het vlak uitgezet, zonder derde
dimensie.

OPTIONELE ARGUMENTEN

-C Stel de x- en y Mercator-coördinaten in ten opzichte van het projectiecentrum [standaard is
ten opzichte van de linkerbenedenhoek van het raster]. Vereist -M.

-D[minlat/maxlat]
Gebruik het uitgebreide breedtegraadbereik -80.738/+80.738. U kunt ook toevoegen
minlat/maxlat als de breedtegraad van het invoer-img-bestand. [Standaard is
-72.006/72.006]. Meestal niet nodig, omdat we de omvang ervan kunnen bepalen
inspectie van de bestandsgrootte.

-E Kan gebruikt worden wanneer -M is niet ingesteld om het uiteindelijke raster te dwingen exact hetzelfde te hebben
regio zoals gevraagd met -R. Standaard is het uiteindelijke gebied een directe projectie van
de oorspronkelijke Mercatorregio en strekt zich doorgaans iets verder uit dan de
het gewenste breedtegraadbereik, en bovendien geldt dat niet voor de verhoging van het raster in de breedtegraad
overeenkomen met de toename van de lengtegraad. De extra herbemonstering introduceert echter kleine
interpolatiefouten en mag alleen worden gebruikt als het uitgangsraster moet overeenkomen met de
gevraagde regio en hebben x_inc = y_inc. In dit geval is de regio ingesteld door -R moet
gegeven in veelvouden van de stapgrootte (bijv. -R0/45/45/72).

-I Aangeven minuten als de breedte van een ingevoerde img-pixel in minuten lengtegraad.
[Standaard is 2.0]. Meestal niet nodig, omdat we de pixelgrootte kunnen bepalen
inspectie van de maat.

-M Voer een bolvormig Mercator-raster uit [standaard is een geografisch lon/lat-raster]. De
De sferische Mercator-projectie van het img-bestand blijft behouden, zodat de regio -R
ingesteld door de gebruiker is enigszins aangepast; het gewijzigde gebied komt overeen met de randen
van pixels [of groepen van navigeren pixels]. De header van het rasterbestand is zo ingesteld dat de x en
De lengtes van de y-as vertegenwoordigen de afstand vanaf de west- en zuidrand van de afbeelding,
gemeten in standaardeenheden van de gebruiker, met -Jm1 en de aangepaste -R. Door het instellen van de
verzuim PROJ_ ELLIPSOÏDE = Bol, de gebruiker kan overlays maken met de aangepaste -R
zodat ze overeenkomen. Zien Voorbeelden onderstaand. De aangepaste -R staat ook in het raster geschreven
headeropmerking, zodat deze later kan worden gevonden. Zien -C coördinaten ten opzichte van instellen
projectie centrum.

-Nnavigeren Gemiddelde van de waarden in de ingevoerde img-pixels navigeren by navigeren vierkanten en creëren
één uitvoerpixel voor elk dergelijk vierkant. Indien gebruikt met -T3 het rapporteert een gemiddelde
beperking tussen 0 en 1. Indien gebruikt met -T2 de uitvoer zal de gemiddelde gegevenswaarde zijn
of NaN, afhankelijk van de vraag of de gemiddelde beperking > 0.5 is. navigeren moet gelijkmatig verdeeld worden
in de afmetingen van het img-bestand in pixels. [Standaard 1 doet geen middeling].

-S[schaal]
Vermenigvuldig de waarden van het img-bestand met schaal voordat u het in een rasterbestand opslaat. [Standaard is
1.0]. Voor recente img-bestanden: img topo-bestanden worden opgeslagen in (gecorrigeerde) meters [-S1];
zwaartekrachtbestanden in de vrije lucht in mGal*10 [-S0.1 om mGal te krijgen]; verticale afbuigingsbestanden in
microradialen*10 [-S0.1 om microradialen te krijgen], verticale zwaartekrachtgradiëntbestanden in
Eotvos*50 [-S0.02 om Eotvos te krijgen, of -S0.002 om mGal/km te krijgen]). Als Nee schaal is gegeven
we proberen de schaal te bepalen door de bestandsnaam te onderzoeken op aanwijzingen.

-T type dan: type dan: verzorgt de codering van beperkingsinformatie. type dan: = 0 geeft aan dat nee
Dergelijke informatie wordt gecodeerd in het img-bestand (gebruikt voor versies van vóór 1995 van de
zwaartekrachtgegevens) en haalt alle gegevens op. type dan: > 0 geeft aan dat er beperkingsinformatie is
gecodeerd (1995 en latere (huidige) versies van de img-bestanden) zodat men dat kan
maak als volgt een rasterbestand: -T1 krijgt gegevenswaarden op alle punten, -T2 krijgt gegevens
waarden op beperkte punten en NaN op geïnterpoleerde punten; -T3 krijgt 1 op
beperkte punten en 0 bij geïnterpoleerde punten [standaard is 1].

-V[niveau] (meer ...)
Selecteer breedsprakigheidsniveau [c]. Hier vooral aanbevolen, omdat het nuttig is om te zien
hoe de coördinaten worden aangepast.

-Wmaxlon
Aangeven maxlon als de maximale lengtegraad van het ingevoerde img-bestand. Versies
sinds 1995 gehad maxlon = 360.0, terwijl sommige eerdere bestanden dat wel hadden maxlon = 390.0.
[Standaard is 360.0].

-^ or voor slechts -
Druk een kort bericht af over de syntaxis van de opdracht en sluit vervolgens af (OPMERKING: op Windows
gebruik gewoon -).

-+ or voor slechts +
Druk een uitgebreid gebruiks(help)bericht af, inclusief de uitleg van eventuele
module-specifieke optie (maar niet de algemene GMT-opties), en wordt vervolgens afgesloten.

-? or geen argumenten
Druk een volledig gebruiks(help)bericht af, inclusief de uitleg van de opties, dan
uitgangen.

--versie
GMT-versie afdrukken en afsluiten.

--show-datadir
Druk het volledige pad af naar de GMT-sharedirectory en sluit af.

GEOGRAFISCH Voorbeelden

De -M optie moet worden uitgesloten als het uitvoerraster geografisch moet zijn
coördinaten. Om gegevens in de regio te extraheren -R-40/40/-70/-30 vanaf wereld_grav.img.7.2 en
opnieuw projecteren om geografische coördinaten op te leveren, u kunt het proberen

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

Omdat de breedtegraad in het img-bestand op gelijke afstand ligt in Mercator-eenheden, wordt de
het resulterende raster komt niet overeen met het opgegeven raster -R precies, en de breedtegraadafstand niet
gelijk aan de lengtegraadafstand. Als u een exacte match met uw nodig heeft -R en dezelfde afstand
voor lengte- en breedtegraad gebruikt u de -E optie:

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

MERCATOR Voorbeelden

Omdat de img-bestanden zich in een Mercator-projectie bevinden, mag u GEEN geografisch raster extraheren
als het uw plan is om een ​​Mercator-kaart te maken. Als je dat hebt gedaan, ben je klaar met projecteren en
herprojectie van het raster, waarbij details op de korte golflengte verloren gaan. Beter te gebruiken -M en teken het raster
met behulp van een lineaire projectie met dezelfde schaal als de gewenste Mercator-projectie (zie GMT
Voorbeeld 29). Om gegevens in de regio te extraheren -R-40/40/-70/-30 vanaf wereld_grav.img.7.2, rennen

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

Merk op dat de -V optie vertelt ons dat het bereik is aangepast
-R-40/40/-70.0004681551/-29.9945810754. Wij kunnen ook gebruiken grdinfo om te vinden dat het rasterbestand
header toont de regio ervan -R0/80/0/67.9666667 Dit is het bereik van x,y dat we zullen krijgen
van een bolvormige Mercator-projectie met behulp van -R-40/40/-70.0004681551/-29.9945810754 en
-Jm1. Als we dus ship.lonlatgrav willen nemen en gebruiken om merc_grav.nc te samplen, kunnen we dit doen:

gmt heeft PROJ_ELLIPSOID Sphere ingesteld

gmt mapproject -R-40/40/-70.0004681551/-29.9945810754 -Jm1i schip.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | gmt kaartproject \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > schip.lonlatgravsat

Het wordt aanbevolen om de bovenstaande methode te gebruiken voor het projecteren en het niet projecteren van de gegevens in dergelijke bestanden
een applicatie, omdat er dan slechts één interpolatiestap is (in grdtrack). Als een
probeert eerst het rasterbestand naar lon,lat te converteren en het vervolgens te samplen, er zijn er twee
interpolatiestappen (bij conversie en bij bemonstering).

Om een ​​lon,lat-raster te maken van het bovenstaande raster kunnen we dit gebruiken

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

In sommige gevallen zal dit niet eenvoudig zijn -R in de twee coördinatensystemen zijn mogelijk niet uitgelijnd
Goed. Wanneer dit gebeurt, kunnen we ook gebruiken (het kan zelfs altijd beter zijn om te gebruiken)

gmt grd2xyz merc_grav.nc | gmt kaartproject \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
gmt oppervlak -R-40/40/-70/70 -I2m -Ggrav.nc

Om een ​​Mercator-kaart van de bovenstaande regio te maken, veronderstellen we onze gmt.conf-waarde voor
PROJ_LENGTH_UNIT is inch. Omdat het bovenstaande merc_grav.nc-bestand wordt geprojecteerd met -Jm1i
het is 80 inch breed. We kunnen een kaart van 8 inch breed maken met behulp van -Jx0.1i op alle kaartprogramma's
toegepast op dit raster (bijv. grdcontour, grdbeeld, grdweergave), en vervolgens voor overlays die
werk in lon,lat (bijv. psxy, pkust) kunnen we bovenstaande aangepast gebruiken -R en -Jm0.1 om te krijgen
de twee systemen op elkaar afstemmen.

We kunnen echter slimmer zijn dan dit. Realiserend dat het invoer-img-bestand pixels 2.0 had
minuten breed (of het controleren van de nx en ny met grdinfo merc_grav.nc) realiseren we ons dat
merc_grav.nc gebruikte de volledige resolutie van het img-bestand en het heeft 2400 bij 2039 pixels, en
bij 8 inch breed is dit 300 pixels per inch. We besluiten dat we er niet zoveel nodig hebben en wij
zal tevreden zijn met 100 pixels per inch, dus we willen de gegevens in 3 bij 3 middelen
vierkanten. (Als we een contourplot willen, zullen we er waarschijnlijk voor kiezen om de gegevens veel meer te middelen
(bijvoorbeeld 6 bij 6) om vloeiende contouren te krijgen.) Omdat 2039 niet deelbaar is door 3, krijgen we een
verschillend aangepast -R deze keer:

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

Deze keer vinden we de aangepaste regio -R-40/40/-70.023256525/-29.9368261101 en de
de uitvoer is 800 bij 601 pixels, een beter formaat voor ons. Nu kunnen we een kunstmatig maken
verlichtingsbestand hiervoor gebruiken grdgradiënt:

gmt gradiënt merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

en als we ook een CPT-bestand hebben met de naam "grav.cpt", kunnen we een in kleur gearceerde reliëfkaart maken
soortgelijk:

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

Stel dat u alleen de beperkte gegevenswaarden uit een img-bestand wilt verkrijgen, in lat/lon
coördinaten. Ren dan img2grd met de -T2 optie, gebruik grd2xyz om de waarden te dumpen, pipe
via grep -v NaN om NaN's te elimineren, en door te voeren kaartproject met het omgekeerde
projectie zoals hierboven.

Gebruik img2grdgmt online met behulp van onworks.net-services