GMT小结.doc
- 文档编号:1398946
- 上传时间:2022-10-22
- 格式:DOC
- 页数:5
- 大小:79.50KB
GMT小结.doc
《GMT小结.doc》由会员分享,可在线阅读,更多相关《GMT小结.doc(5页珍藏版)》请在冰豆网上搜索。
1.psbasemap–B<>-J<>-R<>[-E][-G][-K][-O][-P][-T][-U][-V][-X][-Y][-Z][-c]>map.ps[],<>表示可以缺省
-Rxmin/xmax/ymin/ymaxor-Rxmin/ymin/xmax/ymaxr(加上r的区别)
-X4i『-Y4i』左下角原点右移『上移』4i,画新图时缺省表示离左下角偏移量(1i,1i),覆盖模式时缺省表示离原图左下角偏移(0,0)-Ba20f10g5:
:
WSEN20为刻度说明间隔,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.25p>map.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:
."ConicProjection":
-N1/0.3p-N2/0.25p-A1000-G200-W0.25p/255/0/0-E200/40-I1/0.2p/0/0/255-P>1.ps)
3.echo120.090.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-K>2.ps
echo10040>temp
echo10045>>temp
echo11045>>temp
echo11040>>temp(为举例生成一个temp文件)
psxy-R-JMtemp-W1/0/255/0-O-K>>2.ps绘制直线
psxy-R-JMtemp-W1-G255-O-K>>2.ps绘制方框,-G255为填充颜色
psxy-R-JMtemp-W1t3_3:
0p-O-K>>2.ps绘制虚线,3_3:
0分别为实线部分与空白部分的比例:
初始相位
psxy-R-JMtemp-Sc0.15-W1-G255/0/0-O-K-V>>2.ps用sc命令为绘制圆形,(Sa星号,Sf断层等)
(试着分别运行以上4条psxy语句,看看区别)
-N把那些超出边框的数据画上去。
Ps:
Sa星号,Sb条棒,Sc圆,Sd菱形,Se椭圆,Sf断层,Sh六边形,Si倒三角,Sp点,Ss正方形,St三角形,Sv箭头,Sw扇形,Sx十字.
5.psxyfault1.blnfault2.bln-R-JM-Sf0.1i/0.05iR:
1.2-P-O-W0.5p>>2.ps(正断层)
psxyfault.bln-R-JM-Sf0.1i/0.05iRs:
1.2-P-O-W0.5p>>2.ps(走滑断层)
Sf0.1i/0.05iRs:
1.2i0.1i为断层标记出现的间隔,0.05i为标记的大小,
R(L)表示逆冲方向,l(r)s表示走滑。
:
1.2标记的相位。
一条命令可包含多个信息相同的断层数据文件,每个单独的文件生成一条断
层。
6.psvelo-R-JM1.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.echo11538(经纬度)8(字大小)0.0(逆时针旋转角度)1(字体,1黑体,2斜体)2(位置,2表示所写内容的中心出现在前面经纬度上)WUHAN(内容)>text
pstexttext-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.grdgradientin_grdfile-Gout_grdfile-A0-Nt根据你的网格文件生成一个光照的网格文件,使图产生立体感。
-A0(从北面来的光照)-A0/270(北面照向西面)
命令肯定没错,但是我从来没有成功过,暂时认为是机器原因。
12.grdimage-Iintense.grdtopo_asia.grd-JM8i-R70/150/15/60-Ba10f5-K-Cg.cpt>1.ps
pscoast-R-JM-N1-W1/0/0/0-B-Dl-O-K>>1.ps
intense.grd为光照网格文件,topo_asia.grd为地形网格文件,g.cpt为海拔对应的色谱文件。
13.psmeca-R-JMzycmt2_50.txt-Sm0.1i/4u-K-O-P-T0-N-G255/0/0>cmt.ps
Sm对应的数据类型x,y,depth,mrr,mtt,mff,mrt,mrf,mtf,exp,newX,newY,event_title
其他类型数据对应的命令虽然可以画出beachball,但是画不出颜色(郑勇言).缺乏深度数据时加-o,重要。
14.grdraster
用包含30秒×30秒高程数据的dbase文件覆盖原文件后,可以运行grdraster命令。
即调用高程数据,生成gmt所需的grd格式文件。
例子:
grdraster32-R80/100/25/40-Gtibet_S.grd-I0.5m
由于30秒高程数据被分解成几十个小文件存储,如代号为32的文件包含范围-R60/100/-10/40。
在画地形图时往往需要将好几个高程文件拼接起来,或者只需要截取某个文件中的一部分,前者通过多次调用grdraster命令来实现,后者则通过grdraster命令中的-R范围实现,争取做到运行时间最短。
-Goutfile.grd。
-I后接数据的分辨率。
各个文件的范围,代号及分辨率,可以通过在Dos界面下直接运行grdraster查询得到。
大家可以试着运行:
grdraster32-R80/100/25/40-Gtibet_S.grd-I0.5m
grdraster41-R80/100/40/45-Gtibet_N.grd-I0.5m
grdgradienttibet_S.grd-Gtibet_light1.grd-A180/90-Nt
grdgradienttibet_N.grd-Gtibet_light2.grd
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- GMT 小结
![提示](https://static.bdocx.com/images/bang_tan.gif)