flpd.m

可以直接用

function DetectionMatrix=fLPD(HyperCube,T) 
%LPD函数
%T为主成分含量,应为百分数
HyperCube=double(HyperCube);
[m,n,l]=size(HyperCube);%使用的源数据
Imagery1=reshape(HyperCube,m*n,l);%将3维图象立方体化为2维矩阵
Imagery2D=Imagery1';

BackgroundSignature=LPDBkSignature(HyperCube,T);

P=OrthCompletProjection(BackgroundSignature);%消弱背景的投影算子
ImagerySuppressB=P*Imagery2D;%获取去除背景后的数据，寻找目标过程中的投影均使用此数据

TempImage1=sum(ImagerySuppressB);
DetectionMatrix=reshape(TempImage1',m,n);