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方向的网格间隔
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 …)
指定数据范围
可选选项
- -D[g|n]
利用位模型计算场。-Dg 表示计算引力场，-Dn 表示计算大地水准 [ 同时使用 -E 选项将会计算基于椭球的场 ]。( 目前该选项还没有实现)
- -E
将参考从默认的球改变为椭球。（ 目前该选项还没有实现 ）
- -F[k]filter
根据特定的滤波算法设定滤波参数，默认为球谐系数中的阶数 L ，使用 Fk 时，滤波参数的单位为千米。可选的滤波算法包括：a) 余弦带通滤波：附加四个波长参数 lc/lp/hp/hc 。lc/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_to360.txt 球谐系数文件创建 1 x 1 度全球网格:
gmt sph2grd @EGM96_to_360.txt -GEGM96_to_360.nc -Rg -I1 -V
由于 -D 选项还未实施，该命令结果没有物理意义。EGM96_to_360.txt 远程文件已经失效，可从 https://cddis.nasa.gov/926/egm96/getit.html 下载。
参考文献
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.