Getting started

Looking for an overview? A good place to start is this recorded example (~20 min) that uses our tools to download, process, and interpolate public domain gravity data from the Bushveld Complex in South Africa.

No time for a video? You can also browse the Jupyter notebook used in the example.

New to Python? We provide links to excellent free resources for starting your Scientific Python journey.

Used Fatiando for research? Please cite our tools in your publications.

Source code for this demonstration: leouieda/2021-06-22-gfz

A toolbox for Geophysics

Pooch: Easily download datasets

Pooch is the easiest way to download data files to your computer.

It is used to manage sample data downloads not only by our own tools but also other popular Scientific Python libraries: scikit-image, MetPy, xarray, SHTOOLS, satpy, icepack, histolab, yt, napari, and more.

Download Australian gravity data from figshare
import pooch

# The Digital Object Identifier of the data
doi = "10.6084/m9.figshare.13643837"
# Known MD5 checksum (from figshare)
checksum = "md5:16c94a792003714efee2bdb4f3181d3a"
# Download the netCDF file and check integrity
fname = pooch.retrieve(

# Load the data with xarray
import xarray as xr

# fname is the path to the file
data = xr.load_dataset(fname)

Running this code multiple times will only result in a single download because the data are cached where Pooch can find it.

Verde: Gridding, machine learning style

Verde offers spatial data processing and interpolation (gridding) with a sprinkling of machine learning.

Vertical ground velocity in California interpolated from GPS data with and without weights based on data uncertainty.

Harmonica: All things potential fields

Harmonica is our library for processing, forward modeling, and inversion of gravity and magnetic data. Our goal is to incentivise good practices by carefully designing the software and offering state-of-the-art methods with efficient implementations.

Residual gravity disturbances of the Bushveld Complex, South Africa, gridded to a uniform height with equivalent sources.

Boule: Ellipsoids and normal gravity

Boule defines reference ellipsoids for calculating normal gravity of the Earth and other planetary bodies (Moon, Mars, Venus, Mercury).

Normal gravity of the WGS84 ellipsoid calculated at the Earth's surface using an analytical expression (no free-air correction required).


We are always happy to welcome anyone who is interested in getting involved! Whether it be coding, teaching, designing, or just hanging out. Come say “Hi! 👋🏾” on Slack or at our next online gathering.

Getting involved in open-source can be great way to meet new people, improve your coding skills, and make an impact in your field.

Happy community members at a Fatiando Community Call.


  • 2021/09/01: Documentation for the Python 2.7 fatiando package has been moved to

  • 2021/09/01: Fatiando is now on LinkedIn! Give our page a follow to keep up with the latest releases, events, and other news.

  • 2021/05/20: Santiago, Agustina, and Leo gave a talk to the Geophysical Society of Houston about using Fatiando for potential field data (slides are available).


  • Weekly Fatiando Development Calls: we discuss various aspects of the project. All are welcome, regardless of skill level and prior knowledge! Notes and connection details: fatiando/meeting-notes

  • AGU 2021: Going to the AGU Fall Meeting? Come to our talk! Details at fatiando/agu2021.

Latest releases: Pooch latest version Verde latest version Harmonica latest version Boule latest version