Questo è il comando surfacegmt che può essere eseguito nel provider di hosting gratuito OnWorks utilizzando una delle nostre molteplici workstation online gratuite come Ubuntu Online, Fedora Online, emulatore online Windows o emulatore online MAC OS
PROGRAMMA:
NOME
superficie: dati della tabella della griglia utilizzando spline di curvatura continua a tensione regolabile
SINOSSI
superficie [ tavolo ] outputfile.nc incremento regione [ proporzioni ] [ limite_convergenza ] [
labbassarla ] [ -Lusuperiore ] [ max_iterazioni ] [ ] [ raggio_di ricerca[m|s]] [ [i|b]fattore_tensione
] [ [ livello]] [ fattore_di_relax eccessivo ] [ -a] [ -bi] [ -du] [
-f] [ -h] [ -i] [ -:[i|o]]
Nota: Non è consentito alcuno spazio tra il flag dell'opzione e gli argomenti associati.
DESCRIZIONE
superficie legge triple spaziate in modo casuale (x,y,z) dallo standard input [or tavolo] e produce
un file di griglia binario di valori grigliati z(x,y) risolvendo:
(1 - T) * L (L (z)) + T * L (z) = 0
dove T è un fattore di tensione compreso tra 0 e 1 e L indica l'operatore laplaciano. T = 0
fornisce la soluzione di "curvatura minima" che è equivalente a SuperMISP e ISM
Pacchetti. La curvatura minima può causare oscillazioni indesiderate e falsi massimi locali
minimi (vedi Smith e Wessel, 1990), e potresti voler usare T > 0 per sopprimerli
effetti. L'esperienza suggerisce che T ~ 0.25 di solito sembra buono per potenziali dati sul campo e T
dovrebbe essere maggiore (T ~ 0.35) per i dati di topografia ripida. T = 1 dà una superficie armonica (n
massimi o minimi sono possibili tranne che nei punti dati di controllo). Si consiglia di
l'utente preelabora i dati con blocco, mediana, o modalità di blocco per evitare spaziali
creare alias ed eliminare i dati ridondanti. È possibile imporre limiti inferiori e/o superiori al file
soluzione. Questi possono essere inseriti sotto forma di valore fisso, griglia con valori o semplicemente
essere i valori minimo/massimo dei dati di input. Le condizioni al contorno naturali vengono applicate al
bordi, ad eccezione dei dati geografici con intervallo di 360 gradi in cui applichiamo il confine periodico
condizioni nella direzione della longitudine.
OBBLIGATORIO ARGOMENTI
-Goutputfile.nc
Nome del file di output. L'output è un binario 2-D .nc file. Si noti che la griglia più piccola
la dimensione deve essere almeno 4.
-Ixinc[unità][=|+][/sì[unità][=|+]]
x_inc [e facoltativamente y_inc] è la spaziatura della griglia. Facoltativamente, aggiungi un suffisso
modificatore. geografico (gradi) coordinate: Aggiungi m per indicare i minuti d'arco o s
per indicare i secondi d'arco. Se una delle unità e, f, k, M, n or u è allegato
si assume invece che l'incremento sia espresso in metri, piedi, km, miglia, nautiche
miglio o piede di rilevamento statunitense, rispettivamente, e verrà convertito nell'equivalente
gradi di longitudine alla latitudine media della regione (la conversione dipende da
PROJ_ELLIPSOID). Se /y_inc è dato ma impostato a 0 sarà resettato uguale a x_inc;
altrimenti verrà convertito in gradi di latitudine. Tutti coordinate: Se = is
aggiunto poi il corrispondente max x (est) o puoi y (nord) può essere leggermente regolato
per adattarsi esattamente all'incremento dato [per impostazione predefinita l'incremento può essere regolato
leggermente per adattarsi al dominio dato]. Infine, invece di dare un incremento potresti
specificare la numero of nodi desiderato aggiungendo + all'intero fornito
discussione; l'incremento viene quindi ricalcolato dal numero di nodi e dal
dominio. Il valore di incremento risultante dipende dal fatto che tu abbia selezionato a
griglia registrata sulla griglia o sui pixel; vedere App-file-formati per i dettagli.
Nota: se -Rgrdfile viene utilizzato allora la spaziatura della griglia è già stata inizializzata; utilizzo
-I per sovrascrivere i valori.
-R[unità]xmin/Xmax/ymin/ymax[R] (Di Più ...)
Specificare la regione di interesse.
OPTIONAL ARGOMENTI
tavolo Uno o più ASCII (o binari, vedi -bi[coli][Digitare]) file di tabelle dati contenenti a
numero di colonne di dati. Se non vengono fornite tabelle, leggiamo dallo standard input.
-Aproporzioni
Proporzioni. Se lo si desidera, è possibile aggiungere alle equazioni l'anisotropia della griglia. accedere
proporzioni, dove dy = dx / proporzioni mette in relazione le dimensioni della griglia. [Predefinito =
1 presuppone una griglia isotropa.]
-Climite_convergenza[%]
Limite di convergenza. Si presuppone che l'iterazione converga al massimo assoluto
la variazione in qualsiasi valore della griglia è inferiore a limite_convergenza. (Unità uguali ai dati z
unità). In alternativa, fornire il limite in percentuale della deviazione rms aggiungendo %.
[Il valore predefinito è scalato su 1e-4 della deviazione quadratica media dei dati da a
piano più adatto (dei minimi quadrati).]. Questo è il limite di convergenza finale al
spaziatura della griglia desiderata; per griglie intermedie (più grossolane) la convergenza effettiva
il limite viene scalato dal moltiplicatore della spaziatura della griglia.
-LLabbassarla e -Lusuperiore
Imporre limiti alla soluzione di output. labbassarla imposta il limite inferiore. abbassarla può essere il file
nome di un file griglia con valori limite inferiore, un valore fisso, d impostare al minimo
valore di input, o u per non vincolato [Predefinito]. usuperiore imposta il limite superiore e può
essere il nome di un file griglia con valori limite superiore, un valore fisso, d impostare su
valore massimo di input, o u per non vincolato [Predefinito]. File della griglia utilizzati per impostare il file
i limiti possono contenere NaN. In presenza di NaNs, il limite di un nodo mascherato con
NaN non è vincolato.
-Nmax_iterazioni
Numero di iterazioni. L'iterazione cesserà quando limite_convergenza viene raggiunto o
quando viene raggiunto il numero di iterazioni max_iterazioni. Questa è l'iterazione finale
limitare alla spaziatura della griglia desiderata; per griglie intermedie (più grossolane) l'effettivo
il limite di iterazione viene scalato dal moltiplicatore della spaziatura della griglia. [L'impostazione predefinita è 500.]
-Q Suggerire dimensioni della griglia che abbiano un massimo comun divisore altamente composito. Questo
consente alla superficie di utilizzare diversi passaggi intermedi nella soluzione, cedendo più velocemente
tempi di esecuzione e risultati migliori. Le taglie suggerite da -Q può essere ottenuto alterando
-R e / o -I. Puoi recuperare il -R e -I vuoi più tardi usando grdcampione or
grdcut sull'uscita di superficie.
-Sraggio_di ricerca[m|s]
Raggio di ricerca. accedere raggio_di ricerca nelle stesse unità dei dati x,y; aggiungere m indicare
minuti d'arco o s per secondi d'arco. Viene utilizzato per inizializzare la griglia prima del file
prima iterazione; non ne vale la pena a meno che il reticolo della griglia non sia primo e
non possono avere fasi regionali. [Predefinito = 0.0 e non viene effettuata alcuna ricerca.]
-T[i|b]fattore_tensione
Fattore di tensione[i]. Questi devono essere compresi tra 0 e 1. La tensione può essere utilizzata in
soluzione interna (sopra l'equazione, dove sopprime le oscillazioni spurie) e
nelle condizioni al contorno (dove tende ad appiattire la soluzione avvicinandosi al
bordi). Utilizzando zero per entrambi i valori si ottiene una superficie con curvatura minima con libero
bordi, cioè una spline bicubica naturale. Utilizzo -Tifattore_tensione per impostare l'interno
tensione, e -Tbfattore_tensione per impostare la tensione al contorno. Se non anteponi i or
b, entrambi verranno impostati sullo stesso valore. [Default = 0 per entrambi dà il minimo
soluzione di curvatura.]
-V[ livello] (Di Più ...)
Seleziona il livello di verbosità [c]. -V3 riporterà la convergenza dopo ogni iterazione;
-V riporteranno solo dopo la convergenza di ciascuna griglia regionale.
-Zfattore_di_relax eccessivo
Fattore di rilassamento eccessivo. Questo parametro viene utilizzato per accelerare la convergenza; è
un numero compreso tra 1 e 2. Un valore pari a 1 ripete esattamente le equazioni e lo farà
garantire sempre una convergenza stabile. Valori più grandi sovrastimano l'incrementale
cambiamenti durante la convergenza e raggiungeranno una soluzione più rapidamente, ma potrebbero diventarlo
instabile. Se si utilizza un valore elevato per questo fattore, è una buona idea monitorare
ogni iterazione con il -Vl opzione. [Default = 1.4 converge rapidamente ed è quasi
sempre stabile.]
-acol=Nome[...] (Di Più ...)
Imposta associazioni di colonne spaziali col=Nome.
-bi[coli][T] (Di Più ...)
Seleziona input binario nativo. [Il valore predefinito è 3 colonne di input].
-dunessun dato (Di Più ...)
Sostituisci le colonne di input uguali nessun dato con NaN.
-f[io|o]colinfo (Di Più ...)
Specificare i tipi di dati delle colonne di input e/o output.
-h[io|o][n][+c][+d][+rosservazione][+rtitolo] (Di Più ...)
Salta o produce record di intestazione. Non utilizzato con dati binari.
-icols[l] [sscala][ooffset][,...] (Di Più ...)
Seleziona le colonne di input (0 è la prima colonna).
-:[io|o] (Di Più ...)
Scambia la prima e la seconda colonna su input e/o output.
-^ 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.
GRID VALORI PRECISIONE
Indipendentemente dalla precisione dei dati di input, i programmi GMT che creano file di griglia lo faranno
contenere internamente le griglie in array a virgola mobile a 4 byte. Questo è fatto per conservare la memoria
e inoltre la maggior parte se non tutti i dati reali possono essere archiviati utilizzando la virgola mobile a 4 byte
valori. I dati con maggiore precisione (vale a dire, valori a doppia precisione) perderanno questo
precisione una volta che GMT opera sulla griglia o scrive nuove griglie. Per limitare la perdita di
precisione durante l'elaborazione dei dati si dovrebbe sempre considerare la normalizzazione dei dati prima di
trattamento.
ESEMPI
Per grigliare 5 blocchi di gravità per 5 minuti si intende dai dati ASCII in hawaii_5x5.xyg, utilizzando un
fattore_tensione = 0.25, a limite_convergenza = 0.1 milligal, scrivendo il risultato in un file
chiamato hawaii_grd.nc e monitorando ogni iterazione, prova:
gmt superficie hawaii_5x5.xyg -R198/208/18/25 -I5m -Ghawaii_grd.nc -T0.25 -C0.1 -Vl
Utilizza Surfacegmt online utilizzando i servizi onworks.net