Self Organizing Maps. We are drowning in information and starving for knowledge. A New Approach for Integrated Analysis of Geological Data.

Size: px

Start display at page:

Download "Self Organizing Maps. We are drowning in information and starving for knowledge. A New Approach for Integrated Analysis of Geological Data."

Joel Mosley
5 years ago
Views:

1 Radiometrics Processed Landsat TM Self Organizing Maps A New Approach for Integrated Analysis of Geological Data. We are drowning in information and starving for knowledge. Rutherford D. Roger Stephen.Fraser@csiro.au

2 Background #1 Explorationists/Geoscientists/Miners now gather data faster than it can be interpreted. Traditional multivariate statistical methods are confused by: Non-linear relationships Non-normal data distributions (non-gaussian) Missing or censored data (nulls) Need to analyze disparate or complex data types Categorical (text), continuous, and discontinuous types

3 Background #2 GIS/Mine Planning Packages enable data storage and display; but do not solve the problem of, How do we intelligently analyze and interpret the volumes of data we collect?

4 Clustering & Data Analysis: A brief history Classical Statistical Fisher s Discriminate Analysis, Least-Squares, Principal Components Analysis, Factor Analysis. Modern Statistical - more flexible methods, that estimate within, and between class probabilities: Nearest Neighbour, Projection Pursuit, Causal Networks, CART (Classification & Regression Trees), MARS (Multivariate Adaptive Regression Splines) Machine Learning automatic or logical systems based on logical or binary systems: Artificial Intelligence, Expert Systems, Decision Trees, Neural Nets (most ML is supervised!) Kohonen Nets - (Teuvo Kohonen, 1985) Self Organizing Maps (SOM) - A non-traditional, method of data analysis based on principles of vector quantization and measures of vector similarity. Developed as a computational method to counter neural nets. "Ordered Vector Quantization (Kohonen s preferred name)

5 Scatter Plots: An Effective Tool for Conceptualizing Data Processing methods Consider a grouping of similar/related samples in n-d space Sample point can be considered as a vector based on the contributions of the various constituent components a b Var C c CSIRO

6 Scatter Plots: An Effective Tool for Conceptualizing Data Processing methods a b Var C c CSIRO

7 Background: Self Organizing Maps nd Data Space Vectors A segmentation and visualization technique to explore relationships between diverse data types: Based on principles of vector quantization and measures of vector similarity; Can handle Non-linear relationships and Non-Gaussian data distributions; Can handle categorical (nominal) data and labels ; Can handle Nulls, hence sparse data, can be accommodated. Outputs 2D orderly-maps that represent the nd data structure and maintains the relationships between inputs. 2D Self Organized Map Representation of Samples: Colours indicate similarity or dissimilarity of adjacent nodes CSIRO

quantization and measures of vector similarity; Can handle Non-linear relationships and

Nulls, hence sparse data, can be accommodated.

8 Background: Self Organizing Maps nd Data Space Vectors A segmentation and visualization technique to explore relationships between diverse data types: Based on principles of vector quantization and measures of vector similarity; Can handle Non-linear relationships and Non-Gaussian data distributions; Can handle categorical (nominal) data and labels ; Can handle Nulls, hence sparse data, can be accommodated. Outputs 2D orderly-maps that represent the nd data structure and maintains the relationships between inputs. 2D Self Organized Map Representation of Samples: Colours indicate similarity or dissimilarity of adjacent nodes Variable Contributions for samples shown on the Map - Component Plots : colours indicate spread of values across range of input values. CSIRO

Hodgkinson Input Samples levelled gridded and logtransformed

9 Analysis of Geoscience BC s QUEST Stream & Lake Sediment Geochemical Database Stephen Fraser, Peter Kowalczyk & Jane Hodgkinson Input Samples levelled gridded and logtransformed geochemistry samples x 42 elements: (over ~150,000 km 2 ) CSIRO

10 U-Matrix & K-means2 20 clusters QUEST_all_geochem_levelled_17Apr2009plk_sjf1_35x29_toroid_Km2_20

11 Surficial Geology vs Samples Coded by SOM-derived K-means (20 clusters) Surficial Geology SOM-derived K-means

12 CSIRO K-means Cluster Normalized Elemental Maps

13 Component Plots for each of the 42 elements

14 Cluster-Normalized Elemental Maps Au example Maps for each element available from GBC web site

15 CSIRO Component Plot - Element Distributions

16 Au Component Plot Six groups of BMU's are identified as anomalous in gold These BMU's are selected, and then samples from that BMU are plotted on a map.

17 Sample Sites & Au Component Plot Between Prince George and Quesnel & Mt Polley /Gibraltar

18 Group 1 of the Anomalous Au BMUs spatially plot between Prince George and Quesnel & Mt Polley /Gibraltar

19 Group 2 of the Anomalous Au BMUs

20 Group 3 of the Anomalous Au BMUs

21 Group 4 of the Anomalous Au BMUs

22 Group 5 of the Anomalous Au BMUs

23 Group 6 of the Anomalous Au BMUs Plots around Mt. Milligan

24 CSIRO Cross-plots of Selected Elements

25 Cross-Plot of Au vs Ag SOM Node values Three elevated Au associations evident QUEST_all_geochem_levelled_17Apr2009plk_sjf1_35x29_toroid_Km2_20

26 Spatial Distribution of Selected Ag vs Au Samples High Au

27 Brownfields Regional Geochemical Database -Osborne SOM used to identify populations related to processes SOM used to identify prospective targets Bruce Dickson, Peter Kowalczyk, Gary Sparks ~ 40,000 located (XYZ) geochemical samples with up to 13 elements assayed: ~ 60% of data base is empty CSIRO

$A Demonstration of SOM where a large fraction (60%) of data is missing.$

28 Exploration around the Osborne Cu-Au Deposit, Queensland, Australia. A Demonstration of SOM where a large fraction (60%) of data is missing. CSIRO

29 Component Plots CSIRO Component Plots #1

30 SOM Node - Depth Plot Z 260 AUPM vs depth AUPM Au Vs Depth CSIRO

SOM Node - Depth Plot Z 260 AUPM vs depth 250 240 230 Detrital 220 210 200 Basement/Cover

31 SOM Node - Depth Plot Z 260 AUPM vs depth Detrital Basement/Cover Interface AUPM Au Vs Depth CSIRO SOM can assist in the identification of process

32 Au group I: Detrital Houdini Au component map Kultho r Osborne CSIRO

33 CSIRO Au All Looking ---- North

34 SOM on Voxel Volumes of Petrophysical Data resulting from Unconstrained Geophysical Inversion Seagull (PGE-Ni-Cu) Geophysical Model Data from: Magnetic and Gravity Three-Dimensional (3D) Modelling: Lake Nipigon Region Geoscience Initiative L.E. Reed, D.R.B. Rainsford CSIRO Ontario Ministry of Northern Development and Mines, Ontario Geological Survey MRD report

Seagull Location Lake Nipigon Seagull CSIRO Seagull

[from Magnetic and Gravity Three-Dimensional (3D)

35 Seagull Location Lake Nipigon Seagull CSIRO Seagull Location of Seagull and other magnetic models in area [from Magnetic and Gravity Three-Dimensional (3D) Modelling: Lake Nipigon, Region Geoscience Initiative, L.E. Reed, D.R.B. Rainsford, Ontario, Ministry of Northern Development and Mines, Ontario Geological Survey, MRD report 193, 2006]

Geology draped over Magnetics TMI CSIRO Magnetic model for Seagull. Overhead view of geology draped on amplitude enhanced total field magnetics.

36 Geology draped over Magnetics TMI CSIRO Magnetic model for Seagull. Overhead view of geology draped on amplitude enhanced total field magnetics. [from Magnetic and Gravity Three-Dimensional (3D) Modelling: Lake Nipigon, Region Geoscience Initiative, L.E. Reed, D.R.B. Rainsford, Ontario, Ministry of Northern Development and Mines, Ontario Geological Survey, MRD report 193, 2006]

37 Raw Data Magnetic Susceptibility Density CSIRO

38 SOM on 411,156 samples each: MagSus and Density sf10k_sus_den_12x9_toroid_2_km2_9 CSIRO

39 SOM on 411,156 samples: MagSus and Density CSIRO sf10k_sus_den_12x9_toroid_2_km2_9

40 CSIRO sf10k_sus_den_12x9_toroid_2_km2_9

41 Outlier Analysis Based on measures of how far away a sample is from its SOM Code Vector - QER a b Var C c CSIRO

42 X-Y Plot Quantization Errors CSIRO

43 X-Y Plot Quantization Errors CSIRO

44 X-Y Plot Quantization Errors CSIRO

45 X-Y Plot Quantization Errors CSIRO

46 X-Y Plot Quantization Errors CSIRO

47 X-Y Plot Quantization Errors CSIRO

48 X-Y Plot Quantization Errors CSIRO

49 X-Y Plot Quantization Errors CSIRO

50 X-Y Plot Quantization Errors CSIRO

51 X-Y Plot Quantization Errors CSIRO

52 X-Z Plot Quantization Errors CSIRO

53 SOM on Voxel Volumes of Petrophysical Data resulting from Unconstrained Geophysical Inversion SiroSOM Analysis of UBC-Inverted Geophysical Voxel Data in the Mt Isa Region for Geological Survey of Queensland. Jane H Hodgkinson, Stephen J Fraser January-February 2009 gravgrd to mag list no negs_35x30_toroidkm2_5 CSIRO

54 Location and Input Data Sets Grav Grav Mag Mag CSIRO.

55 U-Matrix and K-means 5 Classes CSIRO.

56 Block Perspectives looking towards North BMUs Clusters Q-errors Grav Mag CSIRO.

57 Block Perspectives looking towards North BMUs Clusters Q-errors Grav Mag CSIRO.

58 Block Perspectives looking towards North BMUs Clusters Q-errors Grav Mag CSIRO.

59 BMUs Clusters Q-errors Grav Mag CSIRO.

60 BMUs Clusters Q-errors Grav Mag CSIRO.

61 BMUs Clusters Q-errors Grav Mag CSIRO.

62 BMUs Clusters Q-errors Grav Mag CSIRO.

63 BMUs Clusters Q-errors Grav Mag CSIRO.

64 BMUs Clusters Q-errors Grav Mag CSIRO.

65 BMUs Clusters Q-errors Grav Mag CSIRO.

66 BMUs Clusters Q-errors Grav Mag CSIRO.

67 BMUs Clusters Q-errors Grav Mag CSIRO.

68 BMUs Clusters Q-errors Grav Mag CSIRO.

69 BMUs Clusters Q-errors Grav Mag CSIRO.

70 BMUs Clusters Q-errors Grav Mag CSIRO.

71 BMUs Clusters Q-errors Grav Mag CSIRO.

72 BMUs Clusters Q-errors Grav Mag CSIRO.

73 BMUs Clusters Q-errors Grav Mag CSIRO.

74 BMUs Clusters Q-errors Grav Mag CSIRO.

75 Conclusions SOM sub-populations can show process or allow targeting of samples with specific characteristics; SOM can assist in identifying: stream and lake sediment geochemical targets; targets and domains in petrophysical volumes. Complex and diverse data types can be analyzed using SOM; Geochemical, Geophysical, Mineralogical (Geotechnical, Geometallurgical.) Any spatial data with complex or disparate inputs

Application and Challenges of Artificial Intelligence in Exploration

Application and Challenges of Artificial Intelligence in Exploration XPLOR 2017 John McGaughey 2017 Mira Geoscience Ltd. Artificial Intelligence Artificial Intelligence is colossally hyped these days,