seisplot

seisplot

Plotting of seismic data using variable-density or variable-area wiggle displays.

Description

The seisplot module provides a single, highly flexible method plot() to display seismic data. Methods toggle() and wipe() can be used to create animations that can be saved as images or movies. There is also a convenience function called spectrum() to display amplitude or phase spectra of seismic data.

The module was designed to get decent displays of seismic data in a way that is more convenient than going through many individual Matplotlib function calls explicitly. However, if required Matplotlib methods can also be called directly. The code is pure Python and kept deliberately simple to get students participating our Geophysics classes and exercises at university going with Python and seismic data.

Key features

• Variable-density image plots.
• Variable-area wiggle plots.
• Highly configurable settings like colors, line widths, colorbars, labels, axis ticks, grid lines, etc.
• Animated toggles between two or more seismic image plots.
• Animated wipes between two seismic image plots.
• Amplitude or phase spectrum plots.

      

Getting Started

Dependencies

Required: numpy, matplotlib

Installation

Install from PyPI:

$> pip install seisplot

Install directly from gitlab:

$> pip install git+https://gitlab.kit.edu/thomas.hertweck/seisplot.git

Editable install from source:

This version is intended for experts who would like to test the latest version or make modifications. Normal users should prefer to install a stable version.

$> git clone https://gitlab.kit.edu/thomas.hertweck/seisplot.git

Once you acquired the source, you can install an editable version of seisplot with:

$> cd seisplot
$> pip install -e .

Brief tutorial

For a demonstration of various features and much more, please visit the "examples" folder in the repository where several Jupyter notebooks (tutorials) are available.

Plotting seismic data (for instance, read with our seisio package and therefore available as Numpy structured array including trace headers) can be as simple as:

import seisplot

fig, ax = seisplot.plot(data, haxis="offset", width=4, height=6,
                        vlabel="Time (s)", hlabel="Offset (m)",
                        vmajorticks=0.2, vminorticks=0.1,
                        hminorticks=500, vgrid="major")

The variables fig and ax are standard Matplotlib figure and axis handles that can be used to tweak the display further. You could also create those first using fig, ax = plt.subplots(1, 1) and pass them to the plot() method. In this way, it is possible to, for instance, create several seismic displays in one figure, or create displays that share the y-axis (usually "time").

A display toggle can, for instance, be produced in the following way:

ani, fig, ax = seisplot.toggle([data_1, data_2, data_diff],
                               interval=1000, repeat_delay=0,
                               hlabel="offset (m)", vlabel="time (s)")

The returned animation-artist object can be used to save an animated image or a movie.

An animated wipe can, for instance, be produced in the following way:

ani, fig, ax = seisplot.wipe(data_1, data_2, blit=True,
                             nwipe=30, wipecolor="blue",
                             interval=5, repeat_delay=0,
                             hlabel="offset (m)", vlabel="time (s)")

Again, the returned animation-artist object can be used to save an animated image or a movie.

Plotting an amplitude spectrum of seismic data (for instance, read with seisio) can be as simple as:

win = np.hamming
fig, ax = seisplot.spectrum(data, amplitude=True, window=win, scale="dB",
                            hlabel="Frequency (Hz)", vlabel="Magnitude (dB)",
                            vgrid="major", hgrid="major")

A phase spectrum could be obtained in the following way:

fig, ax = seisplot.spectrum(data, phase=True, unwrap=True, degree=True,
                            hlabel="Frequency (Hz)", vlabel="Phase (deg)",
                            linewidth=2, title="Unwrapped phase spectrum")

These methods produce similar results to Matplotlib's magnitude_spectrum(), phase_spectrum() and angle_spectrum() functions but have more options. In case of entire gathers, the spectrum is averaged over all traces.

Main author

Dr. Thomas Hertweck, [e-mail]

Citation

If you use the "seisplot" module and you find it useful, getting some feedback would be very much appreciated. If you would like to cite this module, please use, for instance:

Hertweck, T. (2025). seisplot: A Python library for visualisation of seismic data. Version 1.3.0. url: https://gitlab.kit.edu/thomas.hertweck/seisplot/ (visited on 09/18/2025).

Adjust year, version and last visited date as required. Here's a BibTeX entry:

@software{seisplot,
  author  = {Hertweck, Thomas},
  year    = {2025},
  title   = {seisplot: A {P}ython library for visualisation of seismic data},
  url     = {https://gitlab.kit.edu/thomas.hertweck/seisplot/},
  urldate = {2025-09-18},
  version = {1.3.0}
}

License

This project is licensed under the LGPL v3.0 License - see the LICENSE.md file for details