CN104990880A - Sericite mineral relative abundance calculating method based on aviation hyperspectral remote sensing data - Google Patents
Sericite mineral relative abundance calculating method based on aviation hyperspectral remote sensing data Download PDFInfo
- Publication number
- CN104990880A CN104990880A CN201510399242.7A CN201510399242A CN104990880A CN 104990880 A CN104990880 A CN 104990880A CN 201510399242 A CN201510399242 A CN 201510399242A CN 104990880 A CN104990880 A CN 104990880A
- Authority
- CN
- China
- Prior art keywords
- wave band
- sericite
- relative abundance
- absorption position
- sericite mineral
- 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.)
- Granted
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention belongs to the technical field of remote sensing geology survey, and particularly relates to a sericite mineral relative abundance calculating method based on aviation hyperspectral remote sensing data. The method comprises the following steps that data preprocessing is carried out on the aviation hyperspectral remote sensing data; the spectral absorption position of a sericite mineral image is calculated; a mis-extraction area is removed; the absorption depth corresponding to the spectral absorption position of the sericite mineral image is calculated; a sericite mineral relative abundance graph is manufactured. By means of the method, the absorption position of the sericite mineral, the absorption depth and the relative abundance can be accurately calculated.
Description
Technical field
The invention belongs to Remote Sensing Geological Investigation technical field, be specifically related to a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data.
Background technology
High spectrum resolution remote sensing technique is the cutting edge technology of remote sensing fields, can obtain the continuous and complete spectrum of atural object, directly carry out Objects recognition according to curve of spectrum feature.In remote-sensing geology application, adopt high spectrum resolution remote sensing technique can on the basis distinguishing mineral large classes Direct Recognition Within Monominerals.
The Airborne Hyperspectral remote-sensing mineral of prior art extracts the identification that active side overweights mineral species, many employing visual identification methods in the detailed division of dissimilar sericite mineral, be difficult to quantitatively calculate absorption position, and then be difficult to the content of mineral on earth's surface, i.e. relative abundance, accurately calculates.
Summary of the invention
The technical issues that need to address of the present invention are: prior art is difficult to accurately calculate mineral absorption position, relative abundance.
Technical scheme of the present invention is as described below:
Based on sericite mineral relative abundance computing method for Airborne Hyperspectral remotely-sensed data, comprise the following steps: step 1 pair Airborne Hyperspectral remotely-sensed data carries out data prediction; Step 2 calculates sericite mineral image spectral absorption position; Step 3 is removed and is extracted region by mistake; Step 4 calculates the absorption degree of depth corresponding to sericite mineral image spectral absorption position; Step 5 makes sericite mineral relative abundance figure.
Preferably:
In step 1, systems radiate correction and geometry correction are carried out to the Airborne Hyperspectral remotely-sensed data obtained, obtains the spectral radiance data with geographic coordinate; Spectral radiance data, through atmospheric correction and rebuilding spectrum, obtain hyper spectral reflectance data.
Preferably:
In step 2, application hyper spectral reflectance data, according to sericite mica ore thing spectrum characteristic selection high-spectral data characteristic wave bands, cut by the fitting a straight line of computational reflect rate minimum both sides and hand over position, obtain the reflectivity of the sericite mineral position of image spectral absorption accurately and correspondence thereof;
The fitting a straight line of described reflectivity minimum both sides is cut and is handed over the computing method of position as described below:
According to sericite mica ore object light spectrum signature, the wave band b that preliminary selected image spectral absorption position may be in
n; Wavelength is selected to be less than wave band b
n, and with wave band b
nthe most adjacent wave band b
n-1with secondary adjacent wave band b
n-2, select wavelength to be greater than wave band b
n, and and b
nthe wave band b that wave band is the most adjacent
n+1with secondary adjacent wave band b
n+2, adopt following formula to calculate the reflectivity Bx of image spectral absorption position bx and correspondence thereof:
In formula,
B
n-2, b
n-1, b
n+1, b
n+2be followed successively by wave band b
n-2, wave band b
n-1, wave band b
n+1, wave band b
n+2corresponding wavelength;
B
n-2, B
n-1, B
n+1, B
n+2be followed successively by wave band b
n-2, wave band b
n-1, wave band b
n+1, wave band b
n+2corresponding reflectivity.
Preferably:
In step 3, for the image spectral absorption position bx that step 2 calculates, if bx numerical value had both been greater than wave band b
n-1corresponding wavelength b
n-1be less than again wave band b
n+1corresponding wavelength b
n+1, then show that the absorption position bx that step 2 calculates is reasonable value, proceed step 4; Otherwise the absorption position bx that step 2 calculates being rejected as extracting region by mistake, returning the step 2 selected image spectral absorption position wave band that may be in again after rejecting, recalculating sericite mineral image spectral absorption position.
Preferably:
In step 4, according to sericite mica ore object light spectral curve configuration, set left point of shoulder and right point of shoulder, adopt following formula to calculate the absorption degree of depth corresponding to image spectral absorption position:
In formula,
H is the absorption degree of depth corresponding to image spectral absorption position;
B
l, b
rbe followed successively by left point of shoulder, wavelength that right point of shoulder is corresponding;
B
l, B
rbe followed successively by left point of shoulder, reflectivity that right point of shoulder is corresponding.
Preferably:
In step 5, the absorption depth H corresponding according to image spectral absorption position, carries out Threshold segmentation, adopts the color of the different depth to represent relative abundance height, makes sericite mineral relative abundance figure.
Beneficial effect of the present invention is:
(1) a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data of the present invention, the hyper spectral reflectance data that reflectivity values scope is 0 ~ 1 are obtained, for sericite Minerals identification provides remotely-sensed data by the pre-service of Airborne Hyperspectral remotely-sensed data;
(2) a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data of the present invention, Airborne Hyperspectral Remote Sensing Reflectance data decimation characteristic wave bands carries out equations simultaneousness, calculate sericite mineral and absorb position accurately, improve the accuracy absorbing depth calculation;
(3) a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data of the present invention, Airborne Hyperspectral reflectivity data band math obtains the absorption degree of depth, for the accurate calculating of sericite abundance provides foundation.
Accompanying drawing explanation
Fig. 1 is sericite mineral image spectral absorption position view;
Fig. 2 accurately absorbs position view after the matching of sericite mineral spectra;
Fig. 3 is sericite mineral absorption depth calculation schematic diagram.
Embodiment
Below in conjunction with drawings and Examples, a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data of the present invention are described in detail.
A kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data of the present invention, comprise the following steps:
Step 1 pair Airborne Hyperspectral remotely-sensed data carries out data prediction
Systems radiate correction and geometry correction are carried out to the Airborne Hyperspectral remotely-sensed data obtained, obtains the spectral radiance data with geographic coordinate.Spectral radiance data, through atmospheric correction and rebuilding spectrum, obtain the hyper spectral reflectance data of reflectivity values in 0 ~ 1 scope.
Described systems radiate correction, geometry correction, atmospheric correction and rebuilding spectrum are known to the skilled person general knowledge.
Step 2 calculates sericite mineral image spectral absorption position
Application hyper spectral reflectance data, according to sericite mica ore thing spectrum characteristic selection high-spectral data characteristic wave bands, cut by the fitting a straight line of computational reflect rate minimum both sides and hand over position, obtain the reflectivity of the sericite mineral position of image spectral absorption accurately and correspondence thereof.
The fitting a straight line of described reflectivity minimum both sides is cut and is handed over the computing method of position as described below:
According to sericite mica ore object light spectrum signature, the wave band b that preliminary selected image spectral absorption position may be in
n; Wavelength is selected to be less than wave band b
n, and with wave band b
nthe most adjacent wave band b
n-1with secondary adjacent wave band b
n-2, select wavelength to be greater than wave band b
n, and with wave band b
nthe most adjacent wave band b
n+1with secondary adjacent wave band b
n+2, adopt following formula to calculate the reflectivity Bx of image spectral absorption position bx and correspondence thereof:
In formula,
B
n-2, b
n-1, b
n+1, b
n+2be followed successively by wave band b
n-2, wave band b
n-1, wave band b
n+1, wave band b
n+2corresponding wavelength;
B
n-2, B
n-1, B
n+1, B
n+2be followed successively by wave band b
n-2, wave band b
n-1, wave band b
n+1, wave band b
n+2corresponding reflectivity.
In the present embodiment, according to sericite mica ore object light spectrum signature, select high-spectral data characteristic wave bands shown in following table.
Sequence number | Wave band | Wave band corresponding wavelength | Reflectivity |
1 | b1 | 2165nm | B1 |
2 | b2 | 2180nm | B2 |
3 | b3 | 2195nm | B3 |
4 | b4 | 2210nm | B4 |
5 | b5 | 2225nm | B5 |
6 | b6 | 2240nm | B6 |
7 | b7 | 2255nm | B7 |
According to sericite standard mineral spectrum, its image spectrum curve may be arranged in of these three positions of 2195nm, 2210nm and 2225nm in 2165nm ~ 2255nm scope internal reflection rate minimum value.For these three kinds possibilities, carry out following calculating sericite mineral spectral absorption position bx accurately respectively, and the reflectivity Bx of correspondence:
(1) situation of 2195nm is positioned at for image spectrum reflectance minimum
Select b1, b2 wave band adjacent forward with the b3 wave band residing for 2195nm, and b4, b5 wave band adjacent backward, adopt following formula to calculate the reflectivity Bx of spectral absorption position bx and correspondence thereof:
Obtain:
(2) situation of 2210nm is positioned at for image spectrum reflectance minimum
Select b2, b3 wave band adjacent forward with the b4 wave band residing for 2210nm, and b5, b6 wave band adjacent backward, adopt following formula to calculate the reflectivity Bx of spectral absorption position bx and correspondence thereof:
Obtain:
(3) situation of 2225nm is positioned at for image spectrum reflectance minimum
Select b3, b4 wave band adjacent forward with the b5 wave band residing for 2225nm, and b6, b7 wave band adjacent backward, adopt following formula to calculate the reflectivity Bx of spectral absorption position bx and correspondence thereof:
Obtain:
Step 3 is removed and is extracted region by mistake
Adopt band math, ensure that the sericite calculated accurately absorbs position and drops in corresponding wave band interval range, remove the region of by mistake extracting.
Specifically, for the absorption position bx that step 2 calculates, if bx numerical value had both been greater than wave band b
n-1corresponding wavelength b
n-1be less than again wave band b
n+1corresponding wavelength b
n+1, then show that the absorption position bx that step 2 calculates is reasonable value, proceed step 4; Otherwise the absorption position bx that step 2 calculates being rejected as extracting region by mistake, returning the step 2 selected image spectral absorption position wave band that may be in again after rejecting, recalculating sericite mineral image spectral absorption position.
Step 4 calculates the absorption degree of depth corresponding to sericite mineral image spectral absorption position
According to sericite mica ore object light spectral curve configuration, set left point of shoulder and right point of shoulder, left point of shoulder and wavelength corresponding to right point of shoulder are respectively b
land b
r, left point of shoulder and reflectivity corresponding to right point of shoulder are respectively B
land B
r, adopt following formula to calculate the absorption degree of depth corresponding to image spectral absorption position:
In formula,
H is the absorption degree of depth corresponding to image spectral absorption position.
In the present embodiment, left point of shoulder and wavelength corresponding to right point of shoulder are chosen according to data processing experience.Because the wave band of Airborne Hyperspectral data is spaced apart 15nm, the reflection shoulder of absorption peak both sides is often positioned at 2135nm and 2285nm place, experience shows, no matter absorbs position and is positioned at 2195nm, 2210nm or 2225nm, selects these two shoulder wavelength can measure sericite well and absorbs the degree of depth.Therefore, in the present embodiment, left point of shoulder and wavelength corresponding to right point of shoulder are respectively 2135nm and 2285nm.
Step 5 makes sericite mineral relative abundance figure
The absorption depth H corresponding according to image spectral absorption position, carries out Threshold segmentation, adopts the color of the different depth to represent relative abundance height, makes sericite mineral relative abundance figure, for geological analysis provides remote sensing criterion.
Described Threshold segmentation is known to the skilled person general knowledge.
Claims (6)
1., based on sericite mineral relative abundance computing method for Airborne Hyperspectral remotely-sensed data, it is characterized in that: comprise the following steps:
Step 1 pair Airborne Hyperspectral remotely-sensed data carries out data prediction;
Step 2 calculates sericite mineral image spectral absorption position;
Step 3 is removed and is extracted region by mistake;
Step 4 calculates the absorption degree of depth corresponding to sericite mineral image spectral absorption position;
Step 5 makes sericite mineral relative abundance figure.
2. a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data according to claim 1, is characterized in that:
In step 1, systems radiate correction and geometry correction are carried out to the Airborne Hyperspectral remotely-sensed data obtained, obtains the spectral radiance data with geographic coordinate; Spectral radiance data, through atmospheric correction and rebuilding spectrum, obtain hyper spectral reflectance data.
3. a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data according to claim 2, is characterized in that:
In step 2, application hyper spectral reflectance data, according to sericite mica ore thing spectrum characteristic selection high-spectral data characteristic wave bands, cut by the fitting a straight line of computational reflect rate minimum both sides and hand over position, obtain the reflectivity of the sericite mineral position of image spectral absorption accurately and correspondence thereof;
The fitting a straight line of described reflectivity minimum both sides is cut and is handed over the computing method of position as described below:
According to sericite mica ore object light spectrum signature, the wave band b that preliminary selected image spectral absorption position may be in
n; Wavelength is selected to be less than wave band b
n, and with wave band b
nthe most adjacent wave band b
n-1with secondary adjacent wave band b
n-2, select wavelength to be greater than wave band b
n, and and b
nthe wave band b that wave band is the most adjacent
n+1with secondary adjacent wave band b
n+2, adopt following formula to calculate the reflectivity Bx of image spectral absorption position bx and correspondence thereof:
In formula,
B
n-2, b
n-1, b
n+1, b
n+2be followed successively by wave band b
n-2, wave band b
n-1, wave band b
n+1, wave band b
n+2corresponding wavelength;
B
n-2, B
n-1, B
n+1, B
n+2be followed successively by wave band b
n-2, wave band b
n-1, wave band b
n+1, wave band b
n+2corresponding reflectivity.
4. a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data according to claim 3, is characterized in that:
In step 3, for the image spectral absorption position bx that step 2 calculates, if bx numerical value had both been greater than wave band b
n-1corresponding wavelength b
n-1be less than again wave band b
n+1corresponding wavelength b
n+1, then show that the absorption position bx that step 2 calculates is reasonable value, proceed step 4; Otherwise the absorption position bx that step 2 calculates being rejected as extracting region by mistake, returning the step 2 selected image spectral absorption position wave band that may be in again after rejecting, recalculating sericite mineral image spectral absorption position.
5. a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data according to claim 4, is characterized in that:
In step 4, according to sericite mica ore object light spectral curve configuration, set left point of shoulder and right point of shoulder, adopt following formula to calculate the absorption degree of depth corresponding to image spectral absorption position:
In formula,
H is the absorption degree of depth corresponding to image spectral absorption position;
B
l, b
rbe followed successively by left point of shoulder, wavelength that right point of shoulder is corresponding;
B
l, B
rbe followed successively by left point of shoulder, reflectivity that right point of shoulder is corresponding.
6. a kind of sericite mineral relative abundance computing method based on Airborne Hyperspectral remotely-sensed data according to claim 5, is characterized in that:
In step 5, the absorption depth H corresponding according to image spectral absorption position, carries out Threshold segmentation, adopts the color of the different depth to represent relative abundance height, makes sericite mineral relative abundance figure.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510399242.7A CN104990880B (en) | 2015-07-09 | 2015-07-09 | Sericite mineral relative abundance computational methods based on Airborne Hyperspectral remotely-sensed data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510399242.7A CN104990880B (en) | 2015-07-09 | 2015-07-09 | Sericite mineral relative abundance computational methods based on Airborne Hyperspectral remotely-sensed data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104990880A true CN104990880A (en) | 2015-10-21 |
CN104990880B CN104990880B (en) | 2017-10-27 |
Family
ID=54302721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510399242.7A Active CN104990880B (en) | 2015-07-09 | 2015-07-09 | Sericite mineral relative abundance computational methods based on Airborne Hyperspectral remotely-sensed data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104990880B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107389571A (en) * | 2017-07-21 | 2017-11-24 | 中国科学院遥感与数字地球研究所 | A kind of fine discrimination method of mineral constituent high-spectrum remote-sensing |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101871884A (en) * | 2010-06-02 | 2010-10-27 | 中国国土资源航空物探遥感中心 | Atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data |
CN103175801A (en) * | 2013-03-14 | 2013-06-26 | 中国国土资源航空物探遥感中心 | Large-batch automatic hyperspectral remote sensing mineral mapping method |
-
2015
- 2015-07-09 CN CN201510399242.7A patent/CN104990880B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101871884A (en) * | 2010-06-02 | 2010-10-27 | 中国国土资源航空物探遥感中心 | Atmospheric correction and regional mineral map spotting method utilizing multi-scene ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) remote sensing data |
CN103175801A (en) * | 2013-03-14 | 2013-06-26 | 中国国土资源航空物探遥感中心 | Large-batch automatic hyperspectral remote sensing mineral mapping method |
Non-Patent Citations (6)
Title |
---|
GILA NOTESCO ET AL.: "Mineral Classification of Land Surface Using Multispectral LWIR and Hyperspectral SWIR Remote-Sensing Data.A Case Study over the Sokolov Lignite Open-Pit Mines,the Czech Republic", 《REMOTE SENSING》 * |
HONGYUAN HUO ET AL.: "Mineral Mapping and Ore Prospecting with HyMap Data over Eastern Tien Shan,Xinjiang Uyghur Autonomous Regin", 《REMOTE SENSING》 * |
孙雨 等: "相山铀矿岩芯HySpex成像高光谱数据蚀变矿物提取及其地质意义", 《地质与勘探》 * |
孙雨 等: "青海省都兰县阿斯哈金矿区HySpex高光谱矿物填图及其找矿意义", 《地质学报》 * |
杨燕杰 等: "航空成像光谱的蚀变信息提取技术", 《科技导报》 * |
杨燕杰 等: "高光谱影像预处理技术", 《科技导报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107389571A (en) * | 2017-07-21 | 2017-11-24 | 中国科学院遥感与数字地球研究所 | A kind of fine discrimination method of mineral constituent high-spectrum remote-sensing |
CN107389571B (en) * | 2017-07-21 | 2019-06-14 | 中国科学院遥感与数字地球研究所 | A kind of fine discrimination method of mineral constituent high-spectrum remote-sensing |
Also Published As
Publication number | Publication date |
---|---|
CN104990880B (en) | 2017-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104050252B (en) | Hyperspectral remote sensing alteration information extracting method | |
CN103700092B (en) | A kind of forest brulee extraction method based on sequential remote sensing image | |
CN105526874B (en) | A kind of oil film thickness recognition methods based on spectral signature parameter | |
CN103175801B (en) | Large-batch automatic hyperspectral remote sensing mineral mapping method | |
CN109580497B (en) | Hyperspectral mineral abnormal information extraction method based on singularity theory | |
CN104867139B (en) | A kind of remote sensing image clouds and shadow detection method based on radiation field | |
CN107358250A (en) | Body gait recognition methods and system based on the fusion of two waveband radar micro-doppler | |
CN104331880A (en) | Hyper-spectral mixed pixel decomposition method based on geometric spatial spectral structure information | |
CN110188716B (en) | Mineralization information extraction method and system | |
CN103984940A (en) | Method for identifying hematitization based on hyperspectral data | |
WO2024083261A1 (en) | Image spectrum technology-based geological sketching method and system for tunnel face | |
CN103208011A (en) | Hyperspectral image space-spectral domain classification method based on mean value drifting and group sparse coding | |
CN107421892A (en) | A kind of hyperspectral data processing method for water body information | |
Potin et al. | An abrupt change detection algorithm for buried landmines localization | |
CN105528580A (en) | Hyperspectral curve matching method based on absorption peak characteristic | |
CN111709280B (en) | Rapid identification method for hyperspectral gold mineralization alteration scenic area | |
CN103926203B (en) | A kind of for the probabilistic Spectral Angle Mapping method of object spectrum | |
CN101916439B (en) | Space debris hyperspectral sequence detection method based on Hilbert-Huang transform | |
Li et al. | Spectral–spatial depth-based framework for hyperspectral underwater target detection | |
CN104990880A (en) | Sericite mineral relative abundance calculating method based on aviation hyperspectral remote sensing data | |
CN104268896A (en) | Hyper spectrum dimensionality reduction matching method and system based on spectrum sampling histogram | |
CN105259127B (en) | Vegetation index calculation method based on comprehensive absorbability | |
CN104076003A (en) | Extraction method of mineral spectrum absorption characteristic parameters | |
CN103645517A (en) | Comprehensive anomaly extraction method based on blind source separation technology and apparatus thereof | |
CN103927438B (en) | Near-infrared wavelength Variable Selection based on successive projection algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |