grdvolume
- 官方文档:
- 简介:
计算网格数据中某个等值线所包围的表面积和体积
grdvolume 模块读取一个2D网格文件，通过指定某条等值线确定一个Z值平面，并计算由该等值线 约束的区域网格表面积、网格表面到该平面所包围的体积，以及最大平均高度（体积/面积）。 也可以指定一系列等值线，此时该模块会分别计算每个等值线范围内的表面积和体积。
语法
gmt grdvolume
grdfile
[ -Ccval or -Clow/high/delta or -Crlow/high or -Crcval]
[ -Lbase ]
[ -Rregion ]
[ -S[unit] ]
[ -T[c|h] ]
[ -V[level] ]
[ -Zfact[/shift] ]
[ -fflags ]
[ -oflags ]
[ --PAR=value ]
必选选项
- grdfile
输入的2D网格文件名
可选选项
- -C
- -Ccval or -Clow/high/delta or -Crlow/high or -Crcval
指定等值线（Z值平面）以计算由该等值线所包含的表面积、体积以及平均高度（体积/面积）。 若不使用该选项，则返回整个网格文件的表面积、体积和平均高度。
该选项有四种不同的语法：
- -L
- -Lbase
计算体积时加上从 z=base 到等值线的体积
- -S
- -S[unit]
对于地理网格，默认会将度转换为平面距离，默认单位为米。 使用该选项指定其它长度单位，则输出的表面积单位为 unit\(^2\)， 输出的体积单位则是 z_unit * unit\(^2\)。
- -T
- -T[c|h]
找到最大平均高度所对应的等值线
-Th 找到最大平均高度（体积/面积）所对应的等值线
-Tc 找到最大曲率（高度vs等值线值）所对应的等值线
- -R
- -Rxmin/xmax/ymin/ymax[+r][+uunit]
指定数据范围。 (参数详细介绍)
- -V
- -V[level]
设置 verbose 等级 [w]。 (参数详细介绍)
- -Z
- -Zfact[/shift]
将数据减去 shift 再乘以比例因子 fact
- -f
- -f[i|o]colinfo
显式指定当前输入或输出数据中每一列的数据类型。 (参数详细介绍)
- -o
- -ocols[,…][,t[word]]
对输出的数据进行列选择以及简单的代数运算。 (参数详细介绍)
- -^ 或 -
显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）
- -+ 或 +
显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明
- --PAR=value
临时修改GMT参数的值，可重复多次使用。参数列表见 配置参数
示例
计算夏威夷岛中所有陆地部分（高于0等值线）部分的面积（\(km^2\)）、体积（\(km^3\)） 和平均高度（km）:
gmt grdvolume @earth_relief_05m -R190/210/15/25 -C0 -Sk -Z0.001
计算网格表面与等值线 \(Z=250m\) 之间的体积:
gmt grdvolume peaks.nc -Se -C250
在等值线100到300范围内，以10为间隔，计算所有等值线所约束的表面积和体积:
gmt grdvolume peaks.nc -Sk -C100/300/10 > results.d
在等值线100到300范围内，以10为间隔，搜索最大平均高度（即体积与表面积的比）所对应的等值线值:
gmt grdvolume peaks.nc -Sk -C100/300/10 -Th > results.d
计算湖内从表面到300米深度范围内水的体积:
gmt grdvolume lake.nc -Cr-300/0
参考文献
Wessel, P., 1998, An empirical method for optimal robust regional-residual separation of geophysical data, Math. Geol., 30(4), 391-408.