使用GMT绘制花岗岩( Na2O + K2O) -SiO2 分类
2022-04-02 本文已影响0人
cugliming
花岗岩常用的( Na2O + K2O) - SiO2分类图解,原作者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