使用python3读IEEE格式segy文件并绘制变密度剖面图

# code for read IEEE segy format seismic data files 
import matplotlib.pyplot as mp
import numpy as np 
import sys
import struct
import binascii

fileName = "2267x1401.segy"
nTrace = 2207
nSample = 1401
fSegy = open(fileName,"rb")
data = np.zeros((nSample,nTrace))
fSegy.seek(3600,0)  
for iTrace in range(nTrace):
    fSegy.seek(240,1) 
    for iSample in range(nSample):
        tempValue = fSegy.read(4)
        hexValue = tempValue.hex()
        decValue = int(hexValue,16)
        floatValue = struct.unpack('<f',struct.pack('<I', decValue))[0]
        data[iSample][iTrace] = floatValue
fSegy.close()
print(data)
mp.matshow(data, cmap=mp.cm.gray)
mp.show()