Gentle''的个人博客分享 http://blog.sciencenet.cn/u/tao164411096 九分耕耘,一份收获。

博文

GMT绘制站点分布图(详细又杂乱的)

已有 2419 次阅读 2018-7-27 17:02 |个人分类: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中的画笔有三个属性:笔宽、颜色和线型,三者用逗号分隔,即 <width>,<color>,<style> 
绘制底图类型为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 主要由三个属性控制:文字大小、字体、颜色。三个属性之间用逗号分隔,
即<size>,<fonttype>,<fill> 。三者均是可选的,但先后顺序不可乱。若其中任意一个属性被省略,
则使用该属性的默认值或上一次的设置值。字号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 <<END  >> $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 的一些东西(图片)

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备14006957 )

GMT+8, 2019-9-16 16:52

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部