GMT笔记GMT常用画图等命令.docx

GMT笔记GMT常用画图等命令.docx

《GMT笔记GMT常用画图等命令.docx》由会员分享，可在线阅读，更多相关《GMT笔记GMT常用画图等命令.docx（19页珍藏版）》请在冰豆网上搜索。

GMT笔记GMT常用画图等命令

命令：

backtrackerForwardandBackwardflowlinesandhotspottracks

binlegs维护GMT–SYSTEM索引文件

blockmeanL2（x,y,z）tabledatafilter/decimator

blockmedianL1（x,y,z）tabledatafilter/decimator

blockmodeModeestimate（x,y,z）tabledatafilter/decimator

cps把shellscripts和数据文件插入到CompletePostScrip（tCPS）文件

cpsdecode由CompletePostScript（CPS）文件提取shellscripts和数据文件

cpsencode把shellscripts和数据文件插入到CompletePostScript（CPS）文件

dat2gmt把一个ASCII文件转换为二进位gmt文件

filter1dFilter1–Dtabledatasets（timeseries）

fitcircle得到一组数据点的最佳拟合大圆或小圆

gmtGMT数据处理和显示软件包

gmt2bin由gmt文件创建二进制索引文件

gmt2dat二进制gmt文件转换为ASCII文件

gmtconvertASCII和二进制1–D表文件之间转换

gmtconvert数据表文件格式转换

gmtdefaults当前缺省设置列表

gmtinfo得到单独航次的信息

gmtlegs查明一个给定地区的航次

gmtlist由.gmt文件提取数据

gmtmath使用PolishNotation语法对表格数据进行代数操作

gmtpath获得gmt文件的完全路径

gmtselect选择基于多空间标准的数据子集

gmtset改变当前的.gmtdefaults文件所选择的参数

gmttrack一个航迹绘图程序

grd2cpt由一个网格文件创建色谱表

grd2xyz2–D网格文件转换为表格数据

grdclip限定网格数据组的z值范围

grdcontour绘制2–D网格数据组等值线

grdcut由一个网格文件抽取一个子区就是裁剪掉一部分只要感兴趣的部分

grdcutingrid-Goutgrid-Rrang

grdedit改变一个2–D网格文件的标头信息

grdfftPerformoperationsongriddedfilesinthefrequencydomain

grdfilterFilter2–Dgriddeddatasetsinthespacedomain

grdgradient由网格文件计算方向梯度,得到另一个网格文件

grdhisteq网格文件柱状图均化（equalization）

grdimage由一个2–D网格数据组创建影像图

grdinfo得到网格文件的信息

grdlandmaskCreatemaskinggriddedfilesfromshorelinedatabase

grdmaskResetgridnodesin/outsideaclippathtoconstants

grdmath对网格文件进行代数操作

grdpaste沿一条公共边界合并多个网格文件

grdproject把网格数据组投影到一个新的坐标系

grdrasterExtractsubregionfromabinaryrasterandwriteagrdfile

grdreformat把网格文件转换为另一种格式

grdsample对一个2-D网格文件重新采样创建一个新的网格

grdtrack沿1-D轨迹对一个2-D网格数据组重新采样

grdtrend对网格文件进行多项式趋势拟合

grdvector绘制2-D网格矢量区

grdview由一个2-D网格数据组创建3–D透视影像图

grdvolume计算给定等值线以下的体积

hotspotterCreateCVAimagefromseamountflowlines

img2mercgrdExtractregionofimg,preservingMercator,saveasgrd

makecpt创建色谱表文件

makepatternMakeGMTcolorpatternfromb/wpatternoricon

mapprojectTransformationofcoordinatesystemsfortabledata

mgd77togmtConvertanMGD–77asciifiletoabinarygmtfile

minmax报告数据表文件极大极小值

nearneighborNearest–neighborgriddingscheme

originatorAssociateseamountswithhotspotpointsources

projectProjecttabledataontolinesorgreatcircles

psbasemap创建底图

psclipUsepolygonfilestodefineclippingpaths

pscoast在图上绘制（和充填）海岸线、国界线和河流

pscontourContourorimagerawtabledatabytriangulation

pscoupePlotcross–sectionsoffocalmechanisms.

pshistogram绘制柱状图

psimagePlotSunrasterfilesonamap

pslibv3.2APostScriptbasedplottinglibrary

psmaskCreateoverlaytomaskoutregionsonmaps

psmecaPlotfocalmechanismsonmaps

psmegaplotTocreateposter–sizePostScriptplotsfrompage–sizeplot

pspolarPlotpolaritiesontheinferiorfocalhalf–sphereonmaps

psrose绘制扇形或玫瑰花图

psscale在图上绘制灰度图例或彩色图例

pssegyCreateimagemaskedpostscriptfromSEGYfile

pssegyzCreateimagemaskedpostscriptfromSEGYfile

pstext在图上绘制字串

psveloPlotvelocityvectors,crosses,andwedgesonmaps

pswiggle沿地图上一条规机绘制根据时间序列数据绘制曲线

psxy在图上绘制符号、多边形和线条

psxyz在图上绘制3–D符号、多边形和线条

sample1d表格数据组重新采样

spectrum1dComputevariousspectralestimatesfromtime–series

splitxyz把xyz文件分为几段

surface一种连续弯曲的网格算法

trend1dFitspolynomialorFouriertrendstoy=f（x）series

trend2dFitspolynomialtrendstoz=f（x,y）series

triangulatePerformoptimalDelauneytriangulationandgridding

参数：

-B设定图边界尺标的间隔

-H选项让GMT知道输入数据文件具有一个（缺省）或多个头纪录。

如果多于1个头纪录，一定要在-H选项后指定数字。

如-H4。

-O

-K

-P选择肖像模式。

缺省为风景画模式

-R指定绘图范围

-U系统时间标志

-V选择verbose模式，将建立程序运行报告文件stderr[Defaultruns``silently''].

-X移动整个图的水平位置

-Y移动整个图的竖直位置

-:

对于地理数据，通常认为第一列为经度，第二列为纬度值。

如果要调换经纬度顺序，必须使用–:

选项。

1.>psbasemapGMT画底图的指令

 一定要有的参数

   -B 设定图边界尺标的间隔

     a文字批注的间隔

     f边框标记的间隔

     g隔线的间隔

       如果要针对X轴或Y轴做不同的设定，就在参数间加斜线/

     NEWS 每个轴加上标记

     news 每个轴不加标记

     :

.title:

 在图上面加title这个标题

     :

title:

在轴上加title这个标题（只有在-Jx的情况下有用）

:

%:

刻度单位“%”的表示方法，其它亦然

     例如-BNEWsf2a3g2:

xx:

/f3a3g2:

yy:

:

.title:

 -JX（-Jx） 线性投影（LinearProjection）一般绘图用

     -Jx后面接的数字代表每单位的长度（预设长度单位为英吋，加c的话就变公分）

      -JX后面接的数字代表轴的全长

      如果要针对X轴或Y轴做不同的设定，就在数字间加斜线/

 -JM（-Jm） 麦卡托投影（MercatorProjection）地图绘图用

     -Jm后面接的数字代表每度的长度（预设长度单位为英吋，加c的话就变公分）

      -JM后面接的数字代表整张图的宽

      如果要针对X轴或Y轴做不同的设定，就在数字间加斜线/

-JGlon0/lat0/width.lon0/lat0指定投影中心

正形投影（Orthographic:

）

-JBlon0/lat0/lat1/lat2/width.给出投影中心和2条标准平行线

Albers圆锥投影（Albersconic:

）

-JK[f|s]lon0/width.给出中央子午线

Eckert投影（EckertIVandVI:

）

-R 地图的经纬度范围

     西/东/南/北 依序填入经纬度

   >out.ps 把画出来的图输出到out.ps这个档案

 不一定要有的参数

   -E 旋转地图的视野方向

   -G 把图着色，用r/g/b调

　　-K在这一行下面还有后续的指令

   -O 把这一行指令绘制的图迭加到先前画的图上

   -P 把图旋转90度

   -V 画面上会显示本行指令的执行结果

   -X 移动整个图的水平位置

   -Y 移动整个图的垂直位置

 表示方式

psbasemap-R10/70/-3/8-JX4i/3i-B10a4f3g5:

x:

/1a0.6f0.5g1:

y:

:

."Myfirstplot":

EN-P>plot.ps

psbasemap-R10/70/-3/8-JX4i/3i-B10a10f3g5:

cm:

:

x:

/1a1.1f0.5g1:

y:

:

%:

:

."Myfirstplot1":

WS-P>plot.ps

 psbasemap-R0/100/0/100-B20NEWsf3g10:

x:

/20f5g10:

y:

:

.title:

-JX3i-E30/20-G255/155/155-P-V-X2i-Y2i>basemap_string.ps

-E30/20效果图：

2.>pscoastGMT画陆块海洋和海岸线的程序

 一定要有的参数

   -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.>psxyGMT在地图上画点线的程序

 一定要有的参数

-J 投影法，和前面的一样

-R 地图的经纬度范围，和前面的一样

不一定要有的参数

-A 画大圆线

-B 和先前psbasemap的一样

-C 对照颜色，后面要接文件名称,分别设定每个符号的颜色

-Gpdpi/pattern[:

Br/g/b[Fr/g/b]] 或–Gfill：

把符号着色，用r/g/b调色pattern可以是1–90范围的数字，或1–、8–或24–位的Sun光栅文件的名字，如用大写字母–GP代替–Gp，图像将反相，即黑白像素交换（只能用于1–bit图像或预设的位图图案）使用–G的光栅图像大小一定要小于146x146像素。

-K 在这一行下面还有后续的指令

-N 不要跳过超出范围的点

-O 把这一行指令绘制的图迭加到先前画的图上

-P 把图旋转90度

-M 线分段

-Sa在图上画星星，后面要接大小

 -Sc在图上画圆圈，后面要接大小

-Sd在图上画钻石，后面要接大小

-Se在图上画椭圆形，在资料的第3,4,5栏要有方向长轴短轴

-SE椭圆，方位角azimuth（由垂直方向顺时针计），以公里做单位的长轴major、和短轴minor由输入数据文件读取

-Sfgap/tick[l|L|r|R]断层，gap和tick设置短线的长度和间隔，如gap<0即指短线编号；l或r表示短线画在断层线的左侧或右侧[缺省画在线上]，大写L或R表示用三角形代替短线。

-Slsize/string[%font]字符，size为字体大小；增加一个字符或一个字串，字体可任选。

-Sp点，无须设大小（当前分辨率下的1个像素）

-Sh在图上画六角形，后面要接大小

-Ss在图上画正方型，后面要接大小

-St在图上画三角形，后面要接大小

-Si在图上画反三角，后面要接大小

-Sv[thick/length/width][nnorm]在图上画箭头，在数据的第3,4栏要有方向长度

-SV[thick/length/width][nnorm]箭头，用方位角azimuth（由北向东旋转的角度）代替方向direction。

地图的角度根据所选的地图投影来计算

-Sx在图上画叉叉，后面要接大小

-Sw[size扇形，起点和终点的方向（由水平方向逆时针计）由输入数据文件读取。

-V 画面上会显示本行指令的执行结果

-Wwidth[/color][texture][p] 笔宽,注意，根据你的设备的dpi的不同，笔宽为

5可有不同的物理宽度，但是5p则总是5/72英寸，texture属性控制了线的结构

-X 移动整个图的水平位置

-Y 移动整个图的垂直位置

通用的输入数据的格式为（可选项用方括号表示）：

x,y[z][size][σx][σy][symbol]

只有前两列是输入文件所必需的，它们一定是经度和纬度（或x和y）。

其他可选项在遇到以下情况之一时使用：

1.如果你打算分别设定每个符号的颜色，用–C提供一个cpt文件，增加数据文件第3列，让z值在cpt文件所取值范围内。

2.如果你打算分别设定每个符号的大小，在第4列给出符号大小。

3.为绘制误差棒，使用–E选项，并增加1列或2列数据，分别为±dx和±dy值。

–E的形式决定需要1列（–Ex或–Ey）还是2列（–Exy）数据。

4.如果你使用–Sv（或–SV）绘制箭头，则size实际上为2列，包括了每个箭头的方向direction（或方位azimuth）和长度。

5.如果你使用–Se绘制椭圆，则size实际上为3列，包括了每个椭圆的方向direction，以英寸做单位的长轴major和短轴minor（用–SE时，改为方位角azimuth，轴长以公里做单位）.

psxytable_5.11-R0/10/0/10-JX5i-B2/2-G155/255/155-Sa3-V-W1>psxy_xingxing.ps

psxytable_5.11-R0/10/0/10-JX5i-B2/2-G155/255/155-Sb0.5-V-W1>psxy_xingxing.ps

psxytable_5.11-R0/10/0/10-JX5i-B2/2-G155/255/155-V-W1>psxy_xingxing.ps效果图：

psxytable_5.11-R0/10/0/10-JX5i-B2/2–V-W15/255/155/125to>psxy_xingxing.ps

-W15/255/155/125ta

-W15/255/155/125"20_10_5_10:

10"10表示距离原点10个单位

表示方法二：

–W2/255/255/0t20_10_5_10:

10

带“to”点线，

带“ta”虚线，

自定义线：

绘制一个与线型相关的符号：

psxytable_5.11-R0/10/0/10-JX5i-B2/2-V-K-W1/255/155/125>psxy_xingxing.ps

psxytable_5.11-R0/10/0/10-JX5i-B2/2-V-O-Si0.2i-W1/255/155/125>>psxy_xingxing.ps

4.>pstext

–Cdx/dy设定文本和文本框之间的间隔（见–W）

–Ddx/dy设定字串投影位置偏移量

–Gfill设定文本颜色

–L列出字体标识符并退出

–N保持图边界符号完整

–Spen选择空心字体并设置画笔属性

–W[fill][o[pen]]打印文本框，Paintthetextbox;drawtheoutlineifoisappended（alsosee–C）

用pstext输入数据可望含有以下信息：

xysizeanglefontnojustifytext

7个参数的意义为：

xyX坐标，Y坐标；size以点做单位的字体大小；angle文本基线和水平线之间的角度（逆时针测算）；fontno选用字符集代码；justify为字串对齐方式选择，设定文本串相对于x、y坐标的位置。

相应的代码为1~LT，2~CT，3~RT，4~LM，5~CM，6~RM，7~LB，8~CB，9~RB；text为字符串或待标注文本

@+上标使用转换开关

@–下标使用转换开关，:

.title@-2:

@#小型大写字符使用转换开关,:

.titleA@#A:

@!

由后面两个字符创建一个合成符号，:

.titleA@#AE:

@@打印符号@本身

@E@e，:

.title@E@e:

@O@o，:

.title@O@o:

@A@a，:

.title@A@a:

5.>grdimageGMT画二维等高图的指令

   grdgrd档案的档名

-Edpi设置所需要的影像分辨率[缺省为数据分辨率]

   -C 对照颜色，后面要接文件名称

-J 投影方式，和前面一样

-Iintenfile使用intensfile文件的亮度数据建立虚拟照明

-M使用（电视）YIQ转换改为灰度

-B 和先前psbasemap的一样

-K 在这一行下面还有后续的指令

-O 把这一行指令绘制的图迭加到先前画的图上

   -P 把图旋转90度

   -R 地图的经纬度范围，和前面的一样

   -V 画面上会显示本行指令的执行结果

   -X 移动整个图的水平位置

   -Y 移动整个图的垂直位置

 表示方式

grdimagefile.grd-JOc122/21/22/69/4.8i-R117/126/18/26-Cfile.cpt-P-K-V>123.ps

6.>pscontour轮廓表中的数据直接三角剖分[方法]

pscontour-R-Jtable_5.11-B2f1WSne-Wthin-Ctopo.cpt-Lthinnest,--G1i/0>123.ps:

以右上图为相对坐标，绘左下图图框、标注；绘各条边组成三角形网，虚线；根据topo.cpt的间隔画等值线，标注间隔1i。

pscontour-R-Jtable_5.11-B2f1eSnw-Ctopo.cpt-I-X3.25i-O-K>>example_12.ps

以左下图为相对坐标，绘右下图图框、标注；根据topo.cpt的间隔画等值线和填充颜色

pscontour-R-JX6itable_5.11-B2f1WSne-Wthin-Ctopo.cpt-G1i/0-K-U>123.ps

pscontour-R-Jtable_5.11-B2f1WSne-Wthin-Ctopo.cpt-Lthinnest,--G1i/0>12.ps

7.>psscaleGMT画比例尺的指令

  一定要有的参数

   –Dxpos/ypos/length/width[h]设定色谱图例的左侧中心坐标和图例的长／宽。

加h改为水平图例并给

出顶部中心坐标

 不一定要有的参数

   -B 定义比例尺的间隔，和前面一样

   -Ccptfile必需的cpt文件

　-K 在这一行下面还有后续的指令

   -O 把这一行指令绘制的图迭加到先前画的图上

   -P 把图旋转90度

   -V 画面上会显示本行指令的执行结果

   -X 移动整个图的水平位置

   -Y 移动整个图的垂直位置

-W 笔宽

 表示方式