GMT绘制图中图

 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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77


#!/bin/bash

gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF

gmtset ANOT_FONT_SIZE 8p #主注释文字大小

gmtset ANNOT_FONT_SIZE_SECONDARY 6p #次注释大小 

gmtset LABEL_FONT_SIZE 6p

gmtset BASEMAP_TYPE plain #fancy

PS=sichuan.ps
R=102/105/31.5/35
J=X3.5id/4id
PH2=~/gmtdir/jiuzhai
PH=~/gmtdir/grdfile
Rg=96/110/26/36
Jg=X1id

#grdcut $PH/china.grd -Gsichuan.grd -R$R -V

grd2cpt sichuan.grd  -Ctopo -S0/6000/50 -Z > colors.cpt

psxy -R$R -J$J -T -K > $PS

psbasemap -R$R -J$J -B1/1WS -K -Xc -Yc> $PS
psbasemap -R$R -J$J -Ben -K -O >> $PS

grdimage  sichuan.grd  -R$R -J$J -Ccolors.cpt -K -O >> $PS

#sichuan.txt 四川省界 fault_china_2010.dat 中国断层数据
psclip sichuan.txt -R$R -J$J -K -O  -P >> $PS
psxy -R$R -J$J fault_china_2010.dat  -W1.5/purple -M -O -K >> $PS
psclip -C -O -K -P >> $PS

#faultline.txt 发震断层数据
psxy -R$R -J$J faultline.txt  -W4/red -M -O -K >> $PS

echo 103.82 33.20                               | psxy   -R$R -J$J -Sa0.25 -W1/blue -Gred -O -K >> $PS

echo 103.82 33.20 2.0 156 79 -9 3.0 103.5 33.50 | psmeca -R$R -J$J -Sa0.8 -Gred  -Ewhite  -O -K -W1 -L1  -C >> $PS

#mid.dat 断层示意图
gawk '{print $1,$2}'  mid.dat                   | psxy -R$R -J$J -Sc0.05 -W0.1/red  -Gred   -O -K   >> $PS

psxy -R$R -J$J -O -K sichuan.txt  -W2/black -M >> $PS

#地形高度色标
psscale  -Ctmp.cpt   -D2.95i/3.77i/0.9009i/0.08ih    -B50          -O -K  >> $PS
echo 104.44 34.85 8 0 4 1 'KM' | pstext -R$R -J$J -W  -O -K >> $PS

psscale  -Ccolors.cpt -D3.79i/1i/2i/0.1i          -Ba1500f500/:"Elevation(m)": -O -K -Alc >> $PS

echo 103.93 33 6 -65 4 1 Huya Fault       | pstext -R$R -J$J -W  -O -K  -Gred >> $PS
echo 103.47 32.7 6 285 4 1 Minjiang Fault | pstext -R$R -J$J -W  -O -K  -Gred >> $PS
echo 102.7 33.83 6 -8 4 1 Tazang Fault    | pstext -R$R -J$J -W  -O -K  -Gred >> $PS

#station.txt 台站名
gawk '{print $3,$2}' station.txt              |   psxy -R$R -J$J -St0.2 -W1.5/black -Gblack -O -K >> $PS
gawk '{print $3,$2,6,0,4,1,$1}' stationnm.txt | pstext -R$R -J$J -W255/250/250  -O -K >> $PS


pslegend -Dx2.8i/0i/0.7i/0.3i/BL  -R$R -J$J -O -K -L1  -Gwhite << EOF >> $PS
C black
S 0.1c t 0.2 black 0.5p/black 0.2c Station
S 0.1c a 0.2 red 0.5p/red    0.2c Epicenter
EOF

#小图 bound.txt 区域边界(四个角点坐标)
psbasemap -R$Rg -J$Jg -B5/5NE  -Gwhite -K -O --BASEMAP_TYPE=plain >> $PS
grdimage  $PH/china.grd  -R$Rg -J$Jg -Cchina.cpt -O -K>> $PS
psxy -R$Rg -J$Jg -O -K $PH2/sichuan.txt  -W1/black -M >> $PS
echo 103.82 33.20 | psxy -R$Rg -J$Jg  -Sa0.25 -W1/blue -Gred -O -K >> $PS
psxy bound.txt -R$Rg -J$Jg  -W1p,blue -O -L >> $PS

ps2raster $PS -A -P -Tg
Contents
