grdfftgmt - En ligne dans le Cloud

PROGRAMME:

Nom

grdfft - Effectue des opérations mathématiques sur des grilles dans le domaine des nombres d'onde (ou fréquence)

SYNOPSIS

grdfft ingrid [ ingrid2 ] fichier de sortie [ azimut ] [ niveau z ] [ [en échelon|g] ] [ [r|x|y][w[k]] ] [
[r|x|y]params ] [ [en échelon|g] ] [
[f|q|s|nx/ny][+a|d|h|l][+e|n|m][+tlargeur][+w[suffixe]][+z[p]] ] [ en échelon ] [ [niveau] ] [ -fg
]

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

DESCRIPTION

grdfft prendra la transformation de Fourier rapide avant 2D et effectuera une ou plusieurs
opérations mathématiques dans le domaine fréquentiel avant de se retransformer dans l'espace
domaine. Une option est fournie pour mettre à l'échelle les données avant d'écrire les nouvelles valeurs dans une sortie
déposer. Les dimensions horizontales de la grille sont supposées être en mètres. Géographique
les grilles peuvent être utilisées en spécifiant le -fg option qui met à l'échelle les degrés en mètres. Si tu as
grilles avec des dimensions en km, vous pouvez changer cela en mètres en utilisant grdédit ou à l'échelle
sortie avec mathmath.

REQUIS ARGUMENTS

ingrid Fichier de grille binaire 2-D à opérer. (Voir FORMATS DE FICHIERS GRILLE ci-dessous). Pour
opérations interspectrales, donnent également le deuxième fichier de grille ingrd2.

-Gfichier de sortie
Spécifiez le nom du fichier de grille de sortie ou de la table de spectre 1-D (voir -E). (Voir
FORMATS DE FICHIERS GRILLE ci-dessous).

EN OPTION ARGUMENTS

-Aazimut
Prenez la dérivée directionnelle dans le azimut direction mesurée en degrés CW
du nord.

-Cniveau z
Vers le haut (pour niveau z > 0) ou vers le bas (pour niveau z < 0) continuer le champ niveau z
mètres.

-RÉ[en échelon|g]
Différencier le champ, c'est-à-dire prendre d(field)/dz. Cela équivaut à multiplier
par kr dans le domaine fréquentiel (kr est le nombre d'onde radiale). Ajouter une échelle à
multiplier par (kr * en échelon) au lieu. Sinon, ajoutez g pour indiquer que votre
les données sont les hauteurs du géoïde en mètres et les sorties devraient être des anomalies de gravité en mGal.
[La valeur par défaut n'est pas d'échelle].

-E[r|x|y][w[k]]
Estimer le spectre de puissance dans la direction radiale [r]. Endroit x or y juste après
-E pour calculer le spectre dans la direction x ou y à la place. Aucun fichier de grille n'est
créé. Si une grille est donnée alors f (c'est-à-dire, fréquence ou nombre d'onde), puissance[f],
et 1 écart type de puissance[f] sont écrits dans le fichier défini par -G [stdout]. Si
deux grilles sont données nous écrivons f et 8 quantités : Xpuissance[f], Ypuissance[f], cohérente
puissance[f], puissance de bruit[f], phase[f], admittance[f], gain[f], cohérence[f]. Chaque
la quantité est suivie de sa propre estimation d'erreur de développement 1std, donc la sortie est de 17
colonnes de large. Ajouter w pour écrire la longueur d'onde au lieu de la fréquence. Si votre grille est
géographique, vous pouvez ajouter k pour mettre à l'échelle les longueurs d'onde du mètre [Par défaut] à
km.

-F[r|x|y]params
Filtrez les données. Endroit x or y juste après -F à filtrer x or y direction seulement;
la valeur par défaut est isotrope [r]. Choisissez entre un passe-bande conique cosinus, un gaussien
filtre passe-bande ou un filtre passe-bande de Butterworth.

Cosinus-cône :
Spécifiez quatre longueurs d'onde lc/lp/hp/hc dans les bonnes unités (voir -fg) pour concevoir un
filtre passe-bande : longueurs d'onde supérieures à lc ou moins que hc sera coupé,
longueurs d'onde supérieures à lp et moins de hp sera passé, et les longueurs d'onde
entre les deux sera conique en cosinus. Par exemple, -F1000000 / 250000 / 50000 / 10000 -fg
passe-bande, coupant les longueurs d'onde > 1000 km et < 10 km, passant
longueurs d'onde comprises entre 250 km et 50 km. Pour faire un filtre passe-haut ou passe-bas,
donner des tirets (-) pour hp/hc or lc/lp. Par exemple, -Fx-/-/50/10 passe-bas x,
longueurs d'onde de passage > 50 et longueurs d'onde de rejet < 10. -Fy1000/250/-/-
passera haut y, les longueurs d'onde de passage < 250 et les longueurs d'onde de rejet > 1000.

Gaussienne passe-bande:
Ajouter lo/hi, les deux longueurs d'onde dans les bonnes unités (voir -fg) pour concevoir un
filtre passe-bande. Aux longueurs d'onde données, les poids du filtre gaussien seront
être de 0.5. Pour créer un filtre passe-haut ou passe-bas, donnez un trait d'union (-) pour le hi
or lo longueur d'onde, respectivement. Par exemple, -F-/30 passera les données à l'aide d'un
Filtre gaussien à demi-poids à 30, tandis que -F400/- passera le haut
revendre.

Butterworth passe-bande:
Ajouter lo/hi/de commander, les deux longueurs d'onde dans les bonnes unités (voir -fg) Et le
ordre de filtrage (un entier) pour concevoir un filtre passe-bande. Au donné
longueurs d'onde, les poids du filtre de Butterworth seront de 0.5. Pour faire un passe-haut
ou un filtre passe-bas, donnez un trait d'union (-) pour le hi or lo longueur d'onde,
respectivement. Par exemple, -F-/30/2 passera les données à l'aide d'un 2nd ordre
Filtre Butterworth, avec demi-poids à 30, tandis que -F400/-/2 passe-haut
les données.

-JE[en échelon|g]
Intégrez le champ, c'est-à-dire calculez intégral_over_z (champ * dz). C'est
équivalent à diviser par kr dans le domaine fréquentiel (kr est le nombre d'onde radiale).
Ajoutez une échelle pour diviser par (kr * en échelon) au lieu. Sinon, ajoutez g à
indique que votre ensemble de données est constitué d'anomalies de gravité en mGal et que la sortie doit être en géoïde
hauteurs en mètres. [La valeur par défaut n'est pas d'échelle].

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tlargeur][+w[suffixe]][+z[p]]
Choisissez ou renseignez-vous sur les dimensions de grille appropriées pour FFT et réglez en option
paramètres. Contrôlez la dimension FFT :
-Nf forcera la FFT à utiliser les dimensions réelles des données.

-Nq se renseignera sur les dimensions les plus appropriées, les signaler, puis continuer.

-Ns présentera une liste de dimensions facultatives, puis quittera.

-Nnx/ny fera FFT sur la taille du tableau nx/ny (doit être >= taille du fichier de grille). Défaut
choisit les dimensions >= données qui optimisent la vitesse et la précision de la FFT. Si FFT
dimensions > dimensions du fichier de grille, les données sont étendues et réduites à zéro.

Contrôlez la suppression de tendance des données : ajoutez des modificateurs pour supprimer une tendance linéaire :
+d: Supprimer la tendance, c'est-à-dire supprimer la tendance linéaire la mieux ajustée [Par défaut].

+a: ne supprime que la valeur moyenne.

+h: ne supprimer que la valeur médiane, c'est-à-dire 0.5 * (max + min).

+l: Laissez les données tranquilles.

Contrôler l'extension et la réduction des données : utilisez des modificateurs pour contrôler la façon dont l'extension
et tapering sont à effectuer :
+e étend le maillage en imposant une symétrie bord-point [Par défaut],

+m étend la grille en imposant une symétrie miroir de bord

+n désactive l'extension de données.

L'effilage est effectué du bord des données au bord de la grille FFT [100%]. Changer
ce pourcentage via +tlargeur. Quand +n est en vigueur, le tapering est appliqué
à la place des marges de données car aucune extension n'est disponible [0%].

Contrôle de l'écriture des résultats temporaires : pour une enquête détaillée, vous pouvez écrire le
la grille intermédiaire étant transmise à la FFT avant ; cela a probablement été
sans tendance, étendu par symétrie ponctuelle le long de tous les bords et effilé. Ajouter
+w[suffixe] à partir duquel le(s) nom(s) de fichier de sortie seront créés (c'est-à-dire, ingrid_prefix.ext)
[conique], où poste est votre extension de fichier. Enfin, vous pouvez enregistrer la grille complexe
produit par la FFT à terme en ajoutant +z. Par défaut, nous écrivons le réel et
composants imaginaires à ingrid_réel.poste et ingrid_image.poste. Ajouter p pour sauver
à la place la forme polaire de la magnitude et de la phase dans les fichiers ingrid_mag.poste et
ingrid_phase.poste.

-Sen échelon
Multipliez chaque élément par en échelon dans le domaine spatial (après le domaine fréquentiel
opérations). [La valeur par défaut est 1.0].

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

-fg Les grilles géographiques (dimensions de longitude, latitude) seront converties en mètres
via une approximation « Terre plate » en utilisant les paramètres actuels de l'ellipsoïde.

-^ 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.

GRID DISTANCE UNITÉS

Si la grille n'a pas de mètre comme unité horizontale, ajoutez +uunité au fichier d'entrée
nom à convertir de l'unité spécifiée en mètre. Si votre grille est géographique, convertissez
distances aux mètres en fournissant -fg à la place.

CONSIDERATIONS

Les grilles netCDF COARDS seront automatiquement reconnues comme géographiques. Pour les autres grilles
grilles géographiques où vous souhaitez convertir les degrés en mètres, sélectionnez -fg. Si les données
sont proches de l'un ou l'autre pôle, vous devriez envisager de projeter le fichier de grille sur un rectangle
système de coordonnées utilisation grdprojet

EXEMPLES

Continuer vers le haut les anomalies magnétiques du niveau de la mer dans le fichier mag_0.nc jusqu'à un niveau 800 m
au dessus du niveau de la mer:

gmt grdfft mag_0.nc -C800 -V -Gmag_800.nc

Pour transformer les hauteurs du géoïde en m (geoid.nc) sur une grille géographique en gravité à l'air libre
anomalies en mGal :

gmt grdfft geoïde.nc -Dg -V -Ggrav.nc

Transformer les anomalies de gravité en mGal (faa.nc) en déviations de la verticale (en
micro-radians) dans la direction 038, il faut d'abord intégrer la gravité pour obtenir le géoïde, puis
prendre la dérivée directionnelle, et enfin mettre les radians à l'échelle en micro-radians :

gmt grdfft faa.nc -Ig -A38 -S1e6 -V -Gdefl_38.nc

Les secondes dérivées verticales des anomalies gravimétriques sont liées à la courbure de la
champ. Nous pouvons les calculer en mGal/m^2 en différenciant deux fois :

gmt grdfft gravité.nc -D -D -V -Ggrav_2nd_derivative.nc

Pour calculer des estimations interspectrales pour les grilles bathymétriques et gravimétriques co-enregistrées, et
rapporter le résultat en fonction des longueurs d'onde en km, essayer

gmt grdfft bathymetry.nc gravité.grd -Ewk -fg -V > cross_spectra.txt

Pour examiner la grille pré-FFT après détendance, réflexion de symétrie ponctuelle et effilement a
appliqué, ainsi que la sauvegarde des composants réels et imaginaires du spectre brut de
les données dans topo.nc, essayez

gmt grdfft topo.nc -N+w+z -fg -V

Vous pouvez maintenant faire des tracés des données dans topo_taper.nc, topo_real.nc et topo_imag.nc.

Utilisez grdfftgmt en ligne en utilisant les services onworks.net