前段時(shí)間閱讀文獻(xiàn)時(shí)，發(fā)現(xiàn)有一類論文繪制的圖片非常美觀（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

參考文獻(xiàn)

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