CN103984940A - 一种基于高光谱数据识别赤铁矿化的方法 - Google Patents
一种基于高光谱数据识别赤铁矿化的方法 Download PDFInfo
- Publication number
- CN103984940A CN103984940A CN201410242036.0A CN201410242036A CN103984940A CN 103984940 A CN103984940 A CN 103984940A CN 201410242036 A CN201410242036 A CN 201410242036A CN 103984940 A CN103984940 A CN 103984940A
- Authority
- CN
- China
- Prior art keywords
- haematization
- spectral data
- end member
- ferritization
- wave bands
- 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
本发明属于一种赤铁矿化识别方法,具体公开一种基于高光谱数据识别赤铁矿化的方法,该方法包括如下步骤:步骤(1)获取高光谱影像数据和预处理;步骤(2)高光谱影像数据特征波段选择;步骤(3)高光谱数据特征波段图像端元提取;步骤(4)建立光谱特征识别规则区分识别赤铁矿化和褐铁矿化端元;步骤(5)利用混合协调匹配滤波对赤铁矿化和褐铁矿化端元进行填图。本发明的方法能够识别赤铁矿化和褐铁矿化,识别的精度高,矿物的检出限低。
Description
技术领域
本发明属于一种赤铁矿化识别方法,具体涉及一种基于高光谱数据识别赤铁矿化的方法。
背景技术
赤铁矿化又称为红色蚀变,是寻找中、低温热液铀矿床的特殊标志之一,在其他中、低温铜、金矿床中也常有发现,是重要的成矿要素。因此,如何准确识别和提取与成矿密切相关的赤铁矿化分布,对圈定成矿靶区具有重要的现实意义。当前,在遥感蚀变信息提取领域,对铁染信息的提取已开展多年,发展已较为成熟,但鲜有见到进一步对铁染信息进行区分进而识别出赤铁矿化的案例。由于提取铁染信息中同时包含的褐铁矿化,多为其他矿物风化形成,实际找矿意义并不大,影响了遥感信息的实用价值。因此,遥感图像上提取出的铁染信息中究竟有多少是真正与成矿作用密切相关的,成为遥感地质工作者需要进一步探究的问题。
目前对铁染信息的提取多采用多光谱TM、中光谱分辨率Aster数据等,方法归纳主要有以下几种:(1)波段运算法,如波段比值;(2)统计分析法,如主成分分析法、独立成分分析法等;(3)综合法,如将波段比值和主成分分析法相结合。波段运算法通过波段加减组合、比值运算能够增强不同地物之间的差异,但是由于铁化信息的光谱特征较宽缓,特征波长位置较不明确,难以构建准确的运算法则;统计分析法通过数学变换将图像数据进行集中和压缩,将信息集中到少数几个互不相关的新变量中,但是这些新变量不像原始变量含义那么清楚、确切,带有一定的模糊性;综合法尽管将两者进行结合,但是算法中的不确定性仍然存在,并且在运算过程中容易引入其他噪声。总之,这些方法的共同缺陷是:物理意义不明确,且容易引入噪声,最根本的缺陷是难以准确区分识别赤铁矿化和褐铁矿化,因此,有必要开发赤铁矿化识别的新方法。
发明内容
本发明的目的在于提供一种基于高光谱数据识别赤铁矿化的方法,该方法能够识别赤铁矿化和褐铁矿化,识别的精度高,矿物的检出限低。
实现本发明目的的技术方案:一种基于高光谱数据识别赤铁矿化的方法,该方法包括如下步骤:
步骤(1)获取高光谱影像数据和预处理;
步骤(2)高光谱影像数据特征波段选择;
步骤(3)高光谱数据特征波段图像端元提取;
步骤(4)建立光谱特征识别规则区分识别赤铁矿化和褐铁矿化端元;
步骤(5)利用混合协调匹配滤波对赤铁矿化和褐铁矿化端元进行填图。
所述的步骤(1)中对高光谱影像数据进行预处理包括辐射校正、几何校正、影像裁剪、图像掩膜,并对高光谱影像数据进行大气校正和光谱重建,得到高光谱数据的反射率光谱数据。
所述的步骤(2)中采用光谱重采样方法对高光谱数据进行特征波段选择,截取可见光-近红外光谱区间波段作为特征波段。
所述的步骤(3)中对上述步骤(2)中得到的高光谱数据特征波段进行最小噪声分离变换、多维空间旋转、聚类,初步圈定高光谱数据特征波段图像中的端元。
所述的步骤(4)中根据赤铁矿化和褐铁矿化的高光谱数据特征波段区间上的光谱特征差异,建立两者的区分规则,依据规则区分识别赤铁矿化和褐铁矿化端元。
所述的步骤(5)中对述步骤(4)中对识别出的赤铁矿化和褐铁矿化端元进行混合协调匹配滤波处理,分别得到赤铁矿化和褐铁矿化的匹配结果灰度图和不可行性灰度图,实现图像识别填图。
本发明的有益技术效果在于。(1)本发明充分利用高光谱图像图谱合一的特点,从光谱角度建立赤铁矿化和褐铁矿化的光谱特征区分规则,与传统方法相比,通过高光谱遥感填图技术,把遥感信息提取领域引入铁化信息的识别分类中,提高了识别的精度。(2)针对高光谱数据波段较多、信息量巨大的特点,针对识别和提取的铁化信息光谱特征波段区间,对预处理和光谱重建好的高光谱数据进行特征波段选择,达到降维、优化信息量的目的。(3)通过最小噪声分离变换实现信噪分离;通过观察特征值和相关图像,确定数据的内在维数,有利于进一步在多维空间中进行聚类和圈定各类端元。(4)采用线性混合分解与匹配滤波相结合的方法进行混合调制匹配滤波,综合了匹配滤波不需要其他背景端元光谱的优点和线性混合分解的像元中各端元的含量为正且总和为1的约束条件,因而降低矿物的检出限。
附图说明
图1为本发明所提供的一种的基于高光谱数据识别赤铁矿化的方法的流程图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。
如图1所示,一种基于高光谱数据识别赤铁矿化的方法,该方法包括如下步骤:
步骤(1)获取高光谱影像数据和预处理
利用机载或星载成像光谱仪获取高光谱影像数据,成像光谱仪的光谱范围至少覆盖可见光-近红外波段,如机载Casi、Hymap、星载Hyperion等传感器,航空飞行或卫星过境时间选择正午或接近正午,天气晴朗无云。
利用ENVI软件对高光谱影像数据进行预处理,预处理包括辐射校正、几何校正、影像裁剪、图像掩膜等,并对高光谱影像数据进行大气校正和光谱重建,得到高光谱数据的反射率光谱数据。
步骤(2)高光谱影像数据特征波段选择
采用光谱重采样方法对高光谱数据进行特征波段选择,截取可见光-近红外光谱区间400-1000nm之间的波段作为特征波段。
步骤(3)高光谱数据特征波段图像端元提取
对上述步骤(2)中得到的高光谱数据特征波段进行最小噪声分离变换,得到最小噪声分离变换后的各分量;通过观察各分量特征值和相关图像,确定数据的内在维数;将像元投影到在多维空间中旋转、聚类,交互式选择分布于数据云主体的外围或末端的集群,初步圈定高光谱数据特征波段图像中的端元。
步骤(4)建立光谱特征识别规则区分识别赤铁矿化和褐铁矿化端元
根据赤铁矿化和褐铁矿化的高光谱数据特征波段区间上的光谱特征差异,建立两者的区分规则,依据规则区分识别赤铁矿化和褐铁矿化端元。判别规则依据端元光谱的特征波峰、谷的反射率值和波长位置。
获取端元光谱曲线750nm附近处的反射率最大值R1,1000nm附近处的反射率最大值R2;获取750-1000nm间反射率最小值R3,设R3的波长位置为λ3,规则如下:
①若R3<R1<R2,且λ3∈(800,900),判断为赤铁矿化,λ波长单位为纳米;
②若R3<R2<R1,且λ3∈(900,1000),判断为褐铁矿化,λ波长单位为纳米。
步骤(5)利用混合协调匹配滤波对赤铁矿化和褐铁矿化端元进行填图
对述步骤(4)中对识别出的赤铁矿化和褐铁矿化端元进行混合协调匹配滤波处理,分别得到赤铁矿化和褐铁矿化的匹配结果灰度图和不可行性灰度图。由于和端元最匹配的像元有一个大于背景分布值的较高的匹配滤波值和一个较低的不可行性值,通过建立匹配结果灰度和不可行性灰度的二维散点图,圈出图像上赤铁矿化和褐铁矿化最匹配的像元,实现赤铁矿化和褐铁矿化端元图像识别填图。
混合调制匹配滤波采用线性混合分解与匹配滤波相结合的方法。
上面结合附图和实施例对本发明作了详细说明,但是本发明并不限于上述实施例,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可以采用现有技术。
Claims (6)
1.一种基于高光谱数据识别赤铁矿化的方法,其特征在于,该方法包括如下步骤:
步骤(1)获取高光谱影像数据和预处理;
步骤(2)高光谱影像数据特征波段选择;
步骤(3)高光谱数据特征波段图像端元提取;
步骤(4)建立光谱特征识别规则区分识别赤铁矿化和褐铁矿化端元;
步骤(5)利用混合协调匹配滤波对赤铁矿化和褐铁矿化端元进行填图。
2.根据权利要求1所述的一种,其特征在于基于高光谱数据识别赤铁矿化的方法:所述的步骤(1)中对高光谱影像数据进行预处理包括辐射校正、几何校正、影像裁剪、图像掩膜,并对高光谱影像数据进行大气校正和光谱重建,得到高光谱数据的反射率光谱数据。
3.根据权利要求2所述的一种,其特征在于基于高光谱数据识别赤铁矿化的方法:所述的步骤(2)中采用光谱重采样方法对高光谱数据进行特征波段选择,截取可见光-近红外光谱区间波段作为特征波段。
4.根据权利要求3所述的一种,其特征在于基于高光谱数据识别赤铁矿化的方法:所述的步骤(3)中对上述步骤(2)中得到的高光谱数据特征波段进行最小噪声分离变换、多维空间旋转、聚类,初步圈定高光谱数据特征波段图像中的端元。
5.根据权利要求4所述的一种,其特征在于基于高光谱数据识别赤铁矿化的方法:所述的步骤(4)中根据赤铁矿化和褐铁矿化的高光谱数据特征波段区间上的光谱特征差异,建立两者的区分规则,依据规则区分识别赤铁矿化和褐铁矿化端元。
6.根据权利要求5所述的一种,其特征在于基于高光谱数据识别赤铁矿化的方法:所述的步骤(5)中对述步骤(4)中对识别出的赤铁矿化和褐铁矿化端元进行混合协调匹配滤波处理,分别得到赤铁矿化和褐铁矿化的匹配结果灰度图和不可行性灰度图,实现图像识别填图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410242036.0A CN103984940B (zh) | 2014-06-03 | 2014-06-03 | 一种基于高光谱数据识别赤铁矿化的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410242036.0A CN103984940B (zh) | 2014-06-03 | 2014-06-03 | 一种基于高光谱数据识别赤铁矿化的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103984940A true CN103984940A (zh) | 2014-08-13 |
CN103984940B CN103984940B (zh) | 2017-12-26 |
Family
ID=51276902
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410242036.0A Active CN103984940B (zh) | 2014-06-03 | 2014-06-03 | 一种基于高光谱数据识别赤铁矿化的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103984940B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021653A (zh) * | 2016-05-06 | 2016-10-12 | 华南农业大学 | 一种ndvi时间序列重构的方法及系统 |
CN107589094A (zh) * | 2017-09-20 | 2018-01-16 | 鞍钢集团矿业有限公司 | 基于光谱特征的鞍山式铁矿石类型确定方法 |
CN108090434A (zh) * | 2017-12-13 | 2018-05-29 | 赣州好朋友科技有限公司 | 一种矿石快速识别方法 |
CN108333129A (zh) * | 2018-02-12 | 2018-07-27 | 中国科学院遥感与数字地球研究所 | 一种蚀变信息获取方法及系统 |
CN108956482A (zh) * | 2017-05-27 | 2018-12-07 | 核工业北京地质研究院 | 一种火山热液活动中心的高光谱遥感快速识别方法 |
CN108956483A (zh) * | 2017-05-27 | 2018-12-07 | 核工业北京地质研究院 | 一种火山岩型铀矿化地段的航空高光谱遥感快速圈定方法 |
CN109283148A (zh) * | 2018-09-30 | 2019-01-29 | 核工业北京地质研究院 | 一种基于光谱信息自动识别岩石矿物的方法 |
CN109657637A (zh) * | 2018-12-27 | 2019-04-19 | 核工业北京地质研究院 | 一种利用casi图像区别不同岩石中赤铁矿的方法 |
CN109740489A (zh) * | 2018-12-27 | 2019-05-10 | 核工业北京地质研究院 | 一种利用casi图像识别赤铁矿的方法 |
CN113406041A (zh) * | 2021-05-31 | 2021-09-17 | 核工业北京地质研究院 | 一种钠交代岩型铀矿关键蚀变矿物组合获取方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020015151A1 (en) * | 2000-07-27 | 2002-02-07 | Gorin Brian A. | Spectral drift and correction technique for hyperspectral imaging systems |
CN1595203A (zh) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | 高光谱矿物分层谱系识别方法 |
CN101871884A (zh) * | 2010-06-02 | 2010-10-27 | 中国国土资源航空物探遥感中心 | 多景aster遥感数据大气校正与区域性矿物填图方法 |
CN102012528A (zh) * | 2010-11-23 | 2011-04-13 | 北京理工大学 | 一种稀疏植被地区的高光谱遥感油气勘探方法 |
CN103175801A (zh) * | 2013-03-14 | 2013-06-26 | 中国国土资源航空物探遥感中心 | 一种批量化与自动化的高光谱遥感矿物填图方法 |
US20140118738A1 (en) * | 2012-10-31 | 2014-05-01 | Corning Incorporated | Hyperspectral imaging system, monolithic spectrometer and methods for manufacturing the monolithic spectrometer |
-
2014
- 2014-06-03 CN CN201410242036.0A patent/CN103984940B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020015151A1 (en) * | 2000-07-27 | 2002-02-07 | Gorin Brian A. | Spectral drift and correction technique for hyperspectral imaging systems |
CN1595203A (zh) * | 2004-06-29 | 2005-03-16 | 中国国土资源航空物探遥感中心 | 高光谱矿物分层谱系识别方法 |
CN101871884A (zh) * | 2010-06-02 | 2010-10-27 | 中国国土资源航空物探遥感中心 | 多景aster遥感数据大气校正与区域性矿物填图方法 |
CN102012528A (zh) * | 2010-11-23 | 2011-04-13 | 北京理工大学 | 一种稀疏植被地区的高光谱遥感油气勘探方法 |
US20140118738A1 (en) * | 2012-10-31 | 2014-05-01 | Corning Incorporated | Hyperspectral imaging system, monolithic spectrometer and methods for manufacturing the monolithic spectrometer |
CN103175801A (zh) * | 2013-03-14 | 2013-06-26 | 中国国土资源航空物探遥感中心 | 一种批量化与自动化的高光谱遥感矿物填图方法 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021653A (zh) * | 2016-05-06 | 2016-10-12 | 华南农业大学 | 一种ndvi时间序列重构的方法及系统 |
CN108956482A (zh) * | 2017-05-27 | 2018-12-07 | 核工业北京地质研究院 | 一种火山热液活动中心的高光谱遥感快速识别方法 |
CN108956483A (zh) * | 2017-05-27 | 2018-12-07 | 核工业北京地质研究院 | 一种火山岩型铀矿化地段的航空高光谱遥感快速圈定方法 |
CN107589094B (zh) * | 2017-09-20 | 2020-02-21 | 鞍钢集团矿业有限公司 | 基于光谱特征的鞍山式铁矿石类型确定方法 |
CN107589094A (zh) * | 2017-09-20 | 2018-01-16 | 鞍钢集团矿业有限公司 | 基于光谱特征的鞍山式铁矿石类型确定方法 |
CN108090434A (zh) * | 2017-12-13 | 2018-05-29 | 赣州好朋友科技有限公司 | 一种矿石快速识别方法 |
CN108333129A (zh) * | 2018-02-12 | 2018-07-27 | 中国科学院遥感与数字地球研究所 | 一种蚀变信息获取方法及系统 |
CN108333129B (zh) * | 2018-02-12 | 2020-04-21 | 中国科学院遥感与数字地球研究所 | 一种蚀变信息获取方法及系统 |
CN109283148A (zh) * | 2018-09-30 | 2019-01-29 | 核工业北京地质研究院 | 一种基于光谱信息自动识别岩石矿物的方法 |
CN109657637A (zh) * | 2018-12-27 | 2019-04-19 | 核工业北京地质研究院 | 一种利用casi图像区别不同岩石中赤铁矿的方法 |
CN109740489A (zh) * | 2018-12-27 | 2019-05-10 | 核工业北京地质研究院 | 一种利用casi图像识别赤铁矿的方法 |
CN109657637B (zh) * | 2018-12-27 | 2022-07-26 | 核工业北京地质研究院 | 一种利用casi图像区别不同岩石中赤铁矿的方法 |
CN113406041A (zh) * | 2021-05-31 | 2021-09-17 | 核工业北京地质研究院 | 一种钠交代岩型铀矿关键蚀变矿物组合获取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103984940B (zh) | 2017-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103984940A (zh) | 一种基于高光谱数据识别赤铁矿化的方法 | |
Peyghambari et al. | Hyperspectral remote sensing in lithological mapping, mineral exploration, and environmental geology: an updated review | |
van der Meero et al. | Cross correlogram spectral matching: application to surface mineralogical mapping by using AVIRIS data from Cuprite, Nevada | |
CN104050252B (zh) | 一种高光谱遥感蚀变信息提取方法 | |
Yang et al. | TeaNet: Deep learning on Near-Infrared Spectroscopy (NIR) data for the assurance of tea quality | |
Dhumal et al. | Classification of Crops from remotely sensed Images: AnOverview | |
Kumar et al. | Deep remote sensing methods for methane detection in overhead hyperspectral imagery | |
Shahriari et al. | Selection of less biased threshold angles for SAM classification using the real value–area fractal technique | |
Galdames et al. | Rock lithological classification by hyperspectral, range 3D and color images | |
Li et al. | GIS-based detection of grain boundaries | |
Uzar | Automatic building extraction with multi-sensor data using rule-based classification | |
CN114575925B (zh) | 一种基于时空大数据的矿井安全监测和预警的方法 | |
CN108956482A (zh) | 一种火山热液活动中心的高光谱遥感快速识别方法 | |
Yamamoto et al. | An automated method for crater counting using rotational pixel swapping method | |
Ishidoshiro et al. | Geological mapping by combining spectral unmixing and cluster analysis for hyperspectral data | |
Zhang et al. | A novel multitemporal cloud and cloud shadow detection method using the integrated cloud Z-scores model | |
Chen et al. | Improving the accuracy of hydrothermal alteration mapping based on image fusion of ASTER and Sentinel-2A data: a case study of Pulang Cu deposit, Southwest China | |
Tripathi et al. | Comparative evaluation threshold parameters of Spectral Angle Mapper (SAM) for mapping of Chhabadiya Talc Minerals, Jahajpur, Bhilwara, India using Hyperion hyperspectral remote sensing data | |
Ouerghemmi et al. | Urban objects classification by spectral library: Feasibility and applications | |
Sarath et al. | A Study on Hyperspectral Remote Sensing | |
Zhao et al. | SRA-CEM: An improved CEM target detection algorithm for hyperspectral images based on sub-region analysis | |
Wolfe et al. | Hyperspectral analytics in envi target detection and spectral mapping methods | |
CN100561145C (zh) | 一种高光谱矿物填图方法 | |
JP6060377B2 (ja) | 画像クラスタ分析装置 | |
Zhang et al. | Feature Bands Based Unsupervised Hyperspectral Underwater Target Detection Near the Coastline |
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 |