实验四-彩色图像处理.doc

文档介绍：实验四彩***像处理
实验目的
掌握常见的几种颜色模型及各种颜色空间之间的相互转换。
掌握彩***像处理的方法,如彩***像增强、彩***像复原、彩***像检测与分割等。
了解伪彩色处理方法的原理及应用。
实验内容
读入一幅RGB图像,将其分别转换到CMY空间、HSI空间并显示。
利用全彩***像处理的方法,对彩***像进行平滑滤波及锐化处理。
载入一幅清晰彩***像,产生运动模糊并添加一定的高斯噪声,采用维纳滤波复原该图像。
尝试将一幅单***像采用密度分层法和灰度级-彩色变换法转化成伪彩***像。
实验步骤
1 转换到CMY空间
rgb_I=imread('');
cmy_I=plement(rgb_I);
imshow(rgb_I);
figure;
imshow(cmy_I);
HSI空间
rgb=imread('');
imshow(rgb);
rgb=im2double(rgb);
r=rgb(:,:,1);
g=rgb(:,:,2);
b=rgb(:,:,3);
I=(r+g+b)/3;
tmp1=min(min(r,g),b);
tmp2=r+g+b;
tmp2(tmp2==0)=eps;
s=1-3.*tmp1./tmp2;
tmp1=*((r-g)+(r-b));
tmp1=sqrt((r-g).^2+(r-b).*(g-b));
theta=acos(tmp1./(tmp2+eps));
h=theta;
h(b>g)=2*pi-h(b>g);
h=h/(2*pi);
h(s==0)=0;
hsi=cat(3,h,s,I);
figure;
imshow(h);
figure;
imshow(s);
figure;
imshow(I);
2、
rgb=imread('');
figure;
imshow(rgb);
r=rgb(:,:,1); %平滑滤波
g=rgb(:,:,2);
b=rgb(:,:,3);
figure;
imshow(r);
figure;
imshow(g);
figure;
imshow(b);
m=fspecial('average');
r_filtered=imfilter(r,m);
g_filtered=imfilter(g,m);
b_filtered=imfilter(b,m);
rgb_filtered=cat(3,r_filtered,g_filtered,b_filtered);
figure;
imshow(rgb_filtered);
lapMatrix=[1 1 1;1 -8 1;1 1 1]; %锐化
i_tmp=imfilter(rgb,lapMatrix,'replicate');
i_sharped=imsubtract(rgb,i_tmp);
figure;
imshow(i_sharped);
3、
i=imread('');
figure;
imshow(i);
m=fspecial('motion',20,45);
i2=imfilter(i,m,'circular');
noise=imnoise(zeros(size(i)),'gaussian',0,);
figure;
imshow(noise);
i3=double(i2)+noise;
i3=uint8(i3);
figure;
imshow(i3);
i4=deconvwnr(i3,m);
figure;
imshow(i4);
4、
a=imread('');
figure;
imshow(a);
a=double(a);
[m,n]=size(a);
l=256;
for i=1:m
for j=1:n
if a(i,j)<l/4
r(i,j)=0;
g(i,j)=4*a(i,j);
b(i,j)=l;
else if a(i,j)<=l/2
r(i,j)=0;
g(i,j)=l;
b(i,j)=-4*a(i,j)+2*l;
else if a(i,j)<=3*l/4
r(i,j)=4*a(i,j)-2*l;
g(i,j)=l;
b(i,j)=0;
else
r(i,j)=l;
g(i,j)=-4*a(i,j)+4*l;
b(i,j)=0;
end
end
end
end
end
for i=1: