資源描述：
《VC編程實現(xiàn)數(shù)字圖像的邊緣檢測》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、VC編程實現(xiàn)數(shù)字圖像的邊緣檢測數(shù)字圖像的邊緣檢測是圖像分割、目標(biāo)區(qū)域的識別、區(qū)域形狀提取等圖像分析領(lǐng)域十分重要的基礎(chǔ)，圖像理解和分析的第一步往往就是邊緣檢測，目前它以成為機器視覺研究領(lǐng)域最活躍的課題之一，在工程應(yīng)用中占有十分重要的地位。本文向讀者簡單介紹一下這個技術(shù)，并給出了在VisualC++環(huán)境下實現(xiàn)的代碼?！ ∷^邊緣就是指圖像局部亮度變化最顯著的部分，它是檢測圖像局部變化顯著變化的最基本的運算。對于數(shù)字圖像，圖像灰度灰度值的顯著變化可以用梯度來表示，以邊緣檢測Sobel算子為例來講述數(shù)字圖像處理中邊緣檢測的實現(xiàn)： 對于數(shù)字圖像，可以用一階
2、差分代替一階微分； △xf(x,y)=f(x,y)-f(x-1,y); △yf(x,y)=f(x,y)-f(x,y-1) 求梯度時對于平方和運算及開方運算，可以用兩個分量的絕對值之和表示，即： 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,y+1)+2f(x,y+1)+f(x+1,y+1)-f(x-1,y-1)-2f(x,y-1)-f(x+1,y-1); △yf(x,y)=f(x
7、-1,y-1)+2f(x-1,y)+f(x-1,y+1)-f(x+1,y-1)-2f(x+1,y)-f(x+1,y+1); 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(x+1,y-1)f(x-1,y)f(x,y)f(x+1,y)f(x-1,y+1)f(x,y+1)f(x+1,y+1)?我在視圖類中定義了響應(yīng)菜單命令的邊緣檢測Sobel算子實現(xiàn)灰度圖像邊緣檢測的函數(shù)：voidCDibView::OnMENUSobel()//灰度圖像數(shù)據(jù)的獲得參見天極網(wǎng)9.
12、10日發(fā)表的拙作//VC數(shù)字圖像處理一文{ HANDLEdata1handle; LPBITMAPINFOHEADERlpBi; CDibDoc*pDoc=GetDocument(); HDIBhdib; unsignedchar*hData; unsignedchar*data; hdib=pDoc->m_hDIB; BeginWaitCursor(); lpBi=(LPBITMAPINFOHEADER)GlobalLock((HGLOBAL)hdib); hData=lpbi+*(LPDWORD)lpbi+256*sizeof(RGBQUAD)
13、; //得到指向位圖像素值的指針 pDoc->SetModifiedFlag(TRUE);//設(shè)修改標(biāo)志為"TRUE" data1handle=GlobalAlloc(GMEM_SHARE,WIDTHBYTES(lpBi->biWidth*8)*lpBi->biHeight); //申請存放處理后的像素值的緩沖區(qū) data=(unsignedchar*)GlobalLock((HGLOBAL)data1handle); AfxGetApp()->BeginWaitCursor(); inti,j,buf,buf1,buf2; for(j=0;jbi
14、Height;j++)//以下循環(huán)求(x,y)位置的灰度值 for(i=0;ibiWidth;i++) { if(((i-1)>=0)&&((i+1)biWidth)&&((j-1)>=0)&&((j+1)biHeight)) {//對于圖像四周邊界處的向素點不處理 buf1=(int)*(hData+(i+1)*WIDTHBYTES(lpBi->biWidth*8)+(j-1)) +2*(int)*(hData+(i+1)*WIDTHBYTES(lpBi->biWidth*8)+(j)) +(int)(i
15、nt)*(hData+(i+1)*WIDTHBYTES(lpBi->biWidth*8)+(j+1)); 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)+(j+1)); //x方向加權(quán)微分 buf
16、2=(int)(int)*(hData+(i-1)*WIDTHBYTES(lpBi->biWidth*8)+(j+1))
