matlab食物素材去噪代码(毕业设计或大作业个人学习参考代码)

clear all; %clc; indi=load('dat.mat'); indianbe=indi.I0; indianbe=reshape(indianbe,145,145,186); %indianbe=indianbe(1:50,1:50,:); [a,b,c]=size(indianbe); %add noise and calculate the noise image's SNR for i=1:c indian(:,:,i)=imnoise(indianbe(:,:,i),'gaussian',0,0.05); end S=SNR(indianbe,indian); fprintf('The Noise Image SNR is %d.\n',S); indianbet=tensor(indianbe); indiant=tensor(indian); indian1= tensor_unfolding(indian,1); %mode 1 unfolding indian2= tensor_unfolding(indian,2); %mode 2 unfolding indian3= tensor_unfolding(indian,3); %mode 3 unfolding %K1=MWFAIC(indian1);%estimate rank of mode n K1=80; %K2=MWFAIC(indian2); K2=80; %K3=MWFAIC(indian3); K3=50; k=0; thr=1e-9; H10=eye(a);H20=eye(b);H30=eye(c); Rt0=indiant; %tensor Rt1=ttm(indiant,{H10,H20,H30},[1:3]); %dif=norm(Rt1-indianbet); while k<1 %dif>thr && k<1000 Rt11=ttm(indiant,{H20,H30},[2,3]); Rt11d=double(Rt11); Rt11= tensor_unfolding(Rt11d,1); %mode 1 unfolding R1=indian1; %R10=R1*(kron(H20,H30))'; H11=MWFilter(Rt11,R1,K1,a); Rt12=ttm(indiant,{H11,H30},[1,3]); Rt12d=double(Rt12); Rt12= tensor_unfolding(Rt12d,2); %mode 2 unfolding R2=indian2; %R20=R2*(kron(H11,H30))'; H21=MWFilter(Rt12,R2,K2,b); Rt13=ttm(indiant,{H11,H21},[1,2]); Rt13d=double(Rt13); Rt13= tensor_unfolding(Rt13d,3); %mode 3 unfolding R3=indian3; %R30=R3*(kron(H11,H21))'; H31=MWFilter(Rt13,R3,K3,c); % H10=H11; %update % H20=H21; % H30=H31; % Rt0=Rt1; k=k+1; Rt1=ttm(indiant,{H11,H21,H31},[1:3]); %dif=norm(Rt1-Rt0); %fprintf('Iteration No.%d processing done. dif= %f \n',k, dif); Rt1d=double(Rt1); S=SNR(indianbe,Rt1d); fprintf('The Image denoised SNR is %d.\n',S); end fprintf('All iteration done, now get the filtered signal tensor.\n'); estX=ttm(indiant,{H11,H21,H31},[1:3]); filtX=(double(estX)); S=SNR(indianbe,filtX); fprintf('The Image denoised SNR is %d.\n',S); %70, 86 and 136 filtX_rgb(:,:,1)=filtX(:,:,70); filtX_rgb(:,:,2)=filtX(:,:,86); filtX_rgb(:,:,3)=filtX(:,:,136); figure, imshow(uint8(filtX_rgb*255)); filtX_rgb(:,:,1)=indian(:,:,70); filtX_rgb(:,:,2)=indian(:,:,86); filtX_rgb(:,:,3)=indian(:,:,136); figure, imshow(uint8(filtX_rgb*255)); filtX_rgb(:,:,1)=indianbe(:,:,70); filtX_rgb(:,:,2)=indianbe(:,:,86); filtX_rgb(:,:,3)=indianbe(:,:,136); figure, imshow(uint8(filtX_rgb*255));