[发明专利]一种基于交替最小化的多通道盲识别自适应光学图像复原方法

权利要求书

1.一种基于交替最小化的多通道盲识别自适应光学图像复原方法，复原方法包括以下步骤：

步骤(1)：采集至少两幅大气湍流影响下的光学成像；

步骤(2)：基于多通道盲识别原理，通过拉普拉斯方法产生相对稀疏的矩阵R_Δ；

步骤(3)：使用总变分正则化(TV-Total Variation)和增广拉格朗日方法(AL M)估计目标，即U_step；

步骤(4)：在估计的目标的基础上，使用总变分正则化和增广拉格朗日方法估计点扩散函数，即H_step；

步骤(5)：对输出目标进行停止标准判别，若满足则返回求解目标，否则重复步骤(3)和步骤(4)。

2.根据权利要求1所述的一种基于交替最小化的多通道盲识别自适应光学图像复原方法，其特征在于：所述步骤(2)中构建稀疏矩阵R_Δ计算方式为：

R_Δ＝[ΔG₂,-ΔG₁]^T[ΔG₂,-ΔG₁]

其中ΔG₁和ΔG₂为第一幅成像G₁和第二幅成像G₂分别与离散拉普拉斯算子卷积后的矩阵，T表示转置，R_Δ的构建仅取决于成像G。

3.根据权利要求1所述的一种基于交替最小化的多通道盲识别自适应光学图像复原方法，所述步骤(3)中估计目标计算(U_step)步骤为：

步骤(3-1)：使用TV正则化对U_step中的目标u进行估计，即：

其中Q(u)表示求解目标U时采用的正则项，i为像素索引，是对目标u梯度求导，和分别表示目标u的x和y方向的梯度；

步骤(3-2)：矩阵形式表示式为：

其中Dx和Dy是分别是对x和y求导的矩阵，Φ是的大写形式；

步骤(3-3)：对于步骤(3-2)表示式的最小化矩阵形式为：

其中||.||²是L2范数的平方，H是PSF的矩阵形式，g为模糊成像；

步骤(3-4)：应用变量拆分，用v_x替换D_xu，用v_y替换D_yu，得到以下公式：

其中，γ是最小化求解目标u公式中保真项的权重，s.t是subject to的缩写，含义为受约束于；

步骤(3-5)采用增广拉格朗日方法进行处理，同时引入变量a_x和a_y促进解的收敛，其计算公式为：

其中，最终求解目标U的公式使用表示，α为中的惩罚系数。

4.根据权利要求1所述的一种基于交替最小化的多通道盲识别自适应光学图像复原方法，其特征在于：所述步骤(4)中估计点扩散函数(H_step)计算步骤为：

步骤(4-1)：对h_k施加稀疏性，进一步促进算法收敛，为了保证非负性和稀疏性，使用以下公式正则化：

其中δ是控制R_Δ的权重，ψ的定义通过绝对惩罚负值来确保非负性，是估计的PSF大小，h是点扩散函数，Ψ(h)表示正则公式R(h)中确保非负性的项；

步骤(4-2)：求解点扩散函数的最小化问题，求解公式为：

其中，U为目标的矩阵形式；

步骤(4-3)：应用变量拆分w＝h得到以下公式：

步骤(4-4)：采用增广拉格朗日方法，即：

其中，最终求解PSF的公式使用表示，β为中的惩罚系数。