GMT笔记（常用）

葡萄

1.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
2.pscoast时间画陆块海洋和海岸线的程式
一定要有的参数
- 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
3.psxy时间在地图上画点流水生产线的程式
一定要有的参数
- 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
12325
>
121 23
12327
在123.dat的档案里是这样（这是点用的档案）
120 25
12325
121 23
12327
表示方式三
把更多的资料打在资料档案里
psxy 123.dat -R117/126/18/26 -JM7 -Sc -Ccpt >> 123.ps
在123.dat的档案里要有四，栏别名，依序为x，y，颜色代码，圆的大小
4.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
5.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
6.grdcontour GMT画二维等高线图的指令
一定要有的参数
grd grd档案的档名
-C等高线的间距
- J投影方式，和前面一样
不一定要有的参数
- A等高线数值标示的间隔
- B和先前psbasemap的一样
- E视角
- G同一条等高线数值标示的间隔
- K在这一行下面还有后续的指令
-L等高线最低或最高的限制
-N等高线的单位（如果有下 - A的参数才有效值）
-O把这一行指令绘制的图叠加到先前画的图上
-P把图旋转90度
-R地图的经纬度范围，和前面的一样
-T高地（ - T+）或低地（-T）加上标记
-V画面上会显示本行指令的执行结果
-X移动整个图的水平位置
-Y移动整个图的垂直位置
表示方式
grdcontour file.grd -JM7 -R117/126/18/26 -C0.025 -W1 -A1 -P -K -V >> 123.ps
7.grd资料的处理
转换成XYZ值
grd2xyz file.grd -R117/126/18/26 > file.xyz
块状平均（我什么精确度每多少度-报关员资格）
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