Versions Compared

Version Old Version 1 New Version Current
Changes made by

ian macleod (Deactivated)

Tom Popowski

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.
Tip

Here are some tips to get you started. You can edit this page to see how it works!

1. Create a page

  • Click "Create" and select "Blank Page" to create your first page.
  • New pages are created as children of the page you are currently viewing.

2. Add to your page

  • Click "Edit" to enter the Confluence editor and use the page layouts feature to structure your content using sections and columns.
  • Use headings to format your text and drag and drop images into your page to provide visual interest.
  • Click "Insert" and select "Other Macros" to add macros for navigation, special formatting and other media.
Panel

On this page:

Table of Contents

3. Organise your pages

Here are some tips for organising your content.

  • Change the page order

The sidebar on the left displays your pages in a hierarchy. If you have Space Administrator permissions you can click "Space Tools" > "Reorder Pages" to move pages around.

  • Add labels

Labels help keep pages organised and make it easier for you to find the information you need. Click "Labels" at the bottom of a page to add or edit. The "Related pages" section on this page uses labels too!

  • Make templates

Standardise and speed up the page creation process with templates. You can create and format a template with page layouts, standard headings and instructional text for hints and guidelines. Check out our sample page on "Making a template"

Related pages

Filter by label (Content by label)showLabelsfalsespacesGGRsorttitleshowSpacefalsetypepagelabelsdocumentation-space-sample


Page Properties


SummaryA Geosoft Python script that classifies multi-variate data using Self-Organizing Maps (SOM, an AI technique).
Researchers

ian macleod (Deactivated)
Rob Ellis (Deactivated)
telma.aisengart@geosoft.com (Deactivated)

Githubhttps://github.com/GeosoftInc/Geoscience/tree/v9.4.0/Open%20research/self_organizing_maps
Publications

Quantitative Magnetization Inversion

Self-Organizing Maps applied to Magnetization Vector Inversion

Application of SOM technique to the Magnetization vector

See alsoCracknell, 2014, a SOM analysis of Australia



Installation

  1. Install Oasis Montaj, which includes an Anaconda3 Python distribution.

  2. Install the dependent "pyqt5" package as Administrator:

    1. In Start Menu, find and click left mouse button (LMB) on "Anaconda3 (64-bit)"
    2. Click right mouse button (RMB) on "Python Prompt" and select "More" > "Open file location" to show these shortcuts in Windows Explorer

    3. In Windows Explorer, RMB on "Python Prompt" and select "Run as Administrator" to open the Python command prompt

    4. Run the following command 

      Code Block
      pip install pyqt5


  3. Download self_organizing_maps_9_4.zip, extract the content to your Geosoft user\python folder (%USERPROFILE%\Documents\Geosoft\Desktop Applications 9\python)

  4. Open your Oasis montaj project, create a Geosoft database with multivariate data, display the data channels to be analysed, and only those channels.

  5. Run the som_om.py file (Run a GX, change the type to *.py, navigate to C:\Program Files\Geosoft\Desktop Applications 9\user\python and select som_om.py):

    Image Added

  6. Optionally add an entry to your custom menu, so that you don't have to browse to the user\python folder:

    Code Block
    ITEM "SOM Analysis" , som_om.py


If installed correctly you will see the SOM dialog:
Image Added

Usage Notes

This GX analyses multivariate data by grouping data into statistically meaningful groupings using SOM neural-network analysis. The technique is described here:

http://en.wikipedia.org/wiki/Self-organizing_map

In this implementation up to 16 separate data channels can be analysed. The first three channels are used to colour-code the resulting SOM neural network report. The report provides a simple visual impression of the classifications.

The size of the neural network is the square-root of the number of classifications, which can be any of 4,6,16,25,36,49,etc… Two new channels are added (or replaced) in the database:

Class – the classifications, 0 to number of classifications - 1

EuD - Euclidian distance of the point from the assigned neuron class, which is the closest based on Euclidian distance.

If an anomalous percent above 0 is specified, the percent of data that is furthest from the network (by Euclidian distance) is re-classified a second time on its own into a second SOM network of the same dimension such that there will be twice as many classes. For example, choosing 16 classifications and 5% anomalous data will result in 32 classes (16+16) where classes 0 to 15 capture the 95% of data that most closely fits the original network, and classes 16 to 31 are the most extreme 5% reclassified in their own network.

To work with gridded or voxel data:

Export a grid or voxel to a new database.

Sample other grids/voxels into new channels of the database.

Run the SOM analysis.

Save the "Class" channel as a grid or voxel.

Imagination encouraged!

For example, you can classify MVI results based on the x,y,z vector directions and the amplitude, which effectively double-weights the amplitude in the analysis. Because the directions can range + or -, the amplitude is scalar, and all dimensions are likely reasonably and meaningfully scaled relative to each other then no normalization is necessary. The process would be:

Create a database from the amplitude voxel.

Sample the x, Y and Z directions into mx,my,mz channels.

Run a SOM, dimension 9 (or 16), anomalous percent 5.

Convert Class channel to voxel and view.

Clip the data to above 9 (or 16), and this is the SOM of the anomalous locations.