彩色图像处理.doc.doc

实验四 彩色图像处理 实验目的 掌握常见的几种颜色模型及各种颜色空间之间的相互转换。 掌握彩色图像处理的方法，如彩色图像增强、彩色图像复原、彩色图像检测与分割等。 了解伪彩色处理方法的原理及应用。 实验内容 读入一幅RGB图像，将其分别转换到CMY空间、HSI空间并显示。 rgb_I=imread(E:\学习\数字图像\实验结果\b.jpg); cmy_I=imcomplement(rgb_I); imshow(rgb_I); figure; imshow(cmy_I); rgb=imread(E:\学习\数字图像\实验结果\b.jpg); 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=0.5*((r-g)+(r-b)); tmp1=sqrt((r-g).^2+(r-b).*(g-b)); theta=acos(tmp1./(tmp2+eps)); h=theta; h(bg)=2*pi-h(bg); h=h/(2*pi); h(s==0)=0; hsi=cat(3,h,s,I); figure; imshow(h); figure; imshow(s); figure; imshow(I); 利用全彩色图像处理的方法，对彩色图像进行平滑滤波及锐化处理。 rgb=imread(E:\学习\数字图像\实验结果\b.jpg); 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); 载入一幅清晰彩色图像，产生运动模糊并添加一定的高斯噪声，采用维纳滤波复原该图像。 i=imread(E:\学习\数字图像\实验结果\b.jpg); figure; imshow(i); m=fspecial(motion,20,45); i2=imfilter(i,m,circular); noise=imnoise(zeros(size(i)),gaussian,0,0.05); figure; imshow(noise); i3=double(i2)+noise; i3=uint8(i3); figure; imshow(i3); i4=deconvwnr(i3,m); figure; imshow(i4); 尝试将一幅单色图像采用密度分层法和灰度级-彩色变换法转化成伪彩色图像。 a=imread(pillsetc.png); 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;