%?f=dicomread（‘d:\CT1.2.392.200036.9116.2.5.1.16.1613470806.1347435069.475977.dcm‘）;
%?fmin=min（min（f））;
%?if?fmin<0
%?????f=f-fmin;
%?end
%?fmax=max（max（f））;
%?x=double（f）/double（fmax）;

x=imread（‘d:\aa.jpg‘）;
x=rgb2gray（x）;
[MN]=size（x）;
x1=double（x）;
%基本參數(shù)
Fd=0.8;
FD=-1*Fd;
Fe=128;
Xmax=255;
%模糊特征平面
for?i=1:M
????for?j=1:N
????????P（ij）=（1+（Xmax-x1（ij））/Fe）^FD;
????end
end
%模糊增強(qiáng)
times=1;
for?k=1:times
????for?i=1:M
????????for?j=1:N
????????????if?P（ij）<=0.500
????????????????P1（ij）=2*P（ij）^2;
????????????else
????????????????P1（ij）=1-2*（1-P（ij））^2;
????????????end
????????end
????end
????P=P1;
end
%反模糊化
for?i=1:M
????for?j=1:N
????????I（ij）=Xmax-Fe*（（1/P（ij））^（1/Fd）-1）;
????end
end
X=uint8（I）;
figuresubplot（121）imshow（x）;title（‘原始圖像‘）
subplot（122）imshow（X）;title（‘模糊化后的圖像‘）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件?????????874??2012-09-01?11:04??模糊集增強(qiáng)圖像.m