[发明专利]一种基于高光谱数据识别赤铁矿化的方法

权利要求书

1.一种基于高光谱数据识别赤铁矿化的方法，其特征在于，该方法包括如下步骤：

步骤(1)获取高光谱影像数据和预处理；

步骤(2)高光谱影像数据特征波段选择；

所述的步骤(2)中采用光谱重采样方法对高光谱数据进行特征波段选择，截取可见光-近红外光谱区间400-1000nm之间的波段作为特征波段；

步骤(3)高光谱数据特征波段图像端元提取；

步骤(4)建立光谱特征识别规则区分识别赤铁矿化和褐铁矿化端元；

步骤(5)利用混合协调匹配滤波对赤铁矿化和褐铁矿化端元进行填图。

2.根据权利要求1所述的一种，其特征在于基于高光谱数据识别赤铁矿化的方法：所述的步骤(1)中对高光谱影像数据进行预处理包括辐射校正、几何校正、影像裁剪、图像掩膜，并对高光谱影像数据进行大气校正和光谱重建，得到高光谱数据的反射率光谱数据。

3.根据权利要求1所述的一种，其特征在于基于高光谱数据识别赤铁矿化的方法：所述的步骤(3)中对上述步骤(2)中得到的高光谱数据特征波段进行最小噪声分离变换、多维空间旋转、聚类，初步圈定高光谱数据特征波段图像中的端元。

4.根据权利要求3所述的一种，其特征在于基于高光谱数据识别赤铁矿化的方法：所述的步骤(4)中根据赤铁矿化和褐铁矿化的高光谱数据特征波段区间上的光谱特征差异，建立两者的区分规则，依据规则区分识别赤铁矿化和褐铁矿化端元；

获取端元光谱曲线750nm附近处的反射率最大值R1，1000nm附近处的反射率最大值R2；获取750-1000nm间反射率最小值R3，设R3的波长位置为λ3，规则如下：

①若R3＜R1＜R2，且λ3∈(800,900)，判断为赤铁矿化，λ波长单位为纳米；

②若R3＜R2＜R1，且λ3∈(900,1000)，判断为褐铁矿化，λ波长单位为纳米。

5.根据权利要求4所述的一种，其特征在于基于高光谱数据识别赤铁矿化的方法：所述的步骤(5)中对述步骤(4)中对识别出的赤铁矿化和褐铁矿化端元进行混合协调匹配滤波处理，分别得到赤铁矿化和褐铁矿化的匹配结果灰度图和不可行性灰度图，实现图像识别填图。