Gridding Geophysical Aeromagnetic Data from XYZ data in order to obtain a Grid with Magnetic Field, Elevation and other columns using Python or Simpeg


I have Geophysical aeromagnetic data that was collected in 2014 via the Helicopter. The raw data consist of XYZ whereby the X and Y is the location and the Z is the depth. I use Ubuntu Linux and I wanted to use the SIMPEG in order to grid the data but this tool gives me options for analyzing the already gridded data.

I was wondering if you can assist me with any Python library or how can I use SIMPEG in order to grid my data to obtain the Magnetic field in nanotesla, Elevation(m) from the raw XYZ data

Hi mandalfred!

I know that SimPEG has implemented equivalent sources, which would allow you to interpolate scattered magnetic data.

But I can also recommend Harmonica. In there you will find the EquivalentSources class that will allow you to interpolate harmonic data (like gravity and magnetics fields) using equivalent sources.

Hope it’s useful for you!

Hi Santisoler,

This was very helpful and I did my gridding. However, I wanted to work with Simpeg more. In this regard, I have my data consisting of Lat/Long and the Magnetic value in csv format. Going through SIMPEG documentation, I could not use csv format when it comes to simulation and inversion. The data formats that work with SIMPEG are .obs, .msh, .sus and Simpeg PF input file .inp format. I would like to know if from the Mag data that I have, I can produce these formats or else how can I use my data (csv format) with SIMPEG to achieve inversion and simulation.

Thank you for reading through and your response is highly appreciated.