3PSNR.rar_PSNR_PSNR彩色图像_彩色图像psnr_彩色图像评价

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%% 作者：gisplayer %%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%% 邮件地址：[email protected] %%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%% 时间：2007.10 %%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%% 欢迎与大家交流！！ %%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %本程序功能是计算2幅彩色图像的PSNR,以此来评价图像处理的效果%%%%%%%%% % 运行的结果保存在变量PSNR中 clear all close all tu1=imread('原图.bmp'); tu2=imread('处理后图像.bmp'); tu1=im2double(tu1); tu2=im2double(tu2); r1=tu1(:,:,1); %提取图像的R分量 g1=tu1(:,:,2); %提取图像的g分量 b1=tu1(:,:,3); %提取图像的b分量 r2=tu2(:,:,1); %提取图像的R分量 g2=tu2(:,:,2); %提取图像的g分量 b2=tu2(:,:,3); %提取图像的b分量 d1=mean(mean((r1-r2).^2)); m1=max(max(r1(:)),max(r2(:))); m1=m1^2; psnr1=10*log10(m1/d1); d2=mean(mean((g1-g2).^2)); m2=max(max(g1(:)),max(g2(:))); m2=m2^2; psnr2=10*log10(m2/d2); d3=mean(mean((b1-b2).^2)); m3=max(max(b1(:)),max(b2(:))); m3=m3^2; psnr3=10*log10(m3/d3); psnr=(psnr1+psnr2+psnr3)/3