readsegy的python代码

以下是一个简单的示例代码，用于读取SEGY文件的数据：

import segyio

def read_segy(file_path):
    with segyio.open(file_path, "r") as segyfile:
        # 获取文件的相关信息
        print("SEGY file info:")
        print(f" - File format: {segyfile.format}")
        print(f" - Sample count: {segyfile.samples.size}")
        print(f" - Trace count: {segyfile.tracecount}")
        print(f" - Inline count: {segyfile.ilines.size}")
        print(f" - Crossline count: {segyfile.xlines.size}")

        # 读取数据
        traces = segyfile.trace.raw[:]
        return traces

# 用法示例
segy_file_path = "example.segy"
traces = read_segy(segy_file_path)

这段代码使用了`segyio`库来读取SEGY文件。首先，我们通过`segyio.open()`方法打开SEGY文件，并获取了一些相关的信息，如文件格式、样本数、道数、内联数和横断面数。然后，我们使用`segyfile.trace.raw[:]`来读取所有的道数据，并将其存储在一个变量中（这里使用`traces`）。你可以根据自己的需求进行进一步处理和分析。