[发明专利]基于分裂Bregman的干涉高光谱图像快速分解算法

权利要求书

1.一种基于分裂Bregman的干涉高光谱图像快速分解算法，其特征在于：该算法包括以下步骤：

(1)读取一帧待分解的干涉高光谱图像X作为输入参数；

(2)设置X_B表示背景层，设置X_I表示干涉条纹层，二者与X的关系是：X_I＝X-X_B，设置分裂Bregman迭代初始化参数d_x，d_y，b_x，b_y，设置水平方向全变分权值参数λ₁，设置竖直方向全变分权值参数λ₂；

(3)根据分裂Bregman原理利用干涉高光谱图像的单方向性特征进行内层迭代，对X_B，d_x，d_y，b_x，b_y进行寻优操作；

(4)根据分裂Bregman原理进行外层迭代，对变量X进行Bregman更新；

(5)完成Bregman迭代求出背景层最优解X_B，通过X_I＝X-X_B求得干涉条纹层最优解；

(6)完成图像分解。

2.根据权利要求1所述的一种基于分裂Bregman的干涉高光谱图像快速分解算法，其特征在于：步骤(2)所述设置分裂Bregman迭代初始化参数d_x，d_y，b_x，b_y，均初始化为与输入参数X同等大小的矩阵，矩阵中的元素数值全部设置为0。

3.根据权利要求1所述一种基于分裂Bregman的干涉高光谱图像快速分解算法，其特征在于：步骤(3)所述的内层迭代中，用k表示外层迭代次数(内层迭代次数不采用变量表示)，初始化k＝0，令X^k＝X，背景层初始化为X_B^k＝X^k，根据分裂Bregman原理利用干涉高光谱图像的单方向性特征进行内层迭代，对X_B，d_x，d_y，b_x，b_y进行寻优操作：

其中fft2，ifft2分别表示傅里叶变换与逆傅里叶变换，表示卷积，分别表示[0 0 0；0 1 -1；0 0 0]，[0 0 0；-1 1 0；0 0 0]，[0 0 0；0 1 0；0 -1 0]，[0 -1 0；0 1 0；0 0 0]的Sobel卷积核算子。

4.根据权利要求1所述一种基于分裂Bregman的干涉高光谱图像快速分解算法，其特征在于：步骤(4)所述的外层迭代，对输入数据X进行Bregman更新：

(2)k＝k+1。

5.根据权利要求1所述一种基于分裂Bregman的干涉高光谱图像快速分解算法，其特征在于：步骤(5)所述，进行指定次数的外层迭代和内层迭代，即可快速求出背景层最优解X_B，通过X_I＝X-X_B求得干涉条纹层最优解。