科学网

 找回密码
  注册
科学网 标签 GMT

tag 标签: GMT

相关帖子

版块 作者 回复/查看 最后发表

没有相关内容

相关日志

[转载]Linux 下GMT 4.5.9 安装
zhoufcumt 2013-6-7 21:41
转载自: http://seisman.blog.ustc.edu.cn/index.php/archives/515 ,根据自己情况稍微修改了下。 GMT 4.5.9于2013年1月1日正式发布,这个版本依然主要只是bug修订,值得关注的一点在于海岸线数据已经从原来的GSHHS变成了现在的GSHHG,因而安装GMT4.5.9所需要的文件包括以下三个,具体下载以及其他参见原来的文章。 下载安装所需要的三个文件,放在任意文件夹下,比如GMT gmt-4.5.9.tar.bz2 gshhg-gmt-nc3-2.2.2.tar.gz netcdf-3.6.3.tar.gz 注:你可能会看到gshhg-gmt-nc3-2.2.2.tar.gz和gshhg-gmt-nc4-2.2.2.tar.gz两个文件,其分别对应netcdf3和netcdf4,根据自己要安装的netcdf版本确定gshhg版本。 安装netcdf tar -zxvf netcdf-3.6.3.tar.gz //解压netcdf cd netcdf-3.6.3 //进入netcdf文件夹 ./configure --prefix=/opt/netcdf-3.6.3 //configure, prefix指定netcdf的安装路径 make //编译netcdf sudo make install //安装netcdf 安装GMT cd .. //退回到netcdf上级目录 tar -jxvf gmt-4.5.9.tar.bz2 //解压GMT cd GMT4.5.9 //进入GMT文件夹 ./configure --prefix=/opt/GMT-4.5.9 --enable-netcdf=/opt/netcdf-3.6.3/ //配置GMT。prefix为GMT安装路径,--enable-netcdf为前面netcdf的安装路径 make gmt //编译GMT自己的可执行文件及库文件 make suppl //编译其他人提供的GMT扩展包及库文件 sudo make install-gmt //安装GMT自己的可执行文件 sudo make install-suppl //安装GMT扩展包 sudo make install-data //安装GMT数据,其实就是拷贝share目录 sudo make install-man //其实是拷贝man到share sudo make install-doc //其实是拷贝doc到share 安装GMT所需要的海岸线数据 cd .. //返回GMT上级目录 tar -zxvf gshhg-gmt-nc3-2.2.2.tar.gz //解压GMT需要的海岸线数据 sudo cp gshhg-gmt-nc3-2.2.2/* /opt/GMT-4.5.9/share/coast //拷贝数据,将所有文件拷贝到share下的coast目录中 修改环境变量 # for GMT 4.5.9 PATH=/opt/GMT-4.5.9/bin:$PATH MANPATH=/opt/GMT-4.5.9/share/man:$MANPATH LD_LIBRARY_PATH=/opt/GMT-4.5.9/lib:$LD_LIBRARY_PATH
个人分类: GMT|3688 次阅读|0 个评论
GMT第三弹:时间序列的曲线图
热度 2 zhoufcumt 2013-5-31 20:25
本例的脚本和数据文件会在下方提供下载,可根据类似格式和熟读脚本写出你想要的时间序列图,当然matlab可以搞定的,但是既然选择了GMT就要尽量多熟练掌握她的应用吧,GMT == give u surprises in any time, as u know。 NOTE: 如果有更好的建议或者改进请及时通知原作者 zhouforme@gmail.com。 脚本: #!/bin/bash # # Purpose: Plot residuals # Author : Feng Zhou # Date : 2013-05-31 # outps=ats_resi.eps ORI=-P # plot North # plot a line on zero psxy ! -R2001.5/2007.5/-10/10 -JX16.0/7.0 -X2.0 -Y16.7 -K $ORI -W1/200/200/200 -V \ -Ba1.0f0.5/a2.0f1.0sWne $outps 2001.50 0.00 2007.50 0.00 ! # plot ats residual : red line cat $1 | grep -v \* | awk '{print $1, $3, 0.1}' | \ psxy -R -J -K -O $ORI -W2/255/0/0 -V $outps # put text pstext ! -H0 -JX -R -O -K -V $ORI $outps 2001.75 -8 13 0 1 1 N(mm) ! # plot East # plot a line on zero psxy ! -R2001.5/2007.5/-10/10 -JX16.0/7.0 -X0.0 -Y-7.0 -O -K $ORI -W1/200/200/200 -V \ -Ba1.0f0.5/a2.0f1.0sWne $outps 2001.50 0.00 2007.50 0.00 ! # plot ats residual : red line cat $1 | grep -v \* | awk '{print $1, $2, 0.1}' | \ psxy -R -J -K -O $ORI -W2/0/255/0 -V $outps # put text pstext ! -H0 -JX -R -O -K -V $ORI $outps 2001.75 -8 13 0 1 1 E(mm) ! # plot Up # plot a line on zero psxy ! -R2001.5/2007.5/-30/30 -JX16.0/7.0 -X0.0 -Y-7.0 -K -O $ORI -W1/200/200/200 -V \ -Ba1.0f0.5/a5.0f2.5SWne $outps 2001.50 0.00 2007.50 0.00 ! # plot ats residual : red line cat $1 | grep -v \* | awk '{print $1, $7, 0.1}' | \ psxy -R -J -K -O $ORI -W2/0/0/255 -V $outps # put text pstext ! -H0 -JX -R -O -V $ORI $outps 2001.75 -25 13 0 1 1 U(mm) ! # Convert EPS to PNG ps2raster $outps -A -P -Tg # Clean up after ourselves: rm -f .gmt* 效果图: 脚本和数据文件附件: ats_resi.gmt albh_resi.neu
个人分类: GMT|9917 次阅读|5 个评论
咱们的导航定位系统--北斗的轨道图
热度 2 zhoufcumt 2013-5-30 21:42
终于搞定了北斗轨道图,看着感觉不错,GMT的成图效果真好!
个人分类: GMT|4393 次阅读|5 个评论
GMT第二弹:全球站点分布图
热度 2 zhoufcumt 2013-5-29 22:18
全球站点分布图和区域站点分布图类似: NOTE: 如果有更好的建议或者改进请及时通知原作者 zhouforme@gmail.com。 脚本: #!/bin/bash # # Purpose: Plot IGS sites in China # Author : Feng Zhou # Date : 2012-12-10 # ps=igs_global.eps ORI= gmtset ANNOT_FONT_SIZE_PRIMARY 10 PLOT_DEGREE_FORMAT ddd:mm:ssF pscoast -Rg-180/180/-90/90 -JX10i/6i -Ba30f15/a30f15WenS -Dh -A5000 \ -G219/235/199 -S130/199/240 -W0 -K $ORI $ps cat $1 | awk '{print $2, $3}' | psxy -R -J -O -Sc0.15 -Gred -Wthinnest $ps # Clean up after ourselves: rm -f .gmt* 脚本和数据文件下载: igs_global.gmt igs_global.list 效果图:
个人分类: GMT|4324 次阅读|3 个评论
GMT第一弹:区域站点分布图
热度 2 zhoufcumt 2013-5-29 20:40
GMT博大精深,功能强大,我也算是一个刚入门的初学者,希望能将我的学习GMT的经历与大家分享,希望能帮助到没有头绪的初学者,当然对于GMT的高手,我也希望GMTer们能多多交流,共同提高。 NOTE: 如果有更好的建议或者改进请及时通知原作者 zhouforme@gmail.com。 主要针对中国地区IGS站分布图 GMT脚本(igs_cn.gmt): #!/bin/bash # # Purpose: Plot IGS sites in China # Author : Feng Zhou # Date : 2012-12-10 # ps=igs_cn.eps gmtset HEADER_FONT_SIZE 10 X_ORIGIN 4c Y_ORIGIN 8c \ ANNOT_FONT_SIZE_PRIMARY 6 FRAME_WIDTH 0.05c pscoast -R70/140/15/55 -JM5i -Glightyellow -Slightblue -UMade by Feng Zhou -A1000 \ -Ba10f5/5:.IGS Stations in China: -K -Dh -N1/1p -Wthinnest -P $ps cat $1 | awk '{print $1, $2}' | psxy -R -J -O -Sc0.15 -Gred -Wthinnest $ps # Clean up after ourselves: rm -f .gmt* NOTE:$1是可执行脚本后面带的参数,比如本例子./igs_cn.gmt igs_cn.list这样才可以正确运行。 数据文件(igs_cn.list): 115.89 39.61 BJFS LM 121.20 31.10 SHAO LM 114.36 30.53 WUHN LM 87.60 43.81 URUM LM 125.44 43.79 CHAN LM 109.22 34.37 XIAN RM 102.80 25.03 KUNM LM 91.10 29.66 LHAZ LM 87.18 43.47 GUAO RM NOTE:格式为经度,纬度(北纬为正,南纬为负),站点名,站点名与本例子的圆圈表示的相对位置 LM-left middle 效果图如下: 上面的例子画出来的没有站点名(针对很多个站点的时候比较好),如果站点数不多,想加上站点名,数据文件不变,脚本需要改变一下。 改后脚本: #!/bin/bash # # Purpose: Plot IGS sites in China # Author : Feng Zhou # Date : 2012-12-10 # ps=igsname_cn.eps gmtset HEADER_FONT_SIZE 10 X_ORIGIN 4c Y_ORIGIN 8c \ ANNOT_FONT_SIZE_PRIMARY 6 FRAME_WIDTH 0.05c pscoast -R70/140/15/55 -JM5i -Glightyellow -Slightblue -UMade by Feng Zhou -A1000 \ -Ba10f5/5:.IGS Stations in China: -K -Dh -N1/1p -Wthinnest -P $ps cat $1 | awk '{print $1, $2}' | psxy -R -J -O -K -Sc0.15 -Gred -Wthinnest $ps cat $1 | awk '{print $1, $2, 5, 0, 5, $4, $3}' | pstext -R -J -O -Dj0.15/0 -Gred -N $ps # Clean up after ourselves: rm -f .gmt* 效果图: 注明:由于上述脚本出现了问题,因此上传脚本和数据文件,供大家下载。 igs_cn.list igsname_cn.gmt igsnoname_cn.gmt
个人分类: GMT|6530 次阅读|2 个评论
GMT的中文支持
chchao 2013-5-21 11:18
一直以为GMT中不能使用中文,只能用AI等软件进行后期加工。 受下面两篇博文的启发,成功设置GMT中的中文支持,没想到经过简单的设置后也可以直接使用中文。 博文地址: http://xxqhome.blog.163.com/blog/static/1967330202011112810120598/ http://hi.baidu.com/guyueshuiming/item/0052df53852ee4494fff20c3 以下设置在Ubuntu下试验成功,具体方法如下: 1、将字体文件从Windows下拷到Linux下的目录中(如/usr/share/fonts/winFonts) 2、处用gs命令和mkcidfm.ps形成cidfmap 如 sudo gs -q -dBATCH -sFONTDIR=/usr/share/fonts/winFonts -sCIDFMAP=/usr/share/ghostscript/9.04/Resource/Init/FAPIcidfmap /usr/share/ghostscript/9.04/lib/mkcidfm.ps 注意这里的FAPIcidfmap随不同的ghostscript版本,位置有可能不同,文件名也可能会不同(如可能为cidfmap) 3、根据FAPIcidfmap文件中的信息,修改GMT中的字体配置文件 PS_font_info.d(一般在$GMT_SHARE/pslib下),在后面添加 NSimSun--GB-EUC-H 0.700 1 SimHei--GB-EUC-H 0.700 1 KaiTi--GB-EUC-H 0.700 1 LiSu--GB-EUC-H 0.700 1 FangSong--GB-EUC-H 0.700 1 YouYuan--GB-EUC-H 0.700 1 MicrosoftYaHei--GB-EUC-H 0.700 1 MicrosoftYaHeiBold--GB-EUC-H 0.700 1 注意:一定要在字体名称后面加上 --GB-EUC-H ,否则不能正常显示 4、用pstext -L查看字体对应的数字编码,如 35 STSong-Light--GB-EUC-H 36 STFangsong-Light--GB-EUC-H 37 STHeiti-Regular--GB-EUC-H 38 STKaiti-Regular--GB-EUC-H 5、在GMT脚本中,需插入中文的地方用如下命令 @%font_no% 中文 @%% 如需使用楷体,根据上面的字体编码,楷体对应字体编码为38,则形式为 @%38% 我是中文 @%% 注意:文件的编码方式要保存为 gb18030,utf-8编码显示不正常,不知道其它系统上是否是这样?
8296 次阅读|0 个评论
在RedHat Linux下安装GMT
cambaluc 2013-5-10 12:22
再安装了一次GMT。因linux用的少,有的配置还是模糊,所以还是记录一下吧,万一有再装的可参考。 先装NetCDF 先用root用户登录,mkdir /usr/local/netcdf,tar xvfz netcdf-4.4.1.1.tar.gz,cd netcdf-4.2.1.1, ./configure --prefix=/usr/local/netcdf --disable-netcdf-4,make check install,make,make install 再用自己的用户名登录,改文件./home/kj/.bash_profile,(可用ls -a|grep bash找该文件), 加PATH=$PATH:/usr/local/netcdf/bin export NETCDF=/usr/local/netcdf 装GMT,先用root用户 tar -jxvf gmt-4.5.8.tar.bz2,tar -jxvf gmt-4.5.8-non-gpl.tar.bz2,mkdir /usr/local/GMT ./configure --prefix=/usr/local/GMT --enable-netcdf=/usr/local/netcdf, make,make install,cp -r share/ /usr/local/GMT 再用自己的用户, export GMTHOME=/usr/local/GMT, export PATH=/usr/local/GMT/bin:$PATH,source ~/.bashrc 再解压三个海岸线文件 tar -jxvf GSHHS2.2.0_high.tar.bz2及_full,_high,拷到/usr/local/GMT/share/coast/(可用root权限) 试一下:pscoast -R0/360/-90/90 JG280/30/6i -Bg30/g15 -Dc -A5000 -G255/255/255 -S150/50/150 -Pmap.ps 命令来自《能用地学制图工具GMT入门百度文库》 改改参数可绘出下图:
4246 次阅读|0 个评论
Graphic Mapping Tool(GMT)
Master123 2013-4-24 10:06
http://gmt.soest.hawaii.edu/
个人分类: GMT|2340 次阅读|0 个评论
GMT中的project命令
热度 1 zhanggw 2013-2-21 10:26
project命令还是很好用的,比如可以将一个矩形区域内的地震投影到剖面上 project tt.reloc -C${center} -A${azimuth} -Fxyzpqrs -Q -W${width} -L${length} -V ! tang.tmp 只需要设定好中心点,宽度及长度,方位角,-F选项也要注意,设置参数不同,输出的数据格式也不同 任意角度的剖面图投影,就很容易了。
个人分类: SOFT|5679 次阅读|1 个评论
GRDCONTOUR and GRDMATH (GMT)-Plot event map
热度 1 shenxzh 2013-1-27 13:09
gmtset FRAME_PEN 1.0p gmtset FRAME_WIDTH 0.1c gmtset HEADER_FONT 5 gmtset LABEL_FONT 5 gmtset HEADER_FONT_SIZE 10 gmtset LABEL_FONT_SIZE 12 gmtset ANOT_FONT_SIZE 10 gmtset BASEMAP_TYPE plain gmtset TICK_LENGTH -0.1c gmtset TICK_PEN 1.0p gmtset LABEL_OFFSET 0.1c gmtset MAP_TITLE_OFFSET 5p set ps=yuanzhen.ps set scale=2.4e8 set eq=eq.txt set lon=103.6 set lat=36.0 rem grdmath -R-180/180/2/4 -I6/0.1 X 4 MUL PI MUL 180 DIV COS Y 2 POW MUL = dist.nc grdmath -Rg -I1 %lon% %lat% SDIST = dist.nc grdmath -Rg -I1 %lon% %lat% SBAZ = baz.nc rem grdmath -Rg -I1 %lon% %lat% SDIST = dist.nc rem pscoast -R-180/180/-90/90 -Je103.6/36/1:%scale% -B30g30 -Ggray -K -A2500 -Dh -X3.5i -Y1.25i %ps% pscoast -R0/360/-90/90 -Je103.6/36/1:%scale% -I2/0.25p,blue -N1/0.25p -W0.25p,white -Ggreen -Sblue -K -A2500 -Dh -X3.0i -Y1.1i %ps% rem pscoast -R-180/180/-90/90 -JA103.6/36/6i -I2/0.25p,blue -N1/0.25p -W0.25p,white -Ggreen -Sblue -K -A2500 -Dh -X3.0i -Y1.1i %ps% rem pscoast -R-180/180/2/4 -JP6i -I2/0.25p,blue -N1/0.25p -W0.25p,white -Ggreen -Sblue -K -A2500 -Dh -X3.5i -Y1.25i %ps% grdcontour dist.nc -A60 -C30 -S8 -K -O -R -J -Wathick,black -Wcthin,black %ps% rem grdcontour dist.nc -B30Ns -P -C2 -S4 --PLOT_DEGREE_FORMAT=+ddd %ps% grdcontour baz.nc -C45 -S8 -K -O -R -J -Wcthicker,white,- %ps% rem grdcontour dist.nc -A15+u+fblack -S8 -O -K -J -Wathin,black %ps% gawk "{print $1,$2,$3*$3*$3/1000}" %eq% |psxy -: -R -J -O -K -Sc -Gred -Wthinnest %ps% rem psxy -: %eq% -R -J -O -K -Sc0.08 -Gblack -Wthinnest %ps% echo %lon% %lat% | psxy -R -J -O -K -St0.7 -Gblack -Wthin %ps% pause del dist.nc .gmt*
个人分类: Resear note|6936 次阅读|2 个评论
GAMIT 中 error: GMT Warning: Old GMT 3 .gmtdefaults 解决方案
zhangyong1361 2013-1-4 16:26
GAMIT 中 error: GMT Warning: Old GMT 3 .gmtdefaults 解决方案
error: GMT Warning: Old GMT 3 .gmtdefaults file found. May not be fully compatible with GMT 4. GMT Warning: It is recommended that you migrate your GMT 3 settings to GMT 4 settings. 环境说明: 系统:Ubuntu 10.04 LTS GAMIT/GLOBK版本: 10.4(更新到20120824) GMT版本:4.5.1(通过 apt-get install gmt 安装的) 出现的问题文中开头所述,问题原因是存在旧的配置文件,解决思路指定正确的配置文件。 具体解决方法: 1、在GAMIT解算完后生成的DOY目录下,按 “Ctrl+H”显示隐藏文件,找到“.gmtdefaults4”文件。 2、在任何地方启动终端,输入"gmtdefaults -D ~/.gmtdefaults4" 意思是指定GMT默认配置文件为主文件夹下的.gmtdefaults4文件。 3、将第一步中找到的文件复制到主文件夹下,覆盖原文件。 4、打开主文件夹的.gmtdefaults4文件,设置里面的PAPER_MEDIA= a4,保存即可。 配置完后就不会出现之前的错误了,并且绘图成功且美观。如下图: 坐标重复性图: 基线重复性图: wrnm图: 对流层延迟图(总的 Fiexd): 如有问题可以联系交流,E-mail:zhangyong1361@163.com 张勇 2013-01-04
个人分类: GAMIT/GLOBK|5578 次阅读|0 个评论
[转载]GMT绘制震源球:psmeca
热度 1 cgh 2012-12-19 08:43
GMT绘制震源球:psmeca Plot seismological symbols on maps 用法:psmeca infiles -Jparams -Rwest/east/south/north -Sformatscale ] ] ] ] ] x_shift ] x_shift ] ]] 下面对一些特别的选项做下介绍: -C 在绘制震源球时,GMT会首先在震源处(经纬度由infiles给定)绘制一个小圆,然后用直线连接震源与震源球, pen指定线属性,包括宽度、颜色、线型等,pointsize指定小圆的大小。 -D只绘制某一深度范围之内的震源球 -E、-G 指定了震源球压缩、拉伸部分的颜色,默认压缩部分为黑色,拉伸部分为白色 -L 绘制震源球轮廓 -M所有震级的震源球大小相同,大小由-S给出 -N不要跳过那些震源球落在图形边界外的 -r 在文本周围绘制框 -S指定输入文件格式以及震源球大小 (c) Focal mechanisms in Harvard CMT convention X, Y, depth, strike1, dip1, rake1, strike2, dip2, rake2, moment, newX, newY, event_title (a) Focal mechanism in Aki Richard's convention: X, Y, depth, strike, dip, rake, mag, newX, newY, event_title (p) Focal mechanism defined with X, Y, depth, strike1, dip1, strike2, fault, mag, newX, newY, event_title fault = -1/+1 for a normal/inverse fault (m) Sesmic moment tensor (Harvard CMT, with zero trace) X, Y, depth, mrr, mtt, mff, mrt, mrf, mtf, exp, newX, newY, event_title (z) Anisotropic part of seismic moment tensor (Harvard CMT, with zero trace) X, Y, depth, mrr, mtt, mff, mrt, mrf, mtf, exp, event_title (d) Best double couple defined from seismic moment tensor (Harvard CMT, with zero trace) X, Y, depth, mrr, mtt, mff, mrt, mrf, mtf, exp, newX, newY, event_title (x) Principal axis X, Y, depth, T_value, T_azim, T_plunge, N_value, N_azim, N_plunge P_value, P_azim, P_plunge, exp, newX, newY, event_title (t) Zero trace moment tensor defined from principal axis X, Y, depth, T_value, T_azim, T_plunge, N_value, N_azim, N_plunge P_value, P_azim, P_plunge, exp, newX, newY, event_title (y) Best double couple defined from principal axis X, Y, depth, T_value, T_azim, T_plunge, N_value, N_azim, N_plunge P_value, P_azim, P_plunge, exp, newX, newY, event_title -Tn plane n=1: 画第一个节点平面 n=2:画第二个节点平面 n=3:两个都画 ------------------------------------------------------------------------------------------- 实例: #!/bin/bash #This script is used to plot a focal mechanisms on maps PS=meca.ps LON=80 LAT=35 psbasemap -Rg -JE$LON/$LAT/180/6i -B:".GCMT": -K -P ${PS} pscoast -Rg -J -B -Dc -A10000 -Wthin -G200 -K -P -O ${PS} #plot mechanism psmeca -Rg -JE$LON/$LAT/180/6i -Sm1 -T0 -W1p -K -O -P EOF ${PS} -26.55 -60.71 24 1.06 -0.24 -0.82 0.09 0.64 -0.62 26 -26.55 -60.71 200904161457A EOF 效果图 -------------------------------------------------- 震源机制解文件可以以多种方式输入,不做震源所以很多参数都不知道什么意思,自己写震源机制解文件有困难。GlobalCMT可以将catlog导出为GMTpsmeca所需要的格式,对应-Sz选项,可参考之。 参见 http://www.globalcmt.org/CMTsearch.html http://blog.sina.com.cn/s/blog_872ecf9f0100vysp.html
个人分类: Linux|8112 次阅读|1 个评论
[转载]画震源机制剖面的GMT工具:pscoupe
热度 2 cgh 2012-12-19 08:37
画震源机制剖面的GMT工具:pscoupe画震源机制剖面的GMT工具,也可以画其他类似东西的剖面,比如地震深度剖面。 语法: pscoupe files −J parameters −R west / east / south / north −A parameters parameters ] ] ] ] ] ] ] ] ] ] ]] ] ] GMT的震源机制工具,比如psmeca, pscoupe一般默认是绘下半球投影。还要注意的是在利用pscoupe绘震源机制剖面的时候投影面是不同于平面上的。 -J 控制投影方式 -R 坐标范围 -A 选择剖面 -Aa lon1/lat1/lon2/lat2/dip/p_width/dmin/dmax lon1,lat1,lon2,lat2分别是两个端点的经纬度 dip是投影平面的倾角 p_width是投影剖面两侧的宽度 dmin和dmax是距离水平面最小和最大的距离 加上f可以获得框格 −Ab lon1/lat1/strike/p_length/dip/p_width/dmin/dmax −Ac x1/y1/x2/y2/dip/p_width/dmin/dmax 类似于-Aa,但用于笛卡尔坐标系下 −Ad x1/y1/strike/p_length/dip/p_width/dmin/dmax 类似于-Ab,但用于笛卡尔坐标系下 -S 选择数据文件 −Sa scale ]] 用于AKI的震源机制约定,第1~10列分别为经度、纬度、深度、走向、倾角、滑动角、震级、 0、 0、文字说明 −Sc scale 用于哈佛的震源机制约定,第1~14列分别为经度、纬度、深度、节面1的走向、倾角、滑动角、节面2的走向、倾角、滑动角、地震矩的小数和指数部分、0、0、文字说明 −Sp scale ]] 当只有部分震源机制解参数的时候,可用这个。第1~11列分别为经度、纬度、深度、节面1的走向、倾角、节面2走向、+1(逆断层)/-1(正断层)、震级、0、0、文字说明 −Sm|d|z scale ]] 用于绘地震矩张量解。第1~13列分别为经度、纬度、深度、mrr, mtt, mff, mrt, mrf, mtf in 10*exponent dynes-cm、指数部分、0、0、文字说明 −Sx scale ]] 用于已知力轴参数来绘震源机制解。第1~16列分别为经度、纬度、深度、value (in 10*exponent dynes-cm), azimuth, plunge of the T, N, and P axes、指数部分、0、0、文字说明 −s symbol ]] 用其他图形来替代震源机制球。第1~4列分别为经度、纬度、深度、说明。 -B 设置边框刻度 −E fill 膨胀象限的颜色设定 −G fill 压缩象限的颜色 −L 画笔属性 −M 忽略震级大小 −T 控制所绘节面的数量,0-两个节面都绘,1-节面1,2-节面2,默认是0 −a ]],计算并绘P、T轴的位置,符号类型有:( c ) circle, ( d ) diamond, ( h ) hexagon, ( i ) inverse triangle, ( p )point, ( s ) square, ( t ) triangle, ( x ) cross. −e fill 填充T轴 −g fill 填充P轴 http://hi.baidu.com/purplepool/item/4d2a961493945e4ee75e0617
个人分类: Linux|10026 次阅读|10 个评论
GMT 5.0.0b的安装方法
热度 1 Master123 2012-10-20 21:05
在网上能够搜索到的一些同行的博文: (1)GMT5.0在Linux环境下的安装 http://blog.163.com/jey_df/blog/static/18255016120111111104913536/ (2) win7下安装Cygwin/GMT5.0成功 http://xxqhome.blog.163.com/blog/static/1967330202011101573237365/ (3)GMT(the Generic Mapping tools)下载最简单的安装 http://bbs.sciencenet.cn/forum.php?mod=viewthreadtid=535676page=1 这些方法感觉要么程序繁琐,要么文件太多。现在GMT 5.0.0b已经提供了windows版本的下载,为什么不直接下载安装呢?(也许在Linux或其他版本下另有好处,我还不知道) 点击链接 http://www.soest.hawaii.edu/gmt5/ 点击左侧"DOWNLOAD",在右侧高亮显示的字体中选择“ GMT Windows ”,网页下方找到gmt-5.0.0b_install32.exe下载,然后安装后重启。可以在命令行里输入pscoast试一下。
个人分类: GMT|6151 次阅读|1 个评论
简单的GMT安装方式
热度 1 peluo 2012-7-27 14:25
采用sudo apt-get install gmt安装,系统会自动安装其需要的NetCDF库等,默认将NetCDF库安装在usr/lib中,gmt安装在usr/lib/gmt路径下,因此安装完成后可以在.bashrc中这样设置环境变量: export NETCDFHOME="/usr/lib" export GMTHOME="/usr/lib/gmt" export PATH="$PATH:$GMTHOME/bin"
个人分类: Linux|5005 次阅读|1 个评论
GMT 在地图上绘制线条
热度 1 xinxiang 2011-5-18 17:11
曾经因为画断层图而郁闷了很久,最后发现gmt可以干这事。 使用psxy画点,同时将相邻的两点连接成线即可,关键在于如何把很多不同的线条区分开,具体命令如下: psxy fault.dat -JM -R -O -K -Wred -M'370' test.ps 命令解释:fault.dat 包含要画的数据点的坐标文件。其中隐含着绘图的关键:’370‘。每条线之间使用370隔开,如: -115.4 35.3 第一条线的终点坐标 370 -118.4 40.4 第二条线的起始点坐标 如果中间没有370隔开,那么所有的点都会连接在一起。 当然,370是个人定义的标志,也可以换成其他的。 下面是自己画的加利福尼亚地区的断层线图,断层线条和州的边界线条没有区分,都是黑色的。
个人分类: GMT|25126 次阅读|2 个评论
[转载]Windows XP下安装GMT
热度 1 zheyang 2011-3-23 15:53
GMT − The Generic Mapping Tools data processing and display software package GMT is a collection of public-domain Unix tools that allows you to manipulate x,y and x,y,z data sets (filtering, trend fitting, gridding, projecting, etc.) and produce PostScript illustrations ranging from simple x-y plots, via contour maps, to artificially illuminated surfaces and 3-D perspective views in black/white or full color. Linear, log10, and power scaling is supported in addition to 25 common map projections. The processing and display routines within GMT are completely general and will handle any (x,y) or (x,y,z) data as input.Significantly improved versions (3.1-3.3, 3.3.1-6), 3.4, 3.4.1-5, and 4.0-4.5.3 were released between November 1998 and July 2010, culminating in the Mar 2011 introduction of 4.5.6. GMT now is used by 15,000 users worldwide in a broad range of disciplines. The homepage of GMT is http://www.soest.hawaii.edu/gmt/ The installation for windowsis as follows: Reprint from the website 1.准备文件如下 (文件下载网址:http://gmt.soest.hawaii.edu/) netcdf-4.1.1 GMT_basic_install.exe GMT_pdf_install.exe GSHHS_highfull_install.exe 2.执行安装以下文件 (以安装到D盘为例) GMT_basic_install.exe GMT_pdf_install.exe GSHHS_highfull_install.exe 此时会在D:\GMT目录下产生有以下文件夹:bin include lib share 3.把netcdf-4.1.1文件拷到D:\GMT目录下 (解压netcdf 后的文件) 4.把D:\GMT\share下所有 *.d; *.cdf; *.ps 文件及conf文件夹下的所有文件拷到D:\GMT目录下 5.在D:\GMT目录下创建gmt.bat文件 目的:设置GMTHOME,HOME 和 NETCDF,修改PATH,LIB,INCLUDE下面是gmt.bat文件的内容,拷贝下来再,保存为gmt.bat, 然后根据自己的情况加以修改,这里以上面创建的目录为例。 gmt.bat文件的内容如下: ECHO OFF REM ---------------------------------------------------- REM REM @(#)gmtenv.bat 1.10 03/08/00 REM REM Copyright (c) 1991-2000 by P. Wessel and W. H. F. Smith REM See COPYING file for copying and redistribution conditions. REM REM This program is free software; you can redistribute it and/or modify REM it under the terms of the GNU General Public License as published by REM the Free Software Foundation; version 2 of the License. REM REM This program is distributed in the hope that it will be useful, REM but WITHOUT ANY WARRANTY; without even the implied warranty of REM MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the REM GNU General Public License for more details. REM REM Contact info: www.soest.hawaii.edu/gmt REM -------------------------------------------------------------------- REM REM Sets the GMT environment needed under WinXP. REM REM Author: Paul Wessel, 07-MAR-2000 REM REM ---------------------------------------------------- REM Set NETCDF, GMTHOME, and HOME: REM ---------------------------------------------------- SET NETCDF=d:\GMT\netcdf-4.1.1 SET GMTHOME=d:\GMT SET HOME=d:\GMT REM ---------------------------------------------------- REM Must set INCLUDE and LIB if GMT source is to be used REM OR compiled. If not, REM these lines out. REM ---------------------------------------------------- SET INCLUDE=%INCLUDE%;%NETCDF%\INCLUDE SET LIB=%LIB%;%NETCDF%\LIB;%GMTHOME%\LIB REM ---------------------------------------------------- REM STOP HERE, Now appending to PATH: REM ---------------------------------------------------- SET PATH=%PATH%;%GMTHOME%\BIN;%NETCDF%\LIB REM ---------------------------------------------------- ECHO GMT 4.5.2 Environment initialized doskey 6.进入DOS,在DOS下进入D:\GMT执行gmt.bat,显示如下: D:\GMTgmt.bat D:\GMTECHO OFF GMT 4.5.2 Environment initialized 此时可输入一个命令psxy等进行测试,如果出现此命令的说明,则表示成功了。 7.此外gawk也是window下常用程序,可以在下面网站下载 http://gnuwin32.sourceforge.net/packages/gawk.htm 下载后,安装在D:\GMT 目录下(可以安装在任何目录下),只要把D:\GMT\GnuWin32\bin下面的所有文件拷贝到 D:\GMT\bin 下就可以了。 8.PS查看器ghostview的安装 在安装Ghostview前先安装Ghostscript,这两个软件都可以网上自由获取 http://pages.cs.wisc.edu/~ghost/ 按先后顺序依次安装,可以安装在任何一个文件夹。安装后把gs\gs8.63\bin 下的gswin32.exe可执行文件和Ghostgum\gsview下的gsview32.exe可执行文件拷贝到D:\GMT\bin 下即可。 近来初学GMT软件,在Win下先尝试安装了该软件,安装过程中借鉴了很多网上分享的资料,对一些作者表示感谢。
个人分类: 转载|2681 次阅读|0 个评论
Lectures on Varifolds
zjzhang 2009-9-28 22:28
Chapter 4 Theory of Rectifiable n-Varifolds 15 Basic Definitions and Properties 16 First Variation 17 Monotonicity Formula And Basic Consequences 18 Pioncare and Sobolev Inequalities 19 Miscellaneous Additional Consequences Of the Mo This is done when I was suffering lecturing Geometric Measure Theory by Leon Simon. Chapter 5 is devoted to Allard regularity theory, but I've no time to input it in a Latex form.
个人分类: 数学|3624 次阅读|0 个评论
GMT 及Win4GMT简介及其安装技巧
热度 1 kiger 2009-9-19 17:17
GMT-The Generic Mapping Tools ,即通用制图工具,是一个开源的工具集,主要用于 XY 和 XYZ 数据的处理和显示,包括栅格化、过滤和各种地图投影等操作。 GMT 最初是由哥伦比亚大学的两位毕业生 Paul Wessel 和 Walter H. F. Smith 于 1988 年创建的,至今已经发展到版本 4.5.0 。 经过这段时间的了解,感觉该软件主要有以下几个优点:( 1 )完全免费;( 2 )源代码是开放的,可以进行二次开发;( 3 )可以在多种系统下安装和使用,如 UNIX , Linux , MkLinux , Mac OS 9 , Windows 和 OS/2 ,但是每个系统下的安装要求是不一样的;(4) GMT 好比一个工具箱,每个工具之间相对独立,操作起来更为灵活,这也是 GMT 研究小组的一个特别考虑;( 5 )有很多免费数据资源都可以被 GMT 调用,尤其是 NOAA , USGS , NGDC , NEIC 等;( 6 )制作出来的图形图像精美,感觉比 ArcGis 和 MapGis 等软件成图的效果要好。正是由于以上诸多优点,该软件已经被较为广泛应用于医学,工程,物理,数学,地球科学,社会科学等领域。 由于该软件的最初是在 UNIX 系统上发展的,所以对于初学者来说早期的版本安装和使用方面还是有点麻烦的,我也是学习第一次安装该软件,经过摸索,终于在 Windows XP 系统下安装成功。习惯了 Windows 操作的人来说,突然要用 DOS 操作,确实需要一段时间的适应,故可以安装一个 GUI (图形交互界面)软件 Win4GMT 。具体安装过程如下: 1. Install GMT 3 (E g : C : \GMT\ ) . As you will see GMT3 programs will be installed in C : \GMT\ bin and rivers, border, shorelines, patterns etc in C:\GMT\share\coast and gmtenv.bat in C : \GMT\ ( Note: I n GMT 4 .5.0 , gmtenv.bat cant be need any more, and the installer has done the work and set the path. ) 2. S et Home: Modify gmtenv.bat with an ASCII editor in order to set your home dir . SET NETCDF= C : \GMT\ SET GMTHOME= C : \GMT\ SET HOME= C : \ USERS \ USERS \ PROJECTS or C :\ GMT \ TEMP SET PATH=%PATH%; %GMTHOME%\BIN; %NETCDF%\LIB 3. R un cmd command to open this window i n a DOS window . I n GMT 4 .5.0 , .gmtdefaults is created by default in C : \ GMT \bin, you should copy it to C : \GMT\. Then c lose the window . 4. if necessary, install Ghostsview or other postscript viewer or editor . 5. Put netcdf.dll in C : \GMT\. 6. Unzip the Win4GMT and r un setup.exe to install it . To avoid problems , Please install Win4GMT in the GMT Home Directory ( C : \GMT ). 7. Reboot the computer . 8. Configure folders in Preferences menu when starts the program or edit win4gmt.ini in C:\WINDOWS or C:\WINNT folder . Notes: Setup will create a subdirectory ( man ) where will be installed the html version of the GMT manuals. If you have problems with OLEAUT32.DLL (in some cases, problems has been reported) . R ename your oleaut32.dll (just in case) and copy the ftp file in your system directory ( C : \WINDOWS\S ystem) . Then r estart the computer. 其它: GMT下载地址: http://gmt.soest.hawaii.edu/ Win4GMT下载地址: http://www.icm.csic.es/gma
个人分类: 学研心得|11945 次阅读|3 个评论
GMT中的 surface
wangchangzai 2009-8-12 15:31
gmt 中的surface---不完整 2008-05-21 23:34 在gmt中由xyz文件转到grid的文件有两个主要命令:xyz2grd, surface,前面已经提到xyz2grd,xyz2grd要求等间隔采样,有时候显得很不方便。而surface可以提供一个较好的方式 Adjustable tension continuous curvature surface gridding 可调整的曲线表面网格?由manual给出的解释是 A continuous curvature gridding algorithm 应该可以理解为曲率连续的网格插值算法,具体的例子在cookbook的6.14,6.16都有用到 主要选项: surface -Goutput_grdfile_name -Ixinc ] -Rwest/east/south/north ] ] ] ] ] ] ] -G 输出文件名 -xinc/yinc: 给出插值后的网格间隔 -R:指定坐标范围 -A: -Z:是插值的收敛标准 -L:?强行指定解的范围? without any modifiers indicate that x is longitude and periodic in 360: -L constrain the range of output values: -Lllimit specifies lower limit; forces solution to be = limit. -Lulimit specifies upper limit; forces solution to be = limit. limit can be any number, or the letter d for min (or max) input data value, or the filename of a grdfile with bounding values. . Example: -Ll0 gives a non-negative solution. -N sets max n_iterations in each cycle; default = 250. 设置每个圆的最大迭代数 -S 除非网格是病态的 sets search_radius to initialize grid; default = 0 will skip this step. This step is slow and not needed unless grid dimensions are pathological; i.e., have few or no common factors. Append m to give search_radius in minutes. -T adds Tension to the gridding equation; use a value between 0 and 1. default = 0 gives minimum curvature (smoothest; bicubic) solution. 1 gives a harmonic spline solution (local max/min occur only at data points). typically 0.25 or more is good for potential field (smooth) data; 0.75 or so for topography. Experiment. Append B or b to set tension in boundary conditions only; Append I or i to set tension in interior equations only; 关系到插值的结果?详情见cookboo 6.16图 -Z sets over_relaxation parameter. Default = 1.4 Use a value between 1 and 2. Larger number accelerates convergence but can be unstable. Use 1 if you want to be sure to have (slow) stable convergence. 控制收敛速度
个人分类: 生活点滴|14488 次阅读|0 个评论

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-17 12:13

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部