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 中使用 coast 模块绘制海岸线。 coast 模块,顾名思义是海岸线、湖岸线的意思。该模块不仅可以绘制这些岸线,还可以绘制比例尺和行政边界等。
这一节我们将通过绘制一张美国地图来介绍 coast 模块的基本用法。
绘制海岸线和湖岸线
下面的脚本用最简单的命令绘制了一张美国地图。其中 -R 选项设置了绘图区域,-JM15c 设置投影方式和图片尺寸, -Baf 设置底图边框属性,这些在前一节都已经做了介绍。
coast 模块的选项 -W0.5p,black 表示绘制岸线,并设置岸线为 0.5p 宽的黑色线条。
gmt begin coastline gmt coast -R-130/-50/20/60 -JM15c -Baf -W0.5p,black gmt end show
美国海岸线
上图看上去依然有些乱。主要是因为美国有大大小小的很多岛屿和湖泊,GMT 默认会把这些岛屿和湖泊都画出来。如果能够只绘制比较大的岛屿或湖泊就好了,我们可以通过 -A 选项设置要绘制的对象的最小面积(单位为平方千米)。
gmt begin coastline gmt coast -R-130/-50/20/60 -JM15c -Baf -W1/0.5p,black -W2/1p,lightred -A5000 gmt end show
设置要绘制的对象的最小面积
上面介绍了如何使用 -W 绘制海岸线和湖岸线。我们还可以不绘制岸线,而是为陆地和水体设置不同的填充色。
gmt begin coastline gmt coast -R-130/-50/20/60 -JM15c -Baf -A5000 -Gred -Slightblue -Clightred gmt end show
填充颜色
使用 -N 选项可以绘制国界、州界/省界等行政边界。-N1 表示绘制国界线, -N2 表示绘制州界/省界线。该选项绘制的行政边界来自 GSHHG: 全球高分辨率海岸线数据,目前只有美洲各国以及澳大利亚有州界/省界线。
gmt begin coastline gmt coast -R-130/-50/20/60 -JM15c -Baf -A5000 -Gred -Slightblue -Clightred -N1/1p -N2/0.25p gmt end show
绘制国界与州界
最后,我们还需要为地图添加比例尺。为了绘制比例尺,我们需要提供如下参数:
要绘制哪个纬度的比例尺
比例尺在图中的位置
比例尺的长度
在下面的例子中,我们使用了 -Lg-60/25+c25+w1000k+f+u 增加比例尺,其中:
+w1000k 表示比例尺长度为 1000 千米
+c25 表示绘制纬度为北纬 25° 处的比例尺
g-60/25 则表示将比例尺画在北纬 25° 西经 60° 处
+f 表示比例尺的风格为图中所示黑白相间的铁轨形式
+u 表示显示比例尺对应的单位
gmt begin coastline gmt coast -R-130/-50/20/60 -JM15c -Baf -A5000 -Gred -Slightblue -Clightred -Lg-60/25+c25+w1000k+f+u gmt end show
添加比例尺
如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。
绑定邮箱获取回复消息
由于您还没有绑定你的真实邮箱,如果其他用户或者作者回复了您的评论,将不能在第一时间通知您!
发布评论