seismic unix 输入txt文件

1、准备txt文件

复制数据到文本编辑器，保存为input.txt文件

9.5012929e-01
6.1543235e-01 
5.7891305e-02
1.5273927e-02
8.3811845e-01
1.9343116e-01
4.9655245e-01
7.2711322e-01
7.9482108e-01
1.3651874e-01
5.8279168e-01
2.0906940e-01 
4.1537486e-01
2.1396333e-01
6.8333232e-01
2.3113851e-01
7.9193704e-01
3.5286813e-01
7.4678568e-01
1.9639514e-02
6.8222322e-01
8.9976918e-01
3.0929016e-01
9.5684345e-01
1.1756687e-02
4.2349626e-01
3.7981837e-01
3.0499868e-01
6.4349229e-01
2.1255986e-01
6.0684258e-01
9.2181297e-01
8.1316650e-01
4.4509643e-01
6.8127716e-01
3.0276440e-01
8.2162916e-01
8.3849604e-01
5.2259035e-01
8.9389797e-01
5.1551175e-01
7.8332865e-01
8.7436717e-01
3.2003558e-01
8.3923824e-01
4.8598247e-01
7.3820725e-01
9.8613007e-03
9.3181458e-01
3.7948102e-01
5.4167385e-01
6.4491038e-01
5.6807246e-01
8.8014221e-01
1.9913807e-01
3.3395148e-01
6.8084575e-01
1.5009499e-02
9.6009860e-01
6.2878460e-01 
8.9129897e-01
1.7626614e-01
1.3889088e-01
4.6599434e-01
8.3179602e-01
1.5087298e-01
8.1797434e-01
3.7041356e-01
1.7295614e-01
2.9872301e-01
4.3290660e-01
4.6109513e-01
7.6795039e-01
7.2663177e-01
1.3377275e-01
7.6209683e-01
4.0570621e-01
2.0276522e-01
4.1864947e-01
5.0281288e-01
6.9789848e-01
6.6022756e-01
7.0273991e-01
9.7974690e-01
6.6144258e-01
2.2594987e-01
5.6782871e-01
9.7084494e-01
4.1195321e-01
2.0713273e-01
4.5646767e-01
9.3546970e-01
1.9872174e-01
8.4622142e-01
7.0947139e-01
3.7837300e-01
3.4197062e-01
5.4657115e-01
93.9008259e-01
2.7144726e-01
2.8440859e-01
5.7980687e-01
7.9421065e-01
9.9008259e-01
7.4456578e-01
6.0719894e-01
1.8503643e-02
9.1690444e-01
6.0379248e-01
5.2515250e-01
4.2889237e-01
8.6001160e-01
2.8972590e-01
4.4488020e-01
2.5232935e-01
4.6922429e-01
7.6036501e-01
5.9182593e-02
7.8886169e-01
2.6794725e-01
6.2988785e-01
8.2140716e-01
4.1027021e-01
2.7218792e-01
2.0264736e-01
3.0461737e-01
8.5365513e-01
3.4119357e-01
6.9456724e-01
8.7574190e-01
6.4781123e-02
5.2982312e-01
6.0286909e-01
4.3865853e-01
4.3992431e-01
3.7047683e-01
4.4470336e-01
8.9364953e-01
1.9881427e-01
6.7213747e-01
1.8965375e-01
5.9356291e-01
5.3407902e-01
6.2131013e-01
7.3730599e-01
9.8833494e-01
63.9008259e-01
6.4052650e-01
5.0268804e-02
4.9831130e-01
9.3338011e-01
5.7514778e-01
48.7837300e-01
51.7974690e-01 
75.9655245e-01
62.8222322e-01
49.0273991e-01
63.9008259e-01
68.3651874e-01
54.9182593e-02
54.3546970e-01
69.5317

2、转换txt文件为二进式文件

a2b<input.txt>input.bin n1=8 &
读入txt文件，转换为二进制文件，其中转换后每行8个数字

3、准备su的头文件

sustrip<s01.su head=data.headers >s01.bin
在现有文件s01.su中分离头文件和二进制进制数据，头文件保存为data.headers，数据保存为s01.bin，头文件用于与我们自制的何必，数据文件直接丢弃：）

4、生成su文件

supaste<input.bin head=data.headers>input.su *ns=1024
suxwigb<input.su &

4、修改道头信息

surange <input.su&
suchw<input.su key1=tracl a=1000>input_chw.su &
将tracl值改编为1000

![请添加图片描述](https://img-blog.csdnimg.cn/223d59cb3f2145b0a1ddc1ab19ecac72.png

尝试同时修改tracl 和tracr出现错误

% suchw<input.su key1=tracl key2=tracr a=1200 >input_chw.su

本意尝试同时修改tracl和tracr两个参数，但结果发现只有tracl由5变为1205，继续寻找原因。

找到suchw的说明：
有些道头字（如cdp）可以从已有的道头字计算而来，程序suchw就提供了这种功能。
参数有：key1 输出的关键字
key2 输入的关键字
key3 输入的关键字

    a   偏移量
    b   key2关键字的倍数

c key3关键字的倍数
d overall scales
我们可以使用两个道头字的值（key2 和key3），利用下面的方程式计算第三个道头字的值(key1):
val(key1) = (a + b * val(key2) + c * val(key3)) / d
例如：
suchwoutdata

现在的问题时a，b，c，d分别都是多少呢？回到su中的说明文档，找到我们需要的默认值：
Optional parameters:
key1=cdp,… output key(s)
key2=cdp,… input key(s)
key3=cdp,… input key(s)
a=0,… overall shift(s)
b=1,… scale(s) on first input key(s)
c=0,… scale on second input key(s)
d=1,… overall scale(s)
e=1,… exponent on first input key(s)
f=1,… exponent on second input key(s)

The value of header word key1 is computed from the values of
key2 and key3 by:

    val(key1) = (a + b * val(key2)^e + c * val(key3)^f) / d    

由说明可以看出b默认值时1，因此a=1200,val(key2)=5,val(key1)=a+b*val(key2)=1200+5=1205，是不是没毛病？没毛病
这也就就解释了为什么数值跟我们预期的不同。

https://download.csdn.net/download/weixin_56694865/87399529