map=0
[X,map]=imread('lena256x256x8.bmp')%读取图像
p=imfinfo('lena256x256x8.bmp')%读取图像的属性
s=2^p.BitDepth-1%计算图像的最大深度,如8位图像为255.
%-------------------------------------------------------------------------
%简化的图像增强算法:简化Pmn及T变换的算法为:
Y=0%清空Y.
A=0%清空A.
Y=double(X)./s%计算P(m,n),相当于求G(x)
%模糊变换算法
for m=1:p.Height
for n=1:p.Width
if(Y(m,n)<0.5)
A(m,n)=2*Y(m,n)^2
else
A(m,n)=1-2*(1-Y(m,n))^2
end
end
end
B=0%清空B.
B=uint8(A*(2^p.BitDepth-1))%相当于求此贺逆G-1(x')
%%%区域减影增强算法
P=0%清空P.
d=0.00005%%%%设定冗余度
k=2%%%%设定修正参数
for m=1:p.Height
for n=1:p.Width
if(Y(m,n)-A(m,n)>d)
P(m,n)=Y(m,n)+k*(Y(m,n)-A(m,n))
end
if(Y(m,n)-A(m,n)<(-d))
P(m,n)=Y(m,n)+(-k)*(Y(m,n)-A(m,n))
else
P(m,n)=Y(m,n)
end
end
end
B=0%清空B.
B=uint8(P*(2^p.BitDepth-1))%相当于求逆G-1(x')
%%%图像增强几种方法效果对比
I=imread('lena256x256x8.bmp') %读入图像
I=double(I) %变成双精度
K1=filter2(fspecial('average',5),I)/255 %进行滤正大波
K2=imsubtract(I,K1)
K2=immultiply(K2,0.5)
K3=imadd(I,K2)
K2=uint8(K2) %转成8位数据
K3=uint8(K3)
I=uint8(I)
figure
title('图像增强几种方举扒竖法效果对比')
subplot(2,2,1)imshow(I)title('原始图像')
subplot(2,2,2)imshow(K1)title('低通滤波图像')
subplot(2,2,3)imshow(K3)title('线性反锐化掩膜图像')
subplot(2,2,4)imshow(B)title('区域减影增强图像')
从图像可知,这两条曲线分别是sin(x),cos(x)。实现这两笑改条曲线的绘制可以这样编写其代码:
x=0:pi/10:4*pi %生成x向量组,其范围从0到4π,步长为pi/10
y1=sin(x) %计算正弦函数值
y2=cos(x) %计算余弦函数值
plot(x,y1,'--',x,y2,'-') %绘制正弦函数和余弦函数曲线图
legend("sin(x)","cos(x)") %图例标碰含判注老键
xlabel('时间')ylabel('幅值')%标注坐标轴名称
执行代码和运行结果
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)