grdfiltergmt - En ligne dans le Cloud

PROGRAMME:

Nom

grdfilter - Filtrer une grille dans le domaine spatial (ou temporel)

SYNOPSIS

grdfiltre ingrid distance_flag xlargeur[/largeur2][modificateurs] grille [ incrément ] [ i|p|r ]
[ région ] [ ] [ [niveau] ] [ -f]

Remarque : Aucun espace n'est autorisé entre l'indicateur d'option et les arguments associés.

DESCRIPTION

grdfiltre filtrera un fichier de grille dans le domaine temporel en utilisant l'une des convolutions sélectionnées
ou filtres isotropes ou rectangulaires sans convolution et calcul des distances à l'aide de cartes cartésiennes
ou géométries sphériques. Le fichier de grille de sortie peut éventuellement être généré en tant que sous-région
de l'entrée (via -R) et/ou avec un nouvel incrément (via -I) ou l'inscription (via -T). Dans ce
manière, on peut avoir « un espace supplémentaire » dans les données d'entrée afin que les bords ne soient pas utilisés et
la sortie peut être à moins d'une demi-largeur des bords d'entrée. Si le filtre est passe-bas,
alors la sortie peut être moins fréquemment échantillonnée que l'entrée.

REQUIS ARGUMENTS

ingrid Le fichier de grille des points à filtrer. (Voir FORMATS DE FICHIERS GRILLE ci-dessous).

-Ddistance_flag
Tarifs Distance drapeau indique comment la grille (x,y) se rapporte au filtre largeur comme suit:

drapeau = p : grille (px,py) avec largeur un nombre impair de pixels ; Distances cartésiennes.

drapeau = 0 : grille (x,y) mêmes unités que largeur, distances cartésiennes.

drapeau = 1 : grille (x,y) en degrés, largeur en kilomètres, distances cartésiennes.

drapeau = 2 : grille (x,y) en degrés, largeur en km, dx mis à l'échelle par cos(y du milieu), cartésien
les distances.

Les options ci-dessus sont les plus rapides car elles permettent de calculer uniquement la matrice de poids
une fois que. Les trois options suivantes sont plus lentes car elles recalculent les poids pour chaque
latitude.

drapeau = 3 : grille (x,y) en degrés, largeur en km, dx mis à l'échelle par cosinus(y), cartésien
calcul des distances.

drapeau = 4 : grille (x,y) en degrés, largeur en km, Calcul de distance sphérique.

drapeau = 5 : grille (x,y) dans Mercator -Jm1 unités img, largeur en km, Distance sphérique
calcul.

-Fxlargeur[/largeur2][modificateurs]
Définit le type de filtre. Choisissez parmi les filtres de convolution et de non-convolution. Utilisez n'importe quel
code de filtre x (énuméré ci-dessous) suivi du diamètre complet largeur. Cela donne un
filtre isotrope; ajouter /largeur2 pour un filtre rectangulaire (nécessite -Dp or -D0).
Certains filtres autorisent des arguments et des modificateurs facultatifs.

Les filtres de convolution (et leurs codes) sont :

(b) Wagon couvert : Tous les poids sont égaux.

(c) Arc en cosinus : les poids suivent une courbe en arc en cosinus.

(g) Gaussienne : les poids sont donnés par la fonction gaussienne, où largeur est 6 fois
le sigma gaussien conventionnel.

(f) Personnalisé : les poids sont donnés par les valeurs précalculées dans la grille de poids du filtre
filet poids, qui doit avoir des dimensions impaires ; nécessite également -D0 et espacement de sortie
doit correspondre à l'espacement des entrées ou être des multiples entiers.

(o) Opérateur : Les poids sont donnés par les valeurs précalculées dans la grille de poids du filtre
filet poids, qui doit avoir des dimensions impaires ; nécessite également -D0 et espacement de sortie
doit correspondre à l'espacement des entrées ou être des multiples entiers. Les poids sont supposés totaliser
zéro donc aucune accumulation de sommes de poids et de normalisation ne sera effectuée.

Les filtres sans convolution (et leurs codes) sont :

(m) Médiane : renvoie la valeur médiane. Pour sélectionner un autre quantile, ajoutez +qquantile in
la plage 0-1 [La valeur par défaut est 0.5, c'est-à-dire la médiane].

(p) Probabilité de vraisemblance maximale (un estimateur de mode) : renvoie la valeur modale. Si plus
plus d'un mode est trouvé, nous renvoyons leur valeur moyenne. Ajouter +l or +u si tu préfères
voulez retourner la plus basse ou la plus haute des valeurs modales.

(h) Mode histogramme (autre estimateur de mode) : renvoie la valeur modale comme centre
du pic dominant dans un histogramme. Ajouter /largeur de caisse pour spécifier le binning
intervalle. Utiliser le modificateur +c pour centrer les bacs sur des multiples de largeur de caisse [La valeur par défaut a
bords de bac qui sont des multiples de largeur de caisse]. Si plus d'un mode est trouvé, nous
renvoient leur valeur moyenne. Ajouter +l or +u si vous souhaitez plutôt retourner le
la plus basse ou la plus haute des valeurs modales.

(l) Inférieur : renvoie le minimum de toutes les valeurs.

(L) Inférieur : renvoie le minimum de toutes les valeurs positives uniquement.

(u) Upper : renvoie le maximum de toutes les valeurs.

(U) Supérieur : renvoie uniquement les valeurs maximales ou toutes les valeurs négatives.

Dans le cas d' L|U il est possible qu'aucune donnée ne passe le test de signe initial ; dans
dans ce cas, le filtre renverra NaN.

-Ggrille
grille est le fichier de grille de sortie du filtre. (Voir FORMATS DE FICHIERS GRILLE ci-dessous).

EN OPTION ARGUMENTS

-Ixinc[unité][=|+][/yinc[unité][=|+]]
x_inc [et éventuellement y_inc] est l'incrément de sortie. Ajouter m pour indiquer l'arc
minutes, ou s pour indiquer les secondes d'arc. Si le nouveau x_inc, y_inc ne sont PAS des entiers
multiples des anciens (dans les données d'entrée), le filtrage sera considérablement
Ralentissez. [Par défaut : Identique à l'entrée.]

-Ni|p|r
Déterminez comment les valeurs NaN de la grille d'entrée affectent la sortie filtrée : i à
ignorer tous les NaN dans le calcul de la valeur filtrée [Par défaut], r est identique à i
sauf si le nœud d'entrée était NaN, le nœud de sortie sera défini sur NaN (uniquement
s'applique si les deux réseaux sont co-enregistrés), et p qui forcera la valeur filtrée
être NaN si des nœuds de grille avec des valeurs NaN sont trouvés à l'intérieur du cercle de filtre.

-R ouest, est, sud et nord définit la région des points de sortie. [Défaut:
Identique à l'entrée.]

-T Basculez l'enregistrement du nœud pour la grille de sortie de manière à devenir l'opposé de
la grille d'entrée [Par défaut donne le même enregistrement que la grille d'entrée].

-V[niveau] (plus ...)
Sélectionnez le niveau de verbosité [c].

-f[je|o]colinfo (plus ...)
Spécifiez les types de données des colonnes d'entrée et/ou de sortie.

-^ or juste -
Imprime un court message sur la syntaxe de la commande, puis quitte (REMARQUE : sous Windows
utiliser juste -).

-+ or juste +
Imprimez un message d'utilisation détaillé (aide), y compris l'explication de tout
option spécifique au module (mais pas les options communes GMT), puis se ferme.

-? or aucune arguments
Imprimez un message d'utilisation (aide) complet, y compris l'explication des options, puis
sorties.

--version
Imprimer la version GMT et quitter.

--show-datadir
Affichez le chemin complet vers le répertoire de partage GMT et quittez.

GRID DOSSIER FORMATS

Par défaut, GMT écrit la grille sous forme de flottants en simple précision dans un netCDF de plainte COARDS
format de fichier. Cependant, GMT est capable de produire des fichiers de grille dans de nombreux autres grilles couramment utilisées
formats de fichiers et facilite également ce que l'on appelle "l'emballage" des grilles, l'écriture en virgule flottante
les données sous forme d'entiers de 1 ou 2 octets. Pour spécifier la précision, l'échelle et le décalage, l'utilisateur doit
ajouter le suffixe =id[/en échelon/compenser[/nan]], où id est un identifiant à deux lettres de la grille
type et précision, et en échelon et compenser sont un facteur d'échelle et un décalage facultatifs à
appliqué à toutes les valeurs de la grille, et nan est la valeur utilisée pour indiquer les données manquantes. Au cas où
les deux personnages id n'est pas fourni, comme dans =/en échelon qu'un id=nf est assumé. Lorsque
grilles de lecture, le format est généralement reconnu automatiquement. Sinon, le même suffixe
peuvent être ajoutés aux noms de fichiers de grille d'entrée. Voir grdconvertir et Section grid-file-format du
Référence technique GMT et livre de recettes pour plus d'informations.

Lors de la lecture d'un fichier netCDF contenant plusieurs grilles, GMT lira, par défaut, le
première grille bidimensionnelle que l'on peut trouver dans ce fichier. Amener GMT à lire un autre
variable multidimensionnelle dans le fichier de grille, ajoutez ?varname au nom du fichier, où
varname est le nom de la variable. Notez que vous devrez peut-être échapper au sens spécial
of ? dans votre programme shell en plaçant une barre oblique inverse devant celui-ci, ou en plaçant le
nom de fichier et suffixe entre guillemets ou guillemets doubles. Les ?varname le suffixe peut également être utilisé
pour les grilles de sortie, spécifiez un nom de variable différent de la valeur par défaut : "z". Voir
grdconvertir et Sections modifiers-for-CF et grid-file-format du GMT Technical
Référence et livre de recettes pour plus d'informations, en particulier sur la façon de lire les épissures de 3-,
Grilles à 4 ou 5 dimensions.

GÉOGRAPHIQUE ET COORDONNEES

Lorsque le type de grille de sortie est netCDF, les coordonnées seront étiquetées "longitude",
« latitude » ou « heure » ​​en fonction des attributs des données d'entrée ou de la grille (le cas échéant) ou de la
-f or -R option. Par exemple, les deux -f0x -f1t et -R90w/90e/0t/3t se traduira par un
grille de longitude/temps. Lorsque la coordonnée x, y ou z est le temps, elle sera stockée dans la grille
comme temps relatif depuis l'époque comme spécifié par TIME_UNIT et TIME_EPOCH dans le gmt.conf filet
ou sur la ligne de commande. De plus, le unité l'attribut de la variable de temps indiquera
à la fois cette unité et cette époque.

EXEMPLES

Supposons que north_pacific_etopo5.nc soit un fichier de bathymétrie de 5 minutes de 140E à 260E
et 0N à 50N, et vous voulez trouver les médianes des valeurs dans un rayon de 300km (600km
pleine largeur) des points de sortie, que vous choisissez entre 150E à 250E et 10N à 40N,
et vous voulez les valeurs de sortie tous les 0.5 degré. En utilisant des calculs de distance sphérique,
vous avez besoin:

gmt grdfilter north_pacific_etopo5.nc -Gfiltered_pacific.nc -Fm600 \
-D4 -R150/250/10/40 -I0.5 -V

Si nous voulions plutôt un résultat passe-haut, alors on peut effectuer le passe-bas correspondant
filtrer en utilisant un intervalle de grille grossier car grdfilter rééchantillonnera le résultat au même
résolution comme grille d'entrée afin que nous puissions calculer les résidus, par exemple,

gmt grdfilter north_pacific_etopo5.nc -Gresidual_pacific.nc -Fm-600 \
-D4 -R150/250/10/40 -I0.5 -V

Ici, la grille résiduelle_pacific.nc aura la même résolution de 5 minutes que l'original.

Pour filtrer l'ensemble de données dans ripples.nc à l'aide d'un filtre gaussien anisotrope personnalisé exp
(-0.5*r^2) dont les distances r du centre sont données par (2x^2 + y^2 -2xy)/6, avec
axe à un angle de 63 degrés avec l'horizontale, essayez

gmt grdmath -R-10/10/-10/10 -I1 X 2 POW 2 MUL Y 2 POW AJOUTER XY MUL 2 MUL \
SUB 6 DIV NEG 2 DIV EXP DUP SUM DIV = gfilter.nc
gmt grdfilter ondulations.nc -Ffgfilter.nc -D0 -Gsmooth.nc -V

LIMITATIONS

Pour utiliser la -D5 option la grille de Mercator d'entrée doit être créée par img2mercgrd en utilisant le -C
option de sorte que l'origine des valeurs y est l'équateur (c'est-à-dire que x = y = 0 correspond à lon =
latitude = 0).

Utilisez grdfiltergmt en ligne en utilisant les services onworks.net