function [phase_unwrap,N]=Unwrap_TIE_DCT_Iter(phase_wrap)
phi1 = unwrap_TIE(phase_wrap);
phi1=phi1+mean2(phase_wrap)-mean2(phi1); %adjust piston
K1=round((phi1-phase_wrap)/2/pi); %calculate integer K
phase_unwrap=phase_wrap+2*K1*pi;
residue=wrapToPi(phase_unwrap-phi1);
phi1=phi1+unwrap_TIE(residue);
phi1=phi1+mean2(phase_wrap)-mean2(phi1); %adjust piston
K2=round((phi1-phase_wrap)/2/pi); %calculate integer K
phase_unwrap=phase_wrap+2*K2*pi;
residue=wrapToPi(phase_unwrap-phi1);
N=0;
while sum(sum(abs(K2-K1)))>0
K1=K2;
phic=unwrap_TIE(residue);
phi1=phi1+phic;
phi1=phi1+mean2(phase_wrap)-mean2(phi1); %adjust piston
K2=round((phi1-phase_wrap)/2/pi); %calculate integer K
phase_unwrap=phase_wrap+2*K2*pi;
residue=wrapToPi(phase_unwrap-phi1);
N=N+1;
end
end
function
基于强度方程TIE-DCT 算法-matlab
于 2024-05-19 21:18:00 首次发布