一个用matlab实现的基于图像灰度概率分布的图像模糊增强算法

suffer

1. function d = directory(v)
   
2. %求出当前路径上的所有文件列表
   
3. v='IM2';
   
4. d(1).name = '.';
   
5. counter = max(size(d));
   
6. dire = dir(v);
   
7. temp = size(dire,1);
   
8. for loop = 1:temp
   
9.     if dire(loop).isdir == 0
   
10.         counter = counter+1;
    
11.         d(counter).name = fullfile(v,dire(loop).name);  
    
12.     end
    
13.     if strcmp(dire(loop).name,'.') == 0&strcmp(dire(loop).name,'..') == 0&dire(loop).isdir == 1
    
14.         I = directory(fullfile(v,dire(loop).name));
    
15.         ma = max(size(d));
    
16.         for ii = 2:max(size(I));
    
17.             ma = ma+1;
    
18.             d(ma).name = I(ii).name;
    
19.         end
    
20.     end
    
21. end
    
22. for num=2:counter
    
23.     im=imread(d(num).name);
    
24.     im=double(im);
    
25.     mean=mean2(im);
    
26.     [m,n]=size(im);
    
27.     maxGray=max(im(:));
    
28.     minGray=min(im(:));
    
29.     b1=(minGray+mean)/2;
    
30.     b2=(maxGray+mean)/2;
    
31.     for x=1:m
    
32.         for y=1:n
    
33.             if im(x,y)>=minGray & im(x,y)<b1
    
34.                 FuzzyData(x,y)=2*(((im(x,y)-minGray)/(mean-minGray))^2);
    
35.                 im(x,y)=im(x,y)*(FuzzyData(x,y)^2);
    
36.             elseif im(x,y)>=b1 & im(x,y)<mean
    
37.                 FuzzyData(x,y)=1-(2*(((im(x,y)-mean)/(mean-minGray))^2));
    
38.                 im(x,y)=im(x,y)*(FuzzyData(x,y)^2);
    
39.             elseif im(x,y)>=mean & im(x,y)<b2
    
40.                 FuzzyData(x,y)=1-2*(((im(x,y)-mean)/(maxGray-mean))^2);
    
41.                 im(x,y)=im(x,y)*(FuzzyData(x,y)^2);
    
42.             elseif im(x,y)>=b2 & im(x,y)<maxGray
    
43.                 FuzzyData(x,y)=2*(((im(x,y)-mean)/(maxGray-mean))^2);
    
44.                 im(x,y)=im(x,y)*(FuzzyData(x,y)^2);
    
45.             end
    
46.         end
    
47.     end
    
48.     f=fullfile('F:\ExperimentResult\IM2Result1',d(num).name);
    
49.     csvwrite(f,im);
    
50. end
    
51. %csvread

复制代码