/pyatoa

A misfit assessment toolbox for seismic tomography

Primary LanguagePython

Python's Adjoint Tomography Operations Assitant

Misfit assessment for the modern tomographer


Documentation can be found on Read the Docs: https://pyatoa.rtfd.io (in development)


Pyatoa* is a waveform-based misfit quantification package, with additional tools for measurement aggregation, and visualizations of inversion results. It is meant to facilitate the assessment of seismic inversions. Under the hood, Pyatoa is built on, and provides a high-level API for, a few key Python packages:

ObsPy: for seismic data fetching, handling, processing and organization
Pyflex: a Python port of Flexwin, for automatic time window selection
Pyadjoint: evaluating misfit functions and creating adjoint sources
PyASDF: heirarchical, self-describing storage of seismic data
Pandas: large-scale aggregation and manipulation of measurement information

Pyatoa can be used in scripting, interactive Python environments, or written into larger workflow tools. Although applicable in a standalone maner, Pyatoa was designed as a tool to be used in conjunction with SeisFlows, an automated workflow for seismic inversions, and SPECFEM3D Cartesian, a numerical solver for seismic wave propogation.

The design philosophy of Pyatoa is easy-to-use custom-built objects that make tomography research flexible, rapid, and repeatable.

*pronounced Py-uh-toe-uh (ˈpaɪəˈtoʊə), inspired by the famed volcano Krakatoa