一个 GMT 命令由“gmt + 模块 + 选项 + 参数”构成,写成如下形式:
gmt module -Axx+bxxxx -Bxx+axxxx

• gmt 是 GMT 中“唯一”的一个二进制文件,所有 GMT 命令必须以 gmt 开头
• module 是用于完成某个特定操作的模块的名字
• -A 是模块 module 提供的选项,xx 为选项 -A 的参数
• +b 是选项 -A 的子选项,xxxx 为该子选项的参数
一个完整的示例:
gmt pscoast -R0/// -JM6i -Ggray -Wthin -B5 -B+t"Title with spaces" -V -P > map.ps

• 不以 - 开头的参数,都会被当做文件,GMT 会尝试去读取
• 子选项以 + 开头,后接单个字符以及子选项的参数
官方文档中子选项称为 modifier,可以译为“修饰符”,本文档中统一使用“子选项”
子选项以 + 开头,后接单个字符以及子选项的参数
各选项间以空格分隔,选项内部不能有空格
• 选项内部的字符串,若存在空格,应用单引号或双引号括起来
STDOUT)用于显示输出数据,标准错误流(STDERR)用于显示错误消息
. 将标准输出流重定向到数据文件中
. 将标准错误流重定向到日志文件中
. 将标准输出流通过管道传递给下一个命令
对于标准输出流:
• > :将标准输出流重定向到新文件中。若该文件已存在,则覆盖文件中原内容;若该文件不存在,则
创建该文件
• >> :将标准输出流追加到文件中。若文件已存在,则将标准输出流追加到已有文件后面;若文件不
存在,则创建该文件
对于标准错误流,重定向符号是 2> 和 2>> 。这里的 2 表示标准错误流,大于号的含义与标准输出
流相同。

-K  used in the starting plot;-K -O used in processing plots;-O used in the final plot;

管道
举例如下,假设文件 input.dat 中包含了一系列地震的经度、纬度和震级共三列数据,想要在图上
画很多圆表示地震的位置,圆的大小表示震级的大小。可以用类似如下命令:
gmt psxy input.dat -Rxxx ... > test.ps
此时 psxy 模块会读取 input.dat 文件的内容作为其输入。
也可以使用管道:
cat input.dat | gmt psxy ... > test.ps
cat 命令会读取 input.dat 的内容并将其发送到标准输出流,由于使用了管道,标准输出流中的内
容被 gmt psxy 接收作为自己的标准输入流。
当然还可以使用 gawk
gawk '{print $1, $2, $3/10}' input.dat | gmt psxy ... > test.ps
gawk 会读取 input.dat 的内容,并对数据做简单处理并输出。
4. 通过管道输入
管道可以将前一个命令的标准输出作为后一个命令的标准输入:
cat event.loc | gmt psxy -R70/140/20/60 -JM6i -B5 -Sc0.2c -Gred -P > test.ps
5. Here Documents
示例如下,两个 EOF 之间的所有数据都会被传递给 GMT:
gmt psxy -R70/140/20/60 -JM6i -B5 -Sc0.2c -Gred -P > map.ps << EOF
100.0 40.0
110.0 45.0
EOF

Quoting:

想要用 psxy 绘制某数据时,需要提供数据的范围 -R ,而 gmtinfo 模块可以用于计算并输出
数据的范围,即需要将 gmtinfo 的输出作为 psxy 的一个选项。
比如:
$ gmt info in.dat -I1/1
-R0/10/0/10
$ gmt psxy in.dat -JX10c -R0/10/0/10 > map.ps
上面的做法需要人工干预,不适合脚本自动化,可以利用倒引号将 gmtinfo 的输出保存到变量中:
#!/bin/bash
R=`gmt info input -I1/1`
gmt psxy in.dat -JX10c $R > map.ps
上面的例子还可以进一步简化。此处变量 $R 只需要用一次,因而没有必要把 gmtinfo 的输出信息
保存到变量中,可以直接在 psxy 命令中使用倒引号:
$ gmt psxy in.dat -JX10c `gmt info in.dat -I1/1` > map.ps
EOF 可以被替换成其他任意字符(比如 END ),只要保证开始和结束的符
号一致即可

CMT learning的更多相关文章

  1. 【Machine Learning】KNN算法虹膜图片识别

    K-近邻算法虹膜图片识别实战 作者:白宁超 2017年1月3日18:26:33 摘要:随着机器学习和深度学习的热潮,各种图书层出不穷.然而多数是基础理论知识介绍,缺乏实现的深入理解.本系列文章是作者结 ...

  2. 【Machine Learning】Python开发工具:Anaconda+Sublime

    Python开发工具:Anaconda+Sublime 作者:白宁超 2016年12月23日21:24:51 摘要:随着机器学习和深度学习的热潮,各种图书层出不穷.然而多数是基础理论知识介绍,缺乏实现 ...

  3. 【Machine Learning】机器学习及其基础概念简介

    机器学习及其基础概念简介 作者:白宁超 2016年12月23日21:24:51 摘要:随着机器学习和深度学习的热潮,各种图书层出不穷.然而多数是基础理论知识介绍,缺乏实现的深入理解.本系列文章是作者结 ...

  4. 【Machine Learning】决策树案例:基于python的商品购买能力预测系统

    决策树在商品购买能力预测案例中的算法实现 作者:白宁超 2016年12月24日22:05:42 摘要:随着机器学习和深度学习的热潮,各种图书层出不穷.然而多数是基础理论知识介绍,缺乏实现的深入理解.本 ...

  5. Deep learning:五十一(CNN的反向求导及练习)

    前言: CNN作为DL中最成功的模型之一,有必要对其更进一步研究它.虽然在前面的博文Stacked CNN简单介绍中有大概介绍过CNN的使用,不过那是有个前提的:CNN中的参数必须已提前学习好.而本文 ...

  6. Programming Learning - Based on Project

    Today when taking a bath I got a good idea that it is an efficient and interesting way to learn a ne ...

  7. 做中学(Learning by Doing)之背单词-扇贝网推荐

    做中学(Learning by Doing)之背单词-扇贝网推荐 看完杨贵福老师(博客,知乎专栏,豆瓣)的「继续背单词,8个月过去了」,我就有写这篇文章的冲动了,杨老师说: 有时候我会感觉非常后悔,如 ...

  8. 【原】Learning Spark (Python版) 学习笔记(一)----RDD 基本概念与命令

    <Learning Spark>这本书算是Spark入门的必读书了,中文版是<Spark快速大数据分析>,不过豆瓣书评很有意思的是,英文原版评分7.4,评论都说入门而已深入不足 ...

  9. 【原】Learning Spark (Python版) 学习笔记(四)----Spark Sreaming与MLlib机器学习

    本来这篇是准备5.15更的,但是上周一直在忙签证和工作的事,没时间就推迟了,现在终于有时间来写写Learning Spark最后一部分内容了. 第10-11 章主要讲的是Spark Streaming ...

随机推荐

  1. angularjs集成requirejs

    其实说成使用requirejs加载angularjs应用会更贴切一些 <body> <span ng-controller="homeController"> ...

  2. Android自定义ImageView圆形头像

    效果图: 代码如下: RoundImageView.java import cn.comnav.evaluationsystem.R; import android.content.Context; ...

  3. Scalaz(31)- Free :自由数据结构-算式和算法的关注分离

    我们可以通过自由数据结构(Free Structure)实现对程序的算式和算法分离关注(separation of concern).算式(Abstract Syntax Tree, AST)即运算表 ...

  4. 微信公众平台关于fakeid和openid的解析

    今天在开发项目的时候遇到了个问题: 搞不清楚微信官方接口的fromusername(openid)和公众平台内每个粉丝所拥有的fakeid,于是在测试号中开始了对以上两项的研究,结果如下:   1.对 ...

  5. HAXM VT-X (与Hype-V冲突)

    之前一直使用vs emulator. 感觉性能各方面都比较好, 但在我更新完电脑后不知道什么原因各种起不来...  无奈之下想回到Google自带的模拟器. 然后发现intel haxm一直安装失败. ...

  6. 逐行读取txt文件

    header("Content-type:text/html; charset=utf-8"); $handle = fopen('test.php', 'r') or die(' ...

  7. App Store

    App store最新审核标准(2015.3)公布 http://www.cnii.com.cn/mobileinternet/2015-03/24/content_1550301.htm iOS提交 ...

  8. 动手写一个Remoting接口测试工具(附源码下载)

    基于.NET开发分布式系统,经常用到Remoting技术.在测试驱动开发流行的今天,如果针对分布式系统中的每个Remoting接口的每个方法都要写详细的测试脚本,无疑非常浪费时间.所以,我想写一个能自 ...

  9. [Tool] PowerDesigner

    一般项目的生命周期: 1.需求分析 2.需求规格说明书 3.总体设计 4.详细设计 5.编码实现 6.测试,试运行. 7.验收 8.后期维护 PowerDesigner 可以把软件生命周期的每一个阶段 ...

  10. 【原创】记一次Sql2008R2的数据库订阅发布遇到的问题!

    1.首先新建一个空的数据库 CreditTest,里面没有任何表结构,视图,数据等. 2.在已连接的服务器的下方找到“复制”选项,然后兴建发布: ...(新建发布成功) ...(新建订阅成功) 3.问 ...