GMT笔记

来源：互联网发布：淘宝买家秀兼职骗局编辑：程序博客网时间：2024/04/30 18:35

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...图
 综合应用...图
 立体地形图...图

回到目录：回到最前面