python处理地震sac数据_试用python读取地震波形的二进制SAC文件

sac文件比较简单，不过有些设置没有处理，只是单纯的读取并绘图。

代码：

import struct

import pylab

class sacfile_wave:

def read(self,sFile):

#sFile='e:\\10.sac'

f=open(sFile,'rb')

hdrBin=f.read(632)

sfmt='f'*70+'I '*40+'8s '*22+'16s';

hdrFmt=struct.Struct(sfmt)

self.m_header=hdrFmt.unpack(hdrBin)

npts=int(self.m_header[79])

fmt_data='f'*npts

dataFmt=struct.Struct(fmt_data)

dataBin=f.read(4*npts)

f.close()

self.m_data=dataFmt.unpack(dataBin)

print "data len:",len(self.m_data)

def draw(self,sImageFile):

npts=len(self.m_data)

xd=range(1,npts+1)

pylab.figure(1)

pylab.plot(xd,self.m_data,linewidth=0.3)

pylab.savefig(sImageFi