双 Y 轴图

示例贡献者:

陈箫翰

---

GMT 自身并没有绘制双 Y 轴图的功能，读者可以在相同位置绘制两个不同坐标系的图层，叠加在一起成为双 Y 轴图。示例如下：

gmt begin ex005
    # 用 gmt math 命令生成示例数据供接下来使用
    gmt math -T0/30/0.1 T SIN = sin.txt
    gmt math -T0/30/0.1 T SQR = sqr.txt

    # 首先绘制 Sin 函数图层
    # 设置坐标轴边框颜色、标注字体、刻度线、标签字体为蓝色
    gmt set MAP_FRAME_PEN blue MAP_TICK_PEN_PRIMARY blue
    gmt set FONT_ANNOT_PRIMARY blue FONT_LABEL blue
    # 设置x轴，y轴自动间隔，只绘制下边框(S)和左边框(W)
    gmt basemap -R0/30/-1.5/1.5 -JX10c/5c -Bxaf -Byaf+lSin -BWS
    gmt plot sin.txt -W1p,blue

    # 然后绘制 sqr 函数图层，直接在相同位置叠加
    # 设置坐标轴边框颜色、标注字体、刻度线、标签字体为红色
    gmt set MAP_FRAME_PEN red FONT_ANNOT_PRIMARY red
    gmt set MAP_TICK_PEN_PRIMARY red FONT_LABEL red
    # -R选项中改变y轴的范围
    # x轴设置相同，y轴标注间隔100，只绘制右边框(E)
    gmt basemap -R0/30/0/1000 -JX10c/5c -Bxaf -Bya200+lSqr -BE
    gmt plot sqr.txt -W1p,red

    rm sin.txt sqr.txt
gmt end show

双 Y 轴图