区域底图时注意！ -Df 平滑划分区域的线

gmt pscoast -J$J -R$R -Ba0.2f0.1g0.2 -BWSne -A0 -Glightyellow -Slightblue -Wthinnest -I0.1 -W0.1 -Df -X3c -Y12c -P -K > $PS

#!/bin/bash

R=113.6/114.5/22/22.4

J=M105/35/5c

PS=out.ps

gmt set MAP_FRAME_TYPE fancy

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

gmt set FORMAT_GEO_MAP ddd.xxF

gmt set MAP_FRAME_WIDTH 2.5p

gmt set FONT_ANNOT_PRIMARY 6p

gmt set MAP_ANNOT_OFFSET_PRIMARY 0.05c

gmt set MAP_TICK_LENGTH_PRIMARY 4p

#gmt set FONT_TITLE 6p

#gmt set FONT_LABEL 4p,35 MAP_LABEL_OFFSET 1p

# 增加 -Q 可以取消底图

gmt pscoast -J$J -R$R -Ba0.2f0.1g0.2 -BWSne -A0 -Glightyellow -Slightblue -Wthinnest -I0.1 -W0.1 -Df -X3c -Y12c -P -K > $PS

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

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

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

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

 # plot site name  +f size  Dj0.6c   right 0.6cm

gmt pstext site.dat -J$J -R$R -F+f8p,10+j -Dj0.1c/0.1c -O >> $PS

#绘制鹰眼图部分

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

#输出 f-PDF e-emf t-tiff

gmt psconvert -A0.5c -E600 -P -Te -Tt out.ps

rm gmt*