-B 选项

-B 选项用于控制底图边框的绘制。

-B 选项有两套语法，分别用于设置底图边框以及每条边的属性，因而在一个命令中可能需要多次使用 -B 选项。若命令中没有出现 -B 选项，则不绘制底图边框。

边框设置

-B 选项在设置边框属性时的语法为:

-B[axe][+b][+gfill][+n][+olon/lat][+ttitle]

其中：

  • axes 控制显示底图的哪几条边，具体写法在下面会进一步介绍
  • +ttitle 指定当前底图的标题。该标题位于底图上方的中间位置。标题可以是任意字符串，如果是字符串中有空格，则必须用引号将字符串括起来。标题的文本属性由 FONT_TITLE 控制。标题与上边框之间的距离由 MAP_TITLE_OFFSET 控制。
  • +gfill 为底图内部填充背景色，见 填充
  • +n 只绘制边框，而不绘制刻度线和标注
  • +b 仅适用于3D底图，使用该子选项则会根据绘制3D底图的12条边
  • +oilon/lat 指定网格线的参考点。默认情况下，网格线是以北极点作为参考的，如果你想要以另一个点作为参考绘制倾斜的网格线，则可以使用 +o 子选项

通常情况下，只需要使用 axes 和 +ttitle 选项。

axes

axes 用于控制要绘制哪些边以及这些边是否有刻度或标注。axes 的格式为:

WSENZ[1234]wesez[1234]lrbtu

Source Code

../../_images/B-gmtplot-0.png

二维底图（上图左图）有四条边，分别用东西南北（WSEN或wsen）或左右上下（lrtb）的单词首字母表示。每条边都有四种状态：

  1. 不出现某个字母 => 不绘制该字母所对应的边
  2. 出现大写字母WSEN => 绘制某条边，该边有刻度、有标注
  3. 出现小写字母wsen => 绘制某条边，该边有刻度、无标注
  4. 出现小写字母lrtb => 绘制某条边，该边无刻度、无标注

下图给出了不同的 -B 选项绘制出来的效果图。读者可以修改如下命令中的 -B 选项来尝试不同搭配的效果:

gmt basemap -R0/4/0/4 -JX10c -BWS -pdf axes

Source Code

../../_images/B-gmtplot-1.png

3D底图有12条边（上图右图）。对于3D底图而言，Zzu 用于控制Z轴的绘制。

  • Z 表示有刻度和标注
  • z 表示有刻度无标注
  • u 表示无刻度无标注

默认只绘制一条Z轴，可以额外加上 1234 的任意组合来表示要绘制哪些Z轴。其中 1 始终表示位于左下角的Z轴，其他Z轴按逆时针顺序编号。加上 +b 子选项则绘制全部12条边。

下图展示了3D绘图中 -B 选项的不同用法。读者可以修改如下命令中的 -B 选项来实现不同搭配的效果:

gmt basemap -R0/10/0/10/0/10 -JX5c -JZ5c -Bxaf -Byaf -Bzaf -BwesnZ+t'-BwesnZ' -p130/30 -pdf map

Source Code

../../_images/B-gmtplot-2.png

轴设置

X轴、Y轴、Z轴，每条轴都有很多属性，包括刻度间隔、网格线间隔、轴标签以及标注的间隔、前缀和单位。轴属性可以用如下语法控制:

-B[p|s][x|y|z]intervals[+aangle|n|p][+l|Llabel][+slabel][+pprefix][+uunit]

以上语法也可以被拆分为两部分:

-B*[p|s][x|y|z]intervals 和 -B*[p|s][x|y|z][+aangle|n|p][+l|Llabel][+slabel][+pprefix][+uunit]

其中，

  • p|s 表示一级属性（primary）和二级属性(secondary)
  • x|y|z 表示设置哪一条轴的属性
  • interval 设置刻度、网格线、标注的间隔
  • +aangle|n|p 用于设置标注的倾斜角度，其中 angle 是相对于水平方向的旋转角度，取值范围为-90到90。+an 等效于 +a90 即垂直标注， +ap 等效于 +a0 即平行标注。对于Y轴标注而言，不支持任意角度的标注，仅支持 +ap 和 +an。
  • +llabel 用于给指定的轴加标签。默认情况下，X轴标签文字方向平行于X轴，Y轴标签文字方向平行于Y轴。对于Y轴，可以使用 +Llabel 使得Y轴标签文字方向平行于X轴
  • +slabel 与 +llabel 类似，也用于给指定的轴添加标签。当同时使用 +llabel 和 +slabel 时，前者用于指定左轴或下轴的标签，而后者用于指定右轴和上轴的标签。
  • +pprefix 在选中的轴的标注加前缀
  • +uunit 给选中的轴的标注加单位。对于地图而言，标注的单位为度，该符号是自动添加的，由 FORMAT_GEO_MAP 控制

x|y|z

x|y|z 用于指明要设置哪条边的属性，默认值为 xy，即同时设置X轴和Y轴的属性。可以指定单个轴（比如只有 x），也可以同时指定多个轴（比如 xy 和 xyz）。如果想要不同轴有不同的设置，则需要多次使用 -B 选项，每个指定不同的轴。例如:

-Bxaf -Byaf
-Bxyzaf

interval

每个轴都有三个属性，分别是标注（annotation）、刻度（frame）和网格线（grid）。下图展示了这三个名词在绘图时的具体含义。

Source Code

../../_images/B-gmtplot-3.png

interval 用于设置这三个属性的间隔，它是一个或多个 [t]stride[±phase][u] 的组合。

  • t 可以取 a（标注）、f（刻度）、g（网格线），表明要设置轴的哪个属性的间隔
  • stride 用于设置间隔，stride 为0，表示不绘制
  • ±phase 可以用于控制标注、刻度或网格线的起算点
  • u 是间隔的单位，通常只在绘制时间轴时才使用

-B 选项还有一个可以自动计算间隔的功能，-Bafg 会根据当前的区域大小等信息自动计算合适的间隔，-Bxafg -Byafg 则会对X轴和Y轴分别计算合适的间隔。

读者可以将命令:

gmt basemap -JX10c/10c -R0/10/0/10 -Ba2f1g1 -pdf test

中的 -B 选项替换成如下不同的值并查看绘图效果以理解各个参数的含义：

  • -Ba2f1g1
  • -Bxa2 -Bya1
  • -Bxafg -Byafg
  • -Ba2+1f1g1

p|s

对于每个轴来说，都有两个等级的属性可以设置，分别称为p（Primary）和s（Secondary）。

对于地理坐标而言，通常只需要使用默认的Primary属性即可，而Secondary则主要用于坐标轴为时间轴的情况下，此时 p 和 s 分别用于指定不同尺度的时间间隔。在GMT默认的情况下，p 属性的标注比较靠近坐标轴，而 s 属性的标注离坐标轴稍远。p 和 s 的用法与区别，可以参考后面给出的例子。

地理底图

地理底图与一般的坐标轴不同，其底图类型 MAP_FRAME_TYPE 使用 fancy 形式。

Source Code

../../_images/B-gmtplot-4.png

地理底图示例1

-Ba1f15mg5m -BS

下图同时使用了 p 和 s 两级属性。这里 p 属性用于显示弧度，s 属性用于显示弧分。

Source Code

../../_images/B-gmtplot-5.png

地理底图示例2

同时使用P和S两级属性 -Bpa15mf5mg5m -BwSe -Bs1f30mg15m

笛卡尔线性轴

对于一般的线性轴而言，标注的格式由 FORMAT_FLOAT_OUT 决定，其默认值为 %g，即根据数据的大小决定用一般表示还是指数表示，小数位的数目会根据 stride 自动决定。若设置 FORMAT_FLOAT_OUT 为其他值，则会严格使用其定义的格式，比如 %.2f 表示显示两位小数。

Source Code

../../_images/B-gmtplot-6.png

笛卡尔线性轴

-R0/12/0/0.95 -JX3i/0.3i -Ba4f2g1+lFrequency+u” %” -BS

笛卡尔log10轴

由于对数坐标的特殊性，stride 参数具有特殊的含义。下面说明 stride 在对数坐标下的特殊性：

  • stride 必须是1、2、3或负整数-n。
    • 1：每10的指数
    • 2：每10的指数的1、2、5倍
    • 3：每10的指数的0.1倍
    • -n：每10的n次方出现一次
  • 在 stride 后加 l，则标注会以log10的值显示，比如100会显示成2
  • 在 stride 后加 p，则标注会以10的n次方的形式显示，比如10-5

Source Code

../../_images/B-gmtplot-7.png

对数坐标轴

(上) -R1/1000/0/1 -JX3il/0.25i -Ba1f2g3(中) -R1/1000/0/1 -JX3il/0.25i -Ba1f2g3l(下) -R1/1000/0/1 -JX3il/0.25i -Ba1f2g3p

笛卡尔指数轴

正常情况下，stride 用于生成等间隔的标注或刻度，但是由于指数函数的特性，这样的标注会在坐标轴的某一端挤在一起。为了避免这个问题，可以在 stride 后加 p，则标注会按照转换后的值等间隔出现，而标注本身依然使用未转换的值。比如，若stride=1，pow=0.5（即sqrt），则在1、4、处会出现标注。

Source Code

../../_images/B-gmtplot-8.png

指数投影坐标轴

(上) -R0/100/0/0.9 -JX3ip0.5/0.25i -Ba20f10g5 (下) -R0/100/0/0.9 -JX3ip0.5/0.25i -Ba3f2g1p

时间轴

时间轴与其他轴不同的地方在于，时间轴可以有多种不同的标注方式。下面会用一系列示例来演示时间轴的灵活性。在下面的例子中，尽管只绘制了X轴（绘图时使用了 -BS），实际上时间轴标注的各种用法使用于全部轴。

在绘制时间轴时，需要指定时间间隔，时间间隔的单位可以取如下值：

GMT时间单位
Flag Unit Description
Y year Plot using all 4 digits
y year Plot using last 2 digits
O month Format annotation using FORMAT_DATE_MAP
o month Plot as 2-digit integer (1--12)
U ISO week Format annotation using FORMAT_DATE_MAP
u ISO week Plot as 2-digit integer (1--53)
r Gregorian week 7-day stride from start of week (see TIME_WEEK_START)
K ISO weekday Plot name of weekday in selected language
k weekday Plot number of day in the week (1--7) (see TIME_WEEK_START)
D date Format annotation using FORMAT_DATE_MAP
d day Plot day of month (1--31) or day of year (1--366) (FORMAT_DATE_MAP)
R day Same as d; annotations aligned with week (see TIME_WEEK_START)
H hour Format annotation using FORMAT_CLOCK_MAP
h hour Plot as 2-digit integer (0--24)
M minute Format annotation using FORMAT_CLOCK_MAP
m minute Plot as 2-digit integer (0--60)
S seconds Format annotation using FORMAT_CLOCK_MAP
s seconds Plot as 2-digit integer (0--60)

第一个例子展示了2000年春天的两个月，想要将这两个月的每周的第一天的日期标注出来。

gmt begin GMT_-B_time1 pdf,png
gmt set FORMAT_DATE_MAP=-o FONT_ANNOT_PRIMARY +9p
gmt basemap -R2000-4-1T/2000-5-25T/0/1 -JX5i/0.2i -Bpa7Rf1d -Bsa1O -BS
gmt end
../../_images/B-gmtplot-9.png

时间轴示例1

需要注意，-Bsa1O 指定了次级标注的间隔为一个月，由于此处使用的是大写的 O，因而具体的显式方式由 FORMAT_DATE_MAP 决定。根据 FORMAT_DATE_MAP 的说明可知，其值为 -o 表明以月份名格式显式。破折号表示要去掉日期前面的前置零（即02变成2）。

下面的例子用两种不同的方式标注了1969年的两天。图中下面的例子使用周来标注，上面的例子使用日期来标注。

gmt begin GMT_-B_time2 pdf,png
gmt set FORMAT_DATE_MAP "o dd" FORMAT_CLOCK_MAP hh:mm FONT_ANNOT_PRIMARY +9p
gmt basemap -R1969-7-21T/1969-7-23T/0/1 -JX5i/0.2i -Bpa6Hf1h -Bsa1K -BS
gmt basemap -Bpa6Hf1h -Bsa1D -BS -Y0.65i
gmt end
../../_images/B-gmtplot-10.png

时间轴示例2

第三个例子展示了两年的时间，并标注了每年以及每三个月。年标注位于一年间隔的中间，月标注位于对应月的中间而不是三个月间隔的中间。

gmt begin GMT_-B_time3 pdf,png
gmt set FORMAT_DATE_MAP o FORMAT_TIME_PRIMARY_MAP Character FONT_ANNOT_PRIMARY +9p
gmt basemap -R1997T/1999T/0/1 -JX5i/0.2i -Bpa3Of1o -Bsa1Y -BS
gmt end
../../_images/B-gmtplot-11.png

时间示例3

第四个例子展示了一天中的几个小时，通过在R选项中指定 t 来使用相对时间坐标。这里使用了 p 属性和 s 属性，12小时制，时间从右向左增加：

gmt begin GMT_-B_time4 pdf,png
gmt set FORMAT_CLOCK_MAP=-hham FONT_ANNOT_PRIMARY +9p TIME_UNIT d
gmt basemap -R0.2t/0.35t/0/1 -JX-5i/0.2i -Bpa15mf5m -Bsa1H -BS
gmt end
../../_images/B-gmtplot-12.png

时间轴示例4

第五个例子用两种方式展示了几周的时间：

gmt begin GMT_-B_time5 png,pdf
gmt set FORMAT_DATE_MAP u FORMAT_TIME_PRIMARY_MAP Character \
       FORMAT_TIME_SECONDARY_MAP full FONT_ANNOT_PRIMARY +9p
gmt basemap -R1969-7-21T/1969-8-9T/0/1 -JX5i/0.2i -Bpa1K -Bsa1U -BS
gmt set FORMAT_DATE_MAP o TIME_WEEK_START Sunday FORMAT_TIME_SECONDARY_MAP Chararacter
gmt basemap -Bpa3Kf1k -Bsa1r -BS -Y0.65i
gmt end
../../_images/B-gmtplot-13.png

时间轴示例5

第六个例子展示了1996年的前5个月，每个月用月份的简写以及两位年份标注：

gmt begin GMT_-B_time6 pdf,png
gmt set FORMAT_DATE_MAP "o yy" FORMAT_TIME_PRIMARY_MAP Abbreviated
gmt basemap -R1996T/1996-6T/0/1 -JX5i/0.2i -Ba1Of1d -BS
gmt end
../../_images/B-gmtplot-14.png

时间轴示例6

第七个例子：

gmt begin GMT_-B_time7 pdf,png
gmt set FORMAT_DATE_MAP jjj TIME_INTERVAL_FRACTION 0.05 FONT_ANNOT_PRIMARY +9p
gmt basemap -R2000-12-15T/2001-1-15T/0/1 -JX5i/0.2i -Bpa5Df1d -Bsa1Y -BS
gmt end
../../_images/B-gmtplot-15.png

时间轴示例7

弧度轴 \(\pi\) 的标注

如果坐标轴以弧度为单位，用户可以直接指定 \(\pi\) 的整数倍或分数倍作为标注间隔，其格式为 [+|-][s]pi[f] ，其中 s 表示标注间隔是 \(\pi\) 的 s 倍，而 f 表示标注间隔为 \(\pi\) 的 f 分之一。

示例:

gmt basemap -JX10c/5c -R-12pi/12pi/-1/1 -Bxa3pi -pdf test1
gmt basemap -JX10c/5c -R-pi/pi/-1/1 -Bxapi4 -pdf test2

自定义轴

GMT允许用户定义标注来实现不规则间隔的标注，用法是 -Bc 后接标注文件名。

标注文件中以“#”开头的行为注释行，其余为记录行，记录行的格式为:

coord   type   [label]
  • coord 是需要标注、刻度或网格线的位置
  • type 是如下几个字符的组合
    • a 或 i 前者为annotation，后者表示interval annotation
    • 在一个标注文件中，a 和 i 只能出现其中的任意一个
    • f 表示刻度，即frame tick
    • g 表示网格线，即gridline
  • label 默认的标注为 coord 的值，若指定 label，则使用 label 的值

需要注意，coord 必须按递增顺序排列。

下面的例子展示中展示了自定义标注的用法，xannots.txt 和 yannots.txt 分别是X轴和Y轴的标注文件：

cat << EOF > xannots.txt
416.0 ig Devonian
443.7 ig Silurian
488.3 ig Ordovician
542 ig Cambrian
EOF
cat << EOF > yannots.txt
0 a
1 a
2 f
2.71828 ag e
3 f
3.1415926 ag @~p@~
4 f
5 f
6 f
6.2831852 ag 2@~p@~
EOF

gmt begin GMT_-B_custom pdf,png
gmt basemap -R416/542/0/6.2831852 -JX-5i/2.5i -Bpx25f5g25+u" Ma" -Bpycyannots.txt -BWS+glightblue
gmt basemap -R416/542/0/6.2831852 -JX-5i/2.5i -Bsxcxannots.txt -BWS \
              --MAP_ANNOT_OFFSET_SECONDARY=10p --MAP_GRID_PEN_SECONDARY=2p
gmt end
rm -f [xy]annots.txt
../../_images/B-gmtplot-16.png

自定义坐标轴