收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

[Generic Mapping Tools (GMT)] Generic Mapping Tools (GMT) 使用心得和随笔(1.5):散点图之 GMT...

[复制链接]
在前一篇文章中,我们介绍了使用gmt5绘制散点图的过程。大多数代码基于gmt5结构。在新版GMT中,引入了现代模式,因此我们可以更轻松地编写代码。
在本文中,我们仍然以第一篇文章中的散点图为例,重点讨论gmt5和gmt6之间的差异。
正文:(1.5)散点图(Scatter Plot)+拟合 —— GMT6版
依然是先上效果图。
) {& J  _  H1 \& U. c9 x; u
                               
登录/注册后可看大图

0 e6 Y4 C& Z' b1 b+ F7 }; Y/ Q
代码如下:
游客,如果您要查看本帖隐藏内容请回复

- M1 P' m8 ~# b& ?* e. s+ w
在GMT6中,非常重要的改变有两个:
  • 所有绘图以gmt begin开始,gmt end为止
  • 同一张图中,部分重复指令可以省略4 N7 V$ T" R0 `9 A: R! H6 O0 k, j
此外,GMT6还引入了subplot功能,可以方便地绘制子图并调整位置大小。下面,我们开始一一分析。
  • 开始绘图2 v8 C- m1 K# R5 n' ?" J
gmt begin sample_subplot jpeg,pdf) ?- f1 U+ T% T  @5 g/ }, o
与GMT5相比,GMT6开始绘图必须以‘gmt begin’开始。其次,绘图的输出格式可以直接指定。在上述代码中,我们的输出文件名为sample_subplot,并且同时输出了两种格式jpeg和PDF格式。
游客,如果您要查看本帖隐藏内容请回复
) A5 b" f; i) j; w. Y, L. s
之后,我们设定了子图标签的字体和大小,以及颜色条。因为图中的颜色条是公用的,所以我们直接在开始子图前就进行了设置。如果需要不同图使用不同颜色条,则可以分别设置。
此外,这里还有一个关键改变。GMT6的颜色条无需写入一个特定文件(如,之前我们在GMT5篇中设置的compare.cpt)。
  • 定义子图排版布局
    / T8 l# r$ e" t! n
gmt subplot begin 1x2 -Fs$size -A+jTL+o0.2/0.2 -M0.25
6 {: B, ~1 C# \& g
绘制子图,可以通过这条命令来决定整个图片的布局。
gmt subplot begin : 开始绘制子图
1x2: 表示 一行两列。类似地,3x2则表示三行两列。
-Fs$size:表示对每张子图单独设置大小为2.8英寸 (size=2.8 )。而如果换成-Ff$size,则表示整张图的大小为2.8英寸。 注:我们上面定义过默认单位为英寸。
-A:添加自动标签,默认格式是a),b),c)...。
+jTL表示标签定位的锚点在左(L)上(T)
+o0.2/0.2表示对标签的x,y坐标分别进行0.2英寸的微调。
-M0.25:表示每张子图之间的距离为0.25英寸。
以上,我们可以发现,在GMT6中只需简单一行命令,就完成了整张图的排版布局。而在,GMT5中,我们需要不断的对-X,-Y分别调试,来达到最好的排版布局。
  • 绘制子图1
    # @) }8 Q* f7 `5 t
8 D* U8 |/ G& M3 Y0 T# K  k
游客,如果您要查看本帖隐藏内容请回复
7 @4 O) B" Z) M5 Y0 @" v
除了第一行的gmt subplot set 0之外,GMT6的命令与GMT5没有本质区别。这条命令表示以下的绘图命令都在第0号图中进行,即第一张子图中。GMT6中的子图编号,是从0开始,并定位在整幅图的左上。而编号默认是按照先行后列增加的(可以+v修改为先列后行,详见手册)。比如3x2的布局中,0、1、2三张图分别是第一行的左中右,而3、4、5则是第二行。
而在第二行命令里,我们绘制底图所用的psbasemap命令,在GMT6中,简化了拼写,去掉了'ps'这两个字母。在选项中,也有一些变化,列举如下:
  • -JX? : 这里的问号表示让GMT根据之前-R范围和我们在上面定义的子图大小自行计算。这里省去了许多麻烦,保证了每张子图的大小比例都是一样的。 不写问号,则可自行定义每张子图的大小。这在需要绘制不同大小子图时,会用到。
  • -B: 之前的文章中,我们使用的是GMT4的老写法。这里,我就改成了GMT5开始的新写法。其意义首先是,x/y和WSNE之类的设置,可以分别定义。对比见下。
  • -X -K > $out 等命令的消失: 这些命令都被gmt begin/gmt end带来的模块化结构所取代了,所以不再需要。* B( R' }: O' A
GMT4: -Ba200f100g300:"ABC":/a500f100g300:"ABC":WSneGMT5/GMT6: -Bpx200f100g300+l"ABC"  -Bpy500f100g300+l"ABC"  -BWSne #本质上,命令的易读性增加了。
6 k& z/ e  D1 T0 W) s% Y/ Q
  • 子图中画点和线
    9 S/ M$ R9 d* N% R: R2 u% u/ c
gmt plot -Sc0.04 -C4 ^4 T" ]+ I1 O: }# E6 r
除却与之前意义的awk判断句式,在GMT6画点的命令简化了许多。我们把各个版本的命令对比以下:
游客,如果您要查看本帖隐藏内容请回复
7 J4 P' N2 w. m* Z' n" j
我们发现,psxy被plot所取代,而且GMT6省略了-R、-J、-O、-K、>> $out等选项以外,还省略了颜色条的配置文件。这是因为GMT6可以自动选择之前makecpt生成的颜色条配置,无需再指定了。
类似地,画线也进行了简化。
游客,如果您要查看本帖隐藏内容请回复

7 K+ N  D  V  A& P! [
  • 子图中绘制文字, @- `; P. e  Y& T+ j5 [5 s
GMT6与GMT5的文字绘制格式是一样的,仅仅对-R、-J、-O、-K、>> $out等选项进行了省略,而pstext变成了text。
之前,我们的子图标签是用pstext进行绘制的,而如今则是自动添加(subplot begin 中的-A选项)。
  • 子图2的绘制4 n. e* f5 }- n  j7 s; T
gmt subplot set 1- k2 T: q: S8 X6 J+ f
经过了布局排版,画新的子图仅仅需要这样的简单命令。它表示在它之后的命令都用于第二子图(编号1)。
  • 颜色条的绘制
    $ p3 I( [. e% O" o& T
gmt colorbar -DJMR+v+ef -C -Bx+l"Depth" -By+l"m" -Np -B209 V  j) \# x6 r& R! Q; ^( ~
颜色条的绘制与GMT5类似,仅仅是命令名从psscale改为了colorbar。
  • 绘图结束
    / f/ j% \- Q- ]6 L0 Y& J4 N' J
gmt subplot end #结束子图gmt end  #结束所有绘图
- r4 T7 y  J  e
类似于其他编程语言,有begin,就要有end。
3 T' i, D+ h- w, }1 D' W8 R
  • 小结:* f3 M1 C% b! R1 ?% C2 c, g9 G7 M
在本文中,我们只是在第一篇文章中升级了散点图。一般来说,与gmt4和gmt5相比,gmt6消除了重复的无用命令,大大简化了图像转换和多图像设置问题。因此,我最近也升级到了gmt6。
注意:升级到gmt6后,尽管gmt5兼容代码仍在继续(gmt4部分兼容,需要在每个命令之前添加GMT),但gmt4和gmt5兼容格式不再同时出现。出现以下错误:
basemap [ERROR]: Syntax error -B option: Mixing of GMT 4 and 5 level syntax is not possible

# ^1 G- o( O/ F7 y! I: \
下期预告:地图绘制,色块和等值线的绘制。
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
fdfdfdff
活跃在2022-11-6
快速回复 返回顶部 返回列表