一种用高光谱数据计算白云母Al-OH吸收位置的方法技术

本发明专利技术属于高光谱遥感分析技术领域，具体涉及一种用高光谱数据计算白云母Al‑OH吸收位置的方法。包括以下步骤：步骤一、提取出白云母类矿物，结果记作result(Muscovite)；步骤二、对航空高光谱数据进行波段运算，将计算的结果记作result(Al‑OH)；步骤三、求步骤一和步骤二处理的结果的交集为result(mus‑wavelength)；步骤四、对步骤三中的结果进行密度分割。本发明专利技术针对现有光谱分辨率下的航空高光谱数据，提取白云母类矿物单元，通过波段运算计算出白云母类矿物的Al‑OH吸收峰位置，进一步根据吸收峰位置对白云母类矿物进行细分，从而更好地应用于矿产勘查。

【技术实现步骤摘要】
一种用高光谱数据计算白云母Al-OH吸收位置的方法
本专利技术属于高光谱遥感分析
，具体涉及一种用高光谱数据计算白云母Al-OH吸收位置的方法。
技术介绍
航空高光谱数据已经被广泛地应用于地表蚀变矿物填图，目前，人们主要通过航空高光谱成像设备(如澳大利亚的Hymap，加拿大的CASI/SASI/TASI等)获取高光谱数据，并采用一系列技术手段对这些数据进行蚀变矿物信息提取，这些技术手段主要是基于对高光谱图像反射率光谱与光谱库或者端元光谱的特征对比来实现的，可以提取蛇纹石、透闪石、褐铁矿、白云母、蒙脱石、高岭石、绿泥石、绿帘石、方解石等数十种蚀变矿物，识别准确度较高，因而获得较为广泛的应用。在大部分蚀变矿物中，白云母类矿物是分布较为广泛同时也是与成矿关系较为密切的一类蚀变矿物，另一方面，对蚀变矿物的光谱研究发现，白云母类矿物的Al-OH吸收峰位置(位于2200nm附近)与此类矿物中化学组分具有较为密切的联系，吸收峰位置的微量变化常常具有重要的地质指示意义，越来越多的人开始将白云母矿物的Al-OH吸收峰位置值应用于矿产资源潜力评价、成矿流体路径反演和变质岩变质程度上，因而，准确的获取白云母类矿物的Al-OH吸收峰位置值则显得较为重要。传统的高光谱蚀变矿物提取方法可以较为准确地提取白云母类矿物，但受航空高光谱成像设备光谱分辨率(如SASI的光谱分辨率为15nm)的制约，很难对分辨不同白云母类矿物的之间微小变化，如不同Al-OH吸收峰位置(2180nm～2230nm)的白云母类矿物，因而大大制约了航空高光谱数据的进一步应用。因此，如何在现有的航空高光谱数据条件下利用白云母类矿物在不同波段的反射率等参数来计算白云母的Al-OH吸收位置，进而对白云母类矿物进行细分，是本专利技术重点解决的问题。
技术实现思路
本专利技术要解决的技术问题是提供一种用高光谱数据计算白云母Al-OH吸收位置的方法，针对现有光谱分辨率下的航空高光谱数据，提取白云母类矿物单元，通过波段运算计算出白云母类矿物的Al-OH吸收峰位置，进一步根据吸收峰位置对白云母类矿物进行细分，从而更好地应用于矿产勘查。为解决上述技术问题，本专利技术一种用高光谱数据计算白云母Al-OH吸收位置的方法，包括以下步骤：步骤一、利用航空高光谱数据提取白云母类矿物，读取经过大气校正的影像文件，提取出白云母类矿物，结果记作result(Muscovite)，以图像存储，含白云母类矿物的像元值记作1，不含白云母类矿物的像元值记作0；步骤二、对航空高光谱数据进行波段运算，将计算的结果记作result(Al-OH)；步骤三、求步骤一和步骤二获得的结果的交集；求步骤一和步骤二处理的结果的交集为result(mus-wavelength)，将运算结果作为新的像元值在新图像上存储；步骤四、对步骤三中的结果进行密度分割，提高图像的可视化程度，对结果图像进行密度分割处理，0～2180之间的像元统一赋上黑色，而在2180～2230之间，可根据像元值范围进一步细分，进而提高图像的可视化程度，即得到所需的白云母类矿物Al-OH吸收位置图像。所述步骤二中，选择航空高光谱数据中的波段进行运算，主要选择2195nm、2210nm、2225nm三个波段参与运算，结果result(Al-OH)＝2286.96-36.30*b2/b1-39.49*b3/b2，其中b1＝像元光谱在2195nm的反射率值，b2＝像元光谱在2210nm附近的反射率值，b3＝像元光谱在2225nm附近的反射率值，将计算的结果记作result(Al-OH)。所述的步骤一中，读取经过大气校正的影像文件，分别进行最小噪声分离(MNF)、数据维数判断、获取端元光谱、计算纯净像元指数、N维可视化和端元光谱选择、波谱识别的一系列运算过程，得到白云母类矿物。本专利技术的有益技术效果在于：受航空高光谱成像设备的光谱分辨率制约，传统的航空高光谱数据处理方法提取的白云母类矿物单元较为简单，通常根据铝羟基吸收峰位置分为高铝白云母和低铝白云母，而本专利技术通过波段运算对白云母类矿物的铝羟基吸收峰位置进行计算，获得的白云母类矿物Al-OH羟基吸收峰位置分辨率可达1nm，据此可对白云母类矿物单元进行更细致的划分，本专利技术可应用于航空高光谱数据处理领域；本专利技术通过利用航空高光谱数据不同波段的反射率比值进行线性拟合计算白云母类矿物Al-OH吸收峰位置值。附图说明图1为本专利技术一种用高光谱数据计算白云母Al-OH吸收位置的方法流程图。具体实施方式下面结合附图和实施例对本专利技术作进一步详细说明。本专利技术一种用高光谱数据计算白云母Al-OH吸收位置的方法，其具体实现步骤如下：步骤一、利用航空高光谱数据提取白云母类矿物，读取经过大气校正的影像文件，分别进行最小噪声分离(MNF)、数据维数判断、获取端元光谱、计算纯净像元指数、N维可视化和端元光谱选择、波谱识别的一系列运算过程，提取出白云母类矿物，结果记作result(Muscovite)，以图像存储，含白云母类矿物的像元值记作1，不含白云母类矿物的像元值记作0；步骤二、对航空高光谱数据进行波段运算；选择航空高光谱数据中的波段进行运算，主要选择2195nm、2210nm、2225nm三个波段参与运算，结果result(Al-OH)＝2286.96-36.30*b2/b1-39.49*b3/b2，其中b1＝像元光谱在2195nm的反射率值，b2＝像元光谱在2210nm附近的反射率值，b3＝像元光谱在2225nm附近的反射率值，将计算的结果记作result(Al-OH)，以图像形式存储，像元值即为表达式计算的结果值；步骤三、求步骤一和步骤二获得的结果的交集；求步骤一和步骤二处理的结果的交集，表达式为result(mus-wavelength)＝result(Muscovite)*result(Al-OH)，将运算结果作为新的像元值在新图像上存储；步骤四、对步骤三中的结果进行密度分割，提高图像的可视化程度，由于最终的图像像元值存在0值和2180～2230之间的值，而我们关注的是2180～2230之间的像元，因此，对结果图像进行密度分割处理，0～2180之间的像元统一赋上黑色，而在2180～2230之间，可根据像元值范围进一步细分，进而提高图像的可视化程度，即得到所需的白云母类矿物Al-OH吸收位置图像。本文档来自技高网...

【技术保护点】
一种用高光谱数据计算白云母Al‑OH吸收位置的方法，其特征在于：包括以下步骤：步骤一、利用航空高光谱数据提取白云母类矿物，读取经过大气校正的影像文件，提取出白云母类矿物，结果记作result(Muscovite)，以图像存储，含白云母类矿物的像元值记作1，不含白云母类矿物的像元值记作0；步骤二、对航空高光谱数据进行波段运算，将计算的结果记作result(Al‑OH)；步骤三、求步骤一和步骤二获得的结果的交集；求步骤一和步骤二处理的结果的交集为result(mus‑wavelength)，将运算结果作为新的像元值在新图像上存储；步骤四、对步骤三中的结果进行密度分割，提高图像的可视化程度，对结果图像进行密度分割处理，0～2180之间的像元统一赋上黑色，而在2180～2230之间，可根据像元值范围进一步细分，进而提高图像的可视化程度，即得到所需的白云母类矿物Al‑OH吸收位置图像。

【技术特征摘要】
1.一种用高光谱数据计算白云母Al-OH吸收位置的方法，其特征在于：包括以下步骤：步骤一、利用航空高光谱数据提取白云母类矿物，读取经过大气校正的影像文件，提取出白云母类矿物，结果记作result(Muscovite)，以图像存储，含白云母类矿物的像元值记作1，不含白云母类矿物的像元值记作0；步骤二、对航空高光谱数据进行波段运算，将计算的结果记作result(Al-OH)；步骤三、求步骤一和步骤二获得的结果的交集；求步骤一和步骤二处理的结果的交集为result(mus-wavelength)，将运算结果作为新的像元值在新图像上存储；步骤四、对步骤三中的结果进行密度分割，提高图像的可视化程度，对结果图像进行密度分割处理，0～2180之间的像元统一赋上黑色，而在2180～2230之间，可根据像元值范围进一步细分，进而提高图像的可视化程度，即得到所需的白云母类矿物Al-...

【专利技术属性】
技术研发人员：孟树，叶发旺，张川，刘洪成，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京,11