Infrasound array processing for volcano monitoring. Developed for operational use at the Alaska Volcano Observatory and slightly generalized for broader use.
Tested on python 3.7
- obspy
- pandas
- jinja2
Edit config.py. There are
- system level parameters: controlling data source and output location
- processing parameters: controlling data processing details (filters, window length, etc.)
- plotting parameters: slight control on how a few things appear (mostly this allows for differentiating between acoustic and hydroacoustic velocities)
Note that each individual array can have its own unique parameters to selectively override the default (though this wouldn't work well for the output directories)
Run on a cron every 10 minutes:
*/10 * * * * /path_to_file/array_processing.py >> /dev/null 2>&1
The webpage is automatically generated each time using jinja2 to populate the index.template file with the network and array structure from the config file.
(This method was adapted from a forked version of this code by Israel Brewster, and the html and javascript front-end was originally developed by Tom Parker)
Note that if you do not want to process a particular array (e.g. bad/no data), the array must remain defined in config.py in order to remain a viewing option on the webpage (for past data). In this case, you could just comment out the SCNL's for that particular array in config.py
Modify T1 and T2 in back_populate.py and run this script to generate images between time T1 and T2