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

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

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

    
import os
import segyio
import numpy as npfilename = 'testsgy.segy'
with segyio.open(filename) as segyfile:# Memory map file for faster reading (especially if file is big...)segyfile.mmap()# Print binary header infoprint(segyfile.bin)print(segyfile.bin[segyio.BinField.Traces])# Read headerword inline for trace 10print(segyfile.header[10][segyio.TraceField.INLINE_3D])# Print inline and crossline axisprint(segyfile.xlines)print(segyfile.ilines)

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

unable to find sorting #322的办法

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