Open-source software for two-dimensional Fourier processing of gridded magnetic data

Abstract

Magnetic data are widely available and useful in many exploration and applied-geophysics projects. The magnetic data are usually processed, imaged, and interpreted in commercial software packages. The algorithms used in these packages are sometimes difficult to check or tune, and the code is not available for review. However, these packages often have an application programming interface (API) for people to access data and undertake their own processing and data enhancement. In many cases, these APIs use the Python programming language. In the course of developing a new method for transforming magnetic data called reduction to pole and vertical dip (RTPVD), the initial test code was written in Python. This initial code was then rewritten and incorporated into GAMS, an open-source software package capable of using a Python API to read from and then write transformed (or enhanced) data to a commercial database. In addition to RTPVD, the other enhancements GAMS can generate are the zeroth-order analytic-signal amplitude (ASA0), tilt, spatial derivatives of ASA0, the zeroth-order local wavenumber, the first-order analytic-signal amplitude, and the apparent susceptibility. These transformations require that the data be transformed to the wavenumber domain using a fast Fourier transform (FFT), operated on, and then transformed back to the space domain. The FFT and some of the preprocessing steps can be done with a number of built-in Python tools. For the preprocessing steps, some of the available Python options are fast, but they can occasionally introduce unwanted artifacts. Our open-source tool allows users to test the different options and check the intermediate steps to ensure the result is appropriate.

Remanent Magnetization in the Proterozoic Basement Rocks of Iraq: the Southern Desert in Focus

“…For implementing the RTPVD transformation, and subsequently the estimation of the AMS of the Proterozoic basement rocks, we use an open-source Python code called GAMS presented by Smith et al (2022b). GAMS can be applied to any anomalous magnetic field data (Smith et al, 2022b).…”

“…For implementing the RTPVD transformation, and subsequently the estimation of the AMS of the Proterozoic basement rocks, we use an open-source Python code called GAMS presented by Smith et al (2022b). GAMS can be applied to any anomalous magnetic field data (Smith et al, 2022b). Ten grids are calculated when applying GAMS to the anomalous field, and eight of them are pre-processing grids, including derivatives of the field in three directions (two spatial in the x, and y directions, and the third is vertical in the z-direction), ASA, tilt angle, and zero-order local wavenumber (for more detail see Smith et al, 2022a and b), and the other two grids are the RTPVD and AMS.…”

Remanent Magnetization in the Proterozoic Basement Rocks of Iraq: the Southern Desert in Focus

“…For implementing the RTPVD transformation, and subsequently the estimation of the AMS of the Proterozoic basement rocks, we use an open-source Python code called GAMS presented by Smith et al (2022b). GAMS can be applied to any anomalous magnetic eld data (Smith et al, 2022b).…”

“…For implementing the RTPVD transformation, and subsequently the estimation of the AMS of the Proterozoic basement rocks, we use an open-source Python code called GAMS presented by Smith et al (2022b). GAMS can be applied to any anomalous magnetic eld data (Smith et al, 2022b). Ten grids are calculated when applying GAMS to the anomalous eld, and eight of them are pre-processing grids, including derivatives of the eld in three directions (two spatial in the x, and y directions, and the third is vertical in the z-direction), ASA, tilt angle, and zero-order local wavenumber (for more detail see Smith et al, 2022a and b), and the other two grids are the RTPVD and AMS.…”