3-11 GMT绘制色标

范例数据和脚本下载
主要命令：psscale
主要选项：-D
具体用法：在命令行输入：

gmt psscale --help

可以得到具体的用法

目的

1. 设置水平和竖直的色标
2. 确定色标空间位置的方式
3. 色标的标注放置位置
4. 其他

先看结果：

colorbar_ver.jpg

脚本：

#!/bin/bash
RR=R105/120/5/25
JJ=JM8c
JJP=JX8C/4.5C
PS=colorbar_ver.ps
CPT=etopo1.cpt
RANGE=-130/335/5
GRD=nanhai_topo.grd
ILLU=tmp.i

gmt gmtset FORMAT_GEO_MAP ddd:mm:ssF

gmt grdgradient $GRD -A45 -Nt -G$ILLU 
gmt psxy -$JJ -$RR -T -K >$PS
gmt grdimage $GRD -$JJ -$RR -B5 -BWSne -I$ILLU -C$CPT -K -O -P >>$PS 
#gmt psscale -$JJ -$RR -DjRB+jLB+o0.5c/0c+w10c/0.3c -C$CPT -I  -G-4000/1000 -Bx1000 -By+l"m" -K -O >>$PS
gmt psscale -Dx8.5c/0.3c+jLB+w10c/0.3c -C$CPT -I  -G-4000/1000 -Bx1000 -By+l"m" -K -O >>$PS
gmt psxy -$JJ -$RR -T -O >>$PS
gmt psconvert $PS -Tj -A -D../GMT_picture/ 
#rm  *.dat *.i gmt.* 

解释：
-I表示加光照， -G设置范围
-D选项的参数：
x表示按照坐标距离设置色标位置，x8.5c/0.3c意思是距离坐标原点（左下）x方向8.5cm，y方向0.3cm；其他表示方法：
gmt psscale -$JJ -$RR -DjRB+jLB+o0.5c/0c+w10c/0.3c -C$CPT -I -G-4000/1000 -Bx1000 -By+l"m" -K -O >>$PS
如这一行的上一句：jRB, j表示用锚点来表示，RB表示right and bottom，即右下角，即坐标放在图的右下角
+jLb表示色标是以LB，即左下作为起点，然后长宽延伸是10cm和0.3cm，由+w控制
其他可选项：+m表示文字放在另外一边；+e表示色标首尾三角形，一般图上颜色超过色标范围时需要；+h设置水平色标

gmt psscale -$JJ -$RR -DjLT+jLB+o0.5c/1c+w7C/0.3c+h+e+m -C$CPT -I -G-4000/1000 -Bx1000 -By+l"m" -K -O >>$PS
得到水平色标
这里如果吧By+l"m"改为Bx+l"m",那么色标的标识“m”就会出现在最上方。

image.png