基于Visual C++的數(shù)字圖像處理

1、 數(shù)字圖像處理設計1、 設計題目 基于Visual C+的圖像平滑、銳化處理2、 設計要求基于Visual C+的MFC程序設計，實現(xiàn)對灰度圖像簡單的平滑處理與拉普拉斯銳化處理。并能得到處理后所需的圖像。3、 設計原理（一）、圖像平滑 1 圖像平滑的基本概述圖像平滑是指用于突出圖像的寬大區(qū)域、低頻成分、主干部分或抑制圖像噪聲和干擾高頻成分，使圖像亮度平緩漸變，減小突變梯度，改善圖像質量的圖像處理方法。 2 圖像平滑的基本方法 （1）線性濾波（領域平均法）對一些圖像進行線性濾波可以去除圖像中某些類型的噪聲。領域平均法就是一種非常適合去除通過掃描得到的圖像中的噪聲顆粒的線性濾波。領域平均法是空間域

2、平滑噪聲技術。對于給定的圖像中的每個像素點，取其領域S。設S含有M個像素，取其平均值作為處理后所得圖像像素點處的灰度。用一像素領域內(nèi)各像素灰度平均值來代替該像素原來的灰度，即領域平均技術。領域S的形狀和大小根據(jù)圖像特點確定。一般取的形狀是正方形、矩形及十字形等，S的形狀和大小可以在全圖處理過程中保持不變，也可以根據(jù)圖像的局部統(tǒng)計特性而變化，點(m,n)一般位于S的中心。如S為33領域，點(m,n)位于S中心，則 假設噪聲n是加性噪聲，在空間各點互不相關，且期望為0，方差為，圖像g是未受污染的圖像，含有噪聲圖像f經(jīng)過加權平均后為 由上式可知，經(jīng)過平均后，噪聲的均值不變，方差，即方差變小，說明噪聲

3、強度減弱了，抑制了噪聲。 （2）中值濾波 中值濾波是一種非線性濾波，由于它在實際運算過程中并不需要圖像的統(tǒng)計特性，所以比較方便。但是對一些細節(jié)多，特別是點、線、尖頂細節(jié)多的圖像不宜采用中值濾波的方法。中值濾波的基本原理是把數(shù)字圖像或數(shù)字序列中一點的值用該點的一個鄰域中各點值的中值代替。設有一個一維序列，取窗口長度為m(m為奇數(shù))，對此序列進行中值濾波，就是從輸入序列中相繼抽出m個數(shù)，其中為窗口的中心位置，再將這m個點按其數(shù)值大小排列，取其序號為正中間的那作為輸出。用數(shù)學公式表示為： 例如，有一個序列為1，3，4，0，7，則中值濾波為重新排序后的序列0，1，3，4，7中間的值為3。此例若用平均濾

4、波，窗口也是取5，那么平均濾波輸出為。因此平均濾波的一般輸出為: 對于二維序列進行中值濾波時，濾波窗口也是二維的，但這種二維窗口可以有各種不同的形狀，如線狀、方形、圓形、十字形、圓環(huán)形等。二維數(shù)據(jù)的中值濾波可以表示為： 其中，A為濾波窗口。（2） 、圖像銳化 1 圖像銳化的基本概述圖像在傳輸過程中，通常質量都要降低，除了噪聲因素外，圖像一般都要變得模糊。這主要是因為圖像的傳輸或轉換系統(tǒng)的傳遞函數(shù)對高頻成分的衰減作用，造成圖像的細節(jié)輪廓不清晰。圖像銳化就是補償圖像的輪廓，使圖像比較清晰。 2 圖像銳化的基本方法 （1）微分算子法設有一副圖像，它的梯度用數(shù)學概念描述為一個向量，定義為 該向量的模為

5、 上式所表示的稱為圖像的梯度，它是x,y的函數(shù)，實際上也是圖像的梯度圖像。對于數(shù)字圖像來講，微分和只能用查分來近似。上式按差分運算近似后的梯度表達式為 但是，實際計算時，用這種方式計算量較大，為了節(jié)省計算機計算的時間，在計算精度允許的情況下，可蔡勇絕對差法。對應的絕對差算法為 必須注意的是，用梯度計算時，對于一個MN的數(shù)字圖像f(x,y)(x=0,M-1,y=0,N-1)，點(M-1,y)和(X,N-1)的梯度是無法進行計算的。因此，圖像f(x,y)的最后一行（第M-1行）和最后一列（第N-1列）的像素的梯度值應該進行特殊處理。常用的方法是(M-1,y)的梯度值用(M-2,y)的梯度值來代替，

6、(x,N-1)的梯度值用(x,N-2)的梯度值來代替。 （2）拉氏算子拉普拉斯算子是最簡單的各向同性微分算子，具有旋轉不變性，比較適用于改善因為光線的漫反射造成的圖像模糊。其原理是，在攝像記錄圖像的過程中，光點將光漫反射到其周圍區(qū)域，這個過程滿足擴散方程: 經(jīng)過推導，可以發(fā)現(xiàn)當圖像的模糊是由光的漫反射造成時，不模糊圖像等于模糊圖像減去它的拉普拉斯變換的常數(shù)倍。另外，人們還發(fā)現(xiàn)，即使模糊不是由于光的漫反射造成的，對圖像進行拉普拉斯變換也可以使圖像更清晰。拉普拉斯銳化的一維處理表達式是： 在二維情況下，拉普拉斯算子使走向不同的輪廓能夠在垂直的方向上具有類似于一維那樣的銳化效應，其表達式為: 對于離

7、散函數(shù)，拉氏算子定義為 其中 = = = = = 同理 類似的有 所以有 則拉式算子的模板表示為： 4、 源程序清單 #include stdafx.h #include DIPDemo.h #include DIPDemoDoc.h #include DIPDemoView.h #include mainfrm.h #ifdef _DEBUG #define new DEBUG_NEW #undef THIS_FILE static char THIS_FILE = _FILE_;#endif IMPLEMENT_DYNCREATE(CDIPDemoView, CScrollView) BE

8、GIN_MESSAGE_MAP(CDIPDemoView, CScrollView) /AFX_MSG_MAP(CDIPDemoView) ON_COMMAND(ID_ENHA_SMOOTH, OnEnhaSmooth) ON_COMMAND(ID_ENHA_SHARP, OnEnhaSharp) /AFX_MSG_MAP END_MESSAGE_MAP() CDIPDemoView:CDIPDemoView() : m_bLBtnDblClkSrv(0) / TODO: add construction code here CDIPDemoView:CDIPDemoView() void C

9、DIPDemoView:OnDraw(CDC* pDC) BeginWaitCursor(); / 顯示等待光標 CDIPDemoDoc* pDoc = GetDocument(); / 獲取文檔 ASSERT_VALID(pDoc); if(pDoc-m_Image.IsValidate() CPalette* pOldPalette; CPalette* pPalette = pDoc-GetDocPalette(); if(pPalette!=NULL) pOldPalette = pDC-SelectPalette(pPalette, FALSE); pDC-RealizePalett

10、e(); /更新系統(tǒng)調(diào)色板 pDoc-m_Image.Draw(pDC); /繪制圖像 if(pPalette!=NULL) pDC-SelectPalette(pOldPalette, FALSE); EndWaitCursor(); / 恢復正常光標 void CDIPDemoView:OnEnhaSmooth() / 圖像平滑 CDIPDemoDoc* pDoc = GetDocument(); / 獲取文檔 CImgProcess imgInput = pDoc-m_Image; / 輸入對象 if (imgInput.m_pBMIH-biBitCount!=8) AfxMessage

11、Box(不是8-bpp灰度圖像，無法處理！); return; / 檢查圖像是灰度圖 CImgProcess imgOutput = imgInput; / 輸出的臨時對象 intnTempH;/ 模板高度 intnTempW; / 模板寬度 FLOAT fTempC;/ 模板系數(shù) intnTempMX;/ 模板中心元素X坐標 intnTempMY;/ 模板中心元素Y坐標 FLOATaValue25 = 1.0, 1.0, 1.0, 0.0, 0.0, 1.0, 1.0, 1.0, 0.0, 0.0, 1.0, 1.0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0

12、.0, 0.0, 0.0, 0.0, 0.0, 0.0,; BeginWaitCursor(); / 更改光標形狀 nTempH = 3; nTempW = 3; nTempMX = 1; nTempMY = 1; fTempC = (FLOAT) (1.0 / 9.0); / 初始化對話框變量值 imgInput.Template(&imgOutput, nTempH, nTempW, nTempMY, nTempMX, aValue, fTempC); / 調(diào)用Template()函數(shù)實現(xiàn)平滑濾波 pDoc-m_Image = imgOutput; / 將結果返回給文檔類 pDoc-Set

13、ModifiedFlag(TRUE); / 設置臟標記 pDoc-UpdateAllViews(NULL); / 更新視圖 EndWaitCursor(); / 恢復光標 void CDIPDemoView:OnEnhaSharp() / 拉普拉斯銳化 CDIPDemoDoc* pDoc = GetDocument(); / 獲取文檔 CImgProcess imgInput = pDoc-m_Image; / 輸入對象 if (imgInput.m_pBMIH-biBitCount!=8) AfxMessageBox(不是8-bpp灰度圖像，無法處理！); return; / 檢查圖像是灰度

14、圖 CImgProcess imgOutput = imgInput; / 輸出的臨時對象 int nTempH; / 濾波器的高度 int nTempW; / 濾波器的寬度 int nTempMX; / 中心元素的X坐標 int nTempMY; / 中心元素的Y坐標 FLOATfCoef;/ 模板系數(shù) BeginWaitCursor(); nTempW = 3; nTempH = 3; fCoef = 1.0; nTempMX = 1; nTempMY = 1; / 設置拉普拉斯模板參數(shù) imgInput.Template(&imgOutput, nTempH, nTempW, nTemp

15、MY, nTempMX, Template_Laplacian2, fCoef); / 調(diào)用Template()函數(shù)用拉普拉斯模板銳化 pDoc-m_Image = imgOutput; / 將結果返回給文檔類 pDoc-SetModifiedFlag(TRUE); / 設置臟標記 pDoc-UpdateAllViews(NULL); / 更新視圖 EndWaitCursor(); / 恢復光標 5、 設計結果設計程序運行結果如下圖所示，運行后原始圖像如下： 進行圖像平滑處理后如下圖所示：圖像進行拉普拉斯銳化后如下示：6、 收獲和體會在這一個禮拜的圖像處理程序設計中，我經(jīng)過多次的實踐和多次的失

16、敗，終于在最后成功設計好了程序。同時深刻體會到圖像平滑讓圖像中的細節(jié)成分越來越少，而相對的，圖像銳化則是凸顯圖像的細節(jié)部分。圖像平滑處理，在很多噪聲的情況下，圖像平滑能夠有效的減弱噪聲對圖像的影響。無論是高斯噪聲、椒鹽噪聲還是乘性噪聲，都能很好的在圖像平滑中將這幾種噪聲很好的去除。雖然不能完全復原如原圖一樣的情況，但是在圖像平滑過后，噪聲在圖像中的存在確實少了許多。但是對于一些本身比較清晰的圖像，效果就不那么明顯了。圖像平滑的弊端在于會把圖像本身一些很好的細節(jié)部分給模糊掉，從而降低圖像的質量。對一些分辨率很高而且細節(jié)很唯美的圖像來說，圖像平滑來處理噪聲反而效果不是很好。而在圖像的銳化中則是讓圖像變得更為清晰。通過圖像銳化可以使圖像的