GMT系列六 pscoast和psxy命令详解

pscoast和psxy命令详解

pscoast是GMT最为常用的绘制地理底图的命令；psxy是绘制2D多边形、符号或线的命令。以下将详细介绍这两个命令的用法及其常用参数的意义。

1.      pscoast命令

作用：绘制地理底图

常用参数解释：

ü  -Amin_area[/min_level/max_level]        

通常绘制地图都需要用到分层数据，对于数据层级小于min_level或大于max_level的数据，在面积小于min_area的地理特征将不会被绘出。默认选项是0/0/4（绘制所有的地理特征），常用方法如-A500（只设置最小面积）等。

ü  -C[l|r/]fill

设置lakes和river-lakes的填充属性（默认是wet areas,等同于-S选项）。-Clfill，设置laker的填充属性，如果还想填充river，继续使用-C选项，即-Crfill

ü  -Drevolution[+]

设定数据集的分辨率，共有full, high, intermediate, low and crude五个选项。默认选项是low.

+表示如果选定的分辨率达不到的话，自动选择一个较低的分辨率。

ü  -Eazim/elev

设定透视投影的方位角和高程。

ü  -Gfill|c

设定dry areas的填充属性；或者用-Gc设定裁剪区域。

ü  -Iriver_type[/pen]

绘制河流的边界，river_type如下（如果想绘制多种类型的河流，重复-I命令）

River type	Comment
1	Permanent major rivers
2	Additional major rivers
3	Additional rivers
4	Minor rivers
5	Intermittent(间断的) rivers – major
6	Intermittent rivers – additional
7	Intermittent rivers – minor
8	Major canals(运河)
9	Minor canals
10	Irrigation canals (灌溉运河)
a	All rivers and canals(1-10)
r	All permanent rivers (1-4)
i	All intermittent rivers (5-7)
c	All canals(8-10)

ü  -Sfill|c

设定wet areas的填充属性；-Sc用于剪切区域。

ü  -W[level/]pen

绘制海岸线

ü  -L[f][x]lon0/lat0[/slon]/slat/length[m|n|k][+llabel][+jjust][+ppen][+ffill][+u]

在lon0/lat0（实际地理坐标）的位置绘制一个地图比例尺，-Lx用x/y（图形坐标）的位置代替lon0/lat0；

–Lf表示获得的是fancy比例尺；

slat:比例尺是在slat纬度上计算的（对于斜投影，可以提供slon参数）；

length: 绘制多长举例的比例尺（如1500km），默认单位是km，也可以将单位设置成miles或者nautical miles（海里）；

+l设定比例尺的label，与length后的单位对应；

+j设定label的对其方式，可以是left, right, top, bottom.

如果想在比例尺后面绘制一个矩形，可以用+p和+f来设定矩形边框和填充属性。

+u在比例尺所有的annotation后加上单位。

ü  -Nborder_type[/pen]

绘制政治边界。常用的政治边界有四类，如下：

Boder Type	Comment
1	National boundaries
2	State boundaries（仅限美国）
3	Marine boundaries
a	All boundaries

ü  -T[f|m][x]lon0/lat0/size[/info][:w,e,s,n:][+gint[/mint]]

在地图lon0/lat0的位置上绘制一个指北针。-Tx用x/y位置代替lon0/lat0.

size是指北针玫瑰图中玫瑰的直径。

[:w,e,s,n:]在::之间定制各个方向的label.

默认的指北针只绘制北方向，用-Tf可以绘制fancy的指北针，并可在[/info]里定制显示。info包含三种类型，1代表E-W, N-S指向；2增加了NW-SE和NE-SW指向；3又增加了四个方向。

Tm用于绘制磁罗经玫瑰图

ü  -U[just/dx/dy/][c|label]

绘制unix时间标签。just设定时间标签的对其方式（相对图形的左下角），dx/dy表示时间标签相对左下角的位移。如BL/0/0，时间标签就绘制在图形的左下角。

label，在时间后绘制标签。如-UBL/0/0”This is my first example”

ü  -Zzlevel

用于3D投影，设置海岸线的z方向的标签。

ü  -Jz|Zparameter

用于绘制3D图形，设置垂直方向的比例尺

2.      psxy命令

作用：在图形上绘制线、多边形和符号等。

用法：psxy files [options] > output.ps 或 command | psxy [options] > output.ps

详解：

ü  命令的输入：数据可以来源于给定的files（一个或多个文件），如果没有给定数据文件，可以通过管道命令将其他命令的输出作为psxy命令的输入。

ü  -A[m|p]

默认情况下，在绘制线段的时候，用大圆的圆弧来代替；-A选项改变这个默认选项，按直线的方式绘制线段；在绘制地图时，通常有纬圈和经圈两种大选，此时绘制直线时需要加后缀m（经圈）或p（纬圈）来选择按哪个大圆来绘制直线。

ü  -Ccptfile

给出调色板文件，用于绘制符号或多边形的颜色。如果有-S选项，符号的填充色会由输入文件的第三列Z值决定（此时，其他列都将会后移一列，如尺寸大小会从第三列移到到第四列）；如果没有给出-S选项，psxy命令希望输入数据还有多文件头（-m），每个文件头给出这一块数据的Z值，从而确定填充色。

ü  -Ddx[/dy]

在绘制线或符号时，使其位置偏离（offset）实际值dx/dy。默认的偏离值为0. 如果dy没有给出，默认其值等于dx。

ü  -E[x|y|X|Y][n][cap][/[-|+]pen]

绘制误差条，x,y表示绘制x的误差条或者y的误差条。x或者y误差存储在(x,y)或(x,y,size)的后两列。cap表示绘制的误差条终端线条的长度。图中蓝色显示的为误差条，cap控制误差条两端短线的长度。[+-]是在psxy命令调用了cptfile时，用来控制符号和误差条的颜色属性。

ü  -Gfill

填充符号或者多边形的颜色。

ü  -Iintens

模仿自发光的形式，利用提供的强度值（通常在-1~1之间）调制填充色。

ü  -L

绘制封闭的多边形或者将线段连接起来

ü  -N

对于那些超过地图边界的符号，同样要绘制（默认情况下绘制在地图边界上）。该选项对于绘制多边形和直线不起作用。

ü  -S[symbol][size]

绘制符号。如果没有给出符号的尺寸，此时会自动默认输入数据的第三列（有-C选项时，为数据的第四列）为每个符号的尺寸大小（不同数据用不同大小的符号来标示，以显示区别）。如果符号类型没有给出来的话，会默认输入数据的最后一列是符号类型。符号类型有很多种，具体参见man psxy，这里仅给出一些常用的类型。

-S+      +                -Sa          星形         -Sc            圆形

-Sd      菱形         -Sg            八边形    -Sh            六边形

-Ss      正方形    -St             三角形    -Sv            矢量        

-Sx      叉

-Se      椭圆形，需要给出方向（从水平方向按逆时针计算）、长轴和短轴（通常由数据的三四五列给出）

-Sk      绘制一些特殊符号（特殊符号的定制语句放在一个单独的文件里）

ü  -T

忽略所有输入文件或者标准输入。

ü  -W[-|+][pen]

设置绘制直线或者符号轮廓的画笔属性

       注：如果没有给出S的话，将绘制直线。