图像 matlab 峰值信噪比 im2double double

简介
PSNR 是最普遍，最广泛使用的评鉴画质的客观量测法，不过许多实验结果都显示，PSNR 的分数无法和人眼看到的视觉品质完全一致，有可能 PSNR 较高者看起来反而比PSNR 较低者差。
这是因为人眼的视觉对于误差的敏感度并不是绝对的，其感知结果会受到许多因素的影响而产生变化（例如：人眼对空间频率较低的对比差异敏感度较高，人眼对亮度对比差异的敏感度较色度高，人眼对一个区域的感知结果会受到其周围邻近区域的影响）。
（参考百度百科）

计算公式：

2010-12-15 21:45:44 上传

下载附件(10.22 KB)

2010-12-15 21:45:41 上传

下载附件(4.23 KB)

代码

1. % 图像峰值信噪比PSNR、均方根误差MSE
   
2. % By lyqmath
   
3. % Dalian University of Technology
   
4. % School of Mathematical Sciences
   
5. function [PSNR, MSE] = psnr(X, Y)
   
6. % 计算峰值信噪比PSNR、均方根误差MSE
   
7. % 如果输入Y为空，则视为X与其本身来计算PSNR、MSE
   
8. 
   
9. if nargin<2
   
10.     D = X;
    
11. else
    
12.     if any(size(X)~=size(Y))
    
13.         error('The input size is not equal to each other!');
    
14.     end
    
15.     D = X-Y;
    
16. end
    
17. MSE = sum(D(:).*D(:))/prod(size(X));
    
18. PSNR = 10*log10(255^2/MSE);

复制代码

实例

1. function main()
   
2. clc; close all;
   
3. I = imread('rice.png');
   
4. I1 = imnoise(I, 'salt & pepper');
   
5. figure;
   
6. subplot(1, 2, 1); imshow(I); title('原图像');
   
7. subplot(1, 2, 2); imshow(I1); title('加噪声图像');
   
8. [PSNR, MSE] = psnr(I, I1)
   
9. 
   
10. % 图像峰值信噪比PSNR、均方根误差MSE
    
11. % By lyqmath
    
12. % Dalian University of Technology
    
13. % School of Mathematical Sciences
    
14. function [PSNR, MSE] = psnr(X, Y)
    
15. % 计算峰值信噪比PSNR、均方根误差MSE
    
16. % 如果输入Y为空，则视为X与其本身来计算PSNR、MSE
    
17. 
    
18. if nargin<2
    
19.     D = X;
    
20. else
    
21.     if any(size(X)~=size(Y))
    
22.         error('The input size is not equal to each other!');
    
23.     end
    
24.     D = X-Y;
    
25. end
    
26. MSE = sum(D(:).*D(:))/prod(size(X));
    
27. PSNR = 10*log10(255^2/MSE);

复制代码

2010-12-15 21:53:36 上传

下载附件(75.09 KB)

2010-12-15 21:54:31 上传

下载附件(16.56 KB)

总结
有一些方法中采用了分块计算PSNR等参数，用来得到更为有效地信息。也有的方法提出一些新的计算公式，当然这都是为了统计处图像处理前后的差异信息。
类似的，对于相近的公式，也可以如上做编码。

注：部分内容参考网络共享资料。

 

 

坑爹啊有木有！！！

尼玛图像要转换格式有木有！！！

a simple test will show you that
when you use im2double, it will convert your uint8 image from the
range of 0-255 into a double array in the range 0-1. In general, all
Image Processing Toolbox functions that take floating point arrays
want them in the range of 0 to 1. I find it kind of annoying but
that's the way it is. That's basically what "rescaling" means. So
this imdouble() function is there to help you get to what you might
need for later calls to Image Processing Toolbox functions. If you
use double() you would have to divide by 255 yourself if you needed
to. If you're doing your own custom functions then you can leave it
in the original range if you want.

 

Using im2double will rescale the data (if necessary) from your image.
Using double won't.

怪不得总出错！！！

铭记！！！