grd2xyz

官方文档:grd2xyz
简介:将网格文件转换成表数据

grd2xyz 读取一个或多个二进制2D网格文件，并将XYZ数据以ASCII或二进制格式写到标准输出中。ASCII输出的格式由参数 FORMAT_FLOAT_OUT 控制，也可以以单精度或双精度浮点数的形式储存为二进制格式，还可以只输出Z值而不包含 XY坐标数据。

语法

gmt grd2xyz grid [ -C[f|i] ] [ -Rregion ] [ -V[level] ] [ -W[a|weight] ] [ -Z[flags] ] [ -bobinary ] [ -dnodata ] [ -fflags ] [ -ho[n] ] [ -oflags ] [ -sflags ] [ --PAR=value ]

必选选项

grid
要转换的2D网格文件

可选选项

-C[f|i]

输出的XY坐标值用对应的列、行号替代

默认输出的三列数据是：X坐标、Y坐标和Z值。使用该选项，则输出的三列数据为：列号、行号和Z值。其中，行号和列号从0开始算起。使用 -Cf 则行号和列号从1开始算起。若使用 -Ci 会输出两列数据：索引值和Z值。索引值相当于是将二维数组用一维数组表示

-Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)
指定数据范围
使用 -R 选项指定只对网格数据的一个子区域进行操作。若该子区域超过网格边界，则只输出二者共同的区域
-V[level] (more …)
设置 verbose 等级 [c]
-W[a|weight]

输出四列数据XYZW，其中W为 weight [weight默认值为1]

若使用 -Wa 则权重为每个节点所占据的面积。

-Z[flags]

以 ASCII 或二进制形式输出表数据

使用该选项，则输出时只有Z值，没有XY信息。输出Z值的顺序由 flags 决定。若是行优先，flags 的第一个字符可以取：

  • T 表示第一行是y=ymax
  • B 表示第一行是y=ynin

flags 的第二个字符可以取：

  • L 表示每一行的第一个元素是x=xmin
  • R 表示每一行的第一个元素是x=xmax

若是列优先，则 L|R 为第一个字符，B|T 为第二个字符。

对于网格线配准的网格文件而言：

  • 若网格在X方向是周期的，输出数据时不需要包含x=xmax所在的列，则加上 x
  • 若网格在Y方向是周期的，输出数据时不需要包含y=ymax所在的行，则加上 y

若数据需要做字节交换，则加上 w。最后需要指定数据以何种数据类型保存：

  • a ASCII表，每行输出一个Z值
  • c int8_t, signed 1-byte character
  • u uint8_t, unsigned 1-byte character
  • h int16_t, short 2-byte integer
  • H uint16_t, unsigned short 2-byte integer
  • i int32_t, 4-byte integer
  • I uint32_t, unsigned 4-byte integer
  • l int64_t, long (8-byte) integer
  • L uint64_t, unsigned long (8-byte) integer
  • f 4-byte floating point single precision
  • d 8-byte floating point double precision

默认值为 -ZTLa。

-bo[ncols][type] (more …)
设置二进制输出的数据格式
该选项只适用于XYZ输出。若只输出Z值，参考 -Z 选项指定只对网格数据的一个子区域进行操作。若该子区域超过网格边界，
-d[i|o]nodata (more …)
将输入数据中等于 nodata 的记录替换为 NaN，或将输出数据中值为 NaN 的记录替换为 nodata
-f[i|o]colinfo (more …)
指定输入或输出列的数据类型
-ocols[,…][,t[word]] (more …)
设置输出数据列（0表示第一列，t 表示文本列）
-s[cols][+a|+r] (more …)
设置 NaN 记录的处理方式
-^ 或 -
显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）
-+ 或 +
显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明
-? 或无参数
显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明
--PAR=value
临时修改GMT参数的值，可重复多次使用。参数列表见 gmt.conf

时间坐标

GMT可以识别netCDF网格文件中的时间坐标。netCDF变量的 unit 属性会被解析以确定网格文件中时间坐标的起算点和单位。这些时间坐标值会被进一步根据 TIME_UNIT 和 TIME_EPOCH 转换为GMT内部时间值。输出时，默认以相对时间的形式输出，也可以使用 -f 选项指定以绝对时间方式输出。

示例

将一个netCDF网格文件转换为XYZ文件:

gmt grd2xyz @AFR.nc > AFR.xyz

将一个netCDF文件以单精度二进制格式输出其Z值:

gmt grd2xyz @AFR.nc -ZTLf > AFR.b

相关模块

grdedit, grdconvert, xyz2grd