Question

贡献者

周茂

---

官方文档

sph2grd

简介

使用球谐系数计算网格

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

语法

gmt sph2grd [table]

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 表示计算大地水准 [ 同时使用

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

gmt sph2grd @EGM96_to_360.txt -GEGM96_to_360.nc -Rg -I1 -V

由于

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