GMT 入门
GMT 实例
GMT 进阶
- 基础知识
- 表数据
- 网格文件
- 颜色表 CPT
- 标准选项
- 地图投影
- -JX: Linear, logarithmic, power, and time
- -JP: Polar/Cylindrical
- -JA: Lambert azimuthal equal area
- -JB: Albers conic equal area
- -JC: Cassini cylindrical
- -JCyl_stere: Cylindrical stereographic
- -JD: Equidistant conic
- -JE: Azimuthal equidistant
- -JF: Azimuthal gnomonic
- -JG: Azimuthal orthographic/General perspective
- -JH: Hammer equal area
- -JI: Sinusoidal equal area
- -JJ: Miller cylindrical
- -JK: Eckert equal area
- -JL: Lambert conic conformal
- -JM: Mercator cylindrical
- -JN: Robinson
- -JO: Oblique Mercator
- -JPoly: (American) polyconic
- -JQ: Equidistant cylindrical
- -JR: Winkel Tripel
- -JS: General stereographic
- -JT: Transverse Mercator
- -JU: Universal Transverse Mercator (UTM)
- -JV: Van der Grinten
- -JW: Mollweide
- -JY: Cylindrical equal area
- 配置参数
- 地学数据集
- 中国地理空间数据集
- 中文支持
- API
- 模块手册
- basemap
- begin
- clear
- clip
- coast
- colorbar
- contour
- coupe
- end
- figure
- filter1d
- fitcircle
- gmt
- gmt-config
- gmtbinstats
- gmtconnect
- gmtdefaults
- gmtget
- gmtinfo
- gmtlogo
- gmtmath
- gmtselect
- gmtset
- gmtsimplify
- gmtspatial
- gmtsplit
- gmtvector
- gmtwhich
- grd2cpt
- grd2xyz
- grdblend
- grdclip
- grdcontour
- grdconvert
- grdcut
- grdedit
- grdfill
- grdimage
- grdinfo
- grdlandmask
- grdmask
- grdmath
- grdpaste
- grdproject
- grdsample
- grdtrack
- grdtrend
- grdvector
- grdview
- grdvolume
- gshhg
- histogram
- image
- img2google
- img2grd
- kml2gmt
- makecpt
- mapproject
- mask
- meca
- nearneighbor
- plot
- polar
- project
- psconvert
- rose
- sac
- sample1d
- solar
- spectrum1d
- sph2grd
- sphdistance
- sphinterpolate
- sphtriangulate
- surface
- ternary
- text
- triangiulate
- velo
- wiggle
- x2sys_binlist
- x2sys_cross
- x2sys_datalist
- x2sys_get
- x2sys_init
- x2sys_list
- x2sys_merge
- x2sys_put
- x2sys_report
- x2sys_solve
- xyz2grd
- mgd77manage
- mgd77convert
- mgd77header
- mgd77info
- mgd77list
- mgd77magref
- mgd77path
- mgd77sniffer
- mgd77track
- legend
- docs
- subplot
- inset
附录
其他
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_OUT 和 FORMAT_CLOCK_OUT 控制。请确保这些设置可以满足精度需求。
示例
计算
2104.gmt
文件的内部交叉点,TAG 设置为 GMTgmt x2sys_cross c2104.gmt -TGMT > c2104.txt
计算测深文件
A13232.mgd77
和A99938.mgd77
之间的交叉点,TAG 设置为 MGD77gmt 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 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论