pal-King、模糊熵增强算法

自编的pal-King增强算法代码和基于模糊熵的增强算法代码

function y=enhancement（r）
X1=imread（‘lena.jpg‘）;
%X1=imread（‘lirong.jpg‘）;
X1=rgb2gray（X1）;
subplot（221）;
imshow（X1）;
title（‘（a）原始图像‘）;%显示原始图像

%直方图均衡化方法增强的图像
tic
I=histeq（X1）;
subplot（222）;
imshow（I）;
title（‘（b）直方图均衡化增强后的图像‘）
toc
t1=toc%返回程序所用时间

[mn]=size（X1）;
X2=double（X1）;

%Pal-King算法
%图像模糊化
tic
for i=1:m
    for j=1:n             
        f=X2（ij）;           
        g（ij）=（1+（255-f）/309）^（-2）;
    end
end
%模糊增强
for k=1:r
    for i=1:m  
        for j=1:n       
            f=g（ij）;        
            if（0<= f <= 0.5）           
                T（ij）=2*f^2;       
            elseif（0.5<= f <=1）            
                T（ij）=1-2*（1-f）^2;        
            end
            g（ij）=T（ij）;
        end 
    end
end
%进行对比增强算子的逆变换
min=（1+255/309）^（-2）;
for i=1:m  
    for j=1:n
        f=g（ij）;
        if（f            f=min;
        else
            f=g（ij）; 
        end
        T（ij）=255+309*（1-f^（-1/2））;
    end
end
I=Mat2gray（T）;
subplot（223）;
imshow（I）;
title（‘（c）Pal-King算法增强的图像‘）;
toc
t2=toc

%本文算法
%图像模糊化
tic
for i=1:m
    for j=1:n             
        f=X2（ij）;           
        g（ij）=（sin（3.14*（f/255-0.5））+1）/2;       
    end
end
%模糊增强
q=fazhi（X1）;
q=（sin（3.14*（q/255-0.5））+1）/2;  %模糊熵最大化选择的阀值
for k=1:r
    for i=1:m  
        for j=1:n       
            f=g（ij）;        
            if（0<= f <= q）           
                T（ij）=f^2/q;       
            elseif（q<= f <=1）            
                T（ij）=1-（1-f）^2/（1-q）;        
            end
            g（ij）=T（ij）;
        end 
    end  
end
subplot（224）;
imshow（T）;
title（‘（d）本文算法增强的图像‘）;
toc
t3=toc

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1820  2012-06-12 21:34  enhancement.m

----------- ---------  ---------- -----  ----

                 1820                    1