使用GMT绘制花岗岩( Na2O + K2O) -SiO2 分类

花岗岩常用的( Na₂O + K₂O) - SiO₂分类图解，原作者Middlemost（1994）在文献中已经提供了图解各交点的坐标，如下图所示：

花岗岩分类图解坐标数据
所以，可比较方便的直接拿来用GMT进行花岗岩分类图解的绘制。
当然，即使原作者未提供各交点数据，也可以使用一些矢量化软件（如WebPlotDigitizer)获得图解中的坐标数据。

下面为使用GMT的绘制示例：

gmt begin
gmt set MAP_FRAME_PEN 0.5p 
gmt set FONT_ANNOT_PRIMARY 6p,4
gmt set FONT_LABEL 7p,4
gmt figure NaK-Si pdf,png E800,A0.2c
# 底图
gmt basemap -JX8c/6c -R30/90/0/19 -Bxa10+l"SiO@-2@-(%)" -Bya2+l"Na@-2@-O+K@-2@-O(%)"
gmt plot NaK-Si.txt -W0.2p
# 投点
gmt plot -Sc0.1c -Gred -l@~'\147\144'@~Pt@-1@-+jRT+p0.2p << EOF
# γδPt1
# SiO2 Na2O+K2O
72.97   7.92
74.37   8.27
73.9    6.02
73.14   8.54
74.45   7.85
EOF
gmt plot -Sc0.1c -Gblue -l@~'\150\147\142'@~Pt@-1@- << EOF
# ηγβPt1
# SiO2 Na2O+K2O
70.88   8.97 
72.05   7.89 
70.73   8.06 
68.77   8.20 
EOF
gmt plot -Sc0.1c -Ggreen -l@~'\150\147\142'@~mPt@-1@- << EOF
# ηγβmPt1
# SiO2 Na2O+K2O
65.63   7.90 
65.70   8.72 
71.48   7.99 
EOF
gmt plot -Sc0.1c -Gcyan -l@~'\150\147'@~Pt@-1@- << EOF
# ηγPt1
# SiO2 Na2O+K2O
65.64   6.87 
72.40   8.25 
69.88   8.08  
EOF
gmt plot -Sc0.1c -Gmagenta -l@~'\144\150\157'@~Pt@-1@- << EOF
# δηοPt1
# SiO2 Na2O+K2O
54.16   5.39 
54.70   5.22 
54.08   5.72 
EOF
gmt end

NaK-Si.png