13.22 grdinfo¶
| 官方文档: | grdinfo |
|---|---|
| 简介: | 从2D网格文件中提取基本信息 |
能提取的信息包括：
- X、Y、Z的最大和最小值
- 最大/最小Z值所在的位置
- X、Y的网格间隔
- X和Y方向节点数目
- 均值、标准差
- 中位数、绝对中位差（median absolute deviation）
- the mode (LMS), LMS scale of z
- 值为NaN的节点数
- 网格配准方式
13.22.1 必选选项¶
grdfile一个或多个网格文件名
下面展示了一个简单的示例:
$ gmt grdinfo test.nc test.nc: Title: ETOPO5 global topography test.nc: Command: grdraster -R100/150/-30/30 1 -Gtest.nc test.nc: Remark: /opt/GMT-4.5.13/share/dbase/etopo5.i2 test.nc: Gridline node registration used [Geographic grid] test.nc: Grid file format: nf = GMT netCDF format (32-bit float), COARDS, CF-1.5 test.nc: x_min: 100 x_max: 150 x_inc: 0.0833333333333 name: longitude [degrees_east] nx: 601 test.nc: y_min: -30 y_max: 30 y_inc: 0.0833333333333 name: latitude [degrees_north] ny: 721 test.nc: z_min: -10376 z_max: 6096 name: m test.nc: scale_factor: 1 add_offset: 0 test.nc: format: classic
从输出中可以看到很多信息：
- 网格文件中的标题信息；
- 生成该网格文件的命令；
- 网格文件的配准方式，此处为Gridline配准；
- 数据格式为
nf即 netCDF 32位float； - 数据中X维度的最小值x_min、最大值x_max、网格间隔x_inc以及数据点数nx；
- 数据中Y维度的最小值y_min、最大值y_max、网格间隔y_inc以及数据点数ny；
- 数据中Z值的最小值z_min和最大值z_max以及其他信息；
13.22.2 可选选项¶
-C将输出信息以Tab分隔显示在一行中。输出格式为:
name w e s n z0 z1 dx dy nx ny [x0 y0 x1 y1] [med scale] [mean std rms] [n_nan]
默认只输出前11列，其含义看名字就可以猜到了。若使用了
-M、-L1、-L2、-M选项，则会输出对应的方括号内的项。若与
-I选项一起使用，则输出格式为:NF w e s n z0 z1
其中
NF是总网格数目。$ gmt grdinfo test1.nc -C test1.nc 100 150 -30 30 -10376 6096 0.0833333333333 0.0833333333333 601 721
-D[<xoff>[/<yoff>]][+i]将单个网格区域划分为多个特定大小的子区域范围，子区域大小由
-I<dx>/<dy>控制，<xoff>/<yoff>用于指定多个子区域之间的重叠量。+i子选项表明若该子区域内无数据则忽略该区域。$ gmt grdinfo earth_relief_01m.grd -D -I120/90 -R-240/-120/-90/0 -R-120/0/-90/0 -R0/120/-90/0 -R120/240/-90/0 -R-240/-120/0/90 -R-120/0/0/90 -R0/120/0/90 -R120/240/0/90
-F- 以每行输出一个信息的方式显示上例中的输出信息，与
-C选项不兼容。
-I[<dx>[/<dy>]|b|i|r]
-I: 以-I<xinc>/<yinc>的格式输出网格间隔-Ir: 以-R<w>/<e>/<s>/<n>格式输出网格文件的真实范围-Ii: 以-R<w>/<e>/<s>/<n>格式输出 img 补充包生成的网格文件的精确范围-Ib: 输出网格区域范围的四个顶点对应的坐标-I<dx>/<dy>会先获取网格的区域范围，并对该范围做微调使得其是<dx>和<dy>的整数倍，并以-R<w>/<e>/<s>/<n>的格式输出示例:
$ gmt grdinfo test1.nc -I -I0.0833333333333/0.0833333333333 $ gmt grdinfo test1.nc -Ir -R100/150/-30/30 $ gmt grdinfo test1.nc -Ib > Bounding box for test1.nc 100 -30 150 -30 150 30 100 30 $ gmt grdinfo test1.nc -I3/3 -R99/150/-30/30当
-I<dx>/<dy>与-C选项连用时:$ gmt grdinfo test1.nc -I3/3 -C 1 99 150 -30 30 -10376 6096
-L[0|1|2|p|a]报告Z值的其他信息，该选项可多次使用。
-L0: 扫描整个数据并报告Z值的范围，而不仅仅只是从网格的头段中读取Z值范围-L1: 输出中位数以及L1 scale （L1 scale= 1.4826*Median Absolute Deviation）-L2: 输出均值、标准差以及均方根-Lp: 输出 mode (LMS) 和 LMS scale-La: 输出以上全部信息
注意，对于像素配准的地理网格数据而言，每个节点代表的区域面积随着纬度的增加而减小，此时GMT报告的是网格文件在球面平均下的统计值。
-M- 寻找并报告Z最小和最大值所对应的坐标，以及值为NaN的网格点的数目
-R<w>/<e>/<s>/<n>从网格文件中取出一个子区域，并报告该子区域的信息。
若指定的区域范围超过了网格边界，则只提取公共区域内的为网格信息。
-T[<dz>][+a[<alpha>]][+s]以
-T<zmin>/<zmax>/<dz>的格式输出Z值范围-T<dz>：提取Z的最小最大值，并做微调使得最值是<dz>的整数倍，然后以-T<zmin>/<zmax>/<dz>的格式输出+a<alpha>使用该子选项，则会对网格文件中的值进行排序，并排除两端的部分数据。<alpha>代表要排除的数据的百分比，默认值为2，即排除最小的 1%以及最大的1%之后再输出Z值范围+s: 根据Z的绝对最大值，输出一个关于0对称的范围
$ gmt grdinfo test1.nc -T0.1 -T-10376/6096/0.1 $ gmt grdinfo test1.nc -T0.1+s -T-10376/10376/0.1