Skip to main content

A Python Package for Surface Wave Inversion Pre- and Post-Processing

Project description

swprepost - A Python Package for Surface Wave Inversion Pre- and Post-Processing

Joseph P. Vantassel, The University of Texas at Austin

DOI PyPI - License CircleCI Documentation Status Language grade: Python Codacy badge Codecov PyPI - Python Version

Table of Contents


About swprepost


swprepost is an open-source Python package for performing surface wave inversion pre- and post-processing. swprepost was initially developed by Joseph P. Vantassel under the supervision of Professor Brady R. Cox at The University of Texas at Austin. The package continues to be developed by Joseph P. Vantassel.

The package includes multiple class definitions for interacting with the various components required for surface wave inversion. It is designed to integrate seamlessly with the dinver module of the popular open-source software geopsy, however has been written in a general manner to ensure its usefulness with other inversion programs. Furthermore, some of the class definitions provided, such as GroundModel, may even be of use to those working in the geotechnical or geophysical fields, but who do not perform surface wave inversion.

If you use swprepost in your research or consulting we ask you please cite the following:

Joseph Vantassel. (2020). jpvantassel/swprepost: latest (Concept). Zenodo. http://doi.org/10.5281/zenodo.3839998

Vantassel, J.P. and Cox, B.R. (2021). SWinvert: a workflow for performing rigorous 1-D surface wave inversions. Geophysical Journal International 224, 1141-1156. https://doi.org/10.1093/gji/ggaa426

Note: For software, version specific citations should be preferred to general concept citations, such as that listed above. To generate a version specific citation for swprepost, please use the citation tool for that specific version on the swprepost archive.

A Few Examples

All examples presented here can be replicated using the Jupyter notebook titled ReadmeExamples.ipynb in the examples/basic directory.

Import 100 ground models in less than 0.5 seconds

time_start = time.perf_counter()
gm_suite = swprepost.GroundModelSuite.from_geopsy(fname="inputs/from_geopsy_100gm.txt")
time_stop = time.perf_counter()
print(f"Elapsed Time: {np.round(time_stop - time_start)} seconds.")
print(gm_suite)
Elapsed Time: 0.0 seconds.
GroundModelSuite with 100 GroundModels.

Plot the ground models

fig, ax = plt.subplots(figsize=(2,4), dpi=150)
# Plot 100 best
label = "100 Best"
for gm in gm_suite:
    ax.plot(gm.vs2, gm.depth, color="#ababab", label=label)
    label=None
# Plot the single best in different color
ax.plot(gm_suite[0].vs2, gm_suite[0].depth, color="#00ffff", label="1 Best")
ax.set_ylim(50,0)
ax.set_xlabel("Vs (m/s)")
ax.set_ylabel("Depth (m)")
ax.legend()
plt.show()

Plot of 100 best shear wave velocity profiles.

Compute and plot their uncertainty

fig, ax = plt.subplots(figsize=(2,4), dpi=150)
disc_depth, siglnvs = gm_suite.sigma_ln()
ax.plot(siglnvs, disc_depth, color="#00ff00")
ax.set_xlim(0, 0.2)
ax.set_ylim(50,0)
ax.set_xlabel("$\sigma_{ln,Vs}$")
ax.set_ylabel("Depth (m)")
plt.show()

Plot of the lognormal standard deviation of the 100 best shear wave velocity profiles.

Getting Started


Installing or Upgrading swprepost

  1. If you do not have Python 3.6 or later installed, you will need to do so. A detailed set of instructions can be found here.

  2. If you have not installed swprepost previously use pip install swprepost. If you are not familiar with pip, a useful tutorial can be found here. If you have an earlier version and would like to upgrade to the latest version of swprepost use pip install swprepost --upgrade.

  3. Confirm that swprepost has installed/updated successfully by examining the last few lines of text displayed in the console.

Using swprepost

To start learning about swprepost, we recommend walking through the provided examples.

  1. To access the examples you can download the latest release of the project archived on zenodo.

  2. Unzip the project folder titled swprepost-vX.X.X.zip. And move the example directory to any location of your choosing. You can now discard the other files and directories in the .zip.

  3. Explore the Jupyter notebooks in the basic directory for a no-coding-required introduction to the swprepost package. If you have not installed Jupyter, detailed instructions can be found here.

  4. Move to the adv directory and follow the Jupyter notebook title example_swinvert_workflow.ipynb for an example of swprepost applied in the context of the SWinvert workflow (Vantassel and Cox, 2021). This workflow demonstrates how to use swprepost to perform surface wave processing and swbatch for running batch-style surface wave inversion. For more information on swbatch see its GitHub page.

  5. Enjoy!

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

swprepost-2.0.0.tar.gz (73.5 kB view details)

Uploaded Source

Built Distribution

swprepost-2.0.0-py3-none-any.whl (109.9 kB view details)

Uploaded Python 3

File details

Details for the file swprepost-2.0.0.tar.gz.

File metadata

  • Download URL: swprepost-2.0.0.tar.gz
  • Upload date:
  • Size: 73.5 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.0 CPython/3.8.10

File hashes

Hashes for swprepost-2.0.0.tar.gz
Algorithm Hash digest
SHA256 8f854333b081bfcbb7c1926c1ab598d7f3db716f053bcc4599897f472ceb10ab
MD5 685a09a27e298528de6672eef5daab3d
BLAKE2b-256 fb636bfc2ba97049089c7617304f1549d01feed04665ea888dce18d10bd49cf7

See more details on using hashes here.

File details

Details for the file swprepost-2.0.0-py3-none-any.whl.

File metadata

  • Download URL: swprepost-2.0.0-py3-none-any.whl
  • Upload date:
  • Size: 109.9 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.0 CPython/3.8.10

File hashes

Hashes for swprepost-2.0.0-py3-none-any.whl
Algorithm Hash digest
SHA256 9a94adab51262b99580927729b4974d1bd76af2ba4548e14a722d9e451719a92
MD5 890a091072fccaa4a198cfe28999b967
BLAKE2b-256 58eaff093a1ceed96ac5091998f0c2af83f930b1c9f59c7583e4dce3bf7eb786

See more details on using hashes here.

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page