德克萨斯大学埃尔帕索地质研究所GMT实例之三

huolongguo

bat代码：

gmtset PAPER_MEDIA A2
rem Plot Border Around Map
echo -105.7 43.0 25 0 1 2 Topographic Map for CD-ROM 1999 | pstext -R-108.1/-103.5/34/43.2 -Jl-106.25/39/33/44.5/1:2000000 -X1.0 -K -P -V> cdtopo.ps

rem Generate color image
grdimage cdrom30.grd -Ccdrom1.cpt -R-108/-103.75/34.25/43 -Jl-106.25/39/33/44.5/1:2000000 -X1.0 -O -K -V>> cdtopo.ps
rem Put coastlines, national boundaries, rivers and scales and grid
pscoast -R -Jl -Dh -B1 -Ia -Na/5/0/205/96 -Lf-104.5/41/39/100 -S168/255/255 -W3 -K -O -V>> cdtopo.ps
rem Bring in coord files
psxy -: cdllee99.txt -R -Jl -P -Sc0.1/255/0/0 -O -K -V>> cdtopo.ps
psxy -: cdshlle99.txt -R -Jl -P -Sa0.2/0/205/205 -G0/205/205 -O -K -V>> cdtopo.ps
rem Begin plotting legend

echo 0.0 0.0>tmp
echo 0.2 7.7>>tmp
echo 3 7.7>>tmp
echo 3 0.0>>tmp
echo 0 0.0>>tmp
psxy tmp -Jx1.0 -R0/7/0/15 -G255/255/255 -L -W10 -Y0.1 -K -O -V>> cdtopo.ps
del tmp
rem #Plot scale
psscale -Ccdrom1.cpt -D0.8/3.8/7/0.4 -I -Ba600g300 -O -K -V>> cdtopo.ps
rem  #Plot globe in corner
pscoast -R0/180/-90/90 -JA-110/32/5 -B30g30 -Dc -A1000 -G0 -S255/255/255 -X4.5 -O -K -V>> cdtopo.ps
rem Plot yellow area on globe
echo -110 32>tmp
echo -106 32>>tmp
echo -106 49>>tmp
echo -110 49>>tmp
echo -110 32>>tmp
psxy tmp -JA -R -G255/255/0 -L -O >> cdtopo.ps
del tmp
ps2raster cdtopo.ps -A -P -Tg
del .gmt*

效果图：