grdedit

官方文档

grdedit

简介

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

grdedit 模块具有如下功能:

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

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

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

语法

gmt grdedit grid [ -A ] [ -C ] [ -D[+xxname][+yyname][+zzname][+sscale][+ooffset][+ninvalid][+ttitle][+rremark] ] [ -E[a|h|l|r|t|v] ] [ -Goutgrid ] [ -Jparameters ] [ -L[+n|p] ] [ -Ntable ] [ -Rregion ] [ -S ] [ -T ] [ -V[level] ] [ -bibinary ] [ -dinodata[+ccol] ] [ -eregexp ] [ -fflags ] [ -hheaders ] [ -iflags ] [ -:[i|o] ] [ --PAR=value ]

必须选项

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 ),xnameyname 会自动设置

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

对网格做变换。该选项与除 -G 外的其它选项不兼容

  • -Ea 旋转180度

  • -Eh 水平翻转网格(从左到右)

  • -Ev 垂直旋转网格(从上到下)

  • -El 逆时针将网格旋转90度

  • -Er 顺时针将网格旋转90度

  • -Et 对网格进行转置(想象成一个二维矩阵),默认使用该变换

下图展示了不同变换的具体效果:

Source Code

../../_images/607d80ab2318e6694cf022d06e653c15.png
-Goutgrid

默认情况下,grdedit 模块会直接修改并覆盖原始网格文件。使用该选项则将修改后的网格写到新的文件中。

-Jprojection (more …)

设置地图投影方式

使用 -J 选项则将地理相关信息以 CF-1 兼容的元数据形式(可被GDAL识别)保存到 netCDF 文件中。

-L[+n|p]

调整地理网格文件的经度

默认情况下会调整 westeast 使得 west>=-180 或 east <= 180。 +n 则强制经度为负值,+p 则强制经度为正值。

-Ntable

从文件 table 中读入XYZ数据,并用这些XYZ数据替换网格中对应节点的值。

-Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)

指定数据范围

修改网格文件的范围。同时,网格间隔会做相应修改。

-S

将网格沿着经度范围整体偏移,使得其满足 -R 定义的新范围。仅用于全球地理网格数据。

例如,原数据范围是 0/360/-72/72,现将数据整体偏移180度使得数据范围是 -180/180/-72/72:

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

将一个网格线配准的文件变成像素配准的文件,或反之。

使用该选项后,网格线配准的数据的范围将在四个方向上扩大半个网格间隔,像素点配置的数据的范围将在四个方向上缩小半个网格间隔。

注意:修改网格配准方式会导致高频信息丢失,详情见 网格配准

-V[level] (more …)

设置 verbose 等级 [w]

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

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

-dinodata (more …)

将输入数据中值为 nodata 的列替换为 NaN

-e[~]“pattern” | -e[~]/regexp/[i] (more …)

筛选或剔除匹配指定模式的数据记录

-f[i|o]colinfo (more …)

指定输入或输出列的数据类型

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

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

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

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

-^-

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

-++

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

-? 或无参数

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

--PAR=value

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

示例

假定数据文件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