试用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(sImageFile)
		
	def exportAsc(self,sAscFile):
		f2=open(sAscFile,"wt")
		sdataAsc=[str(x) for x in self.m_data]
		sDataAsc='\n'.join(sdataAsc)
		f2.writelines(sDataAsc)
		f2.close()
	
if __name__=="__main__":
	sacfile='d:\\1\\1.sac'
	sac=sacfile_wave()
	sac.read(sacfile)
	sac.draw("d:\\1\\1.png")
	sac.exportAsc("d:\\1\\1.asc")

绘图效果如图：