最近需要處理關(guān)于不同分級(jí)顏色的散點(diǎn)，但是GMT實(shí)現(xiàn)不了，其具體的gmt代碼如下：

gmt?begin?symbols3?png,pdf ????

gmt?makecpt?-Chot?-T0/3/1 ????

gmt?plot?-R0/10/0/10?-JX10c/10c?-Baf?-Sc0.5c?-W1p,black?-C?<<?EOF?

2?3?0

5?6?1?

8?2?2?

EOF?

gmt?end?show

*參考資料：GMT學(xué)習(xí)筆記4-繪制符號(hào) - 知乎 (zhihu.com)

我使用GMT5，但是一直顯示語法錯(cuò)誤！于是我轉(zhuǎn)向matlab，最終完成了達(dá)到了目的，代碼如下：

A=textread('bd.txt');

LLZ.rg = LLZ.rg.*t;

LLZ.rg(LLZ.rg==0)=NaN;

rg_plot(LLZ)

xlim([70,140])

ylim([15,55])

% hold on

x = A(:,1);

y = A(:,2);

sz = linspace(1,269,length(x));

c = linspace(-100,100,length(x));

scatter(x,y,[],c,'filled')

具體實(shí)現(xiàn)效果如下：

歡迎與交流！！ 1565903211