python 使用segyio 读取格式为segy的地震数据 以及解决报错unable to find sorting #322的办法

使用python读取segy格式的地震数据

首先我们需要安装库文件segyio，以及需要的其他库文件numpy(安装方法详见置顶文章，配置pip环境变量以及使用pip安装库文件)

    
import os
import segyio
import numpy as np

filename = 'testsgy.segy'
with segyio.open(filename) as segyfile:

    # Memory map file for faster reading (especially if file is big...)
    segyfile.mmap()

    # Print binary header info
    print(segyfile.bin)
    print(segyfile.bin[segyio.BinField.Traces])

    # Read headerword inline for trace 10
    print(segyfile.header[10][segyio.TraceField.INLINE_3D])

    # Print inline and crossline axis
    print(segyfile.xlines)
    print(segyfile.ilines)

有的同学可能会出现如下错误：

unable to find sorting #322的办法

解决办法其实很简单：你安装的segyio 版本可能高于1.6.0版本，将高版本卸载掉，安装1.6.0版本的segyio库文件，完美解决！