8.3.2. Do an equivalent source inversion through GIFtools
Caution
The equivalent source code is currently under review. Beware of topographic effects in the forward data.
In this recipe, we step through an equivalent source inversion using magnetic data, starting from a TMI data file, a topography file, and a 3D mesh. At the bottom of this page, you will find the files used in this recipe as well as the complete project and output files. The steps for doing an equivalent source inversion for gravity and gravity gradiometry are very similar to this.
This recipe requires the following steps:
STEP 1: Import the necessary files
Start GIFtools or open an already existing project. Because we will be loading in a few files, it is easiest to set the working directory to the folder where the files are located and/or where you would like to do the forward model.
The first step is then to import the magnetic data, the topography file, and the 3D mesh. Importation of data can slightly differ, depending on the file type (GIF format, XYZ, etc). For this example, the magnetic data and the topography file are both in GIF format. Both are imported using the following menus:
For magnetic survey locations: Import → Data → Magnetics → GIF format
For topography: Import → Data → Topography → 3D GIF (3D XYZ)
For the 3D mesh: Import → Mesh → 3D
Note
Note that the mesh does not have any padding in the vertical direction, as required for the equivalent source problem.
Once the data are loaded, GIFtools will have the three items shown in the project tree on the left. It’s a good idea to regularly save your project!
Note
If you have a susceptibility model, you can import the model and skip directly to forward modeling the amplitude data.
Before we can solve for an equivalent source model for this magnetic data, we need to set the I/O for the data and uncertainties. This means we need to assign uncertainties to the data first assignUncert (link needed).
Use the menu structure to set the uncertainty for the magnetic data:
Data Manipulation → Assign uncertainties
A dialog pops up to add a percentage and/or a floor to each header. We completed the dialog by clicking the Assign checkbox for the header “B” and using 0 for percentage and 0.5 nT for the floor. Click “OK” to apply the changes. Now use the menu structure to set the I/O headers for the data:
Data Manipulation → Set I/O headers
For this example, make sure the total-field anomaly is set the “B” and the uncertainty to “B_uncert” by choosing these from the drop-down boxes. Click “OK” to close the dialog and apply the changes. The info panel for the magnetic data item will have updated and all minimum data requirements for the inversion are now set (indicated using ** in the info panel).
STEP 2: Invert for an equivalent source model
Now we can invert for an equivalent source model. To create the equivalent source item, use the menu structure in GIFtools as following:
Create → Processing → Magnetic equivalent source
Choose a directory to write files to and run the equivalent source in. Click “Open”. You will now see the equivalent source item in the GIFtools project tree on the left. For this example, it’s called “MagEquivalentSource (MAGequivalentSource)”. You can always rename the first part of the name.
The info panel shows you which data item, topography, and mesh are being used. In addition, it will show you the sensitivity filename.
The equivalent source item has its own menu which contains the following:
Working directory (to set or view)
Edit options
Copy options
Write
Run
Load results
It also has a visualization menu. These menus are identical to those for an inversion item so we refer you to the inversion recipe for more information on each one.
Use the menu and select Edit options. The dialog shown below pops up:
See the MAGES manual (which can be accessed via the GIF documentation page too) for detailed information about each parameter.
In this example, set the mesh, data, and topography. We will leave the remaining parameters set to the default values. Click “Apply and write files” to set the changes and write the files to the working directory. A message will appear to state that the files have been written. Click “OK”.
Next, use the menu to run the equivalent source. Use the “all” options to run both the sensitivity and the inversion. A command window will pop up and the code will run. Once finished, close the command window. In GIFtools, use the menu to load the last result:
Magnetic equivalent source → Load results → Last result. The equivalent source item will become a folder, with two subfolders (one for predicted data and one for the models).
If you view the model (use the visualization menu for the recovered GIFmodel), you’ll see it is a single layer, just below the topography.
Example
The magnetics data file, topography file, and 3D mesh along with a GIFtools project to repeat the above steps to calculate the equivalent source model: download