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.

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 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.

   

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 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").

Main author

Dr. Thomas Hertweck, [e-mail]

License

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