sample1d

官方文档:

sample1d

简介:

对1D表数据进行重采样

该命令既可以对常规的一维数据(比如时间序列,自变量为时间)进行重采样,也可以对地理坐标下的测线进行重采样。

可选选项

<table>

多列表数据,其中某一列是自变量,其他列为因变量。自变量所在列必须递增或递减。本页面将自变量称为“时间”,因为该命令常用于处理时间序列,实际上自变量可以是任意物理量

-Af|p|m|r|R[+l]

指定插值方式。

  • -Af 保留原始数据点,若有必要,则在原始数据点的中间加上额外的点

  • -Am 对测线进行采样时,先沿着Y方向,再沿着X方向

  • -Ap 对测线进行采样时,先沿着X方向,再沿着Y方向

  • -Ar 等间距采样

  • -AR 等间距采样,但会调整间距以适应自变量的原始长度

  • +l if distances should be measured along rhumb lines (loxodromes)

-Fl|a|c|n[+1|+2]

插值方式

  • l 线性插值

  • a Akima样条插值

  • c natural cubic spline

  • n 不插值,取最近的数据点作为插值后的值

  • +1|+2 插值的同时计算spline的一阶或二阶插值

-I<inc>[<unit>]

默认的等间隔采样间隔是自变量第一个和第二个数据点的间隔,该选项可以自定义采样间隔 <inc>

加上 <unit> 表明数据文件的前两列包含经纬度信息,重采样后的测线的采样间隔的单位是 <units> 。若想要采样笛卡尔坐标下的(x,y),则需要指定单位为 c

-N<knotfile>

<knotfile> 中包含了一系列X坐标轴,使用该选项则会将原始数据插值到这些 X坐标轴数据点上。

-S<start>[/<stop>]

对于等间隔采样而言, <start> 是第一个输出值的X位置, <stop> 是最后一个输出值的X位置。

-T<col>

指定输入数据中的哪列数据是自变量。

示例

输入数据的格式为:

time distance gravity magnetics bathymetry

使用Akima spline插值方式将其采样为1千米等间隔:

gmt sample1d profiles.tdgmb -I1 -Fa -T1 > profiles_equi_d.tdgmb

将0到6之间的数据用cubic spline方式重采样为0.01间隔,不输出数据而是输出一阶偏导(即斜率):

gmt sample1d points.txt -S0/6 -I0.01 -Fc+1 > slopes.txt

测线数据中包含经度、纬度和深度,将其采样为每2海里一个点:

gmt sample1d track.txt -I2n -AR > new_track.dt

同上,但确保包含了原始数据点:

gmt sample1d track.txt -I2n -Af > new_track.dt

To obtain a rhumb line (loxodrome) sampled every 5 km instead:

gmt sample1d track.txt -I5k -AR+l > new_track.dt