x2sys_solve

贡献者:

周茂

官方文档:

x2sys_solve

简介:

使用最小二乘平差计算系统误差改正值，即交叉点平差

x2sys_solve 使用给定的交叉点信息计算系统误差改正值，并应用于每条轨迹以 提高数据质量。本模块可消除 7 种不同的系统误差。

语法

gmt x2sys_solve -Ccolumn -TTAG -Ec|d|g|h|s|y|z [ COE_list.txt ] [ -V[level] ] [ -W[+u] ] [ -bibinary ] [ -dinodata[+ccol] ] [ -x[[-]n] ] [ --PAR=value ]

必选选项

COE_list.txt

由 x2sys_list 模块筛选的某观测量的交叉点文件。 如果使用了 -bi 选项，开始的两列为轨迹的 ID， 否则，就在最后两列给出轨迹名称。 如果不给定该输入文件，就从标准输入读取数据。

-TTAG

指定 x2sys TAG，参见 x2sys_init

-Ccolumn

指定要处理的列，需要和 x2sys_list 中准备的输入数据中的选项一致。

-Ec|d|g|h|s|y|z

设置用来拟合系统误差的模型。模型的类型可以从下面的函数 f(p) 中选择， 其中 p 是每条轨迹所包含的参数，参数个数为 m 个， GMT 将使用最小二乘拟合这些参数。

  • c f(p) = a （a 为常数）；记录中必须包含交叉点不符值（COE）, ID1 和 ID2

  • d f(p) = a + b * d （线性漂移，d 为沿轨距离）； 记录中必须包含：轨迹 1 距离（d1），轨迹 2 距离（d2），COE， ID1 和 ID2

  • g f(p) = a + b sin(y)^2 （1880-1930 重力改正）； 记录中必须包含：交叉点的纬度 y，COE，ID1，ID2

  • h f(p) = a + b cos(h) + c cos(2h) + d sin(h) + e sin(2h) （磁朝向改正，h 为交叉点处的轨迹方位角）；记录中必须包含： 轨迹 1 方位角（h1），轨迹 2 方位角（h2），COE，ID1，ID2

  • s f(p) = a * z （a 是单位缩放改正，z 是交叉点的观测值）； 记录中必须包含：z1，z2，ID1，ID2

  • t f(p) = a + b * (t - t0) （线性漂移， t -t0 是 沿轨累积时间，t 为观测时刻，t0 为测线起始时刻）； 记录中必须包含：t1-t0，t2-t0，COE，ID1，ID2

  • z f(p) = a + b * z （偏移加单位缩放）；记录中必须包含 z1，z2，ID1，ID2

可选选项

-W[+u]

使用该选项则意味着，每种观测值在其前面还有一列，表示其权重， 可用于加权最小二乘平差 [默认没有权重列]。 追加 +u 将输出不含权重的均值和标准差 [默认输出带有权重的统计值]

-V[level] (more …)

设置 verbose 等级 [w]

-bi[ncols][type][w][+l|b] (more …)

设置二进制输入数据的格式

-dinodata (more …)

将输入数据中值为 nodata 的列替换为 NaN

-x[[-]n] (more …)

限制多核算法中能使用的核数（需要GMT开启OpenMP支持）

-^ 或 -

显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）

-+ 或 +

显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明

-? 或无参数

显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明

--PAR=value

临时修改GMT参数的值，可重复多次使用。参数列表见 配置参数

注意事项

-E 选项中的大多数系统误差模型均包含一个常数偏移量。由于交叉点不符值是两个点 的差值，因此，这种常数偏移是相对量，而不是绝对量。为了最终结果可解，GMT 添加 了其他的约束，即所有的常数偏移量的和等于 0。 如果轨迹形成的集合中只包含同类轨迹，与其他轨迹不存在交叉点， 则形成两个轨迹集合，则每个集合的约束是独立的。GMT 会自动确定集合的数量， 并自动添加约束。如果需要设置特定的轨迹具有指定的偏移，则可以从每条轨迹 中减去本身的偏移，然后加上设定的偏移。

示例

以 MGD77 为 TAG 的所有轨道的重力异常均包含一个简单的常数偏移， 对该系统误差进行平差:

gmt x2sys_list COE_data.txt -V -TMGD77 -Cfaa -Fnc > faa_coe.txt
gmt x2sys_solve faa_coe.txt -V -TMGD77 -Cfaa -Ec > coe_table.txt

若重力异常观测数据中存在与时间相关的线性漂移:

gmt x2sys_list COE_data.txt -V -TMGD77 -Cfaa -FnTc > faa_coe.txt
gmt x2sys_solve faa_coe.txt -V -TMGD77 -Cfaa -Et > coe_table.txt

以 MGD77 为 TAG 的轨迹数据生成磁观测交叉点文件 COE_data.txt ， 估计磁朝向改正:

gmt x2sys_list COE_data.txt -V -TMGD77 -Cmag -Fnhc > mag_coe.txt
gmt x2sys_solve mag_coe.txt -V -TMGD77 -Cmag -Eh > coe_table.txt

基于测深交叉点估计单位缩放因子:

gmt x2sys_list COE_data.txt -V -TMGD77 -Cdepth -Fnz > depth_coe.txt
gmt x2sys_solve depth_coe.txt -V -TMGD77 -Cdepth -Es > coe_table.txt

相关模块

x2sys_binlist, x2sys_datalist, x2sys_get, x2sys_list, x2sys_put, x2sys_report, x2sys_init, x2sys_cross, x2sys_merge