Skip to content


Switch branches/tags

Name already in use

A tag already exists with the provided branch name. Many Git commands accept both tag and branch names, so creating this branch may cause unexpected behavior. Are you sure you want to create this branch?

Latest commit


Git stats


Failed to load latest commit information.
Latest commit message
Commit time
June 18, 2019 15:40
July 29, 2019 14:27
November 17, 2020 09:36
March 2, 2020 09:14
July 19, 2019 13:21
July 12, 2019 15:32
June 3, 2019 22:34
October 17, 2019 17:47


Docker Cloud Automated build PyPI DOI Build Status

Magnetic data is one of the most common geophysics datasets available on the surface of the Earth. Curie depth is the depth at which rocks lose their magnetism. The most prevalent magnetic mineral is magnetite, which has a Curie point of 580°C, thus the Curie depth is often interpreted as the 580°C isotherm.

Current methods to derive Curie depth first compute the (fast) Fourier transform over a square window of a magnetic anomaly that has been reduced to the pole. The depth and thickness of magnetic sources is estimated from the slope of the radial power spectrum. pycurious implements the Tanaka et al. (1999) and Bouligand et al. (2009) methods for computing the thickness of a buried magnetic source. pycurious ingests maps of the magnetic anomaly and distributes the computation of Curie depth across multiple CPUs. Common computational workflows and geospatial manipulation of magnetic data are covered in the Jupyter notebooks bundled with this package.


Launch the demonstration at




Mather, B. and Delhaye, R. (2019). PyCurious: A Python module for computing the Curie depth from the magnetic anomaly. Journal of Open Source Software, 4(39), 1544,

Navigation / Notebooks

There are two matching sets of Jupyter notebooks - one set for the Tanaka and one for Bouligand implementations. The Bouligand set of noteboks are a natural choice for Bayesian inference applications.

Note, these examples can be installed from the package itself by running:

import pycurious





You will need Python 2.7 or 3.5+. Also, the following packages are required:

Optional dependencies for mapping module and running the Notebooks:

Installing using pip

You can install pycurious using the pip package manager with either version of Python:

python2 -m pip install pycurious
python3 -m pip install pycurious

All the dependencies will be automatically installed by pip.

Installing with conda

You can install pycurious using the conda package manager. Its required dependencies can be easily installed with:

conda install numpy scipy cython

And the full set of dependencies with:

conda install numpy scipy cython matplotlib pyproj cartopy

Then pycurious can be installed with pip:

pip install pycurious

Conda environment

Alternatively, you can create a custom conda environment where pycurious can be installed along with its dependencies.

Clone the repository:

git clone
cd pycurious

Create the environment from the environment.yml file:

conda env create -f environment.yml

Activate the newly created environment:

conda activate pycurious

And install pycurious with pip:

pip install pycurious

Issue with gcc

If the pycurious installation fails due to an issue with gcc and Anaconda, you just need to install gxx_linux-64 with conda:

conda install gxx_linux-64

And then install pycurious normally.

Installing using Docker

A more straightforward installation for pycurious and all of its dependencies may be deployed with Docker. To install the docker image and start the Jupyter notebook examples:

docker run --name pycurious -p brmather/pycurious:latest


PyCurious consists of 2 classes:

  • CurieGrid: base class that computes radial power spectrum, centroids for processing, decomposition of subgrids.
  • CurieOptimise: optimisation module for fitting the synthetic power spectrum (inherits CurieGrid).

Also included is a mapping module for gridding scattered data points, and converting between coordinate reference systems (CRS).

Below is a simple workflow to calculate the radial power spectrum:

import pycurious

# initialise CurieOptimise object with 2D magnetic anomaly
grid = pycurious.CurieOptimise(mag_anomaly, xmin, xmax, ymin, ymax)

# extract a square window of the magnetic anomaly
subgrid = grid.subgrid(window_size, x, y)

# compute the radial power spectrum
k, Phi, sigma_Phi = grid.radial_spectrum(subgrid)

A series of tests are located in the tests subdirectory. In order to perform these tests, clone the repository and run pytest:

git checkout
cd pycurious
pytest -v

API Documentation

The API for all functions and classes in pycurious can be accessed from


  1. Bouligand, C., Glen, J. M. G., & Blakely, R. J. (2009). Mapping Curie temperature depth in the western United States with a fractal model for crustal magnetization. Journal of Geophysical Research, 114(B11104), 1–25.
  2. Tanaka, A., Okubo, Y., & Matsubayashi, O. (1999). Curie point depth based on spectrum analysis of the magnetic anomaly data in East and Southeast Asia. Tectonophysics, 306(3–4), 461–470.