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

权利要求书

1.一种高光谱金矿化蚀变远景区快速识别方法，其特征在于：包括以下步骤：

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

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

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

所述步骤三还包括：对于二维空间上任一像元(j+1，i+1)，设其吸收峰波长为(b_{abs_peak})nm，反射峰波长为(b_{ref_peak})nm，则吸收峰对应的数值为super_data[b_{abs_peak}，j，i]，反射峰对应的数值为super_data[b_{ref_peak}，j，i]；吸收峰同时满足下面两个条件：

条件一：super_data[b_{abs_peak}，j，i]＜super_data[(b_{abs_peak}+1)，j，i]

条件二：super_data[b_{abs_peak}，j，i]＜super_data[(b_{abs_peak}-1)，j，i]

反射峰同时满足如下两个条件：

条件一：super_data[b_{ref_peak}，j，i]＞super_data[(b_{ref_peak}+1)，j，i]

条件二：super_data[b_{ref_peak}，j，i]＞super_data[(b_{ref_peak}-1)，j，i]

吸收峰和反射峰的数量不止一个，因此使用吸收峰动态列表list_{abs_peak}存储所有吸收峰所在波段的波长数值，使用反射峰动态列表list_{ref_peak}存储所有反射峰所在波段的波长数值；

遍历list_{abs_peak}和list_{ref_peak}，发育绢云母蚀变的像元必须同时满足如下条件：

条件一：有1个吸收峰波长位于(2210-15，2210+15)nm的开区间内

条件二：有1个吸收峰波长位于(2350-10，2350+10)nm的开区间内

条件三：有1个反射峰波长位于(2280-10，2210+10)nm的开区间内

在此基础上，将吸收峰的波长位置小于2205nm的像元厘定为短波云母亚类；

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

所述步骤四还包括：以波长为X轴反射率为Y轴建立平面直角坐标系，用于计算吸收深度的三点坐标分别为(b₂₁₃₅，B₂₁₃₅)、(b_R，B_R)、(b_C，B_C)，利用两点(b₂₁₃₅，B₂₁₃₅)、(b_R，B_R)建立直线方程，求出与垂直y轴直线(x＝b_C)的交点，继而计算出吸收深度线段长度；如方程(1)、(2)、(3)所示：

x＝b_C……………………(2)

Depth_c＝y-B_C……………………(3)

方程(1)(2)(3)联立后，求解出吸收深度

以list_{abs_peak}中位于(2210-15，2210+15)nm开区间内的吸收峰(b_C)nm为中心，两侧各取10nm的波长间隔，得到21个光谱点，在平面直角坐标系中的坐标分别如下：

(b_C-10，B_C-10)，…，(b_C，B_C)，…，(b_C+10，B_C+10)

采用最小二乘法利用抛物线(f(x)＝a0*x²+a1*x+a2)方程拟合这些光谱点，获得抛物线方程参数a0、a1、a2，进而得到拟合后的圆滑曲线方程，随后计算吸收峰b_C处的曲率Curv_C，为地质分析提供热液性质的信息；计算过程如下：

根据高等数学曲率定义，有

将f(x)＝a0*x²+a1*x+a2带入公式(4)，得到

利用公式(5)计算得出Curv_C数值；

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