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 ] [ -eregexp ] [ -fflags ] [ -hheaders ] [ -iflags ] [ -qflags ] [ -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。建议用户事先使用 blockmeanblockmedian,或 blockmode 预处理数据,以避免空间假频与消除冗余数据。解微分方程时,采用自然边界条件。对于经度范围是360度的地理数据,将会采用周期边界条件。

警告

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

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

必选选项

table

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

-Goutputfile.nc

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

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

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

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

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

  • +n 表明 xincyinc 不是网格间隔,而是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