资源简介
自编的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
评论
共有 条评论