|||
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有底图和无底图的效果,如图:
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 21:37
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社