Readers For Laboratory Acoustic Data

ReadBinSeismo.m

function[seismo,size] = ReadBinSeismo(fname)

fid = fopen(fname,'rb');
size = fread(fid,2,'int');
size(2)
seismo = fread(fid,size(2),'double')
fclose(fid)

SeisRead.py

import numpy as np

def ReadBinSeismo(fname):
    
        """
        Takes filename and returns the seismogram data and 
        size of the seismogram.  JRL 8/1/13
        """
    
        f = open(fname,'rb')
        seismo = np.fromfile(f,dtype='<i4',count=-1,sep='')
        
        # For array, size information is encoded in 2 4-byte integers
        size = seismo[0:2]
        
        # After grabbing the size, take the remaining data and re-interpret as
        # doubles
        seismo = seismo[2:]
        seismo.dtype = '<d'
        return seismo.reshape(size[1]),size