基于子空间模型GLRT的高光谱目标检测方法

clear;clc; I= imread('1.tif'); data = im2double(I); [m,n,p]=size(data);%(610,340,103) t=m*n;%207400 data=reshape(data,t,p);%(207400,103),即(样本数，波段数) %% 二进制编码 10.14 % % data = dataT; % filter1 = [-1/2,1,-1/2]; % data = filter2(filter1,dataT,'same') ; %% data2=isinf(data); [inf_r inf_c]=find(data2==1); data(inf_r,:)=[]; [pc,score,latent,tsquare]=pca(data);%pc为主成分系数，score为主成分的结果，latent为方差 feature_after_PCA=score(:,1:20); RES=reshape(feature_after_PCA,m,n,20); save('RES.mat','RES'); red = RES(:,:,1); blue= RES(:,:,2) ; green = RES(:,:,3); red = mat2gray(red); blue = mat2gray(blue); green = mat2gray(green); red = im2uint8(red); blue = im2uint8(blue); green = im2uint8(green); photo=cat(3,red,blue,green); % imshow(photo); % imwrite(RES(:,:,1),'1.jpg'); % imwrite(RES(:,:,2),'2.jpg'); % imwrite(RES(:,:,3),'3.jpg'); % imwrite(RES,'PU_pca_3.jpg') %cumsum(latent)./sum(latent) %ans =0.5832,0.9442,0.9886,0.9916,0.9937