Question

官方文档

gmtsplit

简介

将表数据拆分成单独的段

gmtsplit 将表数据拆分以使每段数据都具有近乎恒定的方位角。该模块可认为是数据提取和 wiggle 模块之间的一个数据过滤模块，也可以用来将一个大的数据集分割成小段。

语法

gmt gmtsplit [ table ] [

table

输入表文件，可以为 2 列，3 列或者 5 列，数据值为 (x,y,[z[,d,h]])。默认情况下认为输入只包含 3 列，若输入数据为 5 列，使用

-Aazimuth/tolerance

计算方位角，只输出和方位角 azimuth 相差小于 tolerance 的段

-Ccourse_change

当检测到航向发生变化且超过 course_change 时，结束该航段

-Dminimum_distance

输出总长度大于 minimum_distance 的段，默认值为 0

-Fxy_filter/z_filter

对数据进行滤波，假定数据为 d 列（x 的增量）的函数。xy_filter 和 z_filter 分别为 x，y 值和 z 值的滤波的长度，单位和距离单位相同（见

使用 NCEI 提供的 JA020015 航次中 -R300/315/12/20 范围内东西方向的测线绘制磁异常的 wiggle 图，并对测线进行 100 km 的低通滤波和 500 km 的高通滤波

gmt mgd77list JA020015 -R300/315/12/20 -Flon,lat,mag,dist,azim |\
    gmt split -A90/15 -F100/-500 -D100 -S -V -fg |\
    gmt wiggle -R300/315/12/20 -Jm0.6i -Baf -B+tJA020015 -T1 \
    -W0.75p -Ggray -Z200 -pdf JA020015_wiggles

原始测量数据为 MGD77 格式，因此建议使用 mgd77list 模块提取 dist 和 azim 而不是使用 gmtsplit 计算。

将一个包含经纬度以及重力观测值的二进制双精度文件，拆分为 survey_xxxx.txt 形式的一系列文件，拆分原则为间隙大于 100 km

gmt split survey.bin -Nsurvey_%03d.txt -V -gd100k -D100 -: -fg -bi3d

相关模块

filter1d, mgd77list, wiggle