WO2024077330A1 - Continuum estimation tool for mining of geological structures - Google Patents
Continuum estimation tool for mining of geological structures Download PDFInfo
- Publication number
- WO2024077330A1 WO2024077330A1 PCT/AU2023/050927 AU2023050927W WO2024077330A1 WO 2024077330 A1 WO2024077330 A1 WO 2024077330A1 AU 2023050927 W AU2023050927 W AU 2023050927W WO 2024077330 A1 WO2024077330 A1 WO 2024077330A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- continuum
- geological
- analysis method
- geological analysis
- abundance
- Prior art date
Links
- 238000005065 mining Methods 0.000 title claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 23
- 239000000463 material Substances 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000000985 reflectance spectrum Methods 0.000 claims description 22
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 10
- 239000011707 mineral Substances 0.000 claims description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 238000000701 chemical imaging Methods 0.000 claims description 2
- 238000005286 illumination Methods 0.000 description 13
- 238000007781 pre-processing Methods 0.000 description 9
- 239000011159 matrix material Substances 0.000 description 8
- 238000013459 approach Methods 0.000 description 6
- 239000000203 mixture Substances 0.000 description 5
- 238000000862 absorption spectrum Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000009472 formulation Methods 0.000 description 3
- 238000010801 machine learning Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 239000013626 chemical specie Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000004800 variational method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/40—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for measuring magnetic field characteristics of the earth
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
- G01J2003/2826—Multispectral imaging, e.g. filter imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/34—Transmitting data to recording or processing apparatus; Recording data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V99/00—Subject matter not provided for in other groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/901—Indexing; Data structures therefor; Storage structures
- G06F16/9017—Indexing; Data structures therefor; Storage structures using directory or table look-up
Definitions
- the present invention generally relates to mining of geological structures.
- a mobile mining spectral scanner can be used to perform hyperspectral scanning of a mine including scanning of mine faces, muck piles, cores, stockpiles and other geological structures.
- the preferred embodiment provides a useful analytical tool for such a purpose.
- a geological analysis method for determining abundance results for one or more valuable materials in a geological structure the method involving: receiving scan data relating to a geological structure; processing the data to determine a continuum; and generating an abundance result of the valuable materials using the continuum.
- the abundance result can be used for mining applications whereby locations of the geological structure having higher abundance are more desirable to mine.
- the method is automated.
- the method may involve capturing the scan data using a hyperspectral imaging device.
- the method may involve forming a hyperspectral data cube.
- the hyperspectral data cube may include reflectance spectra over a spatial period.
- the continuum may be determined using a heuristic, and preferably a greedy heuristic.
- the determined continuum may be optimized based upon minimal curvature and/or minimal offset from reflectance spectra.
- the continuum may be removed from reflectance spectra.
- the continuum may be determined by solving a least squares problem.
- the continuum may be determined by iteratively solving least squares problems.
- the continuum may be determined by splitting a solution vector into a passive and an active set, and iteratively updating the passive set with a highest residual.
- the continuum may be determined by using pre-calculated QR decomposition for at least one least squares problem.
- the abundance result may include a mineral percentage at a location of the geological structure.
- geological analysis system for determining abundance results for one or more valuable materials in a geological structure, the system including: a scanner for capturing scan data relating to a geological structure; a processor for processing the data to determine a continuum; and a generator for generating an abundance result of the valuable materials using the continuum.
- the scanner may include a hyperspectral camera.
- the scanner may be mobile.
- Figure 1 is a schematic diagram of a geological analysis system in accordance with an embodiment of the present invention.
- Figure 2 is a higher level flowchart of a geological analysis method performed using the system of Figure 1 ;
- Figure 3 is a lower level flowchart of the illumination pre-processing step of the higher level flowchart of Figure 2;
- Figure 4 shows outputs of the illumination pre-processing step of Figure 3 at different iterations.
- the system 100 includes a mobile hyperspectral scanner 102, which is an imaging device for capturing scan data 103 relating to a mine 104.
- the mine 104 may include mine faces, muck piles, core, stockpiles and other like geological structures, including minerals or other valuable materials.
- the system 100 also includes a base processor 106 for processing of the hyperspectral scan data 103 to determine a continuum.
- the processor 106 i.e. generator
- the processor 106 also generates an abundance result of the minerals using the determined continuum.
- An electronic display 108 is provided for displaying the abundance result to a geologist 110 who can then determine the best place to mine.
- the base processor 106 is typically in wireless communication with the remote scanner 102, and is also networked to the cloud for enabling the geologist 110 to utilise stored scan data 103 and software via a web-based portal.
- Figure 2 is a flowchart of a geological analysis method 200 performed using the system 100.
- the automated method 200 is used to determine abundance results for one or more minerals in the mine 104.
- the scan data 103 from the scanner 102 and relating to the mine 104, is received by the processor 106.
- the mobile scanner 102 drives around the mine 104 capturing the data at various locations of interest and a hyperspectral data cube is formed.
- the hyperspectral data cube includes reflectance spectra over a spatial period.
- the method 200 involves illumination pre-processing the data 103 to determine a continuum. Illumination pre-processing of the hyperspectral data (reflectance spectra) produces an illumination spectrum using either known methods, or the method detailed below.
- the processor 106 generates an abundance result of the minerals using the continuum.
- the processor 106 continues processing the data 103, using machine learning with the illumination continuum as input, to produce abundance percentages useful for mining applications.
- the processor 106 takes the output of the illumination pre-processing step 204 (in this case the optimal fit of the continuum, removed from reflectance spectra) and generates a result in the form of abundances of minerals.
- the abundance result is displayed on the display 108 or can be printed for the geologist 110.
- the abundance result includes a mineral percentage at a location of the mine 104.
- the abundance result can be used for mining applications whereby locations of the geological structure having higher abundance are more desirable to mine.
- Illumination pre-processing A known step in the processing of hyperspectral data analysis is to estimate and remove the illumination continuum from the reflectance spectra.
- the continuum line connects the local maximas of the reflectance spectra and is frequently estimated using the method developed by Clark & Roush (1984) and summarised by Clark (1999) and improved upon by Clark et al. (2003).
- a first approach to getting reasonable runtimes is using a Greedy Heuristic (see Appendix I below).
- the solution vector is split into a passive and an active set (which is a common approach in Quadratic Programming algorithms).
- the passive set is the set of solution vector elements that are precisely zero and the active set includes the parts that are non-zero.
- the method then involves iterating multiple times based on a heuristic.
- Figure 3 shows the illumination pre-processing step 204.
- setting an element of the solution vector to zero is equivalent to cutting a row out of the A matrix.
- Steps 204.1 and 204.6 involve initialising the working points set to be empty (i.e. the whole reflectance spectra is contained in the active set).
- Steps 204.2 form the problem which is solved using a least squares via a QR factorization. This can then use standard methods, such as case Bjorck (1996), to calculate the solution.
- QRUpdate library (20220)
- steps 204.2 and 204.5 can be used in implementing steps 204.2 and 204.5.
- Steps 204.3 and 204.4 utilise the “greedy heuristic” approach by iteratively updating the passive set with the highest residual to provide speed optimisation.
- Step 204.5 exploits the fact that, in the least squares matrix used to solve step 204.2, setting an element of the solution vector to zero is equivalent to removing a column out of the active matrix. Since the matrix A and the result vector b do not change within iterations (see Appendix I below), the method uses the pre-calculated QR decomposition that is needed for each least squares iteration.
- Figure 4 shows outputs of the illumination pre-processing step of Figure 3 at different iterations.
- the method can pre calculate the QR decomposition that is needed for each least squares iteration. As the passive set is updated, this is equivalent to removing rows from the A matrix.
- These updates are performed directly to the QR factorisation via Bjdrck’s corrected semi-normal equation (CSNE) approach (REF). This greatly speeds up the solution time as it avoids the need for a full matrix solve at every iteration (as the rate limiting factorisation step is avoided).
- CSNE corrected semi-normal equation
- Von Neumann boundary conditions can be used for continuity of derivate at the ends to the signal.
- the regularization parameter can be tuned based on repeated scans of the same sample. Runtime was observed down to approximately 30 seconds on a single thread for a Chip scan. It was implemented to make use of multithreading and saturated at about 4 threads at approximately 2 times quicker. This physics informed normalisation procedure for the samples enables results of this decomposition to be fed into a subsequent machine learning (ML) processing pipeline.
- ML machine learning
- the foregoing method is used for preprocessing the Hyperspectral data that is subsequently used for regressing XRF chemical species percentages from hyperspectral images that is performed.
- the crux problem in turning diffuse reflectance spectra into absorption spectra for mineral species classification is the determination of an unknown “continuum”. This is the remaining unknown after calibration for the known illumination spectra of a scan.
- this is typically estimated via the fitting of a convex hull to the reflectance spectra. There is no solid theoretical foundation for this assumption of convexity or linearity in the continuum.
- the continuum is strictly greater than the reflectance spectrum, i.e. C( ⁇ ) ⁇ p( ⁇ ) ⁇
- the continuum is a smoother function than the absorption spectrum, which is typically composed of narrow absorption bands (where we as- sume that hold the information about material composition).
- equation (5) is simply a second order low pass filter on the reflectance spectrum and may be solved by a variety of simple methods on a discrete wavelength basis.
- the positivity constraint changes the nature of the problem however, it can be shown that this is instead a Quadratic Program via formulation of the dual form (REF), i.e.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Geophysics (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The present invention relates to a geological analysis method for determining abundance results for one or more valuable materials in a geological structure. The method involves receiving scan data relating to a geological structure, and processing the data to determine a continuum. The method further involves generating an abundance result of the valuable materials using the continuum. Advantageously, the abundance result can be used for mining applications whereby locations of the geological structure having higher abundance are more desirable to mine.
Description
CONTINUUM ESTIMATION TOOL FOR MINING OF GEOLOGICAL STRUCTURES
TECHNICAL FIELD
[0001] The present invention generally relates to mining of geological structures.
BACKGROUND
[0002] The reference to any prior art in this specification is not, and should not be taken as an acknowledgement or any form of suggestion that the prior art forms part of the common general knowledge.
[0003] A mobile mining spectral scanner can be used to perform hyperspectral scanning of a mine including scanning of mine faces, muck piles, cores, stockpiles and other geological structures.
[0004] It is desirable to mine geological structures where there is a higher abundance of minerals.
[0005] The preferred embodiment provides a useful analytical tool for such a purpose.
SUMMARY OF THE INVENTION
[0006] According to one aspect of the present invention, there is provided a geological analysis method for determining abundance results for one or more valuable materials in a geological structure, the method involving: receiving scan data relating to a geological structure; processing the data to determine a continuum; and generating an abundance result of the valuable materials using the continuum.
[0007] Advantageously, the abundance result can be used for mining applications whereby locations of the geological structure having higher abundance are more desirable to mine.
[0008] Preferably, the method is automated.
[0009] The method may involve capturing the scan data using a hyperspectral imaging device. The method may involve forming a hyperspectral data cube. The hyperspectral data cube may include reflectance spectra over a spatial period.
[00010] The continuum may be determined using a heuristic, and preferably a greedy heuristic. The determined continuum may be optimized based upon minimal curvature and/or minimal offset from reflectance spectra. The continuum may be removed from reflectance spectra. The continuum may be determined by solving a least squares problem. The continuum may be determined by iteratively solving least squares problems. The continuum may be determined by splitting a solution vector into a passive and an active set, and iteratively updating the passive set with a highest residual.
[00011] The continuum may be determined by using pre-calculated QR decomposition for at least one least squares problem.
[00012] The abundance result may include a mineral percentage at a location of the geological structure.
[00013] According to another aspect of the present invention, there is provided geological analysis system for determining abundance results for one or more valuable materials in a geological structure, the system including: a scanner for capturing scan data relating to a geological structure; a processor for processing the data to determine a continuum; and a generator for generating an abundance result of the valuable materials using the continuum.
[00014] The scanner may include a hyperspectral camera. The scanner may be mobile.
[00015] Any of the features described herein can be combined in any combination with any one or more of the other features described herein within the scope of the invention.
BRIEF DESCRIPTION OF THE DRAWINGS
[00016] Preferred features, embodiments and variations of the invention may be discerned from the following Detailed Description which provides sufficient information for those skilled in the art to perform the invention. The Detailed Description is not to be regarded as limiting the scope of the preceding Summary of the Invention in any way. The Detailed Description will make reference to a number of drawings as follows:
[00017] Figure 1 is a schematic diagram of a geological analysis system in accordance with an embodiment of the present invention;
[00018] Figure 2 is a higher level flowchart of a geological analysis method performed using the system of Figure 1 ;
[00019] Figure 3 is a lower level flowchart of the illumination pre-processing step of the higher level flowchart of Figure 2; and
[00020] Figure 4 shows outputs of the illumination pre-processing step of Figure 3 at different iterations.
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS
[00021] According to an embodiment of the present invention, there is provided a geological analysis system 100 or tool as shown in Figure 1. The system 100 includes a mobile hyperspectral scanner 102, which is an imaging device for capturing scan data 103 relating to a mine 104. The mine 104 may include mine faces, muck piles, core, stockpiles and other like geological structures, including minerals or other valuable materials.
[00022] The system 100 also includes a base processor 106 for processing of the hyperspectral scan data 103 to determine a continuum. The processor 106 (i.e. generator) also generates an abundance result of the minerals using the determined continuum.
[00023] An electronic display 108 is provided for displaying the abundance result to a geologist 110 who can then determine the best place to mine. The base processor 106
is typically in wireless communication with the remote scanner 102, and is also networked to the cloud for enabling the geologist 110 to utilise stored scan data 103 and software via a web-based portal.
[00024] Overview
[00025] Figure 2 is a flowchart of a geological analysis method 200 performed using the system 100. The automated method 200 is used to determine abundance results for one or more minerals in the mine 104.
[00026] At step 202, the scan data 103, from the scanner 102 and relating to the mine 104, is received by the processor 106. The mobile scanner 102 drives around the mine 104 capturing the data at various locations of interest and a hyperspectral data cube is formed. The hyperspectral data cube includes reflectance spectra over a spatial period.
[00027] At step 204, the method 200 involves illumination pre-processing the data 103 to determine a continuum. Illumination pre-processing of the hyperspectral data (reflectance spectra) produces an illumination spectrum using either known methods, or the method detailed below.
[00028] At step 206, the processor 106 generates an abundance result of the minerals using the continuum. The processor 106 continues processing the data 103, using machine learning with the illumination continuum as input, to produce abundance percentages useful for mining applications. The processor 106 takes the output of the illumination pre-processing step 204 (in this case the optimal fit of the continuum, removed from reflectance spectra) and generates a result in the form of abundances of minerals.
[00029] At step 208, the abundance result is displayed on the display 108 or can be printed for the geologist 110. The abundance result includes a mineral percentage at a location of the mine 104. Advantageously, the abundance result can be used for mining applications whereby locations of the geological structure having higher abundance are more desirable to mine.
[00030] Illumination pre-processing
[00031] A known step in the processing of hyperspectral data analysis is to estimate and remove the illumination continuum from the reflectance spectra. The continuum line connects the local maximas of the reflectance spectra and is frequently estimated using the method developed by Clark & Roush (1984) and summarised by Clark (1999) and improved upon by Clark et al. (2003).
[00032] There are other various methods published in the scientific literature, including Mutanga et al. (2004).
[00033] Instead of using the foregoing Convex Hull method, a different method approach is developed which is more complex, aligns with underlying physics and gives better results. A model of the continuum is built and optimized having: (1) minimal curvature, (2) minimal offset from the captured reflectance spectra. In an ideal world this leads to a low pass filter to fit the reflectance spectra using off-the-shelf quadratic program algorithms, this however is computationally slow.
[00034] Part A - Greedy Heuristic to Improve speed
[00035] A first approach to getting reasonable runtimes is using a Greedy Heuristic (see Appendix I below). The solution vector is split into a passive and an active set (which is a common approach in Quadratic Programming algorithms). The passive set is the set of solution vector elements that are precisely zero and the active set includes the parts that are non-zero. The method then involves iterating multiple times based on a heuristic.
[00036] Figure 3 shows the illumination pre-processing step 204. In the least squares formulation, setting an element of the solution vector to zero is equivalent to cutting a row out of the A matrix.
[00037] The greedy heuristic steps are as follows:
204.1. Initialise the passive set to the empty set (i.e. the whole solution vector is active)
204.2. Solve the least squares problem via a QR factorization
204.3. Find the residual with the highest positive value
204.4. Add this residual to the passive set
204.5. Solve the least squares problem on this smaller problem
204.6. Iterate this procedure until there are no more positive residuals (with a tolerance).
[00038] Steps 204.1 and 204.6 involve initialising the working points set to be empty (i.e. the whole reflectance spectra is contained in the active set).
[00039] Steps 204.2 form the problem which is solved using a least squares via a QR factorization. This can then use standard methods, such as case Bjorck (1996), to calculate the solution. A known library, QRUpdate library (20220), can be used in implementing steps 204.2 and 204.5.
[00040] Steps 204.3 and 204.4 utilise the “greedy heuristic” approach by iteratively updating the passive set with the highest residual to provide speed optimisation.
[00041] Part B - Iterative updates to a pre-computed AR Set
[00042] Step 204.5 exploits the fact that, in the least squares matrix used to solve step 204.2, setting an element of the solution vector to zero is equivalent to removing a column out of the active matrix. Since the matrix A and the result vector b do not change within iterations (see Appendix I below), the method uses the pre-calculated QR decomposition that is needed for each least squares iteration.
[00043] These updates are performed directly to the pre-calculated QR factorisation which was original output via step 204.2, but is taken as an input to Bjdrck’s (1996) corrected semi-normal equation (CSNE) approach on each subsequent iteration. This greatly speeds up the solution time as it avoids the need for a full matrix solve at every iteration (as the rate limiting factorisation step is avoided).
[00044] Figure 4 shows outputs of the illumination pre-processing step of Figure 3 at different iterations.
[00045] Since the matrix A and the result vector b do not change within iterations, the method can pre calculate the QR decomposition that is needed for each least squares iteration. As the passive set is updated, this is equivalent to removing rows from the A matrix. These updates are performed directly to the QR factorisation via Bjdrck’s corrected semi-normal equation (CSNE) approach (REF). This greatly speeds up the solution time as it avoids the need for a full matrix solve at every iteration (as the rate limiting factorisation step is avoided).
[00046] Von Neumann boundary conditions can be used for continuity of derivate at the ends to the signal. The regularization parameter can be tuned based on repeated scans of the same sample. Runtime was observed down to approximately 30 seconds on a single thread for a Chip scan. It was implemented to make use of multithreading and saturated at about 4 threads at approximately 2 times quicker. This physics informed normalisation procedure for the samples enables results of this decomposition to be fed into a subsequent machine learning (ML) processing pipeline.
[00047] In conclusion, the foregoing method is used for preprocessing the Hyperspectral data that is subsequently used for regressing XRF chemical species percentages from hyperspectral images that is performed. The crux problem in turning diffuse reflectance spectra into absorption spectra for mineral species classification is the determination of an unknown “continuum”. This is the remaining unknown after calibration for the known illumination spectra of a scan. In Hyperspectral literature this is typically estimated via the fitting of a convex hull to the reflectance spectra. There is no solid theoretical foundation for this assumption of convexity or linearity in the continuum. As a counterpoint to this modelling assumption a basic application Kubelka Munk theory of scattering from infinetesimal layers suggests a model that is the reciprocal of a low order polynomial for the continuum [REF]. More principled ways of estimating the continuum via variational methods are provided herein.
[00048] References
Bjorck, A. (1996). Numerical methods for least squares problems. SIAM.
Clark, R.N., Roush, T.L. (1984). Reflectance spectroscopy: Quantitative analysis techniques for remote sensing applications.: Journal of Geophysical Research: Solid Earth 89, 6329-6340.
Clark, R. N. (1999). Spectroscopy of rocks and minerals and principles of spectro-scopy. In Manual of Remote Sensing, edited by A. N. Rencz, John Wiley, New York.
Clark, R. N., King, T. V. V. and Gorelick, N. S. (1987): Automatic continuum analysis of reflectance spectra. Proceedings of the Third Airborne Imaging Spectrometer Data Analysis Workshop, 30. 138-142.
Friedlander, M., Le Provost, M., Saba, E. (2022). QRUpdate Julia library (v1.0.00 [Computer software]. GitHub. httpsV/github.com/mpf/QRupdate.jl
Mutanga, O. and Skidmore, A. K. (2004). Hyperspectral band depth analysis for a better estimation of grass biomass (Cenchrus ciliaris) measured under controlled laboratory conditions International Journal of applied Earth Observation and Geoinformation, 5, 87-96.
Clark, R. N., G. A. Swayze, K. E. Livo, R. F. Kokaly, S. J. Sutley, J. B. Dalton, R. R. McDougal, and C. A. Gent, Imaging spectroscopy. (2003) Earth and planetary remote sensing with the USGS Tetracorder and expert systems, J. Geophys. Res., 108(E12), 5131 , doi:10.1029/2002JE001847.
[00049] A person skilled in the art will appreciate that many embodiments and variations can be made without departing from the ambit of the present invention.
[00050] In compliance with the statute, the invention has been described in language more or less specific to structural or methodical features. It is to be understood that the invention is not limited to specific features shown or described since the means herein described comprises preferred forms of putting the invention into effect.
[00051] Reference throughout this specification to 'one embodiment’ or 'an embodiment’ means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the present invention. Thus, the appearance of the phrases 'in one embodiment’ or 'in an embodiment’ in various places throughout this specification are not necessarily all referring to the same embodiment. Furthermore, the particular features, structures, or characteristics may be combined in any suitable manner in one or more combinations.
APPENDIX I
Reflectance Spectra Decomposition
A fundamental assumption of this method is that the illumination corrected reflectance spectra that we sample, p(λ), are the product of a wavelength dependent absorption spectrum, K(λ), and a continuum C(λ), i.e. p(λ) = K(λ)C(λ) (1)
The question is then framed as follows,
“how do we estimate the continuum K(λ) for a given reflectance spectra p(λ), without knowing the absorption spectra K(λ)?”.
This is essentially a blind source separation problem where we can use known attributes of the absorption and continuum spectra to pull these two signals
apart. The two crucial known pieces of information that allow for such a de- composition are
1. The continuum is strictly greater than the reflectance spectrum, i.e. C(λ) ≥ p(λ) ∀λ
2. The continuum is a smoother function than the absorption spectrum, which is typically composed of narrow absorption bands (where we as- sume that hold the information about material composition).
To make progress towards this decomposition we can make the following mod- elling assumptions
Using the above desiderata and assuming the 12 norm in (3) and (4) as a linear objective yields the following optimisation problem
where the regularisation parameter vi(A) has been introduced to control the smoothness of the estimated continuum.
Without the constraint that ο(λ) ≥ 0, equation (5) is simply a second order low pass filter on the reflectance spectrum and may be solved by a variety of simple methods on a discrete wavelength basis. The positivity constraint changes the nature of the problem however, it can be shown that this is instead a Quadratic Program via formulation of the dual form (REF), i.e.
Note that in our case on a discrete wavelength basis
where x is the discrete basis o and Ais the finite difference operator.
Naive solution of a Quadratic Program requires exponential runtime, with quadratic complexity for some algorithms for positive definite Q (note we are only guaranteed semi-positive definite for our formulation of the problem). Both
Claims
1. A geological analysis method for determining abundance results for one or more valuable materials in a geological structure, the method involving: receiving scan data relating to a geological structure; processing the data to determine a continuum; and generating an abundance result of the valuable materials using the continuum.
2. A geological analysis method as claimed in claim 1 , wherein the abundance result can be used for mining applications whereby locations of the geological structure having higher abundance are more desirable to mine.
3. A geological analysis method as claimed in claim 1 , involving capturing the scan data using a hyperspectral imaging device.
4. A geological analysis method as claimed in claim 3, involving forming a hyperspectral data cube.
5. A geological analysis method as claimed in claim 4, wherein the hyperspectral data cube includes reflectance spectra over a spatial period.
6. A geological analysis method as claimed in claim 1 , wherein the continuum is determined using a heuristic.
7. A geological analysis method as claimed in claim 6, wherein the heuristic is a greedy heuristic.
8. A geological analysis method as claimed in claim 1 , wherein determined continuum is optimized based upon minimal curvature and/or minimal offset from reflectance spectra.
9. A geological analysis method as claimed in claim 1 , wherein the continuum is removed from reflectance spectra.
10. A geological analysis method as claimed in claim 1 , wherein the continuum is determined by solving a least squares problem.
11. A geological analysis method as claimed in claim 1 , wherein the continuum is determined by iteratively solving least squares problems.
12. A geological analysis method as claimed in claim 1 , wherein the continuum is determined by splitting a solution vector into a passive and an active set, and iteratively updating the passive set with a highest residual.
13. A geological analysis method as claimed in claim 1 , wherein the continuum is determined by using pre-calculated QR decomposition for at least one least squares problem.
14. A geological analysis method as claimed in claim 1 , wherein the abundance result includes a mineral percentage at a location of the geological structure.
15. A geological analysis system for determining abundance results for one or more valuable materials in a geological structure, the system including: a scanner for capturing scan data relating to a geological structure; a processor for processing the data to determine a continuum; and a generator for generating an abundance result of the valuable materials using the continuum.
16. A geological analysis system as claimed in claim 1 , wherein the scanner includes a hyperspectral camera.
17. A geological analysis system as claimed in claim 1 , wherein the scanner is mobile.
18. A geological analysis system as claimed in claim 1 , wherein the geological structure includes any one of a mine face, muck pile, core, or stockpile.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU2022902975 | 2022-10-11 | ||
AU2022902975A AU2022902975A0 (en) | 2022-10-11 | Continuum estimation tool for mining of geological structures |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2024077330A1 true WO2024077330A1 (en) | 2024-04-18 |
Family
ID=90668352
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/AU2023/050927 WO2024077330A1 (en) | 2022-10-11 | 2023-09-25 | Continuum estimation tool for mining of geological structures |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2024077330A1 (en) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103175801A (en) * | 2013-03-14 | 2013-06-26 | 中国国土资源航空物探遥感中心 | Large-batch automatic hyperspectral remote sensing mineral mapping method |
-
2023
- 2023-09-25 WO PCT/AU2023/050927 patent/WO2024077330A1/en unknown
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103175801A (en) * | 2013-03-14 | 2013-06-26 | 中国国土资源航空物探遥感中心 | Large-batch automatic hyperspectral remote sensing mineral mapping method |
Non-Patent Citations (2)
Title |
---|
A.J. BROWN: "Spectral curve fitting for automatic hyperspectral data analysis", IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, IEEE, USA, vol. 44, no. 6, 1 June 2006 (2006-06-01), USA, pages 1601 - 1608, XP093160735, ISSN: 0196-2892, DOI: 10.1109/TGRS.2006.870435 * |
RONAN RIALLAND, CHARLES SOUSSEN, RODOLPHE MARION, VERONIQUE CARRÈRE: "OMP-Based Algorithm for Mineral Reflectance Spectra Deconvolution from Hyperspectral Images", IGARSS 2020, 2020 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM, IEEE, 26 September 2020 (2020-09-26) - 2 October 2020 (2020-10-02), pages 4918 - 4921, XP093160738, ISBN: 978-1-7281-6374-1, DOI: 10.1109/IGARSS39084.2020.9323461 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kruse et al. | Validation of DigitalGlobe WorldView-3 Earth imaging satellite shortwave infrared bands for mineral mapping | |
Li et al. | A compressive sensing and unmixing scheme for hyperspectral data processing | |
Altmann et al. | Supervised nonlinear spectral unmixing using a postnonlinear mixing model for hyperspectral imagery | |
CN104952050B (en) | High spectrum image adaptive de mixing method based on region segmentation | |
Aggarwal et al. | Hyperspectral unmixing in the presence of mixed noise using joint-sparsity and total variation | |
Lorenz et al. | Feature extraction for hyperspectral mineral domain mapping: A test of conventional and innovative methods | |
Zhang et al. | Assessing the impact of endmember variability on linear Spectral Mixture Analysis (LSMA): A theoretical and simulation analysis | |
Uezato et al. | A novel spectral unmixing method incorporating spectral variability within endmember classes | |
Kurz et al. | A review of hyperspectral imaging in close range applications | |
Kirk et al. | Photoclinometry made simple | |
CN104964931B (en) | A kind of martian surface enhydrite quantitative inversion method | |
WO2024077330A1 (en) | Continuum estimation tool for mining of geological structures | |
Nakhostin et al. | Joint anomaly detection and spectral unmixing for planetary hyperspectral images | |
Wolfe et al. | Hyperspectral analytics in envi target detection and spectral mapping methods | |
Chandra et al. | A cognitive method for building detection from high-resolution satellite images | |
Chabrillat et al. | Prediction of common surface soil properties using airborne and simulated EnMAP hyperspectral images: Impact of soil algorithm and sensor characteristic | |
Piroddi et al. | Towards the definition of a low-cost toolbox for qualitative inspection of painted historical vaults by means of modified DSLR cameras, open source programs and signal processing techniques | |
Wolfe et al. | Hyperspectral analytics in ENVI | |
Irgenfried et al. | A framework for storage, visualization and analysis of multispectral data | |
Bui et al. | Hybrid machine learning models for aboveground biomass estimations | |
Johnson et al. | VNIR-SWIR Superspectral Mineral Mapping: An Example from Cuprite, Nevada | |
Basener et al. | A detection-identification process with geometric target detection and subpixel spectral visualization | |
Reyes-Muñoz et al. | Mapping essential vegetation variables over europe using gaussian process regression and sentinel-3 data in Google Earth Engine | |
Brook et al. | Quantitative detection of sediment dust analog over green canopy using airborne hyperspectral imagery | |
Nidamanuri et al. | Understanding the unique spectral signature of winter rape |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 23875946 Country of ref document: EP Kind code of ref document: A1 |