CN113552075A - Gold geochemical element inversion method and system - Google Patents
Gold geochemical element inversion method and system Download PDFInfo
- Publication number
- CN113552075A CN113552075A CN202110856570.0A CN202110856570A CN113552075A CN 113552075 A CN113552075 A CN 113552075A CN 202110856570 A CN202110856570 A CN 202110856570A CN 113552075 A CN113552075 A CN 113552075A
- Authority
- CN
- China
- Prior art keywords
- data
- gold
- geochemical
- transformation
- screening
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
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/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
-
- 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/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
- G01N21/274—Calibration, base line adjustment, drift correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- Immunology (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Pure & Applied Mathematics (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Evolutionary Biology (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention relates to a gold geochemical element inversion method and a gold geochemical element inversion system. The gold geochemical element inversion method comprises the steps of determining screening data according to acquired gold geochemical element data, determining mathematical transformation data according to GF-5 hyperspectral data by adopting a mathematical transformation method, carrying out correlation analysis on the screening data and the mathematical transformation data to obtain correlation data, carrying out regression analysis on the correlation data to obtain a regression equation, carrying out inversion operation according to the regression equation and the mathematical transformation data to accurately obtain gold geochemical element inversion information, and accurately expressing the gold geochemical element inversion information on a remote sensing image according to the obtained gold geochemical element inversion information. And the large-range abnormal information is presumed according to the small-range geochemical data (the acquired gold geochemical element data and GF-5 hyperspectral data), so that the defect that field work is difficult to carry out in a complex area is effectively overcome, and a powerful scientific basis and support are provided for the prospecting.
Description
Technical Field
The invention relates to the field of remote sensing technology and mathematical statistics, in particular to a gold geochemical element inversion method and a gold geochemical element inversion system based on GF-5 hyperspectral data.
Background
The gold geochemical element is an important index for gold mine exploration, and the traditional gold geochemical element research usually obtains the abundance of the gold geochemical element by arranging wild points in a research area, collecting soil samples of the points and analyzing the soil samples in a laboratory. The method has high measurement precision but low efficiency, and simultaneously, the distribution state of the geochemical elements of the soil in the region can be obtained by utilizing an interpolation method, and the influence of abnormal values is large.
Whether the abundance of geochemical elements in a region can be obtained more quickly and efficiently is an important problem for studying by scholars at home and abroad. With the development of the current remote sensing technology towards high spectral resolution, high spatial resolution and high temporal resolution, the hyperspectral remote sensing technology can quickly acquire spectral information in an area range due to the advantages of the platform, and the advantages can fill the defects of the traditional method. The high-resolution five-number (GF-5) is a full-spectrum hyperspectral satellite which is independently researched and developed in China and comprehensively observes the atmosphere and the land in the world, and is an important mark for realizing the ground observation capability of high spectral resolution in China.
In view of this, how to construct a high-precision gold geochemical element inversion method by using GF-5 hyperspectral data becomes an urgent problem to be solved.
Disclosure of Invention
The invention aims to provide a gold geochemical element inversion method and a gold geochemical element inversion system so as to improve the accuracy of a gold geochemical element inversion result.
In order to achieve the purpose, the invention provides the following scheme:
a gold geochemical element inversion method comprising:
acquiring gold geochemical element data, and determining screening data according to the gold geochemical element data;
acquiring GF-5 hyperspectral data of an area to be inverted;
determining mathematical transformation data according to the GF-5 hyperspectral data by adopting a mathematical transformation method;
performing correlation analysis on the screening data and the mathematical transformation data to obtain correlation data;
carrying out regression analysis on the correlation data to obtain a regression equation;
and carrying out inversion operation according to the regression equation and the mathematical transformation data to obtain gold geochemical element inversion information.
Preferably, the obtaining gold geochemical element data and determining screening data according to the gold geochemical element data specifically comprises:
acquiring known gold geochemical element content data in an area to be inverted;
and screening the data of the content of the gold geochemical elements by adopting a hierarchical sampling method to obtain the screened data.
Preferably, the screening of the gold geochemical element content data by using a hierarchical sampling method to obtain the screening data specifically comprises the following steps:
arranging the gold geochemical element content data in a sequence from high metal element content to low metal element content to obtain the total number of samples;
layering the arranged geochemical element content data of the gold according to the total number of the samples;
and screening the data of the content of the gold geochemical elements in each layer according to the data screening mode that the content of the gold geochemical elements is sorted from high to low to obtain the screened data.
Preferably, the determining mathematical transformation data according to the GF-5 hyperspectral data by using a mathematical transformation method specifically includes:
performing radiometric calibration on the GF-5 hyperspectral data;
atmospheric correction is carried out on the GF-5 hyperspectral data after radiometric calibration, and correction data are obtained;
and performing mathematical transformation on the correction data to obtain the mathematical transformation data.
Preferably, the mathematically transforming the correction data to obtain the mathematically transformed data specifically includes:
performing first derivative transformation, second derivative transformation or continuum removal transformation on the correction data to obtain mathematical transformation data; the mathematical transformation data is first derivative transformation data, second derivative transformation data or continuum removal transformation data.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the gold geochemical element inversion method provided by the invention comprises the steps of determining screening data according to the acquired gold geochemical element data, determining mathematical transformation data according to GF-5 hyperspectral data by adopting a mathematical transformation method, carrying out correlation analysis on the screening data and the mathematical transformation data to obtain correlation data, carrying out regression analysis on the correlation data to obtain a regression equation, and finally carrying out inversion operation according to the regression equation and the mathematical transformation data to accurately obtain gold geochemical element inversion information, so that the gold geochemical element inversion information can be accurately expressed on a remote sensing image according to the obtained gold geochemical element inversion information. And the large-range abnormal information is presumed according to the small-range geochemical data (the acquired gold geochemical element data and GF-5 hyperspectral data), so that the defect that field work is difficult to carry out in a complex area is effectively overcome, and a powerful scientific basis and support are provided for the prospecting.
Corresponding to the gold geochemical element inversion method, the invention also provides the following technical scheme:
a gold geochemical element inversion system comprising:
the screening module is used for acquiring gold geochemical element data and determining screening data according to the gold geochemical element data;
the GF-5 hyperspectral data acquisition module is used for acquiring GF-5 hyperspectral data of an area to be inverted;
the mathematical transformation module is used for determining mathematical transformation data according to the GF-5 hyperspectral data by adopting a mathematical transformation method;
the correlation analysis module is used for carrying out correlation analysis on the screening data and the mathematical transformation data to obtain correlation data;
the regression analysis module is used for carrying out regression analysis on the correlation data to obtain a regression equation;
and the gold geochemical element inversion module is used for carrying out inversion operation according to the regression equation and the mathematical transformation data to obtain gold geochemical element inversion information.
Preferably, the screening module comprises:
the element content data acquisition unit is used for acquiring known gold geochemical element content data in the region to be inverted;
and the screening unit is used for screening the gold geochemical element content data by adopting a hierarchical sampling method to obtain the screening data.
Preferably, the screening unit includes:
the sorting subunit is used for sorting the gold geochemical element content data according to the sequence of the metal element content from high to low to obtain the total number of samples;
the layering subunit is used for layering the arranged gold geochemical element content data according to the total number of the samples;
and the screening subunit is used for screening the data of the gold geochemical element content of each layer according to the data of the gold geochemical element content in a mode of sequencing from high to low to obtain the screened data.
Preferably, the mathematical transformation module specifically includes:
the radiometric calibration unit is used for radiometric calibration of the GF-5 hyperspectral data;
the correction unit is used for carrying out atmospheric correction on the GF-5 hyperspectral data subjected to radiometric calibration to obtain correction data;
and the mathematical transformation unit is used for performing mathematical transformation on the correction data to obtain the mathematical transformation data.
Preferably, the mathematical transformation unit includes:
the mathematical transformation subunit is used for carrying out first derivative transformation, second derivative transformation or continuum removal transformation on the correction data to obtain mathematical transformation data; the mathematical transformation data is first derivative transformation data, second derivative transformation data or continuum removal transformation data.
The technical effect achieved by the gold geochemical element inversion system provided by the invention is the same as that achieved by the gold geochemical element inversion method, so that the details are not repeated herein.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a gold geochemical element inversion method provided by the present invention;
FIG. 2 is a flow chart of a first implementation of a gold geochemical element inversion method in accordance with an embodiment of the present invention;
FIG. 3 is a flow chart of a second implementation of a gold geochemical element inversion method in accordance with an embodiment of the present invention;
FIG. 4 is a diagram of a first derivative transformation result of an original image according to an embodiment of the present invention;
FIG. 5 is a diagram of second derivative transformation results of an original image according to an embodiment of the present invention;
FIG. 6 is a diagram illustrating a transformation result of removing original image continuum according to an embodiment of the present invention;
FIG. 7 is a correlation curve of gold geochemical elemental data with remote sensing data according to an embodiment of the invention;
FIG. 8 is a graph of the first derivative inversion results provided by an embodiment of the present invention;
FIG. 9 is a diagram of second derivative inversion results provided by an embodiment of the present invention;
FIG. 10 is a graph of continuum removal inversion results provided by embodiments of the present invention;
FIG. 11 is a schematic structural diagram of a gold geochemical element inversion system provided by the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a gold geochemical element inversion method and a gold geochemical element inversion system, which are used for accurately acquiring gold geochemical element inversion information so as to accurately show gold geochemical elements on a remote sensing image under the condition of a complex geological background.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
As shown in fig. 1, the gold geochemical element inversion method provided by the invention comprises the following steps:
step 100: and acquiring gold geochemical element data, and determining screening data according to the gold geochemical element data.
Step 101: and acquiring GF-5 hyperspectral data of the region to be inverted.
Step 102: and determining mathematical transformation data according to the GF-5 hyperspectral data by adopting a mathematical transformation method.
Step 103: and carrying out correlation analysis on the screening data and the mathematical transformation data to obtain correlation data.
Step 104: and carrying out regression analysis on the correlation data to obtain a regression equation.
Step 105: and carrying out inversion operation according to the regression equation and the mathematical transformation data to obtain gold geochemical element inversion information.
In order to further improve the accuracy of inversion, the process of data screening in step 100 is specifically:
and step 1000, obtaining known gold geochemical element content data in the region to be inverted.
And 1001, screening the data of the content of the gold geochemical elements by adopting a hierarchical sampling method to obtain screened data. The step 1001 specifically includes:
and arranging the data of the gold geochemical element content according to the sequence of the metal element content from high to low to obtain the total number of the samples.
And layering the arranged gold geochemical element content data according to the total number of the samples.
And screening the data of the content of the gold geochemical elements in each layer according to a mode of sorting the content of the gold geochemical elements from high to low to obtain screened data.
Specifically, the gold element content is sorted from high to low, the total number of samples is N, the samples are divided into N layers by taking the N/N number order as a unit, and the number of samples in each layer (A)1,A2,A3,...,AiAnd (i-1, 2, 3.., n) screening data according to a method of ordering gold element content from high to low to obtain screened gold geochemical element data (M-A)1+A2+A3+...+Ai,i＝1,2,3,...,n)。
In order to provide a more accurate data base for the mathematical transformation, step 102 can be implemented in the following specific manner in the present invention:
and step 1020, carrying out radiometric calibration on the GF-5 hyperspectral data.
And step 1021, performing atmospheric correction on the GF-5 hyperspectral data subjected to radiometric calibration to obtain atmospheric correction data.
And step 1022, performing mathematical transformation on the atmosphere correction data to obtain mathematical transformation data.
Based on the specific implementation process of step 102, as shown in fig. 2, another implementation of the gold geochemical element inversion method provided by the invention is as follows:
step 200: collecting data of the known gold geochemical element content of the region.
Step 201: and screening the obtained gold element content data by using a hierarchical sampling method.
Step 202: and acquiring GF-5 hyperspectral data of the region.
Step 203: and carrying out radiometric calibration on the GF-5 hyperspectral data.
Step 204: and performing atmospheric correction on the data subjected to radiation calibration.
Step 205: the data after atmospheric correction is subjected to various mathematical transformations to obtain a special image form.
Step 206: and carrying out correlation analysis and regression analysis on the screened gold geochemical element data and the data subjected to mathematical transformation.
Step 207: and applying a formula obtained after regression analysis to the data subjected to mathematical transformation for inversion to obtain gold geochemical element inversion information.
The two embodiments of the gold geochemical element inversion method have the advantages that the representativeness of the sample is better, and the sampling error is smaller. Then, radiometric calibration and atmospheric correction are carried out on GF-5 hyperspectral data, so that the remote sensing image is closer to reality, mathematical transformation is carried out on the image after atmospheric correction, the relation between the image and gold geochemical element data can be effectively increased, and then correlation analysis and regression analysis are better carried out to obtain a better inversion effect.
The method preferably adopts first derivative transformation, second derivative transformation or continuum removal transformation as a mathematical transformation method to determine mathematical transformation data according to GF-5 hyperspectral data (namely, the data after atmospheric correction is subjected to various mathematical transformations to obtain a special image form). Each mathematical transformation method is implemented as follows:
and performing first derivative transformation on the data after the atmospheric correction to obtain first derivative transformation data.
And performing second derivative transformation on the data after the atmospheric correction to obtain second derivative transformation data.
And carrying out continuum removal transformation on the data after the atmospheric correction to obtain continuum removal transformation data.
The implementation processes of step 103, step 104 and step 206 are different based on the different mathematical transformation methods.
When the first derivative transformation method is adopted to determine the transformation data, correlation analysis needs to be carried out on the gold geochemical element data obtained by screening and the first derivative transformation data subsequently, 4 wave bands with the highest correlation coefficient are selected for carrying out regression analysis, and a formula obtained by the regression analysis is recorded as a first derivative regression equation.
When the second derivative transformation method is adopted to determine the transformation data, correlation analysis needs to be carried out on the gold geochemical element data obtained by screening and the second derivative transformation data subsequently, 4 wave bands with the highest correlation coefficient are selected for carrying out regression analysis, and a formula obtained by the regression analysis is recorded as a second derivative regression equation.
When the continuum removal transformation method is adopted to determine transformation data, correlation analysis needs to be carried out on the gold geochemical element data obtained by screening and the continuum removal transformation data subsequently, 4 wave bands with the highest correlation coefficient are selected for carrying out regression analysis, and a formula obtained by the regression analysis is recorded as a continuum removal regression equation.
Wherein, the invention preferably adopts bivariate correlation analysis as the correlation analysis method, which comprises the following specific steps:
the expression of the multiple linear regression analysis is:
y＝b0+b1·x1+b2·x2+b3·x3+...+bixi(i＝1，2，3，...n；bk(k ═ 1,2,3,.. i) is a regression coefficient)
In the formula, y is the content of gold element, b is a constant, and x is the pixel value of the transformation data.
Further, the regression equations obtained after performing the correlation analysis and the regression analysis are different because the mathematical transformation data used are different. Therefore, the inversion information obtained in the above steps 105 and 207 also corresponds to the mathematical transformation method used.
Specifically, the first derivative regression equation is substituted into the first derivative transformation data, and the obtained result is the first derivative gold geochemical element inversion information.
And substituting the second derivative regression equation into the second derivative transformation data to obtain a result, namely second derivative gold geochemical element inversion information.
And substituting the continuum removal regression equation into the continuum removal transformation data to obtain a result, namely the continuum removal gold geochemical element inversion information.
Based on the above description, a specific process for determining inversion information based on different mathematical transformation methods is shown in fig. 3.
In summary, in the above implementation process, the present invention includes collecting experimental geochemical element data, hierarchically sampling the collected gold geochemical element data, and performing a plurality of mathematical transformations, including a first derivative transformation, on the atmospheric corrected remote sensing image data, describing the rate of change of this function around this point. The second derivative transform, which represents the rate of change of the first derivative. From the figure, it reflects the concave-convex of the function image. Continuum removal transform, also known as envelope removal, geometrically, the envelope of a family of curves is a curve that is tangent to at least one point on each line of the family. And performing correlation analysis on the data after mathematical transformation and the extracted geochemical element data of the gold, wherein the correlation analysis refers to analyzing two or more variable elements with correlation, so that the correlation closeness degree of two variable factors is measured, and several wave bands with the highest correlation are selected for regression analysis according to the magnitude of a correlation coefficient, the regression analysis refers to a statistical analysis method for determining the quantitative relation of mutual dependence between two or more variables, and a formula after the regression analysis is applied to a remote sensing image for inversion to obtain the geochemical element inversion information.
The following describes a specific implementation procedure and technical effects achieved by the gold geochemical element inversion method provided by the present invention, with the implementation of the camellia village gold mine area in which geochemical element data have been acquired. In practical applications, the technical solution of the present invention is not limited thereto.
In the gold mining area of camellia village, the gold element data samples are as many as 20144, so how to perform data screening is of great importance. After analyzing the data composition, the content of the samples is sorted from high to low, the samples are screened by a stratified sampling method, 10 orders of magnitude are stratified, random sampling is carried out in each stratum, and finally 219 samples are selected for experiments, as shown in the following table 1.
Table 1 screening data
An image obtained by preprocessing GF-5 hyperspectral data is used as an original remote sensing image, and in view of the fact that good inversion results cannot be usually obtained by inverting the original image data, the original remote sensing image is subjected to transformation processing, including first derivative transformation, second derivative transformation and continuum removal transformation, as shown in figures 4 to 6, the transformed image is obtained, then gold geochemical element data obtained before is projected onto the transformed image to extract pixel values of different wave bands, the extracted result is combined with an actual sample to perform correlation analysis, and as shown in figure 7, the result is a result after the correlation analysis.
Table 2 first derivative correlation coefficient table
| Wave band | Correlation coefficient |
| B164 | -0.239 |
| B165 | -0.223 |
| B271 | 0.22 |
| B272 | 0.219 |
TABLE 3 second derivative correlation coefficient Table
| Wave band | Correlation coefficient |
| B165 | 0.203 |
| B272 | -0.199 |
| B101 | 0.192 |
| B163 | -0.187 |
Table 4 continuum decorrelation coefficient table
The highest 4 bands were selected for regression analysis, see tables 2-4, to obtain the following formula.
First derivative regression formula:
Y＝-2.478-0.232*XB164-0.268*XB165+0.47*XB288+0.565*XB289
second derivative regression formula:
Y＝14.738+1.663*XB101+0.151*XB163+0.413*XB165-1.472*XB289
continuum removal regression formula:
Y＝205.465-86.624*XB128-200.383*XB165-55.843*XB166+89.358*XB168
the above formulas are respectively substituted into the corresponding processed images to obtain the inversion results as shown in fig. 8 to 10, from which the position information, distribution information, and the like of the gold geochemical elements can be clearly seen.
To sum up, the gold geochemical element inversion method provided by the invention is used under the condition of complex geological background, and specifically comprises the following steps: sampling the collected experimental data by a hierarchical sampling method, then performing various mathematical transformations on the GF-5 hyperspectral data, including first derivative transformation, second derivative transformation and continuum removal transformation, extracting pixel values of the remote sensing image according to the screened data, performing correlation analysis on the screened data and the extracted data, measuring the degree of correlation closeness of two variable factors, selecting several wave bands with the highest correlation to perform regression analysis according to the magnitude of the correlation coefficient, applying a formula after the regression analysis to the remote sensing image to perform inversion to obtain geochemical element inversion information so as to express the geochemical element information on the remote sensing image, the large-range abnormal information is presumed according to the small-range geochemical data, the defect that field work is difficult to carry out in a complex area is overcome, and powerful scientific basis and support are provided for the prospecting of mines.
In addition, corresponding to the gold geochemical element inversion method provided above, the present invention also provides a gold geochemical element inversion system, as shown in fig. 11, which includes: the system comprises a screening module 1, a GF-5 hyperspectral data acquisition module 2, a mathematical transformation module 3, a correlation analysis module 4, a regression analysis module 5 and a gold geochemical element inversion module 6.
The screening module 1 is used for acquiring gold geochemical element data and determining screening data according to the gold geochemical element data.
The GF-5 hyperspectral data acquisition module 2 is used for acquiring GF-5 hyperspectral data of an area to be inverted.
The mathematical transformation module 3 is used for determining mathematical transformation data according to the GF-5 hyperspectral data by adopting a mathematical transformation method.
And the correlation analysis module 4 is used for performing correlation analysis on the screening data and the mathematical transformation data to obtain correlation data.
The regression analysis module 5 is used for performing regression analysis on the correlation data to obtain a regression equation.
And the gold geochemical element inversion module 6 is used for carrying out inversion operation according to the regression equation and the mathematical transformation data to obtain the gold geochemical element inversion information.
Preferably, the screening module 1 comprises: the device comprises an element content data acquisition unit and a screening unit.
The element content data acquisition unit is used for acquiring known gold geochemical element content data in the region to be inverted.
The screening unit is used for screening the data of the content of the gold geochemical elements by adopting a hierarchical sampling method to obtain screening data.
Wherein, above-mentioned screening unit still preferably includes: a ranking subunit, a layering subunit, and a screening subunit.
And the sorting subunit is used for sorting the gold geochemical element content data according to the sequence of the metal element content from high to low to obtain the total number of the samples.
And the layering subunit is used for layering the arranged gold geochemical element content data according to the total number of the samples.
And the screening subunit is used for screening the data of the gold geochemical element content of each layer according to the mode that the gold geochemical element content is sorted from high to low to obtain screening data.
Preferably, the mathematical transformation module 3 specifically includes: a radiometric calibration unit, a correction unit and a mathematical transformation unit.
And the radiometric calibration unit is used for radiometric calibration of the GF-5 hyperspectral data.
And the correction unit is used for performing atmospheric correction on the GF-5 hyperspectral data subjected to radiometric calibration to obtain correction data.
The mathematical transformation unit is used for performing mathematical transformation on the correction data to obtain mathematical transformation data.
Wherein, the mathematical transformation unit further preferably includes: a mathematical transformation subunit.
The mathematical transformation subunit is used for performing first derivative transformation, second derivative transformation or continuum removal transformation on the correction data to obtain mathematical transformation data. The mathematical transformation data is first derivative transformation data, second derivative transformation data, or continuum removal transformation data.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (10)
1. A gold geochemical element inversion method, comprising:
acquiring gold geochemical element data, and determining screening data according to the gold geochemical element data;
acquiring GF-5 hyperspectral data of an area to be inverted;
determining mathematical transformation data according to the GF-5 hyperspectral data by adopting a mathematical transformation method;
performing correlation analysis on the screening data and the mathematical transformation data to obtain correlation data;
carrying out regression analysis on the correlation data to obtain a regression equation;
and carrying out inversion operation according to the regression equation and the mathematical transformation data to obtain gold geochemical element inversion information.
2. The gold geochemical element inversion method of claim 1, wherein the obtaining gold geochemical element data and determining screening data from the gold geochemical element data includes:
acquiring known gold geochemical element content data in an area to be inverted;
and screening the data of the content of the gold geochemical elements by adopting a hierarchical sampling method to obtain the screened data.
3. The gold geochemical element retrieval method according to claim 2, wherein the screening of the gold geochemical element content data by the hierarchical sampling method to obtain the screening data specifically comprises:
arranging the gold geochemical element content data in a sequence from high metal element content to low metal element content to obtain the total number of samples;
layering the arranged geochemical element content data of the gold according to the total number of the samples;
and screening the data of the content of the gold geochemical elements in each layer according to the data screening mode that the content of the gold geochemical elements is sorted from high to low to obtain the screened data.
4. The gold geochemical element inversion method according to claim 1, wherein the determining of the mathematical transformation data from the GF-5 hyperspectral data by the mathematical transformation method specifically comprises:
performing radiometric calibration on the GF-5 hyperspectral data;
atmospheric correction is carried out on the GF-5 hyperspectral data after radiometric calibration, and correction data are obtained;
and performing mathematical transformation on the correction data to obtain the mathematical transformation data.
5. The gold geochemical element inversion method according to claim 4, wherein the mathematically transforming the correction data to obtain the mathematically transformed data specifically comprises:
performing first derivative transformation, second derivative transformation or continuum removal transformation on the correction data to obtain mathematical transformation data; the mathematical transformation data is first derivative transformation data, second derivative transformation data or continuum removal transformation data.
6. A gold geochemical element inversion system, comprising:
the screening module is used for acquiring gold geochemical element data and determining screening data according to the gold geochemical element data;
the GF-5 hyperspectral data acquisition module is used for acquiring GF-5 hyperspectral data of an area to be inverted;
the mathematical transformation module is used for determining mathematical transformation data according to the GF-5 hyperspectral data by adopting a mathematical transformation method;
the correlation analysis module is used for carrying out correlation analysis on the screening data and the mathematical transformation data to obtain correlation data;
the regression analysis module is used for carrying out regression analysis on the correlation data to obtain a regression equation;
and the gold geochemical element inversion module is used for carrying out inversion operation according to the regression equation and the mathematical transformation data to obtain gold geochemical element inversion information.
7. The gold geochemical element inversion system of claim 6 wherein the screening module includes:
the element content data acquisition unit is used for acquiring known gold geochemical element content data in the region to be inverted;
and the screening unit is used for screening the gold geochemical element content data by adopting a hierarchical sampling method to obtain the screening data.
8. The gold geochemical element inversion system of claim 7 wherein the screening unit includes:
the sorting subunit is used for sorting the gold geochemical element content data according to the sequence of the metal element content from high to low to obtain the total number of samples;
the layering subunit is used for layering the arranged gold geochemical element content data according to the total number of the samples;
and the screening subunit is used for screening the data of the gold geochemical element content of each layer according to the data of the gold geochemical element content in a mode of sequencing from high to low to obtain the screened data.
9. The gold geochemical element inversion system of claim 6 wherein the mathematical transformation module includes in particular:
the radiometric calibration unit is used for radiometric calibration of the GF-5 hyperspectral data;
the atmosphere correction unit is used for carrying out atmosphere correction on the GF-5 hyperspectral data subjected to radiometric calibration to obtain correction data;
and the mathematical transformation unit is used for performing mathematical transformation on the correction data to obtain the mathematical transformation data.
10. The gold geochemical element inversion system of claim 9 wherein the mathematical transformation unit includes:
the mathematical transformation subunit is used for carrying out first derivative transformation, second derivative transformation or continuum removal transformation on the correction data to obtain mathematical transformation data; the mathematical transformation data is first derivative transformation data, second derivative transformation data or continuum removal transformation data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110856570.0A CN113552075A (en) | 2021-07-28 | 2021-07-28 | Gold geochemical element inversion method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110856570.0A CN113552075A (en) | 2021-07-28 | 2021-07-28 | Gold geochemical element inversion method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113552075A true CN113552075A (en) | 2021-10-26 |
Family
ID=78104773
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110856570.0A Pending CN113552075A (en) | 2021-07-28 | 2021-07-28 | Gold geochemical element inversion method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113552075A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115101137A (en) * | 2022-06-28 | 2022-09-23 | 贵州逸宸未来科技有限公司 | Geochemical data analysis method based on deep learning |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108663330A (en) * | 2018-04-19 | 2018-10-16 | 中国国土资源航空物探遥感中心 | A kind of vegetation-covered area soil copper inversion method based on blade measured spectra |
| CN109508512A (en) * | 2018-12-28 | 2019-03-22 | 中南大学 | Change the earth's surface constituent content visited and combined with remote sensing and quantifies projectional technique and device |
-
2021
- 2021-07-28 CN CN202110856570.0A patent/CN113552075A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108663330A (en) * | 2018-04-19 | 2018-10-16 | 中国国土资源航空物探遥感中心 | A kind of vegetation-covered area soil copper inversion method based on blade measured spectra |
| CN109508512A (en) * | 2018-12-28 | 2019-03-22 | 中南大学 | Change the earth's surface constituent content visited and combined with remote sensing and quantifies projectional technique and device |
Non-Patent Citations (5)
| Title |
|---|
| XIA ZHANG 等: "INVERSION OF SOIL Cu CONCENTRATION BASED ON BAND SELECTION OF HYPERSPETRAL DATA", 《IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》 * |
| 丁海宁等: "黄土高原土壤铁元素含量遥感反演方法", 《遥感技术与应用》 * |
| 张霞等: "基于铁氧化物特征光谱和改进遗传算法反演土壤Pb含量", 《农业工程学报》 * |
| 王金凤等: "基于高光谱反射率的喀斯特地区土壤重金属锌元素含量反演", 《光谱学与光谱分析》 * |
| 马秀梅等: "高光谱遥感定量反演成矿元素含量的研究进展", 《地质找矿论丛》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115101137A (en) * | 2022-06-28 | 2022-09-23 | 贵州逸宸未来科技有限公司 | Geochemical data analysis method based on deep learning |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109614942B (en) | Forest disturbance long-time sequence monitoring method based on cloud computing platform | |
| Ren et al. | Estimating senesced biomass of desert steppe in Inner Mongolia using field spectrometric data | |
| CN108152235B (en) | Heavy metal content inversion method combining soil indoor and outdoor spectra | |
| CN109799195A (en) | A kind of high-precision fixed analysis method of laser induced breakdown spectroscopy | |
| Wang et al. | Effects of hyperspectral data with different spectral resolutions on the estimation of soil heavy metal content: From ground-based and airborne data to satellite-simulated data | |
| CN110793923A (en) | Hyperspectral soil data acquisition and analysis method based on mobile phone | |
| CN113591700A (en) | Prediction method for potential landslide and river blockage | |
| CN113436153A (en) | Method for predicting carbon components of undisturbed soil profile based on hyperspectral imaging and support vector machine technology | |
| Qu et al. | Estimating vegetation fraction using hyperspectral pixel unmixing method: a case study of a karst area in China | |
| CN113552075A (en) | Gold geochemical element inversion method and system | |
| Singh et al. | Multitemporal spectral analysis for cheatgrass (Bromus tectorum) classification | |
| CN114863294A (en) | Water quality parameter monitoring method and device suitable for brackish water lake | |
| Shamsoddini et al. | Mapping red edge-based vegetation health indicators using Landsat TM data for Australian native vegetation cover | |
| Liu et al. | Comparison of datasets obtained by global-scale geochemical sampling in Australia, China and Europe | |
| Ge et al. | A new perspective to near-infrared reflectance spectroscopy: A wavelet approach | |
| CN110596017B (en) | Hyperspectral image soil heavy metal concentration assessment method based on space weight constraint and variational self-coding feature extraction | |
| CN115825388A (en) | Training method, estimation method, device and equipment of heavy metal estimation model | |
| CN117058522A (en) | Remote sensing earth surface ecological change detection method integrating spectral neural network and leaf area index | |
| CN116773516A (en) | Soil carbon content analysis system based on remote sensing data | |
| Chabrillat et al. | Prediction of common surface soil properties using airborne and simulated EnMAP hyperspectral images: Impact of soil algorithm and sensor characteristic | |
| Li et al. | Comparison of two methods for estimation of soil water content from measured reflectance | |
| CN117129656B (en) | Screening and determining method and device for geochemical investigation index elements | |
| CN113624634B (en) | Method for estimating content of metal elements in buried environment | |
| Wilford et al. | National Geochemical Survey of Australia: comparison of geochemical and airborne radiometric data | |
| AU2021104180A4 (en) | Method and device for quantitatively estimating the content of surface elements by combining geochemical exploration and remote sensing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211026 |