图像峰值信噪比计算

Matlab代码如下

function S=PSNR(s,t)

[m,n,~]=size(s);
s=im2uint8(mat2gray(s));
t=im2uint8(mat2gray(t));
s=double(s);
t=double(t);

sd=0;
mi=m*n*max(max(s.^2));
for u=1:n
    for v=1:n
        sd=sd+(s(u,v)-t(u,v))^2;
    end
end
if sd==0
    sd=1;
end
S=mi/sd;
S=10*log10(S);