Question

官方文档

grdedit

简介

修改网格文件的头段或内容

grdedit 模块具有如下功能：

• 从2D网格文件中读入头段信息，并使用命令行中的值替换头段信息

• 对全球地理网格文件沿着东西方向旋转

• 可以用 x y z 值替换网格文件中特定节点的值

语法

gmt grdedit grid [

grid

要修改的2D网格文件

可选选项

-A

如有必要，则对网格间隔做微调使得其与数据的范围相兼容。仅用于处理 GMT 3.1 之前版本生成的网格文件。

-C

清除网格文件头段区中生成该网格所使用的命令历史

-D[+xxname][+yyname][+zzname][+sscale][+ooffset][+ninvalid][+ttitle][+rremark]

给定网格文件头段中的基本信息：

• +xxname X变量名及其单位，格式为 varname [unit]，比如 “distance [km]”

• +sscale 读入网格数据后要乘以的因子，默认值为 1

• +ooffset 读入数据后并乘以因子后要加的常数，默认值为 0

• +ninvalid 指定值 invalid 用于表示该节点处无有效值，默认为NaN

• +ttitle 网格文件的标题

• +rremark 网格文件的注释信息

其它说明：

• 未指定的项其值保持不变

• 可以给一个空值以重置某一项，比如使用 +t 而不指定标题则设置标题为空

• 若字符串中包含空格则需要用引号括起来

• 若字符串中包含加号 +，需要使用 + 对其进行转义。或者可以使用单引号套双引号的方式，例如 ‘“title with + inside”’

• 若字符串中使用了shell变量，且变量值中包含加号，则需要使用 ${variable/+/\\+}

• 对于地理数据（比如 -fg ），xname 和 yname 会自动设置

-E[a|h|l|r|t|v]

对网格做变换。该选项与除

假定数据文件data.nc的范围为300/310/10/30。下面的命令修改了其数据范围并修改了标题:

gmt grdedit data.nc -R-60/-50/10/30 -D+t"Gravity Anomalies"

数据文件 world.nc 的范围为 0/360/-72/72，下面的命令对数据做了移动，使得数据范围为 -180/180/-72/72:

gmt grdedit world.nc -R-180/180/-72/72 -S

GMT 4.1.3 之前的网格文件不包含足够的信息表明某个网格文件是地理网格。为了添加这一信息，可以使用:

gmt grdedit junk.nc -fg

将网格文件 oblique.nc 逆时针旋转90度，并输出到新文件:

gmt grdedit oblique.nc -El -Goblique_rot.nc

为了确保文件 depths.nc 的经度始终为正值:

gmt grdedit depths.nc -L+p

相关模块

grd2xyz, grdfill, grdinfo xyz2grd