dimfiltergmt - En ligne dans le Cloud

PROGRAMME:

Nom

dimfilter - Filtrage directionnel de fichiers quadrillés 2D dans le domaine spatial (ou temporel)

SYNOPSIS

dimfiltre fichier_entrée.nc distance_flag [mode] fichier_sortie.nc
[ cols ] [ incrément ] [ région ] [ ] [ [niveau] ] [ -f]

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

DESCRIPTION

dimfiltre va filtrer un .Caroline du Nord fichier dans le domaine de l'espace (ou du temps) en divisant le
filtrer le cercle dans n_secteurs, en appliquant l'une des convolutions primaires sélectionnées ou
filtres de non-convolution à chaque secteur, et en choisissant le résultat final en fonction de la
filtre secondaire sélectionné. Il calcule les distances en utilisant des géométries cartésiennes ou sphériques.
La sortie .Caroline du Nord peut éventuellement être généré en tant que sous-région de l'entrée et/ou avec un
nouvelle -Iincrément. De cette façon, on peut avoir un "espace supplémentaire" dans les données d'entrée afin qu'il y ait
n'y aura pas d'effets de bord pour la grille de sortie. Si le filtre est passe-bas, la sortie
peut être moins fréquemment échantillonné que l'entrée. -Q est pour le mode d'analyse d'erreur et seulement
requiert le nombre total de colonnes dans le fichier d'entrée, qui contient les
profondeurs. Enfin, il faut savoir que dimfiltre ne produira pas une sortie lisse comme d'autres
les filtres spatiaux le font car il renvoie une médiane minimale sur N médianes de N secteurs. Les
la sortie peut être approximative à moins que les données d'entrée soient sans bruit. Ainsi, un filtrage supplémentaire
(par exemple, gaussien via grdfiltre) des données filtrées DiM est généralement recommandée.

REQUIS ARGUMENTS

fichier_entrée.nc
La grille de données à filtrer.

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

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 moyen), distances cartésiennes.

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.

-F[mode]
Définit le type de filtre principal. Choisissez parmi les filtres de convolution et de non-convolution.
Ajouter le code du filtre suivi du diamètre complet largeur. Convolution disponible
les filtres 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.

Les filtres sans convolution sont :

(m) Médiane : renvoie la valeur 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 - ou + au filtre
width si vous souhaitez plutôt retourner la plus petite ou la plus grande des valeurs modales.

-N
Définit le type de filtre secondaire et le nombre de secteurs en nœud papillon. n_secteurs doit
être un entier et supérieur à 0. Quand n_secteurs est mis à 1, le filtre secondaire est
Inefficace. Les filtres secondaires disponibles sont :

(l) Inférieur : Renvoie le minimum de toutes les valeurs filtrées.

(u) Supérieur : renvoie le maximum de toutes les valeurs filtrées.

(a) Moyenne : renvoie la moyenne de toutes les valeurs filtrées.

(m) Médiane : renvoie la médiane de toutes les valeurs filtrées.

(p) Mode : renvoie le mode de toutes les valeurs filtrées.

-Gfichier_sortie.nc
fichier_sortie.nc est la sortie du filtre.

EN OPTION ARGUMENTS

-I x_inc [et éventuellement y_inc] est l'incrément de sortie. Ajouter m pour indiquer les minutes,
or c pour indiquer les secondes. Si le nouveau x_inc, y_inc ne sont PAS des multiples entiers de
anciens (dans les données d'entrée), le filtrage sera considérablement plus lent. [Par défaut : Idem
comme entrée.]

-R ouest, est, sudet 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].

-Qcols cols est le nombre total de colonnes dans le fichier de table de texte d'entrée. Pour ce mode,
il s'attend à lire des profondeurs composées de plusieurs colonnes. Chaque colonne représente un
grille filtrée avec une largeur de filtre, qui peut être obtenue par grd2xyz -Z. Le résultat
sera médiane, MAD et moyenne. Ainsi, la colonne avec les médianes est utilisée pour générer
la composante régionale et la colonne avec les MAD sont utilisées pour effectuer l'erreur
analyse.

-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_dbdb5.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 complet
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é. Pour éviter que les médianes ne soient biaisées par
le plan incliné, vous voulez diviser le cercle de filtre en 6 secteurs et choisir le
valeur la plus basse parmi 6 médianes. En utilisant les calculs de distance sphérique, vous avez besoin de :

gmt dimfilter north_pacific_dbdb5.nc -Gfiltered_pacific.nc -Fm600 -D4 \
-Nl6 -R150/250/10/40 -I0.5 -V

Supposons que cape_verde.nc soit un fichier de bathymétrie de 0.5 minute de 32W à 15W et de 8N à
25N, et vous souhaitez supprimer des entités à petite échelle afin de définir une houle dans un
zone s'étendant de 27.5W à 20.5W et 12.5N à 19.5N, et vous voulez la valeur de sortie chaque
Deux minutes. En utilisant les calculs de distance cartésienne, vous avez besoin de :

gmt dimfiltre cape_verde.nc -Gt.nc -Fm220 -Nl8 -D2 -R-27.5/-20.5/12.5/19.5 -I2m -V
gmt grdfilter t.nc -Gcape_swell.nc -Fg50 -D2 -V

Supposons que vous ayez trouvé une plage de largeurs de filtre pour une zone donnée et que vous ayez filtré le
données bathymétriques données en utilisant la gamme de largeurs de filtre (par exemple, f100.nc f110.nc f120.nc
f130.nc), et vous souhaitez définir une tendance régionale à l'aide de la plage de largeurs de filtre, et
vous souhaitez obtenir des estimations de l'écart absolu médian (MAD) à chaque point de données. Alors vous
devra faire :

gmt grd2xyz f100.nc -Z > f100.d
gmt grd2xyz f110.nc -Z > f110.d
gmt grd2xyz f120.nc -Z > f120.d
gmt grd2xyz f130.nc -Z > f130.d
coller f100.d f110.d f120.d f130.d > profondeurs.d
gmt dimfilter depths.d -Q4 > output.z

LIMITATIONS

Lorsque vous travaillez avec des grilles géographiques (lat, lon), les trois filtres de convolution (boxcar,
arc cosinus et gaussien) normalisera correctement les poids de filtre pour la variation de
la taille de la grille avec la latitude, et déterminer correctement quels nœuds sont nécessaires pour le
convolution lorsque le filtre "cercle" traverse une frontière périodique (0-360) ou contient un
pôle géographique. Cependant, les filtres spatiaux, tels que les filtres médians et de mode, n'utilisent pas
poids et ne doit donc être utilisé que sur des grilles cartésiennes (ou à de très basses latitudes).
Si vous souhaitez appliquer de tels filtres spatiaux, vous devez projeter vos données sur une surface égale
projection et exécuter dimfilter sur la grille cartésienne résultante.

SCÉNARIO MODÈLE

Le dim.template.sh est un script shell squelette qui peut être utilisé pour configurer un DiM complet
analyse, y compris l'analyse MAD.

RÉFÉRENCE

Kim, S.-S., et Wessel, P. (2008), Filtrage médian directionnel pour les résidus régionaux
Séparation de la bathymétrie, Géochimie. Géophysique. Géosyste., 9, Q03005,
doi: 10.1029/2007GC001850.

Utilisez dimfiltergmt en ligne en utilisant les services onworks.net