資源描述:
《談vc編程實(shí)現(xiàn)數(shù)字圖像的邊沿檢測(cè)》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在學(xué)術(shù)論文-天天文庫(kù)。
1、談VC編程實(shí)現(xiàn)數(shù)字圖像的邊沿檢測(cè)數(shù)字圖像的邊沿檢測(cè)是圖像分割、目標(biāo)區(qū)域的識(shí)別、區(qū)域外形提取等圖像分析領(lǐng)域十分重要的基礎(chǔ),圖像理解和分析的第一步往往就是邊沿檢測(cè),目前它以成為機(jī)器視覺探究領(lǐng)域最活躍的課題之一,在工程應(yīng)用中占有十分重要的地位。本文向讀者簡(jiǎn)單先容一下這個(gè)技術(shù),并給出了在VisualC環(huán)境下實(shí)現(xiàn)的代碼。 所謂邊沿就是指圖像局部亮度變化最明顯的部分,它是檢測(cè)圖像局部變化明顯變化的最基本的運(yùn)算。對(duì)于數(shù)字圖像,圖像灰度灰度值的明顯變化可以用梯度來(lái)表示,以邊沿檢測(cè)Sobel算子為例來(lái)講述數(shù)字圖像處理中邊沿檢測(cè)的實(shí)現(xiàn): 對(duì)于數(shù)字圖像,可以用一階差分代替
2、一階微分; △xf(x,y)=f(x,y)-f(x-1,y); △yf(x,y)=f(x,y)-f(x,y-1) 求梯度時(shí)對(duì)于平方和運(yùn)算及開方運(yùn)算,可以用兩個(gè)分量的盡對(duì)值之和表示,即: G[f(x,y)]={[△xf(x,y)][△yf(x,y)]}
3、△xf(x,y)
4、
5、△yf(x,y)
6、; Sobel梯度算子是先做成加權(quán)均勻,再微分,然后求梯度,即: △xf(x,y)=f(x-1,y1)2f(x,y1)f(x1,y1)-f(x-1,y-1)-2f(x,y-1)-f(x1,y-1); △yf(x,y)=f(x-1,y-1)2f(x-1,y)
7、f(x-1,y1)-f(x1,y-1)-2f(x1,y)-f(x1,y1); G[f(x,y)]=
8、△xf(x,y)
9、
10、△yf(x,y)
11、; 上述各式中的像素之間的關(guān)系見圖f(x-1,y-1)f(x,y-1)f(x1,y-1)f(x-1,y)f(x,y)f(x1,y)f(x-1,y1)f(x,y1)f(x1,y1)我在視圖類中定義了響應(yīng)菜單命令的邊沿檢測(cè)Sobel算子實(shí)現(xiàn)灰度圖像邊沿檢測(cè)的函數(shù):voidCDibVieent(); HDIBhdib; unsignedchar*hData; unsignedchar*data; hdib=pDoc->m
12、_hDIB; BeginAPINFOHEADER)GlobalLock((HGLOBAL)hdib); hData=lpbi*(LPDodifiedFlag(TRUE);//設(shè)修改標(biāo)志為"TRUE" data1handle=GlobalAlloc(GMEM_SHARE,WIDTHBYTES(lpBi->biWidth*8)*lpBi->biHeight); //申請(qǐng)存放處理后的像素值的緩沖區(qū) data=(unsignedchar*)GlobalLock((HGLOBAL)data1handle); AfxGetApp()->BeginWaitCursor
13、(); inti,j,buf,buf1,buf2; for(j=0;jbiHeight;j)//以下循環(huán)求(x,y)位置的灰度值 for(i=0;ibiWidth;i) { if(((i-1)>=0)&&((i1)biWidth)&&((j-1)>=0)&&((j1)biHeight)) {//對(duì)于圖像四面邊界處的向素點(diǎn)不處理 buf1=(int)*(hData(i1)*WIDTHBYTES(lpBi->biWidth*8)(j-1)) 2*(int)*(hData(i1)*WIDTHBYTES(lpBi->biWidth*
14、8)(j)) (int)(int)*(hData(i1)*WIDTHBYTES(lpBi->biWidth*8)(j1)); buf1=buf1-(int)(int)*(hData(i-1)*WIDTHBYTES(lpBi->biWidth*8)(j-1)) -2*(int)(int)*(hData(i-1)*WIDTHBYTES(lpBi->biWidth*8)(j)) -(int)(int)*(hData(i-1)*WIDTHBYTES(lpBi->biWidth*8)(j1)); //x方向加權(quán)微分
15、 buf2=(int)(int)*(hData(i-1)*WIDTHBYTES(lpBi->biWidth*8)(j1)) 2*(int)(int)*(hData(i)*WIDTHBYTES(lpBi->biWidth*8)(j1)) (int)(int)*(hData(i1)*WIDTHBYTES(lpBi->biWidth*8)(j1)); buf2=buf2-(int)(int)*(hData(i-1)*WIDTHBYTES(lpBi->biWidth*8)(j-1)) -2*(int)(int)*(h
16、Data(i)*WIDTHBYTES(lpBi->biWidth*8)(j-1)