GMT笔记
GMT笔记（周汉强） 出处 http://hcc.earth.sinica.edu.tw/note/gmt/gmt.htm
GMT的网站：
GMT笔记目录检索
Windows版GMT的简易安装使用说明
psbasemap GMT画底图的指令
pscoast GMT画陆块海洋和海岸线的程序
psxy GMT在地图上画点线的程序
grdimage GMT画二维等高图的指令
psscale GMT画比例尺的指令
grdcontour GMT画二维等高线图的指令
grd数据的处理
一些例子
psbasemap GMT画底图的指令
一定要有的参数
-B 设定图边界尺标的间隔
a 文字批注的间隔
f 边框标记的间隔
g 隔线的间隔
如果要针对X轴或Y轴做不同的设定，就在参数间加斜线/
NEWS 每个轴加上标记
news 每个轴不加标记
:.title: 在图上面加title这个标题
:title: 在轴上加title这个标题（只有在-Jx的情况下有用）
例如-BNEWsf2a3g2:xx:/f3a3g2:yy::.title:
-Jx(-JX) 线性投影(Linear Projection)一般绘图用
-Jx后面接的数字代表每单位的长度（预设长度单位为英吋，加c的话就变公分）
-JX后面接的数字代表轴的全长
如果要针对X轴或Y轴做不同的设定，就在数字间加斜线/
-Jm(-JM) 麦卡托投影(Mercator Projection)地图绘图用
-Jm后面接的数字代表每度的长度（预设长度单位为英吋，加c的话就变公分）
-JM后面接的数字代表整张图的宽
如果要针对X轴或Y轴做不同的设定，就在数字间加斜线/
-R 地图的经纬度范围
西/东/南/北 依序填入经纬度
>out.ps 把画出来的图输出到out.ps这个档案
不一定要有的参数
-E 旋转地图的视野方向
-G 把图着色，用r/g/b调色
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
psbasemap -JM7 -R117/126/18/26 -BNEWsf2a3g2:xx:/f3a3g2:y::.test: -P -V -K >123.ps
回到目录：回到最前面
pscoast GMT画陆块海洋和海岸线的程序
一定要有的参数
-J 投影法，和前面的一样
-R 地图的经纬度范围，和前面的一样
不一定要有的参数
-E 旋转地图的视野方向
-G 把干地着色，用r/g/b调色
-S 把湿地着色，用r/g/b调色
-D 图的分辨率(f,h,i,l,c)左往右为分辨率高至低
-I 加河流上去(1-10)大河至小河(a,r,i,c)大小河都画
-W 画海岸线
-N 画边界(1.国界,2.州界,3.海洋,a.全部)
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
pscoast -JM7 -R117/126/18/26 -G50/50/50 -S150/150/150 -W1 -Di -V -O -K >>123.ps
回到目录：回到最前面
psxy GMT在地图上画点线的程序
一定要有的参数
-J 投影法，和前面的一样
-R 地图的经纬度范围，和前面的一样
不一定要有的参数
-A 画大圆线
-B 和先前psbasemap的一样
-C 对照颜色，后面要接文件名称
-G 把符号着色，用r/g/b调色
-K 在这一行下面还有后续的指令
-N 不要跳过超出范围的点
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-M 线分段
-Sa 在图上画星星，后面要接大小
-Sc 在图上画圆圈，后面要接大小
-Sd 在图上画钻石，后面要接大小
-Se 在图上画椭圆形，在资料的第3,4,5栏要有方向长轴短轴
-Sh 在图上画六角形，后面要接大小
-Ss 在图上画正方型，后面要接大小
-St 在图上画三角形，后面要接大小
-Si 在图上画反三角，后面要接大小
-Sv 在图上画向量，在数据的第3,4栏要有方向长度
-Sx 在图上画叉叉，后面要接大小
-V 画面上会显示本行指令的执行结果
-W 笔宽
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式一 把数据直接打在档案里
psxy -JM7 -R117/126/18/26 -W1 -O -K -V‹‹! >>123.ps
120 23
124 27
!
表示方式二 把数据另外打在数据文件里
psxy 123.dat -JM7 -R117/126/18/26 -W1 -M -V -O -K >>123.ps
在123.dat的档案里是这样（这是画线条用的档案）
>
120 25
123 25
>
121 23
123 27
在123.dat的档案里是这样（这是点用的档案）
120 25
123 25
121 23
123 27
表示方式三 把更多的数据打在数据文件里
psxy 123.dat -R117/126/18/26 -JM7 -Sc -Ccpt >> 123.ps
在123.dat的档案里要有四栏数据，依序为X,Y,颜色代码,圆的大小
回到目录：回到最前面
grdimage GMT画二维等高图的指令
一定要有的参数
grd grd档案的档名
-C 对照颜色，后面要接文件名称
-J 投影方式，和前面一样
不一定要有的参数
-B 和先前psbasemap的一样
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-R 地图的经纬度范围，和前面的一样
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
grdimage file.grd -JOc122/21/22/69/4.8i -R117/126/18/26 -Cfile.cpt -P -K -V >> 123.ps
回到目录：回到最前面
psscale GMT画比例尺的指令
一定要有的参数
-D 定义比例尺的位置和大小
如果比例尺是水平的 -D距中央/距顶端/长度/宽度h
如果比例尺是垂直的 -D距左端/距中央/长度/宽度
不一定要有的参数
-B 定义比例尺的间隔，和前面一样
-C 对照颜色，后面要接文件名称
-K 在这一行下面还有后续的指令
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
psscale -Cfile.cpt -D2.8/-1/4/0.22h -Bf1000a2000:"Elevation (m)": -O -K >> 123.ps
回到目录：回到最前面
grdcontour GMT画二维等高线图的指令
一定要有的参数
grd grd档案的档名
-C 等高线的间距
-J 投影方式，和前面一样
不一定要有的参数
-A 等高线数值标示的间隔
-B 和先前psbasemap的一样
-E 视角
-G 同一条等高线数值标示的间隔
-K 在这一行下面还有后续的指令
-L 等高线最低或最高的限制
-N 等高线的单位（如果有下-A的参数才有效）
-O 把这一行指令绘制的图迭加到先前画的图上
-P 把图旋转90度
-R 地图的经纬度范围，和前面的一样
-T 高地（-Ｔ+）或低地（-Ｔ-）加上标记
-V 画面上会显示本行指令的执行结果
-X 移动整个图的水平位置
-Y 移动整个图的垂直位置
表示方式
grdcontour file.grd -JM7 -R117/126/18/26 -C0.025 -W1 -A1 -P -K -V >> 123.ps
回到目录：回到最前面
grd数据的处理
转换成xyz值
grd2xyz file.grd -R117/126/18/26 > file.xyz
块状平均（I是精确度每多少度一格）
blockmean file.xyz -R117/126/18/26 -I0.1 > file2.xyz
转换成grd檔
surface file2.xyz -R117/126/18/26 -I0.1 -Gfile3.grd
查询档案信息
minmax file.xyz
grdinfo file.grd
qaidam
Qdm/user/byk/SUMMER/2003
/home_gobi/byk/SECONDARY/topo
回到目录：回到最前面
一些例子
制作CPT文件的程序 制作CPT文件程序的指令文件 ...跑出来的CPT檔
把地震画到地图上...图
把地震画到地图上3D...图
综合应用...图
立体地形图...图
回到目录：回到最前面