GMT5中地形数据叠加色彩cpt

前段时间阅读文献时，发现有一类论文绘制的图片非常美观（Zhang et al.,2022）：

感谢论文作者Zhang博士提供的帮助，最终完成了类似的效果：

代码如下：

set PS ="NET_grav_vel.ps" 

set R=70/140/15/55    

set J=M14c 

REM gmt grdcut ETOPO1_Bed_g_gdal.grd -R%R% -Gchina.grd

REM gmt grdsample china.grd -Gresamp_china.grd -I0.01/0.01 

REM gmt grdgradient resamp_china.grd -A270 -Ne0.8 -Gchina.int     

REM xyz2grd  mass_g.txt  -R%R% -I0.1  -Gvel_grav.grd

grd2cpt vel_grav.grd -Cjet -S-50/50/0.1> vel_grav.cpt  

grdimage vel_grav.grd  -R%R% -J%J% -B60f60/30f30WenS -Cvel_grav.cpt -Ichina.int -Xc -Yc -K>%ps%   

gmt psbasemap -R%R% -J%J%  -BwesN  -Xc -Yc -O -K >>%PS%

pscoast -R%R% -J%J%  -Df -A300 -S0/191/255 -Lf107.5/30/36/200k+l+jright -I1/0.8p/0/229/238  -N1/0.4p -O -K>>%ps% 

gmt psvelo position.txt -R%R% -J%J% -Se0.1c/0.9/0 -A0.4c+e+p0.1p -Gred -W0.2p,red -X0c -Y0c -O >>%PS%

REM psscale -D3c/10.8c/5.5c/0.2ch -Cvel_grav.cpt -B0.1/:@~m@~Gal/yr: -E -O -K >> %ps%

ps2raster -E50 -P -Tf -A >>%ps%  

欢迎交流：1565903211

参考文献

Zhang, T.; Shen, Z.; He, L.; Shen, W.; Li, W. Strain Field Features and Three-Dimensional Crustal Deformations Constrained by Dense GRACE and GPS Measurements in NE Tibet. Remote Sens. 2022, 14, 2638. https://doi.org/10.3390/rs14112638