Res SU (13)读数据

1. 读.segy生成.su

# segyread tape=文件名.segy verbose=1 ns=采样点数 conv=1>文件名.su
segyread tape=data.sgy verbose=1 ns=50 conv=1>data.su

segyread tape=test.segy verbose=1 ns=100 conv=1>test.su

segyread tape=demo.sgy verbose=1 ns=100 conv=1>demo.su

segyread tape=LX_SEGY005.segy verbose=1 | segyclean > LX_SEGY005.su

2. 画剖面图

# 
suximage < data-tm-flt-dec-bal-cdp.su legend=1 d2=tr.d2 d1=tr.d1 title="test_su" cmap=rgb1 &


suximage < data-tm-flt-ac.su legend=1 d2=tr.d2 d1=tr.d1 title="test_su" cmap=rgb1 &

suximage < modeldata1 title="test_su" cmap=rgb1 &

3. 变密度

# suxwigb < LX_SEGY005.su d2=tr.d2 d1=tr.d1 title="LX_SEGY005"  &
suxwigb < data-tm-flt-dec-bal-cdp.su title="Test of suvlength" &

读数据

# 读取第一道
sustrip < data1.su | xxd -c 4 | sed 1q
# 读取第二道
sustrip < data-tm-flt-dec-bal-cdp.su | xxd -c 4 | sed 2q | tail -1

su文件格式的数据，道头信息如下：
n1=64 n2=32 d1=0.004000
nt=64 ntr=32 dt=0.004000
ns=64
00000000: 0000 0000  ....
现对该数据的第三道中第二个采样点数据进行修改，将其值改为1535。请使用seismic unix和Linux命令进行修改

读取数据存储到txt

• 这样存储的省略的数据

sustrip <  data-tm-flt-dec-bal-cdp.su | xxd -c 4 | sed 2q | tail -1 > data.txt

读取具体的su数据

sushw < data-tm-flt-dec-bal-cdp.su >output14.su key=tracl a=15 b=15
suascii <output14.su >output.txt

sushw < data-tm-flt-dec-bal-cdp.su >output14.su key=tracl,tracr a=24,24
suascii <output14.su >output.txt

• a=2和b=2表示提取第二道数据。

• 第二行命令使用suascii命令将输出SU文件转换为文本文件output.txt

# 使用sushw命令修改地震道集中第100个采样点的数值为新值（假设新值为10）
sushw < data.su key=tracl a=100 b=100 value=10 > modified_data.su