2-16-gmtinfo-与-grdinfo

09 Jul 2018

有时候我们并不需要手动输入-R的范围信息，gmt提供了两个命令，帮助从输入数据中提取范围。
gmtinfo与frdinfo
如对R做变量赋值时，可以：

RR=$(gmt info table.txt -I1/1)

具体使用方法：

gmtinfo 提供2D数据的信息

优点：可以不再手动输入网格数据的信息

1. -R

   gmt info table.txt -I1/1

选项：

-Aa|f|s 输入数据为多段数据或多文件时的处理方式：

1. -Aa ：报告所有文件的总最值
2. -Af ：分别报告各个文件自己的最值
3. -As ：分别报告各段数据自己的最值

-D 对 -I 选项的结果做整体偏移，使得区域范围的中心与原数据的中心对齐。
比如数据的真实范围是 0/3 ，使用 -I2 选项会得到 0/4 ，在此基础上使用 -D 选项，则变成 -0.5/3.5

-EL

l

H

h<col> 返回第 <col> 列最值所在的记录行。

1. -T

   gmt info table.txt -T0.5+c2
   注意：-T与-I只能二选一

grdinfo

1. -R

   -I 后不接任何参数，会以 -I/ 的格式输出网格间隔 -Ir 或 -I- 会以 -R//<\s>/ 的格式输出网格文件的真实范围 -Ib 会输出网格区域范围的四个顶点对应的坐标 -I/ 会先获取网格的区域范围，并对该范围做微调使得其是 dx 和 dy 的整数倍，并 以 -R//<\s>/ 的格式输出

2. -T

gmt grdinfo gmt.grd -T5
-T-5250/3680/5

gmt grdinfo gmt.grd -T5+s
-T-5250/5250/5