img2grdgmt - Online în cloud

PROGRAM:

NUME

img2grd - Extrage subset de fișier img în format Mercator sau Geographic

REZUMAT

img2grd imgfile grdfile regiune tip [ ] [ [minlat/maxlat] ] [ ] [ minute ] [ ] [ navg
] [ [scară] ] [ [nivel] ] [ maxlon ] [ -n]

Notă: Nu este permis niciun spațiu între indicatorul opțiunii și argumentele asociate.

DESCRIERE

img2grd citește un fișier în format img, extrage un subset și îl scrie într-un fișier grilă. The -M
opțiunea dictează dacă proiecția Spherical Mercator a fișierului img este sau nu
păstrat sau dacă o grilă geografică ar trebui scrisă prin anularea proiecției Mercator. Dacă
grila geografică este selectată, puteți solicita și o reeșantionare pe exact -R dat.

NECESARE ARGUMENTE

imgfile
Un fișier în format Mercator img, cum ar fi câmpurile gravitaționale maritime sau topografiile fundului mării
estimat din datele altimetrului satelit de Sandwell și Smith. Dacă utilizatorul a setat
o variabilă de mediu $GMT_DATADIR, Apoi img2grd va încerca să găsească imgfile in
$GMT_DATADIR; altfel va încerca să se deschidă imgfile direct.

-Ggrdfile
grdfile este numele fișierului grilă de ieșire.

-R[unitate]vest/est/sud/nord[/zmin/zmax][r]
vest, est, sud, și nord specificați regiunea de interes și puteți specifica
acestea în grade zecimale sau în format [+-]dd:mm[:ss.xxx][W|E|S|N]. Adăuga r dacă mai jos
Coordonatele hărții din stânga și din dreapta sus sunt date în loc de w/e/s/n. Cei doi
stenografii -Rg și -Rd reprezintă domeniul global (0/360 și -180/+180 în longitudine
respectiv, cu -90/+90 în latitudine). Alternativ, pentru crearea grilei, dați
Rcodelon/lat/nx/ny, În cazul în care cod este o combinație de 2 caractere de L, C, R (pentru stânga,
centru sau dreapta) și T, M, B pentru sus, mijloc sau jos. de exemplu, BL pentru stânga jos.
Aceasta indică ce punct dintr-o regiune dreptunghiulară este lon/lat coordonate se referă
la și dimensiunile grilei nx și ny cu spații ale grilei prin -I este folosit pentru a crea
regiunea corespunzătoare. Ca alternativă, specificați numele unui fișier grilă existent
si -R setările (și spația dintre grilă, dacă este cazul) sunt copiate din grilă.
Utilizarea -Runitate asteapta coordonate proiectate (carteziane) compatibile cu alese -J
și proiectăm invers pentru a determina regiunea geografică dreptunghiulară reală. Pentru
vedere în perspectivă (-p), opțional adăugați /zmin/zmax. În cazul vederii în perspectivă
(-p), un interval z (zmin, zmax) poate fi atașat pentru a indica a treia dimensiune. Acest
trebuie făcut numai atunci când utilizați -Jz opțiunea, nu atunci când utilizați numai -p opțiune.
În acest din urmă caz, este trasată o vedere în perspectivă a planului, fără a treia
dimensiune.

OPTIONAL ARGUMENTE

-C Setați coordonatele Mercator x și y relativ la centrul de proiecție [Implicit este
relativ la colțul din stânga jos al grilei]. Necesită -M.

-D[minlat/maxlat]
Utilizați intervalul extins de latitudine -80.738/+80.738. Alternativ, anexați
minlat/maxlat ca întinderea latitudinii fișierului img de intrare. [Implicit este
-72.006/72.006]. De obicei, nu este necesar, deoarece putem determina măsura de la
verificarea dimensiunii fișierului.

-E Poate fi folosit când -M nu este setat să forțeze grila finală să aibă exact același lucru
regiune conform cererii cu -R. În mod implicit, regiunea finală este o proiecție directă a
regiunea originală Mercator și se va extinde de obicei puțin dincolo de
intervalul de latitudine solicitat și, în plus, creșterea grilei în latitudine nu
potriviți cu incrementul de longitudine. Cu toate acestea, reeșantionarea suplimentară introduce mici
erori de interpolare și ar trebui utilizat numai dacă grila de ieșire trebuie să se potrivească cu
regiunea cerută și au x_inc = y_inc. În acest caz regiunea stabilită de -R trebuie să fie
dat în multipli ai incrementului (de ex., -R0/45/45/72).

-I Indica minute ca lățimea unui pixel img de intrare în minute de longitudine.
[Implicit este 2.0]. De obicei, nu este necesar, deoarece putem determina dimensiunea pixelilor din
verificarea dimensiunii.

-M Ieșiți o grilă sferică Mercator [Implicit este o grilă geografică lon/lat]. The
Se păstrează proiecția sferică Mercator a fișierului img, astfel încât regiunea -R
setat de utilizator este ușor modificat; regiunea modificată corespunde marginilor
de pixeli [sau grupuri de navg pixeli]. Antetul fișierului grilă este setat astfel încât x și
lungimile axei y reprezintă distanța de la marginile de vest și de sud ale imaginii,
măsurată în unități implicite de utilizator, cu -Jm1 și ajustat -R. Prin setarea
lipsă PROJ_ ELIPSOID = Sferă, utilizatorul poate face suprapuneri cu reglajul -R
astfel încât să se potrivească. Vedea EXEMPLE de mai jos. Ajustat -R este scris și în grilă
remarca antet, deci poate fi găsită mai târziu. Vedea -C pentru a seta coordonatele relativ la
centru de proiecție.

-Nnavg Media valorilor din pixelii img de intrare în navg by navg pătrate și creați
un pixel de ieșire pentru fiecare astfel de pătrat. Dacă este folosit cu -T3 va raporta o medie
constrângere între 0 și 1. Dacă este folosit cu -T2 rezultatul va fi valoarea medie a datelor
sau NaN în funcție de faptul dacă constrângerea medie este > 0.5. navg trebuie să se împartă uniform
în dimensiunile fișierului img în pixeli. [Mod implicit 1 nu face o medie].

-S[scară]
Înmulțiți valorile fișierului img cu scară înainte de stocarea în fișierul grid. [Implicit este
1.0]. Pentru fișierele img recente: fișierele img topo sunt stocate în contoare (corectate) [-S1];
fișiere gravitaționale în aer liber în mGal*10 [-S0.1 pentru a obține mGal]; fișiere de deviere verticală în
micro-radiani*10 [-S0.1 pentru a obține micro-radiani], fișiere cu gradient gravitațional vertical
Eotvos*50 [-S0.02 pentru a obține Eotvos, sau -S0.002 pentru a obține mGal/km]). Daca nu scară este dată
încercăm să determinăm scara examinând numele fișierului pentru indicii.

-Ttip tip se ocupă de codificarea informațiilor de constrângere. tip = 0 indică faptul că nu
astfel de informații sunt codificate în fișierul img (utilizat pentru versiunile anterioare 1995 ale fișierului
date gravitaționale) și obține toate datele. tip > 0 indică faptul că informațiile de constrângere sunt
codificat (versiunile 1995 și ulterioare (actuale) ale fișierelor img) astfel încât să se poată
produceți un fișier grilă după cum urmează: -T1 primește valorile datelor în toate punctele, -T2 primește date
valori la punctele constrânse și NaN la punctele interpolate; -T3 primește 1 la
puncte constrânse și 0 la punctele interpolate [Implicit este 1].

-V[nivel] (Mai mult ...)
Selectați nivelul de verbozitate [c]. Recomandat în mod special aici, deoarece este util de văzut
cum sunt reglate coordonatele.

-Wmaxlon
Indica maxlon ca lungime maximă de longitudine a fișierului img de intrare. Versiuni
din 1995 au avut maxlon = 360.0, în timp ce unele fișiere anterioare aveau maxlon = 390.0.
[Implicit este 360.0].

-^ or doar -
Imprimă un mesaj scurt despre sintaxa comenzii, apoi iese (NOTĂ: pe Windows
foloseste doar -).

-+ or doar +
Imprimați un mesaj extins de utilizare (ajutor), inclusiv explicația oricăruia
opțiunea specifică modulului (dar nu opțiunile comune GMT), apoi iese.

-? or Nu. argumente
Apoi imprimați un mesaj complet de utilizare (ajutor), inclusiv explicația opțiunilor
iesirile.

--versiune
Tipăriți versiunea GMT și ieșiți.

--show-datadir
Imprimați calea completă către directorul de partajare GMT și ieșiți.

GEOGRAFIC EXEMPLE

-M opțiunea ar trebui să fie exclusă dacă aveți nevoie ca grila de ieșire să fie geografică
coordonate. Pentru a extrage date din regiune -R-40/40/-70/-30 de la world_grav.img.7.2 și
reproiectați pentru a obține coordonatele geografice, puteți încerca

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

Deoarece spațierea latitudinii din fișierul img este echidistantă în unități Mercator, the
grila rezultată nu se va potrivi cu cea specificată -R exact, iar distanța dintre latitudini nu va fi
egal cu distanța de longitudine. Dacă aveți nevoie de o potrivire exactă cu dvs -R și aceeași distanță
în longitudine și latitudine, utilizați -E opţiune:

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

MERCATOR EXEMPLE

Deoarece fișierele img sunt într-o proiecție Mercator, NU ar trebui să extrageți o grilă geografică
dacă planul tău este să faci o hartă Mercator. Dacă ai făcut asta ai terminat de proiectat și
reproiecția grilei, pierzând detaliile cu lungimea de undă scurtă. Mai bine de folosit -M și trasează grila
folosind o proiecție liniară cu aceeași scară ca și proiecția Mercator dorită (vezi GMT
Exemplul 29). Pentru a extrage date din regiune -R-40/40/-70/-30 de la world_grav.img.7.2, alerga

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

Rețineți că -V opțiunea ne spune că intervalul a fost ajustat la
-R-40/40/-70.0004681551/-29.9945810754. De asemenea, putem folosi grdinfo pentru a găsi că fișierul grid
antetul arată că regiunea sa este -R0/80/0/67.9666667 Acesta este intervalul de x,y pe care îl vom obține
dintr-o proiecție sferică Mercator folosind -R-40/40/-70.0004681551/-29.9945810754 și
-Jm1. Astfel, pentru a lua ship.lonlatgrav și a-l folosi pentru a proba merc_grav.nc, putem face acest lucru:

gmt set PROJ_ELLIPSOID Sphere

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

Se recomandă utilizarea metodei de mai sus de proiectare și dezaproiectare a datelor în astfel de cazuri
o aplicație, deoarece atunci există o singură etapă de interpolare (în grdtrack). Daca unul
mai întâi încearcă să convertească fișierul grid în lon,lat și apoi eșantionează, sunt două
etapele de interpolare (în conversie și în eșantionare).

Pentru a face o grilă lon,lat din grila de mai sus putem folosi

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

În unele cazuri, acest lucru nu va fi ușor ca -R în cele două sisteme de coordonate poate să nu se alinieze
bine. Când se întâmplă acest lucru, putem folosi și (de fapt, poate fi întotdeauna mai bine să folosim)

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

Pentru a face o hartă Mercator a regiunii de mai sus, să presupunem valoarea noastră gmt.conf pentru
PROJ_LENGTH_UNIT este inch. Apoi, deoarece fișierul merc_grav.nc de mai sus este proiectat cu -Jm1i
are o lățime de 80 de inci. Putem realiza o hartă de 8 inci lățime folosind -Jx0.1i pe orice program de hărți
aplicat acestei grile (de ex. grdcontour, grdimage, grdview), iar apoi pentru suprapuneri care
lucrați în lon,lat (de ex., psxy, pscoast) putem folosi ajustarea de mai sus -R și -Jm0.1 pentru a obține
cele două sisteme să se potrivească.

Totuși, putem fi mai deștepți decât asta. Realizând că fișierul img de intrare avea pixeli 2.0
minute lat (sau verificând nx și ny cu grdinfo merc_grav.nc) ne dăm seama că
merc_grav.nc a folosit rezoluția completă a fișierului img și are 2400 x 2039 pixeli și
la 8 inci lățimea este de 300 de pixeli pe inch. Noi decidem că nu avem nevoie de atât de mulți și noi
va fi mulțumit cu 100 de pixeli pe inch, așa că dorim să facem o medie a datelor în 3 pe 3
pătrate. (Dacă dorim o diagramă de contur, probabil că vom alege să facem o medie a datelor mult mai mult
(de exemplu, 6 cu 6) pentru a obține contururi netede.) Deoarece 2039 nu este divizibil cu 3, vom obține un
diferite ajustate -R de data asta:

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

De data aceasta aflăm că regiunea ajustată este -R-40/40/-70.023256525/-29.9368261101 și
ieșirea este de 800 x 601 pixeli, o dimensiune mai bună pentru noi. Acum putem crea un artificial
fișier de iluminare pentru această utilizare grdgradient:

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

iar dacă avem și un fișier CPT numit „grav.cpt” putem crea o hartă în relief umbrită în culori
asa:

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

Să presupunem că doriți să obțineți numai valorile de date constrânse dintr-un fișier img, în lat/lon
coordonate. Apoi fugi img2grd cu -T2 opțiune, utilizare grd2xyz a arunca valorile, pipe
prin grep -v NaN pentru a elimina NaN-urile și trece prin neaproape cu inversul
proiecție ca mai sus.

Utilizați img2grdgmt online folosind serviciile onworks.net