surfacegmt - 云端在线

程序：

您的姓名

表面 - 使用可调节张力连续曲率样条的网格表数据

概要

表面 [ 表 ] 输出文件 增量 地区 [ Aspect_ratio [ 收敛_限制 [
l降低 [ -鲁上 [ 最大迭代次数 ] [ ] [ 搜索半径[m|s] ] [ [i|b]张力因子
] [ [水平] ] [ 过度松弛因子 [ -a] [ -双] [ -du] [
-f] [ -h] [ -i] [ -:[i|o]]

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

商品描述

表面 从标准输入中读取随机间隔的 (x,y,z) 三元组 [或 表] 并产生
通过求解网格值 z(x,y) 的二进制网格文件：
(1 - T) * L (L (z)) + T * L (z) = 0

其中 T 是介于 0 和 1 之间的张力因子，L 表示拉普拉斯算子。 T = 0
给出等效于 SuperMISP 和 ISM 的“最小曲率”解
包。 最小曲率会导致不希望的振荡和错误的局部最大值或
最小值（参见 Smith 和 Wessel，1990），您可能希望使用 T > 0 来抑制这些
效果。 经验表明 T ~ 0.25 通常看起来对潜在场数据和 T
对于陡峭的地形数据，应该更大（T ~ 0.35）。 T = 1 给出一个调和曲面（没有
除控制数据点外，最大值或最小值都是可能的）。 建议将
用户预处理数据 块均值, 块中值或 块模式 避免空间
混叠并消除冗余数据。 您可以对下限和/或上限施加
解决方案。 这些可以以固定值、带值的网格或简单的形式输入
是最小/最大输入数据值。 自然边界条件适用于
边缘，除了我们应用周期性边界的 360 度范围的地理数据
经度方向的条件。

所需 争论

-G输出文件
输出文件名。 输出是一个二进制二维 .NC 文件。 注意最小的网格
维度必须至少为 4。

-I新奇[单元][=|+][/inc[单元][=|+]]
x_公司 [和可选 y_公司] 是网格间距。 可选地，附加一个后缀
修饰符。 地理 （度） 坐标: 追加 m 表示弧分或 s
表示弧秒。 如果其中一个单位 e, f, k, M, n or u 附加
相反，假定增量以米、英尺、公里、英里、航海为单位给出
英里或美国测量英尺，分别，并将转换为等效
该地区中纬度的经度度数（转换取决于
PROJ_ELLIPSOID）。 如果 /y_公司 给出但设置为 0 它将被重置等于 x_公司;
否则将转换为纬度。 全部 坐标：如果 = is
附加然后相应的最大值 x (东），或 y (北) 可能略有调整
完全适合给定的增量 [默认情况下，增量可以调整
稍微适合给定的域]。 最后，您可以不增加增量
指定 数 of 节点 通过附加需要 + 到提供的整数
争论; 然后根据节点数重新计算增量，然后
领域。 产生的增量值取决于您是否选择了一个
网格线注册或像素注册网格； 有关详细信息，请参阅应用程序文件格式。
注意：如果 -R文件 使用则表示网格间距已经初始化； 用
-I 覆盖这些值。

-R[单元]最小/XMAX/YMIN/最大值[R] （更多的 ...）
指定感兴趣的区域。

不是必须的 争论

表 一个或多个 ASCII（或二进制，见 -双[恩科斯][类型]) 数据表文件保存一个
数据列数。 如果没有给出表，那么我们从标准输入中读取。

-AAspect_ratio
纵横比。 如果需要，可以将网格各向异性添加到方程中。 进入
Aspect_ratio, 其中 dy = dx / Aspect_ratio 与网格尺寸相关。 [默认 =
1 假设各向同性网格。]

-C收敛_限制[％]
收敛极限。 当最大绝对值时，假设迭代已经收敛
任何网格值的变化都小于 收敛_限制. （单位与数据 z
单位）。 或者，通过附加 % 给出均方根偏差百分比的限制。
[默认缩放到数据的均方根偏差的 1e-4
最佳拟合（最小二乘法）平面。]。 这是最终收敛极限
所需的网格间距； 对于中间（较粗）网格，有效收敛
限制由网格间距乘数缩放。

-二降低 和 -鲁上
对输出解决方案施加限制。 l降低 设置下限。 降低 可以是
具有下限值、固定值的网格文件的名称， d 设置为最小值
输入值，或 u 对于不受约束的 [默认]。 u上 设置上限并且可以
是具有上限值、固定值的网格文件的名称， d 设置为
最大输入值，或 u 对于不受约束的 [默认]。 用于设置网格文件
限制可能包含 NaN。 在存在 NaN 的情况下，节点的限制被屏蔽
NaN 不受约束。

-N最大迭代次数
迭代次数。 迭代将停止时 收敛_限制 达到或
当迭代次数达到 最大迭代次数. 这是最后的迭代
限制在所需的网格间距； 对于中间（较粗）网格，有效
迭代限制由网格间距乘数缩放。 [默认值为 500。]

-Q 建议具有高度复合的最大公因数的网格尺寸。 这个
允许表面在解决方案中使用几个中间步骤，从而更快地产生
运行时间和更好的结果。 建议的尺寸 -Q 可以通过改变来实现
-R 和 -I. 你可以恢复 -R 和 -I 你想稍后使用 样本 or
切入点 在输出 表面.

-S搜索半径[米|秒]
搜索半径。 进入 搜索半径 与 x,y 数据的单位相同； 附加 m 以表示
弧分或 s 为弧秒。 这用于在开始之前初始化网格
第一次迭代； 除非网格是素数，否则不值得花时间
不能有区域阶段。 [默认值 = 0.0，不进行搜索。]

-T[我|b]张力因子
张力系数[s]。 这些必须介于 0 和 1 之间。张力可以用于
内部解（上面的方程，它抑制了虚假振荡）和
在边界条件下（它趋向于使接近于
边缘）。 对这两个值都使用零会导致最小曲率表面具有自由
边缘，即自然的双三次样条。 用 -钛张力因子 设置室内
紧张，和 -结核病张力因子 设置边界张力。 如果你不预先 i or
b，两者都将设置为相同的值。 [默认 = 0 两者都给出了最小值
曲率解。]

-V[水平] （更多的 ...）
选择详细级别 [c]。 -V3 每次迭代后都会报告收敛情况；
-V 只有在每个区域网格收敛后才会报告。

-Z过度松弛因子
过松系数。 该参数用于加速收敛； 这是
一个介于 1 和 2 之间的数字。 值 1 精确地迭代方程，并且将
始终确保稳定收敛。 较大的值高估了增量
收敛过程中会发生变化，并且会更快地找到解决方案，但可能会变成
不稳定。 如果您对此因子使用较大的值，最好监视
每次迭代 -VI 选项。 [默认 = 1.4 收敛很快，几乎
始终稳定。]

-a山坳=姓名[...] （更多的 ...）
设置非空间列关联 山坳=姓名.

-双[恩科斯][吨] （更多的 ...）
选择本机二进制输入。 [默认为 3 个输入列]。

-du没有数据 （更多的 ...）
替换等于的输入列 没有数据 与 NaN。

-f[我|o]信息 （更多的 ...）
指定输入和/或输出列的数据类型。

-h[我|o][n][+c][+d][+r备注][+r标题] （更多的 ...）
跳过或生成标题记录。 不用于二进制数据。

-i列[升][秒由于平均内核尺寸较大，西米棕榈的加工比类似作物简单。然而，西米棕榈的相对稀缺性降低了潜在的加工规模。][哦抵消][,...] （更多的 ...）
选择输入列（0 是第一列）。

-:[我|o] （更多的 ...）
在输入和/或输出上交换第一列和第二列。

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

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

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

- 版
打印 GMT 版本并退出。

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

电网 王斌金璟閣的價值觀 精准匹配

无论输入数据的精度如何，创建网格文件的 GMT 程序都会
在内部将网格保存在 4 字节浮点数组中。 这样做是为了节省内存
此外，大多数（如果不是全部）真实数据都可以使用 4 字节浮点数存储
值。 精度更高的数据（即双精度值）将丢失
GMT 在网格上运行或写出新网格后的精度。 为了限制损失
在处理数据时，您应该始终考虑在处理数据之前对数据进行标准化
处理。

示例

从 hawaii_5x5.xyg 中的 ASCII 数据中，使用 5 分钟重力块来网格 5
张力因子 = 0.25，一个 收敛_限制 = 0.1 毫克，将结果写入文件
调用 hawaii_grd.nc，并监控每次迭代，尝试：

gmt 表面 hawaii_5x5.xyg -R198/208/18/25 -I5m -Ghawaii_grd.nc -T0.25 -C0.1 -Vl

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