linux下gmt4绘制站点分布图,科学网—GMT绘制站点分布图（整理后的） - 李建涛的博文...

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

#!/bin/bash

R=113.8/114.5/22/22.4

#    -JY//  -Jy//   是中心经线， 是标准纬线

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

#-JU/ Ju/

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/ 表示沿着#   +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有底图和无底图的效果，如图：

转载本文请联系原作者获取授权，同时请注明本文来自李建涛科学网博客。

链接地址：http://wap.sciencenet.cn/blog-3391834-1126509.html

上一篇：gmt set 的一些东西(图片)

下一篇：MATLAB绘图设置