[SACdata,SeisData,filenames]=SACSUNMAT('file1','file2',...,'filen')
reads n SAC files file1, file2, filen with SUN byte order and converts them to MATLAB format. The filenames can contain globbing characters (e.g. * and ?). These are expanded and all matching files loaded.
SACSUNMAT( cellarray ) where cellarray={'file1','file2',...,'filen'} is equivalent to the standard form.
SACdata is an nx1 struct array containing the header variables in the same format as is obtained by using MAT function of SAC2000. SACdata(i).trcLen contains the number of samples.
SeisData is an mxn array (where m=max(npts1, npts2,...)) containing the actual data.
filenames is a nx1 string cell array with the filenames actually read.
Note that writing
[SACdata,SeisData]=sacsun2mat('file1','file2',...,'filen')
is equivalent to the following sequence
sac2000
READ file1 file2 .. filen
MAT
(in fact the failure of above sequence to work properly on my system motivated this script).
SACSUN2MAT was written by F Tilmann (tilmann@esc.cam.ac.uk) based on sac_sun2pc_mat by C. D. Saragiotis (I copied the routines doing the actual work from this code but used a different header structure and made the routine flexible). It was tested on MATLAB5 on a PC but should work on newer versions, too. |