[发明专利]一种高光谱金矿化蚀变远景区快速识别方法

说明书

本发明属于遥感地质调查技术领域，具体涉及一种高光谱金矿化蚀变远景区快速识别方法包括：步骤一：对获取的研究区航天/航空高光谱数据产品开展预处理，得到浮点型高光谱反射率数据；步骤二：利用高光谱反射率数进行光谱增殖处理，得到新高光谱数据；步骤三：利用新高光谱数据，逐像元获取并判别吸收峰和反射峰位置，识别出绢云母矿物并厘定为短波云母亚类；步骤四：采用基于动态位置的解析几何方法计算吸收深度；计算吸收峰bc处的曲率Curvsubgt;C/subgt;；步骤五：选择吸收深度大且曲率大的蚀变短波云母亚类，设置合适阈值进行高端密度切割，得到短波云母蚀变分布图；结合地质综合分析圈定金矿化有利地段，在野外查证基础上优选出金矿化蚀变远景区。

技术领域

本发明属于遥感地质调查技术领域，具体涉及一种高光谱金矿化蚀变远景区快速识别方法。

背景技术

金矿属于贵金属矿产，是关系国民经济运行的重要矿产资源。金矿化的地质过程通常与热液活关系密切，普遍发育远大于矿体平面面积的蚀变带，因而可以被遥感技术探测出来，为圈定找矿靶区提供重要依据。作为遥感前沿之一的高光谱遥感技术，可一次性获取上百个波段的光谱曲线，从而有效的识别地物类型。当前，无论是卫星高光谱数据还是航空高光谱数据，原始数据的波段采样间隔均为6nm～20nm，地物光谱曲线呈折线状，尽管可以有效识别出云母大类矿物，但在云母亚类矿物的详细划分上多采用目视识别方法，精细定量化数值计算工作开展有限，难以更加准确的指导找矿工作。

在高光谱数据处理方面，众多研究集中在光谱相似度填图方法研发和应用方面，在光谱特征参量填图方法方面工作相对较少。由上百个波段组成的高光谱数据必然存在信息冗余，而且针对特定矿物填图时存在很多无光谱特征的波段，如果不加选择的全部纳入处理将增大数据量，增长运算时间。

因此，需要针对上述技术问题设计一种利用高光谱数据处理技术对金矿化蚀变远景区快速识别的方法，以解决上述技术问题。

发明内容

本发明针对上述现有技术的不足，提供一种高光谱金矿化蚀变远景区快速识别方法，用于解决原有技术中高光谱数据处理无法得到精准地理位置，光谱特征参量填图方法较少的技术问题。

本发明的技术方案：

一种高光谱金矿化蚀变远景区快速识别方法，包括以下步骤：

步骤一：对获取的研究区的航天/航空高光谱数据产品开展预处理，得到浮点型高光谱反射率数据；选择包含2130nm～2380nm谱段范围的波段提取绢云母矿物，得到特征优选后的高光谱反射率数据；

步骤二：利用步骤一得到的高光谱反射率数，采用非扭结边界条件的三次样条插值法进行光谱增殖处理，得到波段间隔为1nm的新高光谱数据；

步骤三：利用步骤二得到的新高光谱数据，逐像元获取并判别吸收峰和反射峰位置，识别出绢云母矿物并将吸收峰的波长位置小于2205nm的像元厘定为短波云母亚类；

步骤四：利用步骤三得到的短波绢云母的list_{abs_peak}和list_{ref_peak}，采用基于动态位置的解析几何方法计算吸收深度；采用最小二乘法利用抛物线方程拟合吸收峰附近的21个光谱点，得到拟合后的圆滑曲线方程并计算吸收峰bc处的曲率Curv_C；

步骤五：利用步骤四得到的数据，选择吸收深度大且曲率大的蚀变短波云母亚类，设置合适的阈值进行高端密度切割，得到短波云母蚀变分布图；结合地质综合分析，圈定分布在成矿有利地质体的合适构造部位的短波云母与褐铁矿叠加分布的地段为金矿化有利地段，在野外查证基础上优选出金矿化蚀变远景区。

所述步骤一还包括：对获取的研究区的航天/航空高光谱数据产品开展几何校正，在相对高差大的情况下引入DEM高程数据进行正射校正，在此基础上进行辐射定标，得到具有地理投影坐标的高光谱辐亮度数据；