img2grdgmt - Online nel cloud

PROGRAMMA:

NOME

img2grd - Estrai il sottoinsieme del file img in formato Mercator o Geografico

SINOSSI

img2grd imgfile grdfile regione Digitare [ ] [ [minlat/maxlat]] [] [ verbale ] [ ] [ navigazione
] [ [scala]] [ [livello]] [ maxlon ] [ -n]

Nota: Non è consentito alcuno spazio tra il flag dell'opzione e gli argomenti associati.

DESCRIZIONE

img2grd legge un file in formato img, estrae un sottoinsieme e lo scrive in un file griglia. Il -M
l'opzione determina se la proiezione Spherical Mercator del file img è
conservato o se si deve scrivere una griglia geografica annullando la proiezione di Mercatore. Se
la griglia geografica è selezionata puoi anche richiedere un ricampionamento sull'esatto -R dato.

OBBLIGATORIO ARGOMENTI

imgfile
Un file in formato img Mercator come i campi della gravità marina o della topografia del fondale marino
stimato dai dati dell'altimetro satellitare di Sandwell e Smith. Se l'utente ha impostato
una variabile d'ambiente $GMT_DATADIR, poi img2grd proverò a trovare imgfile in
$GMT_DATADIR; altrimenti proverà ad aprire imgfile direttamente.

-Ggrdfile
grdfile è il nome del file della griglia di output.

-R[unità]ovest/est/sud/nord[/zmin/zmax][R]
ovest, est, sude nord specificare la regione di interesse e si può specificare
in gradi decimali o in formato [+-]dd:mm[:ss.xxx][W|E|S|N]. Aggiungere r se inferiore
vengono fornite le coordinate della mappa a sinistra e in alto a destra invece di w/e/s/n. Il due
scorciatoie -Rg ed -Rod sta per dominio globale (0/360 e -180/+180 in longitudine
rispettivamente, con -90/+90 in latitudine). In alternativa per la creazione della griglia, dai
Rcodelon/lat/nx/ny, Dove codice è una combinazione di 2 caratteri di L, C, R (per sinistra,
al centro o a destra) e T, M, B per l'alto, il centro o il basso. ad esempio, BL per in basso a sinistra.
Questo indica quale punto su una regione rettangolare il lon/lat la coordinata si riferisce
a, e le dimensioni della griglia nx ed ny con spaziatura della griglia via -I è usato per creare
la regione corrispondente. In alternativa, specifica il nome di un file griglia esistente
e la -R le impostazioni (e la spaziatura della griglia, se applicabile) vengono copiate dalla griglia.
utilizzando -Runità si aspetta coordinate proiettate (cartesiane) compatibili con quelle scelte -J
e progettiamo inversamente di determinare l'effettiva regione geografica rettangolare. Per
vista in prospettiva (-p), facoltativamente aggiungi /zmin/zmax. In caso di vista prospettica
(-p), un intervallo z (zmin, zmax) può essere aggiunto per indicare la terza dimensione. Questo
deve essere fatto solo quando si utilizza il -Jz opzione, non quando si utilizza solo il -p opzione.
In quest'ultimo caso viene tracciata una vista prospettica del piano, senza terzo
dimensione.

OPTIONAL ARGOMENTI

-C Imposta le coordinate x e y di Mercatore relative al centro di proiezione [Il valore predefinito è
relativo all'angolo inferiore sinistro della griglia]. Richiede -M.

-D[minlat/maxlat]
Utilizzare l'intervallo di latitudine esteso -80.738/+80.738. In alternativa, append
minlat/maxlat come estensione di latitudine del file img di input. [Il valore predefinito è
-72.006/72.006]. Di solito non richiesto poiché possiamo determinare l'estensione da
controllo della dimensione del file.

-E Può essere utilizzato quando -M non è impostato per forzare la griglia finale ad avere esattamente lo stesso
regione come richiesto con -R. Per impostazione predefinita, la regione finale è una proiezione diretta di
l'originale regione di Mercator e si estenderà tipicamente leggermente oltre il
l'intervallo di latitudine richiesto e inoltre l'incremento della griglia in latitudine non lo fa
corrisponde all'incremento di longitudine. Tuttavia, il ricampionamento extra introduce piccoli
errori di interpolazione e deve essere utilizzato solo se la griglia di uscita deve corrispondere al
regione richiesta e avere x_inc = y_inc. In questo caso la regione impostata da -R deve essere
dato in multipli dell'incremento (.eg, -R0/45/45/72).

-I Indicare verbale come la larghezza di un pixel img di input in minuti di longitudine.
[Il valore predefinito è 2.0]. Di solito non richiesto poiché possiamo determinare la dimensione dei pixel da
controllo delle dimensioni.

-M Output di una griglia Spherical Mercator [L'impostazione predefinita è una griglia geografica lon/lat]. Il
La proiezione sferica di Mercatore del file img viene preservata, in modo che la regione -R
impostato dall'utente viene leggermente modificato; la regione modificata corrisponde ai bordi
di pixel [o gruppi di navigazione pixel]. L'intestazione del file della griglia è impostata in modo che x e
le lunghezze dell'asse y rappresentano la distanza dai bordi ovest e sud dell'immagine,
misurato in unità predefinite dall'utente, con -Jm1 e il regolato -R. Impostando il
difetto PROG_ ELLISSOIDE = Sfera, l'utente può creare sovrapposizioni con la regolazione -R
in modo che corrispondano. Vedere ESEMPI sotto. Il aggiustato -R è scritto anche nella griglia
commento di intestazione, in modo che possa essere trovato in seguito. Vedere -C per impostare le coordinate relative a
centro di proiezione

-Nnavigazione Media i valori nei pixel di input img in navigazione by navigazione quadrati e creare
un pixel di output per ciascuno di questi quadrati. Se usato con -T3 riporterà una media
vincolo tra 0 e 1. Se usato con -T2 l'output sarà il valore medio dei dati
o NaN a seconda che il vincolo medio sia > 0.5. navigazione deve dividere equamente
nelle dimensioni del file img in pixel. [Predefinito 1 non fa la media].

-S[scala]
Moltiplica i valori del file img per scala prima di memorizzare nel file di griglia. [Il valore predefinito è
1.0]. Per i file img recenti: i file img topo sono memorizzati in metri (corretti) [-S1];
file di gravità in aria libera in mgal*10 [-S0.1 per ottenere mgal]; file di deflessione verticale in
microradianti*10 [-S0.1 per ottenere micro-radianti], file di gradiente di gravità verticale in
Éotvos*50 [-S0.02 per ottenere Eotvos, o -S0.002 per ottenere mgal/km]). se no scala viene data
proviamo a determinare la scala esaminando il nome del file per gli indizi.

-TDigitare Digitare gestisce la codifica delle informazioni sui vincoli. Digitare = 0 indica che no
tali informazioni sono codificate nel file img (usato per le versioni precedenti al 1995 del
dati di gravità) e ottiene tutti i dati. Digitare > 0 indica che le informazioni sul vincolo sono
codificati (1995 e versioni successive (attuali) dei file img) in modo che si possa
produrre un file griglia come segue: -T1 ottiene i valori dei dati in tutti i punti, -T2 ottiene dati
valori nei punti vincolati e NaN nei punti interpolati; -T3 ottiene 1 a
punti vincolati e 0 nei punti interpolati [Il valore predefinito è 1].

-V[livello] (Di Più ...)
Seleziona il livello di verbosità [c]. Particolarmente consigliato qui, come è utile vedere
come vengono regolate le coordinate.

-Wmaxlon
Indicare maxlon come estensione massima della longitudine del file img di input. Versioni
dal 1995 hanno avuto maxlon = 360.0, mentre alcuni file precedenti avevano maxlon = 390.0.
[Il valore predefinito è 360.0].

-^ or ad appena -
Stampa un breve messaggio sulla sintassi del comando, quindi esce (NOTA: su Windows
usa solo -).

-+ or ad appena +
Stampa un messaggio di utilizzo esteso (aiuto), inclusa la spiegazione di qualsiasi
opzione specifica del modulo (ma non le opzioni comuni GMT), quindi esce.

-? or no argomenti
Stampare un messaggio di utilizzo completo (aiuto), inclusa la spiegazione delle opzioni, quindi
esce.

--versione
Stampa la versione GMT ed esci.

--show-datadir
Stampa il percorso completo della directory di condivisione GMT ed esci.

GEOGRAFICA ESEMPI

I -M l'opzione dovrebbe essere esclusa se hai bisogno che la griglia di output sia in geografica
coordinate. Per estrarre i dati nella regione -R-40/40/-70/-30 da world_grav.img.7.2 ed
riproiettare per ottenere le coordinate geografiche, puoi provare

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

Poiché la spaziatura della latitudine nel file img è equidistante in unità di Mercatore, il
la griglia risultante non corrisponderà a quella specificata -R esattamente, e la spaziatura della latitudine non lo farà
uguale alla distanza di longitudine. Se hai bisogno di una corrispondenza esatta con il tuo -R e la stessa spaziatura
in longitudine e latitudine, usa il -E opzione:

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

MERCATORE ESEMPI

Poiché i file img sono in una proiezione di Mercatore, NON dovresti estrarre una griglia geografica
se il tuo piano è creare una mappa di Mercatore. Se lo hai fatto finisci di proiettare e
riproiezione della griglia, perdendo i dettagli della lunghezza d'onda corta. Meglio usare -M e traccia la griglia
utilizzando una proiezione lineare con la stessa scala della proiezione di Mercatore desiderata (vedi GMT
Esempio 29). Per estrarre i dati nella regione -R-40/40/-70/-30 da world_grav.img.7.2, correre

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

Notare quello -V l'opzione ci dice che l'intervallo è stato regolato su
-R-40/40/-70.0004681551/-29.9945810754. Possiamo anche usare grdinfo per scoprire che il file della griglia
l'intestazione mostra la sua regione da essere -R0/80/0/67.9666667 Questo è l'intervallo di x,y che otterremo
da una proiezione sferica di Mercatore usando -R-40/40/-70.0004681551/-29.9945810754 e
-Jm1. Quindi, per prendere ship.lonlatgrav e usarlo per campionare merc_grav.nc, possiamo fare questo:

gmt set PROJ_ELLIPSOID Sfera

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

Si consiglia di utilizzare il metodo di cui sopra per proiettare e annullare la proiezione dei dati in tale
un'applicazione, perché allora c'è un solo passaggio di interpolazione (in grdtrack). Se uno
prima prova a convertire il file della griglia in lon,lat e poi campionalo, ce ne sono due
fasi di interpolazione (in conversione e in campionamento).

Per creare una griglia lon,lat dalla griglia sopra possiamo usare

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

In alcuni casi questo non sarà facile in quanto il -R nei due sistemi di coordinate potrebbero non allinearsi
bene. Quando ciò accade, possiamo anche utilizzare (anzi, potrebbe essere sempre meglio usare)

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

Per creare una mappa di Mercator della regione sopra, supponiamo che il nostro valore gmt.conf per
PROJ_LENGTH_UNIT è pollici. Quindi, poiché il file merc_grav.nc di cui sopra è proiettato con -Jm1i
è largo 80 pollici. Possiamo creare una mappa larga 8 pollici usando -Jx0.1i su qualsiasi programma cartografico
applicato a questa griglia (ad es. grdcontorno, grdimage, grdview), e poi per le sovrapposizioni che
lavorare in lon,lat (ad es. psxy, pcoast) possiamo usare quanto sopra aggiustato -R ed -Jm0.1 per ottenere
i due sistemi a combaciare.

Tuttavia, possiamo essere più intelligenti di così. Rendendosi conto che il file img di input aveva pixel 2.0
minuti (o controllando nx e ny con grdinfo merc_grav.nc) ci rendiamo conto che
merc_grav.nc ha utilizzato la piena risoluzione del file img e ha 2400 x 2039 pixel, e
a 8 pollici di larghezza questo è 300 pixel per pollice. Decidiamo che non abbiamo bisogno di così tanti e noi
sarà soddisfatto con 100 pixel per pollice, quindi vogliamo fare la media dei dati in 3 per 3
piazze. (Se vogliamo un diagramma di contorno, probabilmente sceglieremo di fare una media dei dati molto di più
(ad esempio, 6 per 6) per ottenere contorni lisci.) Poiché 2039 non è divisibile per 3 otterremo un
diverso regolato -R questa volta:

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

Questa volta troviamo che la regione modificata è -R-40/40/-70.023256525/-29.9368261101 e il
l'output è 800 per 601 pixel, una dimensione migliore per noi. Ora possiamo creare un artificiale
file di illuminazione per questo utilizzo gradiente:

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

e se abbiamo anche un file CPT chiamato "grav.cpt" possiamo creare una mappa in rilievo colorata
come questo:

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

Supponiamo di voler ottenere solo i valori dei dati vincolati da un file img, in lat/lon
coordinate. Allora corri img2grd con la -T2 opzione, usa grd2xyz per scaricare i valori, pipe
attraverso grep -v NaN per eliminare i NaN e passare attraverso mappaprogetto con l'inverso
proiezione come sopra.

Usa img2grdgmt online utilizzando i servizi onworks.net