數(shù)字圖像處理實驗報告(圖像邊緣檢測)

1、實驗報告 實驗名稱 實驗三圖像邊緣檢測 課程名稱數(shù)字圖像處理 姓名成績 班級學號 日期地點 備注： 1、實驗?zāi)康?(1) 了解并掌握使用微分算子進行圖像邊緣檢測得基本原理； (2) 編寫程序使用 Laplaci a n算子(二階導(dǎo)數(shù)算子)實現(xiàn)圖像銳化，進一步理解圖像 銳化得實質(zhì)； (3) 掌握使用不同梯度算子(一階導(dǎo)數(shù)算子)進行圖像邊緣檢測得原理、方法，根據(jù) 實驗結(jié)果分析各種算子得工作效果； (4) 總結(jié)實驗過程(實驗報告，左側(cè)裝訂)：方案、編程、調(diào)試、結(jié)果、分析、結(jié)論。 2、實驗環(huán)境 (1 ) Win do ws XP / 7 (2) Matla b 7、1/7、14 3、實驗方法 本次實驗

2、要求對2 56X2 56大小，2 56級灰度得數(shù)字圖像lena、img進行處理。 (1)對該圖像進行銳化處理，要求采用La p laci an算子進行銳化，分 1與 2兩種情況， 按如下不同情況進行處理： gi ( m ,n)f(m , n )f g2 (m,n)4f (m , n)f (m1 ,n) f (m 1,n) f (m, n 1) f (m, n 1) I、 要對圖像進行處理，要先讀取該圖像，實驗代碼如下： c los e all; c l ear al l; f i d=f o pen ( le n a、img , r); im a ge=fread(fid , 256,2 5

3、6 ， ui n t8z); fclo se( f i d); II、讀取圖像后，對該圖像得每一像素(不考慮圖像得邊界部分)進行遍歷，根據(jù)公 式(公式相當于做差分)對每一灰度進行計算，將所得得結(jié)果存入一矩陣gl中 (矩陣g1 初始化為該圖像得矩陣)，代碼如下(僅以a=1為例): g 1= ima ge; a = l; x, y= siz e (i m a g e); f o r i=2：( x -1) f o r j =2 :(y 1)?l (i, j )=(1+4 * a)* image(i,j)-a *( image (i + 1,j) +image (i-1,j) +i m age (i

4、,j+1 ) +i ma ge (i , j 1); end end III 、根據(jù)公式對圖像得每一個像素 (不考慮圖像得邊界部分)進行計算，將所得之存 入矩陣2中(g2初始化值為該圖像得矩陣值 )，具體方法與上一步類似，代碼如下(僅以a =l為例)： g2=i ma ge; a=l ; x, y =s ize (ima g e); fo r i = 2： (x-1) for j= 2： (y1) ?g 2 (i, j) = 4 *a * i ma g e( i, j ) a *( ima ge(i+1,j) +ima g e (i -1, j) +im a ge(i ,j+1 )+image

5、( i,j-1 ); e nd end (2 )分別利用 Ro b e r ts、P r ew i tt與Sobel邊緣檢測算子，對原圖像進行邊緣檢測， 顯 示處理前、后圖像。 I 、同(1)中I，不再贅述。 II 、 對圖像進行邊緣檢測,要對圖像得每一像素(不考慮圖像得邊界部分)得灰度進行 遍歷,分別用公式對圖像得水平與垂直方向得邊緣進行檢測,并分別存儲,在進行合成。采 用不同得算子時,變換得公式有所不同。在合成時可以采用三種算法(街區(qū)法、棋盤法、歐 式幾何法 ),本次實驗中均又采用,所以得到三層圖像。代碼如下(僅以街區(qū)法得 Rob er ts 算子為例)： x ,y=size ( i m

6、age); g1 =i mage; gh1 = image; gvl =im a ge; f or i=2:( x 1) fo r j =2： (y1 ) gh1(i ,j)=image( i,j )ima ge(i1, j1); gv1(i, j)=image (i, j1) image(i1, j); gl( i,j)=abs (gh1 (i,j)+ abs(g v1(i,j ); end en d 4、實驗結(jié)果分析 ( 1)、圖像銳化 圖像銳化結(jié)果如圖a與圖b （注:圖a為情況結(jié)果；圖b為情況結(jié)果） 原團像 K6 f | I H fni i-l 9=2時膽的團像 圖b 由圖a與圖b對比可

7、知，圖像得銳化實質(zhì)就就是原圖像與梯度信息進行疊加，相對于原 圖像而言，圖像得邊緣信息得到了加強。由圖b中不同a值所得圖像對比可知，a越大，邊界 越清晰，同時對比圖a中不同a值所得圖像可知，a得大小決定梯度信息與原圖像疊加時所占 得比例 （2）、圖像邊緣檢測 街區(qū)法、棋盤法、歐式距離法得不同檢測算子結(jié)果分別入圖c、圖d與圖e 圖c 圖d -J ll I Prewitt H / 圖e 由圖c、圖d與圖e對比可知，采用街區(qū)法所得得邊界清晰度最高，歐式幾何法其次， 棋 盤法最差。對比同一張圖中得不同子圖可知，S o bel算子得性能更好，Prewitt算子其次，R o Sobel W 子 原團像 Ro

8、berts算子 Roberts算子 Prewitt H F Sobe 子 原凰像 I-1 IIL -lil J 、Z-|. I J- 亠JF1. -1 嚴 A Q I i|-A .予* Af 塔亠b % if * P 丨 1 A Al| V t. L 鳥“ id j J 1 ( J J ； j * 1 / 咲前 /-J jF 1/ 片i r 一二鼻 -J Jl _ 、| | f b erts算子最差。 5、實驗結(jié)論 本次實驗就是對圖像進行銳化與邊緣檢測,通過這次實驗 ,我對數(shù)字圖像處理中圖像得銳 化與邊緣檢測有了更加深刻得了解.加深了圖像銳化與邊緣檢測得原理，掌握了圖像邊緣檢 測得不同方法。學

9、會了使用L a placian算子(二階導(dǎo)數(shù)算子)實現(xiàn)圖像銳化得程序編寫；掌 握了使用不同梯度算子 (一階導(dǎo)數(shù)算子 )進行圖像邊緣檢測得原理、方法。 6、源代碼 (1) 圖像銳化(gl) clo se a ll; c lear a l l; f i d = f o pen (T ena、im g ,rz); im age=fread (fid , 256,2 5 6 , uin t 8); fc lose(fi d); sub plot(2 2 1); ims ho w (ui n t8 (image)，門 ); ti tle (原圖像); 采用拉普拉斯算子 g 1=i ma ge； a=l；

10、 x, y =size(i m age) ; fo r i=2 : (x 1) for j=2 : ( y 1) g 1 (i, j)=( l +4*a ) *imag e( i, j) a* (image (i +1, j) +i m ag e( i-1 , j)+ image (i,j +1) + ima g e(i, j 1); end e nd subpl o t( 22 2 ); i msh ow ( ui nt 8(g1)，口 ); tit l e( a= l 時 g1 得圖像); g1=i mage; a =2 ; x,y =size( im age); f o r i = 2:

11、(x 1) for j=2 ：(y 1 ) g l (i, j )=(1+4 * a) * ima g e (i ,j) a* (image (i+1,j)+ima g e(i-1, j ) +image (i, j+1) + image (i,j 1); end en d subpl ot(2 23) ； i m sh o w ( u in t 8 (g1)，口); title (a=2 時 g1 得圖像); (2) 圖像銳化( g2) cl ose all; cl ear all; f i d=fopen ( I e na、img, rz ); i mage= f r ea d (fid

12、, 256, 2 5 6 , uint 8); f c Io s e (f id); subpIot( 221)； imsh o w ( uint8 (i m age), ); title (原圖像/); 采用拉普拉斯算子 g2 =im a ge; a=l ； x ,y =s i ze( i ma ge ); f or i=2:(x 1) fo r j=2: (y 1) g2(i, j) =4aimage(i,j)a*(image(i1, j) +image(il, j) +ima g e( i , j +1) + i m age (i , j -1); en d end su bp 1 ot

13、 (222); im s how(u i n t 8 ( g 2),); ti t leCa =1 時 g2 得圖像); g2 =im age ; a= 2; x, y= s i z e(i m ag e ); for i=2: ( x 1 ) fo r j=2 ：( y-1) g2(i,j) =4aimage(i ,j)-a*(i mag e( i+1 , j)im age(i1, j) +i mag e(i, jl )+image (i, j-1) ; end e nd subplot ( 223) ; imshow( uin t8(g2), ); tit le (/ a=2 時 g 2得

14、圖像/ ); ( 3)邊緣檢測 clo se all; c1e ar a11; fid = f o p e n (I e na、i mg, rz ); i ma ge=f rea d(fid, 256, 25 6,uint8); fclo se( fi d); fi gu re( Nam ez ,街區(qū)法,Nu mb e r T i tle, o ff ); sub plot (221); im sho w( uint8( image), ) ; titl e (原圖像 )； x,y =size (im age)； gl =image ; gh1=im a ge; gv 1 =ima g e;

15、for i =2 : (x 1) fo r j =2:(y-1 ) gh1 (i, j)=im a ge(i, j ) im ag e (i 1, j-1); gv1(i , j)= ima ge( i ,j 1)-image (i-1,j); g1(i，j )=a bs(gh1(i,j) +abs(gv1(i,j)； end en d subpl ot(22 2 ); im s h ow ( uin t8(g1)，門); ti t le ( R o b erts 算子); g2 =im a ge; gh 2=im a ge; gv2 =image; fo r i =2： (x-1) for

16、j =2:(y1) gh2(i ,j )=(ima ge(i+1,j1)+imag e(i +1,j)+image(i+1 ,j+1 )/3) (image(i1,j-1)+image(i1 ,j)+image(i1，j +1)/3)； gv2(i， j)=(image(i1，j +1)+i mag e(i,j+1)+imag e(i+1,j 1)/3)(ima ge(i -1,j1)+i mage(i,j1)+imag e(i1,j -1)/3); g 2( i, j ) =abs( gh2 (i, j )+abs (gv2( i, j); end end s ubpl ot (223)；

17、ims how( uint8(g2),門); ti tl e( Prewitt 算子 )； g3 =im ag e; g h 3 =image; gv 3 =im a ge; for i =2: (x 1) for j =2 :(y1) gh3(i，j)=(imag e(i+1 ,j-1)+2 *image(i1,j)+image(i+1 ,j+1) -i mage (i 1,j 1) -2*imag e (i 1, j) image(i 1, j+ 1) )/ 4 ; gv3(i， j) =(image(i-1，j +1)+2im age(i， j +1)+i mage(i+1， j +1)

18、im age(i-1， j-1)2ima ge(i ,j 1) -image(i1， j1) /4； g3(i,j) =abs(gh 3(i， j) +abs(gv3(i ,j)； g3(i ， j)=max(gh3(i ，j),gv3(i, j) )； %g 3(i,j) =sqrt(gh3(i， j) *gh 3(i，j) +gv3(i，j )*gv3(i,j) ； end end s ubpl ot(2 24)； ims ho w( u i nt 8(g3), ); tit le( Sobel 算子)； % subplot(23 5 ) % i mshow(ui nt 8(gh 3 ),

19、); title (檢測水平邊界 ); % subpl o t(236); % ims ho w (u int 8(g v3), 口)； % title (檢測豎直邊界 )； f i gu re( Name,棋盤法/ , Num b erT it l e , o ff ; su b p l ot( 221); i m show(uin t8( imag e),); title (原圖像/); x,y=size ( image) ; g 1=im age ; gh 1=im ag e; gv1=i mag e; fo r i=2:(x-1 ) for j=2 ：( y-1) gh1 (i, j)

20、 =i m age( i, j) - i m a ge(i - 1, j 1); gvl( i,j) =image( i,j 1) image (i-1, j); %g 1(i,j) =ab s(gh1(i， j) +ab s(gv1(i， j) ); g 1(i ,j) =ma x(gh1(i， j)， gv 1(i， j) ); end end subpl o t ( 2 2 2); im s how ( ui n t 8( g l),); ti tle( Robe rts 算子) ; g 2 = i ma ge； gh2 =i mage； gv2 =imag e; for i=2：(x1

21、) for j=2 ：( y 1 ) gh 2(i， j) =(i mag e(i1， j-1) +image (i1， j) +image( i+1， j+1)3)( (image( i-1,j-1 )imag e(i-1， j)image(i- 1,j+l ) /3); gv2( i， j)=(i mage( i-1， j+l )+imag e(i,j +1)+image(i+l， j+1) )/3)-(image(i 1， jl) +i mag e(i,j1) +image(i +l， j-1) /3); %g2 (i，j )=abs(gh2(i， j)+abs( gv2(i,j) );

22、 g2(i,j )=max(gh2( i,j)， gv2(i,j) ) ; end end subplot ( 223); imsh ow(uin t8( g2),); t itl eCP rew i t t 算子)； g3 =image; gh3 =image; gv3 =imag e; f or i=2: (x 1) for j=2 ：( y-1) gh3 (i, j) = (i m a ge (i+ 1 , j 1) + 2* im age(i + 1, j) + im ag e (i+1 , j+ 1) im age (i -1, j 1)-2 *i m a ge (i 1, j )-

23、image( i 1, j+1) /4； gv3(i， j)=(image(i 1，j 1)2* image(i， j+1) +image(i +1， j+1)image(i 1,j1) 2*image(i， j1)image (i+1,j1) /4； g3 (i, j ) =abs (g h 3(i, j )+abs (gv 3(i, j); g3( i, j)=max( g h3( i,j) ,gv3 (i ,j); %g3(i ,j )=sqr t (gh3(i, j ) * gh3 (i,j) + g v3(i,j )* g v 3 (i, j ); end end su bp lot

24、 (22 4 ); i ms how(uin t 8 (g3),); title (Sobe l 算子); fi g ure (/ Name , / 歐式幾何法,N umberTitle , of f); subplot ( 221)； i m show (uint 8(i mage), ); t i tl e( 原圖像 ); x, y =size (i m age) ; g1 =i m age; g h1=im a ge; gv 1=i mage; f o r i=2 : (x 1) for j=2 ：(y -1 ) gh 1(i， j)=image( i，j)-i mage(i 1,j1)

25、 ; gv1(i， j)=image (i ,j-1)image(i 1， j); %g 1(i ,j) =ab s(gh1 (i, j )+a b s(g v1 (i, j ); %g1(i,j)=max (g h1 (i ,j)，gv 1(i， j) ); g1(i， j)=sqr t(gh1(i， j)* gh1(i， j) +gv1(i， j) *gv1(i ,j)； end end subplot ( 222) ; i mshow( ui nt 8(g1 ), )； t itle( Robe rts 算子 ); g2=i ma ge； g h2=i m a ge； g v2=i m age； f or i= 2:( x 1 ) f or j=2: (y-1) gh2(i,j) =(image(i+1， j1) +i mage(i +1， j)ima ge( i+ 1 ,j+1 ) )/3) - (image(i 1, j -1) +image (i 1,j)+i m age (i -1, j+1) /3 ); gv2(i,j)=(im age(i-1， j+1)+image(i,j+1) imag e( i