13.20 grdedit

官方文档:grdedit
简介:修改网格文件的头段或内容

该模块具有如下功能:

  1. 从2D网格文件中读入头段信息,并使用命令行中的值替换头段信息
  2. 对全球地理网格文件还可以在东西方向旋转
  3. 可以用 x, y, z 值替换网格文件中特定结点处的值

13.20.1 必选选项

<grid>
要修改的2D网格文件

13.20.2 可选选项

-A
如有必要,则对网格间隔做微调使得其与数据的范围相兼容。仅用于处理 GMT 3.1 之前版本生成的网格文件。
-C
清除网格文件头段区中生成该网格所使用的命令历史
-D[+x<xname>][+y<yname>][+z<zname>][+s<scale>][+o<offset>][+n<invalid>][+t<title>][+r<remark>]

修改网格文件头段中的基本信息:

  • +x<xname> X变量名及其单位,格式为 varname [unit] ,比如 "distance [km]"
  • +s<scale> 读入网格数据后要乘以的因子,默认值为 1
  • +o<offset> 读入数据后并乘以因子后要加入的常数,默认值为 0
  • +n<invalid> 指定特定的数用于表示该节点处无有效值,默认为NaN
  • +t<title> 网格文件的标题
  • +r<remark> 网格文件的注释信息

其他说明:

  1. 未指定的项其值保持不变
  2. 可以给一个空值以重置某一项,比如使用 +t 而不指定标题则设置标题为空
  3. 若文本中包含空格则需要用引号括起来
  4. 对于地理数据(比如 -fg )而言, <xname><yname> 会自动设置

假设数据的范围是 300/310/10/30 ,现修改数据的范围以及标题:

gmt grdedit data.nc -R-60/-50/10/30 -D+t"Gravity Anomalies"
-E[a|h|l|r|t|v]

对网格进行处理,该选项与除 -G 外的其他选项不兼容

  • -Ea rotate the grid around 180 degrees
  • -Eh 水平翻转网格(从左到右)
  • -El 逆时针将网格旋转90度
  • -Er 顺时针将网格旋转90度
  • -Et 对网格进行转置(想象成一个二维矩阵),默认使用该选项
  • -Ev 垂直旋转网格(从上到下)

将网格数据逆时针旋转90度,并将旋转后的网格写到新网格文件中:

gmt grdedit oblique.nc -El -Goblique_rot.nc
-G<outgrid>
默认情况下,该命令会直接修改并覆盖原始网格文件,使用该选项则将修改后的网格写到新的文件中。
-J<params>
使用该选项,则会将地理相关信息以 CF-1 兼容的元数据形式(可被GDAL识别)保存到 netCDF 文件中。
-N<table>
从文件 <table> 中读入XYZ数据,并用这些XYZ数据替换网格中对应节点的值。
-R<w>/<e>/<s>/<n>
修改网格文件的范围。同时,网格间隔会做相应修改。
-S

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

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

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

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

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

13.20.3 示例