Aceasta este comanda grdfiltergmt care poate fi rulată în furnizorul de găzduire gratuit OnWorks folosind una dintre multiplele noastre stații de lucru online gratuite, cum ar fi Ubuntu Online, Fedora Online, emulator online Windows sau emulator online MAC OS
PROGRAM:
NUME
grdfilter - Filtrați o grilă în domeniul spațiu (sau timp).
REZUMAT
grdfilter ingrid steag_distanță xlățime[/lățimea2][modificatori] outgrid [ creştere ] [ i|p|r ]
[ regiune ] [ ] [ [nivel] ] [ -f]
Notă: Nu este permis niciun spațiu între indicatorul opțiunii și argumentele asociate.
DESCRIERE
grdfilter va filtra un fișier grilă în domeniul temporal folosind una dintre convoluția selectată
sau filtre izotrope sau dreptunghiulare fără convoluție și calculează distanțele folosind carteziani
sau Geometrii sferice. Fișierul grilă de ieșire poate fi generat opțional ca subregiune
a intrării (prin -R) și/sau cu un nou increment (prin -I) sau înregistrare (prin -T). In acest
în felul acesta, se poate avea „spațiu suplimentar” în datele de intrare, astfel încât marginile să nu fie folosite și
ieșirea poate fi la o jumătate de lățime a marginilor de intrare. Dacă filtrul este trece jos,
atunci ieșirea poate fi eșantionată mai rar decât intrarea.
NECESARE ARGUMENTE
ingrid Fișierul grilă de puncte care trebuie filtrate. (Vezi mai jos FORMATELE DE FIȘIER GRILĂ).
-Dsteag_distanță
Distanţă pavilion spune modul în care grila (x,y) se referă la filtru lățime după cum urmează:
pavilion = p: grilă (px,py) cu lățime un număr impar de pixeli; distanțe carteziene.
pavilion = 0: grilă (x,y) aceleași unități ca lățime, distante carteziene.
pavilion = 1: grilă (x,y) în grade, lățime in kilometri, distante carteziene.
pavilion = 2: grilă (x,y) în grade, lățime în km, dx scalat de cos(y mijlociu), cartezian
distanțe.
Opțiunile de mai sus sunt cele mai rapide, deoarece permit doar calcularea matricei de greutate
o singura data. Următoarele trei opțiuni sunt mai lente, deoarece recalculează ponderi pentru fiecare
latitudine.
pavilion = 3: grilă (x,y) în grade, lățime în km, dx scalat după cosinus(y), carteziană
calculul distanței.
pavilion = 4: grilă (x,y) în grade, lățime în km, calculul distanței sferice.
pavilion = 5: grilă (x,y) în Mercator -Jm1 unități de imagini, lățime în km, distanță sferică
calcul.
-Fxlățime[/lățimea2][modificatori]
Setează tipul de filtru. Alegeți dintre filtrele de convoluție și non-convoluție. Folosiți oricare
cod de filtru x (enumerate mai jos) urmată de diametrul complet lățime. Aceasta dă un
filtru izotrop; anexează /lățimea2 pentru un filtru dreptunghiular (necesită -Dp or -D0).
Unele filtre permit argumente și modificatori opționali.
Filtrele de convoluție (și codurile acestora) sunt:
(b) Vagon: Toate greutățile sunt egale.
(c) Arcul cosinus: Greutățile urmează o curbă a arcului cosinus.
(g) Gaussian: Greutățile sunt date de funcția Gauss, unde lățime este de 6 de ori
sigma gaussiană convenţională.
(f) Personalizat: Greutățile sunt date de valorile precalculate în grila de greutate a filtrului
fişier greutate, care trebuie să aibă dimensiuni impare; cere de asemenea -D0 și distanța dintre ieșiri
trebuie să se potrivească cu spațierea de intrare sau să fie multipli întregi.
(o) Operator: Greutățile sunt date de valorile precalculate în grila de greutate a filtrului
fişier greutate, care trebuie să aibă dimensiuni impare; cere de asemenea -D0 și distanța dintre ieșiri
trebuie să se potrivească cu spațierea de intrare sau să fie multipli întregi. Se presupune că se însumează greutățile la
zero deci nu se va face nicio acumulare de sume de greutate și normalizare.
Filtrele non-convoluție (și codurile acestora) sunt:
(m) Mediană: returnează valoarea mediană. Pentru a selecta un alt cuantile anexează +qcuantilă in
intervalul 0-1 [Implicit este 0.5, adică mediana].
(p) Probabilitate maximă de probabilitate (un estimator de mod): returnează valoarea modală. Dacă mai mult
se găsește decât un mod, returnăm valoarea lor medie. Adăuga +l or +u dacă mai degrabă
doriți să returnați cea mai inferioară sau cea mai superioară dintre valorile modale.
(h) Mod histogramă (un alt estimator de mod): returnează valoarea modală ca centru
a vârfului dominant într-o histogramă. Adăugați /lățimea binului pentru a specifica binningul
interval. Utilizați modificatorul +c pentru a centra recipientele pe multipli de lățimea binului [Implicit are
marginile bin care sunt multipli ale lățimea binului]. Dacă se găsesc mai multe moduri noi
returnează valoarea lor medie. Adăuga +l or +u dacă mai degrabă doriți să returnați
cea mai inferioară sau cea mai superioară dintre valorile modale.
(l) Inferioară: returnează minimul tuturor valorilor.
(L) Inferioară: returnează doar valoarea minimă a tuturor valorilor pozitive.
(u) Superior: returnează maximul tuturor valorilor.
(U) Superior: returnează numai valorile maxime sau toate negative.
În cazul în care L|U este posibil ca nicio dată să treacă testul semnului inițial; în
în acest caz, filtrul va returna NaN.
-Goutgrid
outgrid este fișierul grilă de ieșire al filtrului. (Vezi mai jos FORMATELE DE FIȘIER GRILĂ).
OPTIONAL ARGUMENTE
-Ixinc[unitate][=|+][/yinc[unitate][=|+]]
x_inc [și opțional y_inc] este Incrementul de ieșire. Adăuga m pentru a indica arcul
minute, sau s pentru a indica secundele de arc. Dacă noul x_inc, y_inc NU sunt întregi
multiplii celor vechi (în datele de intrare), filtrarea va fi considerabil
Mai lent. [Implicit: La fel ca intrarea.]
-Ni|p|r
Determinați modul în care valorile NaN din grila de intrare afectează ieșirea filtrată: Adăugați i la
ignorați toate NaN în calculul valorii filtrate [Implicit], r este la fel ca i
cu excepția cazului în care nodul de intrare a fost NaN, atunci nodul de ieșire va fi setat la NaN (doar
se aplică dacă ambele grile sunt co-înregistrate) și p care va forţa valoarea filtrată
să fie NaN dacă se găsesc noduri-grilă cu valori NaN în interiorul cercului de filtru.
-R vest, est, sud și nord definește Regiunea punctelor de ieșire. [Mod implicit:
La fel ca și intrarea.]
-T Comutați înregistrarea nodului pentru grila de ieșire astfel încât să devină opusul
grila de intrare [Implicit oferă aceeași înregistrare ca și grila de intrare].
-V[nivel] (Mai mult ...)
Selectați nivelul de verbozitate [c].
-f[i|o]colinfo (Mai mult ...)
Specificați tipurile de date ale coloanelor de intrare și/sau de ieșire.
-^ 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.
GRID FILE FORMATE
În mod implicit, GMT scrie grila ca flotoare de precizie unică într-un netCDF de reclamație COARDS
tipul fisierului. Cu toate acestea, GMT este capabil să producă fișiere grilă în multe alte grile utilizate în mod obișnuit
formate de fișiere și facilitează, de asemenea, așa-numita „împachetare” a grilelor, scrierea în virgulă mobilă
date ca numere întregi de 1 sau 2 octeți. Pentru a specifica precizia, scara și decalajul, utilizatorul ar trebui
adăugați sufixul =id[/scară/compensa[/nan]], Unde id este un identificator de două litere al grilei
tipul și precizia și scară si compensa sunt factor de scară opțional și offset să fie
aplicat tuturor valorilor grilei și nan este valoarea folosită pentru a indica datele lipsă. In caz
cele două personaje id nu este prevăzut, ca în =/scară Decât a id=nf este asumat. Cand
citind grile, formatul este, în general, recunoscut automat. Dacă nu, același sufix
poate fi adăugat la numele fișierelor grilei de intrare. Vedea grdconvert și Secțiune grid-file-format al
Referințe tehnice GMT și carte de bucate pentru mai multe informații.
Când citește un fișier netCDF care conține mai multe grile, GMT va citi, în mod implicit, fișierul
prima grilă bidimensională care poate fi găsită în acel fișier. Pentru a convinge GMT să citească altul
variabilă multidimensională în fișierul grilă, anexați ?varname la numele fișierului, unde
varname este numele variabilei. Rețineți că poate fi necesar să scăpați de sensul special
of ? în programul shell, punând o bară oblică inversă în fața acestuia sau plasând
nume de fișier și sufix între ghilimele sau ghilimele duble. The ?varname se poate folosi și sufixul
pentru grilele de ieșire să specifice un nume de variabilă diferit de cel implicit: „z”. Vedea
grdconvert și Modificatori de secțiuni pentru CF și format de fișier grilă din Tehnica GMT
Referințe și Carte de bucate pentru mai multe informații, în special despre cum să citiți îmbinările de 3-,
Grile 4 sau 5-dimensionale.
GEOGRAFICE AND TIMP COORDONATE
Când tipul de grilă de ieșire este netCDF, coordonatele vor fi etichetate „longitudine”,
„latitudine” sau „timp” pe baza atributelor datelor de intrare sau ale grilei (dacă există) sau pe
-f or -R Opțiuni. De exemplu, ambele -f0x -f1t si -R90w/90e/0t/3t va avea ca rezultat a
grila de longitudine/timp. Când coordonatele x, y sau z este timpul, aceasta va fi stocată în grilă
ca timp relativ de la epocă, așa cum este specificat de TIME_UNIT și TIME_EPOCH în gmt.conf fişier
sau pe linia de comandă. In plus unitate va indica atributul variabilei de timp
atât această unitate cât şi epocă.
EXEMPLE
Să presupunem că north_pacific_etopo5.nc este un fișier de batimetrie de 5 minute de la 140E la 260E
și de la 0N la 50N și doriți să găsiți medianele valorilor pe o rază de 300 km (600 km
lățimea completă) a punctelor de ieșire, pe care le alegeți să fie de la 150E la 250E și de la 10N la 40N,
și doriți valorile de ieșire la fiecare 0.5 grade. Folosind calcule de distanță sferică,
ai nevoie:
gmt grdfilter north_pacific_etopo5.nc -Gfiltered_pacific.nc -Fm600 \
-D4 -R150/250/10/40 -I0.5 -V
Dacă am dorit în schimb un rezultat de trecere înaltă, atunci se poate efectua trecerea jos corespunzătoare
filtrați folosind un interval grosier de grilă, deoarece grdfilter va reeșantiona rezultatul la același
rezoluție ca grilă de intrare, astfel încât să putem calcula reziduurile, de exemplu,
gmt grdfilter north_pacific_etopo5.nc -Gresidual_pacific.nc -Fm-600 \
-D4 -R150/250/10/40 -I0.5 -V
Aici, grila residual_pacific.nc va avea aceeași rezoluție de 5 minute ca cea originală.
Pentru a filtra setul de date în ripples.nc folosind un filtru gaussian anizotrop personalizat exp
(-0.5*r^2) ale cărei distanțe r față de centru este dată de (2x^2 + y^2 -2xy)/6, cu major
axa la un unghi de 63 de grade cu orizontala, incercati
gmt grdmath -R-10/10/-10/10 -I1 X 2 POW 2 MUL Y 2 POW ADD XY MUL 2 MUL \
SUB 6 DIV NEG 2 DIV EXP DUP SUM DIV = gfilter.nc
gmt grdfilter ripples.nc -Ffgfilter.nc -D0 -Gsmooth.nc -V
LIMITAREA
Pentru a utiliza -D5, grila Mercator de intrare trebuie creată de img2mercgrd folosind -C
astfel încât originea valorilor y este Ecuatorul (adică, x = y = 0 corespunde lon =
lat = 0).
Utilizați grdfiltergmt online folosind serviciile onworks.net