GMT小结.doc

1. psbasemap B -J -R -E -G -K -O -P -T -U -V -X -Y -Z -c map.ps ，表示可以缺省 -Rxmin/xmax/ymin/ymax or -Rxmin/ymin/xmax/ymaxr (加上r的区别) -X4i-Y4i 左下角原点右移上移4i，画新图时缺省表示离左下角偏移量（1i，1i），覆盖模式时缺省表示离原图左下角偏移（0，0）-Ba20f10g5:WSEN 20为刻度说明间隔，10黑白铁轨间隔，g5为绘制网格线的间隔。Ba20(m/s)缺省时为度，Ba20m表示刻度说明间隔为20分。:里为图框说明。:x: , :y: 形式为坐标轴说明(可以缺省)，:.title:形式为整个图说明，:,%: or :,-%:为刻度后缀，如10%,有-时，10与%间无空格，:=$: or :=-$:刻度的前缀。WSEN表示4个方向上的刻度轴，大写意味着标记该刻度轴。不同刻度轴的说明用 /隔开。(莫卡托投影好像不能标记坐标轴)。2. pscoast -R-90/-70/0/20 -JM6i -P -B5g5 -G0/0/255 -N1/0.5p -W0.25pmap.ps -A 去掉图中细节，比如绘出的一些小湖泊在大尺度图中显示的黑点，数字越大去掉的越多 -Df(h,i,l,c) 海岸线，国界线等的分辨率，从f-c递减 -G 陆地颜色 -Gr/g/b -I 绘制河流，-I1/0.2p/0/0/255,I1表示主要河流，后面是绘制河流的形状颜色参数 -S 水体颜色 -W 海岸线参数 -W0.1p/r/g/b -N 绘制国界线等，-N1/0.1p/r/g/b,N1为国界线，N2本来应该是省界线的，但是中国等国家的信息缺失，能画美国的州界线。 -E200/40 视点转角/倾角 (ps:参考可试运行pscoast -R70/150/15/60 -JB110/35/33/45/6i -Ba10f5:.Conic Projection: -N1/0.3p -N2/0.25p -A1000 -G200 -W0.25p/255/0/0 -E200/40 -I1/0.2p/0/0/255 -P1.ps) 3. echo 120.0 90.0 1.txt 把相关信息写入文本文档，以便直接在gmt命令流中修改，并可在每次命令流的末端删除,免得文件内产生过多文件。 当然，我们也可以直接在文件夹中创建一文件，把内容事先写进去。4. pscoast -R70/150/15/60 -JM6i -Ba10f5 -N1/0.3p -A500 -G200 -E200/40 -W0.25p/255/0/0 -I1/0.2p/0/0/255 -P -K2.ps echo 100 40 temp echo 100 45 temp echo 110 45 temp echo 110 40 temp (为举例生成一个temp文件) psxy -R -JM temp -W1/0/255/0 -O -K2.ps 绘制直线 psxy -R -JM temp -W1 -G255 -O -K2.ps 绘制方框，-G255为填充颜色 psxy -R -JM temp -W1t3_3:0p -O -K2.ps 绘制虚线，3_3:0分别为实线部分与空白部分的比例：初始相位 psxy -R -JM temp -Sc0.15 -W1 -G255/0/0 -O -K -V2.ps 用sc命令为绘制圆形，(Sa星号，Sf断层等) (试着分别运行以上4条psxy语句，看看区别) -N把那些超出边框的数据画上去。 Ps: Sa 星号, Sb 条棒,Sc 圆,Sd 菱形,Se 椭圆,Sf 断层,Sh 六边形,Si 倒三角,Sp 点,Ss 正方形,St 三角形,Sv 箭头,Sw 扇形,Sx 十字.5. psxy fault1.bln fault2.bln -R -JM -Sf0.1i/0.05iR:1.2 -P -O -W0.5p2.ps (正断层) psxy fault.bln -R -JM -Sf0.1i/0.05iRs:1.2 -P -O -W0.5p2.ps (走滑断层) Sf0.1i/0.05iRs:1.2i 0.1i为断层标记出现的间隔，0.05i为标记的大小,R(L)表示逆冲方向，l（r）s表示走滑。:1.2标记的相位。 一条命令可包含多个信息相同的断层数据文件，每个单独的文件生成一条断层。6. psvelo -R -JM 1.dat -Se0.2(箭柄长度)/0.95/0 -G90 -O -K -A0.01(柄宽)/0.1(头长)/0.02(头宽) 1.ps Se对应的数据格式，1，2为起点的经纬度，3，4为箭头经纬度上的分量，5，6为经纬度上的误差，7为误差的相关性，8为起点标注(可缺省)。 -G为箭头颜色，W1p/r/g/b误差椭圆边框信息，-E填充误差椭圆颜色,-L1p/r/g/b给箭头绘制边线。 -N把那些超出边框的数据画上去。在以后的命令中类似。7. echo 115 38(经纬度) 8(字大小) 0.0(逆时针旋转角度) 1(字体，1黑体，2斜体) 2(位置，2表示所写内容的中心出现在前面经纬度上) WUHAN(内容) text pstext text -R -JM -O -N -G0/0/0(字体颜色) 1.ps 给图中添加文字说明8. psscale -Ct2.cpt -D2.5i/-0.25i/4i/0.06ih -O -K -L -B20:x:/:y: -E 1.ps 绘制刻度尺 t2.cpt为色谱文件。2.5i/-0.25i是刻度尺中心位置，4i/0.06i是刻度尺长度和宽度，h表示刻度尺为横着的，缺省时为竖着的。 -L缺省时按照-B后面的数字给出数字说明的间隔，否则按cpt文件给出的间隔标记。 -B后面是刻度说明。 -E自己比较有没有的区别吧。 -A将-B选项中的说明置于刻度尺的另一侧 -M按照你的cpt文件将颜色灰度化 （该命令最好放到所有命令的最后，否则B命令将会影响其后的语句）9. makecpt -Crainbow -T-2/14/2 g.cpt 按照系统给出的色谱文件，从新指定范围，生成新的色谱文件 rainbow为系统自带的色谱，类似的还有cool,gebco等，其他可查询帮助。没事的话，可一个一个试着玩。 -T-2/14/2为确定的范围和间隔。 -I颜色信息反向。最低值对应系统色谱的最高值颜色。 -Z生成的色谱为连续的。 Rainbow: 紫,蓝,青,绿,黄,红(0/300) Cool: 蓝-紫 (0/1) Copper: 黑-红铜色(0/1) Gebco: 海水颜色(青-白)浅色调(-7000/0) Globe: 全球地形色调,高海拔浅色,水平绿色,海蓝紫(-10000/10000) Drywet: 干(褐)-湿(蓝) (0/12) Gray: 灰度(黑-白)(0/1) Haxby: 大地水准面,重力用(高海拔红白,水平青黄,海蓝)(0/32) Hot: 黑-红-黄-白 (0/1) Jet: 深蓝-浅蓝-白-黄-红 (0/1) No_green: 没有绿色,蓝-青-黄-红(-32/32) Ocean: 海洋,蓝-绿-白(-8000/0) Polar: 蓝-白-红 (-1/1) Red2green: 红-白-绿 (-1/1) Relief: 陆地,海洋.(高海拔黄白,水平绿褐,海蓝黑)(-8000/8000) Topo: 地形(绿-红-蓝-青-黄)(-7000/7000) Sealand: (蓝-绿-红)(-7000/3000) Seis: (红-黄-绿-蓝) 地震层析成像图(-1/1) Split: 和polar类似,黑代替白;蓝-黑-红(-1/1) Wysiwyg: (紫-蓝-青-绿 -黄-红) (0/20)10. 创建网格文件命令，xyz2grd，用已知的3D数据生成gmt专用的grd文件，由于机器问题，具体命令我也不懂。11. grdgradient in_grdfile -Gout_grdfile -A0 -Nt 根据你的网格文件生成一个光照的网格文件，使图产生立体感。 -A0(从北面来的光照) -A0/270(北面照向西面) 命令肯定没错，但是我从来没有成功过，暂时认为是机器原因。12. grdimage -Iintense.grd topo_asia.grd -JM8i -R70/150/15/60 -Ba10f5 -K -Cg.cpt 1.ps pscoast -R -JM -N1 -W1/0/0/0 -B -Dl -O -K1.ps intense.grd 为光照网格文件， topo_asia.grd为地形网格文件，g.cpt为海拔对应的色谱文件。13. psmeca -R -JM zycmt2_50.txt -Sm0.1i/4u -K -O -P -T0 -N -G255/0/0 cmt.psSm对应的数据类型x,y,depth,mrr,mtt,mff,mrt,mrf,mtf,exp,newX,newY,event_title 其他类型数据对应的命令虽然可以画出beach ball，但是画不出颜色(郑勇言). 缺乏深度数据时加-o，重要。14. grdraster 用包含30秒30秒高程数据的dbase文件覆盖原文件后，可以运行grdraster命令。即调用高程数据，生成gmt所需的grd格式文件。 例子：grdraster 32 -R80/100/25/40 -Gtibet_S.grd -I0.5m 由于30秒高程数据被分解成几十个小文件存储，如代号为32的文件包含范围R60/100/-10/40。在画地形图时往往需要将好几个高程文件拼接起来，或者只需要截取某个文件中的一部分，前者通过多次调用grdraster命令来实现，后者则通过grdraster命令中的-R范围实现,争取做到运行时间最短。-Goutfile.grd。-I后接数据的分辨率。 各个文件的范围，代号及分辨率，可以通过在Dos界面下直接运行grdraster查询得到。 大家可以试着运行： grdraster 32 -R80/100/25/40 -Gtibet_S.grd -I0.5m grdraster 41 -R80/100/40/45 -Gtibet_N.grd -I0.5m grdgradient tibet_S.grd -Gtibet_light1.grd -A180/90 -Nt grdgradient tibet_N.grd -Gtibet_light2.grd -A180/90 -Nt makecpt -Crelief -T-5000/5000/1000 -Z ti.cpt grdimage -Itibet_light1.grd tibet_S.grd JM5i -R80/