GMT5示例_在南极底图上画点

1.目标

在南极背景地图上绘制点。

注意：需要使用lambertAzimuthal投影。

2.代码

需要准备的数据文件：

文件名：data.txt（经度、纬度）

87.500  -76.200
114.500 -71.350
131.200 -76.570
123.140 -75.180
113.150 -71.260
151.120 -73.010

代码：

# ubuntu
R=0/360/-90/-60         # 经度从0-360,纬度从南极点到南纬60
J=A0/-90/5i     # 投影方式,lambertAzimuthal，0表示经线起点(正上方是0),-90表示极点为南极，5i表示纸张大小？
data=data.txt   # 数据文件
ps=map.ps       # 结果文件

# ps文件头
gmt psxy -J$J -R$R -T -K > $ps

# 海岸线 -W绘制shoreline -G干区域 -S湿区域
gmt pscoast -R$R -J$J -Df -Baf -W1/0.25p,white -Gwhite -Sgrey -K -O >> $ps

# 加点,使用了 -S 选项，则表示要绘制符号，-Sc ：圆， -G 设置符号或多边形的填充色
gmt psxy $data -R$R -J$J -Sc0.08c -Gred -K -O >> $ps

# ps文件尾
gmt psxy -R$R -J$J -T -O >> $ps

# 保存为图片,生成一个与结果文件同名的png
gmt psconvert $ps -Tg -P -E1200 -A

3.结果

posted @ 2017-04-16 19:55  ssh_alitheia  阅读(260)  评论(0)  编辑  收藏  举报