Entre a la pagina que dijo Gabian (http://web.utah.edu/thorne/software.html) y me baje el script rsac.m
Este script tambien funciona con Octave!
La forma de usarlo es muy simple (incluso lo dice en los comentarios del script mismo) y es la siguiente:
octave:49> ls /home/leo/stuff/facu/señales_y_sistemas/Proyecto
Enunciado.pdf my_processing.log scripts señales
Aca scripts y señales son 2 directorios. En scripts tengo los scripts que voy a usar para el TP y en señales todos los archivos que me baje de las estaciones sismograficas
octave:51> ls /home/leo/stuff/facu/señales_y_sistemas/Proyecto/señales
2011.107.01.57.44.2250.GT.CPUP.00.BHN.M.SAC
2011.107.01.57.46.1250.GT.CPUP.00.BHZ.M.SAC
2011.107.01.57.46.6000.GT.CPUP.00.BHE.M.SAC
2011.107.01.57.49.7483.IU.LCO..BHE.M.SAC
2011.107.01.57.55.3483.IU.LCO..BHZ.M.SAC
2011.107.01.58.01.9483.IU.LCO..BHN.M.SAC
2011.107.01.58.45.6195.IU.TRQA.00.BHZ.M.SAC
2011.107.01.58.48.1695.IU.TRQA.10.BHZ.M.SAC
2011.107.01.59.31.5000.GT.PLCA.00.BHN.M.SAC
2011.107.01.59.35.4500.GT.PLCA.00.BHZ.M.SAC
2011.107.01.59.39.2000.GT.PLCA.00.BHE.M.SAC
2011.107.02.04.08.1450.TR.SVB..BHE.M.SAC
2011.107.02.04.09.6350.TR.SVB..BHN.M.SAC
2011.107.02.04.10.6250.CU.GRGR.00.BHZ.M.SAC
2011.107.02.04.12.1650.TR.SVB..BHZ.M.SAC
2011.107.02.04.22.9500.G.FDF.00.BHN.M.SAC
2011.107.02.04.23.6500.G.FDF.00.BHE.M.SAC
2011.107.02.04.27.3500.G.FDF.00.BHZ.M.SAC
Student.sac.tar
octave:53> ls /home/leo/stuff/facu/señales_y_sistemas/Proyecto/scripts
rsac.m
octave:55>
Para poder usar el script rsac.m me tengo que parar en el directorio donde esta:
octave:55> cd /home/leo/stuff/facu/señales_y_sistemas/Proyecto/scripts
octave:57> pwd
ans = /home/leo/stuff/facu/señales_y_sistemas/Proyecto/scripts
octave:59>
Ahora todo lo que hay que hacer es cargar el archivo .SAC en una variable
octave:60> senal = rsac ('../señales/2011.107.01.57.44.2250.GT.CPUP.00.BHN.M.SA
C');
octave:62>
Aca, senal es una variable con 3 columnas con la siguiente informacion (como se ve en los comentarios del script rsac.m):
% Column 1 contains time values.
% Column 2 contains amplitude values.
% Column 3 contains all SAC header information.
Para comprobar que funciono podemos plotear las primeras 2 columnas de la señal y ver como quedo:
octave:62> plot(senal(:,1),senal(:,2))
octave:64>
Les deberia abrir un grafico con la señal que van a analizar.
Espero que sirva para todos los que tenian problemas cargando los archivos SAC.
Saludos