linux下gmt4绘制站点分布图,科学网—GMT绘制站点分布图（详细又杂乱的） - 李建涛的博文...

Linux系统下GMT5.4.4 绘制站点图脚本  主要参考GMT参考手册5.4.2最后设置部分

#!/bin/bash

R=95/115/8/25

EOF 6.5i 控制图宽

-Jm1c 表示使用墨卡托投影，地图上的1度距离投影到画布上为1厘米

-Jm1:10000000 表示使用墨卡托投影，画布上的1 cm代表实际距离中的10000000 cm，即100 km

-JM15c 也表示使用墨卡托投影，整个地图的宽度是15厘米，地图的高度由 -R 和 -J 自动确定

-JX10c/5c 使用线性投影，地图的宽度是10厘米，高度为5厘米

EOF

J=M105/35/6.5i

PS=out.ps

#选择绘制图框类型 plain无铁轨  fancy有铁轨

gmt set MAP_FRAME_TYPE plain

#图中网格线的类型

gmt set MAP_GRID_PEN_PRIMARY 0.25p,gray,2_2:1

#

• ddd:mmF => 35:45W

• ddd:mmG => 35:45 W

• ddd:mm:ss => 40:34:24

• ddd.xxx => 36.250    图框的宽度为3p

gmt set FORMAT_GEO_MAP ddd:mm:ssF MAP_FRAME_WIDTH 3p

# 刻度标注的字体为7p

gmt set FONT_ANNOT_PRIMARY 7p

gmt set MAP_FRAME_PEN

GMT中的画笔有三个属性：笔宽、颜色和线型，三者用逗号分隔，即 ,,

绘制底图类型为plain时边框的画笔属性 [thicker,black]

表 5.1 GMT预定义画笔宽度名

线宽名    线宽     线宽名     线宽

faint     0        thicker    1.5p

default   0.25p    thickest   2p

thinnest  0.25p    fat        3p

thinner   0.50p    fatter     6p

thin      0.75p    fattest    12p

thick     1.0p     obese      18p

#刻度标注与轴之间距离0.1cm

gmt set MAP_ANNOT_OFFSET_PRIMARY 0.1c

#标注刻度线的长度 改为负值刻度线朝内

gmt set MAP_TICK_LENGTH_PRIMARY 3p

#FONT_LABEL 轴标签所使用的字体 [16p,Helvetica,black]，轴标注底部与轴标签顶部间的距离 [4p]

gmt set FONT_LABEL 8p,35 MAP_LABEL_OFFSET 4p

# 增加 -Q 可以取消底图  -L 绘制比例尺``-L 绘制比例尺，其中 g86/18 指比例尺在图中的位置是 86°E、18°N；+c30 指所绘为

30°N，中心经度之处的比例尺；+w1000 指比例尺长度为 1000km；+u 表示在文字标注后

显示长度单位，即 km；+f 表示比例尺样式为黑白相间的 fancy 样式。

#    -Bx0.2f0.1g0.2 -By0.1f0.1g0.2     分别设置刻度

gmt pscoast -J$J -R$R -G244/243/239 -S167/194/223 -B2f5g10 -Lg86/18+c30+w1000+u+f -K > $PS

## 此处使用两次 psxy 命令以实现图中所示的空心圆叠加实心圆的效果

gmt psxy site.dat -J$J -R$R -Sc0.2c -Gblack -K -O >> $PS

#Sc 使用五角星 0.4c 代表大小

gmt psxy site.dat -J$J -R$R -Sc0.4c -W0.5p,black -K -O >> $PS

# plot site name  +f 主要由三个属性控制：文字大小、字体、颜色。三个属性之间用逗号分隔，

即,, 。三者均是可选的，但先后顺序不可乱。若其中任意一个属性被省略，

则使用该属性的默认值或上一次的设置值。字号20p，字体20

# +j 文字对齐方式 (提前输入文档中的)

gmt pstext site.dat -J$J -R$R -F+f20p,20+j -Dj0.6c/0.3c -K -O >> $PS

#绘制鹰眼图部分

gmt psbasemap -R55/150/-4/55 -JM6.5 -B0 -X0 -Y9.8 -K -O >> $PS

gmt pscoast -R -JM6.5 -Df -N1 -W0.2 -A5000 -K -O >> $PS

gmt psxy CN-border-La.dat -J -R -W0.2p -O -K >> $PS

gmt psxy -R -JM6.5 -W1p,red  -O -L <> $PS

95 9

95 22

110 22

110 9

END

#ps转换JPG并且将分辨率改为500 图的边框为2cm

gmt psconvert -A2c -E500 out.ps

gmt psconvert -A -P out.ps

gmt psconvert -A -P -Tf out.ps

输出 f-PDF e-emf t-tiff

gmt psconvert -A2c -E500 -P -Tf -Te -Tt out.ps

转载本文请联系原作者获取授权，同时请注明本文来自李建涛科学网博客。

链接地址：http://blog.sciencenet.cn/blog-3391834-1126229.html

上一篇：注释多行

下一篇：gmt set 的一些东西(图片)