一种绿泥石铁镁含量反演方法、系统、介质、设备及终端

本发明属于地质矿产资源遥感识别，尤其涉及一种绿泥石铁镁含量反演方法、系统、介质、设备及终端。

背景技术：

1、目前，高分五号(gf-5)卫星是世界首颗实现对大气和陆地综合观测的全谱段高光谱卫星，可以获取从390nm到2500nm共330个谱段的高光谱数据，设计寿命8年，2019年3月21日正式投入使用。gf-5首次搭载了大气痕量气体差分吸收光谱仪(emi)、主要温室气体探测仪(gmi)、大气多角度偏振探测仪(dpc)、大气环境红外甚高分辨率探测仪(aius)、可见短波红外高光谱相机(ahsi)和全谱段光谱成像仪(vims)共6台载荷。gf-5卫星对地面物质成分的探测十分精确，可对内陆水体、路表生态环境、蚀变矿物、岩矿类别进行有效探测，为环境监测、资源勘察、防灾减灾等行业提供高质量、高可靠的高光谱数据。全谱段光谱成像仪具有谱段范围宽、空间分辨率高、辐射定标精度高的技术特点。

2、绿泥石是广泛存在于自然界中低级变质作用和热液蚀变作用的产物，常与斑岩铜矿、热液蚀变型金矿、砂岩型铀矿等成矿作用密切相关，研究其在不同地质环境中化学成分的变化，可反映出其形成时的物理化学条件及其光谱特征的变化。绿泥石的晶体化学式一般可表示为y(z4o10)(oh)2y3(oh)6，其中y代表mg2+、fe2+、al3+、fe3+、cr、ni、mn，z代表si、al、fe3+。一般认为绿泥石主要是由热液流体中的铁镁质矿物直接结晶或者蚀变交代围岩又或者沿裂隙沉淀所致。事实上，绿泥石矿物具有独特的诊断性吸收谱带，在特定特征谱带中能够指示离子类矿物、单矿物的存在，通过这些谱带或其组合，能够实现对绿泥石铁镁含量的反演。在不同的成矿环境中，绿泥石的分子晶型和化学成分通常会发生规律性的变化，其化学成分中铁镁含量的特征吸收谱带可以作为指示绿泥石矿物的存在，有着重要的标型意义。绿泥石的主要成分是富铁绿泥石和富镁绿泥石，越靠近热液中心的绿泥石成分上更加富铁，利用两种绿泥石的特征吸收峰位分别不同，以此来指示热液中心，对进一步的矿床勘查工作提供积极的指导。

3、高光谱矿物识别物理方法主要有两大类，分别是以重建光谱与标准光谱相似性度量为基础的光谱匹配方法和以光谱吸收谱带参量为基础的模式识别方法。光谱匹配对矿物光谱的微小差异不够敏感，容易受地形、背景等外界干扰，无法突出光谱谷谱型在相似性计算中重要性，且受图像的信噪比、光谱定标和光谱重建精度等因素的影响比较大。甘甫平等在总结分析常见光谱匹配法的矿物识别能力基础上，提出了一种光谱特征增强的匹配度算法，该算法较其他方法识别矿物信息的检出限更低、正确率更高。但其主要还是对矿物大类或光谱特征具有唯一性时，应用效果较好。甘甫平等分析了矿物光谱的变化规律，评价了光谱参量的敏感性与稳定性，提出了建立矿物识别谱系的思想，参照或借鉴矿物学的分类方法，在可见短波红外光谱区间，分别以主要吸收谱带、谱带组合特征、谱带精细特征和谱带变异特征为基础，对矿物进行“类－族-种－亚种”逐层识别，构成判别决策过程的树状结构。董新丰等依据了甘甫平等提出的矿物分层谱系识别思想，采用了综合光谱特征增强匹配度和特征参量开展了gf-5矿物精细识别应用研究，但仍需提高对gf-5光谱特征理解，进而实现对亚像元、弱信息的识别。绿泥石矿物光谱具有fe-oh和mg-oh的诊断吸收特征，波长位置分别为2245～2260和2320～2360nm，梁树能等对测量的绿泥石岩矿光谱进行分析，获取了岩石样品中绿泥石矿物的特征光谱参量，但并未将此结果结合高分五号短波红外高光谱数据当中，并进一步获取数据中相关化学参数。

4、前人提出了大量基于多光谱遥感数据(landsat tm/etm、worldview-2和aster数据)的矿产快速勘查方法，这些方法主要是针对地质构造中的控矿构造进行解译和图像变换处理，进而突出有利于矿床产出的构造部位。这些方法可以有效识别矿物种类差异大的岩石和明显的地质构造。但是对于矿物亚类的探测和细分效果不佳，这些方法没有专门针对与成矿相关的蚀变矿物的遥感光谱响应提出有效的处理方法。目前多光谱数据源的多波段数据的空间分辨率普遍低，高空间分辨率数据波段少，信噪比较差，有用信息被大量无用信息混杂，严重影响勘查精度。

5、在斑岩铜矿、热液蚀变型金矿、砂岩型铀矿等矿产中，矿体外围富含指示性的绿泥石矿物，同时矿物随分布地段的变化也表现出亚类的变化，这些矿物的亚类变化在国产高分五号数据的短波红外波段具有特征谱带而可以被识别。长期以来，遥感技术领域对矿床的指示矿物亚类变化规律的探测相对比较陌生，针对矿床的绿泥石指示矿物进行遥感识别的尝试更少。利用遥感技术探测构造控矿型的矿产的关键在于全面解译区内的断裂构造，图像变换和处理的方式比较简单直观。但对于挖掘矿床不同部位指示性蚀变矿物亚类的变化，这种方法难以发挥作用。现有技术中没有专门针对与成矿相关的绿泥石指示矿物的遥感光谱响应提出有效的处理方法。

6、目前的勘测方法针对构造容矿型的矿床比较有效，但是对于由于隐伏矿床在地表的指示信息较少，传统的遥感解译和图像变换方法对隐伏矿无效。对于隐伏型矿物，矿物、岩石的光谱特征与其物理化学属性之间的关联分析是高光谱遥感提取岩矿信息的基础。岩石矿物单个诊断性吸收特征可以用吸收位置、吸收深度、吸收宽度、面积、吸收对称性、吸收峰数目和排列次序作完整的表征。在不同的成矿环境中，绿泥石的分子晶型和化学成分通常会发生规律性的变化，其化学成分中铁镁含量的特征吸收谱带可以作为指示绿泥石矿物的存在，有着重要的标型意义。绿泥石的主要成分是富铁绿泥石和富镁绿泥石，越靠近热液中心的绿泥石成分上更加富铁，利用两种绿泥石的特征吸收峰位分别不同，以此来指示热液中心，对进一步的矿床勘查工作提供积极的指导。

7、通过上述分析，现有技术存在的问题及缺陷为：

8、(1)现有高光谱矿物识别物理方法中，光谱匹配对矿物光谱的微小差异不够敏感，容易受地形、背景等外界干扰，无法突出光谱谷谱型在相似性计算中重要性，且受图像的信噪比、光谱定标和光谱重建精度等因素的影响较大。

9、(2)现有基于多光谱遥感数据的矿产勘查方法对于矿物亚类的探测和细分效果不佳，且没有专门针对蚀变矿物的遥感光谱响应提出有效的处理方法。

10、(3)目前多光谱数据源的多波段数据空间分辨率普遍低，高空间分辨率数据波段少，信噪比较差，有用信息被大量无用信息混杂，严重影响勘查精度。

11、(4)目前的勘测方法针对构造容矿型的矿床比较有效，但是对于由于隐伏矿床在地表的指示信息较少，传统的遥感解译和图像变换方法对隐伏矿无效。

技术实现思路

1、针对现有技术存在的问题，本发明提供了一种绿泥石铁镁含量反演方法、系统、介质、设备及终端，尤其涉及一种基于高分五号卫星高光谱数据的绿泥石铁镁含量反演方法、系统、介质、设备及终端。

2、本发明是这样实现的，一种绿泥石铁镁含量反演方法，绿泥石铁镁含量反演方法包括：获取高分五号卫星的高光谱数据集并进行校正运算，得到反映地表反射率值的图像；对反映地表反射率值的图像进行光谱波段裁剪，得到短波红外高光谱数据集并进行savitzky-golay滤波平滑处理；确定参考光谱并与经过滤波平滑处理的短波红外高光谱数据集进行匹配，得到匹配值图像并进行统计分析；对滤波后的短波红外高光谱数据集进行掩膜处理，输出绿泥石矿物分布较纯的区段；对短波红外高光谱数据集进行一阶导数运算，输出光谱谱峰位置图像；计算光谱谱峰位置图像各个像元绿泥石的值，并输出值图像。

3、进一步，绿泥石铁镁含量反演方法包括以下步骤：

4、步骤一，获取高分五号卫星的高光谱数据集，利用偏置和增益数据对地表数字量化值进行辐射校正，并得到卫星视角下观测的地表辐射亮度值；根据大气传输模型对辐射亮度图像进行大气校正，得到反映地表反射率值的图像；

5、步骤二，将大气校正后的反射率影像进行第286～320波段的光谱波段裁剪，得到2142.31～2428.98nm波长范围的短波红外高光谱数据集；对裁剪得到的高分五号短波红外高光谱数据集进行savitzky-golay滤波，削弱光谱数据的噪声；

6、步骤三，在usgs光谱库中选择绿泥石矿物的标准光谱并裁剪至2142.31～2428.98nm波长范围作为参考光谱；采用波谱特征拟合方法对经过滤波平滑的高分五号短波红外高光谱数据集进行匹配，得到匹配值图像；

7、步骤四，对匹配值图像进行统计分析，取均值至均值与方差之和的范围生成掩膜图像；依据掩膜图像对savitzky-golay滤波后2142.31～2428.98nm波长范围的短波红外高光谱数据集进行掩膜处理，输出绿泥石矿物分布较纯的区段；

8、步骤五，选择2184.59～2386.84nm波段作为寻峰范围，对高分五号短波红外高光谱数据集逐个像元进行一阶导数运算；选择导数为0的波段波长值存入图像像元响应的位置，输出光谱谱峰位置图像；

9、步骤六，依据值与2184.59～2386.84nm波长范围内fe-oh吸收峰位置的关系，计算光谱谱峰位置图像各个像元绿泥石的值，并输出值图像。

10、进一步，步骤一中的辐射校正模型为：

11、l＝gain×dn+offset；

12、式中，l为辐射亮度值，gain为增益，offset为偏置，dn为像元亮度值；

13、大气校正模型为：

14、

15、式中，l为传感器像元接收到的总辐射亮度；ρ为像元表面反射率；ρe为像素周围平均表面反射率；s为大气球面反照率；la为大气后向散射辐射率；a和b为取决于大气条件和几何条件的两个系数。

16、进一步，步骤二中的将大气校正后的反射率影像进行第286～320波段的光谱波段裁剪，得到2142.31～2428.98nm波长范围的短波红外高光谱数据集包括：

17、对大气校正后的反射率影像进行光谱波段裁剪，对高分五号高光谱影像数据集390～2500nm波长范围的330个波段，选择第286～320波段的光谱波段进行光谱裁剪，得到2142.31～2428.98nm波长范围的短波红外高光谱数据集；

18、其中，savitzky-golay滤波算法公式为：

19、

20、式中，hi/h为平滑系数，由最小二乘法拟合多项式求得。

21、进一步，步骤三中地在usgs光谱库中选择绿泥石矿物的标准光谱并裁剪至2142.31～2428.98nm波长范围并作为参考光谱；采用波谱特征拟合方法对经过滤波平滑的高分五号短波红外高光谱数据集进行匹配，得到匹配值图像包括：

22、(1)在usgs光谱库中选择绿泥石矿物的标准光谱并裁剪至2142.31～2428.98nm波长范围，并作为影像参考光谱；

23、(2)对处理后的高分五号短波红外高光谱数据集以及usgs光谱库中筛选和处理好的参考光谱进行包络线校正处理，使两者数据归一化，公式如下：

24、

25、

26、式中，rcj为波段j的包络线去除值，rj为波段j的原始光谱反射率，k为吸收曲线中起始点波段和末端点波段间的斜率，λj为第j波段，rend和rstand为吸收曲线中起始点和末端点的原始光谱反射率，λend和λstand为吸收曲线中起始点和末端点的波长；

27、(3)采用光谱特征拟合方法对经过滤波平滑和包络线校正的高分五号短波红外高光谱数据集进行匹配，得到匹配值图像；通过对比增加常量k调整参考光谱的曲线外形，进而实现像元光谱拟合匹配；拟合图像灰度值越大代表匹配度越高，成为目标矿物的可能性越大，计算公式如下：

28、

29、式中，为经过变换的参考光谱特征曲线，p为参考光谱的原光谱特征曲线，和p均经过包络线去除处理，k为常数且k不等于-1。

30、进一步，步骤四中的对匹配值图像进行统计分析，取均值至均值与方差之和的范围生成掩膜图像；依据掩膜图像对savitzky-golay滤波后2142.31～2428.98nm波长范围的短波红外高光谱数据集进行掩膜处理包括：

31、(1)对匹配值图像进行统计分析，取像元统计值均值至均值与方差之和的范围生成掩膜图像；

32、(2)依据掩膜图像对savitzky-golay滤波后2142.31～2428.98nm波长范围的高分五号短波红外高光谱数据集进行掩膜处理，输出绿泥石矿物分布较纯的区段。

33、进一步，步骤五中的选择2184.59～2386.84nm波段作为寻峰范围，对高分五号短波红外高光谱数据集逐个像元进行一阶导数运算；选择导数为0的波段波长值存入图像像元响应的位置，输出光谱谱峰位置图像包括：

34、(1)获取经掩膜处理得到的反映绿泥石矿物分布较纯区段的高分五号短波红外高光谱数据集图像；

35、(2)选择2184.59～2386.84nm波段作为寻峰范围，对逐个像元进行一阶导数运算；选择导数为0的波段波长值存入图像像元响应的位置，输出反映光谱谱峰位置图像，计算公式如下：

36、

37、式中，δλ＝λj-λi，λj>λi。

38、进一步，步骤六中的值与2184.59～2386.84nm波长范围内fe-oh吸收峰位置的关系公式为：

39、y＝102.63-0.023x；

40、式中，y为值；x为fe-oh吸收峰位置，单位cm-1。

41、本发明的另一目的在于提供一种应用所述的绿泥石铁镁含量反演方法的绿泥石铁镁含量反演系统，绿泥石铁镁含量反演系统包括：

42、高光谱数据集获取模块，用于获取高分五号卫星的高光谱数据集并进行辐射定标和大气校正运算，得到反映地表反射率值的图像；

43、图像裁剪滤波模块，用于对反映地表反射率值的图像进行光谱波段裁剪，得到高分五号短波红外高光谱数据集并进行savitzky-golay滤波平滑处理；

44、图像匹配模块，用于确定参考光谱并与经过滤波平滑处理的高分五号短波红外高光谱数据集进行匹配，得到匹配值图像并进行统计分析；

45、数据掩膜处理模块，用于对滤波后的高分五号短波红外高光谱数据集进行掩膜处理，输出绿泥石矿物分布较纯的区段；

46、铁镁含量反演模块，用于对短波红外高光谱数据集进行一阶导数运算，输出光谱谱峰位置图像并计算各个像元绿泥石的值，并输出值图像。

47、本发明的另一目的在于提供一种计算机设备，计算机设备包括存储器和处理器，存储器存储有计算机程序，计算机程序被处理器执行时，使得处理器执行所述的绿泥石铁镁含量反演方法的步骤。

48、本发明的另一目的在于提供一种计算机可读存储介质，存储有计算机程序，计算机程序被处理器执行时，使得处理器执行所述的绿泥石铁镁含量反演方法的步骤。

49、本发明的另一目的在于提供一种信息数据处理终端，信息数据处理终端用于实现所述的绿泥石铁镁含量反演系统。

50、结合上述的技术方案和解决的技术问题，本发明所要保护的技术方案所具备的优点及积极效果为：

51、第一，本发明提供的绿泥石铁镁含量反演方法，是通过获取高分五号卫星的高光谱数据集并进行辐射定标和大气校正，得到反映地表反射率值的图像并进行光谱波段裁剪，得到短波红外高光谱数据集并进行savitzky-golay滤波平滑处理；确定参考光谱并与数据集进行匹配，得到匹配值图像；对滤波后的数据集进行掩膜处理，输出绿泥石矿物分布较纯的区段；对数据集进行一阶导数运算，输出光谱谱峰位置图像并计算各个像元绿泥石的值，输出值图像实现的。本发明能够在较短时间内高效、快速、准确地识别绿泥石矿物富镁-富铁亚类的分布范围，解决了现有技术中蚀变矿物亚类识别的难点问题。

52、与现有技术相比，本发明还具有如下技术效果：

53、(1)本发明是专门针对蚀变绿泥石矿物提出的遥感识别方法，可应用于高寒、高海拔以及人迹罕至的荒漠戈壁的基岩裸露地区的金属矿勘查，这些地区重型勘查设备难以抵达、传统勘查成本很高。因此，本发明对降低这类地区的蚀变矿产找矿勘查成本极为重要。

54、(2)本发明提供的反演方法通过获取与绿泥石蚀变矿物密切相关的矿物高光谱特征，能够在较短时间内高效、快速、准确地识别富铁、富镁绿泥石的分布范围，解决了蚀变矿物亚类识别的难点问题。

55、(3)本发明获取高分五号高光谱影像；高分五号遥感影像为高光谱卫星遥感影像，具有330个波段，尤其是2184.59～2386.84nm波长范围的短波红外部分，反映地表物质的光谱特征更准确，结合fe/(fe+mg)值与2184.59～2386.84nm波长范围内fe-oh吸收峰位置的关系，可以快速区分绿泥石蚀变矿物的亚类。

56、第二，目前应用的各类多光谱数据源，有的空间分辨率低，有的则波段较少，有用信息被大量噪声信息混杂，严重影响勘查的精度。卫星遥感技术的发展和成熟为开发遥感识别方法提供了良好的基础，我国现有的高分五号卫星具有光谱分辨率高、重访周期短、免费使用的优点，因此高分五号高光谱数据可以为矿产资源勘查提供有效的技术方法支撑。高分五号卫星对地面物质成分的探测十分精确，可对内陆水体、路表生态环境、蚀变矿物、岩矿类别进行有效探测，为环境监测、资源勘察、防灾减灾等行业提供高质量、高可靠的高光谱数据。全谱段光谱成像仪具有谱段范围宽、空间分辨率高、辐射定标精度高的技术特点。高光谱遥感识别矿物主要依赖于矿物成分的吸收特征，岩石光谱本质上是矿物的混合光谱，其光谱特征与成分、结构、构造、表面状态等因素有关。本发明基于高分五号卫星短波红外高光谱，使用短波红外数据集进行校正运算、波段裁剪、特征谱带匹配等，能够在人迹罕至的高山高原荒漠区域内高效、快速、准确地识别与矿床相关的绿泥石矿物富镁－富铁亚类分布。本发明是针对蚀变矿物提出的遥感识别方法，可应用于高寒、高海拔以及人迹罕至的荒漠戈壁的基岩裸露地区的金属矿勘查，这些地区重型勘查设备难以抵达、传统勘查成本很高。因此，本发明通过对降低这类地区的蚀变矿产找矿密切相关的矿物高光谱特征，能够在短时间内获取绿泥石矿物亚类分布范围，解决了蚀变矿物亚类识别的难点问题。本发明提供的方法前人尚未提出，本发明提出的绿泥石铁镁含量反演方法能够填补相关领域的空白，对指导找矿工作具有一定意义。

57、第三，作为本发明的权利要求的创造性辅助证据，还体现在以下几个重要方面：

58、本发明的技术方案解决了人们一直渴望解决、但始终未能获得成功的技术难题：

59、通过各种蚀变矿物寻找隐伏型矿床对于国民经济发展和国防安全具有重要意义。新疆、西藏、青海、甘肃等地区矿产资源丰富，但是自然环境恶劣、交通条件极差、地形起伏大，导致该地区地质工作程度低，矿产勘查工作程度低，效果差。目前应用的各类多光谱数据源，有的空间分辨率低，有的则波段较少，有用信息被大量噪声信息混杂，严重影像勘查的精度。因此，亟须一种方法能够在人迹罕至的高原、高山区内高效、快速、准确地识别绿泥石矿物亚类分布规律。矿物识别是高光谱遥感技术优势之一，已在地质矿产领域取得了显著应用效果。随着光谱分辨率的逐渐提高，高光谱遥感应用于矿物识别逐渐从识别矿物种类向矿物亚类、矿物成分、相关化学参数等精细信息的识别方向发展，且随着应用实践的深入，对蚀变矿物的信息需求也越来越大。高分五号(gf-5)卫星的可见短波红外高光谱相机为绿泥石矿物铁镁含量的反演提供了可能。

60、为了探索矿物精细光谱特征与其矿物化学参数间的内在关联性，为利用高光谱技术探测矿物微观成分信息，本研究通过高分五号卫星短波红外高光谱数据中的绿泥石矿物光谱特征不同实现来探讨富铁、富镁绿泥石的分布范围，以弥补多光谱或全色成像遥感定量应用的局限和不足。本发明提供的技术方案是在以高分五号高光谱影像为基础，依据usgs光谱库标准谱识别绿泥石信息，探测绿泥石分布区域，并结合值与2184.59～2386.84nm波长范围内fe-oh吸收峰位置的关系，进而计算光谱谱峰位置图像各个像元绿泥石的值，并输出值图像，然后根据值识别富镁、富铁绿泥石分布地段，最后通过野外调查验证，寻找富镁、富铁绿泥石，不仅避免了仅仅使用landsat、aster等多光谱遥感数据，导致光谱分辨率低，异常识别不够准确，且混杂噪声，严重影响提取精度的问题，而且还实现高海拔、艰险地区蚀变矿床绿泥石辅助的快速找矿，为矿产勘查工作服务，提高工作效率，节约国家和社会投入的资金。所以本发明提供的方法或者系统，能够在较短时间内高效、快速、准确地识别富镁、富铁绿泥石蚀变矿物分布范围，且前人尚未提出，本发明提供的绿泥石铁镁含量反演方法能够填补相关领域的空白，对指导找矿工作具有一定意义。