sph2grd

贡献者:

周茂


官方文档:

sph2grd

简介:

使用球谐系数计算网格

sph2grd 模块利用球谐系数表文件计算特定的网格值。球谐系数模型的格式为 L,M,C[L,M],S[L,M]

语法

gmt sph2grd [table] -Ggrdfile -Iincrement -Rregionregion [ -D[g|n] ] [ -E] [ -F[k]filter ] [ -N[norm] ] [ -Q] [ -V[level] ] [ -bibinary ] [ -hheaders ] [ -iflags ] [ -rreg ] [ -x[-]n ] [ --PAR=value ]

必选选项

table

一个或多个球谐系数文本文件 [ 或者二进制文件,见 -bi]。文件预期 的格式为四列,分别比表示阶 L,度 M,然后是余弦和正弦系数。

-Ggrdfile

grdfile 为输出网格文件名。

-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 …)

指定数据范围

可选选项

-D[g|n]

利用位模型计算场。-Dg 表示计算引力场,-Dn 表示计算大地水准 [ 同时使用 -E 选项将会计算基于椭球的场 ]。( 目前该选项还没有实现)

-E

将参考从默认的球改变为椭球。( 目前该选项还没有实现 )

-F[k]filter

根据特定的滤波算法设定滤波参数,默认为球谐系数中的阶数 L ,使用 Fk 时, 滤波参数的单位为千米。可选的滤波算法包括: a) 余弦带通滤波:附加四个波长参数 lc/lp/hp/hclc/hc 之外的系数被截断, lp/hp 范围内的则可以通过,其余部分将 tapering 。将波长替换为 - 可以跳过设置, 例如: -F-/-/50/75 为低通滤波。 b) 高斯带通滤波:附加两个波长参数 lo/hi,滤波 amplitudes 为 0.5 。将波长替换为 - 可以跳过设置,例如: -F70/- 为高通高斯滤波。

-N[norm]

对球谐系数使用规格化。可选的规格化方法包括: m 数学正规格化化 - 表面内积和为 1 (默认)。 g 大地测量规格化 - 表面内积和为 4pi 。 s 施密特规格化 - 和地磁学中相同。

-V[level] (more …)

设置 verbose 等级 [w]

-bi[ncols][type][w][+l|b] (more …)

设置二进制输入数据的格式

-h[i|o][n][+c][+d][+msegheader][+rremark][+ttitle] (more …)

跳过或生成指定数目的头段记录

-icols[+l][+sscale][+ooffset][,][,t[word]] (more …)

设置输入数据列及简单变换(0表示第一列,t 表示文本列)

-r[g|p] (more …)

设置网格配置方式 [默认为网格线配准]

-x[[-]n] (more …)

限制多核算法中能使用的核数(需要GMT开启OpenMP支持)

-^-

显示简短的帮助信息,包括模块简介和基本语法信息(Windows下只能使用 -

-++

显示帮助信息,包括模块简介、基本语法以及模块特有选项的说明

-? 或无参数

显示完整的帮助信息,包括模块简介、基本语法以及所有选项的说明

--PAR=value

临时修改GMT参数的值,可重复多次使用。参数列表见 配置参数

示例

使用 egm96_to36.txt 球谐系数文件创建 1 x 1 度全球网格:

gmt sph2grd @EGM96_to_36.txt -GEGM96_to_36.nc -Rg -I1 -V

由于 -D 选项还未实施,该命令结果没有物理意义。

参考文献

Holmes, S. A., and Featherstone, W. E., 2002, A unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalized associated Legendre functions: J. Geodesy, v. 76, p. 279-299.

相关模块

gmt, grdfft, grdmath