CN103207415A - Method for extracting mineralization-alteration information - Google Patents
Method for extracting mineralization-alteration information Download PDFInfo
- Publication number
- CN103207415A CN103207415A CN2012105868558A CN201210586855A CN103207415A CN 103207415 A CN103207415 A CN 103207415A CN 2012105868558 A CN2012105868558 A CN 2012105868558A CN 201210586855 A CN201210586855 A CN 201210586855A CN 103207415 A CN103207415 A CN 103207415A
- Authority
- CN
- China
- Prior art keywords
- pixel
- component
- alteration
- remote sensing
- image
- 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
Abstract
The invention discloses a method for extracting mineralization-alteration information. The method comprises the steps of effective background pixel identifications, selection of local variable windows, in-window main component analysis, construction of abnormal variables, abnormal recognition and assignment and alteration abnormal output. Interference information such as water, cloud, ice and snow is eliminated through identifications and statistics of effective background pixels; and statistics variables are limited within the windows rather than the whole image range during main component analysis through local variable windows, therefore further reduction of interference of ambient environment noise is facilitated, and the extraction capacity of alteration information is improved. By means of the method, the potential of mineralization-alteration charting of thematic mapper (TM) and advanced spaceborne thermal emission and reflection radiometer (ASTER) remotely-sensed can be further mined, and the indication capacity in geological exploration can be improved.
Description
Technical field
The present invention relates to a kind of extracting method of Extract Mineralized Alteration Information, belong to the geologic prospect field.
Background technology
The mineral deposit is in the overall process that generates, develops, and its minerogentic condition shows as the geologic anomaly event in the geology evolutionary process, thereby the searching geologic anomaly becomes one of the necessary condition in prediction mineral deposit.Geochemical anomaly, geophysical anomaly and Extract Mineralized Alteration remote sensing abnormal have become indicator for deposit parameter commonly used in the exploration prospecting.The principal component analysis (PCA) of feature guiding is that the Crosta technology has been brought into play the important ore deposit indicative function of looking for, the information of mineral deposit alterated rocks strengthen handle and mineral exploration in obtained application widely.But traditional Crosta technology is based on the view picture image carries out the proper vector value that each component is calculated in statistical study, cause having comprised in the statistical computation process interfere informations such as water body, cloud, ice and snow, thereby make resulting abnormal results also comprise much noise, influenced the precision that Extract Mineralized Alteration Information is extracted to a certain extent.
Summary of the invention
The object of the invention is to provide a kind of method that Extract Mineralized Alteration Information is extracted precision that improves.
Technical scheme of the present invention may further comprise the steps:
A. import a width of cloth remote sensing image A;
B. the remote sensing image A of input is carried out differentiation and the sign of effective background pixel, wherein effectively the background pixel refers to be non-corrupt data pixel, non-water body pixel, non-cloud pixel, non-ice and snow pixel and non-vegetation pixel;
C. the remote sensing image that carries out effective background pixel sign through the b step active window with N * N pixel is scanned and adds up each pixel, wherein N is the odd number between 11 to 31, this statistics is centered by the pixel of analyzing spot place, search for since 11 * 11 windows, enlarge search window successively, largest extension to 31 * 31 windows, the ratio that accounts for whole window pixel quantity up to the quantity of effective background pixel reaches 51%;
D. the remote sensing image that carries out effective background pixel sign through the b step is selected to reflect 4 remote sensing wave bands of Extract Mineralized Alteration according to Kroes tower method, to in N * N window of selecting among the step c all effective background pixels being carried out principal component analysis (PCA) based on feature guiding, obtain PC1, PC2, PC3 and PC4 component then;
E. according to strong absorption and the high reflectance signature of Extract Mineralized Alteration spectrum 4 components that obtain in the steps d are judged, obtain the component at the unusual place of Extract Mineralized Alteration, the method for judgement is the component that just is judged as the unusual place of Extract Mineralized Alteration according to the component of the opposite in sign of the proper vector of the low wave band of the high wave band of reflectivity and reflectivity;
F. the component at the unusual place of Extract Mineralized Alteration that obtains through step e is carried out the judgement of light and shade pixel, the method for judgement be when the high wave band characteristic of correspondence vector value of reflectivity on the occasion of, then show as bright image unit unusually, otherwise then show as dark pixel;
G. unusual light and shade pixel conversion and assignment: the arbitrary wave band with the remote sensing image A of input is that template copies, and makes up an abnormal image B; Unusually, then the value of the proper vector of this component center pixel is directly composed the pixel of giving abnormal image B institute correspondence position during for bright image unit when what judge through the f step, when unusual on duty of judging through the f step with-1 with the proper vector of this component center pixel then during for dark pixel, dark pixel is converted to composes the pixel of giving abnormal image B institute correspondence position after the bright image unit again, obtain final abnormal image C.
H. export final abnormal image C.
Method of the present invention utilizes the locally variable window that the view picture image is divided into several independent statistical analysis unit, and at the inner interfere informations such as water body, cloud, ice and snow of getting rid of of each separate unit, and then carry out judgement and the selection of principal component analysis (PCA) calculating and the unusual component of Extract Mineralized Alteration, can effectively reduce neighbourhood noise, improve the precision that Extract Mineralized Alteration Information is extracted.More than 620 mineral deposit (point) on ground such as the Kunlun, east, statistics Qinghai, north, metallogenic belt, three rivers section, the Qinling Mountains, east, Henan, East Tianshan of Xinjiang and Altay mountain, East Africa, Australia, the correlation ratio of finding the unusual and iron of Extract Mineralized Alteration that conventional Crosta method extracts, plumbous zinc, gold, copper, molybdenum deposit (point) is about 45%, and the Extract Mineralized Alteration that the inventive method is extracted correlation ratio unusual and with iron, plumbous zinc, gold, copper, molybdenum deposit (point) is about 51%, can effectively improve the precision that Extract Mineralized Alteration Information is extracted.The inventive method not only has stronger screening and filtering function to interfere information than conventional Crosta method, and weak alteration Information is also had stronger recognition capability.The inventive method can also further be excavated the potentiality of remotely-sensed data Extract Mineralized Alteration drawing such as TM and ASTER, improves its indication ability in geologic prospecting.
Embodiment
For clear more technical scheme of the present invention, specific embodiment is described in detail as follows:
The invention process case is selected TM and the concrete object of implementing of ASTER remote sensing image conduct in certain plumbous zinc polymetallic ore district of section, north, San Jiang metallogenic belt, Qinghai.
Implementation step is as follows:
A. import a width of cloth remote sensing image A, select a width of cloth TM and ASTER remote sensing image as input one width of cloth remote sensing image A respectively.
B. the remote sensing image A of input is carried out differentiation and the sign of effective background pixel, wherein effectively the background pixel refers to be non-corrupt data pixel, non-water body pixel, non-cloud pixel, non-ice and snow pixel and non-vegetation pixel.Can utilize the normalized differential vegetation index NDVI of TM and ASTER and the vegetation of normalization water body index NDWI and water body to identify (as table 1)
The vegetation index of table 1 TM and ASTER and water body index
C. the locally variable window is selected.The remote sensing image that carries out effective background pixel sign through the b step active window with N * N pixel is scanned and adds up each pixel, wherein N is the odd number between 11 to 31, this statistics is centered by the pixel of analyzing spot place, search for since 11 * 11 windows, enlarge search window successively, largest extension to 31 * 31 windows, the ratio that accounts for whole window pixel quantity up to the quantity of effective background pixel reaches 51%;
The characteristic wave bands combination that the main Extract Mineralized Alteration Crosta of table 2 analyzes
Annotate: ++ be the height reflection,--be strong absorption
D. the remote sensing image that carries out effective background pixel sign through the b step is selected to reflect 4 remote sensing wave bands (as table 2) of Extract Mineralized Alteration according to Kroes tower method, to in N * N window of selecting among the step c all effective background pixels being carried out principal component analysis (PCA) based on feature guiding, obtain PC1, PC2, PC3 and PC4 component then;
E. according to strong absorption and the high reflectance signature of Extract Mineralized Alteration spectrum 4 components that obtain in the steps d are judged, obtain the component at the unusual place of Extract Mineralized Alteration, the method for judgement is the component that just is judged as the unusual place of Extract Mineralized Alteration according to the component of the opposite in sign of the proper vector of the low wave band of the high wave band of reflectivity and reflectivity.For example when the opposite in sign of the proper vector of TM1 and TM3 wave band place component, then this component is exactly the component that TM iron dyes unusual place; When the opposite in sign of the proper vector of ASTER6 and ASTER7 wave band place component, then this component is exactly the component at the unusual place of ASTER kaolinization.
F. the component at the unusual place of Extract Mineralized Alteration that obtains through step e is carried out the judgement of light and shade pixel, the method for judgement be when the high wave band characteristic of correspondence vector value of reflectivity on the occasion of, then show as bright image unit unusually, otherwise then show as dark pixel.For example the proper vector of dying TM3 wave band in the component of unusual place when TM iron on the occasion of the time, then TM iron dyes unusually and showing as bright image unit in unusual component gray scale image, otherwise then shows as dark pixel.
G. unusual light and shade pixel is changed and assignment.Arbitrary wave band of remote sensing image A with input is that template copies, and makes up an abnormal image B; Unusually, then the value of the proper vector of this component center pixel is directly composed the pixel of giving abnormal image B institute correspondence position during for bright image unit when what judge through the f step, when unusual on duty of judging through the f step with-1 with the proper vector of this component center pixel then during for dark pixel, dark pixel is converted to composes the pixel of giving abnormal image B institute correspondence position after the bright image unit again, obtain final abnormal image C.
H. export final abnormal image C.
The TM iron that utilizes the inventive method and conventional Crosta method that above-mentioned case study on implementation is extracted dyes unusual and ASTER chapeau de fer abnormal results compares analysis, statistics such as table 3.The TM iron that the inventive method and conventional Crosta method are extracted dyes unusual average and is respectively 0.2004 and 0.1179, standard deviation is respectively 3.1319 and 3.5238, the TM iron that this method is extracted dyes the average of the conventional Crosta method of anomaly ratio and wants big by 69.7%, it is little by 12.5% that standard deviation is wanted, and illustrates that this method extraction TM iron dyes strong unusually more and concentrated.The unusual average of ASTER chapeau de fer that the inventive method and conventional Crosta method are extracted is respectively 0.1252 and 0.0476, standard deviation is respectively 3.0795 and 3.6237, the average of the conventional Crosta method of the ASTER chapeau de fer anomaly ratio that this method is extracted wants big by 163.0%, it is little by 54.4% that standard deviation is wanted, and illustrates that also the ASTER chapeau de fer of this method extraction is strong unusually more and concentrated.Find that by the analysis to above-mentioned case study on implementation this method has stronger screening and filtering function to interfere information, simultaneously weak alteration Information is had stronger recognition capability.
The alteration abnormal results statistics that table 3 the inventive method and conventional Crosta method are extracted
Claims (1)
1. the extracting method of an Extract Mineralized Alteration Information is characterized in that, may further comprise the steps:
A. import a width of cloth remote sensing image A;
B. the remote sensing image A of input is carried out differentiation and the sign of effective background pixel, wherein effectively the background pixel refers to be non-corrupt data pixel, non-water body pixel, non-cloud pixel, non-ice and snow pixel and non-vegetation pixel;
C. the locally variable window is selected.The remote sensing image that carries out effective background pixel sign through the b step active window with N * N pixel is scanned and adds up each pixel, wherein N is the odd number between 11 to 31, this statistics is centered by the pixel of analyzing spot place, search for since 11 * 11 windows, enlarge search window successively, largest extension to 31 * 31 windows, the ratio that accounts for whole window pixel quantity up to the quantity of effective background pixel reaches 51%;
D. the remote sensing image that carries out effective background pixel sign through the b step is selected to reflect 4 remote sensing wave bands of Extract Mineralized Alteration according to Kroes tower method, to in N * N window of selecting among the step c all effective background pixels being carried out principal component analysis (PCA) based on feature guiding, obtain PC1, PC2, PC3 and PC4 component then;
E. according to strong absorption and the high reflectance signature of Extract Mineralized Alteration spectrum 4 components that obtain in the steps d are judged, obtain the component at the unusual place of Extract Mineralized Alteration, the method for judgement is the component that just is judged as the unusual place of Extract Mineralized Alteration according to the component of the opposite in sign of the proper vector of the low wave band of the high wave band of reflectivity and reflectivity;
F. the component at the unusual place of Extract Mineralized Alteration that obtains through step e is carried out the judgement of light and shade pixel, the method for judgement be when the high wave band characteristic of correspondence vector value of reflectivity on the occasion of, then show as bright image unit unusually, otherwise then show as dark pixel;
G. unusual light and shade pixel conversion and assignment: the arbitrary wave band with the remote sensing image A of input is that template copies, and makes up an abnormal image B; Unusually, then the value of the proper vector of this component center pixel is directly composed the pixel of giving abnormal image B institute correspondence position during for bright image unit when what judge through the f step, when unusual on duty of judging through the f step with-1 with the proper vector of this component center pixel then during for dark pixel, dark pixel is converted to composes the pixel of giving abnormal image B institute correspondence position after the bright image unit again, obtain final abnormal image C;
H. export final abnormal image C.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210586855.8A CN103207415B (en) | 2012-12-30 | 2012-12-30 | A kind of extracting method of Extract Mineralized Alteration Information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210586855.8A CN103207415B (en) | 2012-12-30 | 2012-12-30 | A kind of extracting method of Extract Mineralized Alteration Information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103207415A true CN103207415A (en) | 2013-07-17 |
| CN103207415B CN103207415B (en) | 2015-08-26 |
Family
ID=48754690
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210586855.8A Expired - Fee Related CN103207415B (en) | 2012-12-30 | 2012-12-30 | A kind of extracting method of Extract Mineralized Alteration Information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103207415B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103852807A (en) * | 2014-03-20 | 2014-06-11 | 青海省地质矿产研究所 | High-altitude permafrost region porphyry type copper polymetallic mine exploration technology combination method |
| CN105068136A (en) * | 2015-07-27 | 2015-11-18 | 中国地质调查局武汉地质调查中心 | Potential positioning evaluation method for copper and gold mine of Indo-China peninsula demonstration zone based on multi-source information |
| CN105574621A (en) * | 2016-01-18 | 2016-05-11 | 中国地质科学院矿产资源研究所 | Porphyry copper ore prediction system and method based on remote sensing alteration abnormity |
| CN107192673A (en) * | 2017-05-18 | 2017-09-22 | 成都理工大学 | Integrated geological mapping method based on ASTER and underground core spectral measurement technology |
| CN109406405A (en) * | 2018-10-11 | 2019-03-01 | 核工业北京地质研究院 | A kind of high-definition remote sensing quantitative estimation method suitable for salt lake water body salinity |
| CN113406041A (en) * | 2021-05-31 | 2021-09-17 | 核工业北京地质研究院 | Method for obtaining key altered mineral combination of sodium-intercrossed rock type uranium ore |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101539998A (en) * | 2009-04-29 | 2009-09-23 | 中国地质科学院矿产资源研究所 | Alteration remote sensing abnormity extraction method and system |
| CN101788685A (en) * | 2010-02-11 | 2010-07-28 | 中国土地勘测规划院 | Remote sensing earthquake damage information extracting and digging method based on pixels |
-
2012
- 2012-12-30 CN CN201210586855.8A patent/CN103207415B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101539998A (en) * | 2009-04-29 | 2009-09-23 | 中国地质科学院矿产资源研究所 | Alteration remote sensing abnormity extraction method and system |
| CN101788685A (en) * | 2010-02-11 | 2010-07-28 | 中国土地勘测规划院 | Remote sensing earthquake damage information extracting and digging method based on pixels |
Non-Patent Citations (8)
| Title |
|---|
| 张远飞等人: "基于遥感信息图像提取矿化蚀变信息", 《有色金属矿产与勘查》, vol. 8, no. 6, 1 December 1999 (1999-12-01), pages 604 - 606 * |
| 李建国等人: "基于ETM与ASTER数据的矿化蚀变信息提取方法研究", 《地质调查与研究》, vol. 30, no. 3, 1 September 2007 (2007-09-01), pages 234 - 240 * |
| 沈利霞等人: "不同植被覆盖率地区遥感矿化蚀变提取研究", 《现代地质》, vol. 22, no. 2, 1 April 2008 (2008-04-01), pages 293 - 298 * |
| 王亚红等人: "利用数据进行矿化蚀变信息提取的方法研究", 《遥感信息》, 1 June 2009 (2009-06-01), pages 68 - 72 * |
| 王蔚等人: "应用ASTER数据提取矽卡岩型锡矿蚀变遥感异常信息", 《国土资源遥感》, no. 92, 15 March 2012 (2012-03-15), pages 155 - 159 * |
| 薛红等人: "基于光谱和纹理的SVM矿化蚀变信息提取研究", 《地质找矿论丛》, vol. 23, no. 3, 1 September 2008 (2008-09-01), pages 254 - 259 * |
| 陈劲松等人: "高寒山区矿化蚀变信息提取应用研究", 《遥感信息》, 1 February 2012 (2012-02-01), pages 77 - 80 * |
| 韩玲等人: "利用遥感技术队新疆西天上地区矿化蚀变信息的提取1", 《遥感信息》, 1 June 2007 (2007-06-01), pages 49 - 51 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103852807A (en) * | 2014-03-20 | 2014-06-11 | 青海省地质矿产研究所 | High-altitude permafrost region porphyry type copper polymetallic mine exploration technology combination method |
| CN105068136A (en) * | 2015-07-27 | 2015-11-18 | 中国地质调查局武汉地质调查中心 | Potential positioning evaluation method for copper and gold mine of Indo-China peninsula demonstration zone based on multi-source information |
| CN105574621A (en) * | 2016-01-18 | 2016-05-11 | 中国地质科学院矿产资源研究所 | Porphyry copper ore prediction system and method based on remote sensing alteration abnormity |
| CN107192673A (en) * | 2017-05-18 | 2017-09-22 | 成都理工大学 | Integrated geological mapping method based on ASTER and underground core spectral measurement technology |
| CN107192673B (en) * | 2017-05-18 | 2020-11-06 | 成都理工大学 | Integrated geological mapping method based on ASTER and underground core spectral measurement technology |
| CN109406405A (en) * | 2018-10-11 | 2019-03-01 | 核工业北京地质研究院 | A kind of high-definition remote sensing quantitative estimation method suitable for salt lake water body salinity |
| CN113406041A (en) * | 2021-05-31 | 2021-09-17 | 核工业北京地质研究院 | Method for obtaining key altered mineral combination of sodium-intercrossed rock type uranium ore |
| CN113406041B (en) * | 2021-05-31 | 2023-03-17 | 核工业北京地质研究院 | Method for obtaining key altered mineral combination of sodium-substituted rock type uranium ore |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103207415B (en) | 2015-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103207415A (en) | Method for extracting mineralization-alteration information | |
| Zhang et al. | Lithologic and mineral information extraction for gold exploration using ASTER data in the south Chocolate Mountains (California) | |
| Rajesh | Application of remote sensing and GIS in mineral resource mapping-An overview | |
| Shen et al. | Water body extraction from Landsat ETM+ imagery using adaboost algorithm | |
| CN110806605B (en) | Remote sensing ore searching method for rare earth-uranium ore in high-latitude and high-cold area | |
| CN111141698A (en) | Lithology classification method based on thermal infrared emissivity | |
| CN102426625A (en) | Method for extracting porphyry copper typical altered zone mineral information based on ASTER satellite data | |
| CN108573206A (en) | Metallogenic prognosis method based on analytic hierarchy process (AHP) | |
| Vural et al. | Remote sensing technique for capturing and exploration of mineral deposit sites in Gumushane metallogenic province, NE Turkey | |
| CN110147778A (en) | Rare Earth Mine exploits recognition methods, device, equipment and storage medium | |
| Alasta | Using Remote Sensing data to identify iron deposits in central western Libya | |
| CN109359621B (en) | Pegmatite type lithium ore prospecting method based on multi-source remote sensing data | |
| Joseph et al. | Application of remote sensing method for geological interpretation of Sokoto Plain, Nigeria | |
| Pour et al. | Evaluation of earth observing-1 (EO1) data for lithological and hydrothermal alteration mapping: A case study from Urumieh-Dokhtar volcanic belt, SE Iran | |
| Soe et al. | Remote sensing and GIS studies of alteration and predictive mineral exploration in the Central Volcanic Arc, Myanmar | |
| Abhary et al. | Mapping hydrothermal mineral deposits using PCA and BR methods in BAFT 1: 100000 geological sheet, Iran | |
| Wu et al. | A fully automatic method to extract rare earth mining areas from Landsat Images | |
| CN114386498A (en) | Method for extracting iron dyeing information by using WorldView-3 remote sensing data | |
| Kujjo | Application of remote sensing for gold exploration in the Nuba Mountains, Sudan | |
| Zhao et al. | Ore-and bio-geochemical survey based on the landsat remotely sensed data in and around the Dexing porphyry copper-polymetal ore-field, southeastern China | |
| CN111160182A (en) | Method for identifying surface calcium-bearing rock of ancient river | |
| Gloaguen et al. | Towards 4D Virtual Outcrops with Hyperspectral Imaging | |
| Elhamdouni et al. | RETRACTED ARTICLE: Automatic geological mapping using remote sensing data: case of the Zgounder deposit (Anti-Atlas, Morocco) | |
| Saeedi et al. | Prospecting for Gold Mineralization Using Geochemical, Mineralogical, and WorldView-2 Data: Siyah Jangal Area Case Study, Northern Taftan Volcano, SE Iran | |
| Aisabokhae et al. | Supervised classification of Landsat-8 band ratio images for geological interpretation of Sokoto, Nigeria |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150826 Termination date: 20161230 |