grdvolume¶
官方文档: | 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网格文件名
可选选项¶
- -Ccval or -Clow/high/delta or -Crlow/high or -Crcval
指定等值线(Z值平面)以计算由该等值线所包含的表面积、体积以及平均高度(体积/面积)。若不使用该选项,则返回整个网格文件的表面积、体积和平均高度。
该选项有四种不同的语法:
- -Ccval 指定单个等值线并计算等值线内的区域面积、体积和平均高度
- -Clow/high/delta 指定多条等距等值线并计算每个等值线所包含的区域面积、体积和平均高度
- -Crlow/high 计算两个Z值平面之间的体积
- -Crcval 计算 cval 到网格最小值范围内的体积
- -Lbase
- 计算体积时加上从 z=base 到等值线的体积
- -S[unit]
- 对于地理网格,默认会将弧度转换为 “Flat Earth” 下的距离,默认单位为米,使用该选项指定其它长度单位,则输出的表面积单位为 unit^2,输出的体积单位则是 unit^2 * z_unit。
- -T[c|h]
找到最大平均高度所对应的等值线
- -Th 找到最大平均高度(体积/面积)所对应的等值线
- -Tc 找到最大曲率(高度vs等值线值)所对应的等值线
- -Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)
- 指定数据范围
- -V[level] (more …)
- 设置 verbose 等级 [c]
- -Zfact[/shift]
- 将数据减去 shift 再乘以比例因子 fact
- -f[i|o]colinfo (more …)
- 指定输入或输出列的数据类型
- -ocols[,…][,t[word]] (more …)
- 设置输出数据列(0表示第一列,t 表示文本列)
- -^ 或 -
- 显示简短的帮助信息,包括模块简介和基本语法信息(Windows下只能使用 -)
- -+ 或 +
- 显示帮助信息,包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
- 显示完整的帮助信息,包括模块简介、基本语法以及所有选项的说明
- --PAR=value
- 临时修改GMT参数的值,可重复多次使用。参数列表见 gmt.conf
示例¶
计算夏威夷岛中所有陆地部分(高于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.