python读取修改segy文件

读取和修改SEGY文件可以使用Python中的segyio库。首先需要安装segyio库，可以使用pip命令进行安装：

pip install segyio

读取SEGY文件的示例代码：

import segyio

filename = 'example.sgy'
with segyio.open(filename, 'r', ignore_geometry=True) as segyfile:
    # 获取文件头信息
    print(segyfile.text[0])
    # 获取道头信息
    print(segyfile.header[0])
    # 获取数据
    data = segyfile.trace.raw[:]

修改SEGY文件的示例代码：

import segyio

# 读取SEGY文件并进行修改
filename = 'example.sgy'
with segyio.open(filename, 'r+', ignore_geometry=True) as segyfile:
    # 修改文件头信息
    segyfile.text[0] = 'New file header information'
    # 修改道头信息
    for i in range(segyfile.tracecount):
        segyfile.header[i][segyio.TraceField.SourceX] += 100
        segyfile.header[i][segyio.TraceField.SourceY] += 100
        segyfile.header[i][segyio.TraceField.ReceiverX] += 100
        segyfile.header[i][segyio.TraceField.ReceiverY] += 100
    # 修改数据
    data = segyfile.trace.raw[:]
    data *= 2
    segyfile.trace.raw[:] = data