1.一种适合于热红外高光谱遥感的SiO2含量定量计算方法,其特征在于:包括以下步骤:
步骤1.热红外高光谱遥感数据预处理;
步骤2.热红外高光谱遥感数据大气校正;
步骤3.热红外高光谱数据温度/发射率分离;
步骤4.SiO2含量定量计算经验公式确定;
步骤5.SiO2含量图获取。
2.根据权利要求1所述的一种适合于热红外高光谱遥感的SiO2含量定量计算方法,其特征在于:
步骤1中,包括辐射校正、几何校正、图像镶嵌流程,用于将数个条带拼接成整个研究区遥感图像;
步骤2包括以下步骤:
步骤2.1
对热红外高光谱遥感数据进行大气透过率曲线、上行辐射和下行辐射模拟,得到研究区大气透过率模拟曲线及上下行辐射模拟曲线;
步骤2.2
按照步骤1所采用的热红外高光谱遥感数据波段进行重采样,得到重采样后大气透过率曲线,进而得到热红外高光谱遥感数据大气校正图像;
步骤3中,利用归一化法进行温度/发射率分离,得到发射率图像;
步骤4包括以下步骤:
步骤4.1
通过JHU光谱库中各种岩性的热红外光谱曲线的观察分析,找出与SiO2 含量相关的波段范围;
步骤4.2
采集不同岩性的岩石标本进行热红外光谱数据采集;
步骤4.3
将步骤4.2采集到的热红外光谱数据进行重采样,使其与步骤1采用的热红外高光谱遥感数据相匹配,并结合步骤4.1找出SiO2含量密切相关的波段组合;
步骤4.4
对步骤4.2所述岩石标本进行地球化学测试,获取相应的SiO2含量数据;
步骤4.5
将步骤4.3得到的与SiO2含量密切相关的波段组合在不同组合方式下,进行四则运算、对数运算及指数运算,然后结合步骤4.4得到的SiO2含量进行回归分析,获取相关性最高的波段组合形式,即
|ln[B1×B2÷(B3×B4]|,
其中,B1=B8.60μm,B2=B9.81μm,B3=B8.45-8.61μm或B8.63-8.75μm,B4=B9.25-9.81μm;
步骤4.6
使用最小二乘法对公式系数a、b进行统计分析,得出SiO2含量定量计算经验公式,即SiO2%=a|ln[B1×B2÷(B3×B4]|+b;
步骤5中,将步骤4.6所述经验公式代入发射率图像,进行波段运算,获取SiO2含量图。
3.根据权利要求2所述的一种适合于热红外高光谱遥感的SiO2含量定量计算方法,其特征在于:步骤1中,采用航空热红外高光谱遥感TASI数据。
4.根据权利要求2所述的一种适合于热红外高光谱遥感的SiO2含量定 量计算方法,其特征在于:步骤4.1中,与SiO2含量相关的波段范围为8.60μm的发射峰和9.25-9.81μm波段范围,其中,8.60μm的发射峰的相对高度与SiO2含量呈正相关关系,9.25-9.81μm波段范围内的曲线上升趋势的大小与SiO2含量呈正相关关系。