array_processing

Various array processing tools for infrasound and seismic data. By default uses least-squares to determine the trace velocity and back-azimuth of a plane wave crossing an array in sliding time windows. More advanced processing (such as least-trimmed squares) is easily integrated. Also provides tools to characterize the array response, uncertainty, source-location of a spherical wave crossing the array, etc. See documentation and example.py for more info.

General References and Suggested Citations

Least squares and array uncertainty:

Szuberla, C. A. L., & Olson, J. V. (2004). Uncertainties associated with parameter estimation in atmospheric infrasound arrays, J. Acoust. Soc. Am., 115(1), 253–258. https://doi.org/doi:10.1121/1.1635407

Least-trimmed squares:

Bishop, J. W., Fee, D., & Szuberla, C. A. L. (2020). Improved infrasound array processing with robust estimators, Geophys. J. Int., vol(num), ###–###. https://doi.org/10.1093/gji/ggaa110

Installation

We recommend that you install this package into a new or pre-existing conda environment. (Ensure that the environment contains all of the packages listed in the Dependencies section.)

To create a new conda environment for use with this and other uafgeotools packages, execute the following terminal command:

$ conda create --name uafinfra --channel conda-forge obspy

This creates a new environment called uafinfra with ObsPy and its dependencies installed.

Install the required packages from Github. Execute the following terminal commands to install waveform_collection:

$ conda activate uafinfra
$ git clone https://github.com/uafgeotools/waveform_collection.git
$ cd waveform_collection
$ pip install --editable .

The final command installs the package in "editable" mode, which enables the installed package to be updated with a git pull in the local repository. This install command only needs to be run once.

Next, install lts_array in a similar manner to waveform_collection.

$ cd ..
$ git clone https://github.com/uafgeotools/lts_array.git
$ cd lts_array
$ pip install --editable .

Finally, install array_processing.

$ cd ..
$ git clone https://github.com/uafgeotools/array_processing.git
$ cd array_processing
$ pip install --editable .

Dependencies

uafgeotools packages:

• waveform_collection
• lts_array

Python packages:

• ObsPy

Usage

Import the package like any other Python package, ensuring the correct environment is active. For example,

$ conda activate uafinfra
$ python
>>> import array_processing

Documentation is available online here. For a usage example, see example.py.

Authors

(Alphabetical order by last name.)

Jordan Bishop
David Fee
Curt Szuberla
Liam Toney
Andrew Winkelman