1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
REM GMT EXAMPLE 40
REM $Id$
REM
REM Purpose: Illustrate line simplification and area calculations
REM GMT progs: psbasemap, pstext, psxy, gmtsimplify, gmtspatial
REM DOS calls: gawk, del
REM
echo GMT EXAMPLE 40
set ps=example_40.ps
gmt spatial GSHHS_h_Australia.txt -fg -Qk > centroid.txt
gmt psbasemap -R112/154/-40/-10 -JM5.5i -P -K -B20 -BWSne+g240/255/240 -Xc > %ps%
gmt psxy GSHHS_h_Australia.txt -R -J -O -Wfaint -G240/240/255 -K >> %ps%
gmt psxy GSHHS_h_Australia.txt -R -J -O -Sc0.01c -Gred -K >> %ps%
gmt simplify GSHHS_h_Australia.txt -T500k > T500k.txt
gmt spatial GSHHS_h_Australia.txt -fg -Qk | gawk "{printf 'Full area = %.0f km@+2@+\n', $3}" > area.txt
gmt spatial T500k.txt -fg -Qk | awk '{printf "Reduced area = %.0f km@+2@+\n", $3}' > area_T500k.txt
gmt psxy -R -J -O -K -W1p,blue T500k.txt >> %ps%
gmt psxy -R -J -O -K -Sx0.3i -W3p centroid.txt >> %ps%
echo 112 -10 T = 500 km | gmt pstext -R -J -O -K -Dj0.1i/0.1i -F+jTL+f18p >> %ps%
gmt pstext -R -J -O -K area.txt -F+14p+cCM >> %ps%
gmt pstext -R -J -O -K area_T500k.txt -F+14p+cLB -Dj0.2i >> %ps%
gmt psbasemap -R -J -O -K -B20+lightgray -BWsne+g240/255/240 -Y4.7i >> %ps%
gmt psxy GSHHS_h_Australia.txt -R -J -O -Wfaint -G240/240/255 -K >> %ps%
gmt psxy GSHHS_h_Australia.txt -R -J -O -Sc0.01c -Gred -K >> %ps%
gmt simplify GSHHS_h_Australia.txt -T100k > T100k.txt
gmt spatial T100k.txt -fg -Qk | gawk "{printf 'Reduced area = %.0f km@+2@+\n', $3}" > area_T100k.txt
gmt psxy -R -J -O -K -W1p,blue T100k.txt >> %ps%
gmt psxy -R -J -O -K -Sx0.3i -W3p centroid.txt >> %ps%
echo 112 -10 T = 100 km | gmt pstext -R -J -O -K -Dj0.1i/0.1i -F+jTL+f18p >> %ps%
gmt pstext -R -J -O -K area.txt -F+14p+cCM >> %ps%
gmt pstext -R -J -O -K area_T100k.txt -F+14p+cLB -Dj0.2i >> %ps%
gmt psxy -R -J -O -T >> %ps%
REM Clean up
del centroid.txt area*.txt T*.txt
|
