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