GEO2000 B
GEO2000 B
GEOMETRICS INC.
Magmap 2000
User Guide
24891-01 Rev. H
Magmap 2000
Magmap 2000 4.0, User Guide
©Geometrics Inc.
2190 Fortune Dr
San Jose, CA 95131
Phone 408.954.0522 • Fax 408.954-0902
[email protected]
2
NEW MAGMAP 2000 FEATURES
25 Oct 2002
Manual Addendum
There are several exciting new features in the October 2002 release of MagMap 2000. We list them here
as a manual addendum and as a PDF for website distribution.
I. The Problem:
In the contour map generation routine, the Legend was always visible and if you made the map with many
contour interval colors, the legend became bigger than the map!
The Solution:
Now by right clicking on the map you can select “Plot Color Legend “ box or not. If you turn it off, then
you can get a larger picture of the map on the screen for zooming and further analysis. Remember that
when in PLOT DATA display (profile display) you can select anomalies for flagging by right clicking on
the profile and selecting SET FLAG. When you do that, the flag appears on the profile (number 1 below),
on the position map and also on the color contour map showing the position of the anomaly. Put the cursor
on the anomaly and read the X, Y or Latitude Longitude of the cursor on the info bar at the bottom
3
4
II. The Problem:
Sometimes there may be spikes in the data that are hard to filter using the Despike feature of MagMap.
We have now added a much more powerful despike editor called the Range Despike editor.
The Solution:
With the profile plot on the screen, go to Filter then Range Despike. You will see an icon appear that
will allow you to draw a box on the plot with your mouse. All data that exceeds the minimum or
maximum boundaries of the box will be cut from the data set, but only with in the lateral boundaries of
the box. This gives fine control as well and sweeping editing abilities to the user.
5
After Despike:
After the Range Despike editor has been applied the data is now
modified and the results will be shown in the Contour Map and when
the data is exported to Surfer, MagPick or Geosoft for additional
processing. This is a quick and easy way to clean up the occasional
spike or dropout.
The Solution:
Bring in the basestation file and then go to Filter, Smooth Basestation Data. You can select a value from
1 to 10 for the smoothing, 10 being the highest filtering. Note that the filter curve is shown prior to your
acceptance. If you go back to Filter and Smooth Basestation Data, you can choose Accept which will
modify the data in memory. Then you can export it or plot it. Remember that your original data is never
changed (*.bin files) but you can save a *.stn file over the original *.stn file.
6
Basestation data 4 hours
7
Basestation filter 4 preview, further zoom
NOTE: This feature will be implemented on regular data by November 2002. Then you can apply
smoothing to both basestation and field data!
8
IV. The Problem
Zooming in to see the data had a maximum limit using the circle with a + in it (or an X+ or Y+)
The Solution:
Our programmers have come up with a magnifying glass zoom tool like those found in graphics
programs. The Icon is a circle with a small box in it. Simply select this icon and then draw a box around
the portion of data that you want to zoom.
When you have finished zooming and wish to return to the original coordinate display, simply click on
the green icon with the arrows pointing outward .
V. The Problem
Sometimes we may want to apply a “frequency filter” to our data that becomes a spatial frequency, i.e.,
we are looking for something specific, a source of a given size. Or we may want to eliminate a certain
high or low frequency component in the data so that we can more clearly see the “targets” of interest.
The Solution:
For instance, let us say that we are interested in the geological signal because we are looking for
magnetite concentrations associated with placer gold deposits. However, there are also magnetic
materials right on the surface due to weathering of the host rock. We can raise the sensor to “look”
deeper (essentially applying a low pass distance filter), but we can also filter out the high frequencies in
the following manner.
9
Bring in data with low and high frequency components and contour
Now look at the data plot and apply a factor 7 spline to it, press review
If it looks sufficiently “low passed” then accept the data and generate a contour map. Compare with the
raw data contour map above.
10
Or one can “look” even deep with higher spline factors.
What if we want to look at only the surface features, for instance if we are doing an archeological or UXO
survey, we may want to accentuate the higher frequencies and remove the deeper geologic signal. In this
case we apply a high order (high number) spline filter to the data and THEN CLICK THE DIFFERENCE
button so that the data that is presented when we press review is the difference between the raw data and
the splined data. Then in this fashion, only the higher data is saved into memory when we press the
Accept button.
Then when we contour the data we get high frequency (small target) data display.
11
VI. The Problem
Noise spikes due to operation of the sensor close to the dead zone (see dead zone sensor orientation
discussion in technical report MTR-120) can be troublesome to remove. However with our new range
despike and spline filter modules, the data clean up is simple and quick.
The Solution:
12
First we use the Range Despike tool to eliminate the obvious spikes:
Next we can do a quick spline fit for REVIEW using the Smooth Readings dialog box under Filter. If we
wish we can do some additional despiking after applying the filter to see the results:
13
And now apply a factor 3 spline filter (Review) again from the Filter dialog and see results.
14
Now compare with and without despiking and filter
For more information, please refer to the manual and visit our website often for more updates.
15
Table of Contents
1. WHAT IS MAGMAP? 23
3. IMPORTING A SURVEY 32
16
3.3.3. OPENING A G-856 SURVEY ASCII FILE 39
3.3.4. OPENING A MAGLOG PROJECT FILE 39
3.3.5. OPENING AN INDIVIDUAL MAGLOG FILE 40
3.3.6. OPENING A G-856 MAGMAP FOR DOS FILE 40
3.3.7. OPENING A MAGSEA PROJECT 40
3.3.8. OPENING A MAGSEA INDIVIDUAL FILE 40
3.3.9. OPENING A MGS FILE (MARINE GRADIOMETER SOFTWARE) 40
3.3.10. OPENING A G-858 BIN FILE AND CORRESPONDING .CSV FILE 40
3.3.11. OPENING SURFER™ DAT FILE. 40
3.3.12. SHOWING ALL FILES (*.*) 46
3.4. DEFINING THE GRID (G-858 SIMPLE SURVEY MODE OR G−856 ONLY) 46
3.5. CONCLUSION 47
17
5.1.1. CREATING A 2D OR 3D PLOT 88
5.1.2. ‘EQUALIZING COLORS’ TECHNIQUE FOR CONTOUR PLOT 90
5.1.3. PLOT OPTIONS FOR 2D/3D PLOT 91
5.1.4. USER INTERACTIONS FOR 2D/3D PLOT 93
5.2. SHADED RELIEF PLOT 94
5.3. STACK PROFILE PLOT 95
5.3.1. WHAT IS A STACK PROFILE? 95
5.3.2. HOW TO PLOT A STACK PROFILE 98
5.3.3. SCALING A STACKED PROFILE 100
5.3.4. HOW TO TURN “AUTO-SCALING” OFF 100
5.4. BASIC SCALING OPERATIONS 100
5.5. TRUE MAP COORDINATES 102
5.6. USING GRID LINES 102
5.7. MARKING ANOMALIES FOR FURTHER STUDY 107
5.7.1. SETTING A FLAG 108
5.7.2. REMOVING A FLAG 110
5.7.3. LOCATING A POSITION 110
8. GPS 133
18
8.4. GPS OFFSET 138
8.5. SMOOTHING POSITIONS 140
8.5.1. RECOMMENDED SMOOTHING PROCEDURE 145
8.5.2. WARNINGS ABOUT SMOOTHING 145
8.5.3. REFERENCES ON SPLINES 146
8.6. TROUBLESHOOTING 147
8.6.1. MAP APPEARS INCORRECT WHEN USING GEOGRAPHICAL COORDINATES.147
8.6.2. MAP APPEARS INCORRECT WHEN USING UTM COORDINATES. 147
8.6.3. CANNOT GET ORIGIN OF GRID TO CORRESPOND TO (0,0). 147
19
14. PRINTING 177
20
18.2. MANIPULATING THE G-858 DATA 221
18.3. THE BASE STATION DATA 223
18.4. ANALYSIS OF G-858 DATA FROM WITHIN MAGMAP 223
18.4.1. SETTING A PROPER SCALE 223
18.4.2. USING THE STACK PROFILE 225
18.5. EXPORTING THE DATA 227
18.6. FURTHER DATA MANIPULATION 228
18.7. ANALYSIS IN SURFER 230
18.7.1. USING THE CONTOUR AND SHADED RELIEF PLOTS 230
18.7.2. GETTING A POSITION OF THE ANOMALY FROM THE CONTOUR PLOT 232
18.7.3. ADDITIONAL OBSERVATIONS 233
21
Note to new G-858 Users.
This software replaces the DOS version of software also included with the G−858. In
some cases, G−858 and G−856 magnetometers were shipped with an early MSDOS
based version of MagMap. We strongly recommend that users upgrade to Windows 95,
98, 2000, or NT in order to employ the full power and utility of Magmap. In keeping
with the common naming conventions, this manual refers to Magmap and the newly
released MagMap2000. New features will be added to the MagMap2000 version in our
continual pursuit to deliver the highest levels of support to our customers in their quest
to “see beneath the surface of the earth”.
22
1. What is Magmap?
M
agmap is a post acquisition processing software used to analyze and process data
from a number of different instruments. It allows the user to download data from
the G−858, and modify the positions that were entered while doing a survey. It also
supports other instruments such as the G−856, G−880, G−881, G-877 the
Ohmmapper, MetalMapper and the EM-61. The resulting data grid is graphically displayed on
the PC screen as a Quality Control measure and to assure coherence between the original
acquisition parameters and the final locations. Editing functions include moving lines and
positions, changing the direction of lines, rotating lines, and deleting parts of lines or entire lines.
It also provides a means of performing diurnal corrections, plotting the acquired magnetic field
values and filtering them for spikes, dropouts or leading errors. Vertical and horizontal
(longitudinal and transverse) gradiometer data is supported. In cases where there is GPS data,
MagMap2000 has facilities to smooth this data and use these coordinates in place of mapped
coordinates. Finally, the program will export data for surface plotting programs, such as
Golden Software’s Surfer and GeoSoft’s suite of analysis plotting programs.
Binary Data
Files 1 - 5
File / Import
IBM PC
Binary Data
.BIN extension
ASCII Data
.STN extension
File / Export
23
1.1.1. Downloading Data from the G-858
The time, position and magnetic field data are stored in the G-858 console
in one of 5 files numbered 1 through 5. The data are stored in a compressed binary
format for maximum use of the available RAM. While the user may select ASCII
download via RS-232 with a program of his choosing (i.e. Procom), this can take up
to 20 times longer than downloading the compressed binary file with the Magmap
Import program and is not recommended. Magmap Import provides an
“Autospeed” function that will make RS-232 connections at up to 115 KBaud on
compatible systems. EM-61 data is treated like Magnetometer data, OhmMapper
data is processes within Magmap but is described in a separate manual, PN ------
Finally, after positioning the data, Magmap will write an ASCII file suitable
for input to a 3D surface plotting program, such as Surfer for Windows™ or
Geosoft™ Analysis and Plotting Software. Such files consist of x and y coordinates,
field values, difference between values (in the case of gradiometer data ), and time
and date for each magnetometer reading. These files are typically about 3 MB per
hour of data collected.
After data corrections are applied, the corrected raw data may be stored in
ASCII format with an .STN extension. This will be necessary if you have not
completed all the editing you wish to do, and wish to continue at a later time.
1.2. Requirements
The G-858 can collect an enormous amount of data. Marine Mags create
even larger files. Even modest size surveys will require large computer resources in
terms of RAM and hard disk space. You will need a minimum of 16 MB of RAM
24
installed on your computer. Plenty of hard disk space will be necessary as well, if
you want to store lots of data. Typical Windows95, Windows 98, Windows NT, and
Windows 2000 installations will have enough resources.
This program will not run under MS-Windows 3.1. It will run under
Windows NT 3.51, but some of the dialog boxes and windows will not look as
pictured in this document.
1.3. Installation
1. First, decide to which drive and directory you wish to install the Magmap
program. It is prudent to keep it in the default directory.
(C:\program_files\magmap2k)
3. Run the SETUP program found on that disk. There are several ways to do
this. One is to double click on My Computer, double click on the floppy
drive, then double click on SETUP. Or, click on START, then select RUN.
Type A:\SETUP and press ENTER.
1.4. Usage
The three modes of operation of the G-858 correspond to different ways of
manipulating the data inside the Magmap program:
You can determine that you have a valid base station file that overlaps with
the field data by going to your export menu. If you do not have the message
“No valid base station file has been opened” at the bottom of your export
menu, this means you have a valid base station file.
25
1.4.2. Simple Survey
You must describe your survey to the Magmap program so that the
program can properly locate the data points. You will manually input the starting
line direction, the line and mark spacing and whether you took data bi-directionally.
Note: You are able to edit and reposition data in both simple
survey and mapped survey modes. Simple survey requires some
positioning information to be entered. In mapped survey, the
positioning information is entered in the field and you will most
likely not need to edit any positions, or enter any further
information.
26
2. How To Quickly Get Started
This is section is intended for the user who wishes to quickly start downloading and
editing a G−856 or G−858 survey. Much more detail about the various functions is
given in later sections. Please read those sections for more useful information about
the capabilities of this program when you have time.
• To correct a walking deviation error (for example you did not walk
straight and end up at another line, select “Add Line Deviation”.
27
2.2.3. Selecting and editing segments
To select a segment click the right mouse button on the direction symbol (arrow)
located between two positions. Segment operations are only available for
continuous surveys.
• To correct a walking deviation error (for example you did not walk
straight between 2 positions and end up at another line, select “Add
Segment Deviation”.
• To select more positions press the control key (CTRL) and repeat
the previous operation.
• To select a single position press the CTRL key and click the left
mouse button inside the corresponding position.
28
• If the positions were already selected they will be deselected, if they
were not selected they will be selected (pressing the CTRL key
toggles the selection.)
• To deselect all the positions left click the mouse button anywhere.
2.2.5.2. EDITING THE GROUP OF POSITIONS
Once you have finished selecting the positions, right click the mouse button on the
green area of the map.
• To reverse and shift the selected segments select “Reverse And Shift
Selected Segment (s) by”.
29
2.5. How to view a profile of the readings
• To plot a profile of only the readings between two positions click
the right mouse button on the direction symbol (arrow) located
between the two positions, then select “Plot Mag Field”.
• To plot a profile of only the readings for a specific line, right click
the mouse button on the red or green square of the corresponding
line, then select “Plot Mag Field”.
• To plot a profile of the readings for the whole survey, right click the
mouse button on the green area of the map, then select “Plot Mag
Field”. Individual anomalies may be flagged using the right mouse
click. Select fag: dentifier flag will appear on position and contour
maps.
Position the cursor on one of the lines of the map and press your right mouse
button. Select “Plot Stack Profile on Top of Map”. Then select “Plot Sensors Setup” and
choose which sensor you want to plot. You should then see a series of profiles
appear over your lines. Plus Minus Keys control the amplitude of the stacked
profiles.
30
You can zoom in and out by using the zoom buttons (a toolbar) located at the top
of your screen. If they are grayed out, make sure that you have disabled “Always Fit
in Window” from the Scale Menu.
Despike will only remove spike(s) above 1000 nT. (Use despike with caution in
discrete survey mode, as the readings may change rapidly near an anomaly, and thus
the anomaly may be interpreted as bad data and filtered out.).
The diurnal correction will be done automatically if the time of any opened base
station file matches the time of the survey.
31
3. Importing a survey
3.1. Downloading a survey from a G−856 or G−858
Magnetometer
3.1.1. Select import from the menu
Select “File” from the menu. Choose “Import G−856 Data” if you have a G−856 magnetometer or
“Import G−858 Data” if you have a G−858 data.
Figure 1
Note: Most of the following examples use the G−858 magnetometer but the
procedure for downloading data from the G−856 magnetometer is the same.
In the Data Transfer Menu select “TRANSFER UNDER CONTROL OF PC” and press enter.
The G-858 download setup dialog box should now appear on the PC screen as shown in Figure
2.
32
3.1.3.3 3.1.3.1 3.1.3.2
217BSelectin Setting the 216BSetting
g the data set Communication the baud
port rate
Figure 2
3.1.3.5 3.1.3.4
Setting “download 218BSetting 3.1.3.6
only decompress “file name to 220BRea
later” save as” dy to
Make sure you did not select the communication port where the mouse or external modem is
connected.
33
3.1.3.2. SETTING THE BAUD RATE
The G-858 supports different baud rates ranging from 1200 BD to 115 kBD and an auto mode
“Auto Detect” where the PC will select the best baud rate. For most cases setting the Baud rate to
115kb should work just fine. If you are having problems downloading your survey at this baud
rate with lots of retries, lower the value to “19200” or smaller.
Once you exit Magmap your settings will be saved. You do not need to set these values again,
unless you connect the G-858 to a different port or are having problems at high baud rates.
3.1.3.3. SELECTING THE DATA SET
In the G-858 you cannot store more than 5 surveys at any time. Each survey is represented by a
number from 1 to 5 that correspond to the data set number.
By default every data set is selected to be downloaded (check marks near the data set number). To
prevent a data set from being downloaded click the left mouse button on the check mark of the
data set you do not want to download. The check mark should now be removed and the “file
name to save as” and “BROWSE” button of the corresponding data set grayed out.
Left clicking the mouse button again on the data set number will enable this data set to be
downloaded.
3.1.3.4. SETTING “FILE NAME TO SAVE AS”
Next to the data set number you should type the drive, directory name and the name of the file
you wish to save as on the PC.
“C:\Program Files\Magmap\dataset1.bin”.
Where:
You will see later that you can save your survey as an ASCII (*.STN) file. In any case, you should
always keep this binary file. You cannot edit this file. By caution, after downloading your survey,
you should back up this binary file. Given the cost of your survey it is worth doing this extra step.
If you wish to browse your computer to select another drive or directory, left click the mouse
button on the browse button. A standard Windows 95 “Save Dialog Box” will appear on the
screen.
3.1.3.5. SETTING “DOWNLOAD ONLY DECOMPRESS LATER”
The data downloaded from the G-858 are in a compressed format. If you wish only to download
your survey and not process the survey data now, you can keep your survey in the compressed
34
format to save some space on your hard drive. By default this feature is disabled. (Current hard
drive sizes may make this unnecessary.)
3.1.3.6. READY TO DOWNLOAD
Once you have defined all the data sets you wish to download, click the left mouse button on the
“DOWNLOAD NOW” button or “Cancel” to abort the download operation and close this
dialog box.
You should see the following screen while the downloading is happening.
This screen allows you to see the status of the downloading (whether a connection has been
established, etc.). You can also see the COM port selected and the baud rate.
• For every survey you download a new DOS window will appear on the screen
reporting the status of the downloading process.
Once your survey has been downloaded, if the option “Decompress later” has not been set, your
survey file will be opened automatically.
Note: If your downloading does not work, you can have Magmap use the DOS
utility, binxfer.exe. To enable this, you need to edit the Magmap.ini file by
going to the Windows system directory and double clicking on “Magmap.ini”.
Find the string, “BINXFER = 0” and replace it with “BINXFER = 1”.
Magmap will now use the DOS utility to download.
If you do not need to open the file manually, you can skip the next paragraph that explains how
to open an already downloaded file and go directly to paragraph 3.4, “Defining the grid (G-858
simple survey mode or G−856 only)”.
• G-858 ASCII file (Simple Survey, Mapped Survey and Base Station).
35
• G-856 ASCII file (Survey and Base Station).
• Magsea projects
• G-858 bin files with a .csv file that gives corrected GPS positions
Note: Take care in selecting the proper category. The program acts
differently for different types of projects.
And there are several ways to open these files.
Figure 3
Once the open dialog box appears on the screen, select the desired file then click on the Open
button.
36
3.2.2. Opening a file using Drag & Drop Windows Explorer
You can use Explorer to import several survey files.
Select every file you wish to import by pressing the CTRL key and the left mouse button. Then
drag these files while pressing the left mouse button onto the MAGMAP window.
MagMap will sometimes recognize the type of file that is being loaded, but occasionally it will not.
Each option in the “Open” dialog box is specifically designed for a particular type of file which
will be discussed below.
37
3.3.1. G-858 Compressed files
This option will allow you to open 858 compressed files.
If you are opening a mapped G-858 binary compressed file you do not need to enter any survey
information as it has been stored with your data in the file. Simple survey files require some
information, such as, line spacing, covered in Section 3.4.
A dialog box indicating decompression and status information will be shown on the screen
(Figure 4).
Figure 4
Otherwise if you are opening a G-858 ASCII file that has been decompressed using the windows
version of MagMap called Magmap you will not be prompted for this dialog box as the survey
information has been stored in the ASCII file.
Figure 5
Once you have selected the type of survey you did, click the left mouse button on the “OK”
button and your survey file will be read as shown on the next dialog box (Figure 6)
38
Figure 6
The line number of this dialog box indicates a count of lines of text of the file being read in and
not the survey line number.
You will also need to set the year in which the survey was started. This is because the G−856
stores the date as the Julian day – e.g. Day 221 out of 365 is recorded as 221.
Figure 7
Once you have selected the type of survey you did, click the left mouse button on the “OK”
button and your survey file will be read.
Note: There are some files that Maglog NT produces that are not used in
Magmap (e.g. comment files). Magmap will warn you that it cannot read these
due to an unknown file type. In this case, it is safe to ignore this message.
39
3.3.5. Opening an individual Maglog file
You can open an individual Maglog file by selecting the file type: “Maglog Files (*.*)”. You can
then individually open each file associated with your project.
No user input is required when opening Maglog files except for Interpolator files using open as
sufer.dat, etc.
You may get a warning message such as the one below warning you of duplicate time stamps.
You can choose to smooth the readings immediately by pressing “Yes” or you can leave the data
as it is and correct it later.
If you have readings with duplicate time stamps, you should get an error message such as
described in the previous section.
Depending on the user selection, the following information can be loaded from a SURFER data
file:
40
1. Positions and Data. In this case it is assumed that each reading in the SURFER file has
an associated position.
2. Positions ONLY. In this mode the user loads only positions, i.e., the file is used in the
same way as a GPS file in case of MagLog NT Survey files. Magnetic (or other) data
should be loaded separately and may have a different sample interval than the position
sample rate. Data can loaded as a SURFER file as well (see below) or as MagLog NT
device file(s). While the user can load multiple data files, only one position file is allowed.
Note that the Line Number file can also be loaded in this case.
3. Data ONLY. In this mode the user loads only data. Multiple SURFER files can be
loaded in this mode. Positions and line numbers are loaded separately as MagLog NT
GPS and Line Number file. Alternatively positions can be loaded as a SURFER file (see
above).
4. Use as a Base Station. In addition, any file of the structure with magnetometer and
time information can be used as a base station file for diurnal correction as long as the
dates and times of the survey and base file intersect (base data must overlap rover data).
This means that virtually any base station magnetometer can be used with MagMap to
diurnally correct any other magnetometer data file.
Go to File/Open dialog and select SURFER *.dat files (Interpolator *.INT, Text *.txt) files as the
file type. Then select your file of interest and press “Open”. The program makes no assumptions
about the possible file format and therefore a “Data file definition” dialog is required:
41
Use this dialog to specify the type of data provided in the separate file columns. MagMap uses the
first file line (header) to name the data columns. If there is no header information, then the first
line of data is used to define each column.
• Position information. Specify the columns with position information. For example,
1:Mag_X means that 1st column has name “mag_X” and is selected as X position. Y
position is selected as “2:mag_Y”. For latitude and longitude data the check box below
must be checked. For UTM or local coordinate systems do not check this box.
NOTE: If the data is Lat Long you MUST check the box or your resultant map will
show no positions (green screen, one red and one green position box, nothing
more.)
• Time information. To display data correctly MagMap needs time and date information.
It is looking for key words “time” and “date” in the file header and will automatically
select them. However you can always select other columns as appropriate. For a data file
without time information, date and time values will be generated internally because
MagMap needs these values to be able display the data in a linear fashion. Check “File
42
has no time stamps, assign time by default” box to use default time stamps. Keep in
mind that in this case the time displayed by MagMap is synthetic, and so no time related
operations (like base station correction) can be used.
o US time format as Month / Day / Year. Fields are typically separated by forward
slash; year has short notation. (like 07 stands for 2007).
o European time format. The same as US, but day goes first.
Note that the user can use any characters to separate fields in the date string, except digits,
commas, tabs, colons and spaces. For instance 10—12—1998 is perfectly correct time
stamp. Year can be recorded as full year (as 2007) or as short year (07). Years less then
fifty are to be considered of 21st century, years more than 50 belongs to 20th century.
Examples: 07 is 2007, 77 is 1977.
• Line information. MagMap searches the file header for the word “line” to use that
column as the line number. However you can always select a different column for line
number as appropriate. If no line numbers are available, you can display the data as one
line (check “Don’t use lines”) or you can automatically assign line numbers based on the
distance between neighboring positions (check “Split data file into lines by distance” and
enter “Split distance”). When distance between two sequential reading positions exceeds
the “split distance” MagMap automatically starts a new line and increments the line
number.
• Data channels. The last step is to specify which columns contain the magnetometer
data. The “Data channels” list will list the columns selected for import. Press “Add data
channel” and the following dialog box appears with all data columns listed:
43
Select the magnetometer data channels (for example could include magnetometer field
data, signal strength, depth, altitude, etc.) and press “Ok”. Hold the “Crtl” key to select
multiple channels at once. The “Select channel to add” dialog box will be closed and the
channels added to the data channels list. To remove an already selected channel, highlight
it and press the “Remove” button. To remove all channels press the “Remove all
channels” button.
• How to load. This dialog box controls and constrains the type of data that can be loaded
into MagMap. Typically the user needs to load both data and positions; however it is
possible to load disparate data from separate files, as it described at the beginning of this
section. Select appropriate type here. Based on this selection, some input fields in the
dialog are grayed out: for instance, if only positions are loaded the data section is disabled,
if only data is to be loaded then the position section is disabled. Several “data only”
and/or base station files can be opened simultaneously to over different time spans. This
would normally be used when you have a multiple-day G-858 or MagLog survey with
base station files that begin and stop each day. Then you would load the entire field
survey and several base station files to apply diurnal correction on a day-to-day basis.
Keep in mind however that multiple base stations which were positioned at different
locations will have different base offsets and your final corrected data set under these
conditions may need to be block leveled. Observing the base station data files
concurrently on one plot will enable the user to insert an offset to the base station files to
enable leveling of the data upon export.
• Use Advanced Settings. This selection allows tuning of the data import routines. This
provides for removal of multi-string headers from a data files or use of a special character
denoting comment lines which the program will not load. In addition you can scale the
data while loading. (For instance in some systems a magnetic field value of 40000.000 nT
could be recorded as 40000000). Note that all previous selections (such as the
information denoting position and data columns, etc) are erased when this button is
checked, and the “Advanced” button becomes available. When clicked, it shows the
following dialog:
44
Here the following fields are available to the user:
• Header lines: Number of lines of text to skip from the beginning of the file.
• Comment characters: If strings begin with any of the characters in that field, it is assumed
to be comment line and is omitted.
• X, Y scale factors: Multiply positions by these numbers while reading the data.
• Data scale factor: the same as above, but for data fields. All data fields are multiplied by
this coefficient while reading from the files.
• Field separator: Describes how data columns are separated in the file. Possible values are
“SPACE” and “COMMA or SPACE”. In latter case either space or comma are
considered to separate fields.
Here is an example how INTERMAGNET IAGA-2002 1 format can be loaded using the MagMap
advanced features:
1
International Real-time Magnetic Observatory Network (INTERMAGNET) data is
available from http://www.intermagnet.org/Welcom_e.html; see “Conditions of
use” on INTERMAGNET web site. Note that data is recorded in UTC time.
45
Here the program skips the first 12 lines; in addition, each line starting with “#” is
considered to be a comment.
3.4. Defining the grid (G-858 simple survey mode or G−856 only)
When you open a G−858 simple survey mode or G-856 file you will need to define the
parameters of the grid. This does not apply to the G-858’s mapped survey mode because the
magnetometer has already stored the coordinates of the positions.
In the Grid Orientation dialog box (Figure 8) you will need to enter several parameters:
• First position origin: this is the origin of the survey’s first position.
• Survey direction: this is how you walked relative to the first position.
You cannot cancel the grid settings. The “Cancel” button is grayed out. You must set all the
parameters to plot the map.
46
You will see later that you can change these parameters, and even use the grid setting for mapped
survey mode.
Figure 8
3.5. Conclusion
Regardless of the way you import your data (by downloading your survey, or by using the open
dialog box Figure 3, or using Explorer) the following window (Figure 9) will appear on the screen.
Some surveys will only have the lower window (mapped) initially open.
47
This window will This window will
appear if you have appear if you have
imported a base imported a survey
station file file
Figure 9
Note: If no profile is initially open, you can create one by right clicking on the
mapped screen and selecting “Plot Mag Field”. If you click on a position
point, it will the profile of the line corresponding to the position point clicked.
If you click on another part of the map, you will get a profile of all the data.
48
4. Editing of the map
4.1. Map user interface description and definitions
A start of line is the first position of a line and is represented by a green square (Figure 10).
An end of line is the last position of a line and is represented by a red square (Figure 10).
A position (called Mark in G-858 and Station in G-856) is represented by a blue square (Figure
10). GPS positions also show as Blue Squares.
Duplicate positions are positions acquired at the same location (same X and Y coordinates) with
a different mark or line number. You get duplicate positions when you survey the same line or
segment twice or when you shift a line or segment on top of another line or segment. They are
represented by a double inlaid square.
A segment is represented by an arrow, that also indicates the direction the user was walking
between two positions (Figure 10).
A grid is represented by a white dotted rectangle (Figure 10). Gridlines can be configured using
the Scale menu.
Every magnetic reading is represented by a black point (Figure 10). If there are many readings
between two positions, it will look like a black line drawn between the two positions. By zooming
in, every reading will be represented by a single point.
When the mouse cursor is placed on top of a position, the status bar will indicate information
about the position such as time, mark, line, X and Y location.
• Zoom in: zoom in the map in the X and Y directions with the same zoom factor.
• Zoom out: zoom out the map in the X and Y directions with the same zoom
factor.
• Zoom in X: zoom in the map in the X direction only. The scale in the Y
direction remains unchanged.
• Zoom out X: zoom out the map in the X direction only. The scale in the Y
direction remains unchanged.
• Zoom in Y: zoom in the map in the Y direction only. The scale in the X
direction remains unchanged.
49
• Zoom out Y: zoom out the map in the Y direction only. The scale in the X
direction remains unchanged.
• Window Zoom: zoom in map in the user selected window. To select zoom
area, point mouse into one of the corners, click and drag mouse to the opposite
corner. Zoom rectangle is shown with “rubber band” After releasing of mouse
button selected rectangle is zoomed to extend s of the window. Areas outside of
selected rectangle are not shown. To bring picture back to full extend press “Fit
in window” button.
• Fit in window: the map is scaled to fit exactly inside the window. X and Y do
not have the same scaling factor. If the window is resized, the map will not be
scaled again to fit inside the window.
• Always fit in window: the map is scaled to always fit exactly inside the window.
X and Y do not have the same scaling factor. If the window is resized, the map
will be scaled again to fit inside the window.
50
ZOOM IN MAP TOOL BAR
X
ZOOM OUT X
SEGMENT
represented by
an arrow.
POSITION
represented by a
blue square.
Mark 1 X=0
Line 0 Y=100
07/26/95 13:28:03
START OF LINE
Figure 10
represented by a DATE AND TIME
green square. of position X LOCATION
Y LOCATION
Mark 0 X=0 of position.
of position.
Line 0 Y=0
STATUS BAR MARK NUMBER
indicating position of position.
MAG READINGS information
represented by a black LINE NUMBER
line or black dot. of position.
51
4.2. Selecting and editing a single position
To select a single position right click the mouse button on the position symbol (small blue square)
(Figure 11).
Figure 11
To perform the desired operation “Move Position To”, “Shift Position By” or “Delete
Position” move the mouse cursor on top of the corresponding item, then left click the mouse
button or use the keyboard arrow keys to select the item and press the “ENTER” Key.
Figure 12
Type in the X and Y field the coordinates of the new location, then click the left mouse button on
the “OK” button or pressed the “ENTER” key.
To cancel the operation left click on the CANCEL button or press the “ESC” key.
To switch between the X, Y, OK, and CANCEL controls, use the mouse or press the “TAB”
key.
52
You can copy and paste values between every edit control.
The following commands are common to all edit controls in this program and Windows 95/98
and will not be described for the other dialog boxes:
• To copy: press “CTRL” and “C” keys or “CTRL” and “INSERT” keys.
• To paste: press “CTRL” and “V” keys or “SHIFT” and “INSERT” keys.
• To undo press “CTRL” and “Z” keys or “ALT” and “BACKSPACE” keys.
• To redo press “CTRL” and “Z” keys or “ALT” and “BACKSPACE” keys.
• Click the right mouse button and a pop up menu will appear.
Figure 13
53
4.2.2. Shifting a position
Figure 14
Type in the X and Y field the offset to the new location, then click the left mouse button on the
“OK” button or pressed the “ENTER” key.
For example , in the survey plotted in Figure 10, if you shift by X=-2 and Y=-50 the position at
MARK 1 LINE 0 located at X=0 Y=100 you will obtain the picture in Figure 13.
Figure 15
If the information reported by this dialog box confirms that this is the position you want to delete
click on the “YES” button, otherwise click on the “NO” button to abort.
Deleting a position does not delete the magnetic field readings acquired during the survey.
If a position is deleted, magnetic field readings are linearly interpolated between the previous and
next position of the deleted position.
For example , in the survey plotted in Figure 10, if you delete the position at MARK 1 LINE 0
located at X=0 Y=100 you will obtain the picture in Figure 16.
54
Figure 16
Deleting a position does delete the position and the corresponding magnetic field acquired at
that location.
Figure 17
• To correct a walking deviation error (for example you did not walk straight and
end up at another line, select “Add Line Deviation”.
55
• To move only the start of line or end of line position, select “Move Position
To”.
• To plot every magnetic field reading of this line, select “Plot Mag Field”.
Figure 18
Type in the X and Y field the offset to the new location, then click the left mouse button on the
“OK” button or pressed the “ENTER” key.
Every position of the corresponding line will be shifted by the X and Y offset.
For example , in the survey plotted in Figure 10, if you shift by X=2.5 and Y=-50 LINE 0 you
will obtain the picture in Figure 19.
Figure 19
A line can be reversed (mirrored) in the X and Y directions, in the X direction only or the Y
direction only. The last two options should not be used for horizontal or vertical lines.
56
For most cases select “Reverse In Both Directions”.
Figure 20
For example, in the survey plotted in Figure 10, if you reverse LINE 0 you will obtain the picture
in Figure 21.
Figure 21
If your line is vertical and you did not walk straight and end up off line at the end of the line,
select “deviate in the X direction”. If your line was horizontal, select “deviate in the Y direction”.
57
Figure 22
For example, let us assume that in Figure 10 you walk off line on LINE 0 and really end up at the
end of line at X=2.5 and y=600. To fix this problem select “deviate in the X direction”. When
the following dialog box (Figure 23) appears on the screen type 2.5 and click on “OK”.
Figure 23
You should now have fixed your walking deviation as shown in Figure 24.
Figure 24
If by mistake you start line 0 at X=2.5 instead of X=0 and end up at X=0 and Y=600. First
reverse the line, then deviate the line in the X direction by 2.5 and reverse the line again.
58
4.3.4. Rotating a line
Figure 25
Type in the rotation field the angle of the rotation in degrees (from -360 to 360 degrees).
Type in the X and Y field the coordinates of the rotation origin. Once you enter all the fields, left
click the mouse button on the “OK” button or pressed the “ENTER” key.
For example , in the survey plotted in Figure 10, if you rotate LINE 0 by 180 degrees at the
origin X=0 and Y=0 you will obtain the picture in Figure 26.
Figure 26
59
4.3.6. Deleting a line
Figure 27
If the information reported by this dialog box confirms that this is the line you want to delete
click on the “YES” button, otherwise click on the “NO” button to abort.
Deleting a line will delete all the positions and magnetic field readings acquired on this
line.
For example , in the survey plotted in Figure 10, if you delete LINE 0, you will obtain the picture
in Figure 28.
Figure 28
60
In discrete survey mode there is only one reading per position and located at the position itself
therefore there is no segment.
In continuous survey mode readings are acquired between positions and linearly interpolated
between them.
To select a segment right click the mouse button between two positions on the direction symbol
(arrow). Segment operations are only available for continuous survey.
Figure 29
• To correct a walking deviation error (for example you did not walk straight
between 2 positions and end up at an other line, select “Add Segment
Deviation”.
• To plot the magnetic field between 2 positions, select “Plot Mag Field”.
• To delete a segment select “Delete Segment”. (This will delete the readings
between the two positions).
Every operation on a segment is basically the same as operating on a line. You can think of a
segment as being a line with only two positions.
61
4.4.1. Shifting a segment
Figure 30
Type in the X and Y field the distance to the new location, then click the left mouse button on the
“OK” button or pressed the “ENTER” key.
The two positions of the corresponding segment will be shifted by the X and Y offset and all the
magnetic readings inside that segment will be linearly interpolated between the two new positions.
For example , in the survey plotted in Figure 10, if you shift by X=-2.5 and Y=-50 the second
segment in line 0 you will obtain the picture in Figure 31.
Figure 31
If you shift that segment again Mark 1 and Mark2 will not be duplicated.
A segment can be reversed in X and Y direction, in X direction only and Y direction only. The
last two options should not be used for horizontal or vertical segments.
62
Figure 32
For example , in the survey plotted in Figure 10, if you reverse the second segment of line 0 you
will obtain the picture in Figure 33.
Figure 33
The positions at each end of the rotated segment have been duplicated in order to preserve the
interpolation of the readings in the segments on either side of the rotated segment.
If your segment is vertical and you did not walk straight and end up off line at the end of the
segment, select “deviate in the X direction”. If your segment was horizontal, select “deviate in the
Y direction”.
63
Figure 34
For example, let us assume that in Figure 10 you walk off segment on the second segment of line
0 and really end up at the end of the segment at X=2.5 and Y=200. To fix this problem select
“deviate in the X direction”. When the following dialog box (Figure 35) appears on the screen
type 2.5 and click on “OK”.
Figure 35
You should now have fixed your walking deviation as shown in Figure 36.
Figure 36
64
Mark2 has been duplicated.
Figure 37
Type in the rotation field the angle of the rotation in degrees (from -360 to 360 degrees).
Type in the X and Y fields the coordinates of the rotation origin. Once you enter all the fields, left
click the mouse button on the “OK” button or press the “ENTER” key.
For example , in the survey plotted in Figure 10, if you rotate the second segment of line 0 by 90
degrees at the origin X=0 and Y=100 you will obtain the picture in Figure 38.
Figure 38
65
Figure 39
For more information on this window, see the section about plotting readings (section 0 page
116). Readings are not shown in their correct linear position, only in the sequence that they were
taken.
This plot is used as a quality control device to insure that the data is free from excessive noise and
dropouts. It is very important to immediately remove any zero values or spikes from the data (see
section 7). After applying various filters, the data should be re-plotted for confirmation of results.
Figure 40
If the information reported by this dialog box confirms that this is the segment you want to delete
click on the “YES” button, otherwise click on the “NO” button to abort.
Deleting a segment will delete all the magnetic field readings acquired for this segment.
It may also delete the end positions if there are no other readings attached to these positions.
66
For example , in the survey plotted in Figure 10, if you delete the second segment in line 0, you
will obtain the picture in Figure 41.
Figure 41
• To select a group of positions click the left mouse button and draw a rectangle
around the positions to select (Figure 42). When you have finished drawing the
rectangle, release the left button. Every position inside that rectangle will be
selected (indicated by a white square) as shown in Figure 43.
Figure 42
Figure 43
67
• To select more positions press the control key (CTRL) and repeat previous
operation (Figure 44 and Figure 45).
Figure 44
Figure 45
• To select a single position press the CTRL key and click the left mouse button
inside the corresponding position.
• When you select positions, if the positions were already selected they will be
deselected (Figure 46 and Figure 47). If they were not selected they will be
selected (pressing the CTRL key toggles the selection.)
68
Figure 46
Figure 47
• To deselect all the positions click the left mouse button anywhere.
4.5.2. Editing a group of positions
Once you have finished selecting the positions click the right mouse button on the green area of
the map.
69
Figure 48
• To line up the selected positions along a slope select “Line Up Selected Positions
Along A Slope”.
The following segment operations are only available for continuous surveys.
• To reverse and shift the selected segments select “Reverse And Shift Selected
Segment(s) by”.
• To delete a segment select “Delete Segment”. (This will delete the readings of
every selected segment).
Once you have selected the operation to perform, one of the following dialog boxes will appear.
70
4.5.2.1. SHIFTING A GROUP OF POSITIONS
Figure 49
Type in the X and Y fields the distance for every position to the new location, then click the left
mouse button on the “OK” button or pressed the “ENTER” key.
For example , in the survey plotted in Figure 50, if you shift by X=10 and Y=0 all the selected
positions you will obtain the picture in Figure 51.
Figure 50
Figure 51
71
Figure 52
The behavior of the program differs depending upon whether you have performed a continuous
or discrete survey.
Deleting a group of positions does not delete any magnetic field readings acquired during the
survey.
If a group of positions is deleted, magnetic field readings are linearly interpolated between the
previous and next positions outside of the selected group of positions.
Deleting a group of positions does delete every position and the corresponding magnetic field
acquired at their location.
72
Figure 53
Click the right mouse button and select “Line Up Selected Positions) Along A Slope” from the
pop up menu.
The dialog box as shown in Figure 54 should now appear on the screen.
Figure 54
Type the coordinate of the top left corner of the map in the “point 1 coordinates” field.
Type the coordinate of the top right corner of the map in the “point 2 coordinates” field.
You should now get the right shape of the map as shown below in Figure 55.
73
Figure 55
74
Figure 56
• To reverse the direction of every line or every other line select “Reverse Grid”.
• To export your survey to Surfer or Geosoft Analysis and Plotting Software select
“Export”.
• To plot every magnetic field reading of the survey, select “Plot Mag Field”.
Once you have selected the operation to perform, one of the following dialog boxes or windows
will appear:
Figure 57
Type in the X and Y fields the offset to the new location, then click the left mouse button on the
“OK” button or pressed the “ENTER” key.
75
Every position of the survey will be shifted by the X and Y offset.
For example , in the survey plotted in Figure 10, if you shift the grid by X=5 and Y=100 you will
obtain the picture in Figure 58.
Figure 58
To reverse every line of the survey select “Reverse Every Line” and “Reverse In Both
Directions”.
To reverse every even line of the survey select “Reverse Every Even Line” and “Reverse In
Both Directions”. Every other line, starting with the first line acquired in the survey, will be
reversed.
To reverse every odd line of the survey select “Reverse Every Odd Line” and “Reverse In Both
Directions”. Every other line, starting with the second line acquired in the survey, will be reversed.
Figure 59
For example, in the survey plotted in Figure 10, if you reverse every line in both directions you
will obtain the picture in Figure 60.
76
Figure 60
Figure 61
Type in the rotation field the angle of the rotation in degrees (from -360 to 360 degree).
Type in the X and Y field the coordinates of the rotation origin. Once you enter all the fields, click
the left mouse button on the “OK” button or press the “ENTER” key.
For example , in the survey plotted in Figure 10, if you rotate the grid by 90 degrees at the origin
X=0 and Y=0 you will obtain the picture in Figure 62.
77
Figure 62
In the case of mapped survey mode, you can redefine your grid. This will change the coordinates
of every position acquired in your survey. You should only use this feature if you set up invalid
grid information in the G-858 magnetometer before doing your survey, such as invalid line or
mark spacings. Potentially, every position in the survey may be changed.
78
Figure 63
For more information on this window, see the section about plotting readings (Section 0 page
116). Readings are not shown in their correct linear position, only in the sequence that they were
taken.
• How do we transform local 858 coordinates into State Plane coordinates or some other
coordinate system? If both coordinate systems were defined completely without errors it
would be sufficient to measure only two corresponding points in each system.
• How do we take into account possible irregularities in the 858 coordinate system such as
non-parallel lines, curved lines or irregular mark spacing along the line? In this case it is
not enough to define the transform using only two points. It could be that a new
coordinate for every positional ”mark” would be required.
Keep in mind that methods described below also can be used for direct transformation
between Latitude / Longitude pairs, State Plane coordinates and vise versa.
The transform functions will typically be of most value where the survey is conducted
without GPS but the final map product requires GPS positioning. The process is as follows:
After the survey is completed the positions of the four corners of the survey area are
measured with the geodetic grade GPS, and coordinates are recorded as Latitude Longitude
or State Plane coordinates. (Alternatively, the positions of all survey marks (including end of
line marks) could be measured with the GPS). Note: State Plane coordinates can be easily
obtained with from Lat/Lon pairs using CORPSCON software program which is freely
79
available from the Army Corps of Engineers at:
(http://crunch.tec.army.mil/software/corpscon/corpscon.html).
Example 1: Point for Point Mark replacement. You have State Plane positions for all
survey “marks” (including end of the lines), you want the survey grid in Latitude/Longitude
or State Plane coordinates.
The first string is the header and is not used. The second line has two columns labeled
XGlobal and YGlobal, which contain State Plane (UTM or Lon/Lat) coordinates, two
columns with local (X-Y 858 coordinates) and columns for the line and mark numbers. The
program actually uses only the State Plane and Line Mark pairs. All marks must be present in
the file; omitted pairs are removed from the data in memory (a warning message is
displayed). Integer numbers must be used for Line/Mark pairs.
Note: When data is exported, magnetometer data is interpolated between existing marks.
To perform the transformation, load the G-858 STN or BIN file and select “Survey Setup
/ Replace positions” from the menu:
80
Click Yes and a dialog box opens for you to specify the replacement file name. Make sure
you check the “New coordinates are Longitude and Latitude” box at the bottom of the
dialog if you want to output position data in Lat/Lon, i.e., this box should be checked only if
the new coordinates are to be geographical.
The position screen will now be updated with the new positions and may look like this:
81
If there is an obvious misplacement, use “Undo” to return to the previous position data
display.
In this example all data positions are now in State Plane coordinate system.
Example 2: Local 858 X-Y positions and Global State Plane or GPS coordinates are
known for 4 or more corresponding points.
This is a more common case where the corners of the area are measured with GPS. Then
the new data position coordinates are simply recomputed from one coordinate system to
another. The program does not try to account for any 858 system positioning irregularities.
Four points known in both coordinate systems are sufficient to derive the transformation.
The new points need not necessarily coincide with any of 858 survey marks. The
requirements are that no more then two points can be connected with a straight line (we are
defining a new surface).
X a 1 xy a 2 x a3 y a 4
Y a 5 xy a 6 x a 7 y a 8
For example X, Y would be Global State Plane (or UTM, or new Lat/Lon positions) and x,
y would be in the local 858 X-Y coordinate system. If both coordinate systems are
accurately defined, a1 a5 are equal to zero and a2 a7 are equal to one. Coefficients a2 a7
describe scaling transformation; a6 a3 describe rotation transformation and a4 a8 describe
82
translation transformation between two coordinate systems. Non-linear terms a1 a5 describe
inconsistencies in the application of the model after the transformations have been
computed. For instance, if user defines local coordinate grid as 10 x 10 meter, and wants to
transfer it into State Plane, which is defined in feet, a2 a7 would be 3.28 (scale coefficient
between meters and feet). If he or she converts to UTM (defined in meters) then a2 a7 would
be 1. Furthermore, if in reality 10 x 10 meter grid defined by the user is not square and for
example is 9.6 x 10.2 meters, other terms in the equation for X, Y may also be used.
MagMap can use any number of points to compute this transformation. The recipe is to use
well-defined points with minimal positional errors. For instance in the picture below you
would not use points 1,2,3,4 if you know that point 3 is misplaced.
The local coordinates of point 3 (in x,y) will not show this misplacement. The correct
solution would be to use points 1,2,4,5. Point 3will still be misplaced after transformation,
and although nothing can be done here to correct its position, the rest of the grid will be
correctly located.
If the quality of reference points is not known, you can use some number of them to obtain
a least squares fit solution to minimize the errors. The fitting process is done automatically
when more than 4 points are used in the solution.
83
As before, open the Survey Setup and select Local System Transformation. The following
dialog is displayed:
At the beginning, this dialog box will be blank. To add points interactively, press the “Add
point…” button and fill in fields in the New Reference Point dialog:
In the Global section type geographical (State plane, UTM, etc) coordinates in the first two
fields. Type the local X-Y coordinates below in the Local dialog box. Press Ok. As soon as
four or more points are added, the program computes a solution. Solution coefficients are
displayed in the dialog box shown on previous page, along with computed errors (DX, DY
84
and Distance). Note that point 1 in the dialog box has a Distance mismatch of 21.86 m.
Clearly this point will cause an incorrect repositioning of other points and should be
removed from the computation. Points in the list can be deleted or modified by selecting the
appropriate strings and pressing, “Delete” or “Edit” buttons.
The program treats the entered positions in differing manners depending on whether the
“New coordinates are Longitude and Latitude” box is checked or not. If this box is left
unchecked, the values DX, DY, Distance are computed as simple Cartesian distances
between points in whatever units (feet or meters) the Global system employed. If the box is
checked, the program uses an ellipsoidal approximation and computes distances in meters
using Lat/Lon pairs. It is the user responsibility to select correct type of coordinates.
Selecting “Lat/Lon” when the actual data is in UTM or State Plane would result in huge
DX., DY and Distance error values.
If there are only four points, all errors are zero. Each time you add a point, modify or
remove a point, the transformation coefficients are recomputed, new solution coefficients
are displayed and DX, DY and Distance are updated. It is advisable to save this table into a
file, (use Save Points button) for future reference. File format is shown below:
The first line of the record is the header and is omitted. Note that the file format is
compatible with the format used to replace all the marks (Example 1). These two files are
interchangeable.
After the transformation is completed press “Ok” and the transformation information is
read into MagMap. All points on the screen are replaced with new positions. If there is an
obvious error, use “Undo” to revert to the earlier screen.
The general recipe in this case is that if a few points are known to have accurate positions,
use them. If you are not sure of the accuracy of your positions, use as many points as you
can to minimize errors through least square fit. Note that irregularities in 858 coordinates
cannot be removed with this method.
In this case there is a net of points where both local 858 and State Plane coordinates are
known. Reference points may or may not coincide with Line / marks. See example below:
85
The program computes the Affine Transform in each of the triangles:
X a1 x a2 y a3
Y a4 x a 5 y a6
Because the triangle defines a plane surface, the corners of the triangle (reference points
above) will be computed with zero error. However note that if the reference net is sparse
then the curvature of the survey lines will not be taken into account. If there are errors in the
reference points the survey lines could be curved.
If the State plane coordinates are known for each start and stop mark of each line, it is
possible to obtain the transformation by defining only a few reference points. Here is an
example:
86
If it is known that line 3 positions are suspect, you can use points 1A-1B, 2A-2B, 3B, 4A-
4B, 5A-5A. A good result will not be achieved using points 1A-1B, 3A-3B, 5A-5B.
MagMap2000 implements this function the same way as in “replacement of all marks”.
Select “Survey Setup / Coordinates Network”
If some of the marks are outside of the polygon defined by the triangle mesh, they will be
removed.
If the transformation appears to be miscomputed, use the “Undo” button to restore original
positions in case of error and check your reference point file for errors.
87
5. Other views of Readings
Magmap has the ability to display readings in different ways, depending upon what kind of survey
you have.
3-D plotting capability is included in Magmap. The purpose of the 3D displays is to give a quick
overview at the data. They are not intended for in depth data analysis or presentation but can be
printed.
Note: “Plot Stack Profile On Top Of Map” may not be available for some
types of surveys, specifically GPS surveys.
Selecting “Plot Mag Field in 2D Contours” and “Plot Mag Field in 3D Surface” will give a 2D contour
plot and 3D surface plot, respectively. If this is done from the survey map, both plots use the
current number of lines as X grids and a fixed 100 Y grids to create the mesh points. If this is
done from the GPS map, both plots create a fixed 100 x 100 mesh points. The grid value (Z) is
interpolated by two magnetic field readings that are closest to the grid point. The following
pictures are the 2D/3D examples for the survey on Stanford Test site:
88
2D Contour Plot:
3D Surface Plot:
By default, it will create 10 contour levels using maximum and minimum grid values. The value
and color represented at each level are listed to the right.
These plots use non-offset (GPS antenna to sensor) coordinates. In order to get a plot
that takes offsets, heading error removal, or de-striping operations into account, you will
need to export the file and then either view the output in a plotting program such as
Surfer or reload the file as a Surfer (*.dot) file into Magmap 2000.
89
Note: You can choose to plot any sensor (or pseudo gradient if available) by
right clicking on the map and selecting the option: “Plot sensors setup”. You
can then choose the sensor you wish to plot.
As it can be seen, most of this graph is dominated by one color. The other colors are sparingly
represented and this occurs because each color represents a roughly equal part of the total
magnetometer range. Since the range is divided up into equal parts, this technique will show a
bad contour plot if a drop-out exists, or some small noise exists, because then the range is a lot
larger and each color must span a much larger range. For instance, if we graphed the above and
had one drop-out in the data (that we didn’t filter), we would probably see an entirely red graph
because blue would represent a range of about 0 – 10000, green a range of about 20000 – 30000,
and red, most of the valid data taken.
We can solve this problem by using a technique of equalizing colors. This means that each color
is represented in the contour plot the same amount. In the graph below, we see this at work:
90
Note: The aspect ratio has not been adjusted.
This contour plot was done before any data was filtered or drop-outs removed. Thus, using
conventional techniques, this plot should have been entirely one color using the number of colors
available (for the same reason as discussed earlier). However, in this plot, blue spans about 400
nT whereas green spans about 8 nT – still, each color appears about the same amount on the
graph.
Thus, the technique can be summarized as: the range for each color is picked in such a way that
each color is represented the same amount on the graph.
Selecting “Set Min and Max Field” will allow the user to change the maximum and minimum Z
axis scaling. Type in the range of values you are interested in and press OK button when done.
This will repaint the window using the new defined maximum and minimum fields. The legend
to the right will be updated at the same time. Clicking on “Restore Default” button will reset
maximum and minimum fields to be default values.
Selecting “Plot Sensors Setup” will let user select the readings for a specific sensor to be plotted.
For the 2D/3D plot, it only allows one sensor to be plotted at a time. “Plot All Sensors” option
is inactive for the same reason. After selecting the new sensor and clicking on OK button, it will
repaint the window using readings from selected sensor. However, if user accesses “Plot Sensors
Setup” from Profile Plot Mag Field Plot, it allows multiple sensors selection. To do this, hold
down the CTRL key, and then click each sensor you want to select.
91
5.1.3.2. SET GRID SIZE
You can also set the number of gridding points used by clicking on the “Set Grids” option on the
menu (obtained by right clicking on the contour plot).
In this case, 100 grid points are used to define X and Y. (This is the default).
The X and Y width are calculated by taking the total X and Y distance (respectively), and dividing
by the number of grid points desired. Hence, if you want to have a different width, the number
of grid points will change accordingly. There is a maximum of 300 grid points in a given
direction.
One important tip to remember is that it will take longer to generate the graph if you select more
grid points. However, the resolution might be better if you have more grid points in the direction
you took your measurements. Faster computers are always better for tasks such as these.
92
5.1.4. User Interactions for 2D/3D plot
In 2D/3D plot, user can rotate, translate, scale, and zoom the graph . The procedure for each
operation is explained as follows:
Note: Some mouse drivers interfere with these dual button controls.
• Translation (2D/3D)
• Scaling (2D/3D)
• Zooming (2D/3D)
◊ Press “r”
Moving mouse cursor insider the contour/surface map will display the line number (LN), mark
number (MK), x and y coordinates, representative date/time, grid value (FIELD), and
representative reading number of the closest grid point at status bar.
93
5.2. Shaded Relief plot
The shaded relief plot is used to give a different view of the data, similar to what you might get via
a contour plot, but more customizable. It is produced by assuming there is a light source above a
3-D topographical map of the data. The location of the light source is user defined, and the
resulting shadows are graphed on a 2-D plot. You can decrease the effects that huge spikes have
in your graphed data by using scaling factors that work to reduce the influence of the spike. You
can also adjust the viewpoint by changing the relative location of the light source.
To produce a shaded relief plot, you need to first obtain a contour map. (See the section on
contour maps for instructions on how to do this). From here, right click on the mouse as shown
below.
After you select “Set Min and Max Field”, you should obtain the following dialog.
From here, press the “Shaded Relief” button to get a dialog that will let you set the parameters for
your shaded relief plot.
94
You need to select the option “Plot shaded relief map” to plot the map. Your map will be plotted
after you press “OK” and exit out of the ‘Vertical Scaling’ dialog.
The sun azimuth and elevation parameters allow you to set the position of the light source in
relation to the topographical map. The data scale setting allows you to vary the height
differences of the different parts of the topographical map. A larger value of the data scale will
favor large peaks, whereas smaller values will allow you to see smaller variations better.
An example of a shaded relief plot is shown below. Note, this is obtained from the same data
that produced the contour map shown earlier in this shaded relief discussion.
For certain surveys, it is useful to view the magnetometer data of different lines at the same
time in order to look for trends that occur between different lines. Since an anomaly
95
usually can be seen from adjacent lines, we can distinguish anomalies from noise by looking
at several lines. We can also locate objects such as pipelines that extend over many lines.
Magmap’s stack profile presents a way of looking at the magnetometer data of several lines
at the same time.
It scales the magnetometer data over each individual mapped survey line as seen in the
example below:
Note the similarity of the second line (X = 74) to the plot of our magnetometer readings:
96
We have zoomed in on the left line (X = 74) so that we can compare it with the magnetometer readings
below.
We can see that this early peak also appears above on the stack profile.
Currently Magmap supports handheld land-based surveys such as surveys done with the G-
856 and G-858 magnetometers.
97
Note: Stacked Profile is not available for GPS date.
From your mapped screen, right click on your mouse and select: “Plot Stack profile”.
• Plot sensors setup – allows you to choose which sensor(s) you would like to plot,
including the pseudo gradient, if available.
• Plot All Sensors – will automatically plot all sensors
• Don’t plot any Sensors – will remove all stack plots from the screen (good for getting
rid of your plotted profile).
• Zoom in Stack Profile – allows you to increase the scale of the magnetometer readings
(the positions will remain the same but the variations in the readings will be easier to
see). You can also press the “+” key to do this.
• Zoom out Stack Profile – allows you to decrease the scale of the magnetometer
readings. You can also press the “−“ key to zoom out.
In the case of the above survey, there is one sensor, so we choose “Plot Sensors Setup”. We
need to select the single sensor.
98
“Plot Sensors Setup” menu with single sensor selected
If we instead choose the option “Plot All sensors”, Magmap will attempt to graph all the
sensor readings on a single scale (e.g. the first sensor has readings averaging around 50000
gamma, whereas the second (non-existent) sensor has readings that are 0.) As a result, it
will be difficult to see any small features on the stack profile since there is now a scale that
has to span 50000 gammas compared to 5000 gammas when only one sensor is graphed.
Note: Due to the above, it is useful to de-spike the data and remove drop-outs
from the data before you try to create a stacked profile. Otherwise it will be
difficult to make out small features in your profile.
Your final result should look similar to the stack profile below:
99
5.3.3. Scaling a stacked profile
The scaling for a stacked profile works the same as for all windows. This will be discussed in
more detail in the following section.
The menu options that will immediately change the scale are
• Zoom in X and Y
• Zoom in X
100
• Zoom in Y
• Zoom out X
• Zoom out Y
• Window zoom. zoom in graph in the user selected window. To select zoom area, point
mouse into one of the corners, click and drag mouse to the opposite corner. Zoom rectangle
is shown with “rubber band” After releasing of mouse button selected rectangle is zoomed
to extend s of the window. Areas outside of selected rectangle are not shown. To bring
picture back to full extend press “Fit in window” button.
Once you select one of these menu options, your cursor will change to reflect the type of
zooming operation you are doing. You can then zoom in on a specific point by pressing the left
mouse button. This zooming feature will remain turned on until you select the menu option
again, or deselect it from the toolbar.
You can go back to a view of the full map by selecting “Fit In Window”.
You can also get rid of all scaling permanently by selecting “Always fit in Window”.
You can also zoom by using the toolbar buttons (pictured below)
• Zoom in X and Y
• Zoom in X
• Zoom out X
• Zoom in Y
• Zoom out Y
• Window Zoom: zoom in graph in the user selected window. To select zoom area, point
mouse into one of the corners, click and drag mouse to the opposite corner. Zoom rectangle
is shown with “rubber band” After releasing of mouse button selected rectangle is zoomed
to extend s of the window. Areas outside of selected rectangle are not shown. To bring
picture back to full extend press “Fit in window” button.
101
• Fit full map in window (but keep scaling options enabled)
If you want the scaling of X and Y to be the same, you can select the option: “Use true map
coordinates”. In this case, our original grid that was 30 x 100 would appear like the grid below:
102
You can enable grid lines by choosing “Grid Lines” from the “Scale” menu. By default, these are
not enabled. It should bring up the following dialog box:
Dialog box for “Grid lines” menu option (time series plot)
Both dialog boxes behave in a similar fashion so we will discuss the first dialog box -- this
dialog box allows you to set grid lines for X and Y and choose the scaling method. The
two methods available are:
• Auto grid interval (automatic scaling) -- This assumes you want a fixed number of grid
lines and will pick the scaling between each grid line so that you always see the same
number of lines on your grid, no matter how much you zoom.
103
• Manual interval – you pick a step that you want. The grid lines will remain the same
through all zoom factors. (e.g. the distance between them will scale so that they always
represent a constant value.)
Note: The time series plot mesh step is set in minutes. The time axis can also
have an irregular mesh because different survey lines were taken at different
times (you may not have a continuous set of readings over the entire time).
Here are examples that illustrate the two different types of scaling:
In each example, we show the initial graph, and then zoom in on the lower left hand corner
(corresponding to the coordinate 0,0). Both examples show the initial shot and two
additional shots corresponding to two successive zooms.
The first example deals with automatic scaling. When you look at all three shots, notice
that the number of scale lines never changes.
104
Initial map with “Auto scale grid lines” selected.
105
From here, we can observe that the number of white grid lines always remains the same, no
matter how many times we zoom. The spacing between the grid lines does change.
This is different when we manually scale. The parameters are shown below:
106
Graph obtained from zooming in on (0,0). (green starting square)
Note: As one can see, the total number of lines always remains the same
when there is a manual scale (grid is fixed) and the number of lines shown on
a given view remains the same when there is an automatic scale.
It is recommended that you use grid lines before printing. Note that the appearance on the paper greatly
differs from what you see on the screen.
5.7. Marking anomalies for further study
In the course of a typical Magmap analysis, you usually will use at least two windows – one
showing the positions involved in your survey, and the second showing a profile of the readings
over time. It is often useful to be able to correlate an interesting feature seen on the reading
profile with the positions given in another window. This can be done by the use of flags.
107
The flags placed on the profile will simultaneously appear on the position map and on the
contour map, allowing the user to place the cursor on the flag on the map, for instance, and
easily identify the anomaly’s position. This feature works in either X,Y or Lat/Long spatial
mapping configurations.
To set a flag, you need to have an X,Y map (default screen for G-858 surveys), or a map
with GPS positions. From here, you may click in an empty area (but inside the survey
grid area) and select “Plot Mag Field”.
This will bring up a profile plot of the readings.
Once you have the profile plot, you may right click on part of the map that shows the
reading that you would like to locate on the position map.
The above screen shows the selection of the flag using the menu obtained by right
clicking on the spot where you would like to place the flag. Once you have placed a
flag, you should see the following.
108
Figure 3 – position map with corresponding flag
Note that the ‘1’ shown on the magnetic profile corresponds to the ‘1’ shown on the
position plot.
If you want to make a flag on the GPS plot, you need to bring up a magnetic profile
corresponding to the GPS map. In the above example, we do this by selecting the GPS
menu and then selecting “Draw new map using GPS data and the features of a Maglog
NT survey”. You should then see a GPS map and the corresponding magnetometer data.
Note: The magnetometer data associated with the GPS positions is treated
differently than the magnetometer data associated with the G-858 flags and
positions. Thus, you need to choose the appropriate magnetic profile that
corresponds to the map you want to flag (if you want to use the G-858 map,
then click on that map to get the appropriate data, and if you want to use the
GPS map, use the data associated with that map). This also applies to any
109
data editing operations – if you edit magnetometer data associated with the
standard G-858 position map, you will need to repeat the editing operations if
you wish to use the GPS map with the data. Since most users only use one
type of map, this shouldn’t cause any problems.
You can have as many flags as you wish. The flags should also be visible on any
contour plots you generate. It is necessary to close a contour plot and reopen it to see
the added flags.
Future additions to this program may include a table output of the anomaly positions
110
6. Base station and survey readings
The Earth’s magnetic field varies over time periods on the order of seconds to hours,
similar to survey periods. The variations are caused by internal changes in the Earth’s
magnetic field and the distorting effects of the solar wind. The distortions result in
micropulsations and magnetic storms that occur at any time, and during the daylight
hours, there are diurnal or daily magnetic field changes.
This section describes the MagMap2000 procedure for making diurnal and
micropulsation corrections to data collected with the Geometrics G-858 or G-856AX
magnetometers, with either of the same magnetometers set-up as the base station (all
combinations are acceptable).
The plotting of readings for a segment, line, or grid is similar to the plotting of base
station readings. Therefore, only the base station case will be described in the following
paragraphs.
From these readings, MagMap will determine which readings it can use. In the best case, the
base station takes readings at the same rate as the magnetometer. This means that there will be a
base station reading that corresponds to every roving magnetometer reading. In this case, the
diurnally corrected reading can be calculated by the following formula:
In other cases where fewer base station readings are taken, Magmap will linearly time interpolate
between known readings for the intermediate readings. The diurnally corrected reading will then
be calculated in much the same way as above.
111
Note: In cases where the two magnetometers are not set to the same time, you
can add a given time to the base station clock to make it correspond to the
roving magnetometer clock.
If you do not have a proper base station file open, Magmap will inform you of this at export
time. You will see the message: “No valid base station files are open…” This message will not appear
if you have a valid base station file open – e.g. a base station file where the roving magnetometer
and base station times overlap.
The quality of the correction may be degraded by a mismatch in times between the base station
and survey data, a time mismatch combined with highly variable base station data, and noisy or
spiky base station data. Although it is optimal to synchronize the instrument clocks and place
the base station in a magnetically quiet location for the survey, this is not always logistically
feasible. MagMap2000 provides methods to deal with of these potential problems – see the
following section.
The resulting screen allows you to adjust the base station’s clock to correspond with the
magnetometer clock. Additionally you can choose to add a constant value to each outputted
reading.
Figure 64
In the “bias to add to diurnal magnetic values” just type any positive or negative value you want
to add to every base station magnetic field.
The “Diurnal” column in the XYZ file consists of the difference between the sensor reading and
the base station reading, plus the offset entered in this dialog box. This is in case you wish the
diurnal values to be in the 50,000 nT range, rather than relatively small positive and negative
values.
112
6.2.2. Adding a bias to diurnal clock values
A time mismatch between the base station and survey data will cause misalignment of the data
values and cause subtraction of the wrong base station value during the correction. One can
imagine that, if the base station data curve has large variations, a mismatch would yield even
larger magnitude inaccuracies than if a smoothly varying curve was used. In turn, if the base
station data curve is essentially flat (constant) during the survey period, there would be little
concern about the effects of a time mismatch.
One may use the bias shift option to apply the shift needed to match the start times of the base
station and survey data. The period of the base station data file should be equal to or greater
than that of the survey data.
In the “bias to add to diurnal clock” just type any positive or negative value in seconds that you
want to add to the time of all acquired magnetic field readings of the base station. This will allow
for you to correct for errors due to the clock of the base station being different from the clock of
the magnetometer.
Smoothing of the base station data with third degree spline filter before applying the diurnal
correction routine can alleviate this problem. Replacing the original base station readings with
their continuous equivalent effectively moves the noise into a very low frequency domain. In this
manner local magnetic anomalies are not distorted with noisy base station readings. It is
recommended that the user despike the data prior to using the spline-smoothing filter.
This option is available under Filter / Smooth base station readings menu:
113
After this option is selected, the following dialog appears on the screen:
A larger “smoothing degree” value provides more smoothing. The user should try different values by
pressing the Preview button and examining the results as it is shown on example below:
The black line represents a preview of the smoothed values. Note that the minimum amount of
smoothing that gives a reasonable result should be used. Employing too high a smoothing factor will
tend to remove small micropulsations and other rapid field variations that may be real. Some
experimentation will be required to get a good result, but in general a spline degree setting of 5 to 10
will be sufficient.
114
After a satisfactory result has been reached, the “Smooth base station data” dialog box should be
recalled and the Accept button clicked. The original base station readings are then replaced with the
smoothed values:
Now the user can proceed with the field reading export as usual. Smoothed base station data is used
instead of original base station readings to perform diurnal correction.
“Remove” button erases results of smoothing from the plot without accepting them.
Figure 65
In the previous graph (Figure 65), the Y axis represents the amplitude of the magnetic field in nT
and the X axis the reading number.
The time of the first and last reading is shown on the X axes.
115
When you move the mouse cursor in this window additional information will be shown on the
status bar.
The first item shows the magnetic field amplitude, the second item is the reading number in the
survey, and the last one is the time of the reading.
The base station scale tool bar is similar to the map scale tool bar described in paragraph 4.1,
page 49 except for the last button, “True map coordinates,” which is not available and grayed
out.
Figure 66
By default, only the readings for the first sensor are plotted. You can plot if you wish the
readings of the second sensors or simultaneously the readings of both sensors.
You can also change the scaling parameters by selecting “Set Min and Max” as shown below in
Figure 67.
Figure 67
You can set the maximum and minimum magnetic field value to be plotted in the first two items
of the dialog box.
You can restore the default value by clicking the left mouse button on “Restore Default”.
116
6.5. Base station file formats
The base station file may have any name with the .BIN extension for a binary file and a .STN
extension for an ASCII file. The .BIN file is a raw data dump and not stored in a user-readable
format; do not modify files of this type. Depending on the magnetometer used for the base
station, the .STN data formats would appear as follows.
The values are line number, Julian day, time, station number, and magnetic field reading to one-
tenth (decimal not shown).
The values are an internal tag number indicating type of data, first magnetic field reading, second
magnetic field reading (zero unless using a second sensor), time, date, and a second internal tag
number.
The base station data file time period must be equal to or greater than that of the survey data file.
At a minimum, the file must contain two data points at times equal to or greater than the start
and stop times of the survey data file.
It is possible to use SURFER ™ file for base station readings. File should consist of at least three
columns: readings, date and time.
117
7. Processing magnetic field readings
7.1. Introduction
Despiking is a technique used for removing readings known to be incorrect. These spikes are usually
identified because they differ significantly from the adjacent readings. They can be caused by operating
the sensor near the dead zone, under high tension lines, and or over sites where there is a high
accumulation of small magnetic surface clutter and the operator is carrying the sensor just a few inches
above the ground surface.
However, often, anomalies can be mistaken for spikes, so we need to develop a scheme that allows us to
distinguish anomalies from spikes.
One criterion that we can use to distinguish multi-point spikes from anomalies is the number of points
they span. It is reasonable that we will not see a 10000 gamma anomaly for .1 seconds (e.g. about one
point for a magnetometer taking readings at a rate of ten times per second). However, if we saw this jump
last for a long time (hundreds of readings), we would probably believe that it is an anomaly rather than an
incidental spike. A lot of this is also dependent on the speed of the operator.
We can also use the height of the anomaly to characterize an anomaly. A jump of one gamma might be
considered to be noise, whereas a jump of 5000 gamma might be an anomaly or a spike.
• identifying and removing spikes (non zero anomalous readings over a short duration)
When we remove drop-outs, we try to remove the effects due to being near a dead-zone or near a very
large ferrous mass. These are zero readings that can show up in the data for short to long periods of time
– usually showing that the magnetometer was operating under conditions where good readings were not
possible.
Since these readings are not valid in our data and can overwhelm our efforts to see smaller anomalies, we
remove them. We can choose between linearly interpolating between bad readings and removing the
readings from the file without linearly interpolating them. Because most gridding programs will grid over
the missing data, we recommend removing them. For more information, please see the following
discussion on removing drop-outs.
To identify spikes, we suggest zooming in on potential anomalies and then counting the number of
readings to see if they are really anomalies. To zoom in on a section of the map, de-select the “AUTO”
button and click on the “X+” tool-bar button to allow zooming in the X direction. Once you are at a
scale that allows you to distinguish the anomaly, place the cursor on the left and right at the beginning and
end of the anomaly and note the reading numbers on the information bar at the bottom of the plot
display. You can get the number of readings involved in the anomaly by subtracting the two numbers.
In general, we would think that any anomaly greater than one second (10 readings or less at 0.1 second per
118
reading) or perhaps less than one second depending on walking or traverse speeds, would be determined
to be a real magnetic field shift and not some artifact of the survey process.
Once a true “spike” is identified, it becomes necessary to define the limits of the spike editing process
such that actual data is not modified in any way. We do this by setting the spike edit parameters by
clicking on Filter and then Despike Setup. Default values are 5000 nT and 1 reading – a spike by anyone’s
definition!
Depending upon the type of spike identified, some iterative trial and error process may be required to
remove the data spikes. Using the number of readings identified in the first part of the investigation
process above, and using the maximum amplitude of the spike are good values with which to start. If the
application of the spike filter has no effect on the profile plot, then a larger number of readings and/or a
lower maximum anomaly level may have to be applied. Just click UNDO to remove the last spike edit
process and then apply the new parameters.
When spike editing, keep in mind that all data will be edited that meet the criteria, not just the anomaly
you are investigating. Sometimes removing the largest, very sharp anomalies can enhance your ability to
see smaller objects in the survey area. The process of removing spikes is discussed in more detail later this
chapter.
“Remove Drop Outs” will automatically delete those readings with exactly zero value and adjacent readings
around those zero values (+/- one reading). This will cause no 1D plotting for those readings and no
output for them when exporting.
The difference between “Linearly Interpolate Drop Outs” and “Remove Drop Outs” is “Linearly
Interpolate Drop Outs” will interpolate those zero readings by good adjacent readings. In other words,
those zero readings and adjacent +/- one reading are replaced by interpolated values. Therefore, those
readings will still be plotted or exported by the new values after filtering. “Remove Drop Outs” is usually
the preferred method because most gridding programs simply grid over the missing data.
Please note that “Linear Interpolate Drop Outs” cannot be run if “Remove Drop Outs” has been run
first because there is no meaning to do so.
119
7.4. Removing spikes
Magmap allows the user to define the maximum number of readings allowed in a spike and the height of
the spike.
This allows the user to enter the peak threshold and select the maximum number of readings (width of an
anomaly, in other words) which could be considered to represent a spike. When done, clicking on the
Apply button will automatically find those spikes that match the criteria and filter them out. If the reading
window is open, it will update the display automatically. Therefore, the user can easily check how well the
despike filter works using the current settings. If the effect was not what was desired, the user can click
on the Undo button to reverse the filtering. Clicking on the Close button will close this setup dialog box
with the new settings saved in memory. Therefore, next time user selects “Despike” from Filter menu, it
will always use the current despike setting to filter readings. The default value for peak threshold is 5000
and maximum spike readings is 1.
Use “despiking” carefully for (non-continuous) discrete surveys. Readings may have been acquired far
apart and therefore may look like spikes to the algorithm.
Note: High amplitude noise variations usually occur over a very short time span. Thus, if you have an
anomalous set (more than ten) readings, there is a good chance that this is not noise.
120
7.5. Destriping
7.5.1. Introduction
Destriping is a new feature offered in MagMap 2000 which offers the user the ability to remove
coherent heading error offsets from data collected in parallel line mode. In addition, the program
can perform a “smart” high frequency filter along the line to remove the small sinusoidal noise
associated with “stride” or quantization processes inherent in the collection of the data.
We define heading errors as errors that are dependent on the direction of travel. When a survey
is performed as a series of bi-directional parallel lines, these errors can be very noticeable and
significant. Heading errors can be created by any number of sources, but the primary sources are
heading error of the sensing system itself (often in the 1 to 2 nT range), offsets due to carrying
some ferromagnetic object on the operator, or from the logging console or some other sensor
device such as a GPS receiver or antenna being too close to the sensor. Only since the advent of
very high sensitivity optically pumped Cesium vapor magnetometers have these errors been
visible and the necessity to remove them become apparent.
Heading errors are traditionally removed by using a technique that attempts to find constant
magnetic offsets for each line traveled. Since the causes of heading errors are usually magnetic
objects that change orientation as the user changes orientation, in cases of heading error, there
should be a constant magnetic field offset added to or subtracted from the magnetic field for each
direction of the user. Thus, we attempt to find a technique that removes this constant magnetic
field offset from each line. The line-to-line offset minimization uses a local min/max technique
that searches down every line for local minima (or maxima) and seeks through an iterative
process to minimize this offset on a line-to-line basis. The technique we use therefore has a
minimum distortion effect on large anomalies while performing quite well in the flat areas of the
field where heading error offsets are most apparent. It assumes that the lines are run straight and
parallel.
The input to the Destiping program is the *.dat file exported from MagMap after all filtering and data positioning
edits have been made. The output of Destriping is a *.dat file and therefore can be called back into MagMap for
viewing. However, better visual results will result if the data is gridded and displayed in third party program like
Surfer or Geosoft. When the routine runs, it gives diagnostic data about the destriping process which will give the
user important information about the success of the process and whether an additional iteration or varied parameter
is required.
Destriping like GPS Interpolator (for large marine or airborne files), is accessed through the File
menu selection. Click on File, then Destripe Data to access a dialogue box (shown below). The
pertinent parts are:
121
• Name of input file. This is probably a dat file that was created using Magmap previously. It
should have space separated columns with X, Y, and field values (at minimum).
Additionally lines numbers are useful, but not mandatory.
• Name of output file. The default name is output.dat, so it is advised to select a new name
that is more related to the name of your file.
• Name for statistics file – this is enabled if you have selected the option, “Output Statistics
File”. The statistics file will report the statistical variance and the maximum offsets of the
destriping process.
This screen allows you to specify how to read the input file. Since destriping uses a columnar
approach to interpreting the data format, you need to tell the program which columns hold the
pertinent data. Usually the X and Y columns are listed, as well as one of the magnetometer or
other data file columns, and the line number column. If you have collected multiple sensor data,
then you will have the option of destriping either sensor or the gradient columns.
122
Input file parameters screen
By default, the X, Y, Field, and line columns are set to the first four columns present in the file.
You will probably need to change the settings to reflect the correct columns used. Click on the
down arrow to select the proper data column.
123
• Line column - Destriping uses the line number column to determine which lines it is
working on. Select the proper column for the line number. If the output file does not have a
line number included, you will need to select the option, “Use distance interval to split the
lines”.
• ‘Use Distance Interval to split the lines’ checkbox – This is used when there is no line number stored in the
input file. In this case, you need to specify a distance interval to specify when a new line is started. This
distance interval is defined as the maximum distance between the positions of two adjacent readings. E.g., if
you do a land survey and have line spacings of 5 meters, it is probably safe to put a value of 2 or higher for this
reading. If the value is too large, it will not notice when you have shifted lines, and if the value is too small, it
may assume you have started a new line, when in fact you have just walked a little faster, or had a GPS error
that made the distance larger.
•
7.5.2.3. DESTRIPING PARAMETERS DIALOG
Next, Magmap will bring the data into the destriping program and report the statistical findings.
You should see a screen such as the one below:
• X min, X max, Y min, Y max – These are respectively, the minimum and maximum X and Y
values found during a search of the file
• X cells, Y cells – these are the number of gridding cells calculated based on the interval
specified below.
• Interval – this should approximate the average distance between adjacent points. For a series
of data points equidistant, there should be no performance or accuracy gain by specifying a
smaller interval. However, for points which are not equidistant, a smaller interval may help
with accuracy. Usually the default value works well for most situations.
• Destriping azimuth: This is the angle between ‘north’ in your gridding system, and a line
drawn perpendicular to your line. Usually the default works well in this case.
124
• Iterations: This is the number of times the program will increase/decrease the iteration step.
• Iteration step: this is amount the program will try to vary the field for a particular line on a
given iteration.
• Max change: This is the maximum magnetic field amount that a line will be changed.
The two parameters, ‘iterations’ and ‘iteration step’ are best explained by understanding the
method by which the destriping routine works. During the destriping process, we compare the
readings across a number of lines that are taken perpendicular to the direction of travel. If we
compare all the magnetic readings across a given perpendicular line, there will be a point on a
survey line at which a minimum value occurs, and a point on a survey line at which a maximum
value occurs. For each survey line, we make a total of all the minimums and maximums
occurring on a perpendicular line. In destriping, usually one or more lines will have an abnormal
amount of maximums, and another one will have an abnormal amount of minimums. By
changing the field values of each line by a small constant value, we try to have an even
distribution of maximums and minimums for each line, and hence destripe the data. We take
small steps (iteration steps) because we don’t want to overshoot – (e.g. add too much to a
minimum or subtract too much from a maximum) but we also need to take enough steps to make
a difference (iterations). Hence, a good result involves picking a small enough iteration step to
not overshoot the correct value, and additionally doing enough iterations.
This dialog is used to help us smooth the data to remove high frequency noise (“stride noise”).
Typically, stride noise is very small (on the order of 0.2nt) and we do not wish to effect any
actual anomaly data, so we allow an upper limit on what we will consider stride noise. The
default for this parameter is 1nT. We consider any anomaly of greater than 1nT that is coherent
to be real and we do not want to remove it. We recommend initially leaving this setting at 1nT.
If for instance the stride noise is larger due to some magnetic object on the survey personnel,
then this value should be increased. This dialog box is shown below.
125
A low pass filter is used to filter the stride noise. Additionally, to keep from filtering out actual anomalies, we use a
setting “Dynamic range for smoothing”. This setting is used as a guideline when trying to distinguish anomalies
from noise.
The ‘Run Destriping’ dialog box gives us to options -- either running Destriping in a test mode
or to actually apply the filters and write out the files. Clicking on “Run Destriping” will run the
program without writing out the file so that the users can see the results of the parameter
selection. In particular, we want to see the maximum offsets corrected and the average and
maximum disbalance after the line-to-line balancing is applied. If the values are still too high,
then you are given the option to go back to the earlier screens to select different parameters.
Typically we like to see less than 5 nT maximum correction and a maximum disbalance after
destriping of perhaps 5%. When you are happy with the results, press Finish and the actual
destriping process with file output will occur.
126
Bring the resultant file into Surfer or Geosoft and grid and display the data (we recommend
shaded relief for maximum contrast and dynamic range display) and observe the results
compared to the original *.dat file. If the results are insufficient, we recommend running the
previous destriping output file through the process a second time. This often generates a
smoother result with the downside of removing more of the high frequency content of the data
file.
We believe that judicious use of the spike editing and destriping routines in MagMap 2000 will
afford the user powerful data correction capabilities. Some experimentation will be necessary to
get optimum results.
Shown below is data taken and mapped without the destriping routine and with the destriping
routine run once. Note that the stripes on the original data are much more pronounced.
127
Shaded relief after one pass Destriping higher amplitude Original Data at higher amplitude
Shown below is data taken after destriping has been done twice.
128
Data twice destriped with higher amplitude settings Original data at higher amplitude
129
7.6. Pseudo gradient
During loading of G-858 compressed or ASCII files, Magmap computes a psuedo gradient – the
difference between sensors one and two. This appears in the program as a third sensor (e.g. a G-858
magnetometer would appear to have three sensors as opposed to two). All kinds of processing available
for total magnetic field are therefore available for the pseudo gradient, including 2D and 3D maps, line
graphs, and stack profiles. For example, the plot selection plot sensors setup dialog looks like the
following:
This ‘pseudo gradient’ option is also available when selecting what to plot for the contour plot and the 3D
plot.
130
Contour map after one pass Destriping Original contoured data
131
Shaded relief after one pass Destriping higher amplitude Original Data at higher amplitude
132
8. GPS
8.1. Overview
Magmap will only import GPS data that has been stored with the following format:
NMEA output message format with a GGA message type.
• Draw new map using GPS and the features of a regular 858 survey
• Draw new map using GPS and the features of an airborne survey
If you select either item, a new map will be drawn (as shown in Figure 69) using only
GPS information. Each flag is a GPS position fix.
Every latitude and longitude in the data stream is converted from degrees, minutes,
seconds to decimal degrees within a range of values from -180.00000 degrees to
180.00000 degrees.
When you export your survey to Surfer or Geosoft Analysis and Plotting Software,
magnetic readings will be linearly interpolated between GPS locations in latitude and
longitude decimal degrees. UTM conversion to meters is optional see section 8.3..
The user interface for a GPS map window is about the same as the regular map
window using flag locations. Many of the editing functions that change the
coordinates of a location are not available for a GPS map. You may move a
(presumably bad) GPS position to a new location. You also can delete bad GPS
positions, delete readings and plot readings.
133
degrees to be faithfully represented in the Surfer map grid
annotation (in previous versions you need to use UTM coordinates
to get accurate positioning).
Figure 68
There are small differences in the way the program will behave depending upon with
which option you choose to plot the GPS window. You should choose the first
option (“Draw new map using GPS and the features of a regular 858 survey”) if you don’t
need to specify additional offset (antenna to sensor) position points to be calculated.
It also has an export menu that is simpler to use. The second option (“Draw new map
using GPS and the features of an airborne survey”) will allow you to specify extra positions
relative to the GPS antenna to be calculated. Calculating additional positions will be
discussed later in this chapter.
134
8.3. Enabling/Disabling UTM Coordinates
Magmap has the facility to do a transformation from geographical to UTM coordinates. This can be done
from within the menu option : “UTM Setup”. You should see the following dialog box that will allow you
to set specific parameters related to your survey.
• Enable UTM Transformation: This allows you to turn the UTM transformation on and
off. When it is off, Magmap will use geographical coordinates. Magmap will use the given
parameters to calculate UTM parameters when it is checked.
• Ellipsoid name: This is the type of ellipsoid you use to approximate the earth. Since the
earth is not a perfect sphere some ellipses may better approximate the region you are
surveying than others.
• Autoset projection parameters: This is useful when you don’t know what the prime
meridian you are using is, or when you data where your prime meridian may change. This will
set the prime meridian and scale factor automatically for you.
• Prime Meridian: If you don’t press “Autoset projection parameters”, you should fill
this in. This value has a large impact on the UTM coordinates generated due to the
non-uniformity of an ellipse. It can be quickly estimated as the value of the longitude
but you should look it up below for a better calculation of your positions.
135
The Universal Transverse Mercator (UTM) Coordinate System uses zone codes instead
of specific projection parameters. The table that follows lists UTM zone codes as used
by GCTPc Projection Transformation Package.
136
• Northing / Easting: These are constant values that you would like to add to your
UTM coordinates. UTM coordinates are usually large numbers on the order of a
million meters but your survey may only be looking at thousands of meters. This
allows you to simplify your final numbers by adding or subtracting a large offset.
Occasionally you may see a warning such as the one below when you enter the UTM setup screen:
In most cases, it is a good idea to have the program automatically set the UTM parameters since
a bad value for the prime meridian can seriously affect how your data appears on the screen.
137
8.4. GPS Offset
This feature is available for surveys where you earlier selected “Draw new map using GPS data and the features
of an airborne survey”. It allows you to specify one or more offsets in relation to the GPS that you want to
be calculated and exported. This is needed when you have a GPS but your magnetometers are not
located at the same place as your position point was taken. In the case of a handheld magnetometer
(G−858) survey, the GPS antenna is always a fixed offset away from your magnetic sensor. In the case of
a marine magnetometer (G-880 or G-881), the sensor is towed on a cable of fixed length. In order to
produce a final data set without positional errors, distance offsets from the GPS antenna to the sensor
must be calculated. These values are then used to recalculate the actual sensor position.
The position points that you specify in the offset screen will only be calculated when you export your data.
The recalculated sensor positions can be viewed in Magmap, but must be exported and read back in.
You can set offsets by pressing “GPS Offset”. This should bring up the following dialog box:
• Points to look forward: Magmap needs to establish the direction the ship or person is
travelling. This tells it how many points forward from the current point it should use to
establish the direction of travel.
138
• Points to look backward: This tells Magmap how many points backward it should use to
establish the direction of travel.
Note: The number of points you use to calculate the heading can be very important.
We recommend that you choose between 2 and 5 points for both options. If you are
moving in a straight line, more points will tend to give a better answer. If you are
doing marine survey and making a lot of turns, it is better to choose a small number of
points.
• GPS clock bias: This is a number that you would like added to your GPS clock time. This
is one way you can fix your positions if you have a constant offset. In one case, if the GPS is
located 1.6 meters behind the magnetometer and we are walking at 1 m/s, we could use a
GPS bias of –1.6 seconds to fix the positions. (Distance = Velocity * Time).
• Multiple Offsets setup: This button allows you to set more position points. The dialog box
is shown below:
From this dialog box, you can add and delete extra position points to be calculated. You can also specify
whether to use a dragging algorithm to calculate positions or not.
• X: This is the horizontal offset in meters (see diagram in “GPS Offset Setup”) of your
position from your GPS. Assume the GPS is at the coordinates (0,0). Then, from here, you
139
need to specify how far right (positive X) or left (negative X) your position is from the GPS.
Assume that you are facing in your direction of travel.
• Y: This is the offset in meters parallel to the direction of travel of your position from the
GPS. If your position is in front of the GPS, this will be represented by a positive Y. If it is
located behind the GPS, this offset should be negative.
• Z: This is the vertical offset in meters of your position from the GPS. If the GPS is above
the position point you would like Magmap to calculate, this offset should be negative. If it is
below the point you would like to calculate, this offset should be positive.
Note: GPS of vertical positions are typically two or three times less accurate
than X, Y accuracies.
• Name: This is an identifier by which this position will be known. When you are exporting,
you can choose to export the position by selecting this name. Additionally, this name will be
printed in the heading of the export file, should you choose to export the position.
• “Use Dragging”: This is used for marine surveys to give better position calculations. If you
check this and select a positive dragging distance in meters, Magmap will calculate the
position you specify using previous positions of the boat to give a more accurate boat
positions. Typically this is not used in land surveys.
• If you select this option, make sure that you don’t include your cable length when you specify
your distance, “Y”. E.g. If you have a single 880 magnetometer that you tow with fifty
meters of cable, you should specify a dragging distance of 50 and a Y offset of 0, assuming
that the GPS is located at your tow point at the rear of the vessel. If your GPS is not at the
same position as the tow point, you would enter the Y offset for the GPS with respect to the
tow point.
• Add: This will add a position point to the list. It will use the current parameters filled in
under X, Y, Z, and Name.
• Delete: This will delete the selected position point from the list.
Note: You cannot edit a position point once it has been created. If you want to
modify an existing position point, select it, delete it, and type in your correct offsets.
Then press “Add”.
140
This feature allows you to smooth your GPS data and get rid of spikes in your position. After selecting
the menu option : “Smooth positions”, you should see a menu such as the one below that lists all of the lines
you have in your survey. You will also see series of parameters you can set to manipulate the magnitude
and type of smoothing.
This dialog box gives us a lot of information. It lists all of the lines Magmap has identified and allows you
to set individual smoothing parameters for each line.
You can click on any of the parameters (except for “Line #”) and change the value. The parameters are:
141
• Line # -- this is the number of the line you will be smoothing.
• Spl. Smooth 1: This gives you the option of turning smoothing on or off for a given line.
• Despike: This allows you to choose whether you would like to despike before smoothing or
not.
• Jump: This controls the distance that is considered a “spike”. The program will calculate a
smoothing curve. It will then calculate the distance between these smoothed points and your
original points. If any of these distances are greater than the distance you specified in “jump”,
the point will be replaced by the point on the smoothed curve. All other points falling within
the “jump” distance of the smoothed curve will be left alone.
• Parameter: This controls the degree of smoothing for the second smoothing.
Note: You apply any change in parameters (except Line Number) by right clicking on
the gray heading buttons. E.g. to change the Despike option, we would press
“Despike”. This will then give a menu that will allow you to change the value and
apply the option to all lines.
You can also change the parameters for an individual line by selecting the line and pressing “Change Line
Parameters”. Below we have chosen to edit the parameters for line 4.
142
We have changed the first smoothing parameter from four to six (higher degree of smoothing) and turned
the second smoothing with the parameter to four. After pressing “OK”, the changes are reflected in the
dialog box.
To smooth the data, we press “Smooth”. This should draw a series of yellow lines that will show what your
smoothed data will look like if you apply the changes.
Note: No changes will be applied until you press “Close and Apply”. If you don’t
like the results of your smoothing, you can exit the screen by pressing “Close”, or you
can change the smoothing parameters and press “Smooth” again.
For the example data shown before, pressing “Smooth” results in the following screen:
143
As you can see, line four (fifth from left) appears straighter than the other lines because it has a higher
smoothing parameter and the second degree of smoothing is turned on.
Once applied, these changes cannot be undone. You will have to reload your file if you want to
undo your changes.
144
8.5.1. Recommended Smoothing Procedure
We recommend the following for optimum results:
1) Turn off despiking and your second smoothing. Figure out your first smoothing parameter by trying
lower numbers first, and higher numbers later. You want your smoothing to be strong enough to
minimize any effects of spikes. From here, there are two things you need to think about:
a) Your data has obvious spikes that you would like to take out.
• In this case, it is okay if your curve doesn’t approximate most of the actual data. You should
aim to get a curve that takes out the effects of the spike. The data that falls within the
acceptable “jump” limits will not be changed. This data will only be changed on the second
smoothing. Go to step 2.
• Here you should try to make your smoothed curve approximate your data as much as you
can. Try to determine what is due to error and what is true positional reporting, and use small
smoothing parameters to make your final path look reasonable. You probably won’t need to
do any additional despiking or smoothing (e.g. ignore Steps 2 and 3 if you get a reasonable
curve).
2) Despike your data: After you have your first smoothed curve, choose a “jump” parameter that will
remove the spikes.
3) Smooth your data a second time: If your path still doesn’t look acceptable, apply smoothing a second
time … probably to a smaller degree than the first. In most cases you will be doing this after
despiking so you will want to use a smaller degree of smoothing because your first smoothing had to
take out the effects of the spike.
If the technique described above doesn’t work, experiment with the parameters and see if you can get a
path that looks reasonable.
If all else fails, you can always manually delete certain problem points and proceed from there.
145
8.5.3. References on splines
For more information about splines, you can refer to the following:
The Theory of Splines and their Applications: Ahlberg, J.H., E.N.Nilson, Walsh, J.L, Academic Press,
New York, 1967
• This book is out of print but can be found in some online used book stores (e.g.
BarnesAndNoble.Com)
146
8.6. Troubleshooting
8.6.1. Map appears incorrect when using geographical coordinates.
Occasionally there is a problem where an offset was put into the UTM setup, UTM was disabled, but the
offset still remains. Sometimes this will affect the map that is created. The easiest way to get rid of this
problem is to select UTM, set both offsets to zero and select “OK”. This will redraw the map in UTM
coordinates. You can then go back into the UTM setup and deselect this item and the map should be
redrawn correctly. If this doesn’t work, try exiting the program and reloading your file.
In UTM setup, set your northing and easting offsets to 0. Press “OK”.
On your GPS map, look for the two numbers on the bottom left of your screen and write them down.
Then go into your UTM setup again and input the negative of the first number for your “northing”
parameter, and enter the negative of the second number for your “easting” parameter. In the above case,
I would use northing = -500000, easting = -6000000. Press “OK” to accept the changes.
147
9. Other devices
9.1. Metal Mapper
The metal mapper is a device which can be connected to a G-858 console instead of the magnetic sensor.
From a technical point of view, there is no difference between data collected with a true magnetic sensor,
and data collected with an EM device – this means that Magmap will not be able to distinguish between
the two different types of data and it is up to the user to tell the program what to do.
1) Download data from the G-858 console as it is described for the magnetometer. However, check
the button “Download only, decompress later.” The result will be a binary file on your hard disk.
2) Go to “Open” (on the file menu). Select “Metal mapper binary files” in the dialog box. If you
do not do it, Magmap will assume that your file is a magnetometer file.
3) Load the file as usual. Note that the menu item, “Metal mapper” under “Survey setup” is available.
Like the magnetometer, the metal mapper has two sensors, but it does not have a pseudo gradient. As
soon as you have worked once with metal mapper files, binary files that are in your history list are treated
as metal mapper files. This changes after you work with the magnetometer again.
1) For each line, calculate the median value for each sensor and subtract this from the data.
2) Create a set of new readings by summing the absolute values of both metal mapper channels.
These new values will replace the readings for sensor one and the readings for sensor two will be
equal to 0.
To do these transformations, go to “Survey setup / Metal mapper”. The following dialog should appear on
the screen:
148
Choose “Remove median & abs(ch1) + abs(ch2)” and press OK. The data will be transformed and all of your
views will be updated accordingly. If you save your data now, the transformed data will be saved.
If you want to undo this action, go to “Survey setup” and choose “Original Readings (intact)” as shown in the
dialog box below:
After you have pressed OK, all data will be restored and views updated.
Note: This option is not available if you load an ASCII .stn file. Hence, you will not
be able to reverse any changes once you save to a .stn file, and then load this file.
149
10. Using Multi-Telegram messages with MagMap2000
Some devices used during marine or airborne surveys can produce multi-telegram messages. Most
common example is the GPS with its variety of NMEA strings. Depending on GPS model and
configuration, it can produce GGA, VTG, GLL, etc messages.
These messages could not be used with MagMap2000 in the past. However newest software version
allows handling of such devices. Here we provide and example how GPS stream can be used to extract
VTG messages and include results in the MagMap200 export output file.
It’s important to understand that MagMap2000 uses prefix to select strings of interest. If beginning of data
string matches to the user-provided prefix, data is accepted (or rejected, depending on policy). Then data
is parsed the same way as for Generic Serial Device. Channel names can be assigned by default (in this
case program counts all available channels) or can be defined by the user. In latter case some of the data
channels could be omitted (for instance, only speed could be extracted from the VTG message and
meaningful channel name can be assigned).
Below this kind of data is referred as “prefixed data”. GPS VTG message is used below as an example,
but procedure works the same way for any kind of prefixed data.
To use this feature (sequence of operation is not important) either of following conditions should be met:
• At least one GPS or SURFER file should be loaded. Those files should refer to the same time
interval as prefixed data.
• MagLog Survey file should be loaded. This includes GPS, line number file, magnetometer (or
other devices). After all this information is loaded, prefixed messages could be loaded from the
same GPS file.
• Prefixed data should be acquired with MagLog program, or at least have time stamps at the end of
the message. Time should sequentially increase.
As a result, user is able to export speed values with its positions, pretty much the same as it is done for
magnetic field.
To start loading prefixed data, go to File / Open menu and select file type as Generic Prefixed Serial Data as it
is shown on the following screen shot:
150
Then select file you wish to load, for instance GPS file (it also contains VTG messages). Program will
show preview of the file (first 20 strings) and allow prefix definition as it shown here:
151
To define prefix, type “$GPVTG” and select Policy as Accept data WITH prefix. It means all other strings in
the file are ignored.
As defined by NMEA VTG messages has the following type of information: heading (degrees from true
or magnetic North) speed over the ground (knots), speed over the grounds (kilometers per hour). As an
NMEA message, it also has checksum. It means by parsing this message program will produce for
channels. Last channel (checksum) is meaningless. Be default channels have names CH1, CH2, CH3 etc.
Those default names can be re-assigned by pressing Assign channel names button. Assignment dialog is
displayed as the following:
152
Use this dialog to re-assign channel order and names. If any of the assigned channels are defined (right
pane of the dialog is not blank) they are used in place of defaults. In no channel is specifically assigned
(right pane of this dialog blank) then default names from the left pane are used. In this example we only
want to get speed in knots (channel #2) so select CH2 on the left and press Add button. Type desired
name as shown (no spaces please):
153
Press OK to complete assignment, and then press OK one more time to load actual file. Progress dialog is
shown on the screen:
File is loaded the same way as MagLog NT Generic serial device file. When load is completed, graph window
is shown on the screen (in addition to other data windows). Note that you should not close this window
till you complete data processing. You can perform all kinds of operations with these data. The window is
shown here:
154
Now, when all other data loaded and is ready for export, you can see speed in the export list:
If data is exported into separate files with default file names, prefix is used as part of this filename (except
non-alphabetical characters, such as “$” in this example).
This example (speed values) is mostly for QC purposes (for instance, one could check what was actual
speed during the survey). More practical example would be to use height measurements from GGA
155
message. Those are taken in already, but only if other data (such as magnetic field) is present. Let’s assume
the survey where GPS is the only device, and the goal is to plot an elevation map. Loading GPS data into
MagMap as usual, and then loading GPS again as prefixed data can achieve this.
Note: This may not work as expected, if GPS during the survey loses some of the
fields in its GGA message (for instance, due to poor satellite coverage). This will shift
channels in the message and some other channel can substitute height, which is not
desirable.
156
11. Other utilities
11.1. BintoAsc utility
G-858 files have a binary format and therefore it is impossible for the user to modify them. However,
sometimes it is necessary to access the data directly (this might be due to a hardware problem), or if the
user wants to do data manipulation outside of Magmap.
The BintoASC utility allows the user to convert the binary data file to an ASCII file and save it.
Note: If the user uses Magmap to save ASCII data, often data like GPS data or serial
data logged onto the console will not be saved in the same file.
To use this feature, go to the File menu and select Binary 858 ->ASCII. You should see the following
dialog:
Here you can select the input and output file names. Additionally, you can select if you want to output all
your data, or output specific types of readings.
When you press “OK”, the program will run a conversion routine and create your output file. You are
free to load this file into Magmap or use it for other purposes. Please see the G-858 manual for the
description of the formats of the above types of entries.
157
11.2. GPS Interpolator
11.2.1. Introduction
This utility is used to assign locations to large marine or airborne Maglog NT data files using
interpolation of GPS positions. Such large files are not well handled in MagMap and should
therefore be processed in Surfer or Geosoft type programs.
Features:
Magmap provides a wizard to simplify use of this utility. It can be accessed by selecting the
menu item “Interpolate GPS wizard” located in the File menu.
158
11.2.2. Select Input and Output Files
After you have selected “Interpolate GPS wizard”, you should see the first panel of the
wizard that allows selection of input files, input file format, output file, and output file
format.
• Data file: name of the MagLog NT data file to be processed. Most of the data files are
in ASCII (text) format. Data files from the GR-820 spectrometer are in binary format.
These files will contain GR-820 as part of the filename. Data file records must end
with a MagLog NT time stamp formatted as _MM/DD/YY_HH:MM:SS.sss.
159
• Gps_file: name of the MagLog NT location data file to be processed. This file will
usually contain ‘GPS’ as part of the filename. GPS file records must contain NMEA
$GPGGA data and end with a MagLog NT time stamp formatted as
_MM/DD/YY_HH:MM:SS.sss. If GPS records contain ‘*’, the GPS checksum will
be verified.
160
A sample acceptable GPS reading is shown below:
$GPGGA,114804.00,1838.5075,N,11228.2062,E,1,8,01.0,+00058,M,+00052.9,M
08/10/99 09:05:02.530
• Line_file: name of the MagLog NT line file to be processed. This file will
contain the line number as part of the filename.
This file is created by Maglog NT and should list the times that each line
starts and stops. A sample line file is shown below:
161
11.2.3. Interpolating Other Data Sets
It is possible to interpolate additional data sets to the magnetometer time using the MagMap2000
GPS Interpolator Wizard. This involves a simple linear interpolation based on time stamps in the
data files to be interpolated and a magnetometer data set with time stamps (time stamps means
that every reading has a time associated with it when it arrived at the computer). For example, one
could interpolate fluxgate magnetometer data acquired during an airborne survey to the main
cesium magnetometer times. This interpolated data set could then be used for Compensation
processing.
If there is no need to interpolate additional files, simply press “Next” to advance to the next
screen.
The middle of the dialog box is initially blank. Files are added by clicking on “Add data
file…” and removed by using the “Remove data file(s)…” button. Press “Add data file”
button to add new data file to the list. First a standard “Open file” dialog appears where
the name of an existing file will be selected. Then a second dialog box appears where
162
fields in the file will be specified.
The data file name is listed in the upper box and the beginning of the file is listed in “File
beginning” box so the user can make decisions about which parts of the messageto use
for interpolation.
Use the channel numbers dialog box to specify which “channels” (columns) are to be
interpolated. It is possible to specify more then one channel using a comma separated list.
Magmap2000 uses the same parsing rule as MagLogNT, i.e., any character that is not a
digit, dot or sign is interpreted as field separator.
Names should be assigned to the channels using the channel names enter field. Use a
comma-separated list to specify the names for more then one channel. No spaces are
allowed in the channel name designator box. These names will appear in the headers of
Surfer or GeoSoft interpolated output files.
The use of the prefix designator allows filtering of multi-message files. When the prefix is
left blank, all strings in the file are used for interpolation. A prefix causes the program to
accept only strings that start with the designated prefix.
The task is to have both fields (HEHDT and HEROT) interpolated to the magnetometer time
and appear in the interpolated output file as two columns.
Procedure: Press “Add data file” on the Wizard’s screen and select GYRO compass file. Enter
prefix “$HEHDT” in the prefix field, enter the channel number as 1 and type the channel name
as “Heading”. Then press “Add data file” again and select the same file. Now use “$HEROT” as
prefix, set the channel as 1 and the name as “Rotation”. With this setup the program will read the
input file twice: once taking only heading strings and a second time taking only rotation strings.
Both these variables are interpolated to the magnetometer time and appear in the resultant file as
two separate columns with the names “Heading” and “Rotation”.
The ‘Transfer Lat/Lon -> UTM on input’ check box must be checked to enable the other
items on this panel. If this box is not checked, GPS_Interpolator operates with unchanged
geographical coordinates which results in less accurate positions.
The ‘Transfer UTM -> Lat/Lon on output’ check box us used to generate accurate
interpolated locations in geographical coordinates.
164
Ellipsoid parameters section:
The ‘Ellipsoid name’ list box allows choice of 12 pre-defined ellipsoids or a user defined
ellipsoid. Pre-defined ellipsoids include WGS 84, GRS 80, WGS 72, Australian 1965,
Krasovsky 1940, International 1924, Hayford 1909, Clarke 1880, Clarke 1866, Airy
1830, Bessel 1841, and Everest 1830.
‘Major axis’ and ‘Flattening’ text edit boxes are enabled only for user defined ellipsoids.
• Central meridian (degrees) - text edit box is used for manual entry of the
Central meridian to be used for location conversion. This box is disabled
when ‘Choose central from data’ check box is checked.
• Scale factor - should usually be left at 0.9996. The text edit box is provided
for the rare case where a different scale factor can be useful.
For more discussion on individual UTM parameters, please refer to the UTM
discussion in the Magmap2000 manual.
165
The third panel is used to start the interpolation process using parameters set on
the first two panels.
166
12. Exporting to Geosoft or Surfer
Before you export, you should be aware of which type of survey you are dealing
with. This is important because different options are available for different types of
surveys and the export menus and procedures appear quite different.
Each of the following sections will define the type of surveys that are supported.
This includes
• Surveys done with GPS where you drew the GPS map using the option
“Draw new GPS map using the features of a regular 858 survey”.
To export your survey to Geosoft Analysis and Plotting Software or Surfer, you
must have an open survey file.
When several survey files are open, you must select the survey you want to export by
selecting the window (just click the mouse button anywhere in the window).
If you used a base station and want to do a diurnal correction to your data, you must
also open every base station file which covers the range of times desired.
Once you have an active survey, select “Export” from the File menu item or right
click the mouse button inside the grid and select “Export” from the pop up menu.
The following export dialog box (Figure 69) will appear on the screen.
167
Figure 69
In the case of mapped survey mode, you can redefine your grid. This will change the
coordinates of every position acquired in your survey. You should only use this
feature if you set up invalid grid information in the G-858 magnetometer before
doing your survey, such as invalid line or mark spacing.
In the “Sensor type” combo box choose the configuration of your sensor(s):
12.1.1.1.ONE SENSOR (1 OR LEFT OR TOP OR
FRONT)
You had one sensor only, and it was connected to the first connector.
12.1.1.2.ONE SENSOR (2 OR RIGHT OR BOTTOM OR
BACK)
You had one sensor only, and it was connected to the second connector.
12.1.1.3.HORIZONTAL GRADIENT (RIGHT - LEFT
WHEN WALKING ACCORDING TO
“FORWARD” DIRECTION)
You had two sensors (second sensor on the right of the first sensor).
168
You want the horizontal gradient with the left sensor magnetic field subtracted from
the right sensor magnetic field when you walked according to the picture defined in
“FORWARD DIRECTION”.
For example, if you set the “FORWARD DIRECTION” to the first item (line
going up in the increasing Y direction):
Every survey line going in that direction will have the gradient computed as the
following: right -left.
Every survey line going in the opposite direction will have the gradient computed as
the following: left -right.
12.1.1.4.HORIZONTAL GRADIENT (LEFT - RIGHT
WHEN WALKING ACCORDING TO
“FORWARD DIRECTION”)
You had two sensors (second sensor on the right of the first sensor)
You want the horizontal gradient with the right sensor magnetic field subtracted
from the left sensor magnetic field when you walked according to the picture
defined in “FORWARD DIRECTION”.
For example if you set the “FORWARD DIRECTION” to the first item (line going
up in the increasing Y direction):
Every survey line going in that direction will have the gradient computed as the
following: left - right.
Every survey line going in the opposite direction will have the gradient computed as
the following: right -left.
12.1.1.5.VERTICAL GRADIENT (TOP - BOTTOM)
You had two sensors, one on top of the other.
You want the vertical gradient with the bottom sensor magnetic field subtracted
from the top sensor magnetic field.
12.1.1.6.VERTICAL GRADIENT (BOTTOM - TOP)
You had two sensors, one on top of the other.
You want the vertical gradient with the top sensor magnetic field subtracted from
the bottom sensor magnetic field.
12.1.1.7.LONGITUDINAL GRADIENT (REAR -
FRONT WHEN WALKING ACCORDING TO
“FORWARD DIRECTION”)
You had two sensors,second sensor in front of the first sensor.
169
You want the longitudinal gradient with the front sensor magnetic field subtracted
from the rear sensor magnetic field when you walked according to the picture
defined in “FORWARD DIRECTION”.
For example if you set the “FORWARD DIRECTION” to the first item (line going
up in the increasing Y direction):
Every survey line going in that direction will have the gradient computed as the
following: rear - front.
Every survey line going in the opposite direction will have the gradient computed as
the following: front -rear.
12.1.1.8.LONGITUDINAL GRADIENT (FRONT -
REAR WHEN WALKING ACCORDING TO
“FORWARD DIRECTION”)
You had two sensors, second sensor in front of the first sensor.
You want the longitudinal gradient with the rear sensor magnetic field subtracted
from the front sensor magnetic field when you walked according to the picture
defined in “FORWARD DIRECTION”.
For example if you set the “FORWARD DIRECTION” to the first item (line going
up in the increasing Y direction):
Every survey line going in that direction will have the gradient computed as the
following: front -rear.
Every survey line going in the opposite direction will have the gradient computed as
the following: rear - front.
12.1.1.9.2 SENSORS TO COVER 2 LINES
You had two sensors, with the second sensor on the right of the first sensor.
You do not want the horizontal gradient, but instead assume that each sensor covers
one line (to get a better land surface coverage). Often used for archeology and UXO
surveys.
In the “bias value for magnetics” just type any positive or negative value you want to
add to the gradient.
170
12.1.4. File format
You can choose between Geosoft or Surfer output format.
• Sensor 1 reading
• Sensor 2 reading **
• Time, Date, Line, and Mark (all features in land-based surveys done
with G-858 magnetometer)
171
12.1.8. Base station
If no base station files are currently opened, a blinking warning message will remind
you that if you wish to export your data with diurnal correction you must first open
the base station file or files. If you have a proper base station file open, no message
will be displayed and your data will be diurnally corrected.
For example, when you start a line at X=0 and Y=0 and end a line at X=0 and
Y=100 and the magnetometer cycle time is set to one second, you will only get the
reading only a second later and the magnetometer will automatically record readings
every second. By definition when you press the end of line key, you stop recording,
and therefore you may or may not record the last reading at exactly X=0 and
Y=100.
• Airborne survey
• Land based survey where GPS was loaded using the option “Draw new map using GPS
and the features of a Maglog NT survey”.
172
Export menu for Maglog NT survey
This menu is quite different from the one that was discussed in the previous section. One of the
most noticeable changes is that you can output the additional positions specified in your GPS
Offset menu.
• Export data separately to default filenames – This will export your data to separate
files. You should get one file per device.
Be careful if you choose this option. This option will ignore any filename you type in!
• Select Output Path/Filename – choose any filename and path you wish to store your
data in.
• Select Fields to Export – Highlight all fields you wish to export. The highlighted
fields will be exported with a heading identifying each column. You can select more
than one field by holding down your CTRL key and selecting objects with your left
mouse button.
• Export MagPick profile – This will export a profile that you can use with the
MagPick program.
This outputs a file in Surfer format. You can select any combination of fields to export and this
will export them (with headings) in a format compatible with Surfer.
Here we chose five fields to output – two position pairs that we specified, and the first field
reading.
Note: Surfer will only read data files with the extension “.dat”, “.slk”, “.csv”
or “.txt”.
173
12.2.2.2.GEOSOFT
This outputs a file in Geosoft format. A Geosoft Analysis and Plotting Software file has the
extension “.xyz” and it has headings and a separate line that lists the line number. A typical
Geosoft file is shown below:
mag1_X mag1_Y mag2_X mag2_Y G-858_1
LINE 0
-122.1849910 37.4302225 -122.1849192 37.4301530 50199.07600
-122.1849918 37.4302233 -122.1849200 37.4301537 50096.44800
-122.1849926 37.4302240 -122.1849208 37.4301544 49944.51800
-122.1849933 37.4302247 -122.1849215 37.4301552 49810.18600
Here, the fields outputted and the format look very similar to the Surfer file. The only noticeable
difference is the separate line that lists the line number. This is updated later in the file.
12.2.2.3.TEXT (SPACE DELIMITED)
This outputs all the fields with one space between each reading.
12.2.2.4.SURFER – 3 COLUMN
This export procedure is particularly useful when you have two sensors with corrected positions
for each sensor and you would like to export into one Surfer compatible file. Normally, the
easiest thing to do would be to export a file two times – one file for the first set of readings and
positions, and the second for the second set of readings and positions. You would then need to
concatenate these two files. We tried to save some work by creating this file format. For each set
of positions that you select, you may select one reading. This assumes that the first positions you
specified will go with your first reading, and the second positions you specified will go with your
second reading, etc. If you don’t have a position set to go with a reading, or vice versa, you
will get an error message such as the one below:
In one case, we chose to fill out the output screen like the following:
174
Here is a portion of the output file we obtained:
X Y Data
-122.1849910 37.4302225 50199.07600
-122.1849192 37.4301530 0.00000
-122.1849918 37.4302233 50096.44800
-122.1849200 37.4301537 0.00000
Note: In this example, the first magnetometer had valid readings and the
second had a series of “0” readings. This is evident when we look at the output
file – one line corresponds to the first magnetometer and the second
corresponds to the second magnetometer. Three fields are always outputted.
175
13. Saving your editing
Magmap allows you to save your editing only in the G-858 ASCII format, so can use your favorite
editor or spreadsheet to view or edit your data.
The “Save” option is not available with any Maglog NT files. This applies to the following:
• Any G-858 files that stored GPS data that was then opened as a Maglog NT file. If
you want to view this file in its ASCII form, you need to use the program
bintoasc.exe as described below.
To save a file, select “File” from the menu then “Save” or “Save as”.
If you save a GPS map window, only GPS locations and magnetic readings will be saved.
Positions, field notes, and pause events will not be saved.
If you save a regular map window (using mark location), only position locations and magnetic
readings will be saved. GPS data, field notes, and pause events will not be saved.
You should never save a file with the same file name as your original file.
If you wish to process the data yourself or if you write a program that reads G-858 ASCII format
you can use the DOS utility program “bintoasc.exe” that is shipped with Magmap. It allows you
to decompress a G-858 binary compressed file to G-858 ASCII format.
176
14. Printing
Magmap allows you to print any Magmap window. There are also a few options available that
will help you print.
• You can change the title of a plot by clicking on the middle of the title on a window.
This is useful when you need to distinguish different plots.
• You can see a preview of your printing by choosing “Print Preview” from your file
menu. It will immediately show you a preview of how your hard copy should differ
from your screen layout. Namely, Magmap will set a white background and scale the
picture to occupy most of the available space. This will allow you to have a
reasonable layout, but it will prevent you from using exact scaling. Additionally, for
contour maps, the coordinate mesh will be plotted.
You can also choose the type of printer you want to use and the parameters by selecting “Print
setup” from the file menu. We recommend that you use a landscape paper orientation. The
printers available are set up using the standard Windows printer setup.
177
15. Example : Manipulating a Survey with a .csv file
The option “G-858 bin files and .csv file” currently supports files that are created from Trimble’s
Pathfinder software. The “.csv” file is a file that has differentially corrected positions that were
produced from the Pathfinder program. The 858 binary file should have GPGGA data strings
recorded as serial data and there should be a one to one correspondence between readings in the
.csv file and GPGGA strings in the binary file.
This section describes the operations you need to go through to load these files and some basic
data operations you can do with this.
Make sure that: 1) You select files of type “G-858 bin files & csv file”
2) You select both files to be opened. You can do this by holding the “CTRL”
button and left clicking on the applicable files.
After this, you should see your normal position screen without your GPS data. In this case, we
got the following which reflected the manual position points we took.
178
Opening screen after loading .csv files
We will be prompted to add some time corrections to the data. This allows us to correlate the
time reported from the GPS (UTC) with the time that the correction file is given in (Greenwich
time).
179
Our GPS data is shown below:
It is important to set these offsets correctly. If they are set wrong (e.g. we use 8 hours, -20
seconds instead of the above), we can get a dramatically different picture, such as the one shown
below:
180
Map plotted using incorrect GPS time corrections
From here, you can despike your data and do your normal GPS smoothing operations.
181
16. Sample Analysis from a Maglog File
This section discusses how to analyze data from Maglog NT files. It specifically covers operations
that you might need to do, such as calculating additional magnetometer position points,
smoothing GPS data, and removing spikes from your data.
You need to make sure that you use one of the two options available in the “open” dialog:
• Maglog projects (*.Survey) – this is used to load a survey and all associated files
182
After this is done, you should see your two files appear in Magmap as follows:
GPS file
183
16.3.1. Rescaling the data
We would first like to re-scale the data so that we can see the source and extent of the dropouts.
We can re-scale by clicking on the right mouse button. This will bring up the following pop-up
menu:
From here, we should select “Set Min and Max Field”. The following dialog box will then come
up:
By default, the minimum field is set to the minimum reading received in the survey. Since the
survey has dropouts, the minimum field will be set to 0. In this case, we need to change it to
49900 (or some other value that will allow us to see the data) so that we can see some features
besides the dropouts. If we press “OK”, our G-880 screen changes to reflect the scale change.
184
880 screen after scale change
We can plot this by right clicking on the mouse while the magnetometer window is active. We
will get the same pop-up menu that we saw before:
185
In this example, we choose the options “Analog_Device_1” and “Analog_Device_5” because we
know that we have two 880 magnetometers, each with signal, depth, and altitude measurements.
The data that is logged will appear as the following to Magmap:
Thus, if we wanted to graph the two signals instead, we would choose to display analog channels
2 and 6.
After we have pressed “OK”, the magnetometer screen should look like the following:
186
Scaled magnetometer data
The generic option: “Linearly Interpolate Drop-outs” will remove drop-outs with varying numbers of
points involved. It will replace the drop-outs points with valid points nearby. This doesn’t
involve much conscious effort on behalf of the user. If we select this option through the despike
menu, we will get the following:
187
In this case, the drop-outs have been removed. However, there are still spikes that need to be
removed. The readings that are left after removing drop-outs are not zero readings – but they still
are significantly different from normal readings to count as spikes.
We can remove spikes using the “Despike Setup” also available through the “Filtering” menu.
In this case, these options work well. If we used this to filter our drop-outs, we would need to
use more than one point for our “Maximum Spike Readings”. This is how the program determines
it is dealing with a spike. Since spikes that are due to error tend to be short in duration (a few
readings), if you have any more readings composing your spike than the “Maximum Spike
Readings”, Magmap will assume that this is not a spike and it will leave it alone.
Note: Magmap keeps a history of the changes you have done to your data. If
you don’t like the changes you made, you can press “Undo” to get back your
original data. In a few cases, (e.g. where you remove the drop-outs, or linearly
interpolate them) you will not be able to undo your changes after your first
despiking. Be careful and undo any changes you don’t want to keep before
despiking.
You should get the following:
188
Scaled magnetometer data with drop-outs removed
If you want, you can re-scale this again to fit the bounds of the data.
You can also remove all your drop-outs by despiking, rather than by selecting “Remove Drop-outs”.
Be careful doing this though because some drop-outs that last for a long time may not get
removed if they are longer than the maximum number of points involved allowed in a spike.
189
This section will show how to change from geographical to UTM coordinates, and how to
smooth the data.
190
In this dialog box, we select Enable UTM Transformation, which will cause the coordinates to
be output in UTM format. Additionally, we set the Northing and Easting parameters equal to
0. This will make it easier to add offsets later. We also allow the computer to choose the prime
meridian.
Now we would like to graph this in a local coordinate system – e.g. instead of using the huge
numbers that are given from absolute UTM coordinates, we would like to shift all points by a
given amount so that the final numbers are more manageable. To do this, we can set offsets.
Again, we go to the “UTM Setup” menu and fill in the offsets as follows:
191
Note: Here we select the values for Northing and Easting by taking the
negative of the two coordinates numbers in the bottom left corner of the GPS
graph. (See graph).
Note that the bottom left corner now corresponds to (0,0). If you export data, it will also be in
this coordinate system (with the offsets you selected).
192
16.4.2. Smoothing GPS Data
We would now like to smooth the data. We can do this by selecting “Smooth Position” from the
GPS menu.
Since this data was done as one line, we only see one set of parameters we can change (this may
be different for Maglog surveys taken using several lines!)
193
GPS data with smoothed path drawn
We accept this smoothing by pressing “Close and Apply”. The final map shows our smoothed
positions:
194
16.5. Other operations
16.5.1. Setting Offsets
In this survey (as in most surveys), the positions of your magnetometers will not be the same as
the GPS positions reported. This can be easily seen by looking at a contour map generated by
Magmap for the example data set.
This should appear to be a grid of 9 circular anomalies. However, due to the difference between
the GPS position and the magnetometer position, the actual map doesn’t come out that way.
We can correct this by setting offsets. Select the “GPS Offset” menu option (under GPS).
195
This will give you the following menu:
Recall that we have two magnetometers: one is dragged 40 meters behind the boat, and the
second is dragged 50 meters behind the boat. In the above dialog box, we have specified a
dragging distance between the two magnetometers (45 meters) and we have specified two
positions we would like to calculate with relation to this dragging distance (mag1 and mag2).
You have the opportunity to export these positions and use them in gridding programs such as
Surfer.
Note: These positions are calculated later. They are not available from within
Magmap.
16.5.2. Exporting to a file
We have now smoothed and despiked our data. We would like to export it so that we can make a
grid in a 3D plotting program such as Surfer.
196
To do this, select “Export”. You should see the following dialog box:
Additionally, we do not want to export data separately (this would keep the GPS and 880 data in
separate files), so we have unselected “Export Data Separately to Default Filenames”.
By default, all fields are highlighted and exported. If you don’t want to export all of the fields, you
can hold down the “Ctrl” key and select the fields you would like to export (as seen above).
Note: In the output format that we used in the example, we would have to look
at the data separately (one magnetometer at a time) or work with it further to
concatenate all the readings. You could instead choose to output the file in
“Surfer 3 Column” format which would create three columns and would
combine the readings for both magnetometers.
197
17. Analyzing an 858 Land Survey with GPS Positions
This section covers a sample analysis of an 858 land-based survey that has GPS input. In this
survey, we will look at the results obtained using the GPS readings as position points, and also at
the results obtained from doing a time interpolation between the endpoints of the lines.
This sample survey also has some common problems that might be encountered during a survey.
We will need to delete certain lines, and de-spike the data. There is also base station data.
• 0-70.bin -- This file has all measurements for the first 36 lines (up until X = 70).
• 72-100.bin – This file has the conclusion of the survey (until X = 100)
• There are drop-outs at the beginning due to a the magnetometer warming up.
198
17.3. Analysis of magnetometer files
17.3.1. Loading the file
In this section we will manipulate the two magnetometer files to get the correct mapped
information.
We begin by opening the file. The file type is “G-858 Binary File” as shown below:
Note: This menu will default to the last file type you selected. If you do not
select “G-858 Compress Files” as shown above, you will still get your data,
but you might not get your survey map.
199
Initial map for first G-858 file
If you wish to read more about editing the map, see the following chapter.
We picked this option because we will need to add a position offset to our data (our GPS and
magnetometer are not located in the same place). This option will allow us to do this.
200
GPS data
Magnetometer data
201
Additionally, we will want to add a position offset to account for the magnetometer being located
1.6 meters in front of the GPS. We would also like our final map to be in UTM coordinates.
17.3.4.1.DELETING A GPS LINE
In this file, we can safely ignore the line that was done twice because the position data is still okay.
For the GPS data, if a line was done twice, the GPS will report the same positions twice and this
should not affect any attempts to graph the data.
However, we still need to delete the last part of the line where the data could be bad due to the
loss of battery power.
To do this, we select the end of the line and then right click on the mouse as shown below:
In this case, there are two bad segments to delete. When you have finished deleting the bad
segments, your final map should look like the following:
Note: The sample file, 0-70.bin was improperly terminated and can cause
problems during export. If you have problems exporting this file, reload the
GPS file and don’t delete the lines above.
202
GPS map with segments deleted
17.3.4.2. CHANGING TO UTM COORDINATES
We intend to later import this data into Surfer so it would be helpful to have the coordinate
system in a more manageable format (UTM).
To do this, we need to first display the UTM coordinates with no offsets, and then change it into
a local grid that has some more manageable numbers.
We can change to UTM by selecting our GPS menu. Make sure that your GPS screen is the
active one – otherwise you will not be able to find this menu.
203
Note: Here we let Magmap automatically set the central meridian.
After we are done and press “OK”, we should get the following map:
204
As can be seen, this map is in UTM coordinates but they are quite large and difficult to deal with.
We would like to have this in a local coordinate system, where the bottom left-hand corner
corresponds to the coordinates (0,0). To do this, we need to add offsets.
We need to look at the GPS screen in UTM coordinates to get these offset values. In the above
case, we need to have a northing offset of –4142942.21 and an easting offset of –572012.08 (we
took the negative of the numbers located in the bottom left hand screen). – Remember that we
can find our offsets by using the formula:
Here, “Value” (for northing) was 4142942.21 and NewValue was the value we wanted to have –
namely 0.
After we have set these offsets, our map should appear as follows:
205
Magmap can correct for this by outputting an extra position point. (See Chapter 7 for more
discussion).
We would like to have two position points outputted – one labeled “GPS” and one labeled
“Magnetometer”. We can do this by selecting the “GPS” menu, and then selecting “GPS
Offset”.
From here, press “Multiple Offsets Setup”. Then fill out the dialog box as follows:
1) Make sure “Use Dragging” is disabled. This is mainly used in marine surveys where the
distance can be variable.
2) Delete the initial position point (you cannot edit and existing position point).
After you have finished, you are done with what you need to do with the GPS data for this file. It
is useful to save your work.
Note: The GPS position is included so that we can output our original
positions and compare the results later. It is not usually necessary to keep this
position.
206
17.3.5. Editing Magnetometer data that goes with the GPS
We also need to edit our magnetometer data that will be used with the GPS. As seen earlier, we
had data with drop-outs.
The first thing that we can do is to remove the drop-outs. We will choose to linearly interpolate
the drop-outs with good readings (see Chapter 6 for more of a discussion on filtering data).
To remove the drop-outs, make sure that the magnetometer screen is selected. Then, from your
“Filter” menu, select “Linearly Interpolate Drop Outs”.
Your new magnetometer screen should re-scale to reflect the filtered data.
207
Magnetometer data with drop outs linearly interpolated
There is also a downward spike in the middle of the screen that appears suspect. We can choose
to despike the data using a spike value of around 500 (spikes of 50 are common) but this won’t do
anything. If we look closer (e.g. run our cursor on the readings and see different values for
different times), we notice that this “spike” actually consists of several changing readings – in
other words, it is not a spike at all.
208
To open the file, you can select it as follows:
Note: You can open other files you saved earlier in the same way.
From here, you should see your base station data.
209
We can re-scale the data by pressing the right mouse button. We fill out the following:
Our new map shows us that the data is okay, except for the very beginning.
In this case, we don’t need to do anything more to the data because these spikes occurred a long
time before the survey began. Thus, we have done all we need to do with our base station file. (It
is always a good idea to save it though).
210
To export this data, we select the screen with the GPS data. This is important because certain
export options are only available through the GPS screen.
We can then go to the file menu and select “Export”. This should bring up the following menu:
• Time
• DIURNAL_3 – the diurnally corrected first sensor reading (Base Station Reading is
subtracted from sensor reading)
Note: We have not selected G-858_4 and DIRURNAL_4 because we only had
valid data on one sensor.
211
We also need to make sure that we pick a suitable name for our output file and uncheck “Export
Data Separately to Default File Names”.
When we are finished, we press “Export Now” and the data is exported. Some sample data from
the output file appears below:
magnetometer_X magnetometer_Y GPS_X GPS_Y QUAL_IND DOP G-858_3 DIURNAL_3
TIME
67.48 1.17 68.57 0.01 2.00 1.00 50199.07600 -469.64900 11:26:36.900
67.41 1.26 68.50 0.09 2.00 1.00 50096.44800 -367.06300 11:26:37.000
67.34 1.34 68.43 0.16 2.00 1.00 49944.51800 -215.13900 11:26:37.100
67.27 1.42 68.35 0.24 2.00 1.00 49810.18600 -80.81000 11:26:37.200
67.20 1.50 68.28 0.31 2.00 1.00 49719.18700 10.18600 11:26:37.300
67.13 1.58 68.21 0.39 2.00 1.00 49658.80000 70.56500
11:26:37.
As can be seen, we have exported all of the fields that we wanted into a Surfer
compatible format.
Note: Some versions of Magmap2000 have problems exporting this file. This
is due to the fact that the file 0-70.bin was improperly terminated. If you get a
message that prevents you from exporting such as: “One of the lines is shorter
than Look forward points + Look backward points. Export Aborted”, one
work around is to not delete any lines on the GPS map. It should work fine
then.
212
17.5. Further Data manipulation
17.5.1. Processing the second file
The second file also needs to be processed. This data has no noticeable spikes and the map is
correct, so the processing should be easy (refer to previous sections for operations you might
choose to do).
In this case, we loaded the file, drew a map using the GPS readings, and exported the data (no
special manipulation required).
Make sure you select the same export fields as you did in the first file. You should also
pick a new filename to export to. Make sure your base station file is open.
Note: Many of the settings you selected in the first file will be valid for the
second. These settings remain (e.g. GPS offsets, UTM settings) until you
change them.
If you have difficulty graphing the GPS data (e.g. the option is grayed out), make sure that you
have closed the first file. Magmap cannot have two GPS maps of the same type open at one
time. You should also check and make sure you have the map, rather than the base station file
selected.
We can choose to combine these files either in Surfer (refer to Surfer documentation), or do it
through DOS. To do it through DOS, you need to do the following steps:
1) Remove the heading from your second file. This is important because Surfer recognizes
headings only from the first line.
• You can do this easily from a text editor by deleting the first line.
2) Concatenate the two files with a command at the DOS prompt such as :
Copy a + b c
• Here, “a” and “b” are your first and second files and “c” is the name of the new
output file.
In this case, we named our two output files 0-70.dat and 72-100.dat. We wanted to
output to a file called 0-100.dat, so we typed:
213
These two steps are all that is necessary to concatenate a file for use in Surfer.
We find the shaded relief plot to provide the most useful view of our data. The contour plot
doesn’t work as well here because we have a large range of data, and some smaller anomalies will
not show up when graphed.
Here, we used
• X = Mag_X
• Y = Mag_Y
• Z = DIURNAL_3
We can also see several anomalies from looking at the map. There is a long pipe that runs across
the top of the screen, and another pipe that connects and runs parallel to the Y axis (X = 40).
Also note that these coordinates are graphed using the offsets we specified in the UTM screen.
We can compare this plot with one we obtain using the non shifted positions (replace mag_X and
214
140.00
120.00
100.00
80.00
60.00
40.00
20.00
As you can see, here we can see errors introduced because the magnetometer and the GPS are
not in the same place.
We can also look at the data that has not been diurnally corrected (but the positions have been
corrected)
215
18. G-858 Survey with no GPS data
Magnetometer surveys were typically done using a pre-defined, measured grid that would allow
the surveyor to establish a few position points, and time interpolate the rest. Recent
developments have allowed the use of GPS positioning but there are still many surveys done
without GPS positioning.
The sample survey discussed in Chapter 12 had GPS positioning but we also had a mapped
survey setup. We will use the same files but will approach the analysis from a mapped survey
perspective.
We begin by opening the file. The file type is “G-858 Binary File” as shown below:
Note: This menu will default to the last file type you selected. If you do not
select “G-858 Compress Files” as shown above, you will still get your data,
but you might not get your survey map.
216
Initial map for first G-858 file
2) The last line is bad because the magnetometer battery was drained. We need to delete
this line.
We will start by deleting the line at X = 52. To do this, we need to select the three position
segments, and then right click to select “Delete Selected Segment(s)” as shown below.
You can select more than one position point by holding down “CTRL”
Map with position points selected
217
Here, we have selected the three position points (we can pick these out by looking at the statistics
on the bottom of the screen. We search for X = 52, or Line = 26.
After we have deleted the line, the map should appear as follows:
We now need to shift all of the other segments over so that our other map coordinates will be
correct. To do this, we need to select all the segments and then press the right mouse button to
bring up the pop-up menu. From here we can select “Shift Selected Segment(s) by” as shown below:
218
Since we want to move all of the segments 2 meters to the left and we don’t want to change the Y
offsets at all, we should use an X offset of –2 as shown below.
In the new map, we see that the lines we selected have shifted over by 2 meters.
219
Map with segments shifted
Now we need to delete the last line where the battery went dead (X = 72). In order to do this, we
do the same delete operation that we did earlier. We select the two position points and then
delete the line.
220
18.1.3. Saving current work
At this point, it is advisable to save the current work in case there are changes to make in the
future.
To do this, press “Save”. Since we have not specified a save file before for this, we will be
prompted to type a file name.
In the above example, there are several “.stn” files already present. It will be helpful to pick a
naming convention so that we can easily identify the files associated with a given open file. In this
case, we chose to add “Map” to the initial file name. We can save this by pressing “OK”.
Note: The magnetometer data associated with the GPS and the magnetometer
data associated with the map are treated separately. Any changes (e.g.
despiking, removing drop-outs) that are done to one data set are not applied to
the other.
221
You can look at the G-858 data by positioning the cursor on a point that is not the endpoint of a
line (e.g. in this case, we positioned it on one of the black lines with arrows). Right click on it and
you should get the menu seen below:
After selecting “Plot Mag Field”, you should see the following magnetometer data.
Magnetometer data
The discussion of how to remove drop-outs and despike this data is given in section 0 above.
222
18.3. The base station data
The base station data is treated the same as in the previous chapter – section 0 above
To turn on the grid lines, go to “Scale” and select “Grid Lines”. You should get the following
menu:
Here we have chosen to manually specify our grid lines. We have chosen intervals that are
pertinent to our survey – e.g. the survey is 100 meters wide, so it makes sense to scale and try to
locate anomalies on the scale of 1-5 meters.
We have disabled Auto gridding because it is more difficult to find the values that will be assigned
to the intervals.
223
Map with grid
Note: Magmap will apply the grid options to the active window. Make sure
that you have selected the window you want to grid.
224
18.4.2. Using the Stack profile
The stack profile is particularly useful in a survey such as this, where the lines are straight and
parallel. We can plot a stack profile on our map by selecting an endpoint and right clicking to get
the menu:
This should bring up a menu that will allow us to select the sensor we want to plot – Sensor 1.
At this scale, we can see some anomalies but we would like to zoom in on them to characterize
them more accurately.
225
We can right click on our map and select “Plot Stack Profile on Top of Map”. From here, select
“Zoom In Stack Profile” as shown below:
Note: We can also change the height of the stack profile peaks by using the
“+” and “-“ keys.
After zooming a few times, we get the stack profile shown below:
226
18.4.2.1. CHARACTERIZING AN ANOMALY THROUGH THE
STACK PROFILE
There are a few things we can immediately tell from looking at the stack profile. One thing is that
we have a few objects that look linear.
We can tell this because there are several anomalies on the stack profile that have about the same
height and width and can be arranged in a line. This is characteristic of objects such as pipes.
These particular anomalies came from a wooden pipe wrapped in steel wire that was buried at the
survey site in the 1800s. The wood had long since rotted away but the steel wire surrounding it
was picked up by the magnetometer.
You can also look at other anomalies and try to get some basic characterizations such as depth, X
and Y coordinates, and width. You can even begin to guess what type of object it is. For more
discussion, read “Applications Manual for Portable Magnetometers” by S. Breiner. A thorough
discussion on interpretation of anomalies is given in Chapter 5.
227
This data can now be exported for further analysis. The technique is a little different than the
technique we saw in Chapter 17. For one thing, we have no GPS data to export.
We can export our data by selecting “Export” from the file menu.
In this case, we have chosen to export the position, sensor reading, and diurnally corrected sensor
position to a file named 0-70.dat.
Note: Even if there is GPS data present, it will be ignored in this case. You
need to select an active GPS window in order to export GPS data. GPS and
mapped data are always treated separately.
228
Make sure that you export this file with the same output parameters!
Refer to section 17.5.2 for information on how to concatenate the output files.
229
18.7. Analysis in Surfer
This section assumes you are familiar with basic operations in Surfer. If you are not, consult the
Surfer manual for discussions on the plots used.
We can easily make out the pipe features on this map. We can also see several other anomalies on
the map.
This map is useful because we can get a good view of where all the anomalies are. We can also
get a good idea of how deep they are by seeing how wide they are. A large peak (very dark spot
on the contour graph or a large peak on the stack plot) can be indicative of a large object, or
something near the surface. The width of the anomaly will distinguish this. An object buried
close to the surface will tend to have a high peak but narrow width, whereas the same object
buried deeper will give a larger width and shorter peak.
This is discussed in more detail in the Applications Manual for Portable Magnetometers by S.
Breiner.
You can also start to make a guess about the nature of the object by looking at the contour plot.
Naturally, objects such as pipes that lie horizontally are quite prominent on the graph. You can
make other general identifications by looking at whether the anomaly appears to be a magnetic
monopole or a magnetic dipole. For instance, in the contour plot below, we can see a number of
dipoles and a number of monopoles.
230
Contour map of data
These anomalies are dipoles that were produced by buried drums lying horizontally. In each case,
the number and orientation of the drums differed, as well as the depth. In the lower right hand
anomaly (3), 3 drums were spread out and buried shallowly (about 1.8 meters). In the lower left
hand dipole (2), two drums were buried at right angles at about 2.0 meters depth. One can see the
slight differences in the shape of the dipoles.
231
The drums are buried in the same orientation in the upper right hand dipole (1), but they are
deeper. You can see that the anomaly appears weaker and harder to see.
Additionally, since the magnetic field is induced, one can see that the dipoles point in a north-
south direction. The drums themselves point along the Y axis of the plot but the induced dipole
points in the direction of the magnetic field.
This one is from a drum buried vertically, the top of it about 2.0 meters from the depth. Other
monopoles are characterized by an object’s long axis being perpendicular to the direction of the
magnetic field. In this case, the magnetic field was equatorial (horizontal to the ground) and the
long axis was vertical. Note the converse relation – a dipole is seen when the long axis is parallel
to the magnetic field, as in the case of long sheets of metal being buried parallel to the ground, or
drums being buried parallel to the ground.
As you can see, we can start to make guesses about the orientation and depth of things from a
basic understanding of the contour plot.
There is a more in depth discussion on estimating depth given in the Applications Manual for
Portable Magnetometers -–Chapter 5.
232
This figure shows how you can find the location of the anomaly by finding the minimum and
maximum. Namely:
This is only a rule-of-thumb but it is useful when you have a large area and don’t know where to
start digging.
Keep these factors in mind when you are planning your survey.
233
19. Index table
Ellipsoid name, 136
Export, 31, 78, 169
A
ASCII file, 38 F
File
B
export, 78, 169
Base station saving, 177
plotting, 112 Filtering readings, 119
Base station magnetometer, 116 FIltering readings, 31
Baud rate, 34
Binary file, 39 G
BintoASC, 158
Geosoft, 31, 78, 169
GPS clock bias, 140
C
GPS Interpolator, 122, 159
contour plot, 88, 90, 91, 92, 94, 111, 131, 215, GPS Offset, 139, 141, 174, 196, 207
231, 233 Grid
Coordinate Transformation, 79 editing, 28, 74
csv, 36, 40, 175, 179, 180 Reversing, 76
Rotating, 77
settings, 78
D Shifting, 75
Group
DAT file, 78, 172 editing, 29, 67
data scale, 96 selecting, 29, 67
Decompress, 35
Deleting
line, 60 H
segment, 66
Despiking, 119 Host software, 23
destriping, 122
Deviation I
line, 57
segment, 63 Importing a survey, 32
Diurnal correction, 114, 168 Installation, 25
Diurnal parameters, 113
Downloading, 32
J
Downloading a survey, 27, 32
Dragging, 141, 207 Jump, 143
Dropouts, 121
L
E
Line
Easting, 138, 192, 193 deleting, 60
Editing deviation, 57
grid, 28, 74 editing, 27, 55
lines, 27, 55 moving, 59
positions, 27, 52 reversing, 56
segments, 28, 60 rotating, 59
editing positions, 27 shifting, 56
elevation, 96
234
M Rotation
segment, 65
MagMap
Data flow, 23 S
metal mapper, 149
Moving Saving to a file, 177
line, 59 Scale factor, 138
Multiple Offsets setup, 140 scaling, 49, 50, 91, 94, 101, 102, 103, 105, 117
Segment
N Definition, 60
deleting, 66
Northing, 138, 192, 193 deviation, 63
editing, 28, 60
Editing, 60
O Operations, 61
reversing, 62
offset, 27, 28, 29, 52, 54, 55, 56, 61, 62, 70, 75, rotating, 65
76, 89, 113, 138, 139, 140, 141, 148, 201, 203, Selecting, 60
206, 220 shifting, 62
Opening a file, 35 Selecting a group, 29, 67
Output file, 78, 172 Settings
grid, 78
P shaded relief plot, 94
Shifting
Parameter, 143 grid, 75
Plotting readings, 30 line, 56
Position segment, 62
deleting, 54 Sloping edges, 30
editing, 27 Smoothing, 142, 146, 194
moving, 52 Spl. Smooth, 143
selecting, 27, 52 splines, 147
shifting, 54 Stack profile, 96, 99, 226
prefixed data, 151 State Plane, 79
Prime Meridian, 136 sun azimuth, 96
printing, 178 Surfer, 31, 78, 169
psuedo gradient, 131 Survey
plotting readings, 78, 117
R
T
Readings
despiking, 119 True map coordinates, 50, 103, 117
filtering, 119
Readings, plotting
U
entire survey, 78, 117
parameters, 117 Undo, 30
Reversing Uploading, 32
grid, 76 UTM, 79, 134, 136, 137, 138, 148, 191, 192,
line, 56 203, 204, 205, 206, 214, 215
segment, 62
Rotating
grid, 77 X
line, 59
XYZ file, 78, 172
235