CN103175801A - 一种批量化与自动化的高光谱遥感矿物填图方法 - Google Patents
一种批量化与自动化的高光谱遥感矿物填图方法 Download PDFInfo
- Publication number
- CN103175801A CN103175801A CN2013100819668A CN201310081966A CN103175801A CN 103175801 A CN103175801 A CN 103175801A CN 2013100819668 A CN2013100819668 A CN 2013100819668A CN 201310081966 A CN201310081966 A CN 201310081966A CN 103175801 A CN103175801 A CN 103175801A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- continuum
- image
- mineral
- reference spectra
- 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
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种批量化与自动化的高光谱遥感矿物填图方法,它有七大步骤:步骤一、光谱数据读入;步骤二、光谱波段截取;步骤三、光谱连续统去除;步骤四、影像光谱吸收深度计算;步骤五、影像光谱的标准化;步骤六、影像光谱与参考光谱相似度计算;步骤七、根据步骤五与步骤六的计算结果提取矿物异常。该方法无需从高光谱遥感数据中提取参考光谱并进行人工判读,所有数据的处理采用相同的参考光谱,便于实现数据的批量化与自动化处理,矿物异常提取时综合利用了光谱相似度与特征参量,提高了矿物识别的可靠性与准确度。本发明在高光谱遥感地质勘查工程化应用领域里具有实用价值和广阔的应用前景。
Description
技术领域
本发明涉及一种批量化与自动化的高光谱遥感矿物填图方法,属于高光谱遥感矿物识别技术领域,它适用于对高光谱遥感地质勘查中大量数据进行批量化与自动化处理,并提取矿物种类,半定量估算矿物含量。
背景技术
高光谱遥感具有图谱合一的特点,在地质勘查中应用广泛,可根据矿物的精细光谱特征进行矿物信息的识别。常用的识别方法有三类:基于光谱专家知识的方法、基于光谱相似度的方法、基于光谱分解的方法。
基于光谱专家知识的方法的基本思路是,综合各种光谱专家知识(光谱吸收特征位置、深度、面积等),建立矿物识别的决策树从而完成矿物的识别,代表性的方法是“高光谱矿物分层谱系识别方法”(已获得国家发明专利,专利公开号为CN1595203);基于光谱相似度的方法的基本思路是,基于各种光谱相似性算法(光谱角、匹配滤波等),将待识别矿物的光谱与参考光谱进行直接对比,相似度最高的矿物即被认为是待识别的矿物,参考光谱可取自标准矿物光谱库(常用的为美国地质调查局-USGS发布的矿物光谱库),但由于影像光谱与标准矿物光谱库测试条件(光照等)差异甚大,识别准确度往往较低,实际应用中参考光谱一般采用基于循环迭代方法的提取算法(常用的有纯像元指数法)从影像光谱中提取,并对提取的光谱进行人工判读。上述两类方法实现过程中需要较多的人工参与,难以支撑工程化应用中大量数据的批量化与自动化处理。基于光谱分解的方法中的代表性的方法是“一种高光谱遥感矿物含量定量反演方法”(已获得国家发明专利,专利号为ZL200910143190.1),实现过程中人工参与较少,适合于批量化与自动化处理,同时提高了矿物识别的定量化程度,但运算速度较慢,且参考光谱数目不可大于波段数,应用中存在一定的局限。
发明内容
本发明的目的是,针对现有方法人工参与多、运算步骤多,难以适应工程化应用的问题,提出了一种批量化与自动化的高光谱遥感矿物填图方法,它是一种综合运用光谱相似度与特征参量的高光谱遥感矿物填图方法。光谱相似度计算前首先对影像光谱进行去壳、标准化,抑制影像光谱与参考光谱测试条件差异对光谱特征的影响,参考光谱直接取自标准矿物光谱库,无需利用循环迭代方法从影像中提取,减少了运算步骤与人工的参与,从而实现了大量大量数据的批量化与自动化处理。
技术解决方案
本发明的技术解决方案见附图1,包括以下七个步骤:
(1)光谱数据读入;
(2)光谱波段截取;
(3)光谱连续统去除;
(4)影像光谱吸收深度计算;
(5)影像光谱的标准化;
(6)影像光谱与参考光谱相似度计算;
(7)根据步骤五与步骤六的计算结果提取矿物异常。
步骤一:光谱数据读入
这些光谱数据包括高光谱遥感数据与参考光谱数据,其中高光谱遥感数据是通过机载或星载高光谱传感器获取的,参考光谱数据来自于标准矿物光谱库(常用的为美国地质调查局-USGS发布的矿物光谱库)。
步骤二:光谱波段截取
其截取方法为,截取接近700—1300nm谱段的光谱识别褐铁矿、针铁矿、赤铁矿、黄钾铁钒,截取接近2100—2320nm谱段的光谱识别白云母、蒙脱石、高岭石、地开石、石膏、明矾石、黄钾铁钒,叶腊石,截取接近2200—2400m谱段的光谱识别绿泥石、绿帘石、角闪石、滑石、蛇纹石、方解石、白云石。
步骤三:光谱连续统去除
其实现方法为,光谱连续统去除(Continum Removal)是一种用于分离光谱吸收特征的光谱分析方法,广泛应用高光谱遥感中去除背景影响并分离某些特定物质的吸收特征。连续统被定义为反射光谱曲线中反射峰之间的线性连接部分,连接折线在反射光谱峰值处的外角大于180度。光谱连续统去除即以反射光谱除以连续统光谱,如下公式所示。
其中,Rcr为连续统去除后的吸收光谱,R吸收光谱,Rc为吸收光谱的连续统。经过连续统去除之后,端点处反射率为1,端点之间反射率均小于1。影像光谱与参考光谱均需进行连续统去除。
步骤四:影像光谱吸收深度计算
其计算方法为,在经过波段截取、连续统去除的光谱上,寻找最小值,并利用数字“1”减去该“最小值”。
步骤五:影像光谱的标准化
其计算方法为,以参考光谱标准,对影像光谱进行标准化,以抑制因测试条件差异对光谱相似度计算造成的影响。标准化系数计算公式如下所示。
步骤六:影像光谱与参考光谱相似度计算
其计算方法为,步骤四计算的影像光谱吸收深度除以标准化后的影像光谱与去连续统后的参考光谱。计算公式如下所示。
步骤七:根据步骤五与步骤六的计算结果提取矿物异常;
设定相似度阈值,认为大于阈值的像元即为目标矿物,大于阈值的像元的影像光谱吸收深度(步骤四计算结果)即为目标矿物的相对丰度。若一个工作区有多景/航带数据,可在镶嵌后设定阈值,也可之前设定阈值。
(四)附图说明
附图1为本发明的流程图;
(五)具体实施方式
见附图1,为了更好的说明本发明的方法与步骤,利用某地区AVIRIS机载高光谱数据为例,进行了矿物填图试验。
(1)试验所用的设备为图形工作站,规格型号为Dell Precision4700,操作系统为WindowsXP(64位),CPU为2.66GHz,内容为32GB,硬盘为1T。
(2)具体步骤如下:
步骤一,读入美国Cuprite地区AVIRIS高光谱遥感数据以及参考光谱数据。AVIRIS高光谱遥感数据长×宽为400×350个像元,谱段范围为1990.8—2479nm,波段数为50。参考光谱数据取自USGS矿物光谱库,包括高铝白云母(吸收位置为2200.8nm)、中铝白云母(吸收位置为2210.8nm)、低铝白云母(吸收位置为2220.8nm)、高岭石、埃洛石、地开石、蒙脱石、石膏、明矾石、黄钾铁钒、绿泥石、绿帘石、方解石、白云石、蛇纹石、角闪石;
步骤二,截取2101—2320.3nm波段数据识别高铝白云母(吸收位置为2200.8nm)、中铝白云母(吸收位置为2210.8nm)、低铝白云母(吸收位置为2220.8nm)、高岭石、埃洛石、地开石、蒙脱石、石膏、明矾石、黄钾铁钒,截取2200.8—2409.6nm波段数据识别绿泥石、绿帘石、方解石、白云石、蛇纹石、角闪石;
步骤三,光谱连续统去除;
连续统被定义为反射光谱曲线中反射峰之间的线性连接部分,连接折线在反射光谱峰值处的外角大于180度。光谱连续统去除即以反射光谱除以连续统光谱,如下公式所示。
其中,Rcr为连续统去除后的吸收光谱,R吸收光谱,Rc为吸收光谱的连续统。经过连续统去除之后,端点处反射率为1,端点之间反射率均小于1。影像光谱与参考光谱均需进行连续统去除。
步骤四,影像光谱吸收深度计算;
计算方法为,在经过波段截取、连续统去除的光谱上,寻找最小值,并利用数字“1”减去该“最小值”。
步骤五,影像光谱的标准化;
计算方法为,以参考光谱标准,对影像光谱进行标准化,以抑制因测试条件差异对光谱相似度计算造成的影响。标准化系数计算公式如下所示。
步骤六,影像光谱与参考光谱相似度计算;
计算方法为,步骤四计算的影像光谱吸收深度除以标准化后的影像光谱与去连续统后的参考光谱。计算公式如下所示。
步骤七,设定相似度阈值,认为大于阈值的像元即为目标矿物,大于阈值的像元的影像光谱吸收深度(步骤四计算结果)即为目标矿物的相对丰度。在试验区,识别出了高铝白云母(吸收位置为2200.8nm)、中铝白云母(吸收位置为2210.8nm)、低铝白云母(吸收位置为2220.8nm)、高岭石、埃洛石、地开石、明矾石、黄钾铁钒、方解石等八种矿物。
Claims (1)
1.一种批量化与自动化的高光谱遥感矿物填图方法,其特征在于:该方法具体步骤如下:
步骤一:光谱数据读入
这些光谱数据包括高光谱遥感数据与参考光谱数据,其中高光谱遥感数据是通过机载或星载高光谱传感器获取,参考光谱数据来自于标准矿物光谱库;
步骤二:光谱波段截取
其截取方法为,截取接近700—1300nm谱段的光谱识别褐铁矿、针铁矿、赤铁矿、黄钾铁钒,截取接近2100—2320nm谱段的光谱识别白云母、蒙脱石、高岭石、地开石、石膏、明矾石、黄钾铁钒,叶腊石,截取接近2200—2400m谱段的光谱识别绿泥石、绿帘石、角闪石、滑石、蛇纹石、方解石、白云石;
步骤三:光谱连续统去除
其实现方法为,光谱连续统去除即Continum Removal是一种用于分离光谱吸收特征的光谱分析方法,广泛应用高光谱遥感中去除背景影响并分离某些特定物质的吸收特征;连续统被定义为反射光谱曲线中反射峰之间的线性连接部分,连接折线在反射光谱峰值处的外角大于180度,光谱连续统去除即以反射光谱除以连续统光谱,如下公式所示:
其中,Rcr为连续统去除后的吸收光谱,R吸收光谱,Rc为吸收光谱的连续统;经过连续统去除之后,端点处反射率为1,端点之间反射率均小于1,影像光谱与参考光谱均需进行连续统去除;
步骤四:影像光谱吸收深度计算
其计算方法为,在经过波段截取、连续统去除的光谱上,寻找最小值,并利用数字“1”减去该“最小值”;
步骤五:影像光谱的标准化
其计算方法为,以参考光谱标准,对影像光谱进行标准化,以抑制因测试条件差异对光谱相似度计算造成的影响,标准化系数计算公式如下所示:
步骤六:影像光谱与参考光谱相似度计算
其计算方法为,步骤四计算的影像光谱吸收深度除以标准化后的影像光谱与去连续统后的参考光谱,计算公式如下所示:
步骤七:根据步骤五与步骤六的计算结果提取矿物异常;
设定相似度阈值,认为大于阈值的像元即为目标矿物,大于阈值的像元的影像光谱吸收深度即为目标矿物的相对丰度;若一个工作区有多景/航带数据,在镶嵌后设定阈值,也可之前设定阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310081966.8A CN103175801B (zh) | 2013-03-14 | 2013-03-14 | 一种批量化与自动化的高光谱遥感矿物填图方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310081966.8A CN103175801B (zh) | 2013-03-14 | 2013-03-14 | 一种批量化与自动化的高光谱遥感矿物填图方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103175801A true CN103175801A (zh) | 2013-06-26 |
CN103175801B CN103175801B (zh) | 2015-06-24 |
Family
ID=48635806
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310081966.8A Expired - Fee Related CN103175801B (zh) | 2013-03-14 | 2013-03-14 | 一种批量化与自动化的高光谱遥感矿物填图方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103175801B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103698312A (zh) * | 2013-12-31 | 2014-04-02 | 中国科学院地质与地球物理研究所 | 一种基于激光拉曼测定白云石有序度的方法 |
CN103954567A (zh) * | 2014-04-15 | 2014-07-30 | 塔里木大学 | 基于连续统去除法的土壤盐分测定方法 |
CN103984940A (zh) * | 2014-06-03 | 2014-08-13 | 核工业北京地质研究院 | 一种基于高光谱数据识别赤铁矿化的方法 |
CN104574283A (zh) * | 2013-10-23 | 2015-04-29 | 核工业北京地质研究院 | 一种用于叶蜡石信息提取的高光谱影像处理方法 |
CN104990880A (zh) * | 2015-07-09 | 2015-10-21 | 核工业北京地质研究院 | 基于航空高光谱遥感数据的绢云母矿物相对丰度计算方法 |
CN109580497A (zh) * | 2018-12-13 | 2019-04-05 | 中国国土资源航空物探遥感中心 | 一种基于奇异性理论的高光谱矿物异常信息提取方法 |
CN110702618A (zh) * | 2019-10-12 | 2020-01-17 | 东华理工大学 | 一种金属硫化物矿物显微高光谱图像的自动识别方法 |
CN111267984A (zh) * | 2020-01-21 | 2020-06-12 | 山东大学 | 基于高光谱技术分析的隧道内不良地质体识别系统与方法 |
CN111521579A (zh) * | 2020-05-07 | 2020-08-11 | 中国自然资源航空物探遥感中心 | 一种基于高光谱遥感数据的白云母及其光谱端元识别方法 |
WO2024077330A1 (en) * | 2022-10-11 | 2024-04-18 | Plotlogic Pty Ltd | Continuum estimation tool for mining of geological structures |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5606413A (en) * | 1995-01-19 | 1997-02-25 | Northrop Grumman Corporation | Real time spectroscopic imaging system and method |
CN101551471A (zh) * | 2009-05-19 | 2009-10-07 | 中国国土资源航空物探遥感中心 | 一种高光谱遥感矿物含量定量反演方法 |
-
2013
- 2013-03-14 CN CN201310081966.8A patent/CN103175801B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5606413A (en) * | 1995-01-19 | 1997-02-25 | Northrop Grumman Corporation | Real time spectroscopic imaging system and method |
CN101551471A (zh) * | 2009-05-19 | 2009-10-07 | 中国国土资源航空物探遥感中心 | 一种高光谱遥感矿物含量定量反演方法 |
Non-Patent Citations (3)
Title |
---|
YAN BOKUN ET AL.: ""Minerals mapping of the lunar surface with clementine UVVIS/NIR data based on spectra unmixing method and Hapke model"", 《ICARUS》, no. 208, 11 February 2010 (2010-02-11), pages 11 - 19 * |
刘圣伟等: ""成像光谱技术在典型蚀变矿物识别和填图中的应用"", 《中国地质》, vol. 33, no. 1, 28 February 2006 (2006-02-28), pages 178 - 186 * |
甘甫平等: "《高光谱遥感技术在地质领域中的应用》", 《国土资源遥感》, no. 4, 15 December 2007 (2007-12-15), pages 57 - 60 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104574283A (zh) * | 2013-10-23 | 2015-04-29 | 核工业北京地质研究院 | 一种用于叶蜡石信息提取的高光谱影像处理方法 |
CN104574283B (zh) * | 2013-10-23 | 2018-07-27 | 核工业北京地质研究院 | 一种用于叶蜡石信息提取的高光谱影像处理方法 |
CN103698312A (zh) * | 2013-12-31 | 2014-04-02 | 中国科学院地质与地球物理研究所 | 一种基于激光拉曼测定白云石有序度的方法 |
CN103954567A (zh) * | 2014-04-15 | 2014-07-30 | 塔里木大学 | 基于连续统去除法的土壤盐分测定方法 |
CN103954567B (zh) * | 2014-04-15 | 2016-12-07 | 塔里木大学 | 基于连续统去除法的土壤盐分测定方法 |
CN103984940B (zh) * | 2014-06-03 | 2017-12-26 | 核工业北京地质研究院 | 一种基于高光谱数据识别赤铁矿化的方法 |
CN103984940A (zh) * | 2014-06-03 | 2014-08-13 | 核工业北京地质研究院 | 一种基于高光谱数据识别赤铁矿化的方法 |
CN104990880A (zh) * | 2015-07-09 | 2015-10-21 | 核工业北京地质研究院 | 基于航空高光谱遥感数据的绢云母矿物相对丰度计算方法 |
CN104990880B (zh) * | 2015-07-09 | 2017-10-27 | 核工业北京地质研究院 | 基于航空高光谱遥感数据的绢云母矿物相对丰度计算方法 |
CN109580497A (zh) * | 2018-12-13 | 2019-04-05 | 中国国土资源航空物探遥感中心 | 一种基于奇异性理论的高光谱矿物异常信息提取方法 |
CN109580497B (zh) * | 2018-12-13 | 2020-01-03 | 中国自然资源航空物探遥感中心 | 一种基于奇异性理论的高光谱矿物异常信息提取方法 |
CN110702618A (zh) * | 2019-10-12 | 2020-01-17 | 东华理工大学 | 一种金属硫化物矿物显微高光谱图像的自动识别方法 |
CN111267984A (zh) * | 2020-01-21 | 2020-06-12 | 山东大学 | 基于高光谱技术分析的隧道内不良地质体识别系统与方法 |
CN111521579A (zh) * | 2020-05-07 | 2020-08-11 | 中国自然资源航空物探遥感中心 | 一种基于高光谱遥感数据的白云母及其光谱端元识别方法 |
CN111521579B (zh) * | 2020-05-07 | 2023-02-24 | 中国自然资源航空物探遥感中心 | 一种基于高光谱遥感数据的白云母及其光谱端元识别方法 |
WO2024077330A1 (en) * | 2022-10-11 | 2024-04-18 | Plotlogic Pty Ltd | Continuum estimation tool for mining of geological structures |
Also Published As
Publication number | Publication date |
---|---|
CN103175801B (zh) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103175801B (zh) | 一种批量化与自动化的高光谱遥感矿物填图方法 | |
Rogge et al. | Integration of spatial–spectral information for the improved extraction of endmembers | |
CN112990313B (zh) | 高光谱图像异常检测方法、装置、计算机设备及存储介质 | |
Koch et al. | Automated drill core mineralogical characterization method for texture classification and modal mineralogy estimation for geometallurgy | |
CN109580497B (zh) | 一种基于奇异性理论的高光谱矿物异常信息提取方法 | |
CN105069309A (zh) | 一种识别水文时间序列非线性趋势的方法 | |
CN103984940A (zh) | 一种基于高光谱数据识别赤铁矿化的方法 | |
CN101916440B (zh) | 基于数字图像形态学理论的高光谱异常检测方法 | |
CN109472287A (zh) | 基于二维Gabor小波的三维荧光光谱特征提取方法 | |
CN110188716A (zh) | 一种矿化信息提取方法及系统 | |
KR101750520B1 (ko) | 표적신호 분리를 이용하여 탐지성능이 향상된 초분광 영상의 표적물질 탐지방법 | |
CN110763913A (zh) | 一种基于信号分段分类的导数谱平滑处理方法 | |
CN112485203A (zh) | 一种基于高光谱成像分析的重金属污染分析方法 | |
CN115343685A (zh) | 应用于病害识别的多维度探地雷达检测方法、装置、设备 | |
Ishidoshiro et al. | Geological mapping by combining spectral unmixing and cluster analysis for hyperspectral data | |
CN114820351A (zh) | 高光谱噪声滤波方法、装置、设备和存储介质 | |
CN111898681B (zh) | 一种基于遥感云平台与决策树的火烧迹地精细提取方法 | |
CA3013020A1 (en) | Method and apparatus for non-intrusive program tracing with bandwidth reduction for embedded computing systems | |
CN104268896A (zh) | 基于光谱抽样直方图的超光谱降维匹配方法及系统 | |
CN102880861A (zh) | 基于线性预测倒谱系数的高光谱图像分类方法 | |
Chang et al. | Field-programmable gate array design of implementing simplex growing algorithm for hyperspectral endmember extraction | |
CN108007902B (zh) | 一种用高光谱数据计算白云母Al-OH吸收位置的方法 | |
CN111814727B (zh) | 一种基于高光谱的潜在指纹检测方法 | |
CN104036509A (zh) | 一种基于压缩感知的高光谱混合像元分解的方法 | |
Reddy et al. | Geospatial Approach for the Analysis of Forest Cover Change Detection using Machine Learning |
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 |
Granted publication date: 20150624 Termination date: 20170314 |
|
CF01 | Termination of patent right due to non-payment of annual fee |