试用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,sImag