数字图像实验四

1、数字图像处理实验 实验四班级：08电子A 姓名 陈欣 学号 0815211010CHAPTER11:5. ANSWER:A图用程序实现计算如下：四分量：A=magic(10)50;al=bwlabel(A);bwlabel(A,4)ans =1 1 0 0 0 4 4 4 4 01 1 0 0 0 4 4 4 4 00 1 1 0 0 4 4 4 4 01 1 0 0 0 4 4 4 4 01 1 0 0 0 4 4 4 4 00 0 3 3 3 0 0 0 0 50 0 3 3 3 0 0 0 0 52 0 0 3 3 0 0 0 0 50 0 3 3 3 0 0 0 0 50 0 3 3

2、3 0 0 0 0 5即做多有5个4分量八分量：A=magic(10)50;al=bwlabel(A);bwlabel(A,8)ans =1 1 0 0 0 1 1 1 1 01 1 0 0 0 1 1 1 1 00 1 1 0 0 1 1 1 1 01 1 0 0 0 1 1 1 1 01 1 0 0 0 1 1 1 1 00 0 1 1 1 0 0 0 0 10 0 1 1 1 0 0 0 0 12 0 0 1 1 0 0 0 0 10 0 1 1 1 0 0 0 0 10 0 1 1 1 0 0 0 0 1即做多有2个8分量B图用程序实现计算如下：四分量：B=magic(10)60;bl

3、=bwlabel(B);bwlabel(B,4)ans =1 1 0 0 0 4 4 0 0 01 1 0 0 0 4 0 0 5 00 1 1 0 0 0 0 5 5 01 1 0 0 0 0 5 5 5 01 1 0 0 0 5 5 5 0 00 0 3 3 3 0 0 0 0 60 0 3 3 3 0 0 0 0 62 0 0 3 3 0 0 0 0 60 0 3 3 3 0 0 0 0 00 0 3 3 3 0 0 0 0 0即做多有6个4分量八分量：B=magic(10)60;bl=bwlabel(B);bwlabel(B,8)ans =1 1 0 0 0 3 3 0 0 01 1

4、0 0 0 3 0 0 1 00 1 1 0 0 0 0 1 1 01 1 0 0 0 0 1 1 1 01 1 0 0 0 1 1 1 0 00 0 1 1 1 0 0 0 0 4 0 0 1 1 1 0 0 0 0 4 2 0 0 1 1 0 0 0 0 40 0 1 1 1 0 0 0 0 00 0 1 1 1 0 0 0 0 0即做多有4个8分量12.ANSWER:四分量实验程序如下：n=imread(nodules1.tif);nt=n100;nl=bwlabel(nt,4); max(nl(:)ans =134即有134个四分量标记八分量标记实验程序如下：n=imread(nodu

5、les1.tif);nt=n100; nl=bwlabel(nt,8); max(nl(:)ans = 132即有132个八分量标记14. ANSWER:4分量程序：C=zeros(11,11);C(2,5:7)=1;C(3:4,3:9)=1;C(5:7,2:10)=1;C(8:9,3:9)=1;C(10,5:7)=1;f4=(x(5)&(x(2)*x(4)*x(6)*x(8);lut=makelut(f4,3);Cw=applylut(C,lut);subplot(121),imshow(C);subplot(122),imshow(Cw)所的图像对比如右所示：8边界程序：clear all

6、C=zeros(11,11);C(2,5:7)=1;C(3:4,3:9)=1;C(5:7,2:10)=1;C(8:9,3:9)=1;C(10,5:7)=1;f8=(x(5)&(x(1)*x(2)*x(3)*x(4)*x(6)*x(7)*x(8)*x(9);lut=makelut(f8,3);Cw=applylut(C,lut);subplot(121),imshow(C);subplot(122),imshow(Cw) 所的图像对比如右所示：15. ANSWER:四边界：x,y=meshgrid(-5:5,-5:5);C=(x.2+y.2)20f4=(x(5)&(x(2)*x(4)*x(6)*

7、x(8);lut=makelut(f4,3);Cw=applylut(C,lut);subplot(121),imshow(C);subplot(122),imshow(Cw)八边界：x,y=meshgrid(-5:5,-5:5);C=(x.2+y.2)20;f8=(x(5)&(x(1)*x(2)*x(3)*x(4)*x(6)*x(7)*x(8)*x(9);lut=makelut(f8,3);Cw=applylut(C,lut);subplot(121),imshow(C);subplot(122),imshow(Cw) 则经验证均正确！CHAPTER124. ANSWER:1-1：A=zer

8、os(6,6);A(2,2:4)=1;A(3:5,2:5)=1A %图像矩阵 c=chaincode4(A) %链码数nc=normalize(c) %标准化链码数lc=length(c);cl=c(2:lc) c(1);mc=mod(cl-c,4) %形状数A =0 0 0 0 0 00 1 1 1 0 00 1 1 1 1 00 1 1 1 1 00 1 1 1 1 00 0 0 0 0 0c = 3 3 3 0 0 0 1 1 2 1 2 2nc = 0 0 0 1 1 2 1 2 2 3 3 3mc = 0 0 1 0 0 1 0 1 3 1 0 11-2：B=zeros(7,7);B

9、(2,2:4)=1;B(3,2:5)=1;B(4:5,2:6)=1;B(6,3:6)=1;B %图像矩阵 c=chaincode4(B) %链码数nc=normalize(c) %标准化链码数%A2=rot90(A);%c2=chaincode4(A2);%normalize(c2);lc=length(c);cl=c(2:lc) c(1);mc=mod(cl-c,4) %形状数B = 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 1 1 1 1 0 0 0 1 1 1 1 1 0 0 1 1 1 1 1 0 0 0 1 1 1 1 0 0 0 0 0 0 0 0c =3 3 3

10、 0 3 0 0 0 1 1 2 1 2 1 2 2nc =0 0 0 1 1 2 1 2 1 2 2 3 3 3 0 3mc =0 0 1 3 1 0 0 1 0 1 3 1 3 1 0 13-1：A=zeros(9,8);A(2,4:5)=1;A(3,3:6)=1;A(4,2:5)=1;A(5,3:5)=1;A(6,3:6)=1;A(7,3:7)=1;A(8,5:6)=1;A %图像矩阵 c=chaincode8(A) %链码数nc=normalize(c) %标准化链码数lc=length(c);cl=c(2:lc) c(1);mc=mod(cl-c,8) %形状数A =0 0 0 0

11、0 0 0 00 0 0 1 1 0 0 00 0 1 1 1 1 0 00 1 1 1 1 0 0 00 0 1 1 1 0 0 00 0 1 1 1 1 0 00 0 1 1 1 1 1 00 0 0 0 1 1 0 00 0 0 0 0 0 0 0c = 5 0 4 5 0 4 7 2 1 0 3nc =0 3 5 0 4 5 0 4 7 2 1mc = 3 4 1 3 4 3 3 7 7 3 23-2：B=zeros(8,9);B(2,4)=1;B(2,6)=1;B(3,3:7)=1;B(4,2:7)=1;B(5:6,2:8)=1;B(7,3:4)=1;B(7,6:7)=1;B %图像

12、矩阵 c=chaincode4(B) %链码数nc=normalize(c) %标准化链码数lc=length(c);cl=c(2:lc) c(1);mc=mod(cl-c,4) %形状数B =0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 1 1 1 1 1 0 0 0 1 1 1 1 1 1 0 0 0 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 0 0 0 1 1 0 1 1 0 0 0 0 0 0 0 0 0 0 0c =3 2 3 2 3 3 0 3 0 1 0 0 3 0 1 0 1 2 1 1 2 1 3 2 2 1nc =0

13、0 3 0 1 0 1 2 1 1 2 1 3 2 2 1 3 2 3 2 3 3 0 3 0 1mc =3 1 3 1 0 1 3 1 1 3 0 3 11 3 1 1 3 0 1 3 2 3 0 3 211.ANSWER:6X4：A=zeros(6,8);A(2:5,2:7)=1;A %图像形状b=boundary4(A) %边界c=complex(b(:,1),b(:,2);subplot(221); plot(c,o),axis(1,6,1,8),axis square;title(原始边界图);f=fft(c);fl=zeros(size(f);fl(1:2)=f(1:2);subp

14、lot(222); plot(ifft(fl),o),axis(1,6,1,8),axis square;title(傅里叶算子边界2);fl(1:8)=f(1:8);subplot(223); plot(ifft(fl),o),axis(1,6,1,8),axis square;title(傅里叶算子边界8);fl(1:16)=f(1:16);subplot(224); plot(ifft(fl),o),axis(1,6,1,8),axis square;title(傅里叶算子边界16);A =0 0 0 0 0 0 0 00 1 1 1 1 1 1 00 1 1 1 1 1 1 00 1

15、1 1 1 1 1 00 1 1 1 1 1 1 00 0 0 0 0 0 0 0b =3 24 25 25 35 45 55 65 74 73 72 72 62 52 42 3结果如右图所示12X8：clear all B=zeros(10,14);B(2:9,2:13)=1;B %图像形状b=boundary4(B) %边界c=complex(b(:,1),b(:,2);subplot(221); plot(c,o),axis(1,10,1,14),axis square;title(原始边界图);f=fft(c);fl=zeros(size(f);fl(1:8)=f(1:8);subpl

16、ot(222); plot(ifft(fl),o),axis(1,10,1,14),axis square;title(傅里叶算子边界8);fl(1:16)=f(1:16);subplot(223); plot(ifft(fl),o),axis(1,10,1,14),axis square;title(傅里叶算子边界16);fl(1:36)=f(1:36);subplot(224); plot(ifft(fl),o),axis(1,10,1,14),axis square;title(傅里叶算子边界36);B =0 0 0 0 0 0 0 0 0 0 0 0 0 00 1 1 1 1 1 1

17、1 1 1 1 1 1 00 1 1 1 1 1 1 1 1 1 1 1 1 00 1 1 1 1 1 1 1 1 1 1 1 1 00 1 1 1 1 1 1 1 1 1 1 1 1 00 1 1 1 1 1 1 1 1 1 1 1 1 00 1 1 1 1 1 1 1 1 1 1 1 1 00 1 1 1 1 1 1 1 1 1 1 1 1 00 1 1 1 1 1 1 1 1 1 1 1 1 00 0 0 0 0 0 0 0 0 0 0 0 0 0b = 3 24 2 5 2 6 2 7 28 29 29 39 49 59 69 79 89 99 109 119 129 138 137

18、136 135 134 133 132 132 122 112 102 92 82 72 62 52 42 32 22 2 CHAPTER138.ANSWER:HSV:a=imread(autumn.tif);figure,subplot(211),imshow(a); title(原始图像);ah=rgb2hsv(a);subplot(212),imshow(ah);title(hsv);ah(:,:,1)=histeq(ah(:,:,1);a2=ah(:,:,1);figure,subplot(221),imshow(a2);title(H);ah(:,:,2)=histeq(ah(:,:

19、,2);a3=ah(:,:,2);subplot(222),imshow(a3); title(S);ah(:,:,3)=histeq(ah(:,:,3); a4=ah(:,:,3);subplot(223),imshow(a4);title(V); YIQ: clear alla=imread(autumn.tif);figure,subplot(211),imshow(a); （ HSV的增强效果图）title(原始图像);an=rgb2ntsc(a);an(:,:,1)=histeq(an(:,:,1); a2=ntsc2rgb(an);subplot(212),imshow(a2);t

20、itle(yiq); ay=histeq(an(:,:,1);figure,subplot(311),imshow(ay);title(Y);ai=histeq(an(:,:,2);subplot(312),imshow(ai);title(I);aq=histeq(an(:,:,3);subplot(313),imshow(aq);title(Q); 由上两组图比较可以看出用HSV的 （YIQ的增强效果图）增强比较好！ 10.ANSWER:对RGB三分量检测边界： RGB边界检测f=imread(flowers.tif);f1=edge(f(:,:,1); f2=edge(f(:,:,2);f3=edge(f(:,:,3); fe2=f1|f2|f3;figure,imshow(fe2)右图所示，为对RGB的三分量进行边界检测，可以靠近如13.18所示边界图。对RGB三分量检测边界：f=imread(flowers.tif);fh=rgb2hsv(f);f1=edge(fh(:,:,1); f2=edge(fh(:,:,2); f3=edge(f