GMT绘制产状极坐标投影图

线产状数据投影到吴氏网上是一个点，用Stereonet可以很好的完成。
GMT的-JP极坐标线性投影也可以很方便的完成，如下：

data=sigma_lon_lat.txt
gmt begin
gmt figure polar png,pdf A0.1c,E600
gmt set FORMAT_GEO_MAP=+D
gmt set FONT_TAG 8p,4
gmt subplot begin 1x3 -Ff15c/6c -A'(a)' -M0.1c/0.1c/0.2c/0.2c # -M左/右/下/上 的空白距离
gmt subplot set 0,0 -A'@~\163@~1' 
    # +a表示自北方向顺时针旋转,+f反转径向即圆点处为90度
    gmt basemap -R/0/360/0/90 -JP?+a+f -B+gwhite -Bxg90 -Bya30g30 --MAP_FRAME_PEN=0.5p 
    gmt plot $data -i2,3 -Sc0.01c -W0.1p
gmt subplot set 0,1 -A'@~\163@~2' 
    gmt basemap -R/0/360/0/90 -JP?+a+f -B+gwhite -Bxg90 -Bya30g30 --MAP_FRAME_PEN=0.5p
    gmt plot $data -i4,5 -Sc0.01c -W0.1p
gmt subplot set 0,2 -A'@~\163@~3' 
    gmt basemap -R/0/360/0/90 -JP?+a+f -B+gwhite -Bxg90 -Bya30g30 --MAP_FRAME_PEN=0.5p
    gmt plot $data -i6,7 -Sc0.01c -W0.1p
gmt subplot end
gmt end

效果图：