SU学习笔记（二）

四、编辑SU数据

一旦SU格式的数据读入和道头设置正确后，常常要对数据进行操作和编辑。SU提供的数据处理和编辑的命令有：

SUWIND – 根据关键字开时窗

SUSORT – 基于segy道头关键字进行道分选

SURAMP – 从道起点到道终点进行线性斜坡化(Linearly taper)至零值

SUTAPER – 对一个数据窗进行斜坡化(taper)至零值

SUNULL – 产生空道（都为零值）

SUZERO – 在一时间窗内数据充零

SUKILL – 地震道充零

SUNUTE – 根据关键道头字确定的距曲线的距离，去除用户指定的多边形曲线上边或下边的值

SUVELENGTH – 调整变长度地震道到相同长度地震道

SUVCAT – 将一个数据附加到另一个数据（一道接一道进行）

1.        suwind – 用关键字选定一定时窗地震道

非常常见的操作就是我们时常需要浏览或处理地震数据中的一小块数据。suwind命令使得我们可以根据大量的参数信息设置窗口来选定感兴趣的数据。

通过道头字来设定数据窗口的大小

suwind最简单的用法就是，通过用户设定道头关键字的最小和最大值来选定数据：

key=tracl           设置窗口的道头字(see segy.h)

min=LONG_MIN    数据关键道头字的最小值

max=LONG_MAX   数据关键道头字的最大值

例如，用道数作关键字对命令suplane产生的数据进行窗口大小设置：

suplane | suwind key=tracl min=5 max=10 | sugethw key=tracl | more

对于一个大的数据体，应该使用计数（count）参数，而不是最大值。如果直接设定最大值（max），命令suwind将查询全部数据后，再选择在最小值和最大值之间的道集，这是因为程序认为道的标志（trace labeling）可能多次出现。例如，可以比较下面的两个命令的结果：

suplane ntr=100000 | suwind key=tracl min=5 max=10 | sugethw tracl | more

suplane ntr=100000 | suwind key=tracl min=5 count=5 | sugethw tracl | more

更复杂的窗口大小选择如下(例如，使用分数，decimating data)：

j=1     每个数据的第j道不选…

s=0     以s为基数(如果((key - s)%j) == 0)

下面的例子用suplane生成的数据，每两道抽取一道（即抽取道2、4、6…）：

suplane | suwind key=tracl j=2 | sugethw key=tracl | more

或者每隔2道抽取1道，基数为1（即抽取道1、3、5…）：

suplane | suwind key=tracl j=2 s=1 | sugethw key=tracl | more

也可以用suwind接受或拒绝某些地震道：

reject=none      按指定的关键字跳过的地震道

accept=none      用指定的关键字选择地震道

例如下面的命令所示，地震道中的第3、8、9道将不显示：

suplane | suwind key=tracl reject=3,8,9 | sugethw key=tracl | more

参数accept选项有点特殊，它表示接收这些道，即使这些道被拒绝过。例如：

suplane | suwind key=tracl reject=3,8,9 accept=8 | sugethw key=tracl | more

如果你只想接收列表中的道，这时需要设置参数“max=0”,下列所示只显示出第8道:

suplane | suwind key=tracl accept=8 max=0 | sugethw key=tracl | more

计数（count）参数覆盖接收（accept）参数，所以如果你想真正无条件接受道，就不能指定计数(count)参数。

选取时间窗口（Time gating）

有关窗口的第二个问题就是时间窗的设定。垂直时间窗（time gating）的选项为：

tmin = 0.0                选取最小时间

tmax = (from header)    选取最大时间

itmin = 0                 选取最小时间采样

itmax = (from header)    选取最大时间采样

nt = itmax-itmin+1        时间采样点数

设置的结果要么是选取最小时间采样和选取最大时间采样（itmin和itmax），要么选取最小时间和选取最大时间（tmin和tmax），这样就将产生一个最接近采样的时窗。

2.        susort– 根据segy道头字分选道

susort命令利用Unix系统sort命令，通过道头关键字对地震道进行分选。

例如，根据两个关键字(cdp和offset)分选数据（数值按升序），可用下面的命令：

susort <indata.su >outdata.su cdp offset

如果offset关键字按降序，cdp关键字按升序进行分选，可用下面的命令：

susort <indata.su >outdata.su cdp -offset

3.        suramp和sutaper– 数据斜坡化（Tapering）

很多地震处理算法因为数据的突变边界而显示假象，在数据边界对振幅进行斜坡化(tapering)，就是压制这种假象的最容易的方法。在SU中，我们可以使用sutaper命令对数据边界进行斜坡化。

例如，对每5道数据从头到尾进行斜坡化：

sutaper <diskfile >stdout ntaper=5

suramp可以平滑从起始道和(或)末道。例如，斜坡化从0到tmin=0.05s,然后再向下从tmax=1.15s到末道：

suramp <diskfile tmin=.05 tmax=1.15 >stdout

4.        sukill、suzero、sunull和sumute– 数据充零

对有噪声的道、位于数据边界的道进行充零是有益的。或者是产生一些空道，在绘图时最为联系数据集中间的分割标志。

sukill – 道充零

对一批地震道充零键入：

sukill <stdin >stdout min=MIN_TRACE count=COUNT

这里参数count是要充零的总道数，参数min是要充零的这些道中的最小道号。

sunull – 产生空道数据

有时需要生成道值为零的地震数据块。下面的命令产生共有NTR道，时间采样数为NT的地震数据：

sunull nt=NT ntr=NTR <stdin >stdout min=MIN_TRACE count=COUNT

suzero – 在一段时间窗内对数据充零

suzero itmin=MIN_TIME_SAMPLE itmax=MAX_TIME_SAMPLE <indata.su > outdata.su

sumute – 数据去除

要进行高精度的压制操作，可用命令sumute对SU数据进行处理。

sumute <indata.su >outdata.su key=KEYWORD xmute=x1,x2,x3,... tmute=t1,t2,t3,...

下面用suplane生成数据，然后做去除处理与原始生成数据比较：

suplane | suxwigb &

suplane | sumute key=tracl xmute=1,10,12 tmute=.06,.1,.11 | suxwigb &

去除根据xmute=和tmute=参数确定的多边形曲线内每一个初至。

5.        suvcat和cat – 数据合并

有两种方法可以将一种数据附加到另一个上（合并）。第一种使用Unix命令cat，简单的将第二个文件的数据放到第一个文件中去。

cat data1.su data2.su > data3.su

另外，可能需要对地震道重编号：

cat data1.su data2.su | sushw key=tracl a=1 > data3.su

第二种是将第二个数据集中的每一道垂直的附加到第一个数据集中每一道的末尾。这就需要使用suvcat命令：

suvcat data1.su data2.su > data3.su

在这个例子中，就不需要修改道头字。

6.        suvlength – 调整可变长度道到相同采样点数

有时数据中每道含有不同的采样点数。下面我们对用suplane命令生成的数据，用命令suvlength处理的道相同的采样点数：

suplane nt=64 > data1.su

suplane nt=32 > data2.su

cat data1.su data2.su > data3.su

对于上面合并生成的数据data3.su，如果想用SU程序处理该数据文件将会失败，因为大多说SU程序要求数据块含有相同的采样点数。使用命令suvlength可以解决该问题：

suvlength ns=64 < data3.su > data4.su

suxwigb < data4.su title="Test of suvlength" &

这将使所有的道含有相同的长度。

五、SU数据常用操作

1.    suaddnoise – 对地震道加噪声

下面是两个使用suaddnoise命令的例子：

suplane | suxwigb title="no noise" &

suplane | suaddnoise | suxwigb title="noise added" &

suplane | suaddnoise sn=2 | suxwigb title="noise added" &

2.    sugain – su数据增益

增益命令含有多个选项：

scaling the data, 数据比例伸缩

multiplying the data by a power of time, 数据与时间幂相乘

taking the power of the data, 数据取幂

automatic gain control, 自动增益控制

trapping noise spiked traces, 含噪声脉冲的地震道陷波滤波

clipping specified amplitudes or quantiles, 指定的振幅或分量裁减限制

balancing traces by quantile clip, rms value, or mean, 根据等分法、均方根值或均值道均衡

biasing or debiasing the data. 偏离或去偏离数据

操作优先等级如下面的方程所示：

out(t) = scale * BAL{CLIP[AGC{[t^tpow * exp(epow * t) * ( in(t)-bias )]^gpow}]}

实例：

suplane | suaddnoise > data.su

suxwigb < data.su title="Ungained Data" &

sugain < data.su scale=5.0 | suxwigb title="Scaled data" &

sugain < data.su agc=1 wagc=.01 | suxwigb title="AGC=1 WAGC=.01 sec &

sugain < data.su agc=1 wagc=.2 | suxwigb title="AGC=1 WAGC=.1 sec &

sugain < data.su pbal=1 | suxwigb title="traces balanced by rms" &

sugain < data.su qbal=1 | suxwigb title="traces balanced by quantile" &

sugain < data.su mbal=1 | suxwigb title="traces balanced by mean" &

sugain < data.su tpow=2 | suxwigb title="t squared factor applied" &

sugain < data.su tpow=.5 | suxwigb title="square root t factor applied" &

3.    suop – su数据的一元操作

suop对SU数据进行一元函数运算和操作，主要功能包括：

absolute value,              取绝对值

signed square root,           带符号开方

square,                     平方

signed square,               带符号位平方

signum function,             正负号函数

exponential,                 指数

natural logarithm,             取自然对数

signed common logarithm,     带符号位常用对数

cosine,                     余弦函数

sine,                       正弦函数

tangent,                     正切函数

hyperbolic cosine,            双曲余弦函数

hyperbolic sine,              双曲正弦函数

hyperbolic tangent,           双曲正切函数

divide trace by Max. Value,     地震道数据用最大值除

express trace values in decibels: 20 * slog10 (data)  地震道数据用分贝表示

negate values,                数据取反

pass only positive values,       只选正值

pass only negative values.       只选负值

例子：

suplane | suaddnoise > data.su

suop < data.su op=abs | suxwigb title="absolute value" &

suop < data.su op=ssqrt | suxwigb title="signed square root" &

suop < data.su op=sqr | suxwigb title="signed square" &

4.    suop2 – su数据的二元操作

程序suop2用来对两个SU数据进行操作，该命令支持的计算有：

difference,     相减

sum,           相加

product,       相乘

quotient,       相除

前面的选项假定每个SU数据道采样数相同，后面4个选项假定第二个文件只有一道。

difference of a panel and a single trace,  一块数据和某一道数据相减

sum of a panel and a single trace,        一块数据和某一道数据相加

product of a panel and a single trace,    一块数据和某一道数据相乘

quotient of a panel and a single trace.   一块数据和某一道数据相除

并且有8个对等shell文本命令做这些操作：

susum file1 file2 == suop2 file1 file2 op=sum

sudiff file1 file2 == suop2 file1 file2 op=diff

suprod file1 file2 == suop2 file1 file2 op=prod

suquo file1 file2 == suop2 file1 file2 op=quo

For: panel "op" trace operations:

suptsum file1 file2 == suop2 file1 file2 op=ptsum

suptdiff file1 file2 == suop2 file1 file2 op=ptdiff

suptprod file1 file2 == suop2 file1 file2 op=ptprod

suptquo file1 file2 == suop2 file1 file2 op=ptquo

所有这些操作都是调用suop2进行计算。

例子：

suplane > junk1.su

suxwigb < junk1.su | suxwigb title="Data without noise" &

suplane | suaddnoise > junk2.su

suxwigb < junk2.su | suxwigb title="Data with noise added" &

suop2 junk2.su junk1.su op=diff | suxwigb title="difference" &

注意，文件名应出现在操作参数”op=”之前。

六、SU信号处理

1.        Fourier变换

在SU软件包中有1D和2D Fourier变换程序。2D变换包括地震F-K变换，假定输入数据的块空间是时间的，第二维是空间的；对非地震K1-K2变换，假定输入是纯二维空间（x1, x2）数据。

1D Fourier变换

SUFFT – 从实数时间道变换到复数频率道（正变换）

SUIFFT – 从复数频率道变换到实数时间道（逆变换）

SUAMP – 输出振幅、相位、实部和虚部（从frequency, x）

SUSPECFX – 地震数据Fourier频谱分析（时间T到频率F）

例如：

suplane | suxwigb title="Original Data" &

suplane | sufft | suifft | sushw key=d1,dt a=0,4000 | suxwigb &

结果与输入完全一样，除了地震道结果中有更多的采样点，因为变换时要补零。

要浏览sufft程序输出的振幅谱和相位谱、实部和虚部，操作如下：

suplane | sufft | suamp mode=amp | suxwigb title="amplitude" &

suplane | sufft | suamp mode=phase | suxwigb title="phases" &

suplane | sufft | suamp mode=real | suxwigb title="real" &

suplane | sufft | suamp mode=imag | suxwigb title="imaginary" &

SU数据格式可以存储复数数据的实部和虚部。键入下面的命令可以查看格式的道头设置位置：

suplane | sufft | surange

显示出道头字trid=11，显示出FFT变换输出数据是如何排列的.

当然，绝大多数时候，我们只想快速看看一个地震道或一块地震数据的振幅谱，这时可以用suspecfx命令：

suplane | suspecfx | suximage title="F-X Amplitude Spectrum" &

这将直接显示输入SU数据每道的振幅谱。

2D Fourier变换

如果数据实际上为（时间、空间）坐标数据，那么2D Fourier变换输出的就是F-K(频率、波数)域数据。如果考虑两个空间坐标（x1,x2）数据，进行2D Fourier变换，这时输出就是（k1,k2）2D波数域数据。

SUSPECFK – 数据的F-K域Fourier频谱分析

SUSPECK1K2 - （x1,x2）数据2D（K1,K2）Fourier频谱分析

例子：

suplane | suspecfk | suximage title="F-K Amplitude Spectrum" &

suplane | suspeck1k2 | suximage title="K1-K2 Amplitude Spectrum" &

2.        Hilbert变换、道属性和时间-频率域

SUHILB – Hilbert变换

SUATTRIBUTES – 道属性（瞬时振幅、相位和频率）

SUGABOR – 通过Gabor变换（类似多参数滤波分析），输出地震数据时间-频率结果

进行Hilbert变换的例子如下：

suplane | suhilb | suxwigb title="Hilbert Transform" &

查看suattributes命令生成道属性，将相应的显示出瞬时振幅、相位和频率：

suvibro | suxgraph title="Vibroseis sweep" &

suvibro | suattributes mode=amp | suxgraph title="Inst. amplitude" &

suvibro | suattributes mode=phase unwrap=1.0 | suxgraph title="Inst. phase" &

suvibro | suattributes mode=freq | suxgraph title="Inst. frequency" &

使用sugabor命令可以进行时间-频率域分析：

suvibro | sugabor | suximage title="time frequency plot" &

3.        Radon变换 – Tau P滤波

Radon变换或“Tau_P”变换时压制多次波和进行数据改造的有效方法，涉及的程序有:

SUTAUP - 正向和逆向（forwared and inverse）T-X和F-K全局倾斜叠加

SUHARLAN - 信号噪声分离，应用Harlan可逆线性变换方法（1984）

SURADON - 计算正或逆Radon变换，或用抛物线Radon变换估计并去除多次波

SUINTERP - 使用自动同相轴拾取方法进行道插值

相关使用方法，请参阅：$CWPROOT/src/demos/Tau P下的例子

4.        1D滤波操作

地震处理中很大一部分工作可以称为滤波。在SU中1D滤波应用程序可以进行简单的滤波操作，和更复杂的反褶积和子波整形处理。这些操作成为1D，是因为这些操作处理数据是一道一道进行的。

SU包含的命令有：

SUFILTER - 应用零相位正弦平方斜坡滤波

SUBFILT - 应用Butterworth带通滤波

SUACOR - 自相关

SUCONV, SUXCOR - 利用用户提供的滤波器进行褶积和相关

SUPEF - Wiener预测误差滤波

SUSHAPE - Wiener整形滤波

SURESAMP - 时间域重采样

SUFRAC - 采用一般的时间微分、积分、相位移动，注：时间域数据

SUMEDIAN - 沿着用户确定的多边形曲线，根据道头字指定的与曲线的距离，进行中值滤波

SUTVBAND - 时变带通滤波（正弦平方滤波）

1）sufilter - 应用零相位正弦平方斜坡滤波

sufiletr提供了可用于一般的带通、带陷、低通、高通和陷波滤波，使用sufilter进行滤波的例子如下：

suplane | sufilter f=10,20,30,60 amps=0,1,1,0 | suxwigb title="10,20,30,60 hz bandpass" &

suplane | sufilter f=10,20,30,60 amps=1,0,0,1 | suxwigb title="10,20,30,60 hz bandreject" &

suplane | sufilter f=10,20,30,60 amps=1,1,0,0 | suxwigb title="10,20 hz lowpass" &

suplane | sufilter f=50,60,70 amps=1,0,1 | suxwigb title="60 hz notch" &

2）subfilt - 应用Butterworth带通滤波

subfilt与程序sufilter工程相似，它使用Butterworth滤波器进行滤波：

suplane | subfilt fstoplo=10 fpasslo=20 fpasshi=30 fstophi=60 | suxwigb title="10,20,30,60 hz bandpass” &

3）suconv, suxcor - 使用用户设计的滤波器进行褶积和相关

褶积和互相关可分别使用suconv和suxcor命令实现。滤波器可以通过在命令行输入向量、或者通过含有某个道的SU格式文件提供。

下面是一个可控联系扫描信号进行相关处理的例子。通过联合使用suvibro、suplane和suconv，产生类似可控连续震动的数据。生成可控连续震动的命令如下：

suvibro > junk.vib.su

suplane | suconv sufile=junk.vib.su > plane.vib.su

使用surange命令查看junk.vib.su数据的信息：

surange < junk.vib.su

显示可控震源扫描信号有2500个采样点，因此做下面所示的相关分析：

suxcor < plane.vib.su sufile=junk.vib.su | suwind itmin=2500 itmax=2563 | sushw key=delrt a=0.0 > data.su

参数itmin=sweeplength和itmax=sweeplength+nsout，参数nsout是期望输出的采样点数，最后使用sushw设置道延迟为0。选择参数itmin=sweeplength将保证数据从正确的值起始。

4）supef - Wiener预测误差滤波

5）sushape - Wiener整形滤波

5.        2D滤波操作

二维空间（k1,k2）和（F,K）域滤波有助于改变数据的倾角信息。下面的这些程序初步的提供了一套波数(K)域和频率-波数(F，K)域滤波操作命令：

SUKFILTER - radially symmetric K-domain, sin^2-tapered, polygonal filter

SUK1K2FILTER - symmetric box-like K-domain filter defined by the cartesian product of two sin^2-tapered polygonal filters defined in k1 and k2

SUKFRAC - apply FRACtional powers of i-k to data, with phase shift

SUDIPFILT - DIP-or better-SLOPE Filter in f-k domain

6.        suresamp - 数据时间域重采样

七、SU数据绘图及显示

Seismic Unix软件包包括一部分（普通C语言类型的浮点格式和SU格式）数据绘图的图形工具，均可在X-Windows环境下，在屏幕上显示或生成PostScript格式的图形文件。SU中能够绘制的图形类型有：

contour plots, 等值线绘图

gray or colorscale image plots, 灰度或彩色图像

wiggle trace plots, 地震道波形图

line or symbol graphs, 测线和符号绘图

movies, 电影动画

3D cube plots (PostScript only). 3D立体图

1.    X-Windows中绘制一般浮点型数据

在X-Windows环境下，用于绘制浮点型数据（没有SU道头）的程序有：

XCONTOUR - 调用矢量绘图来绘制f(x1,x2)的等值线图

XIMAGE - 绘制灰度或彩色图像

XWIGB - 绘制变面积地震道波形图

XGRAPH - 绘制曲线

XMOVIE - 绘制动画电影

例子：

suplane | sustrip > data.bin

xcontour < data.bin n1=64 n2=32 title="contour" &

ximage < data.bin n1=64 n2=32 title="image" &

xwigb < data.bin n1=64 n2=32 title="wiggle trace" &

xmovie < data.bin n1=64 n2=32 title="movie" &

上面命令中n1=64,n2=32表示数据大小为：46*32。

xgraph < data.b n=5

完成上述操作后，要消去窗口，可以点击窗口然后键入“q”即可退出。

2.    X-Windows中SU格式数据绘图

相应，在X-Windows环境下，用于绘制SU格式数据的程序有：

SUXCONTOUR - 调用矢量绘图来绘制f(x1,x2)的等值线图

SUXIMAGE - 绘制灰度或彩色图像

SUXWIGB - 绘制变面积地震道波形图

SUXGRAPH - 绘制曲线

SUXMOVIE - 绘制动画电影

SUXMAX - SU数据每道最大值、最小值和最大绝对值X-Windwos图

例子：

suplane | suxcontour title="contour" &

suplane | suximage title="image" &

suplane | suxwigb title="wiggle trace" &

suplane | suxgraph title="graph" &

suplane | suxmovie title="movie" &

suplane | suxmax title="max" &

绘制地震道时，使用真偏移距，需要使用“key=”关键字参数。例如：

suplane | suchw key1=offset key2=tracl a=0 b=100 | suxwigb key=offset &

用suxmovie命令制作电影动画：

suplane > junk1.su

suplane | suaddnoise sn=20 >> junk1.su

suplane | suaddnoise sn=15 >> junk1.su

suplane | suaddnoise sn=10 >> junk1.su

suplane | suaddnoise sn=5 >> junk1.su

suplane | suaddnoise sn=3 >> junk1.su

suplane | suaddnoise sn=2 >> junk1.su

suplane | suaddnoise sn=1 >> junk1.su

suxmovie < junk1.su n2=32 title="frame=%g" loop=1 &

这里参数n2=32表示数据每幅有32道。“%g”用来在图形标题上显示出图形的桢数，参数“loop=1”表示用连续循环方式运行电影。

要加快或减慢电影的放映速度，只需要在图形的右下角点击并拖动窗口来放大或缩小图像。按在鼠标最右边的键一次将暂停，按第二次将重新放映。

3.    一般浮点型数据的PostScript绘图

可用于一般浮点型数据的PostScript绘图的命令有：

PSCONTOUR - PostScript CONTOURing of a two-dimensional function f(x1,x2),

PSIMAGE - PostScript IMAGE plot of a uniformly-sampled function f(x1,x2),

PSCUBE - PostScript image plot of a data CUBE,

PSGRAPH - PostScript GRAPHer Graphs n[i] pairs of (x,y) coordinates,

PSMOVIE - PostScript MOVIE plot of a uniformly-sampled function f(x1,x2,x3),

PSWIGB - PostScript WIGgle-trace plot of f(x1,x2) via Bitmap, 用位图绘制地震道

PSWIGP - PSWIGP - PostScript WIGgle-trace plot of f(x1,x2) via Polygons. 用多边形绘制地震道

例子：

suplane | sustrip > data.bin

pscontour < data.bin n1=64 n2=32 title="contour" > data1.eps

psimage < data.bin n1=64 n2=32 title="image" > data2.eps

pscube < data.bin n1=64 n2=32 title="cube plot" > data4.eps

pswigb < data.bin n1=64 n2=32 title="bitmap wiggle trace" > data3.eps

pswigp < data.bin n1=64 n2=32 title="wiggle trace" > data4.eps

psmovie < data.bin n1=64 n2=32 title="movie" > data5.eps

a2b < data.ascii n1=2 > data.bin

n=5

psgraph < data.bin n=5 > data6.eps

4.    SU数据的PostScript绘图

相应，可用于SU数据的PostScript绘图命令有：

SUPSCONTOUR - PostScript CONTOUR plot of an SU data set

SUPSIMAGE - PostScript IMAGE plot of an SU data set

SUPSCUBE - PostScript CUBE plot of an SU data set

SUPSGRAPH - PostScript GRAPH plot of an SU data set

SUPSWIGB - PostScript Bit-mapped WIGgle plot of an SU data set

SUPSWIGP - PostScript Polygon-filled WIGgle plot of an SU data set

SUPSMAX - PostScript of the MAX, min, or absolute max value on each trace of a SU data set

例子：

suplane > junk.su

supscontour < junk.su title="contour" > data1.eps

supsimage < junk.su title="image" label1="sec" label2="trace number" > data2.eps

supscube < junk.su title="cube plot" > data4.eps

supswigb < junk.su title="bitmap wiggle trace" > data3.eps

supswigp < junk.su title="wiggle trace" > data4.eps

supsmovie < junk.su title="movie" > data5.eps

supsmax < junk.su title="max" > data5.eps

5.    其它PostScript绘图工具

PSBBOX - change BoundingBOX of existing PostScript file 改变PostScript文件四周的框架

PSMERGE - MERGE PostScript files 合并PostScript文件

MERGE2 - MERGE2 PostScript figures onto one page 合并两幅图像

MERGE4 - MERGE4 figures onto one page 合并四幅图像

PSLABEL - output PostScript file consisting of a single TEXT string on a specified background.(Use with psmerge to label plots.)

PSMANAGER - printer MANAGER for HP 4MV and HP 5Si Mx Laserjet PostScript printing

PSEPSI - add an EPSI formatted preview bitmap to an EPS file

例子：

suplane > junk.su

suplane | sufilter > junk1.su

supswigb < junk.su title="Wiggle trace" label1="sec" label2="trace number" > junk1.eps

supsimage < junk.su title="Image Plot" label1="sec" label2="trace number" > junk2.eps

supscontour < junk.su title="Contour Plot" label1="sec" label2="trace number" > junk3.eps

supswigb < junk1.su title="Filtered" label1="sec" label2="trace number" > junk4.eps

合并两幅图到一个文件中去：

merge2 junk1.eps junk2.eps > junk.m2.eps

合并4幅图到一个文件中去：

merge2 junk1.eps junk2.eps junk3.eps junk4.eps > junk.m4.eps

6.    道拾取工具(Trace Picking Utilities)

X-Windows地震道、图像和等值线绘图程序中都具有如下特点：将光标放在图形上某点，然后键入字母“s”，那么该点的坐标就保存在内存中；当键入字母“q”时，这些值就存储到用户指定的一个文件“mpicks”中。

SU中另外还有两个程序用于拾取：

XPICKER - X wiggle-trace plot of f(x1,x2) via Bitmap with PICKing

SUXPICKER - X-windows WIGgle plot PICKER of an SU data set

SUPICKAMP - pick amplitudes within user de_ned and resampled window