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

权利要求书

1.一种基于航空高光谱遥感数据的绢云母矿物相对丰度计算方法，其特征在于：包括以下步骤：

步骤1对航空高光谱遥感数据进行数据预处理；

步骤2计算绢云母矿物影像光谱吸收位置；

步骤3去除误提取区域；

步骤4计算绢云母矿物影像光谱吸收位置对应的吸收深度；

步骤5制作绢云母矿物相对丰度图；

步骤1中，对获取的航空高光谱遥感数据进行系统辐射校正和几何校正，得到具有地理坐标的光谱辐亮度数据；光谱辐亮度数据经过大气校正和光谱重建，得到高光谱反射率数据；

步骤2中，应用高光谱反射率数据，依据绢云母矿物光谱特征选择高光谱数据特征波段，通过计算反射率最低值两侧的拟合直线截交位置，获得绢云母矿物准确的影像光谱吸收位置及其对应的反射率；

所述反射率最低值两侧的拟合直线截交位置的计算方法如下所述：

依据绢云母矿物光谱特征，初步选定影像光谱吸收位置可能处于的波段b_N；选择波长小于波段b_N、且与波段b_N最相邻的波段b_N-1和次相邻的波段b_N-2，选择波长大于波段b_N、且与b_N波段最相邻的波段b_N+1和次相邻的波段b_N+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)]]>

式中，

b_n-2、b_n-1、b_n+1、b_n+2依次为波段b_N-2、波段b_N-1、波段b_N+1、波段b_N+2对应的波长；

B_n-2、B_n-1、B_n+1、B_n+2依次为波段b_N-2、波段b_N-1、波段b_N+1、波段b_N+2对应的反射率。