http://blog.csdn.net/liuuze5/article/details/50773160

1.灰度差

1. function out_val=difference_absolute(img);
2. I=rgb2gray(img);
3. [m,n]=size(I);
4. f=0.0;
5. I=double(I);
6. for x=1:m-1;
7. for y=1:n-1;
8. Ix=I(x+1,y)-I(x,y);
9. Iy=I(x,y+1)-I(x,y);
10. f=f+abs(Ix)+abs(Iy);
11. end
12. end
13. out_val=f/m/n;

2.邊緣強度

1. function outval = edge_intensity(img)
2. % OUTVAL = EDGE_INTENSITY(IMG)
4. if nargin == 1
5. img=rgb2gray(img);
6. img = double(img);
7. % Create horizontal sobel matrix
8. w = fspecial(‘sobel‘);
10. % Get the size of img
11. [r c ] = size(img);
13. gx = imfilter(img,w,‘replicate‘);
14. gy = imfilter(img,w‘,‘replicate‘);
16. for m = 1 : r
17. for n = 1 : c
19. g(m,n) = sqrt(gx(m,n)*gx(m,n) + gy(m,n)*gy(m,n));
21. end
22. end
23. outval = mean(mean(g));
24. else
25. error(‘Wrong number of input!‘);
26. end

3.main

1. clear all;
2. close all;
3. clc;
4. PathName = ‘.\2\‘;
5. FileName = ‘.jpg‘ ;
7. n=20
8. tic
10. for ii=1:n
11. imageName=strcat(PathName,num2str(ii),FileName);
12. im = imread(imageName);
13. a(ii)=difference_absolute(im);
14. b(ii)=edge_intensity(im);
15. end
16. toc;
17. a=a/max(a);
18. b=b/max(b);
19. figure(1);plot(a,‘-r+‘);hold on
20. plot(b,‘-b*‘);hold on
21. xlabel(‘len.s position‘);
22. ylabel(‘value‘);
23. title(‘清晰度評價函數‘);

4.輸入圖片

5.對兩組結果歸一化後如下所示：

第14張圖片最清晰，與人眼觀察結果一致。

對另一組圖片進行測試，結果如下：

處理效果上邊緣算子明顯優於灰度差，但時間上消耗很大。

邊緣算子：96.872847 seconds.
灰度差： 4.645993 seconds.

圖像清晰度評價