數(shù)字圖像實(shí)驗(yàn)四

1、數(shù)字圖像處理實(shí)驗(yàn) 實(shí)驗(yàn)四班級：08電子A 姓名 陳欣 學(xué)號 0815211010CHAPTER11:5. ANSWER:A圖用程序?qū)崿F(xiàn)計(jì)算如下：四分量：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個(gè)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個(gè)8分量B圖用程序?qū)崿F(xiàn)計(jì)算如下：四分量：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個(gè)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個(gè)8分量12.ANSWER:四分量實(shí)驗(yàn)程序如下：n=imread(nodules1.tif);nt=n100;nl=bwlabel(nt,4); max(nl(:)ans =134即有134個(gè)四分量標(biāo)記八分量標(biāo)記實(shí)驗(yàn)程序如下：n=imread(nodu

5、les1.tif);nt=n100; nl=bwlabel(nt,8); max(nl(:)ans = 132即有132個(gè)八分量標(biāo)記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) 則經(jīng)驗(yàn)證均正確！CHAPTER124. ANSWER:1-1：A=zer

8、os(6,6);A(2,2:4)=1;A(3:5,2:5)=1A %圖像矩陣 c=chaincode4(A) %鏈碼數(shù)nc=normalize(c) %標(biāo)準(zhǔn)化鏈碼數(shù)lc=length(c);cl=c(2:lc) c(1);mc=mod(cl-c,4) %形狀數(shù)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) %鏈碼數(shù)nc=normalize(c) %標(biāo)準(zhǔn)化鏈碼數(shù)%A2=rot90(A);%c2=chaincode4(A2);%normalize(c2);lc=length(c);cl=c(2:lc) c(1);mc=mod(cl-c,4) %形狀數(shù)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) %鏈碼數(shù)nc=normalize(c) %標(biāo)準(zhǔn)化鏈碼數(shù)lc=length(c);cl=c(2:lc) c(1);mc=mod(cl-c,8) %形狀數(shù)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) %鏈碼數(shù)nc=normalize(c) %標(biāo)準(zhǔn)化鏈碼數(shù)lc=length(c);cl=c(2:lc) c(1);mc=mod(cl-c,4) %形狀數(shù)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結(jié)果如右圖所示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的增強(qiáng)效果圖）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的增強(qiáng)效果圖）增強(qiáng)比較好！ 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的三分量進(jìn)行邊界檢測，可以靠近如13.18所示邊界圖。對RGB三分量檢測邊界：f=imread(flowers.tif);fh=rgb2hsv(f);f1=edge(fh(:,:,1); f2=edge(fh(:,:,2); f3=edge(f