返回介绍

GMT 入门

GMT 实例

GMT 进阶

附录

其他

x2sys_cross

发布于 2023-08-10 21:46:56 字数 3403 浏览 0 评论 0 收藏 0

贡献者

周茂


官方文档

x2sys_cross

简介

计算轨迹之间的交叉点

x2sys_cross 用来确定轨迹之间所有的交叉点,并输出时间,位置,沿轨距离,朝向,速度以及交叉点不符值和所有观测的平均值。默认地,x2sys_cross 会同时计算内部和外部的交叉点。

语法

gmt x2sys_cross track(s)

tracks

可以是一个或者多个 ASCII,二进制或 COARDS NetCDF 1-D 沿轨数据文件。如果输入的沿轨数据文件很多,可以把所有的文件名放在一个轨迹列表文件中,然后使用 = 后紧跟轨迹列表文件名的形式,例如,如果轨迹列表文件名为 tracks.lis,即为 =tracks.lis 。如果 tracks.lis 中的轨迹文件名不包含后缀, gmt 将自动根据 TAG 添加后缀。gmt 将首先在当前文件夹下搜索这些轨迹文件,如果设置了 $X2SYS_HOME/TAG/TAG_paths.txt 的话,然后在该文件中包含的路径中查找。如果没有设置 $X2SYS_HOME ,默认设置为 $GMT_SHAREDIR/x2sys。

-TTAG

指定 x2sys TAG,参见 x2sys_init

可选选项

-Apairs

只处理 pairs 文件中包含的轨迹对 [默认对所有的轨迹对计算交叉点]。 pairs 可以使用 x2sys_get

交叉点不符值会以 ASCII 格式(除非设置了 -bo )输出到标准输出。输出数据的第一行为 TAG,第二行为运行的命令,第三行为每列的名称。这三个记录作为头部信息,前面都有 # 。对于每对轨迹,将另有注释信息用来记录两个轨迹文件的名称,起始/结束/距离信息,如果观测文件中有时间,起始/结束为绝对时间,如果没有,设置为 NaN;距离为整个轨迹的长度。下一行为交叉点数据记录,包括两个轨迹的 x 坐标,y 坐标,时间,距离,朝向,速度,接下来是选择的观测量。内插算法建议使用 Akima 样条内插,这种算法对离群值敏感度较小。

符号传统

如果 track_a 和 track_b 按此顺序被输入到命令中,则交叉点不符值为 Value(track_a) - Value(track_b)。

精度和格式

每列的输出格式通过 FORMAT_FLOAT_OUT 控制,如果是地理坐标,则为 FORMAT_GEO_OUT ,时间通过 FORMAT_DATE_OUTFORMAT_CLOCK_OUT 控制。请确保这些设置可以满足精度需求。

示例

计算 2104.gmt 文件的内部交叉点,TAG 设置为 GMT

gmt x2sys_cross c2104.gmt -TGMT > c2104.txt

计算测深文件 A13232.mgd77A99938.mgd77 之间的交叉点,TAG 设置为 MGD77

gmt x2sys_cross A13232.mgd77 A99938.mgd77 -Qe -TMGD77 > crossovers.txt

参考文献

Wessel, P. (2010), Tools for analyzing intersecting tracks: the x2sys package. Computers and Geosciences, 36, 348-354.

Wessel, P. (1989), XOVER: A cross-over error detector for track data, Computers and Geosciences, 15(3), 333-346.

相关模块

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

如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。

扫码二维码加入Web技术交流群

发布评论

需要 登录 才能够评论, 你可以免费 注册 一个本站的账号。
列表为空,暂无数据
    我们使用 Cookies 和其他技术来定制您的体验包括您的登录状态等。通过阅读我们的 隐私政策 了解更多相关信息。 单击 接受 或继续使用网站,即表示您同意使用 Cookies 和您的相关数据。
    原文