Reading events from a Seisan database - geoscience-community-codes/GISMO GitHub Wiki

INCLUDE DEMO DATASET FROM MVOE_ !!!!!!!!!!!!!!!!!!!

Here we load events from a Seisan catalog. A Seisan "Sfile" contains all the metadata for 1 event. These Sfiles are stored in a flat-file database structure the path to which is: $SEISAN_TOP/REA/databaseName. Sfiles are organized in year/month subdirectories under this path.

The following will navigate this where in this case $SEISAN_TOP = '/raid/data/seisan' and the databaseName is MVOE_ which stands for the Montserrat Volcano Observatory Event database. (In Seisan, databaseName is limited to exactly 5 characters).

This example will load all Sfiles for January 1st, 2000. This is a slow function to run as MATLAB is slow at parsing text files, and there are many events per day in this particular database.

>> catalogObject = readEvents('seisan', ...
        'dbpath', fullfile('/raid','data','seisan','REA','MVOE_'), ...
        'starttime', '2000/01/01', ....
        'endtime', '2000/01/02');

This Events object can now be explored using eev, which is modelled on the program of the same name in Seisan. This is an interactive function.

>> catalogObject.eev()

You will see:

2000-01-01 00:51:06 l   NaN:  ?

From there type h to see the list of options. For example, f takes you to the next event, s will show you the S-file, p will plot the corresponding waveforms providing the waveform data are available at the path indicated in the S-file.