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
附录
其他
DCW: 世界数字图表
DCW: 世界数字图表
为了绘制某个特定行政区域的边界,首先需要知道这些行政区域的代码。
洲代码
七大洲都有各自的代码,其代码分别为:
AF: 非洲(Africa)
AN: 南极洲(Antarctica)
AS: 亚洲(Asia)
EU: 欧洲(Europe)
OC: 大洋洲(Oceania)
NA: 北美洲(North America)
SA: 南美洲(South America)
国家代码
每个国家都有一个国家代码,国家代码可以从如下三个途径查找:
在维基百科页面 ISO Country Codes 中查找
使用命令
gmt coast -E+l
查看国家代码列表从DCW辅助文件
dcw-countries.txt
中查找,其文件格式为:洲代码 国家代码 国家名
该文件共计约 250 个国家。文件内容大致如下:
AS BH Bahrain AS BN Brunei AS BT Bhutan AS CN China AS CX Christmas Island AS GE Georgia AS HK Hong Kong AS HM Heard Island and McDonald Islands AS ID Indonesia AS IL Israel AS IN India
其中可以看到,中国的国家代码为 CN。
省/州代码
目前有如下八个国家的省界/州界数据:
AR: 阿根廷
AU: 澳大利亚
BR: 巴西
CA: 加拿大
US: 美国
CN: 中国
IN: 印度
RU: 俄罗斯
省代码可以从 DCW 辅助文件 dcw-states.txt
中查找到,其文件格式为:
国家代码 省代码 省名
使用 gmt coast -E+L
命令可以列出 DCW 数据中的所有省代码。如果想进一步筛选出某个国家(以中国 CN 为例)的省代码,则可以使用 gmt coast -E+L | grep CN
命令。
以中国的数据为例,其包括全部 34 个省级行政区域:23 个省(包括台湾省),5 个自治区、4 个直辖市以及 2 个特别行政区(香港、澳门)。每个省用两个字符表示(如 AH 表示安徽):
CN.AH Anhui CN.BJ Beijing CN.CQ Chongqing CN.FJ Fujian CN.GD Guangdong CN.GS Gansu CN.GX Guangxi CN.GZ Guizhou CN.HA Henan CN.HB Hubei CN.HE Hebei CN.HI Hainan CN.HK Xianggang (Hong Kong) CN.HL Heilongjiang CN.HN Hunan CN.JL Jilin CN.JS Jiangsu CN.JX Jiangxi CN.LN Liaoning CN.MO Aomen (Macao) CN.NM Nei Mongol CN.NX Ningxia CN.QH Qinghai CN.SC Sichuan CN.SD Shandong CN.SH Shanghai CN.SN Shaanxi CN.SX Shanxi CN.TJ Tianjin CN.TW Taiwan CN.XJ Xinjiang CN.XZ Xizang CN.YN Yunnan CN.ZJ Zhejiang
使用说明
GMT中至少有两处会使用DCW数据:
洲代码、国家代码和省代码都由两个字符构成,为了避免可能的冲突,GMT通过如下方式区分:
在洲代码前加上 = 号表示某个大洲,比如 =AS 表示亚洲
国家代码不需要做任何处理格式,比如 GB 表示英国
省代码的格式为 country.state,即必须在省代码前加上国家代码才可以,比如 US.TX 表示美国 Texas 州
使用示例
绘制洲界
绘制主要大洋洲国家的边界:
gmt coast -R100/190/-50/10 -JM12c -Baf -E=OC+p0.25p,red -png dataset_dcw_01
绘制国界
绘制澳大利亚边界:
gmt coast -JM12c -Baf -EAU+p0.25p,red -png dataset_dcw_02
绘制省/洲界
绘制澳大利亚昆士兰州(Queensland),并设置边界颜色和填充颜色。其中 -R 选项后跟区域代码 AU.QLD 可间接指定该区域范围, +R2 表示在原有范围外扩大2度:
gmt coast -RAU.QLD+R2 -JM12c -Baf -EAU.QLD+p1p,blue+gred -png dataset_dcw_03
导出省/洲界数据
导出昆士兰州的边界数据:
gmt coast -EAU.QLD -M > Queensland.dat
这里只需要使用 -M 选项即可。
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论