Question

官方文档

grdinfo

简介

提取网格文件的基本信息

grdinfo 模块读取一个2D网格文件并报告网格文件的相关信息。能提取的信息包括：

• X、Y、Z的最大和最小值

• 最大/最小Z值所在的位置

• X、Y的网格间隔

• X和Y方向节点数目

• 均值、标准差

• 中位数、绝对中位差（median absolute deviation）

• the mode (Least Median of Squares; LMS), LMS scale of z

• 值为NaN的节点数

• 网格配准方式

• 网格类型（笛卡尔数据或地理数据）

语法

gmt grdinfo grdfiles [

grdfile

一个或多个网格文件名

可选选项

-C[n|t]

将输出信息以Tab分隔显示在一行中。输出格式为:

name w e s n z0 z1 dx dy nx ny [x0 y0 x1 y1] [med scale] [mean std rms] [n_nan] registration gtype

默认只输出前13列。方括号中的信息仅当使用 -M、-L1、-L2、-M 选项时才会输出。

使用 -Ct 则将文件名 name 放在最后一列；使用 -Cn 则只输出数值列。

registration 为网格的配准方式，0表示网格线配准，1表示像素配准； gtype 表示数据类型，0表示笛卡尔数据，1表示地理数据。

若与

$ gmt grdinfo @earth_relief_30m.grd
earth_relief_30m.grd: Title: Earth Relief at 30 arc minutes
earth_relief_30m.grd: Command: grdfilter SRTM15+V2.nc -Fg55.6 -D1 -I30m -rg -Gearth_relief_30m.grd=ns --IO_NC4_DEFLATION_LEVEL=9 --PROJ_ELLIPSOID=Sphere
earth_relief_30m.grd: Remark: Obtained by Gaussian Cartesian filtering (55.6 km fullwidth) from SRTM15+V2.nc [Tozer et al., 2019; https://doi.org/10.1029/2019EA000658]
earth_relief_30m.grd: Gridline node registration used [Geographic grid]
earth_relief_30m.grd: Grid file format: ns = GMT netCDF format (16-bit integer), CF-1.7
earth_relief_30m.grd: x_min: -180 x_max: 180 x_inc: 0.5 (30 min) name: longitude n_columns: 721
earth_relief_30m.grd: y_min: -90 y_max: 90 y_inc: 0.5 (30 min) name: latitude n_rows: 361
earth_relief_30m.grd: z_min: -9458 z_max: 5888 name: elevation (m)
earth_relief_30m.grd: scale_factor: 1 add_offset: 0
earth_relief_30m.grd: format: netCDF-4 chunk_size: 145,181 shuffle: on deflation_level: 9

从输出中可以看到很多信息：

• 网格文件中的标题信息；

• 生成该网格文件的命令；

• 网格文件的配准方式，此处为Gridline配准；

• 数据格式为 ns 即16位整型；

• 数据中X维度的最小值x_min、最大值x_max、网格间隔x_inc以及数据点数nx；

• 数据中Y维度的最小值y_min、最大值y_max、网格间隔y_inc以及数据点数ny；

• 数据中Z值的最小值z_min和最大值z_max以及其他信息；

相关模块

grd2cpt, grd2xyz, grdedit