上图为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

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

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

1）H font|- header

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

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

2）D [offset] pen [-|+|=]

• 绘制一条水平线

• 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_LABEL。justification 为对齐方式，可以取 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,