surface

官方文档:

surface

简介:

使用可调节张量连续曲率样条插值法对数据进行网格化

语法

gmt surface [ table ] -Goutputfile.nc -Iincrement -Rregion [ -Aaspect_ratio|m ] [ -Cconvergence_limit[%] ] [ -Jparameters ] [ -Dbreakline_file[+z[level]] ] [ -Lllower ] [ -Luupper ] [ -Mmax_radius ] [ -Nmax_iterations ] [ -Q ] [ -Ssearch_radius[m|s] ] [ -T[i|b]tension_factor ] [ -V[level] ] [ -Zover-relaxation_factor ] [ -aflags ] [ -bibinary ] [ -dinodata[+ccol] ] [ -eregexp ] [ -fflags ] [ -hheaders ] [ -iflags ] [ -qiflags ] [ -rreg ] [ -:[i|o] ] [ --PAR=value ]

描述

surface 从标准输入或ASCII格式的文件 table 中读取离散的数据点 (x,y,z)，通过解微分方程，得到网格化的数据 z(x,y) 并生成grid网格文件

\[(1 - t) \nabla ^2(z) + t \nabla (z) = 0,\]

其中 t 是0到1之间的张量因子， \(\nabla\) 是拉普拉斯算子。对于势能场数据，建议 t ~ 0.25 。对于陡峭的地形数据，则建议 t ~ 0.35。建议用户事先使用 blockmean， blockmedian，或 blockmode 预处理数据，以避免空间假频与消除冗余数据。解微分方程时，采用自然边界条件。对于经度范围是360度的地理数据，将会采用周期边界条件。

警告

在每个维度上至少必须有4个点，否则无法使用 surface 。

如果你的数据是等间隔的 (x,y,z) 数据，请不要使用 surface 生成网格文件。正确的做法是使用 xyz2grd 生成网格文件。

必选选项

table

一个或多个ASCII或二进制表数据。若不提供表数据，则会从标准输入中读取。

-Goutputfile.nc

输出的文件名。输出文件的格式参见《 网格文件格式 》。

-Ixinc[+e|n][/yinc[+e|n]]

指定X和Y方向的网格间隔

  • xinc 和 yinc 为 X 和 Y 方向的网格间隔。对于地理坐标，可以指定网格间隔单位 [默认单位为度]

  • +e 微调X和Y方向范围的最大值，使得其是网格间隔的整数倍（默认会微调网格间隔以适应给定的数据范围）

  • +n 表明 xinc 和 yinc 不是网格间隔，而是X和Y方向的节点数。此时会根据节点数、网格区域范围以及网格配准方式重新计算网格间隔。

注意：

  • 若 yinc 设置为0，则表示其与 xinc 相同

  • 若使用 -Rgrdfile 选项，则网格间隔和配准方式已经根据网格文件自动初始化，此时依然可以使用 -I 和 -r 覆盖相应的值

-Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)

指定数据范围

可选选项

-Aaspect_ratio|m

设置网格纵横比 aspect_ratio ，其中 dy = dx / aspect_ratio 。对于地理数据，可以使用 -Am 设置网格纵横比为纬度平均值的余弦。不设置本项时，默认网格纵横比为1，即各向同性网格。

-Lllower and -Luupper

对微分方程的解加上限制。 -Lllower 设置下限，其中 lower 可以是一个grid网格文件的文件名、或一个固定值。 -Luupper 设置上限，同样 upper 可以是一个grid网格文件的文件名、或一个固定值。

-T[i|b]tension_factor

设置张量因子，必须在0到1之间。 不设置时默认为0。

-V

显示多项式系数等详细信息。

示例

将ASCII数据 hawaii_5x5.xyg 网格化为5分，设置 tension_factor = 0.25 :

gmt surface hawaii_5x5.xyg -R198/208/18/25 -I5m -Ghawaii_grd.nc -T0.25