GMT绘制三维柱状图

huolongguo

bat代码：

REM Purpose: Make a 3-D bar plot
REM GMT progs: blockmedian, pstext, psxyz
REM DOS calls: echo, del
REM

echo Make a 3-D bar plot
blockmedian topo.dat -R90/100/35/40 -I10m -V > tmp.10m

数据太密，使用blockmedian按10分的间隔抽稀
echo 2 16.5 15 23 1 1 This is the surface of cube | pstext -R0/12/0/20 -Jx1 -K -V -P>cdm_3d.ps
psxyz tmp.10m -B2/2/1000:"Topography(m)"::".Qaidam Basin Topo":wsNEZ+ -R90/100/35/40/0/7000 -JM12 -JZ15 -E30/45 -So0.1ub0 -P -U"xiaoqiang" -Wthinnest -Glightblue -O -V >> cdm_3d.ps

wsNEZ后面加+是为了显示立方体轮廓线
del .gmt* tmp.*

效果图：