之前我講述了matlab中如何對散點進(jìn)行分級顯示matlab中對散點進(jìn)行分級顯示 - 嗶哩嗶哩 (bilibili.com)。現(xiàn)在我來介紹如何在gmt5中繪制散點圖。

此處用到的函數(shù)模塊為

gmt psxy -R%R% -J%J% -O -K point_label.txt -Sc0.15c -W0.5p,red -Cpoint.cpt -X0c -Y0c >>%ps%

其中，使用-G選項可以為符號填充單一顏色。如果要要符號顏色根據(jù)某個數(shù)值不同而變化，需要使用-C選項，-C后加cpt文件。其中的point_label.txt文件為三列屬性，前兩列為經(jīng)緯度信息，而第三列為顏色對應(yīng)的cpt值

2 3 0?

5 6 1?

8 2 2

下面為全部的代碼：

gmtset PLOT_DEGREE_FORMAT? ? ? ? ? ?ddd:mm:ssF

gmtset MAP_FRAME_TYPE plain

gmtset FRAM_WIDTH 0.04c

gmtset LABEL_FONT 2

gmtset LABEL_FONT_SIZE 8p

gmtset ANNOT_FONT_SIZE_PRIMARY 8p

gmtset MAP_FRAME_WIDTH 1p

gmtset MAP_ANNOT_OFFSET_PRIMARY 0.1c

set J=Cyl_stere/0/0/10c

set R=70/140/15/55

set PS="hydro.ps"

gmt makecpt -Cgray -T1000/10000/100 -I >2.cpt

gmt makecpt -Cpolar -T-100/100/25 >point.cpt

grdimage 1.grd? -R%R% -J%J% -B10f10/10f10WENS? -C2.cpt -Ichina.int -X2c -Y8c -K>%ps%??

pscoast -R%R% -J%J% -Dh -N1 -A1000 -I1 -Wthinnest,black -X0c -Y0c -O -K >>%ps%

gmt psxy -R%R% -J%J% -O -K point_label.txt -Sc0.15c -W0.5p,red -Cpoint.cpt -X0c -Y0c >>%ps%

gmt psxy tianshan.txt -J%J% -R%R%? -W0.5,red,- -X0c -Y0c -O -K >> %ps%

gmt psxy haihe.txt -J%J% -R%R% -W0.3,red,- -X0c -Y0c -O -K >> %ps%

echo 75 51 (a)| gmt pstext -J%J% -R%R% -F+f10p,5,black? -X0c -Y0c -N -O -K>>%PS%

gmt psscale -Dx0.5c/-01.2c+jTC+w2.5c/0.2c+h+e --FONT_ANNOT_PRIMARY=8p,Helvetica,black -X4.5c -Y7c -O -Cpoint.cpt -Ba50>>%ps%

ps2raster %PS% -A? -Tg? -E512