Abstract:
A method for generating, displaying and manipulating the orientation of three-dimensional vectors on a two-dimensional surface for visualization, statistical analysis and correlation. The method involves collecting vector data, transforming each vector into spherical co-ordinates, then assigning, a unique color related to a position on a spherical color model. A two-dimensional raster image filled with this unique color is centered at the point of measurement for each vector orientation as obtained from calculations or instrumentation capable of measuring this orientation. This methodology offers the user the ability to discriminate the location of specific orientations as well as the ability to define an enhanced full color gamut gradation for a specific range of orientation. This rendering of color-coded vector orientation enables easier understanding by the viewer of large data sets.

Description:
FIELD OF THE INVENTION  
       [0001]     The present invention relates, in general, to data processing and color-coded raster rendering of information, and more specifically, to a method for displaying the orientation of three-dimensional vector data on a two dimensional display or surface. This methodology has specific applications in the field of geophysical surveys.  
       BACKGROUND OF THE INVENTION  
       [0002]     The term “ferromagnetic” used hereinafter shall constitute the inclusion of antiferromagnetic and ferromagnetic materials.  
         [0003]     The term “Currie temperature” or “Currie Point” used hereinafter shall also encompass the Ne&#39;el temperature as would be used in the context of antiferromagnetic and ferromagnetic materials.  
         [0004]     In the study of vectors within vector fields it is often useful to correlate information that relates to the magnitude, orientation and spatial relationship of the individual vectors within said field.  
         [0005]     While numerous effective methodologies exist for displaying the magnitude of individual vectors, the ability to observe commonality of individual vector orientation within vector fields containing thousands of individual vectors on a surface or within a volume can be difficult to discern.  
         [0006]     The applicant is aware of several attempts in prior art, which relate particularly to the visualisation of three-dimensional vector fields. Reference may be had by turning to the following publications hereinafter further referred to in  superscript ; 
        Reference 1 : “Visualizing the Magnitude and Direction of Three-Dimensional Vector Fields” by Timothy J. Peter in AP-S International Symposium Digest, IEEE Antennas and Propagation Society, Volume 3, 1993, pages 1118 to 1121.     Reference 2 : “Visualization of Three-Dimensional Vector Fields and Functions” by J. Simkin, IEEE on Visual 3D Fields Digest, 1995/205, pages 1/1-2; and     Reference 3 : “Vector Field Animation with Maps” by Yamrom, IEEE Computer Graphics and Applications, Volume 15, Part 2, pages 22 to 24, March, 1995;     Reference 4 : “Combined 3D Interpretation of Airborne, Surface and Borehole Vector Magnetics” of Mueller, E. L., Morris, W. A., Killeen, P. G., and Balch, S. at the McConnell Nickel Deposit, Proceedings of Exploration 97: Fourth Decennial International Conference on Mineral Exploration, edited by A. G. Gubins, 1997, p. 657-666.     References 5 : “Mapping Remnant Magnetism Using the Local Phase” of Thruston, J. B., in Geophysics, Volume 66, No. 4, July-August, 2001, p. 1082-1089.        
 
         [0012]     Current methodologies for indicating vector orientation consist of streamlines, hedgehogs, stream-ribbons, stream-polygons and textures maps 3 . In all these methods there are difficulties in observing the vector orientation orthogonal to the observation surface or in selection of unique desired orientations.  FIG. 5  of the drawing sheets illustrates a vector grid that shows this effect. A novel texture mapping procedure is hereby disclosed that overcomes many of the shortfalls of the earlier methods.  
         [0013]     Many applications exist for the use of an effective display of vector orientation including studies in fluidics, electromagnetic fields, finite element analysis, magneto hydrodynamics, plasma physics, astrophysics magnetic and digital elevation models (DEMs). We will present an embodiment in which the application of the present method is illustrated in the field of geophysical geomagnetic mineral exploration. It is to be recognized that this embodiment is but only one example of the usefulness of this Vector Orientation Visualization Method.  
         [0014]     Geomagnetic surveys have many uses and are of particular interest in mineral exploration, as they may identify areas likely to contain some types of valuable deposits. They can also identify various lithologies, and provide some types of structural information.  
         [0015]     A common method of performing such a survey is to measure the magnetic field at various points along the Earth&#39;s surface while recording the magnitude variations in the magnetic field strength. Although dip needle surveys have been used in the past, until recently, the orientation of the magnetic field was not generally measured. The magnitude of variations in the readings was primarily due to variations in the magnetic susceptibility of the constituent minerals within the rock, and to a much lesser degree, the variation of the orientation of the remnant field. Such surveys, while giving useful data, which related primarily to lithological effects, disregarded the information of the orientation of the magnetic field.  
         [0016]     Recently, due to the introduction of tri-axial gradient magnetometer surveys, it has become possible to retrieve vector oriented gradient data. In addition, tri-axial fluxgate total magnetometer surveys can also be used to measure the orientation of the total magnetic field vector. Tri-axial fluxgate magnetometers have also been used in along a borehole to effectively delineate the partial outline of a dipole about a ferromagnetic ore body 4 .  
         [0017]     The ability to effectively display the orientation of the magnetic vectors over a surface or within a block of rock as by borehole surveys can greatly add to the geologic understanding of the survey area and can lead to a better understanding of the location and extent of several varieties of ore deposits.  
         [0018]     To understand the benefit of knowing the orientation of the vectors within the survey area it is important to understand the means by which variations in the orientation of the vectors can occur. Afier a high temperature thermal event, certain minerals in the earth will heat enough beyond their Currie Point Ne&#39;el temperature to lose their magnetism. When the temperature lowers, these minerals regain their magnetism and will now have a thermal remnant magnetic field aligned with that of the earth at the time of going through their Currie temperature. Since the earth&#39;s magnetic field is known to change directions and flip at various times, the differently oriented remnant magnetic vectors of these remagnetized minerals can indicate spatially segregated intrusives of contemporaneous age, multiple pulsed intrusives within a similar rock type, the extent of their thermal aureoles and the presence of folds and faults. Similarly alteration due to chemical—known as chemical remnant magnetism, CRM—and pressure changes—known as piezoremnant magnetism—can result in the formation of new minerals, which will retain a remnant orientation of the earth&#39;s field orientation at that time—known as detrital remnant magnetism (DRM) caused by the quiescent sedimentation of fine magnetic material aligning with the earth&#39;s magnetic field at the time of sedimentation.—DRM is another remnance whose variation can be used to recognize later faults and folds.  
         [0019]     For the purpose of this embodiment a magnetic survey can be performed using either a tri-axial fluxgate magnetometer or tri-axial gradient magnetometer. The vertical gradient can also be calculated if both horizontal gradients and the total magnetic field are known. These types of surveys produce the magnitude of the magnetic field for each of the three-orthogonal directions. The data is collected for process and analysis over the survey area by establishing a two-dimensional grid over the area of interest. Measurements of the magnetic field are taken at specific points on the grid either by aircraft or by an operator on the ground working along a cut or demarked grid. Each collection point is identified as a point on the two-dimensional grid, and will have an associated three-dimensional vector representing the magnetic field at that point. Alternatively, if the data should be collected along a borehole as a function of its distance from the collar, the location of the three-dimensional vector can be represented from the surveyed borehole information as a point in 3D space with x, y, z co-ordinates. The x,y co-ordinates can match the two dimensional grid of the surveyed area of interest mentioned earlier.  
         [0020]     Peters Reference 1 : This reference describes a method of representing an electromagnetic vector field defined over a three-dimensional surface. A triangular grid represents the surface, and the method uses color shaded magnitude contours and directed vector lines to characterize the vector field. Referring to the Figures, it may be seen that in this reference, the direction of a vector field is represented using arrows, and the magnitude of the vector field is displayed using colors. Such visualization is suited to vector fields for which the direction changes gradually and continuously. The magnitude of the vector field can then be visually identified by the intensity of the colours used.  
         [0021]     Simkin Reference 2 : This reference provides an overview of three-dimensional visualization, and only contemplates displaying vectors using arrows. The direction of the vector is shown using arrows while the length of the arrow visualizes the magnitude, In other words, the longer the arrow, and the higher the magnitude of the vector.  
         [0022]     Yamrom Reference 3 : This reference describes a method of animating vector fields and discusses the subject matter of displaying a 3D-vector field by changing the texture of the viewed surface with respect to time.  
       SUMMARY OF THE INVENTION  
       [0023]     According to an aspect of the present invention, there is provided a method of displaying three-dimensional vector field information on a two-dimensional display. After collecting the information from an area of interest, the Cartesian vectors are transformed into spherical co-ordinates. A colour model such as Munsell, HSV, HLS, or YIQ, is modified into a spherical colour model. In the embodiment of the modified HSV model, the spherical model is developed such that Hue is distributed evenly along the equator while the Value varies from white at the North Pole to black at the South Pole. Each of the vectors are plotted emanating from the origin of the coloured sphere using their spherical co-ordinates. An HSV colour is assigned to the area of the vector on the map based on the colour occurring at the corresponding piercing point of the vector with the designated HSV spherical colour model surface. This HSV colour is then used to fill the area or volume of influence of the vector. In the case of a monitor display the HSV values are converted to their corresponding RGB values. Ability exists using this model to define a particular orientation range of azimuth and dip such that the whole colour gamut may be displayed within a narrow range and hence further delineate variations in orientation.  
         [0024]     The utility of this methodology includes but is not limited to geophysical surveys.  
         [0025]     For purposes of clarification and to provide a fill disclosure, it was deemed absolute that color drawing be included in the present application in order to fully understand the scope of inventiveness.  
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0026]     These and other advantages of the invention will become apparent upon reading the following detailed description and upon referring to the drawings in which:  
         [0027]      FIG. 1  is a flowchart outlining an embodiment of a method of data collection and display of three-dimensional information using a two-dimensional display.  
         [0028]      FIG. 1   a  is a flowchart outlining an embodiment of a method of data collection and display of three-dimensional information using a two-dimensional display and indicating the variables used to describe the data and its manipulation.  
         [0029]      FIG. 2  is a schematic diagram of a grid of irregular outline established to define the collection sites of the three-dimensional information.  
         [0030]      FIG. 3  is an example of a standard base station calibration.  
         [0031]      FIG. 4  is an illustration of the definition of a spherical co-ordinate system based on Mathematical Conventions.  
         [0032]      FIG. 4   a  is an illustration of the definition of a spherical co-ordinate system based on Geological Conventions.  
         [0033]      FIG. 5  is a schematic diagram of a two-dimensional grid surface defined for locating the collection sites of the three-dimensional information (3D vectors).  
         [0034]      FIG. 6  illustrates the notation used in the orientation of the axis and coordinate systems.  
         [0035]      FIG. 7  shows an HSV colour coded sphere sectored in planes of equal area, each denoting a specific HSV colour with the 3D vectors being oriented so as to emanate from the origin to intersect a uniquely coloured surface of the sphere. It should be noted that this sectored sphere is presented in this illustration for the basis of clarity of methodology. In reality a continuum of millions of colours will be more frequently used.  
         [0036]      FIG. 8  is a table illustrates azimuth and dip angles of the 3D vectors as per the grid in  FIG. 2 .  
         [0037]      FIG. 9  is a table illustrating assigned HSV colour values to the 3D vectors as per the grid in  FIG. 2 .  
         [0038]      FIG. 10  illustrates a further endowment of this methodology that allows for the ability to illustrate a scattergram distribution of vector orientation on a projection.  
         [0039]      FIG. 10   a  illustrates a section taken from  FIG. 10  further expanding the details of the orientations of the vectors.  
         [0040]      FIGS. 11, 11   a  and  11   b  illustrate an endowment of this methodology that shows how a range of vector orientations may be sectored out and given higher visual discrimination to further emphasis subtle changes in vector orientation.  
         [0041]      FIG. 12  illustrates the described methodology of colour coded vector orientation as being applied to voxels located at sampled points along the 3D survey trajectory of a curvilinear borehole.  
         [0042]      FIG. 13  HSV colour assignment illustrates the relationship of the observed geological spherical co-ordinates and the development of the modified HSV spherical colour model. 
     
    
       [0043]     While the invention will be described in conjunction with illustrated embodiments, it will be understood that it is not intended to limit the invention to such embodiments. On the contrary, it is intended to cover all alternatives, modifications and equivalents as may be included within the spirit and scope of the invention as defined by the appended claims.  
       DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT  
       [0044]     In the following description, similar features in the drawings have been given similar reference numerals.  
         [0045]     Turning to  FIG. 1 , a flowchart outlining an embodiment of a method of displaying three-dimensional information, such as three-dimensional vectors, using a two-dimensional display, is shown. More specifically, the present embodiment is directed at the collection of three-dimensional magnetic field data  14  in an area of interest  10  on the surface of the Earth from a plurality of mobile stations to extract the residual magnetic field orientation of the underlying rocks. Said data is then calibrated  16  at a base station  12 , converted to spherical co-ordinates  18 , assigned colour values  20  then plotted  22  to a two-dimensional surface. A further detail of this flowchart and the variable names used to describe this data is shown in  FIG. 1   a.  Although the present embodiment is directed at a method of displaying three-dimensional magnetic field vectors for use in the study of geophysics, thermodynamics, fluidics, plasma physics, magneto hydrodynamics and the like, the method may also be used to display vectors representing any three-dimensional information in a two-dimensional display.  
         [0046]     Initially, in the case of a surface ground based geomagnetic survey the area of interest is selected and defined by setting up a grid over the area form which the magnetic field data is to be collected  10 , in the case of an aircraft survey, location of the flight lines would be obtained by utilizing rapid precision Global Positioning System (GPS) equipment.  FIG. 2  illustrates an example of a grid  100  including a carefully laid out baseline  102  which is located centrally within the grid and more or less parallel to the feature desiring investigation and dividing the grid area into two parts  104  and  106 . This baseline  102  is cut free of brush and trees so as to allow for a line of sight to measure out the locations of the lines  112 . Due to the trend of the baseline in  FIG. 2  being in a more or less South-North direction, each side of the baseline is segregated in an Eastern and Western section respectively. Lines  112  perpendicular to the baseline  102  further define the grid  100  and are demarked from the baseline via surveying or by the use of a turning board. Said lines  112  are then blazed or line cut and picketed at regular station intervals for later measurement with the instrument. In the exemplary illustration presented, the lines  112  are shown 100-meters apart with measuring stations  116  at 25-meter intervals along the line emanating from the baseline  102 . For computational purposes, an outline of rectangular area  114  is defined around each of the measuring stations  116  as follows: In the direction of the lines  112 , the opposite sides of the rectangle are located at the midpoint of each of the consecutive measuring stations  116 , which in the example case would be 12.5 meters. In our example, 12.5 meters would be used to define the on-line rectangle dimension of the measuring stations at the end of each of the lines  112 . Orthogonal to this, the distance between sides of the rectangle can be arbitrarily set between the previously determined online side-dimension as a minimum and the distance between lines  112  as a maximum. Various algorithms are known in the art for calculation of values between lines  112 .  
         [0047]     At each of the measuring stations  116 , the magnetic field data, the time of the measurement and the station location is recorded within the grid  100 .  
         [0048]     In order to distinguish measurements from the various measuring stations  116 , the measuring stations are denoted as LN(orientation)+(distance)(section). This is a standard notation in the exploration industry in North America and is further explained below. The ‘LN’ values represents the measuring station&#39;s line location with respect to a selected origin of the baseline  102  with L0+00 representing the measuring station at the origin on L00+00. In our example, since the lines  112  are orthogonal to the baseline and 100-meters apart, n=BLdistance/100, ‘distance’ being the distance along the baseline  102  from the selected origin. As one progresses in a northward ‘orientation’, the lines  112  are sequentially labelled L00, L01N, L02N . . . LN with (N+1) being the number of lines  112 . The ‘distance’ value represents the distance the measuring station  116  is located away from the baseline  102  and the ‘section’ value represents the section  104  or  106  in which the measuring station is located. Therefore, for a measuring station in the Eastern section  104  located 50-meters northward away from the baseline  102  on the perpendicular line corresponding to L04N, this measuring station may be denoted as L04N+50E. In the grid of  FIG. 2 , the number of baseline measuring stations  116  is indicated by the tick marks along the lines  112 . Further reference to  FIG. 2  illustrates the designation of the measuring stations  116  along the baseline. In respect to those measuring stations  116  occurring along the baseline the notation BL(orientation)+(BLdistance)(direction) is used.  
         [0049]     After the grid  100  has been defined, a magnetically quiet area, close to the grid, is selected so that a tri-axial magnetometer base station  118  may be located (step  12 ) and appropriately positioned to collect signals which are used to assist in the calibration of the signals collected at the measuring stations  116  The base station  118  instrument&#39;s clock is set to match the clock of the portable mobile tri-axial instrument. This will be described in more detail below.  
         [0050]     Prior to collecting the data (step  14 ) the operator initially records, from a distance, x,y,z readings in a magnetically quiet area with the tri-axial magnetometer, this could be a tri-axial fluxgate instrument suspended at the proper height and orientation by a remote nonmagnetic device. A second set of readings is then taken in the same location with the operator holding the instrument in the normal manner to take a reading. These two readings are recorded so that any magnetic effect from the presence of the operator can be removed. This procedure may be undertaken to confirm a consistent operator null correction  
         [0051]     Having undertaken the above operator nulling procedure, the operator walks with the instrument along the baseline and lines such that the individual bypasses each measuring station at least once. At each measuring station, the individual records and stores magnetic field data in the form of magnetic field vector information along with the station location and the time the measurement is recorded. Care is taken in the positioning and orientation of the instrument, aided by a sighting and levelling device to maintain a consistent x,y,z orientation of the probes. The magnetic field vector information collected includes information corresponding to the magnitude and direction of the magnetic field at the measuring station location. This may be achieved by simply writing down the magnetic field data on a notepad or by entering such information into a data collector for later dumping on a computer. The recorded signals X 1 ,Y 1  and Z 1  may be represented as X m , Y m  and Z m  with “ m ” representing the measuring station designation.  
         [0052]     At the end of each days work after the individual has collected the magnetic field data from each of the measuring stations visited, the individual retrieves calibration information from the base station corresponding to the time period over which the magnetic field data was being recorded. After retrieving the calibration information from the base station, the calibration information and the recorded signals are stored in a computer and the recorded signals are calibrated. This is a standard procedure known to those skilled in the art.  
         [0053]     The calibration process allows for the magnetic field data retrieved from the measuring stations to be a truer indication of their actual magnitude and direction. It will be understood that the magnetic field data may be affected by external conditions such as solar storms or diurnal variations.  
         [0054]     Turning to  FIG. 3 , a sample display of calibration information is shown. Base station averages, X bavg , Y bavg , and Z bavg , are calculated by averaging the magnitude of the stable magnetic field vector, in the X, Y and Z directions at the base station over the time period of the survey. Since the base station is assumed to be located in a magnetically quiet location away from vehicular traffic, power lines, transmitter sites, etc. it may be assumed that any recorded deviations from the base station stable average indicates that any measurements taken at that instance in time have been affected by the same external condition. For instance, at time=6:15, as indicated on  FIG. 3 , the X value measured at the base station is seen as a spurious X b  positive peak above the pre-calculated average, X b  and therefore a base station calibration correction value X basestation  equals the absolute value between X b  and the pre-calculated average, X bavg .  
         [0055]     After selecting a first-recorded signal X m , Y m , Z m  from the mobile station, the time at which the signal was recorded is retrieved. The calibration values for this time is then calculated as explained above. If the reading was taken at 6:15, the base-station calibration values X basestation6:15 , Y basestation6:15 , and Z basestation6:15  are retrieved and subtracted from the measuring station values X m , Y m  and Z m  corresponding to the first recorded signal in order to perform a first step of calibrating the measuring station value. After the first reading is calibrated, the magnetic field data from a second reading is selected and calibrated by retrieving the time measurement from the magnetic field data and calculating and retrieving the calibration values corresponding to that time measurement. This step is repeated until all of the magnetic field data from each of the measuring stations is calibrated. Therefore X calibrated  equals X m −X basestation(time)  which is equal to X m −|X bavg −X b(time) |. Values for Y calibrated  and Z calibrated  are calculated in a similar manner.  
         [0056]     After the values X calibrated , Y calibrated  and Z calibrated  for each measuring station have been calculated. A second step may be performed to remove the geomagnetic field vector from these values in order to obtain a residual magnetic value from the magnetic field data.  
         [0057]     This is achieved by calculating the Earth&#39;s magnetic field using the standard International Geomagnetic Reference Field (IGRF) procedures to determine the theoretical values with respect to the area of interest at the time of the survey These theoretical values obtained may be represented as X igrf , Y igrf  and Z igrf .  
         [0058]     The residual magnetic data values X residual , Y residual  and Z residual  may then be calculated from the equations X calibrated −X igrf , Y calibrated −Y igrf , and Z calibrated −Z igrf , respectively.  
         [0059]     After calculating the residual magnetic values, further methods known to the trade 5  may be performed which causes the induced magnetic field value to be removed from the residual magnetic data values to obtain remnant field values at the measuring stations.  
         [0060]     In order to calculate the remnant values X remnant , Y remnant  and Z remnant , the equations X residual −X induced , Y residual −Y induced , and Z residual −Z induced , are respectively determined. The remnant value represent the magnetic field data values were frequently acquired at the time the ferromagnetic rock forming minerals were formed or went through a Currie temperature episode of its ferromagnetic minerals. This remnant vector data can allow for the discrimination as outlined earlier in the claims  79 ( c ). After calculating either the residual or remnant data these x, y, z Cartesian values are transformed to mathematical spherical co-ordinates using the following equations: 
        r math =squareroot(X 2 +Y 2 +Z 2 ) where r is the radius of the unit sphere used to determine spherical co-ordinates;     Θ(theta) math =arctan (Y/X) where theta represents the azimuth angle, which is measured from the positive x-axis toward the positive y-axis theta and is a value between 0° and 360°;     Φ(phi) math =arctan ((squareroot(X 2 +Y 2 )/Z) where phi is the angle measured from the vertical positive z axis counterclockwise toward the x,y plane and is a value between 0° to 180°.        
 
         [0064]     These mathematical spherical co-ordinates are illustrated clearly in  FIG. 4 .  
         [0065]     The mathematical spherical co-ordinates of the magnetic vector at each measuring station may then be represented by (r math , Θ math , Φ math ).  
         [0066]     A conversion of the mathematical spherical co-ordinates (r math , Θ math , Φ math ) is then translated into the geological co-ordinates (r g , Θ g , Φ g ) by the following equation: (r g , Θ g , Φ g )=[(r math ), (Θ math −270), (90−Φ math )] as clearly illustrated in  FIG. 4   a.  These values are tabled in  FIG. 8 .  
         [0067]     After calculating the geological spherical co-ordinates for each measuring station, a unique colour based legend may be determined by assigning colour or patterns to the spherical co-ordinates. Several colour models or notations exist which have direct symmetry with respect to spherical co-ordinates such as red, green, blue (RGB); Commission Internationale de l&#39;Eclairage (CIE), hue, saturation, value (HSV); hue, saturation, luminance (HSL), CIE XYZ, YIQ, Munsell, TekHVC, CIE LUV etc . . . The last three of these models allow for a more continuous perception of colour to the human eye. Although any of these colour models may be used in the preparation of a suitable colour model via the proper mathematical translation. In this embodiment a modified spherical HSV colour model is used, where hue variation is assigned the azimuth angle. Axial symmetry about dipole axis may also be illustrated by having hue vary as dip.  
         [0068]     One example of assigning spherical co-ordinates to an HSV model may be realized by the following two equations, the use of each being dependent on the value of phi:  
         [0069]     When Φ&gt;0[H, S, V]=[(Θ g ), (2Φ g /π), (1)], when Φ&lt;0[H, S, V]=[(Θ g) , (1), (2-2Φ g /π)],  FIG. 6  illustrates this colour assignment.  
         [0070]     The above assignment relates a unique HSV colour value to the direction of a three dimensional vector corresponding to measured and calculated values from the measuring stations which can be applied to the area of interest.  
         [0071]      FIG. 9  graphically illustrates this principle. The particular range of values for the HSV colour values defines the limits of the above mathematical expressions. However, other parametric vector assignments can be used for different colour models or for the creation of a new unique colour model or even textural patterns.  
         [0072]     After each of the spherical co-ordinate vectors have been associated with an HSV value, see the table of  FIG. 10 , the colour representing the HSV value is then applied as a fill colour to a rectangular area centred at each of the corresponding measuring station locations along the lines with respect to the overall grid.  FIG. 5  shows the vectors as arrows occurring at the centre of a rectangular cell around the location of their measurement point illustrating the difficulty in determining change in magnitude, verses a change in foreshortened length due to a change in dip.  
         [0073]      FIGS. 8 and 9  provide two-dimensional displays of the magnetic field data in accordance with an embodiment of the present method. For enhanced clarity, in  FIGS. 5, 8  and  9 , the lines of the grid are not shown but lie along the centre of each grid cell.  
         [0074]     These coloured areas located about each of the measuring stations may also be input into any of the many standard contour algorithms that exist for plotting of these coloured pixelated areas to produce a colour map with gradational interstitial colours applied between the lines and covering the whole area. These interstitial colours could represent, along observed trends, the shortest distance between the two colours as illustrated on the spherical colour model.  
         [0075]     With the advent of computer processing capabilities, the unique display of this vector orientation data lends itself to a number of useful statistical presentations. In  FIG. 10  we see that an area of the colour coded data grid area can be selected ( FIG. 10   a ) and the orientation of all vectors within this area can be plotted on an orthogonal grid as a function of dip angle versus Azimuth ( FIG. 10   b ). This can allow for the fine cognitive discrimination of very subtle orientation regimes.  
         [0076]      FIG. 11  shows another attribute of this vector orientation presentation, which allows for the enhanced discretion of specific selected orientation. In this scenario a specific range of dip and azimuth angles are selected. This selection is illustrated in the highlighted quadrant of the sphere ( FIG. 11   a ) or the rectangular area expressed on the orthogonal dip verses azimuth plot ( FIG. 11   b ). A new colour chart is then created that contains the full complement of colours as the original model ( FIG. 11   c ). The appropriate cells that contain the dip and azimuth values within this range are displayed on the plan view of the grid area utilizing the full gamut of colours available.  
         [0077]     Additionally the magnitude of the vectors may also be shown on the coloured illustration either as dashed black and white line contours or as a Digital Elevation Model.  
         [0078]     Alternatively to the ground acquisition method described earlier, an aircraft flying over the grid may retrieve the above information. After the grid has been determined, the aircraft flies over a prescribed calculated grid lines to record the magnetic vector information at regular time intervals corresponding to a specified distance given the airspeed flown. In order to ensure that the aircraft follows the calculated grid lines, a GPS is be used to guide the aircraft. In this aerial data collection no baseline or ground demarcation of lines is required.  
         [0079]     The magnetic field vector information is collected via instruments on the aircraft as they fly over the measuring station and stored on record keeping equipment, such as a data logger, located on the aircraft. Calibration and corrections to the data are carried out later after each day&#39;s flying.  
         [0080]     An example of the 3D versatility of this colour coded vector orientation methodology is also given in  FIG. 12  where the method illustrates vector orientation sample points collected along a borehole within a given volume by the use of coloured voxels.