# COLOR_MODEL = HSV
0 300 1 1 1 299 1 1 1 299 1 1 2 298 1 1 2 298 1 1 3 297 1 1 ???????
例中设z值由–2到14，分别把–2、14与rainbow.cpt文件的0、300的颜色对应，并将HSV变换为RGB格式，得g.cpt内容如下：
–2 255 0 255 0 255 0 255（由–2到0间隔，RGB=255/0/255） 0 73 0 255 2 73 0 255 ???????
12 255 0 0 14 255 0 0
B 0 0 0 （背景，即<–2，RGB=0/0/0）
F 255 255 255 （前景，即>14，RGB=255/255/255） N 128 128 128 （空白， RGB=128/128/128）
用grdxyz读出的geosat5.grd的内容如下，为xyz格式的网格文件 164.333 47.5833 –1.433 164.417 47.5833 –1.44 164.5 47.5833 –1.451 ???????
grdimage命令用geosat5.grd 和g.cpt来绘制彩色充填的等值区图。 t.cpt内容如下
–7.438 255 0 255 –6.08677 127 0 255 （z由–7.438到–6.08677，RGB由255/0/255
渐变到127/0/255）
–6.08677 127 0 255 –5.75589 0 0 255 ??????? –4.15904 255 127 0 2.804 255 0 0 ???????
用grdxyz读出的gtopo5.grd的内容如下，为xyz格式的网格文件 159.833 47.5833 –5.244 159.917 47.5833 –5.249 160 47.5833 –5.25 ???????
用grdxyz读出的etopo5_int.rgd的内容如下，为xyz格式的网格文件，z值为由etopo5.grd计算的向北面的方向梯度。物理意义为太阳在正北方向时，不同的地形的亮度。 159.833 47.5833 0.0519122 159.917 47.5833 0.0118648 160 47.5833 –0.000182332 ???????
––––––––––––––––––––––––––––––––––––––– Figure 6.4: 3–D透视网线图 批处理文件如下
echo –10 255 0 255 > zero.cpt echo 0 100 10 100 >> zero.cpt
grdcontour geoid.grd –Jm0.45i –E60/30 –R195/210/18/25 –C1 –A5 –G4i –K –P –X1.5i –Y1.5i –U/–1.25i/ –1.25i/\
pscoast –Jm –E60/30 –R –B2/2NEsw –G0 –O –K >> example_04.ps
26
echo 205 26 0 0 1.1 | psxyz –Jm –E60/30 –R –SV0.2i/0.5i/0.4ii –W1p –O –K –N >> example_04.ps echo 205 29.2 36 –90 1 5 N | pstext –Jm –E60/30 –R –O –K –N >> example_04.ps
grdview topo.grd –Jm –Jz0.34i –Czero.cpt –E60/30 –R195/210/18/25/–6/4 –N–6/200/200/200 –Qsm –O –K –B2/2/2:\（km）\–Y2.2i >> example_04.ps
echo 3.25 5.75 60 0.0 33 2 H@#awaiian@# R@#idge | pstext –R0/10/0/10 –Jx1i –O >> example_04.ps del zero.cpt
部分命令注释
echo 用echo建立一个灰度.cpt文件zero.cpt grdcontour –E60/30：视点转角/水平倾角 echo 205 26??
205 26 0 0 1.1作为psxyz的输入值，即由标准输入给出画符号（SV）的起点坐标（205,26,0），方向（0），长度（1.1i）。–Jm –E –R都同上，–SV箭体宽箭/头长箭/头宽，–N符号可画在图框外；–E60/30：视点转角/倾角。
echo 205 29.2 ??
205 29.2 36 –90 1 5 N 作为pstext的标准输入，依次为X坐标，Y坐标，字体大小（p），角度，字符集编号，字符对齐选择（5为中中对齐），字符内容。
grdview 设Z轴的刻度为每单位0.34i；–R195/210/18/25/–6/4：x1/x2/y1/y2/z1/z2；–N–6/200/200/200：z=–6为网格图底面，灰色；–Qsm：设置3D图类型，s表面图（类似Surfer软件的最小曲率算法），m网格图，i影象图；–B2/2/2:\（km）\：x,y,z刻度间隔均为2，只标注Z轴轴名 echo 3.25 5.75 ?? 3.25 5.75 60 0.0 33 2 H@#awaiian@# R@#idge作为pstext的标准输入，@#字母大小写转换开关 del zero.cpt 删除色谱文件
––––––––––––––––––––––––––––––––––––––––
Figure 6.5: 3–D 照明表面图 批处理文件如下
grdmath –R–15/15/–15/15 –I0.3 X Y HYPOT DUP 2 MUL PI MUL 8 DIV COS EXCH NEG 10 DIV EXP MUL = sombrero.grd
echo –5 128 5 128 >gray.cpt
grdgradient sombrero.grd –A225 –Gintensity.grd –Nt0.75
grdview sombrero.grd –JX6i –JZ2i –B5/5/0.5SEwnZ –N–1/255/255/255 –Qs –Iintensity.grd –X1.5i –K –Cgray.cpt –R–15/15/–15/15/–1/1 –E120/30 –U/–1.25i/–0.75i/\echo 4.1 5.5 50 0 33 2 z（r）= cos（2@~p@~r/8）* e@+–r/10@+ | pstext –R0/11/0/8.5 –Jx1i –O >> example_05.ps
27
部分命令注释
grdmath 把方程式Z（γ）=cos（2πγ/8）
γ/10
×e–，用RPN语法写出（[（HYPOT DUP 2 MUL PI MUL 8 DIV）COS）（（（EXCH NEG）10 DIV）EXP）MUL ]，计算.grd文件sombrero.grd
echo –5 128??建立灰度文件gray.cpt grdgradient 由sombrero.grd 计算225°的方向梯度文件intensity.grd。–Nt 见 Figure 6.2注释。
grdview 设Z轴高度为6i；xyz坐标的刻度单位分别为5、5、0.5；–N–1/255/255/255：z= –1为网格图底面，白色；
–Qs：设置3D图类型为表面图（类似Surfer软件的最小曲率算法）；光照强度由文件intensity.grd给出； 灰度由灰度文件gray.cpt给出；–E120/30：视点转角/倾角
γ
echo 4.1 5.5?? 作为pstext的标准输入，@#字母大小写转换开关。写图名“Z(γ)=cos(2πγ/8)*e–/10 ”
–––––––––––––––––––––––––––––––––––––––– Figure 6.6: 两类直方图 批处理文件如下
psrose fractures.d –A10r –S1.8in –U/–2.25i/–0.75i/\6 in Cookbook\–P –G0 –R0/1/0/360 –X2.5i –K –B0.2g0.2/30g30 > example_06.ps pshistogram –Ba2000f1000:\（m）\ v3206.t –R–6000/0/0/30 –JX4.8i/2.4i –G200 –O –Y5.5i –X–0.5i –L0.5p –Z1 –W250 >> example_06.ps del .gmt*
部分命令注释
psrose –A10r：玫瑰花图分区为10度；–S1.8in：半径1.8i，投影标准化为0～1； –R0/1/0/360：极坐标，半径方向0～1，圆周方向0～360度
pshistogram 直方图，–B：x轴标注间隔、刻度、轴名，y轴标注间隔、刻度、轴名，图名
––––––––––––––––––––––––––––––––––––––– Figure 6.7: 一个典型的位置图
pscoast –R–50/0/–10/20 –JM9i –K –GP300/26 –Dl –W0.25p –B10 –U\7 in Cookbook\> example_07.ps
psxy1 –R –JM –O –K –M fz.xy –W0.5pta >> example_07.ps
gawk \$1–360.0, $2, $3*0.01}\quakes.xym | psxy –R –JM –O –K –H1 –Sci –G255 –W0.25p >> example _07.ps
psxy2 –R –JM –O –K –M isochron.xy –W0.75p >> example_07.ps psxy3 –R –JM –O –K –M ridge.xy –W1.75p >> example_07.ps
28
echo –14.5 15.2 > tmp echo –2 15.2 >> tmp echo –2 17.8 >> tmp echo –14.5 17.8 >> tmp
psxy4 –R –JM –O –K –G255 –W1p –A tmp >> example_07.ps echo –14.35 15.35 > tmp echo –2.15 15.35 >> tmp echo –2.15 17.65 >> tmp echo –14.35 17.65 >> tmp
psxy5 –R –JM –O –K –G255 –W0.5p –A tmp >> example_07.ps
echo –13.5 16.5 | psxy –R –JM –O –K –Sc0.08i –G255 –W0.5p >> example_07.ps
echo –12.5 16.5 18 0 6 5 ISC Earthquakes | pstext –R –JM –O –K >> example_07.ps echo –43 –5 30 0 1 6 SOUTH > tmp echo –43 –8 30 0 1 6 AMERICA >> tmp echo –7 11 30 0 1 6 AFRICA >> tmp
pstext –R –JM –O –S0.75p –G255 tmp >> example_07.ps del .gmt*
部分命令注释
pscoast –GP300/26：陆地填充用模板26，反相，深灰（300）底色。 文件fz.xy内容如下：
> 629 OCEAN: 2 FEATURE: 6 NAME: KANE 284.74011 35.04543 284.82080 34.99857 ???????
> 445 OCEAN: 2 FEATURE: 6 NAME: FIFTEEN TWENTY 298.96030 18.06335 299.04568 18.01421 ???????
文件fz.xy为区内大洋转换断层的数据，为多段数据格式（–M）。每一段数据由一行“>”开头的文字和一系列x、y坐标组成，用来绘制一条断层，文字行包含了断层名字和x、y坐标组数等信息。 Psxy1 用文件fz.xy的数据来绘制多条转换断层 文件quakes.xym内容如下 lon lat mag1 338.80 –8.4 4.8 346.68 –8.36 4.4 ???????
gawk命令由文件quakes.xym中读取大洋中脊地震的x、y坐标和震级，作为psxy绘制地震分布图的输入数据。
文件isochron.xy和ridge.xy分别为区内“等时”断层和大洋中脊断层的数据，数据格式与fz.xy相同。 psxy2 用文件isochron.xy 的数据来绘制“等时”断层 psxy3 用文件ridge.xy 的数据来绘制大洋中脊断层
29
echo –14.5 15.2 ??4行echo语句给文件tmp赋值 psxy4 用tmp数据绘制矩形内框
echo –14.35 15.35 ??4行echo语句给文件tmp重新赋值 psxy5 用tmp数据绘制矩形外框
echo –13.5 16.5 | psxy ??在矩形框内绘制地震震中图例（圆圈），echo给出圆心位置
echo –12.5 16.5 18 0 6 5 ISC Earthquakes | pstext?? 在图例（圆圈）右面写“ISC Earthquakes”，echo后面参数依次为x坐标、y坐标、字符大小（字号）、角度（0）、字体（6）和对齐方式，5表示坐标值为字串的中心。
echo –43 –5??3行echo语句给文件tmp重新赋值，pstext 用tmp数据分别写“SOUTH”，“AMERICA”，“AFRICA” 3个字串。对齐方式6表示坐标值为字串右边中点。
–––––––––––––––––––––––––––––––––
Figure 6.8: 3–D 柱状图 grd2xyz topo.grd > tmp.d
psxyz tmp.d –B1/1/1000:\（m）\–R–0.1/5.1/–0.1/5.1/–5000/0 –JM5i –JZ6i –E200/30 –So0.0833333ub–5000 –P –U\–W0.25p –G240 –K > example_08.ps echo 0.1 4.9 24 0 1 9 This is the surface of cube | pstext –R –JM –JZ –Z0 –E200/30 –O >> example_08.ps
部分命令注释
grd2xyz 把2进位的.grd文件topo.grd转变为ASCII文件（x、y、z）tmp.d tmp.d文件内容如下： 0 5 –2823 0.0833333 5 –3027 0.166667 5 –3195 0.25 5 –3384
psxyz 读数据文件tmp.d，设定x、y、z轴刻度标注的间隔，z轴标注轴名“Topography（m）”，图名“ETOPO5”，标注WSZ；
设定x、y、z的范围，墨卡托投影边长5i，Z轴长6i，–So0.0833333ub–5000：画由基点到z的3D柱，柱宽0.083个x轴单位，基点设为–50000；–E200/30 视点转角/倾角。
echo 0.1 4.9 ?? 在立方体上表面标注文本
––––––––––––––––––––––––––––––––––––––––
Figure 6.9: 沿测线变化的时间–序列图Time–series as ``wiggles'' along a track
pswiggle all.xys –R185/250/–68/–42 –U\9 in Cookbook\–K –Jm0.13i –Ba10f5 –G0 –Z2000 –W0.25p –M –S240/–67/500/@~m@~rad > example_09.ps psxy1 –R –Jm –O –K ridge.xy –W1.75p >> example_09.ps psxy2 –R –Jm –O –K –M fz.xy –W0.5pta >> example_09.ps REM Plot labels
gawk \（NF = 5）print $3, $4, 10, 50, 1, 7, $2}\all.xys | pstext –R –Jm –D–0.05i/–0.05i –O >>
30