[发明专利]基于航空高光谱遥感数据的绢云母矿物相对丰度计算方法

摘要

本发明属于遥感地质调查技术领域，具体涉及一种基于航空高光谱遥感数据的绢云母矿物相对丰度计算方法。本发明的方法包括以下步骤对航空高光谱遥感数据进行数据预处理；计算绢云母矿物影像光谱吸收位置；去除误提取区域；计算绢云母矿物影像光谱吸收位置对应的吸收深度；制作绢云母矿物相对丰度图。本发明能够准确计算绢云母矿物吸收位置、吸收深度和相对丰度。

主权项

一种基于航空高光谱遥感数据的绢云母矿物相对丰度计算方法，其特征在于：包括以下步骤：步骤1对航空高光谱遥感数据进行数据预处理；步骤2计算绢云母矿物影像光谱吸收位置；步骤3去除误提取区域；步骤4计算绢云母矿物影像光谱吸收位置对应的吸收深度；步骤5制作绢云母矿物相对丰度图；步骤1中，对获取的航空高光谱遥感数据进行系统辐射校正和几何校正，得到具有地理坐标的光谱辐亮度数据；光谱辐亮度数据经过大气校正和光谱重建，得到高光谱反射率数据；步骤2中，应用高光谱反射率数据，依据绢云母矿物光谱特征选择高光谱数据特征波段，通过计算反射率最低值两侧的拟合直线截交位置，获得绢云母矿物准确的影像光谱吸收位置及其对应的反射率；所述反射率最低值两侧的拟合直线截交位置的计算方法如下所述：依据绢云母矿物光谱特征，初步选定影像光谱吸收位置可能处于的波段bN；选择波长小于波段bN、且与波段bN最相邻的波段bN‑1和次相邻的波段bN‑2，选择波长大于波段bN、且与bN波段最相邻的波段bN+1和次相邻的波段bN+2，采用下式计算影像光谱吸收位置bx及其对应的反射率Bx：(Bn-2-Bn-1)/(Bn-2-Bx)=(bn-1-bn-2)/(bx-bn-2)(Bn+2-Bn+1)/(Bn+2-Bx)=(bn+2-bn+1)/(bn+2-bx)]]>式中，bn‑2、bn‑1、bn+1、bn+2依次为波段bN‑2、波段bN‑1、波段bN+1、波段bN+2对应的波长；Bn‑2、Bn‑1、Bn+1、Bn+2依次为波段bN‑2、波段bN‑1、波段bN+1、波段bN+2对应的反射率。