(4)
式中,矩阵ΩZ中的σi(i=1, 2, …, 2κ 1)为计算得到的奇异值;UZ和VZ分别表示Z(m1,n1)的左奇异矩阵和右奇异矩阵。
对矩阵奇异值分解后,舍去较小奇异值可达到去噪的效果,据此构造降噪之后的复干涉相位矩阵Z(m1,n1)*如下[18]
(5)
式中,[0]2κ×(2κ 1)为2κ×(2κ 1)的零矩阵;σ*为奇异值的主要部分。
从Z(m1,n1)*中创建子矩阵Y0、Y1和Y2
(6)
式中,Y0表示取Z(m1,n1)*前2κ行和列;Y1表示取Z(m1,n1)*后2κ行和前2κ列;Y2表示取Z(m1,n1)*前2κ行和后2κ列,即Y0和Y1以及Y2等3子矩阵尺寸为2κ×2κ。对Y0做奇异值分解,并构建Y0*、Y1*和Y2*
(7)
式中,U00和V00分别表示Y0的左奇异矢量和右奇异矢量,包含了奇异矩阵的主要元素;U01和V01则为余下的奇异矩阵。
计算相位梯度估计值g(m1,n1)|(m2,n2)