grdgradientgmt - En ligne dans le Cloud

PROGRAMME:

Nom

grdgradient - Calcule la dérivée directionnelle ou le gradient à partir d'une grille

SYNOPSIS

dégradé in_grdfile fichier_out_grd [ azim[/azim2] ] [ [a][c][o][n] ] [
[s|p]azim/élévation[/ambiant/diffuser/spéculaire/briller] ] [ drapeau ] [ [e][t][amp][/sigma[/compenser]]
] [ région ] [ fichier de pente ] [ [niveau] ] [ -fg ] [ -n]

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

DESCRIPTION

dégradé peut être utilisé pour calculer la dérivée directionnelle dans une direction donnée (-A),
ou la direction (-S) [et l'ampleur (-D)] du gradient vectoriel des données.

Les valeurs estimées dans la première/dernière ligne/colonne de sortie dépendent des conditions aux limites (voir
-L).

REQUIS ARGUMENTS

in_grdfile
Fichier de grille 2D à partir duquel calculer la dérivée directionnelle. (Voir FORMATS DE FICHIERS DE GRILLE
ci-dessous).

-Gfichier_out_grd
Nom du fichier de grille de sortie pour la dérivée directionnelle. (Voir FORMATS DE FICHIERS DE GRILLE
ci-dessous).

EN OPTION ARGUMENTS

-Aazim[/azim2]
Direction azimutale pour une dérivée directionnelle ; azim est l'angle dans les x,y
plan mesuré en degrés positifs dans le sens des aiguilles d'une montre du nord (la direction +y) vers
est (la direction +x). Le négatif de la dérivée directionnelle,
-[dz/dx*sin(azim) + dz/dy*cos(azim)], est trouvé; la négation donne des valeurs positives
lorsque la pente de z(x,y) est descendante dans le azim direction, le bon sens pour
ombrager l'éclairage d'une image (voir grdimage et vue grd) par une source lumineuse
au-dessus du plan x,y brillant de la azim direction. En option, fournissez deux
azimuts, -Aazim/azim2, auquel cas les gradients dans chacune de ces directions sont
calculé et le plus grand est retenu ; c'est utile pour
l'éclairage des données avec deux directions de structures alignées, par exemple, -A0/270
s'illumine du nord (en haut) et de l'ouest (à gauche).

-D[a][c][o][n]
Trouvez la direction du gradient positif (pente ascendante) des données. à la place
trouver l'aspect (la direction de la pente descendante), utiliser -Oui. Par défaut, les directions sont
mesuré dans le sens des aiguilles d'une montre à partir du nord, comme azim in -A dessus. Ajouter c utiliser le conventionnel
Angles cartésiens mesurés dans le sens inverse des aiguilles d'une montre à partir de la direction x positive (est).
Ajouter o pour signaler des orientations (0-180) plutôt que des directions (0-360). Ajouter n
pour ajouter 90 degrés à tous les angles (par exemple, pour donner des grèves locales de la surface).

-E[s|p]azim/élévation[/ambiant/diffuser/spéculaire/briller]
Calculez la radiance lambertienne appropriée à utiliser avec grdimage et vue grdL’
La réflexion lambertienne suppose une surface idéale qui reflète toute la lumière qui
le frappe et la surface apparaît également brillante dans toutes les directions d'observation. azim
et élever sont l'azimut et l'élévation du vecteur lumière. En option, fournir ambiant
diffuser spéculaire briller qui sont des paramètres qui contrôlent les propriétés de réflectance
de la superficie. Les valeurs par défaut sont : 0.55/0.6/0.4/10 Pour laisser certaines des valeurs
intact, spécifiez = comme nouvelle valeur. Par exemple -E60/30/=/0.5 définit le azim élever
et diffuser à 60, 30 et 0.5 et laisse les autres paramètres de réflectance
intact. Ajouter s d'utiliser un algorithme lambertien plus simple. Notez qu'avec ce formulaire
vous n'avez qu'à fournir les paramètres d'azimut et d'élévation. Ajouter p d'utiliser le
Approximation linéaire par morceaux de Peucker (algorithme plus simple mais plus rapide ; dans ce cas
le azim et élever sont câblés à 315 et 45 degrés. Cela signifie que même si vous
fournir d'autres valeurs, elles seront ignorées.)

-Ldrapeau Condition limite drapeau peut être x or y or xy indiquant que les données sont périodiques dans la plage de
x ou y ou les deux, ou drapeau peut être g indiquant les conditions géographiques (x et y sont
lon et lat). [La valeur par défaut utilise des conditions "naturelles" (seconde dérivée partielle normale
au bord est zéro).]

-Rapporter][amp][/sigma[/compenser]]
Normalisation. [Par défaut : aucune normalisation.] Les gradients réels g sont décalés et
mis à l'échelle pour produire des dégradés normalisés gn avec une amplitude de sortie maximale de amp.
If amp n'est pas donné, par défaut amp = 1. Si compenser n'est pas donné, il est mis à
moyenne de g. -N rendements gn = amp * (g - compenser)/max(abdos(g - compenser)). -Ne
normalise à l'aide d'une distribution de Laplace cumulative donnant gn = amp * (1.0 -
exp (sqrt(2) * (g - compenser)/ sigma)) où sigma est estimée en utilisant la norme L1 de
(g - compenser) s'il n'est pas indiqué. -NT normalise à l'aide d'un Cauchy cumulatif
rendement de distribution gn = (2* amp / PI) * atan( (g - compenser)/ sigma) où sigma
est estimée en utilisant la norme L2 de (g - compenser) s'il n'est pas indiqué.

-R[unité]x min/xmax/ymin/ymax[r] (plus ...)
Spécifiez la région d'intérêt. En utilisant le -R L'option sélectionnera une sous-section de
in_grdfile la grille. Si cette sous-section dépasse les limites de la grille, seule la
région commune sera extraite.

-Sfichier de pente
Nom du fichier de grille de sortie avec les magnitudes scalaires des vecteurs de gradient. A besoin -D
mais fait -G optionnel.

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

-n[b|c|l|n][+a][+bBC][+c][+tпорог] (plus ...)
Sélectionnez le mode d'interpolation pour les grilles.

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

ASTUCES

Si vous ne savez pas quoi -N options à utiliser pour créer un fichier d'intensité pour grdimage or
vue grd, un bon premier essai est -Ne0.6.

Habituellement, 255 nuances sont plus que suffisantes à des fins de visualisation. Vous pouvez économiser 75 % du disque
espace en ajoutant =nb/a au nom du fichier de sortie fichier_out_grd.

Si vous souhaitez créer plusieurs cartes illuminées de sous-régions d'un grand ensemble de données et que vous
besoin que les effets d'éclairage soient cohérents sur toutes les cartes, utilisez le -N option et
fournir la même valeur de sigma et compenser à dégradé pour chaque carte. Une bonne supposition est
compenser = 0 et sigma Trouvé par grdinfo -L2 or -L1 appliqué à un gradient grd non normalisé.

Si vous avez simplement besoin du x- ou y-dérivés de la grille, utiliser mathmath.

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.

EXEMPLES

Pour faire un fichier pour éclairer les données dans geoid.nc en utilisant des gradients exponés dans
la plage [-0.6,0.6] imitant les sources lumineuses dans les directions nord et ouest :

gmt grdgradient geoïde.nc -A0/270 -Ggradients.nc=nb/a -Ne0.6 -V

Pour trouver les orientations azimutales du tissu du fond marin dans le fichier topo.nc :

gmt grdgradient topo.nc -Dno -Gazimuths.nc -V

Références

Horn, BKP, Hill-Shading et la carte de réflectance, Actes de l'IEEE, Vol. 69, non.
1, janvier 1981, p. 14-47. (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)

Utilisez grdgradientgmt en ligne à l'aide des services onworks.net