基于主成分分析PCA相位提取算法程序

基于主成分分析（PCA）相位提取算法，可以用于大于等于3幅干涉图的相位提取。

function phi_cor = PCA_J（I_input）
[mnk] = size（I_input）;
X = reshape（I_inputm*nk）;
Xf = zeros（m*nk）;
for i = 1:k
	temp = X（:i）;
	%Compound image after filtering the mask
	Xf（:i）= temp;
end

[M N] = size（Xf）;
Xm = mean（Xf2）;
Xd = Xf-repmat（Xm1N）;
if （N < M） %less images than image length
	C = Xd‘*Xd;
	[V D Vt] = svd（C）;
	U = Xd*V;
	U = U./repmat（sqrt（diag（D）‘）M1）;
else %more images than image length
	C = Xd*Xd‘;
	[U D Ut] = svd（C）;
end
    U（:1） = max（U（:2））.*（U（:1）./max（U（:1）））;
    U1 = X（:1）.*0;
    U2 = X（:2）.*0;
    
    U1 = U（:1）;
    U2 = U（:2）;
    
    phi_cor = reshape（atan2（-U2U1）mn）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         663  2017-08-16 15:00  PCA_J.m
     文件        1015  2017-07-25 11:16  PCA_test_demo.m