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.结果