psnr.rar_PSNR_信噪比计算_峰值信噪比PSNR

%Project: Calculates the PSNR (Peak Signal to Noise Ratio) % of images A and A', both of size MxN function [A] = psnr(base,result) % convert to doubles % base=double(base); % result=double(result); base=imread('D:\超分辨率重建\超分辨率重建\result_KE_PO_1.tif'); [M,N]=size(base); b_base=block_im(base,8); %b_basew为原始图像 b_result=block_im(result,8); %b_result为处理结果 for i=1:8 for j=1:8 % m_base(i,j)=max(max(double(b_base{i,j}).^2)); m_base(i,j)=max(max(im2double(b_base{i,j}).^2)); end end % psnr_num=double(M*N*max(m_base(:))); psnr_num=M*N*max(m_base(:)); for i=1:8 for j=1:8 % s_dif(i,j)=sum(sum(double(b_dif{i,j}).^2)); s_dif(i,j)=sum(sum(double(im2double(b_base{i,j})-im2double(b_result{i,j})).^2)); end end % psnr_den=double(sum(s_dif(:))); psnr_den=sum(s_dif(:)); % psnr_num=M*N*max(max(base.^2)); % calculate numerator用信号功率最大值乘以象素点个数)); % calculate den % psnr_den=sum(sum((base-result).^2ominator 噪声的功率之和 A=10*log10(psnr_num/psnr_den); % calculate PSNR