|||
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
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-23 10:09
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社