SSIM与PSNR的计算方式

有关SSIM的计算：

function re=SSIM(X,Y)　　%返回值在0-1之间，越大质量越好。X=double(X);Y=double(Y);ux=mean(mean(X));uy=mean(mean(Y));sigma2x=mean(mean((X-ux).^2));sigma2y=mean(mean((Y-uy).^2));sigmaxy=mean(mean((X-ux).*(Y-uy)));k1=0.01;k2=0.03;L=255;c1=(k1*L)^2;c2=(k2*L)^2;c3=c2/2;l=(2*ux*uy+c1)/(ux*ux+uy*uy+c1);c=(2*sqrt(sigma2x)*sqrt(sigma2y)+c2)/(sigma2x+sigma2y+c2);s=(sigmaxy+c3)/(sqrt(sigma2x)*sqrt(sigma2y)+c3);re=l*c*s;end

有关PSNR的计算：

K(i,j)为整个矩阵的均值

function re=PSNR(X,Y)　　

	bit=8;maxnum=2^bit-1;          %图像有多少灰度级

	[h,w]=size(X);temp=sum(sum((X-Y).^2))/(h*w);re=20*log10(maxnum/sqrt(temp));