干涉极值法提取相位,干涉法与相位法相比哪种测量方法更好,matlab

%干涉极值相位提取算法% http://www.doc88.com/p-3877508231882.html function phi_cor = EVI(I) % 函数说明：输入的I为MxNx2的矩阵，两幅干涉图的大小为MxN，输入时合并为一个参数I, 输出的phi_cor表示通过EVI提取出的二维包裹相位。 [I_mx ,I_my] = find(I(:,:,1) == max(max(I(:,:,1)))); delta_E_max = acos(I(I_mx,I_my,2)./I(I_mx,I_my,1)); [m1,n1] = size(delta_E_max); delta_E1 = sum(sum(delta_E_max)); [I_ix ,I_iy] = find(I(:,:,1) == min(min(I(:,:,1)))); delta_E_min = acos(I(I_ix,I_iy,2)./I(I_ix,I_iy,1)); [m2,n2] = size(delta_E_min); delta_E2 = sum(sum(delta_E_min)); delta_E = (delta_E1+ delta_E2)/(m1*n1+m2*n2); Fz = I(:,:,1)*cos(delta_E) - I(:,:,2); Fm = I(:,:,1)*sin(delta_E); phi_cor = atan2(Fz,Fm);