Skip to main content
This is a pre-production deployment of Warehouse. Changes made here affect the production instance of PyPI (
Help us improve Python packaging - Donate today!

Reads and (partially) writes seismic data in Landmark's Geoprobe format

Project Description

python-geoprobe is a python module to read and write geoprobe horizons, volumes, and faults.


This implementation is based on reverse-engineering the file formats, and as such, is certainly not complete. However, things seem to work.


python-geoprobe supports basic reading and writing of geoprobe-formatted volumes, horizons, fault sticks, 2D seismic, colormaps, and TSurfs.

As an example of reading and display a slice from a volume:

import matplotlib.pyplot as plt
import geoprobe

vol = geoprobe.volume('examples/data/Volumes/example.vol')

# By default, this will not load the data into RAM. Instead,
# ```` will be a memmapped numpy array

# Now let's load everythign into RAM. Note that ``vol.load()`` returns the
# array, but ```` will be the same in-memory array after this, as well

# Indexing the volume works in "model" (inline/crossline/z) coordinates.
# We could also use ``vol.XSlice(2300)``, but indexing gives more flexibility.
data = vol[2300, :, :].T

# Display stretched. Note that for no vertical exaggeration, we'd use
# ``aspect=1/vol.dyW``
fig, ax = plt.subplots()
ax.imshow(data, cmap='gray_r', aspect='auto',
          extent=[vol.ymin, vol.ymax, vol.zmax, vol.zmin])
ax.set(title='Inline 2300', xlabel='Crossline', ylabel='Depth (m)')

We can also read/write geoprobe-formatted binary horizons (ascii horizons currently not supported):

import matplotlib.pyplot as plt
import geoprobe

hor = geoprobe.horizon('examples/data/Horizons/channels.hzn')

# Some basic information about the horizon, to show useful attributes
print ('The horizon has a total of %i points, %i of which are'
       ' auto-tracked' % (, hor.surface.size))
print 'The horizon has %i manually picked lines' % len(hor.lines)
print 'The inline coordinates range from', hor.xmin, 'to', hor.xmax
print 'The crossline coordinates range from', hor.ymin, 'to', hor.ymax
print 'The depth/time coordinates range from', hor.zmin, 'to', hor.zmax

# Display the horizon:
fig, ax = plt.subplots()

# hor.grid is a 2D array of the Z-values stored in the horizon
im = ax.imshow(hor.grid, cmap='gist_earth_r', origin='lower',
               extent=(hor.xmin, hor.xmax, hor.ymin, hor.ymax))

# Plot the manual picks
# Here, "line" is a numpy structured array with fields 'x', 'y', 'z', etc.
# "line_info" is a 4-tuple of (xdir, ydir, zdir, ID) (and is unused here)
for line_info, line in hor.lines:
   ax.plot(line['x'], line['y'], color='gray')

# Other information...
cb = fig.colorbar(im, orientation='horizontal')
cb.set_label('Depth in meters below sea level')
ax.set(title='An example horizon file', xlabel='Inline', ylabel='Crossline')


Joe Kington <>


Requires python >= 2.5 and numpy. Some functions (e.g. horizon.toGeotiff) require gdal and its python bindings. The plotting functions in utilities (e.g. utilities.wiggles) and some swfault functionality requires matplotlib.


Installation should be as simple as “python install”

Release History

Release History

This version
History Node


History Node


Download Files

Download Files

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

File Name & Checksum SHA256 Checksum Help Version File Type Upload Date
geoprobe-0.3.1-py2-none-any.whl (40.7 kB) Copy SHA256 Checksum SHA256 py2 Wheel Sep 26, 2015
geoprobe-0.3.1.tar.gz (36.1 kB) Copy SHA256 Checksum SHA256 Source Sep 26, 2015

Supported By

WebFaction WebFaction Technical Writing Elastic Elastic Search Pingdom Pingdom Monitoring Dyn Dyn DNS Sentry Sentry Error Logging CloudAMQP CloudAMQP RabbitMQ Heroku Heroku PaaS Kabu Creative Kabu Creative UX & Design Fastly Fastly CDN DigiCert DigiCert EV Certificate Rackspace Rackspace Cloud Servers DreamHost DreamHost Log Hosting