GMT绘图笔记

（1）GMT6绘制海洋深度分布，陆地白色地形

gmt begin csdn png,pdf E600

rem 绘制底图

rem 边框参数

gmt set MAP_FRAME_TYPE fancy

rem 标注参数

gmt set MAP_FRAME_PEN 0.8p,black

gmt set FONT_ANNOT_PRIMARY 8p,Helvetica,black

gmt set FONT_ANNOT_SECONDARY 8p,Times-Roman,black

gmt set MAP_TICK_PEN 0.8p,black

gmt set MAP_TICK_LENGTH 1.5p

gmt set LABEL_FONT 8p

gmt set FONT_LABEL 10p

gmt set LABEL_FONT_SIZE 8p

gmt set ANNOT_FONT_SIZE_PRIMARY 8p

gmt set MAP_FRAME_WIDTH 1.5p

gmt set FRAM_WIDTH 0.5c

gmt set MAP_ANNOT_OFFSET_PRIMARY 0.15c

gmt set MAP_TICK_LENGTH_PRIMARY 0.15c

gmt set FORMAT_ANNOT_PRIMARY 0.3p,Times-Roman,black

gmt set FORMAT_GEO_MAP ddd:mm:ssf

rem 绘制地形图

    gmt basemap -JX8c/5c -R110/140/25/40 -BWSen -Ba5f5 -Xc -Yc -t10    

    xyz2grd garce_djf.txt -JX8c/5c -R110/150/-25/-5 -I1 -Gm1.grd

    gmt grdsample m1.grd -Gm11.grd -I0.25

    gmt grdimage m11.grd -Chaxby

    gmt grdimage earth_relief_04m.grd -Cgray2 -I+d

    gmt grdimage earth_relief_04m.grd -Ctest1.cpt -I+d

    gmt coast -JX8c/5c -R110/140/25/40 -A500 -W0.1p,black

gmt end

结果如下，具体的cpt文件可以私信本人获取。

（2）GMT6绘制误差棒

gmt begin csdn_err png,pdf E600

    gmt set MAP_FRAME_PEN 1.2p,black

    gmt set FONT_ANNOT_PRIMARY 16p,Helvetica,black

    gmt set FONT_ANNOT_SECONDARY 16p,Times-Roman,black

    gmt set MAP_TICK_PEN 0.4p,black

    gmt set MAP_TICK_LENGTH 1p

    gmt set LABEL_FONT 7p

    gmt set FONT_LABEL 16p,Helvetica,black

    gmt set LABEL_FONT_SIZE 4p

    gmt set ANNOT_FONT_SIZE_PRIMARY 4p

    gmt set MAP_FRAME_WIDTH 3p

    gmt set MAP_FRAME_TYPE plain

    gmt set FRAM_WIDTH 2c

    gmt set MAP_ANNOT_OFFSET_PRIMARY 0.15c

    gmt set MAP_TICK_LENGTH_PRIMARY 0.15c

    gmt set FORMAT_ANNOT_PRIMARY 0.5p,Times-Roman,black

    gmt set FORMAT_GEO_MAP ddd:mm:ssf

    gmt basemap -R0.5/12.5/-6/6 -JX15c/8c -Bxf1a1 -Byf2a2 -BWSne -Byaf+l"Sea level anomaly (cm)" -Bxaf+l"Month"  -Xc -Yc

    gmt plot -R0.5/12.5/-6/6 -JX15c/8c -Sc0.2c -Gred -W0.1p,red season1.txt -Ey+w6p+p0.6,red

    gmt plot -R0.5/12.5/-6/6 -JX15c/8c -W0.5p,red season1.txt

    gmt plot -R0.5/12.5/-6/6 -JX15c/8c -Sc0.2c -Gblue -W0.1p,blue season2.txt -Ey+w6p+p0.6,blue

    gmt plot -R0.5/12.5/-6/6 -JX15c/8c -W0.5p,blue season2.txt

    echo 1 -5 (d) | gmt text -R0.5/12.5/-6/6 -JX15c/8c -F+f24p,4,black

gmt end

结果图

代码中的season.txt文件格式，分别是x、y、误差

1    0.892910462431733    2.18657827720253
2    0.816539654771009    2.18565002057906
3    1.55476239053796    2.02515516262056
4    3.01085605197287    1.72926006293462
5    2.37523950477438    1.63946766895842
6    0.281308538521274    1.57955715541389
7    -1.99860585858323    1.66805777837017
8    -3.43648315123605    1.65395955833412
9    -3.13051359594750    1.70276453214295
10    -1.63360010526121    1.84872922921405
11    0.103412043161049    1.91850200877146
12    1.15651664142905    1.91247118928509

♥欢迎点赞收藏♥