[发明专利]一种适合于热红外高光谱遥感的石英含量定量计算方法

权利要求书

1.一种适合于热红外高光谱遥感的石英矿物含量定量计算方法，其特征在于：具体包括以下步骤：

步骤一、对研究区的热红外高光谱遥感数据预处理；

步骤二、热红外高光谱遥感数据大气校正，具体包括以下子步骤：

步骤2.1热红外高光谱遥感数据进行大气透过率曲线、上行辐射和下行辐射模拟，得到研究区大气透过率模拟曲线及上下行辐射模拟曲线；

步骤2.2按照步骤一所采用的热红外高光谱遥感数据波段进行重采样，得到重采样后大气透过率曲线，进而得到热红外高光谱遥感数据大气校正图像；

步骤三、热红外高光谱数据温度/发射率分离，利用归一化法进行温度/发射率分离，得到发射率图像；

步骤四、石英矿物含量定量计算经验公式确定，具体包括以下子步骤：

步骤4.1热红外光谱曲线的观察分析，找出与石英矿物含量相关的波段范围；

步骤4.2采集不同岩性的岩浆岩标本进行热红外光谱数据采集；

步骤4.3将步骤4.2采集到的热红外光谱数据进行重采样，使其与步骤一采用的热红外高光谱遥感数据相匹配，并结合步骤4.1找出石英矿物含量密切相关的波段组合；

步骤4.4对步骤4.2所述岩石标本进行地球化学测试，获取全量分析数据；

步骤4.5根据各样品全量分析数据，使用CIPW标准矿物计算方法计算出各样品的石英矿物含量；

步骤4.6将步骤4.3得到的与石英含量密切相关的波段组合在不同组合方式下，进行四则运算、对数运算及指数运算，然后结合步骤4.4得到的石英含量进行回归分析，获取相关性最高的波段组合形式；

步骤4.7根据石英矿物含量及其对应的发射率光谱曲线中相关性最高的波段组合，得出石英矿物含量定量计算经验公式；

步骤五、石英矿物含量图获取，将步骤4.7所述经验公式代入发射率图像，进行波段运算，获取石英矿物含量图。

2.根据权利要求1所述的一种适合于热红外高光谱遥感的石英矿物含量定量计算方法，其特征在于：所述的步骤三中，石英含量进行回归分析，相关性最高的波段组合形式为ln[B2^2/(B1*B3)]；其中，B1＝B_8.49，B2＝B_8.60μm，B3＝B_8.71μm，B_8.82可作为B_8.49的备用波段，下标表示该波段的波长。

3.根据权利要求2所述的一种适合于热红外高光谱遥感的石英矿物含量定量计算方法，其特征在于：得出石英矿物含量定量计算经验公式，y＝8208.9x³-5525x²+1228.2x+0.3139，x＝ln[B2^2/(B1*B3]。

4.根据权利要求3所述的一种适合于热红外高光谱遥感的石英矿物含量定量计算方法，其特征在于：结合了CIPW标准矿物计算方法，以获取SiO₂、TiO₂、Al₂O₃、TFe₂O₃、MnO、MgO、CaO、Na₂O、K₂O及P₂O₅的全量分析数据计算出各样品石英矿物含量。

5.根据权利要求4所述的一种适合于热红外高光谱遥感的石英矿物含量定量计算方法，其特征在于：与石英矿物含量相关的波段范围具体为8.60μm的发射峰至两侧发射谷范围的光谱特征。