GMT 入门
GMT 实例
GMT 进阶
- 基础知识
- 表数据
- 网格文件
- 颜色表 CPT
- 标准选项
- 地图投影
- -JX: Linear, logarithmic, power, and time
- -JP: Polar/Cylindrical
- -JA: Lambert azimuthal equal area
- -JB: Albers conic equal area
- -JC: Cassini cylindrical
- -JCyl_stere: Cylindrical stereographic
- -JD: Equidistant conic
- -JE: Azimuthal equidistant
- -JF: Azimuthal gnomonic
- -JG: Azimuthal orthographic/General perspective
- -JH: Hammer equal area
- -JI: Sinusoidal equal area
- -JJ: Miller cylindrical
- -JK: Eckert equal area
- -JL: Lambert conic conformal
- -JM: Mercator cylindrical
- -JN: Robinson
- -JO: Oblique Mercator
- -JPoly: (American) polyconic
- -JQ: Equidistant cylindrical
- -JR: Winkel Tripel
- -JS: General stereographic
- -JT: Transverse Mercator
- -JU: Universal Transverse Mercator (UTM)
- -JV: Van der Grinten
- -JW: Mollweide
- -JY: Cylindrical equal area
- 配置参数
- 地学数据集
- 中国地理空间数据集
- 中文支持
- API
- 模块手册
- basemap
- begin
- clear
- clip
- coast
- colorbar
- contour
- coupe
- end
- figure
- filter1d
- fitcircle
- gmt
- gmt-config
- gmtbinstats
- gmtconnect
- gmtdefaults
- gmtget
- gmtinfo
- gmtlogo
- gmtmath
- gmtselect
- gmtset
- gmtsimplify
- gmtspatial
- gmtsplit
- gmtvector
- gmtwhich
- grd2cpt
- grd2xyz
- grdblend
- grdclip
- grdcontour
- grdconvert
- grdcut
- grdedit
- grdfill
- grdimage
- grdinfo
- grdlandmask
- grdmask
- grdmath
- grdpaste
- grdproject
- grdsample
- grdtrack
- grdtrend
- grdvector
- grdview
- grdvolume
- gshhg
- histogram
- image
- img2google
- img2grd
- kml2gmt
- makecpt
- mapproject
- mask
- meca
- nearneighbor
- plot
- polar
- project
- psconvert
- rose
- sac
- sample1d
- solar
- spectrum1d
- sph2grd
- sphdistance
- sphinterpolate
- sphtriangulate
- surface
- ternary
- text
- triangiulate
- velo
- wiggle
- x2sys_binlist
- x2sys_cross
- x2sys_datalist
- x2sys_get
- x2sys_init
- x2sys_list
- x2sys_merge
- x2sys_put
- x2sys_report
- x2sys_solve
- xyz2grd
- mgd77manage
- mgd77convert
- mgd77header
- mgd77info
- mgd77list
- mgd77magref
- mgd77path
- mgd77sniffer
- mgd77track
- legend
- docs
- subplot
- inset
附录
其他
grdconvert
- 官方文档
- 简介
将网格文件转换为其它网格文件格式
语法
gmt grdconvert ingrdfile
- ingrdfile[=id[+sscale][+ooffset][+ninvalid]]
要读入的网格文件。
若读入的网格文件不是标准的netCDF格式文件,则需要加上 =id 以指定网格文件格式(见 网格文件格式)。此外,
+sscale 对数据做比例缩放,即将数据乘以 scale
+ooffset 对数据做偏移,即将数据加上 offset
+ninvalid 数据中哪个值表示无效值
需要注意的是,在读入网格文件时,总是先缩放再偏移。
若 id=gd,则使用GDAL库检测数据格式并读入数据。实际上,当GMT遇到其无法识别的文件格式时,总是自动使用GDAL库读入数据,但可能会遇到问题,此时可以设置 id=gd强制使用GDAL库读取。
- -Goutgrdfile[=id[+sscale][+ooffset][+ninvalid]][:driver[/datatype]]]
要写入的网格文件。
若要写的网格文件格式不是标准的netCDF格式,则需要加上 =id 以指定网格文件格式(见 网格文件格式)。此外:
+sscale 对数据做比例缩放,即将数据乘以 scale
+ooffset 对数据做偏移,即将数据加上 offset
+ninvalid 数据中哪个值表示无效值
需要注意的是,在写网格文件时,总是先偏移再缩放。若想要将数据以整型保存以减小文件大小,子选项 +s 和 +o 经常会用到。此外,还可以使用 +sa+oa 让GMT自动选择合适的比例因子和偏移量以生成整型网格文件。
当 id= gd时,网格文件将使用GDAL库写入。此时可以进一步指定 driver 和 datatype。 driver 由GDAL提供(如 netCDF, GTiff 等),datatype 则可以取 u8|u16|i16|u32|i32|float32,其中 i 和 u 分别表示有符号和无符号整型。datatype 默认值为 float32。
写网格文件时,可以考虑设置 IO_NC4_DEFLATION_LEVEL 以减小生成的文件大小,并进一步优化读写性能。
可选选项
- -N
在生成native二进制文件时,不将GMT网格文件头段写到文件中。
- -Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)
指定数据范围
- -V[level] (more …)
设置 verbose 等级 [w]
- -Z[+sfactor][+ooffset]
在写网格文件前,从数据中减去 offset 并将结果乘以 factor。
注意:该选项功能与 +s 和 +o 子选项的功能不同。使用 +s 和 +o 子选项,数据也会进行缩放和偏移,修改后的数据会被写入到网格文件中,但同时缩放因子和偏移量会写入到网格文件的元信息中。因而,在读取这些新生成的网格时,会自动根据元信息中的值恢复数据的原始值。而使用该选项则永久修改了网格文件中的数据值。
- -f[i|o]colinfo (more …)
指定输入或输出列的数据类型
- -^ 或 -
显示简短的帮助信息,包括模块简介和基本语法信息(Windows下只能使用 -)
- -+ 或 +
显示帮助信息,包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
显示完整的帮助信息,包括模块简介、基本语法以及所有选项的说明
- --PAR=value
临时修改GMT参数的值,可重复多次使用。参数列表见 配置参数
注意事项
GMT默认只能读取并处理2D单变量网格。对于多变量、多维度网格文件,需要使用额外的语法指定要读取的变量或维度,详情见 读 netCDF 文件。
示例
将网格文件转换成四字节native浮点型网格:
gmt grdconvert data.nc ras_data.b4=bf -V
将网格文件转换成二字节短整型文件,将其乘以10并减去32000,并设置无数据节点的值为-9999:
gmt grdconvert values.nc shorts.i2=bs/10/-32000/-9999 -V
从一个三维网格文件中提取第二层数据:
gmt grdconvert climate.nc?temp[1] temp.nc -V
相关模块
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论