chaoshu的个人博客分享 http://blog.sciencenet.cn/u/chaoshu

博文

科研绘图(1):大小、颜色和符号类型-图例绘制

已有 2730 次阅读 2022-10-31 13:34 |个人分类:GMT|系统分类:科研笔记

Figure 3

上图为Nature文章Figure3:Seasonal peak uplift from GPS.

Uplift and seismicity driven by groundwater depletion in central California

1、图形大小变化;

2、图形颜色变化;

3、图例同时展示;

绘制符号 — GMT 中文手册 (gmt-china.org)

大小、颜色和符号类型都变化的符号

前面介绍了如何绘制大小变化或颜色变化或符号类型变化的符号。那么如何实现颜色、大小以及符号三者都变化的符号呢?可以使用 -S 但不指定符号类型以及符号大小,同时使用 -C 选项。

关键问题在于,此时的数据数据格式是怎样的。根据 GMT 的规定,此时输入数据的格式为:

X   Y   Z   size    symbol

即第三列控制符号填充色,第四列为符号大小,最后一列为符号类型代码。

gmt begin symbols
gmt makecpt -Chot -T0/3/1
gmt plot -R0/10/0/10 -JX10c/10c -Baf -S -W1p,black -C << EOF2   3   0   0.3 c5   6   1   0.8 t8   2   2   0.5 iEOFgmt end show


参照上面方法绘图:gmt plot GNSS_Pha_Amp.txt -R%R% -J%J% -Scp -C

fig9_10_GNSS_Pha_Amp.png

接下来重要的是编写图例!

legend — GMT 中文手册 (gmt-china.org)

1)H font|- header

  • 为图例指定一个居中的标题。

    header 为标题,font 为文字属性。若字体为 - 则使用默认字体 FONT_TITLE

2)D [offsetpen [-|+|=]

  • 绘制一条水平线

    • offset 为线条左右顶端与图例边框的空白距离 [默认为0]

    • pen 为线条的画笔属性。若未指定 pen,则使用 MAP_GRID_PEN_PRIMARY。若 pen 设置为 -,则绘制一条不可见的线(供 V 记录使用)

    • 默认情况下,线条上下各留出四分之一的行间距,-|+|= 分别表示线条上方无空白、线条下方无空白和线条上下均无空白。

3)N [ncolumns or relwidth1 relwidth2 … relwidthn]

  • 修改图例中的列数 [默认为1列]

    该记录仅对 S 和 L 记录有效。该记录指定的列数会一直有效直到再次使用 N 记录。 ncolumns 用于指定若干个等宽的列,relwidth1 relwidth2 … relwidthn 用于指定每列所占的相对宽度,所有宽度的和应等于 -D 选项所设置的宽度相等。

4)S [dx1 symbol size fill pen [ dx2 text ]]

  • 在图例中绘制符号或线段

    若只有 S 而不接其它任何信息,则直接跳至下一列。若 symbol 取 f q 或 v,可以在符号后加上更多的子选项,详情见 plot 模块 -S 选项。某些符号(例如矩形、椭圆等)需要指定多个 size,应将多个 size 用逗号分隔作为 size 即可。如果只给了一个 size,则其余 size 由GMT默认值决定。

    • dx1 符号中心与左边框的距离。若为 - 则自动设置为最大的符号大小的一半。 dx1 除了可以指定距离,还可以使用 L|C|R 表示符号在当前列的对齐方式

    • symbol 指定要绘制的符号类型,见 plot 命令的 -S 选项。symbol 为 - 表示绘制线段

    • size 符号大小

    • fill 符号的填充色。使用 - 表示不填充。fill 也可以用 z=val 的形式从CPT文件中根据Z值查找颜色

    • pen 对于符号设置其轮廓属性,对于线段设置其画笔属性。使用 - 表示不绘制轮廓

    • dx2 是 text 与左边框的距离。使用 - 则自动设置为最大符号大小的1.5倍

    • text 是符号的文字说明,字体由 FONT_ANNOT_PRIMARY 控制

5)L font|- justification label

  • 在某一列增加指定的文字

    label 为显示的文本,font 为字体。若 font 为 - 则使用默认字体 FONT_LABELjustification 为对齐方式,可以取 L|C|R,分别表示左对齐、居中对齐和右对齐


    rem 边框的画笔属性

    gmt set MAP_FRAME_PEN faint,black

    rem 主刻度和次刻度的长度

    gmt set MAP_TICK_LENGTH 0.8p

    rem 刻度的画笔属性

    gmt set MAP_TICK_PEN faint,black

    rem 标注的开始位置与刻度尾端间的距离

    gmt set MAP_ANNOT_OFFSET 2p


    echo H 4p,Times-Roman Peak vertical position > legends.txt

    echo D 0.2p 0.2p >> legends.txt

    echo N 2 >> legends.txt

    echo L 3p,Times-Roman C Phase  >> legends.txt

    echo L 3p,Times-Roman C Amplitude >> legends.txt

    echo L 3p,Times-Roman C (Julian day) >> legends.txt

    echo L 3p,Times-Roman C (mm) >> legends.txt

    gmt legend legends.txt -DjBR+w1.5c/1.5c -F+pfaint,255/255/255,+g255

    gmt colorbar -Ccc.cpt -DjBR+w0.8c/0.1c+o1.15c/0.08c+v -Bxa30f30

    echo N 1 > legends.txt

    echo S 0.1p c 1p gray - 6p 0-3.0 >> legends.txt

    echo G 1.5p >> legends.txt

    echo S 0.1p c 3p gray - 6p 3.0-6.0 >> legends.txt

    echo G 3p >> legends.txt

    echo S 0.1p c 6p gray - 6p 6.0-9.0 >> legends.txt

    echo G 6p >> legends.txt

    echo S 0.1p c 9p gray - 6p 9.0-12.0 >> legends.txt

    gmt legend legends.txt -DjBR+w0.5c/0.5c+o0.2c/0.5c -F+pfaint,255/255/255,


fig9_10_GNSS_Pha_Amp.png



https://m.sciencenet.cn/blog-858128-1361657.html

上一篇:势(potential)
下一篇:科研绘图(2):地理数据预处理和颜色cpt

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-28 21:19

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部