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?	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