Python中的地震SEGY数据操作工具
import segyio
三维速度sgy数据读取:
def get_segy_cube(fname,ilinenum,xlinenum):
if ilinenum *xlinenum ==0:
f = segyio.open(fname,"r")
else:
f = segyio.open(fname,"r", iline=ilinenum, xline=xlinenum)
v = f.trace.raw[:]
nil = len(f.iline)
nxl = len(f.xline)
nz = len(f.samples)
v.shape = (nil, nxl, nz)
f.close()
return v
默认字段为189和193,等价于使用77和73(注意顺序!)
复制道头: