峰值信噪比PSNR

Halcom

1. %% 压缩比函数
   
2. function psnr=PSNR(A,B)
   
3. sizeA=size(A);
   
4. sizeB=size(B);
   
5. if sizeA~=sizeB
   
6.     error('图像尺寸不一致')
   
7. end
   
8. if A==B
   
9.     error('两幅图像一样')
   
10. end
    
11. max2_A=max(max(A));
    
12. max2_B=max(max(B));
    
13. min2_A=min(min(A));
    
14. min2_B=min(min(B));
    
15. if max2_A>255 || max2_B>255 || min2_A<0 || min2_B<0  % 灰度值约束
    
16.     error('输入的灰度值范围应为[0,255]')
    
17. end
    
18. error_diff=A-B;
    
19. decibels=20*log10(255/(sqrt(mean(mean(error_diff.^2)))));
    
20. psnr=decibels;
    
21. % disp(sprintf('PSNR=+%5.2fdB',decibels))
    

复制代码