pan11的个人博客分享 http://blog.sciencenet.cn/u/pan11

博文

gmt4的一些简单命令

已有 1898 次阅读 2018-1-14 14:24 |个人分类:程序技巧|系统分类:科研笔记| gps, GMT

感谢群内吧友 地学gamit群:518498135

自己习惯了gmt4的用法,一直使用,仅仅介绍一些命令的简单作用,但现在陆续开始应用GMT5,一些用法稍微改进,陆续会有更新gmt5的用法,详情见gmt主页http://gmt.soest.hawaii.edu/gmt4/gmt/html/gmt_services.html

以及主页http://gmt.soest.hawaii.edu/


-R  设定范围 xmin/xmax/ymin/ymax 定义直角坐标地理坐标系

         xleft/yleft/xuright/yurightr 后面加r 这个表示倾斜的投影直角地理 描述最下面左边的点和最右边上边的点坐标

         grdfile  描述

         -Rd代表 -180/180/-90/90 -Rg表示0/360/-90/90

         -R2009-08-01T/2009-09-30T/25/40

-B  设定坐标轴间距a标注间隔f为框架的间隔g网格化

   WSne为大写为出现方位轴 小写为隐藏该坐标轴

-X  移动图的左右位置

-Y  移动图的上下位置

-O  链接上图

-K  链接下图

-H  头文件的个数

-P  选择portrait plot模式输出

-U  画出系统时间图附在下面

-V  描述完成过程在显示屏中

注意事项:

1、不要省略参数,在每一行参数使用变量pscoast -J$J -R$R -B$B -W1p -A1000 -K > $PS

                              psxy -J$J -R$R -Sa0.5c -Gred -O >> $PS << EOF

2、不要省略单位,一般默认Jx10 为c厘米,添加上

  1. 尽量不要在$HOME下执行GMT命令(可能会产生临时文件,难以清理)

  2. 尽量不要使用命令行执行GMT命令(因为你很可能会忘记你刚刚执行过哪些命令)

  3. 使用gmtset的脚本,最后一定要记得删除.gmtdefaults4

3、对于-B选项,合理的用法是仅在第一个命令中使用







1、 blockmean      L2(x,y,z)值的filter或者是decimator(每十个删去一个)

2、 blockmedian    L1(x,y,z)tables值的filer或是decimator(每十个删去一个)

3、 blockmode      众数估计(x,y,z)table值的filter或者是(每十个删去一个)

4、 fiter1d        filter 一维table值(时间序列)

5、 fitcircle      对于一系列的点找到最佳适合的大的圆或者小的圆

6、 gmt2rgb        将sun raster或者grid的文件转换为红蓝绿的grid

7、 gmtconvert     将数据table 从一个格式到另一个格式

8、 gmtdefaults    罗列目前默认的设置

9、 gmtmath        对于table数据进行的数学计算

                  对于ascll数据进行第三列计算

10、gmtselect      基于多个空间标准选择table数据的子集

11、gmtset         改变在.gmtdefaults文件中的参数

12、grd2cpt        从grd文件中制作cpt文件

13、grd2xyz        将grd文件转换为table文件

14、grdblend       混杂合成几个部分重复重叠的grd一个新的grd文件

15、grdclip        在已经grdding的数据束中限制z的范围

                  对grd文件进行修补 约束 -Sa 大于的值设置 -Sb 小于的值设置

16、grdcontour     将二维网格化数据等高线

17、grdcut         从一个grd文件中生成子grd文件

18、grdedit        修改一个二维grd文件的头信息

19、grdfft         在频率域执行操作在grd文件

20、grdfilter      在空间域中filter的grd文件

21、grdgradient    在grd文件中计算方向导数

22、grdhisteq      均等化柱状图grd文件

23、grdimage       从一个二维grd文件中生成图像

24、grdinfo        得到grd信息

25、grdlandmask    从shoreline 数据包中生成masking数据

26、grdmask        重新设置grd的node输入或者输出clip路径到常量中

                  主要用于将与想要的数据 距离 或者用度数 设置为 多少 之远 范围外进行屏蔽  想要的数据进行成图,不想要的数据不显示

27、grdmath        对于grd文件的数学计算

                 主要是第三列进行计算科学运算

28、grdpaste       沿着公共的边缘粘帖grd文件

29、grdproject     投射grd数据到一个新的坐标系统

30、grdreformat    将grd格式转变为其他格式

31、grdsample      重新取样一个二维grd文件到一个新的grd文件

32、grdtrack       沿着一维的轨迹进行对二维grd文件取样 沿着一系列的数据进行投影

                  grdtrack 沿着xy的剖面进行对grd文件的取样,test文件的前两列必须与grd文件的前两列性质一样,

                  grdtrack  test -Gmi.grd > junk.dat 得到的junk文件是源文件后面加一列取样值

33、grdtrend       多项式拟合趋势grd文件

34、grdvector      画出二维grd文件的矢量图

35、grdview        从3维的角度画出二维grd文件

36、grdvolume      计算在一个特定等高线的表面的体积

37、greenspline    在一维二维三维 的数据进行格林函数的插值

                 最小曲率样条插值 -Sc -D1  笛卡尔坐标  -sp -D3 球面坐标

                 连续曲率样条插值 -St -D1  笛卡尔坐标  -sq -D3 球面坐标

                 

                 

38、makecpt        制作cpt文件

39、mapproject     转换table数据的坐标系统


40、minmax         显示table数据的特殊值

                  -I 显示R范围

41、nearneighbor   near-neighbor网格化

                 插值方法

42、project        投射table数据到一条线或者大圆中

                  三种方式 1 用-A 和-C 方位角和中心位置经纬度

                         2 用-C 和-E 起点和终点位置

                         3用-C和-T 中心位置和旋转轴位置

                  -C 起点x/y

                  -E 投影的终点x/y

                  -G 生成的间隔

                  -A 方位角 从北向开始 为o度

                  -L 为投影长度

                  -W 为投影宽度

                  -F 输出格式 xyzpqrs 分别为 xy为点的坐标,z为剩余不变的东西 pq为平行和垂直于投影面的距离 rs为投影到新的截面投影点的对应的原图的坐标

43、ps2raster      修剪转换ps文件为raster图像,eps格式,或者pdf格式

                  转换为jpeg的格式 ps2raster file.ps -A -Tj

44、psbasemap      生成底图

45、psclip         用多边形文件定义修剪路径

46、pscoast        画出海岸线、边界和河流

47、pscontour      等高线画图或者三角剖分table数据

48、pshistogram    画直方图

49、psimage        画sun raster文件

50、pslegend       画图例

51、psmask         建立一些覆盖遮住图像的一些区域

52、psrose         画扇形或者玫瑰图表

53、psscale        画灰色比例或者颜色比例

                  -D定义中心/上方(水平),中心/左侧(垂直),然后是图例的大小比例长/宽,默认是垂直,后加h为水平

                   -D6.5i/2i/7.5c/1.25c 7.5长1.25宽,离画图6i宽,4i高

54、pstext         画字符串

                  "4.5 -1.2 27 0 1 LM <-- y=sin(x)+cos(x)"

                  x坐标位置 y坐标位置 字体大小 旋转角度 字型 表示位置LM为左边的中间  表示内容

                  L为左边 C为中间 CM为中间的中心 R为右边 TR为右边的上边 M为中心 B为下边 LB为左边的下边

55、pswiggle       画table数据时间序列沿着轨迹

56、psxy           画符号多边形线

                  psxy -R0/241/0/6000 -JX6.3i/1.4i track4.d -W1p -X3.5 -Y6 -Ba40f10:"km":/a1000f500:."Profile_DD"::"Elev(m)":WnSe -K >profile_D.ps

                  -R范围,-J投影方式,文件可以为grd文件,-W画笔,-B

                  在文件中拥有> 然后加数据 在psxy中加 -M 为画多条线段,若在文件中没有> 则是仅仅画点。

                  有三列,psxy将解释第三列作为符号大小,如果没有符号说明,用了-S,则解释数据的最后一列,

57、psxyz          画符号多边形线在三维

58、sample1d       重新取样一维数据

59、spectrum1d     从时间序列中计算不同的谱估计

60、splitxyz       拆分xyz数据文件几部分

61、surface        连续曲率grid算法

62、trend1d        拟合多项式或者傅立叶趋势到y=f(x)

                  trend1d ex1.dat -Nf10 -Fxmr > ex1.fit

                  #-F为x=x,y=y,m=modelf(x),r=residual y-m 输出残差到新的数据为后面最后一列  -w在拟合中用到加权为在原数据中第三列为

                   -N f为傅立叶拟合一般用于又cossin项的  加r为强制约束拟合,数字2为最小二乘,多项式拟合。

                  去趋势-Fxr  不加m,拟合用m 后面加r输出残差(拟合的数据与原先数据的差值)                

63、trend2d        拟合趋势二维z=f(x,y)

64、triangulate    执行最佳的delaunay插值网格化

65、xyz2grd        转变xyz数据为grd文件

66、grdraster      从一个二进制raster文件中提取子区域并写入grd文件

67、psmeca         画震源机制图

                 -R 范围 -B设置边界间隔 -J 设置投影 infile输入文件为多个或者一个

                 -Ddmin/dmax 只绘制某个深度范围之内的地震

                 -E 震源求拉伸象限的颜色 -e指定T轴的颜色  -p指定p轴线条属性

                 -G 震源球压缩象限的颜色 -g指定p轴颜色   -t指定t轴线条属性

                 -L 绘制震源求外部轮廓的属性            -r在文本周围加上方框

                 -M  默认震源求大小与震级大小成正比,这个代表所有震源求一样大小

                 -Z  指定cpt文件 根据数据文件的第三列值以及cpt文件指定深度

                 -C  将最后的两列变成新的经纬度,然后用c进行线条颜色属性偏移

                 -a  计算并绘制p轴与t轴 可以指定大小与P轴t轴的符号

                 -Sc global cmt格式 地震矩单位dyne-cm  scale指定的为震级为5的震源球大小相当于4.0e23 dyne.cm 尾数为4 指数为23

                     x y depth strike1 dip1 strike2 dip2 rake2 mantissa(地震矩的尾数) exponent(地震矩的指数) newx newy title

                 -Sa 按照Aki&Richards的约定

                     x  y depth strike dip rake mag newx newy title

                 -Sm global cmt格式 默认为这个格式

                    x   y depth mrr mtt mff mrt mrf mtf exp newx newy title

68、pscoupe        画震源机制的十字剖面   命令基本与pemeca一样

                  -Aalon1/lat1/lon2/lat2/dip/p_width/dmin/dmax

                   两个断点的经纬度,dip为投影平面的倾角,p_width为投影平面两侧的宽度 后面为距离水平面的最小和最大距离

                  -Ablon1/lat1/strike/P_length/dip/p_width/dmin/dmax

69、pspolar        画震源机制下半球的极性

70、psvelo         画速度矢量,十字画以及契形图

71、nc2xy          转化netcdf格式文件为ascii xy文件

72、sphdistance    计算最近的距离来自空间球面数据的voronoi

73、sphtriangulate 对球面空间数据进行最佳delaunay插值或者voronoi构建

74、sphinterpolate 张力网格化数据对于球面空间数据

75、grdrotater     对grd文件进行欧拉极旋转




http://blog.sciencenet.cn/blog-651121-1094803.html

上一篇:QOCA命令简介以及简单使用
下一篇:GLOBK——正确求解坐标点位

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备14006957 )

GMT+8, 2018-9-22 03:58

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部