segysak 0.5.1

SEG-Y Seismic Data Inspection and Manipulation Tools using Xarray

Access the full documentation for SEGY-SAK from Gh-pages

SEGY-SAK aims to be your Python Swiss Army Knife for Seismic Data.

To do this SEGY-SAK offers two things; a commandline interface (CLI) for inspecting and converting SEG-Y data to a more friendly format called NETCDF4, and by providing convenience functions for the data using xarray. This includes lazy loading and streaming of data via the Xarray backend engine provided by SEGY-SAK.

We try hard to load the data the same way every time so your functions will work no-matter which cube/line you load. The xarray conventions we use are outlined in the documentation.

Why NETCDF4? Well, NETCDF4 is a fancy type of enclosed file binary format that allows for faster indexing and data retrieval than SEGY. We try our best to scan in the header information and to make it easy (or easier) to load SEGY in different formats, different configuration (2D, 2D gathers, 3D, 3D gathers). We do all this with the help of segyio which is a lower level interface to SEGY. If you stick to our xarray format of files we also offer utility functions to return to SEG-Y so you can export to other software.

Current Capabilities

• CLI:

  • Convert 2D, 3D and gathers type SEG-Y to NETCDF4 and back. The NETCDF4 files are one line open with xarray.open_dataset.

  • Extract sub-volumes via cropping operations.

  • Read the EBCIDC/text header of the seismic file.

  • Perform a limited header scan.

• Xarray and Python API:

  • Load SEG-Y data to a xarray.Dataset, dimensions are created by analysing trace header byte locations provided by the user.

  • Access header information and text headers in Python with convenience functions.

  • Select traces by UTM X and Y coordinates.

Installation

SEGY-SAK can be installed by using pip. Dependencies are specified in the pyproject.toml.

Python Package Index via pip

From the command line run the pip package manager

python -m pip install segysak

Install from source

Clone the SEGY-SAK Github repository and in the top level directory run

python -m pip install

CLI Quick Start

The command line interface (CLI) provides an easy tool to convert or manipulate SEG-Y data. In your Python command-line environment it can be accessed by calling segysak.

For a full list of options run

segysak --help

Any SEG-Y files converted using the convert command. For example

segysak convert test.segy

Can be loaded into Python using xarray.

test = xarray.open_dataset('test.SEISNC')

Complete Documentation

The complete documentation for SEGY-SAK can be found at Gh-Pages