实验四-彩色图像处理

实验四彩色图像处理实验目的掌握常见的几种颜色模型及各种颜色空间之间的相互转换。掌握彩色图像处理的方法，如彩色图像增强、彩色图像复原、彩色图像检测与分割等。了解伪彩色处理方法的原理及应用。实验内容读入一幅RGB图像，将其分别转换到CMY空间、HSI空间并显示。利用全彩色图像处理的方法，对彩色图像进行平滑滤涉及锐化处理。载入一幅清晰彩色图像，产生运动模糊并添加一定的高斯噪声，采用维纳滤波复原该图像。尝试将一幅单色图像采用密度分层法和灰度级-彩色变换法转化成伪彩色图像。实验步骤1转换到CMY空间rgb_I=imread('peppers.png');cmy_I=imcomplement(rgb_I);imshow(rgb_I);figure;imshow(cmy_I);HSI空间rgb=imread('peppers.png');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(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('peppers.png');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=[111;1-81;111];%锐化i_tmp=imfilter(rgb,lapMatrix,'replicate');i_sharped=imsubtract(rgb,i_tmp);figure;imshow(i_sharped);3、i=imread('peppers.png');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);4、a=imread('liftingbody.png');figure;imshow(a);a=double(a);[m,n]=size(a);l=256;fori=1:mforj=1:nifa(i,j)<l/4r(i,j)=0;g(i,j)=4*a(i,j);b(i,j)=l;elseifa(i,j)<=l/2r(i,j)=0;g(i,j)=l;b(i,j)=-4*a(i,j)+2*l;elseifa(i,j)<=3*l/4r(i,j)=4*a(i,j)-2*l;g(i,j)=l;b(i,j)=0;elser(i,j)=l;g(i,j)=-4*a(i,j)+4*l;b(i,j)=0;endendendendendfori=1:mforj=1:ng2c(i,