Who am I

我的照片
Hefei, Anhui, China
Research Fields: Satellite Geodesy

2012年2月3日星期五

GMT系列一 GMT简介

GMT简介
GMT是什么?
GMT是一款地学绘图软件,由美国哥伦比亚大学地球观测所的两名研究生Paul WesselWalter H. F. Smith1987年创立。后经不断发展与改进,现在该软件已经在地理、地球物理等学科中得到广泛应用。

GMT能做什么?
GMT的主要功能是在各种各样的地图上绘制地震、等高线、地形图等等,当然一般的线性函数或对数、指数函数图都可以绘制。此外该软件还提供了大量的预处理和分析功能,能对数据进行滤波、拟合、内插等工作,不同数据间的格式进行转换,还可进行简单的数学分析甚至谱分析等等。

为什么GMT如此普及?
价格决定一切。GMT是一款完全免费的自由软件,任何用户都可以免费得到它。GMT可以搭载各种操作系统平台上,无论是windows用户还是unix用户都可以安装该软件。此外,高质量的PostScript输出也是众多用户喜欢的原因之一。

GMT的安装(以Ubuntu系统为例)
GMT软件的各个系统安装包都可以在其官方网站上下载,但是如果使用的是Ubuntu系统,新立得软件包里自带了GMT。因此,在Ubuntu系统中安装GMT的步骤如下:
1.      Ubuntu软件中心,下载并安装新立得软件包(synaptic package manager;
2.      sudo apt-get install GMT         安装GMT软件
至此,GMT软件已经安装成功了。但是,如果此时使用GMT的话,所有命令之前都需要加上“GMT”,例如pscoast命令,要改成GMT pscoast,非常的不方便。我们可以通过更改$PATH环境变量来解决。添加环境变量步骤如下:在主目录下~下找到.bashrc文件,vim .bashrc; 添加环境变量,export PATH=”$PATH:/usr/lib/gmt/bin”;用source .bashrc命令运行.bashrc文件。此时,新打开的terminal在使用GMT命令时就无需再加”GMT”了。
注:GMT官方网站上提供的30个例子都可以从新立得软件里找到。
打开新立得软件,搜索GMT-examples,标记——重加载——安装即可。安装完成后,可以在/usr/share/doc/gmt-examples/examples找到所有的例子。

第一个GMT例子
打开terminal终端,让我们开始GMT之旅。
首先,新建一个文本文件,另文件名为plot.job
vim plot.job
然后,在plot.job文件里输入以下代码。这些代码的含义将会在以下内容里详细讲解。
#! /bin/bash
# Purpose: create on the same page 4 maps of North America
ps=example_1.ps
gmtdefaults -D > ./.gmtdefaults4
gmtset ANNOT_FONT_SIZE_PRIMARY 8p \
HEADER_FONT_SIZE 16p            \
 LABEL_FONT_SIZE 10p                \
 LABEL_OFFSET 0.15c
R="-R-140/-50/20/65"

# Plot the 1st map
pscoast ${R} -JM2.75i -Ba20f5/a20f5WSNE -Ggrey -Swhite -Wblack -Dc -K -P > $ps
# Plot the 2nd map
pscoast -R -J -B20f5/20f5WSNE -Glightbrown -Slightblue -Lf-63.5/30/40/1500k+lkm \
    -Di -X3.5i -O -K >> $ps
# Plot the 3rd map
pscoast -R -J -O -K -B -Glightbrown -Slightblue -I1,2/blue -Wdarkblue -N2,- \
    -N1/red -Di -X-3.5i -Y2.75i >> $ps
# Plot the 4th map
pscoast -R -JL-95/42.5/33/45/2.75i -O -B20f5/20f5:."Lambert projection":ESwn -Wdarkblue -Glightbrown -Slightblue -X3.5i -Di -N2,- -N1/red -I1,2/blue >> $ps
第三步:增加文件的“执行”权限,并运行该文件。
chmod +x plot.job           !增加“执行权限”
bash plot.job                    !执行该文件
最后:查看结果(.ps文件即为输出结果)。
通常查看结果又很多方法,这里我们主要介绍三种。
最简单的方法,通常的pdf阅读器就可以查看该文件。
第二种方法:gs example_1.ps
第三种方法:gv example_1.psgv是一个专门查看.ps文件的一个程序,因此在使用该命令时,需要先安装。sudo apt-get install gv。安装完成后,即可通过gv命令查看文件了。

没有评论:

发表评论