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
附录
其他
使用方法
GMT 内置数据
GMT 内置数据在安装 GMT 时就已经安装,可以直接使用。
GMT 远程数据
工作原理及技术细节
经过一些数据处理,GMT官方把一些公开数据处理成不同分辨率、不同配准方式的数据,并保存在GMT专门的数据服务器上。当用户第一次使用某个分辨率的某个数据时,GMT会自动从服务器上下载该数据文件,并保存到GMT用户目录下的特定子目录下(默认为 ~/.gmt/server
目录的子目录),然后再读取该文件。当用户再次使用该数据时,GMT会自动从本地GMT用户目录下读取该数据文件,而无需再次从服务器下载。
数据分辨率
通常,公开数据提供者只提供一个原始分辨率版本。但有些时候用户不需要原始分辨率。GMT对原始数据做笛卡尔高斯滤波以进行减采样,提供了多种不同分辨率以及不同配准方式的版本,以满足不同的用户需求。
数据分块
通常,一个分辨率、一个配准方式的数据以单个文件的形式保存在GMT服务器中。但,对于高分辨率的数据,这会导致单个文件的大小变得很大,使得网络下载变得困难。为了解决这一问题,对于高分辨率的数据,GMT官方将其分割成若干个小区块,每个区块以单独的文件保存在服务器上,通常每个区块对应的文件大小约5 MB或更小。当用户需要某个区域的高分辨率数据时,可以通过 -R 选项指定区域,GMT会自动下载该区域的所有区块的数据,然后自动合并成单个网格数据供用户使用。这样可以极大减小数据下载量。例如,想要用3秒的数据绘制2度x2度的区域,GMT只需要下载该区域几MB的网格数据,而无需下载整个7 GB的3秒数据,极大的节省了下载时间。
数据压缩
GMT以无损且高压缩率的netCDF-4短整型格式保存单个网格数据,使得GMT提供的数据比原始数据占据更小的硬盘空间,而不丢失任何分辨率信息。对于分区块的高分辨率数据,GMT将数据以 JPEG2000 格式保存在GMT服务器中。该格式具有超高的压缩率,因而可以极大减小数据下载量。当将这些 JPEG2000 格式的数据下载到本地后,GMT会自动将其转换为短整型 netCDF4 格式。转换过程中使用了GDAL的功能,因而要求GMT支持GDAL。
数据更新
当原始数据发布了新版本,或GMT提供的数据被发现存在问题时,GMT会对数据服务器上的数据进行不定时的更新。GMT只提供数据的最新版本。如果你需要老版本的数据,需要自行到数据提供者的官方网站下载数据并自行处理。
GMT检查数据更新的机制是:
第一次使用任一远程数据时,GMT会从数据服务器下载一个几十 KB 的文本文件
gmt_data_server.txt
并保存到用户的~/.gmt/server
目录下当再次使用任一远程数据时,GMT会检查本地文件
gmt_data_server.txt
的日期,若其日期比当前时间差一天(该时间阈值用户可自行控制)以上,则再次从GMT服务器下载文件gmt_data_server.txt
并替换本地的老文件读取
gmt_data_server.txt
文件的内容,根据其内容确定GMT数据服务器上的远程数据是否有更新。若本地的数据文件比GMT数据服务器上的数据老,则 GMT会自动将本地已下载的数据文件删除,并重新从GMT服务器下载。由此使用远程数据的更新。
使用方式
GMT为每种远程数据都提供了不同的分辨率和配准方式的数据,具体信息在各个数据的小节会详细介绍。所有GMT远程数据都可以通过指定特殊文件名的方式来使用: @remote_name_rru_reg
@ 符号用于告诉GMT这是一个远程数据
remote_data 是远程数据的具体名字
rr 是一个两位整数,用于指定要使用的数据分辨率,单位由 u 指定。单位 u 可以取为 d、m 或 s 分别表示度、分和秒
reg 指定了数据的网格配准方式,可以取 g 或 p,分别表示网格线配准或像素配准。该项是可选的,若不指定 reg,则默认返回网格线配准的网格数据(除非不存在网格线配准的数据)。
除此之外,某些远程数据还有一个默认的CPT文件,若用户不指定CPT文件,则在绘制数据时会使用该默认CPT文件。
备注
关于每个数据的具体名称、分辨率、默认配准方式,在各个数据的小节中会详细列出。
例如,若需要使用30分的网格线配准的地形数据,则可以直接指定文件名:
gmt grdimage @earth_relief_30m_g -pdf map
对于分区块的高分辨率数据,如果想要将某个区域的数据以单个文件的形式保存,并供其它软件或命令使用,可以使用 grdcut 命令实现。例如:
gmt grdcut @earth_relief_02m_g -Gearth_at_2m.grd
参数设置
你可以通过多个参数和命令进一步控制远程数据的获取过程以及数据在本地所占据的硬盘空间:
通过参数 GMT_DATA_SERVER 设置使用离你最近的GMT数据服务器以加快下载速度
通过参数 GMT_DATA_SERVER_LIMIT 设置允许下载的远程数据的文件大小,默认无限制
使用参数 GMT_AUTO_DOWNLOAD 彻底关闭数据自动下载功能。若关闭,你将无法再下载或更新远程数据,但已下载到本地的数据依然可以使用。
通过参数 GMT_DATA_UPDATE_INTERVAL 控制检查远程数据是否更新的频率 [默认为每天检查一次]
通过 clear 命令清理整个 server 目录下的远程数据或部分数据
备注
GMT数据服务器目前在全球范围内有多个镜像。详细的镜像列表见 https://www.generic-mapping-tools.org/mirrors 页面。
对于国内用户,建议直接使用中科大LUG提供的国内镜像。修改方式为:
gmt set GMT_DATA_SERVER http://china.generic-mapping-tools.org
然后将生成的 gmt.conf
文件复制到GMT用户目录 ~/.gmt
(Linux/macOS) 或 C:\Users\XXX\.gmt
(Windows)下。
离线使用
如果你即将有一段无网络或低网速的时间,或者你不希望在第一次需要某个数据等待数据的下载,你可以使用 gmtget 模块将数据提前下载到本地。
你可以下载所有远程数据:
gmt get -Ddata
也可以只下载地球相关数据:
gmt get -Ddata=earth
或者只下载某一类或某几类数据:
gmt get -Ddata=earth_relief,earth_mask
高分辨率数据的数据量大,你也可以只下载低分辨率的数据。比如只需要低于1分分辨率的地形数据:
gmt get -Ddata=earth_relief -I1m
更多功能,见 gmtget 模块文档。
自定义数据
使用自定义数据(即建立 GMT 数据库)的方法很简单,只要将数据集中放在某个特定目录下,然后修改环境变量 GMT_DATADIR 使得其包含该目录即可。在使用数据库里的数据时,就可以直接指定数据文件名,GMT 会自动到 GMT_DATADIR 所指定的目录中寻找该数据,而无需指定数据文件的完整路径。
Linux 和 macOS 用户可以把数据库放在 ~/GMTDB
目录下,然后在 ~/.bashrc
中添加如下语句:
export GMT_DATADIR=~/GMTDB
Windows 用户可以把数据放在 D:\GMTDB (路径最好不要有空格)目录下,然后打开“我的电脑”->“属性”->“高级”->“环境变量”,添加环境变量,变量名为 GMT_DATADIR,值为 D:\GMTDB,最后可能需要重启电脑使得环境变量生效。
如果喜欢将不同的数据分类放在不同的目录下,则可以向 GMT_DATADIR 添加多个目录。多个目录之间用逗号 , 分隔。例如:
export GMT_DATADIR=~/GMTDB/data1,~/GMTDB/data2
Linux 和 macOS 用户,可以进一步简化为:
export GMT_DATADIR=~/GMTDB/
GMT_DATADIR 中的目录若以 / 结尾,GMT 则会在 ~/GMTDB
及其子目录下递归地寻找文件。注意,Windows 不支持这一功能。
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论