rotsmoother
- 贡献者:
- 最近更新日期:
2025-10-17
- 官方文档:
- 简介:
从有限旋转集中获取平均旋转和协方差矩阵
rotsmoother 读取一个总重建表格，并根据旋转年龄（或所有旋转）将其划分为若干子组， 计算每个子组的平均旋转（以及可选的协方差矩阵）。
语法
gmt rotsmoother [ table ] [ -A ] [ -C ] [ -N ] [ -S ] [ -Tages ] [ -V[level] ] [ -W ] [ -Z ] [ -bbinary ] [ -dnodata[+ccol] ] [ -eregexp ] [ -fflags ] [ -hheaders ] [ -iflags ] [ -oflags ] [ -sflags ] [ -:[i|o] ] [ --PAR=value ]
注意：选项标志和相关参数之间不允许有任何空格。
必需选项
- table
包含 (lon lat time angle [weight]) 值的旋转表名称。
可选选项
- -A
使用张开角（opening angles）作为时间的代理。当无法为旋转指定时间时适用。 在这种情况下，输入应包含 lon lat angle [weight] 记录，且 -T 的设置将对应于角度而非时间。 默认情况下，输入应为 lon lat time angle [weight]，且 -T 表示时间。
- -C
为每个平均旋转计算协方差矩阵。 该过程将每个有限旋转转换为四元数（quaternion），求取平均四元数（即平均旋转）， 并将所有旋转视为相对于平均旋转的扰动。 然后根据这些扰动计算出协方差矩阵。
- -N
确保所有旋转极点位于北半球。 [默认确保张开角为正值。]
- -S
确保所有旋转极点位于南半球。 [默认确保张开角为正值。]
- -Tages
设置所需的年龄分组。 若只给定一个时间点，则附加所需的时间值，所有旋转将用于计算单个平均旋转。 若需按等间隔划分的重建年龄区间，请使用 -Tstart/stop/inc， 或 -Tstart/stop/npoints+n。 若为非等间隔的重建年龄集合，请通过文件的第一列传入，例如 -Tagefile。 我们读取或生成的年龄将作为分箱边界（bin boundaries），输出每个箱内所有旋转的平均时间。
- -V[level] (more …)
设置 verbose 等级 [w]
- -W
在输入表的最后一列中指定权重，以执行加权平均旋转 [默认不使用权重]。
- -Z
输出负的张开角 [默认输出正值]。
- -bi[ncols][type][w][+l|b] (more …)
设置二进制输入数据的格式
- -bo[ncols][type][w][+l|b] (more …)
设置二进制输出的数据格式
- -d[i|o]nodata (more …)
将输入数据中等于 nodata 的记录替换为 NaN，或将输出数据中值为 NaN 的记录替换为 nodata
- -e[~]“pattern” | -e[~]/regexp/[i] (more …)
筛选或剔除匹配指定模式的数据记录
- -f[i|o]colinfo (more …)
指定输入或输出列的数据类型
- -h[i|o][n][+c][+d][+msegheader][+rremark][+ttitle] (more …)
跳过或生成指定数目的头段记录
- -icols[+l][+sscale][+ooffset][,…][,t[word]] (more …)
设置输入数据列及简单变换（0表示第一列，t 表示文本列）
- -ocols[,…][,t[word]] (more …)
设置输出数据列（0表示第一列，t 表示文本列）
- -s[cols][+a|+r] (more …)
设置 NaN 记录的处理方式
- -:[i|o] (more …)
交换输入或输出中的第一和第二列
- -n[b|c|l|n][+a][+bBC][+c][+tthreshold] (more …)
设置网格文件的插值方式
- -^ 或 -
显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）
- -+ 或 +
显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明
- --PAR=value
临时修改GMT参数的值，可重复多次使用。参数列表见 配置参数
大地坐标与地心坐标
所有球面旋转均作用于地心坐标。 这意味着输入的数据点和网格被视为采用大地坐标表示， 在进行旋转之前，必须先将其转换为地心坐标。 旋转操作完成后，最终重建得到的点将再转换回大地坐标。 如果将椭球设置参数 {term}`PROJ_ELLIPSOID` 更改为 Sphere， 则可跳过此默认的转换过程。
示例
要以 3 Myr 的时间步长平滑旋转组，并确保极点位于北半球，执行:
gmt rotsmoother rotations.txt -N -T3/3/30 -V > rot_means.txt
要平滑所有旋转并计算单个平均旋转（分配时间为 5），同时输出对应的协方差矩阵，执行:
gmt rotsmoother rotations.txt -C -T5 -V > rot_means_cov.txt
相关模块
backtracker, gmtpmodeler, grdpmodeler, grdspotter, hotspotter, originater, rotconverter