[发明专利]低秩联合协同表示的高光谱遥感影像异常探测方法

权利要求书

1.一种低秩联合协同表示的高光谱遥感影像异常探测方法，其特征在于：包括以下步骤：

步骤1，将进行异常探测的高光谱遥感影像转换为二维矩阵其中m为波段数，n为像元数；

步骤2，采用均值漂移算法对X进行聚类，获得C个类簇，搜索半径为r，每一类簇的聚类中心为

步骤3，对于每一个类簇，执行以下步骤：

(a)统计类簇内像元个数

(b)若m_i≥a，其中a为类簇大小阈值，继续执行步骤(c)；若m_ia,执行下一次循环的步骤(a)；

(c)计算类簇内各个像元与聚类中心c_i的欧氏距离，按照升序排列并选取前b个对应的像元参与构建字典D；

(d)重复执行步骤(a)-(c)，直到对所有类簇循环完毕，获取字典其中s为字典的原子数；

步骤4，计算D中每一个原子和X中每一个像元之间的欧氏距离，确定权重矩阵W＝(w_ij)_s×n，其中w_ij为D中的第i个原子与X中第j个像元之间的欧氏距离；

步骤5，将构建的目标函数写成增广的拉格朗日函数形式，并采用LADMAP方法迭代求解，具体执行以下步骤：

(a)初始化所有变量α₀＝J₀＝Y_2,0＝0，E₀＝Y_1,0＝0，μ₀＝0.01，μ_max＝10¹⁰，ρ₀＝1.1，ε₁＝10^-6，ε₂＝10^-2，k＝0，其中α₀,J₀为初始化的系数矩阵，Y_1,0和Y_2,0为初始的拉格朗日乘子，E₀为初始的异常目标矩阵；

(b)更新其中Θ为奇异值阈值；

(c)更新J_k+1＝(μ_kα_k+1+Y_2,k)./(μ_k+2βW)，其中β为惩罚系数；

(d)更新其中λ为惩罚系数，Ω为l_2,1最小化算子；

(e)更新拉格朗日乘子Y_1,k+1＝Y_1,k+μ_k(X-Dα_k+1-E_k+1)，Y_2,k+1＝Y_2,k+μ_k(α_k+1-J_k+1)；

(f)计算

(g)更新μ_k+1＝min(μ_max,ρμ_k)；

(h)判断迭代是否收敛，若收敛则输出α_k+1,J_k+1,E_k+1；若不收敛则令k＝k+1,重复执行步骤(b)-(g)直至迭代收敛；

步骤6，输出异常图像。

2.根据权利要求1所述低秩联合协同表示的高光谱遥感影像异常探测方法，其特征在于：所述步骤4中建立的采用欧氏距离为权重的权重矩阵W＝(w_ij)_s×n，其具体实现公式为：

w_ij＝||D_i-X_j||₂

式中，D_i为D中的第i个字典原子，X_j为X中的第j个像元。