Linux系统下GMT5.4.4绘制站点图

#!/bin/bash

R=113.8/114.5/22/22.4

#    -JY<lon>/<lat>/<width>  -Jy<lon>/<lat>/<scale>  <lon> 是中心经线，<lat> 是标准纬线 

#-jm6c 一度为6cm ，-JM6c 宽度为6cm  ，Ju：通用横向 Mercator(UTM) 投影  

#-JU<zone>/<width> Ju<zone>/<scale

J=m6c

PS=out.ps

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

gmt set MAP_FRAME_TYPE fancy

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

gmt set MAP_FRAME_PEN 0.1p,solid

#图中网格线的类型

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

gmt set FORMAT_GEO_MAP ddd.xxF

#图框的宽度为3p

gmt set MAP_FRAME_WIDTH 3p

# 刻度标注的字体为7p

gmt set FONT_ANNOT_PRIMARY 7p

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

gmt set MAP_ANNOT_OFFSET_PRIMARY 0.1c

#标注刻度线的长度 -3p 改为负的刻度朝内

gmt set MAP_TICK_LENGTH_PRIMARY -3p

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

gmt set FONT_LABEL 8p,35 MAP_LABEL_OFFSET 4p

#使用A4纸绘图

gmt set PS_MEDIA=A4 

# a 刻度标注值得间隔，f 刻度线间隔 g 网格线间隔 -Bx -By 表示分别对xy轴进行控制 WSEN 大写表示显示刻度和标注，小写表示不显示刻度标注   -Q 控制不显示底图，若想显示底图去掉 -Q 并加上 -Df 即可

gmt pscoast -J$J -R$R -Bx0.2f0.1g0.2 -By0.2f0.1g0.2 -BWSen -Q -K  > $PS

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

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

gmt psxy site.dat -J$J -R$R -Sc0.15c -W1p,black -K -O >> $PS

# 绘制站点名，-F控制 其中+f主要有三个属性控制，文字<大小>，<字体>，<颜色>，三个属性之间用逗号  #隔开，三者均可选，但先后顺序不能乱，若其中一个被省略则使用该属性的默认值或上一次设置值

#  +j 表示对齐方式（在输入文档中设置），具体参考博客《gmt set的一些东西》  -Dj<dx>/<dy> 表示沿着#   +j对齐的方向便宜

gmt pstext site.dat -J$J -R$R -F+f7p,2+j -Dj0.15c/0.15c -O >> $PS

#输出图片的边缘留0.01cm的白 ， dip为700 的jpg格式

gmt psconvert -A0.01c -E700 out.ps

#输出j-jpg f-pdf e-emp t-tiff 

gmt psconvert -A0.01c -E600 -P -Tj -Tf -Te -Tt out.ps

输入文件site.date建立模板：

# 经度    纬度   对齐方式  城市名

113.894   22.114   ML   HKNP

114.228   22.108   ML   HKOH

若命令中直接使用-F+f+j 则需要加入一列字体大小，字号参数

# 经度    纬度    字体    对齐方式  城市名

113.894   22.114   7p,2   ML   HKNP

绘制了fancy和plain有底图和无底图的效果，如图：