CN103900964A - Hyperspectral image processing method used for extracting muscovite information - Google Patents
Hyperspectral image processing method used for extracting muscovite information Download PDFInfo
- Publication number
- CN103900964A CN103900964A CN201210572981.8A CN201210572981A CN103900964A CN 103900964 A CN103900964 A CN 103900964A CN 201210572981 A CN201210572981 A CN 201210572981A CN 103900964 A CN103900964 A CN 103900964A
- Authority
- CN
- China
- Prior art keywords
- wave band
- band image
- image
- pixel value
- deducts
- 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
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
The invention relates to a hyperspectral image processing method used for extracting muscovite information. The method includes preserving the wavebands having significant spectral signatures of muscovite, extracting images the wavebands of which are at 1355 nm, 1415 nm, 1445 nm, 1820 nm, 1940 nm, 2090 nm, 2120 nm, 2165 nm, 2195 nm, 2225 nm, 2270 nm, 2300 nm and 2345 nm, performing a series of determination and calculation, and calculating abundance values of the muscovite in different zones in the image scope. The method can remove the wavebands without significant signatures, thus highlighting the spectral signatures of the muscovite in an information extraction process, reducing influences of other ground features or noises, and reducing the amount of data to be processed. Through an IDL program, the method can achieve an objective of information extraction of final results by less manual operation, thus increasing the precision and the speed of muscovite information extraction.
Description
Technical field
The present invention relates to a kind of Technique for Hyper-spectral Images Classification for white mica information extraction, particularly relate to a kind of impact that reduces other atural objects or noise, reduced the Technique for Hyper-spectral Images Classification for white mica information extraction of the data volume of processing.
Background technology
The white mica information extracting method of current target in hyperspectral remotely sensed image is mainly spectrum all band coupling or the Spectral matching of partial continuous wave band, specific algorithm has spectrum angle, mixes demodulation filtering etc., because the material on earth's surface forms seldom by single mineral composition, these methods are subject to the impact of other ground-object spectrums or noise in the process of information extraction, information extraction precision is relatively low.Secondly existing spectrum extracting method manual steps is many, has increased artificial error in judgement.The 3rd is that high-spectral data wave band is many, and data volume is large, and the existing method processing time is long, has reduced speed and the application scale of data processing.Therefore, how in the process of white mica information extraction, to reduce impact, manual steps and the deal with data amount of other atural objects or noise, become one of forward position of current target in hyperspectral remotely sensed image processing.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of impact that reduces other atural objects or noise, has reduced the Technique for Hyper-spectral Images Classification for white mica information extraction of the data volume of processing.
For solving the problems of the technologies described above, a kind of Technique for Hyper-spectral Images Classification for white mica information extraction of the present invention, comprises successively:
The first step, obtains Hyperspectral imaging; To the pre-service of carrying out of Hyperspectral imaging, carry out atmospheric correction, obtain the image data of ground surface reflectance;
Second step, carries out resampling to the image data of ground surface reflectance, extracts wave band at 1355nm, 1415nm, 1445nm, 1820nm, 1940nm, 2090nm, 2120nm, 2165nm, 2195nm, 2225nm, 2270nm, 2300nm, the image of 2345nm;
The 3rd step, when the pixel value of 1355nm wave band image is greater than the corresponding pixel value of 1415nm wave band image, judges that this pixel satisfies condition one; Obtain the result image H1 that 1355nm wave band image deducts 1415nm wave band image;
The 4th step, when the pixel value of 1415nm wave band image is less than the corresponding pixel value of 1445nm wave band image, this pixel satisfies condition two; Obtain the result image H2 that 1445nm wave band image deducts 1415nm wave band image;
The 5th step, when the pixel value of 1820nm wave band image is greater than the corresponding pixel value of 1940nm wave band image, judges that this pixel satisfies condition three; Obtain the result image H3 that 1820nm wave band image deducts 1940nm wave band image;
The 6th step, when the pixel value of 2090nm wave band image is greater than the corresponding pixel value of 2120nm wave band image, judges that this pixel satisfies condition four; Obtain the result image H4 that 2090nm wave band image deducts 2120nm wave band image;
The 7th step, when the pixel value of 2165nm wave band image is greater than the corresponding pixel value of 2195nm wave band image, judges that this pixel satisfies condition five; Obtain the result image H5 that 2165nm wave band image deducts 2195nm wave band image;
The 8th step, when the pixel value of 2225nm wave band image is less than the corresponding pixel value of 2270nm wave band image, judges that this pixel satisfies condition six; Obtain the result image H6 that 2270nm wave band image deducts 2225nm wave band image;
The 9th step, when the pixel value of 2300nm wave band image is greater than the corresponding pixel value of 2345nm wave band image, judges that this pixel satisfies condition seven; Obtain the result image H7 that 2300nm wave band image deducts 2345nm wave band image;
The tenth step, selects the pixel to condition seven that simultaneously satisfies condition, and obtains H1, H2, H3, H4, H5, H6, H7 within the scope of above-mentioned pixel is added and H.
This method of the present invention has only been used 13 wave bands, and data volume has reduced 87%, and owing to being computing machine from moving a step extraction, has reduced the operation steps such as selection of principal component transform, end member wave spectrum, and arithmetic speed can improve more than 8 times.Owing to having removed most of wave band little to information extraction relation, reduce the interference to its spectrum of other materials or noise, improve the precision of information extraction.Rapid extraction to white mica information in airborne-remote sensing has good function and significance.
Embodiment
The present invention, by Hyperspectral imaging resampling, extracts specific band, carries out a series of judgements and calculating, calculates the Abundances of the zones of different white mica in image capturing range.
Specifically comprise successively:
The first step, obtains Hyperspectral imaging; To the pre-service of carrying out of Hyperspectral imaging, carry out atmospheric correction, obtain the image data of ground surface reflectance;
Second step, carries out resampling to the image data of ground surface reflectance, extracts wave band at 1355nm, 1415nm, 1445nm, 1820nm, 1940nm, 2090nm, 2120nm, 2165nm, 2195nm, 2225nm, 2270nm, 2300nm, the image of 2345nm;
The 3rd step, when the pixel value of 1355nm wave band image is greater than the corresponding pixel value of 1415nm wave band image, judges that this pixel satisfies condition one; Obtain the result image H1 that 1355nm wave band image deducts 1415nm wave band image;
The 4th step, when the pixel value of 1415nm wave band image is less than the corresponding pixel value of 1445nm wave band image, this pixel satisfies condition two; Obtain the result image H2 that 1445nm wave band image deducts 1415nm wave band image;
The 5th step, when the pixel value of 1820nm wave band image is greater than the corresponding pixel value of 1940nm wave band image, judges that this pixel satisfies condition three; Obtain the result image H3 that 1820nm wave band image deducts 1940nm wave band image;
The 6th step, when the pixel value of 2090nm wave band image is greater than the corresponding pixel value of 2120nm wave band image, judges that this pixel satisfies condition four; Obtain the result image H4 that 2090nm wave band image deducts 2120nm wave band image;
The 7th step, when the pixel value of 2165nm wave band image is greater than the corresponding pixel value of 2195nm wave band image, judges that this pixel satisfies condition five; Obtain the result image H5 that 2165nm wave band image deducts 2195nm wave band image;
The 8th step, when the pixel value of 2225nm wave band image is less than the corresponding pixel value of 2270nm wave band image, judges that this pixel satisfies condition six; Obtain the result image H6 that 2270nm wave band image deducts 2225nm wave band image;
The 9th step, when the pixel value of 2300nm wave band image is greater than the corresponding pixel value of 2345nm wave band image, judges that this pixel satisfies condition seven; Obtain the result image H7 that 2300nm wave band image deducts 2345nm wave band image;
The tenth step, selects the pixel to condition seven that simultaneously satisfies condition, and obtains H1, H2, H3, H4, H5, H6, H7 within the scope of above-mentioned pixel is added and H.
The value of H represents the Abundances of pixel white mica, and H value is larger, represents that the abundance of white mica in pixel is larger, and content is higher.
Claims (1)
1. for a Technique for Hyper-spectral Images Classification for white mica information extraction, comprise successively:
The first step, obtains Hyperspectral imaging; To the pre-service of carrying out of Hyperspectral imaging, carry out atmospheric correction, obtain the image data of ground surface reflectance;
Second step, carries out resampling to the image data of ground surface reflectance, extracts wave band at 1355nm, 1415nm, 1445nm, 1820nm, 1940nm, 2090nm, 2120nm, 2165nm, 2195nm, 2225nm, 2270nm, 2300nm, the image of 2345nm;
The 3rd step, when the pixel value of 1355nm wave band image is greater than the corresponding pixel value of 1415nm wave band image, judges that this pixel satisfies condition one; Obtain the result image H1 that 1355nm wave band image deducts 1415nm wave band image;
The 4th step, when the pixel value of 1415nm wave band image is less than the corresponding pixel value of 1445nm wave band image, this pixel satisfies condition two; Obtain the result image H2 that 1445nm wave band image deducts 1415nm wave band image;
The 5th step, when the pixel value of 1820nm wave band image is greater than the corresponding pixel value of 1940nm wave band image, judges that this pixel satisfies condition three; Obtain the result image H3 that 1820nm wave band image deducts 1940nm wave band image;
The 6th step, when the pixel value of 2090nm wave band image is greater than the corresponding pixel value of 2120nm wave band image, judges that this pixel satisfies condition four; Obtain the result image H4 that 2090nm wave band image deducts 2120nm wave band image;
The 7th step, when the pixel value of 2165nm wave band image is greater than the corresponding pixel value of 2195nm wave band image, judges that this pixel satisfies condition five; Obtain the result image H5 that 2165nm wave band image deducts 2195nm wave band image;
The 8th step, when the pixel value of 2225nm wave band image is less than the corresponding pixel value of 2270nm wave band image, judges that this pixel satisfies condition six; Obtain the result image H6 that 2270nm wave band image deducts 2225nm wave band image;
The 9th step, when the pixel value of 2300nm wave band image is greater than the corresponding pixel value of 2345nm wave band image, judges that this pixel satisfies condition seven; Obtain the result image H7 that 2300nm wave band image deducts 2345nm wave band image;
The tenth step, selects the pixel to condition seven that simultaneously satisfies condition, and obtains H1, H2, H3, H4, H5, H6, H7 within the scope of above-mentioned pixel is added and H.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210572981.8A CN103900964B (en) | 2012-12-25 | 2012-12-25 | A kind of Technique for Hyper-spectral Images Classification for muscovite information extraction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210572981.8A CN103900964B (en) | 2012-12-25 | 2012-12-25 | A kind of Technique for Hyper-spectral Images Classification for muscovite information extraction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103900964A true CN103900964A (en) | 2014-07-02 |
| CN103900964B CN103900964B (en) | 2016-10-05 |
Family
ID=50992425
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210572981.8A Active CN103900964B (en) | 2012-12-25 | 2012-12-25 | A kind of Technique for Hyper-spectral Images Classification for muscovite information extraction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103900964B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105787916A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | Hyperspectral image processing method for extracting information of dickite |
| CN105787914A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | Hyperspectral image processing method for alunite information extraction |
| CN105787483A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | Hyperspectral image processing method for extracting information of lizardite |
| CN105784602A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | High spectral image processing method for extracting information of rhodochrosite |
| CN105787915A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | Hyperspectral image processing method for extracting information of jarosite |
| CN106845326A (en) * | 2015-12-04 | 2017-06-13 | 核工业北京地质研究院 | A kind of glacier recognition methods based on Airborne Hyperspectral remotely-sensed data |
| CN108007902A (en) * | 2016-10-27 | 2018-05-08 | 核工业北京地质研究院 | A kind of method that muscovite Al-OH absorptions position is calculated with high-spectral data |
| CN109596535A (en) * | 2018-12-20 | 2019-04-09 | 核工业北京地质研究院 | A kind of spectral manipulation method suitable for extracting ferric ion abundance messages |
| CN109738371A (en) * | 2018-12-20 | 2019-05-10 | 核工业北京地质研究院 | A kind of spectral manipulation method suitable for extracting ferrous ion abundance messages |
| CN111044464A (en) * | 2019-12-19 | 2020-04-21 | 核工业北京地质研究院 | Data processing method suitable for extracting vegetation abundance information |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102609711A (en) * | 2012-02-21 | 2012-07-25 | 核工业北京地质研究院 | Information extraction method applicable to hyperspectral image |
| CN102636778A (en) * | 2012-02-21 | 2012-08-15 | 核工业北京地质研究院 | Information extracting method suitable for high-spectrum image |
-
2012
- 2012-12-25 CN CN201210572981.8A patent/CN103900964B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102609711A (en) * | 2012-02-21 | 2012-07-25 | 核工业北京地质研究院 | Information extraction method applicable to hyperspectral image |
| CN102636778A (en) * | 2012-02-21 | 2012-08-15 | 核工业北京地质研究院 | Information extracting method suitable for high-spectrum image |
Non-Patent Citations (5)
| Title |
|---|
| KUMAR S.ET AL.: "Best-bases feature extraction algorithms for classfication of hyperspectral data", 《GEOSCIENCE AND REMOTE SENSING》, vol. 39, no. 7, 31 July 2001 (2001-07-31), pages 1368 - 1379, XP011021771 * |
| WINTER M E，ET AL.: "comparison of Approaches for determining End-members in hyperspectral data", 《IEEE AEROSPACE CONFERENCE PROCEEDINGS》》, 31 December 2000 (2000-12-31), pages 305 - 313 * |
| 杨燕杰，赵英俊: "航空成像光谱的蚀变信息提取技术", 《科技导报》, vol. 29, no. 23, 31 December 2011 (2011-12-31), pages 57 - 61 * |
| 杨燕杰，赵英俊: "高光谱在油气勘探中的国内外研究现状", 《科学技术与工程》, vol. 11, no. 6, 28 February 2011 (2011-02-28), pages 1671 - 1815 * |
| 杨燕杰等: "富二价铁岩石信息提取及在铀矿勘查中的应用", 《测绘科学》, vol. 36, no. 5, 30 September 2011 (2011-09-30), pages 77 - 78 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105787916A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | Hyperspectral image processing method for extracting information of dickite |
| CN105787914A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | Hyperspectral image processing method for alunite information extraction |
| CN105787483A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | Hyperspectral image processing method for extracting information of lizardite |
| CN105784602A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | High spectral image processing method for extracting information of rhodochrosite |
| CN105787915A (en) * | 2014-12-22 | 2016-07-20 | 核工业北京地质研究院 | Hyperspectral image processing method for extracting information of jarosite |
| CN106845326A (en) * | 2015-12-04 | 2017-06-13 | 核工业北京地质研究院 | A kind of glacier recognition methods based on Airborne Hyperspectral remotely-sensed data |
| CN106845326B (en) * | 2015-12-04 | 2020-10-23 | 核工业北京地质研究院 | Glacier identification method based on aviation hyperspectral remote sensing data |
| CN108007902A (en) * | 2016-10-27 | 2018-05-08 | 核工业北京地质研究院 | A kind of method that muscovite Al-OH absorptions position is calculated with high-spectral data |
| CN108007902B (en) * | 2016-10-27 | 2020-06-19 | 核工业北京地质研究院 | Method for calculating Al-OH absorption position of muscovite by using hyperspectral data |
| CN109596535A (en) * | 2018-12-20 | 2019-04-09 | 核工业北京地质研究院 | A kind of spectral manipulation method suitable for extracting ferric ion abundance messages |
| CN109738371A (en) * | 2018-12-20 | 2019-05-10 | 核工业北京地质研究院 | A kind of spectral manipulation method suitable for extracting ferrous ion abundance messages |
| CN111044464A (en) * | 2019-12-19 | 2020-04-21 | 核工业北京地质研究院 | Data processing method suitable for extracting vegetation abundance information |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103900964B (en) | 2016-10-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103900964A (en) | Hyperspectral image processing method used for extracting muscovite information | |
| CN103900965A (en) | Hyperspectral image processing method used for extracting calcite information | |
| CN103903225A (en) | Hyperspectral image processing method for dolomite information extraction | |
| CN105787915A (en) | Hyperspectral image processing method for extracting information of jarosite | |
| CN107220988B (en) | Part image edge extraction method based on improved canny operator | |
| CN105891230B (en) | Fruit appearance detection method based on spectral image analysis | |
| CN106845326B (en) | Glacier identification method based on aviation hyperspectral remote sensing data | |
| CN108444954B (en) | Spectral signal peak detection method, device and system | |
| CN102636778A (en) | Information extracting method suitable for high-spectrum image | |
| CN103902998A (en) | High-spectral image processing method for chlorite information extraction | |
| Demir et al. | Hyperspectral image classification using denoising of intrinsic mode functions | |
| CN103901497A (en) | High-spectral image processing method for illite information extraction | |
| CN103900967B (en) | A kind of Technique for Hyper-spectral Images Classification for kaolin information extraction | |
| Qu et al. | Dimensionality reduction and derivative spectral feature optimization for hyperspectral target recognition | |
| CN103900966B (en) | A kind of Technique for Hyper-spectral Images Classification for allochite information extraction | |
| CN103902999B (en) | A kind of Technique for Hyper-spectral Images Classification for montmorillonite information extraction | |
| CN105784602A (en) | High spectral image processing method for extracting information of rhodochrosite | |
| CN105405102A (en) | High-spectral image processing method for gibbsite information extraction | |
| CN105787483A (en) | Hyperspectral image processing method for extracting information of lizardite | |
| CN104574283B (en) | A kind of Technique for Hyper-spectral Images Classification for pyrophyllite information extraction | |
| CN104573690B (en) | A kind of Technique for Hyper-spectral Images Classification for gypsum information extraction | |
| CN104732488B (en) | A kind of Technique for Hyper-spectral Images Classification for actinolite information extraction | |
| CN105787914A (en) | Hyperspectral image processing method for alunite information extraction | |
| CN105787916A (en) | Hyperspectral image processing method for extracting information of dickite | |
| CN102819750A (en) | Hyperspectral image classification method based on F-EEMD (fast ensemble empirical mode decomposition) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |