compare.m

用matlab实现的JPEG算法

%  compare(name,'sortie.bmp') ;


function compare(FID,name_in,num1,name_out,num2,taille) 

name_in
name_out

[M_in,nb_comp1]=init_image(FID,name_in);
[M_out,nb_comp2]=init_image(FID,name_out);

size(M_in);
size(M_out);
n1 = double(nb_comp1) ;
n2 = double(nb_comp2) ;

%aff_image(num1,M_in); title(name_in);
aff_image(num2,M_out); title(name_out);

%return 

if n1 ~= n2 
   disp('erreur de taille') 
   return 
end   

fprintf(FID,'***** image in = %s / image out = %s *****\n',name_in, name_out);
for i=1:nb_comp1
   
	[A(i),B(i),C(i)] = mesure_PSNR(M_in(:,:,i),M_out(:,:,i));
	[A(i),B(i),C(i)] = mesure_PSNR(M_in(:,:,i),M_out(:,:,i));
	[A(i),B(i),C(i)] = mesure_PSNR(M_in(:,:,i),M_out(:,:,i));
	fprintf(FID,'composante %d : EQM = %f, SNR = %f, PSNR = %f \n',i,A(i),B(i),C(i));

end 
[M,N] = size(M_in(:,:,i)) ;

taux = nb_comp1*M*N/taille ;
st_taux = num2str(taux) ;

st_rate = num2str(24/taux) ;
st_taux = strcat('rate=',st_rate,'-->facteur X ',st_taux);

title(st_taux);

fprintf(FID,'taille =%d -> (%g bpp) taux de %g \n',taille,taille*8/(M*N),taux);


%============================================================