数字图像实验六彩色图像处理Word格式文档下载.docx

数字图像实验六彩色图像处理Word格式文档下载.docx

《数字图像实验六彩色图像处理Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《数字图像实验六彩色图像处理Word格式文档下载.docx（10页珍藏版）》请在冰豆网上搜索。

）;

cmy=imcomplement（rgb）;

subplot（231）;

imshow（rgb）;

title（'

原图'

subplot（232）;

imshow（cmy）;

CMY中的图'

cmy=im2double（cmy）;

r=cmy（:

:

1）;

g=cmy（:

2）;

b=cmy（:

3）;

i=（r+g+b）/3;

m1=min（min（r,g）,b）;

m2=r+g+b;

m2（m2==0）=eps;

s=1-3.*m1./m2;

m1=0.5.*（（r-g）+（r-b））;

m1=sqrt（（r-g）.^2+（r-b）.*（g-b））;

n=acos（m1./（m1+eps））;

h=n;

h（b>

g）=2*pi-h（b>

g）;

h=h/（2*pi）;

h（s==0）=0;

hsi=cat（3,h,s,i）;

subplot（233）;

imshow（hsi）;

HSI中的图'

subplot（234）;

imshow（h）;

H分量'

subplot（235）;

imshow（s）;

S分量'

subplot（236）;

imshow（i）;

I分量'

2读入一幅图片，先对其进行平滑，再进行锐化处理，并显示。

注意请分别在GRB空间和HSI空间进行，然后在对两个空间中的结果进行比较，给出结论。

r=rgb（:

g=rgb（:

b=rgb（:

%subplot（231）;

%imshow（r）;

%subplot（232）;

%imshow（g）;

%subplot（233）;

%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）;

subplot（221）;

imshow（rgb_filtered）;

RGB空间平滑'

lapMatrix=[111;

1-81;

111];

i_tmp=imfilter（rgb_filtered,lapMatrix,'

replicate'

i_sharped=imsubtract（rgb,i_tmp）;

subplot（222）;

imshow（i_sharped）;

RGB空间平滑+锐化'

rgb=im2double（rgb）;

r_filtered=imfilter（h,m）;

g_filtered=imfilter（s,m）;

b_filtered=imfilter（i,m）;

subplot（223）;

HSI空间平滑'

i_tmp=imfilter（rgb,lapMatrix,'

i_sharped=imsubtract（rgb_filtered,i_tmp）;

subplot（224）;

HSI空间平滑+锐化'

3载入一幅清晰彩色图像，产生运动模糊并添加一定的高斯噪声，采用维纳滤波复原该图像。

4读入工具箱中others\rose.jpg图片，进行边缘检测，并将花朵分割出来

3、

i=imread（'

figure;

motion'

20,45）;

i2=imfilter（i,m,'

circular'

noise=imnoise（zeros（size（i））,'

gaussian'

0,0.05）;

imshow（noise）;

i3=double（i2）+noise;

i3=uint8（i3）;

imshow（i3）;

i4=deconvwnr（i3,m）;

imshow（i4）;

4、

a=imread（'

liftingbody.png'

imshow（a）;

a=double（a）;

[m,n]=size（a）;

l=256;

fori=1:

m

forj=1:

n

ifa（i,j）<

l/4

r（i,j）=0;

g（i,j）=4*a（i,j）;

b（i,j）=l;

elseifa（i,j）<

=l/2

g（i,j）=l;

b（i,j）=-4*a（i,j）+2*l;

=3*l/4

r（i,j）=4*a（i,j）-2*l;

b（i,j）=0;

else

r（i,j）=l;

g（i,j）=-4*a（i,j）+4*l;

end

end

g2c（i,j,1）=r（i,j）;

g2c（i,j,2）=g（i,j）;

g2c（i,j,3）=b（i,j）;

g2c=g2c/256;

imshow（g2c）;