grdgradientgmt - 云端在线

程序：

您的姓名

grdgradient - 计算网格的方向导数或梯度

概要

梯度 in_grd 文件 out_grd 文件 [ 阿齐姆[/阿齐姆2] ] [ [a][c][o][n] ] [
[s|p]阿齐姆/埃莱夫[/环境/扩散/镜面/闪耀] ] [ 旗 ] [ [e][t][amp] [/西格玛[/抵消]]
[ 地区 [ 斜率文件 ] [ [水平] ] [ -fg [ -n]

请注意： 选项标志和相关参数之间不允许有空格。

商品描述

梯度 可用于计算给定方向上的方向导数（-A),
或方向（-S) [和幅度 (-D)] 数据的向量梯度。

输出的第一行/最后一行/列中的估计值取决于边界条件（请参阅
-L).

所需 争论

in_grd 文件
用于计算方向导数的二维网格文件。 （参见网格文件格式
下文）。

-Gout_grd 文件
方向导数的输出网格文件的名称。 （参见网格文件格式
下文）。

不是必须的 争论

-A阿齐姆[/阿齐姆2]
方向导数的方位角方向； 阿齐姆 是 x,y 方向的角度
平面从北（+y 方向）向正顺时针方向测量，以度为单位
东（+x 方向）。 方向导数的负值，
-[dz/dx*sin(阿齐姆) + dz/dy*cos(阿齐姆)]，被发现； 否定产生正值
当 z(x,y) 的斜率是下坡时 阿齐姆 方向，正确的感觉
对图像的照明进行着色（参见 格图 和 视图）通过光源
在 x,y 平面上方，从 阿齐姆 方向。 （可选）提供两个
方位角， -A阿齐姆/阿齐姆2，在这种情况下，每个方向上的梯度是
计算并保留较大的一个； 这对于
用两个方向的线状结构照亮数据，例如， -A0/270
从北（上）和西（左）照亮。

-D[a][c][o][n]
找出数据的正（上坡）梯度方向。 改为
找到坡向（下坡方向），使用 -大。 默认情况下，方向是
从北向顺时针测量，如 阿齐姆 in -A 多于。 附加 c 使用常规的
从正 x（东）方向逆时针测量的笛卡尔角。
附加 o 报告方向 (0-180) 而不是方向 (0-360)。 附加 n
为所有角度添加 90 度（例如，给予表面局部打击）。

-E[s|p]阿齐姆/埃莱夫[/环境/扩散/镜面/闪耀]
计算适合使用的朗伯辐射率 格图 和 视图。 该
朗伯反射假设有一个理想的表面，可以反射所有光线
撞击它，从各个方向观察，表面都显得同样明亮。 阿齐姆
和 提升 是光矢量的方位角和仰角。 可选地，供应 环境
扩散 镜面 闪耀 这是控制反射率属性的参数
的表面。 默认值为： 0.55/0.6/0.4/10 留下一些值
保持不变，指定 = 作为新值。 例如 -E60/30/=/0.5 设置 阿齐姆 提升
和 扩散 为 60、30 和 0.5，保留其他反射率参数
未受影响。 附加 s 使用更简单的朗伯算法。 请注意，使用此表格
您只需提供方位角和仰角参数。 附加 p 使用
Peucker 分段线性近似（更简单但更快的算法；在这种情况下
这些因素包括原料奶的可用性以及达到必要粉末质量水平所需的工艺。 阿齐姆 和 提升 硬连线为 315 度和 45 度。 这意味着即使您
提供其他值，它们将被忽略。）

-L旗 边界条件 旗 也许 x or y or xy 表明数据是周期性的
x 或 y 或两者，或 旗 也许 g 表示地理条件（x 和 y 是
经度和纬度）。 [默认使用“自然”条件（二阶偏导数正态
到边缘为零）。]

-网][amp] [/西格玛[/抵消]]
正常化。 [默认值：无归一化。] 实际梯度 g 是偏移量并且
缩放以产生归一化梯度 gn 最大输出幅度为 amp.
If amp 没有给出，默认 amp = 1。 如果 抵消 没有给出，它被设置为
平均值 g. -N 产量 gn = amp *（g - 抵消)/最大(绝对值(g - 抵消））。 -氖
使用累积拉普拉斯分布进行归一化 gn = amp * (1.0 -
exp（开方(2) * (g - 抵消)/ 西格玛）） 在哪里 西格玛 使用 L1 范数估计
(g - 抵消）如果没有给出。 -Nt 使用累积柯西标准化
分配收益率 gn = (2 * amp / PI) * 阿坦( (g - 抵消)/ 西格玛） 在哪里 西格玛
使用 (g - 抵消）如果没有给出。

-R[单元]最小/XMAX/YMIN/最大值[R] （更多的 ...）
指定感兴趣的区域。 使用 -R 选项将选择一个小节
in_grd 文件 网格。 如果该分段超出网格边界，则只有
公共区域将被提取。

-S斜率文件
具有梯度向量标量大小的输出网格文件的名称。 需要 -D
但使 -G 可选的。

-V[水平] （更多的 ...）
选择详细级别 [c]。

-fg 地理网格（经度、纬度的尺寸）将转换为米
通过使用当前椭球参数的“平坦地球”近似值。

-n[b|c|l|n][+a][+bBC][+c][+t门槛] （更多的 ...）
为网格选择插值模式。

-^ or 只是 -
打印一条关于命令语法的短消息，然后退出（注意：在 Windows 上
只用 -).

-+ or 只是 +
打印广泛的使用（帮助）消息，包括对任何
模块特定选项（但不是 GMT 通用选项），然后退出。

-? or 没有 参数
打印完整的使用（帮助）消息，包括选项的解释，然后
退出。

- 版
打印 GMT 版本并退出。

--显示数据目录
打印 GMT 共享目录的完整路径并退出。

电网 距离 单位

如果网格没有米作为水平单位，则附加 +u单元 到输入文件
要从指定单位转换为米的名称。 如果您的网格是地理网格，请转换
通过提供距离到米 -fg 代替。

HINTS

如果你不知道什么 -N 用于制作强度文件的选项 格图 or
视图，一个好的第一次尝试是 -氖0.6.

通常 255 种色调对于可视化目的来说绰绰有余。 您可以节省 75% 的磁盘空间
通过将 =nb/a 附加到输出文件名来空格 out_grd 文件.

如果您想制作大型数据集子区域的多个照明地图，并且您
需要所有地图上的照明效果保持一致，请使用 -N 选项和
提供相同的值 西格玛 和 抵消 至 梯度 对于每张地图。 一个好的猜测是
抵消 = 0和 西格玛 被发现 数据信息 -L2 or -L1 应用于非标准化梯度 grd。

如果您只是需要 x- 要么 y-网格的导数，使用 数学.

电网 文件 FORMATS

默认情况下，GMT 在 COARDS-complaint netCDF 中将网格写出为单精度浮点数
文件格式。 但是，GMT 能够在许多其他常用网格中生成网格文件
文件格式，也有助于所谓的“打包”网格，写出浮点数
数据为 1 或 2 字节整数。 要指定精度、比例和偏移量，用户应该
添加后缀 =id[/由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。/抵消[/南]]， 在哪里 id 是网格的两个字母的标识符
类型和精度，以及 由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。 和 抵消 是可选的比例因子和偏移量
应用于所有网格值，以及 南 是用于指示缺失数据的值。 如果
两个字符 id 不提供，如 =/由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。 比一个 id=nf 假设。 什么时候
阅读格子，格式一般是自动识别的。 如果不是，相同的后缀
可以添加到输入网格文件名中。 看 转换 和节网格文件格式
GMT 技术参考和说明书以获取更多信息。

当读取包含多个网格的 netCDF 文件时，GMT 将默认读取
可以在该文件中找到的第一个二维网格。 哄GMT阅读另一个
网格文件中的多维变量，追加 ?变量名 到文件名，其中
变量名 是变量的名称。 请注意，您可能需要转义特殊含义
of ? 在 shell 程序中通过在它前面放置一个反斜杠，或者通过放置
引号或双引号之间的文件名和后缀。 这 ?变量名 也可以使用后缀
为输出网格指定一个与默认值不同的变量名：“z”。 看
转换 和 GMT Technical 的 Sections modifiers-for-CF 和 grid-file-format
参考和食谱以获取更多信息，特别是关于如何阅读 3-、
4 维或 5 维网格。

示例

使用标准化梯度创建一个文件来照亮 geoid.nc 中的数据
范围 [-0.6,0.6] 模拟北向和西向的光源：

gmt grdgradient geoid.nc -A0/270 -Ggradients.nc=nb/a -Ne0.6 -V

要在文件 topo.nc 中查找海底结构的方位角方向：

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

参考文献：

Horn, BKP，山体阴影和反射图，IEEE 会议录，卷。 69、没有。
1，1981 年 14 月，第 47-XNUMX 页。 （http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)

使用 onworks.net 服务在线使用 grdgradientgmt