matlab 地形图 .grd,利用全球地形数据绘制EPR的地形图并添加测线

bat脚本

set PS="topo.ps"

set R=104:30W/104:00W/9:20N/10:00N

set topodata="topo15.grd"

gmt grdcut %topodata% -R%R% -GcutTopo.grd

gmt grdgradient cutTopo.grd -Nt -A50 -GcutTopo_i.grd

gmt makecpt -Crainbow  -T-3200/-2500 > colorTopo.cpt

REM 绘制底图

gmt set FORMAT_GEO_MAP=ddd:mm:ss

gmt psbasemap -R%R% -JM10c -Bxafg -Byafg -BWesN -Xc -Yc -K -P> %PS%

gmt grdimage cutTopo.grd -IcutTopo_i.grd -R -J -CcolorTopo.cpt -Q -O -K >>%PS%

REM 绘制colorbar

gmt psscale -DjCB+w10c/0.15i+o0/-0.6i+h -B100 -R -J -CcolorTopo.cpt  -G-3200/-2500 -O -K>>%PS%

REM gmtconvert 11.txt -o1,2,3,4 >trackline.dat

REM 绘制地震测线

gmt psxy EPR_Trackline_multi.dat -R -J -W0.8p -O -K >> %PS%

REM gmt psxy -R -J -T -O >>%PS%

REM gmt psconvert %PS% -A -Tg -P -Z

del .gmt* gmt.* cutTopo*.grd colorTopo.cpt tmp legend.txt