[发明专利]基于同质图斑的高光谱影像异常探测方法

权利要求书

1.一种基于同质图斑的高光谱异常探测方法，其特征在于，包括步骤：

S1、采用光谱降维法对原始高光谱数据X进行变换，得到变换后的特征向量Y；

S2、采用多光谱影像分割方法对降维后的影像Y进行分割，得到标记影像L；

S3、搜索标记影像L中标记号相同的像素，形成像素集，每一个像素集对应着一个同质斑块；

S4、采用加权法计算每个同质斑块的均值及方差；

S5、以同质斑块为单位，依次计算每个像素的异质度；

S6、确定异质度的阈值，依次判断每个像素是否为异常点。

2.如权利要求1所述的方法，其特征在于，在步骤S1中变换方法具体为：

Y=W_PCAX    (1)

式中:W_PCA为PCA的投影矩阵，其由光谱数据协方差矩阵的特征值按大小排列后，前m个特征值所对应的列特征向量组成的矩阵。

3.如权利要求1所述的方法，其特征在于，在步骤S2中，采用基于背景梯度分布模型的分水岭影像分割方法对多光谱影像Y进行分割，得到标记图像L。

4.如权利要求3所述的方法，其特征在于，所述步骤S2具体包括：

其中，影像的梯度影像IG具体为：

IG=maxi{IGi}---(2)]]>

式中，为每一个波段的梯度影像，其中水平方向的梯度算子及竖直方向梯度算子分别为GDx=-xexp(x2+y22σ2)/(πσ2)]]>及GDy=-y·exp(x2+y22·σ2)/(π·σ2),]]>σ为高斯函数的方差，x、y为影像的行列坐标值；

IGB=IG⊗LD---(4)]]>

式中，为高斯低通滤波算子；

标记影像定义为

IM=1，if(IG>IG_B)    (5)

其中，IM为标记影像，当梯度值大于梯度趋势影像值时为1，否则为0，IG为梯度影像，IG_B为梯度趋势影像，采用分水岭影像分割方法进行影像分割，得到最终的标记影像L。