Introduction Mapping the depth of the oceans globally was one of the greatest successes of geophysics in the 20th century. Without bathymetry, we would not know the locations of mid-ocean ridges, volcanic seamounts, and transform faults, to name a few, all of which are key elements of plate tectonics. But […]

How to Load Sandwell & Smith Gravity Data in Python?

Introduction I was in Paris last week for the 2017 Subsurface Hackathon. It was superbly organised by Agile Scientific and Total. Together, they managed to draw a large (about 60 persons) and diverse crowd to participate to the challenge of producing some working and usable software in only two days. […]

Three notebooks to jump start a data science project

Introduction This post is associated with a Jupyter notebook available on GitHub. In a previous post, I explained how to add colour maps to the 3D environment of OpendTect. The method is simply to convert the RGB image of the map into an indexed colour image. The resulting grid can […]

How to add cross-sections to OpendTect?

Introduction This post is associated with a couple of Jupyter notebooks available on GitHub. OpendTect is a great piece of software that allows you to load, process and interpret seismic data. In OpendTect, 2D lines and 3D volumes are displayed in a nice 3D environment that is easy to manipulate. Horizons, […]

How to add maps to OpendTect?

Summary Hillshading simulates the variable illumination of a surface by a directional light source. It is a great method to represent relief on a map and it works very well with potential field data too, not only with topographic data. Creating shaded maps in Python with matplotlib is easy and a few examples are provided […]

Hillshading with matplotlib

Summary When displaying geophysical data in a map, one may find it difficult to show both the presence of extremes and the subtle variations in the background signal. Histogram equalization is there to help, as it redistributes intensities and increases the contrast. In this new implementation for Python and matplotlib, the equalization is applied […]

Histogram Equalization in Python and matplotlib

Summary This post might not have much to do with geophysics but it is another way to share information and ideas about my interests. I am comparing in  a basic and simple way the results of a 5 km run at several locations in and around Leeds. In the process, I have […]

Basic statistics with Parkrun data

Working with seismic data is great but sometimes technical issues such as loading and converting files come in the way of the more interesting part of the job. One of the most common problems I have encountered in my career is the mismatch between the coordinate system of seismic data and […]

Manipulate SEGY files with segy2segy