使用SCons自动执行Madagascar命令

Madagascar是著名的开源的地震资料处理工具包，基本可以代替古老的Seismic Unix用于学习和科研任务。Madagascar所有工具包都通过命令行界面交互，而完成一个研究项目，往往需要诸多步骤，执行大量命令。虽然也可以考虑将命令写入Shell脚本，但Madagascar提供了一种更加完善解决办法，那就是将操作写入SCons脚本，并使用SCons软件自动执行。

SCons is an Open Source software construction tool—that is, a next-generation build tool. Think of SCons as an improved, cross-platform substitute for the classic Make utility with integrated functionality similar to autoconf/automake and compiler caches such as ccache. In short, SCons is an easier, more reliable and faster way to build software.

执行SCons脚本十分简单，只要建立名为SConstruct的字符文件，在该文件所在的目录中执行scons命令，就能够让SCons代替人工自动执行处理操作了。

运行截图：windows上XTerminal的界面是原始了一点

基本语法

SCons使用Python作为其脚本语言，SConstruct脚本文件实际上相当于一个Python脚本，由于Python的基本语法简单易学，因而不需要专门学习Python编程语言也可以轻松使用，同时，Python能够实现复杂的逻辑，使得使用SCons脚本这种方式具有处理复杂任务的能力。

四种SCons脚本命令

为了让SCons执行Madagascar命令，专门为其设计了四个基本函数：

• Fetch （获取数据）
• Flow （执行处理）
• Plot （绘制图鉴）
• Result （输出图鉴）

这些函数定义了，文件与Madagascar命令的关系。

Flow

显然，Flow函数是其中最为核心的部分。其语法可以表示为：

Flow(output file, input file, command)

三个参数的类型都是字符串，最后的command字符中包含要执行的Madagascar命令和其相应的参数。
例如：

Flow("spike1","spike","scale dscale=4.0")

一个Flow命令就构建了输出文件、输入文件、所执行命令之间的关系，定义了一个操作步骤。相对应的Madagascar命令行为：

< spike.rsf sfscale dscale=4.0 > spike1.rsf

注意：Flow函数的参数中，文件名不需要加后缀

Plot和Result

Plot和Result两个函数负责执行绘图命令，它们的区别在于，Plot将生成的vpl文件保存在当前工作路径下，而Result则会将vpl文件保存在当前路径下的一个Fig文件夹中，以便使用Latex创建论文时使用。
它们的语法类似：

Plot(input file, command)

Result(input file, command)

这两个函数接受指定的输入文件，执行指定绘图命令，例如如果要生成一个数据文件的图鉴，可以使用：

Plot("spike1", "sfgraph pclip=100")

这样SCons就会知道我们要利用spike1.rsf数据文件使用sfgraph创建图鉴spike1.vpl

SConstruct基本结构

from rsf.proj import * # Remember, this statement comes first... ALWAYS
Flow("spike",None,"sfspike n1=100 k1=50") # None is a trick, see Advanced SCons for more information
Flow("spike1","spike","sfadd scale=4.0")
Flow("noise","spike1","sfnoise")
Plot("spike",'sfgraph title="spike" ' ) # Note string nesting
Plot("spike1",'sfgraph title="spike1" ') 
Plot("noise",'sfgraph title="noisy" ')
Result("noise",'sfgraph title="noisy" pclip=75 ')
End() # Remember, this always ends the script.