[发明专利]一种提取岩芯氧化还原过渡带信息的光谱数据处理方法

说明书

本发明属于沉积型铀矿床相关信息提取技术领域，具体涉及一种适用于提取岩芯氧化还原过渡带信息的光谱数据处理方法。本发明包括如下步骤：步骤一、获取岩芯样品光谱数据，并按照自上而下的顺序记录光谱测量的位置；按照光谱自上而下的顺序对光谱数据时行编号；步骤二、对光谱数据进行预处理，进行大气校正，获取样品或地物光谱的反射率数据；步骤三、将批量光谱数据转换为高光谱栅格影像数据，每个像元对应一条光谱数据；步骤四、对高光谱数据进行重采样，提取波段在423nm，514nm，557nm，612nm，753nm，833nm，887nm，985nm，的影像；步骤五、对步骤四中的影像进行判断，确定铀矿的有利成矿带。本发明提高了数据处理的速度和信息提取的精度，减少了人为因素的影响。

技术领域

本发明属于沉积型铀矿床相关信息提取技术领域，具体涉及一种适用于提取岩芯氧化还原过渡带信息的光谱数据处理方法。

背景技术

氧化还原过渡带是沉积型铀矿勘查的重要目标。通过对光谱数据处理提取岩石中不同类型矿物的丰度，是目前矿物光谱分析常用的方法。当前的光谱处理方法主要是光谱全波段匹配或是部分连续波段的光谱匹配，具体算法有光谱角、混合解调滤波等，由于地表的物质组成很少是由单一矿物组成的，这些方法在信息提取的过程中易受其他矿物光谱或噪声的影响，提取信息精度较低。其次现有的光谱信息提取方法人工操作步骤多，增加了人为判断误差。因此，如何在提取岩芯氧化还原过渡带信息提取过程中减少其他地物或噪声的影响，人工操作误差等影响，成为当前光谱数据处理的前沿之一。

发明内容

本发明解决的技术问题：本发明提供一种适用于提取岩芯氧化还原过渡带信息的光谱数据处理方法，简化了数据处理的流程，可以批量处理光谱数据，集中提取目标信息，提高了数据处理的速度和信息提取的精度，减少了人为因素的影响。

本发明采用的技术方案：

一种提取岩芯氧化还原过渡带信息的光谱数据处理方法，包括如下步骤：步骤一、获取岩芯样品光谱数据，并按照自上而下的顺序记录光谱测量的位置；按照光谱自上而下的顺序对光谱数据时行编号；步骤二、对光谱数据进行预处理，进行大气校正，获取样品或地物光谱的反射率数据；步骤三、将批量光谱数据转换为高光谱栅格影像数据，每个像元对应一条光谱数据；步骤四、对高光谱数据进行重采样，提取波段在423nm，514nm，557nm，612nm，753nm，833nm，887nm，985nm，的影像；步骤五、对步骤四中的影像进行判断，确定铀矿的有利成矿带。

所述步骤五中的判断方法为：当423nm波段影像的像元值小于514nm波段影像相对应的像元值，判断该像元满足条件一，获得514nm波段影像像元值减去423nm波段影像的像元值的结果影像像元值H1；当887nm波段影像的像元值大于985nm波段影像相对应的像元值，判断该像元满足条件二，获得887nm波段影像像元值减去985nm波段影像的像元值的结果影像像元值H2；当557nm波段影像的像元值小于612nm波段影像相对应的像元值，判断该像元满足条件三，获得612nm波段影像像元值减去557nm波段影像的像元值的结果影像像元值H3；当753nm波段影像的像元值大于833nm波段影像相对应的像元值，判断该像元满足条件四。

选择同时满足条件一、二的像元，获得上述像元范围内的H1和H2相加的和H；选择同时满足条件三、四的像元，获得上述像元范围内的H3。

H的值代表像元二价铁的丰度值，H值越大，代表像元中二价铁的丰度越大，即含量相对越高。

H3的值代表像元三价铁的丰度值，H3值越大，代表像元中三价铁的丰度越大，即含量相对越高。

按照光谱记录自上而下的顺序，选择三价铁丰度减少，二价铁丰度增强的区域，为氧化还原过渡带，这些区域容易形成铀矿，为铀矿的有利成矿带。

提取波段的范围为423±5nm，514±5nm，557±5nm，612±5nm，753±5nm，833±5nm，887±5nm，985±5nm。