【正文】 } } //顯示圖像 UpdateAllViews(NULL)。 //刪除臨時的圖像數(shù)據(jù)空間 delete ImgData。}運(yùn)行結(jié)果如圖221：圖221 圖像的低通濾波 指紋圖像的高通濾波及編程實(shí)現(xiàn)和運(yùn)行結(jié)果 圖像的高通濾波圖像的高通濾波類似于圖像的低通濾波，通過對低頻帶小波系數(shù)的衰減處理，將圖像的邊緣信息突出。如果完全去除了圖像LL頻帶的小波系數(shù)，那么將進(jìn)行圖像的邊緣提取操作，如圖222所示。 圖222 去除LL小波系數(shù)的高通濾波處理利用小波變換處理分頻帶的圖像信息處理是十分容易的，通過改變各頻帶小波系數(shù)的幅值，或者是小波系數(shù)的相位信息均能完成很好的處理。這種邊緣的處理只需要通過去除LL頻帶的小波系數(shù)，而保留高頻帶的小波系數(shù)，然后將處理后的結(jié)果進(jìn)行逆變換便能完成。其思路與低通濾波非常相似。CWFilter::HPass_Filter()實(shí)現(xiàn)了這種邊緣提取的濾波操作，CWvltDoc將圖像數(shù)據(jù)傳遞給它，函數(shù)將處理后的逆變換結(jié)果返回。完成整個濾波處理的CWvltDoc函數(shù)是CWvltDoc::OnFilterSharpness()，其對應(yīng)的菜單是“濾波處理→高通濾波1”高通濾波的實(shí)現(xiàn)及程序、結(jié)果 如果將LL的小波系數(shù)進(jìn)行有效的衰減，那么便起到圖像銳化的效果。至于衰減的系數(shù)可以隨實(shí)際應(yīng)用進(jìn)行調(diào)整，如果為了保證圖像信息的不失真，也可以通過增強(qiáng)高頻帶小波系數(shù)的方法來實(shí)現(xiàn)。 下面的函數(shù)代碼實(shí)現(xiàn)了圖像銳化處理。（1）函數(shù)描述： HPass_Filter完成圖像的高通濾波處理，得到圖像的邊緣信息。（2）函數(shù)參數(shù)： short **pData：二維指針，存放正則化后的小波系數(shù) int nHeight：圖像屬性參數(shù)，數(shù)值為原始圖像的高度值 int nHeight_H：圖像屬性參數(shù)，數(shù)值為原始圖像高度值的一半 int nWidth：圖像屬性參數(shù)，數(shù)值為原始圖像的寬度值 int nLayer：小波變換的層數(shù) （3）函數(shù)返回值：函數(shù)無返回值，小波系數(shù)參數(shù)是利用指針進(jìn)行調(diào)用（4）函數(shù)代碼void CWFilter::HPass_Filter2(short **pData, int nHeight, int nWidth,int nLayer){ short **spOriginData, **spTransData0, **spTransData1。 int iHeight = nHeight, iWidth = nWidth。 spOriginData = pData。 spTransData0 = new short *[nHeight]。 spTransData1 = new short *[nHeight]。 m_FilterData = new short * [nHeight]。 for(int i = 0。 i nHeight。 i ++) { spTransData0[i] = new short [nWidth]。 spTransData1[i] = new short [nWidth]。 m_FilterData[i] = new short [nWidth]。 } CWvltTrans *Trans。TransDWT_TriLayers(spOriginData,spTransData0,spTransData1,nHeight,nHeight/2, nWidth, nWidth/2, nLayer, )。 for(i = 0。 i nHeight。 i ++) { for(int j = 0。 j nWidth 。 j ++) { m_FilterData[i][j] = spTransData1[i][j]。 if(i nHeight / 8 amp。amp。 j nWidth / 8) m_FilterData[i][j] /= 2。 }} delete spTransData0。 delete spTransData1。 DWT_Inverse(iHeight, iWidth, 3)。} 該函數(shù)衰減了LL頻帶的小波系數(shù)：衰減為原來的1/2，而高頻頻帶的小波系數(shù)不衰減，這樣逆變換后，使得圖像的邊緣信息十分明顯。CWFilter::HPass_Filter2（）將逆變換后的數(shù)據(jù)傳遞給CWvltDoc,再通過CWvltDoc::OnFilterSharpness()來完成全部的操作，它所對應(yīng)的菜單如圖223所示，消息ID的定義如圖224所示。圖223 菜單選項(xiàng)圖224 消息ID的定義該函數(shù)代碼如下：void CWvltDoc::OnFilterSharpness2() { // TODO: Add your mand handler code here //讀取數(shù)字圖像的文件頭，獲取圖像的屬性參數(shù)LPBITMAPINFOHEADER lpBitmapInfoHeader = (LPBITMAPINFOHEADER)(m_pBitmap+14)。LPBITMAPFILEHEADER lpBitmapFileHeader = (LPBITMAPFILEHEADER)m_pBitmap。 unsigned char *lpData = m_pBitmap + lpBitmapFileHeaderbfOffBits。 unsigned long biHeight = lpBitmapInfoHeaderbiHeight。 unsigned long biWidth = lpBitmapInfoHeaderbiWidth。 unsigned long biAlign = (biWidth*3+3)/4 *4。 unsigned long bmSize = biHeight * biAlign。 if(m_pTransfered==NULL)m_pTransfered=(unsigned char*) malloc (bmSize)。 if(m_pTransfered==NULL)return。 //圖像小波系數(shù)的高通濾波層數(shù)為3層 int MaxLayer = 3。 short **ImgData。 float fTempBufforDisp。 //圖像矩陣坐標(biāo)與圖像數(shù)據(jù) int tempR, tempG, tempB, x, y, cur。 //分配圖像數(shù)據(jù)的內(nèi)存空間 ImgData = new short * [biHeight]。 for(int i = 0。 i biHeight。 i ++) { ImgData[i] = new short [biWidth]。} //獲取顯示緩存中的原始圖像數(shù)據(jù) for(y=0。 y(int)biHeight。 y++) { for(x=0。 x(int)biWidth。 x++) { cur = y*biAlign+3*x。 tempB=lpData[cur]。 tempG=lpData[cur+1]。 tempR=lpData[cur+2]。ImgData[biHeight1y][x]=(short)(*tempR+*tempG+*tempB)。 //no problem }} //創(chuàng)建小波濾波器類 CWFilter Filter。 //小波高通濾波處理 (ImgData, biHeight, biWidth, MaxLayer)。 //屏蔽圖像復(fù)原操作標(biāo)志 m_bFilter = TRUE。 //將處理后的圖像數(shù)據(jù)放入顯示緩存中 for(y=0。 y(int)biHeight。 y++) { for(x=0。 x(int)biWidth。 x++) { cur= y*biAlign+3*x。 //當(dāng)前像素的位置 m_pTransfered[cur] = (unsigned char)[biHeight 1 y][x]。m_pTransfered[cur+1]= (unsigned char)[biHeight 1 y][x]。m_pTransfered[cur+2]= (unsigned char)[biHeight 1 y][x]。 }} //顯示圖像 UpdateAllViews(NULL)。 //刪除臨時的圖像數(shù)據(jù)空間 delete ImgData。}運(yùn)行結(jié)果如圖225： 圖225 指紋圖像的高通濾波 指紋灰度圖像濾波去噪的目的是去除灰度指紋圖像中的叉連、斷點(diǎn)及模糊不清的部分．傳統(tǒng)的灰度圖濾波方法，如均值濾波、中值濾波、NN 最頻值濾波等，僅僅簡單地把指紋圖當(dāng)作灰度圖來處理，忽視了指紋圖的一個重要特性：紋線的方向性，所以效果均不理想．采用方向圖濾波可克服傳統(tǒng)灰度圖濾波算法的這一缺陷，充分利用指紋圖中紋線的方向性，可以取得很好的濾波效果．方向圖濾波由點(diǎn)方向圖求取、點(diǎn)方向圖去噪聲、塊方向圖求取、塊方向圖去噪以及方向圖濾波5個基本步驟組成．點(diǎn)方向圖表示指紋圖像中每一像素點(diǎn)灰度變化最小的方向．分別求各個方向上的灰度變化，即 (218)其中：是方向d上的第k個點(diǎn)的灰度值；N 為方向數(shù)；n為鄰點(diǎn)數(shù)．如圖211所示，取8個方向，6個鄰點(diǎn)．由于數(shù)字圖像由離散的點(diǎn)組成，計(jì)算Sd 時方向上的取樣點(diǎn)未必與圖像上的點(diǎn)一一對應(yīng)，因此，直接定義了各個方向上的取樣點(diǎn)(如圖212所示)，無需要按方向計(jì)算，從而大大提高了點(diǎn)方向圖求取的速度． 圖210 點(diǎn)方向圖 圖211 取樣點(diǎn)示意圖而相當(dāng)于在求點(diǎn)方向圖之前先進(jìn)行了一次均值濾波操作．dl是與d垂直的方向，即dl=(d+4)mod8，點(diǎn)(i,j)的方向應(yīng)該是S[d]取值最小而S[dl]取值最大的方向．這樣就不僅考慮了指紋紋線的切線方向灰度變化最小，同時考慮它的法線方向應(yīng)是灰度變化最大的方向,求得min(Sd／Sdl)，則點(diǎn)(i,j)的方向?yàn)閐．為提高點(diǎn)方向圖的精度，對其進(jìn)行一次去噪操作：對于每一點(diǎn)，若其方向與它周圍8鄰域的方向都不相同，則將此點(diǎn)視為噪音點(diǎn)，將其方向改為周圍8鄰域的主要方向．考慮到指紋紋線方向的連續(xù)性，在一個小區(qū)域內(nèi)，紋線的方向一般是一致的，所以用區(qū)域的方向去代替每一點(diǎn)的方向更有意義．將點(diǎn)方向圖分成1515大小的塊，對每一塊計(jì)算方向直方圖N(d)，求得直方圖峰值對應(yīng)的方向d，若 (219)則d即為該塊的方向，并以此方向作為該子塊內(nèi)每一點(diǎn)的方向，否則該子塊無塊方向．這樣不僅考慮了最大方向，還考慮了最大方向的兩個相鄰方向，從而在一定程度上彌補(bǔ)了用數(shù)字化方向表示連續(xù)變化方向的缺陷．塊方向圖濾波去噪思想與點(diǎn)方向圖濾波去噪思想一致．采用方向中值濾波方法對指紋圖像進(jìn)行方向?yàn)V波．經(jīng)比較，77模板的濾波效果最好，模板如圖212所示．圖212 方向中值濾波模板 方向加權(quán)中值濾波算法 在指紋圖象處理中，由于指紋圖象二值化后，在指紋的脊線上仍然存在著或多或少，大大小小的孔洞，且脊線邊緣并不光滑．即有不少的缺口和突出物．若直接對其進(jìn)行細(xì)化，則這些缺口、突出物就會使細(xì)化線變得參差不平，并出現(xiàn)一些短小的分叉，即毛刺；同時，處于紋線上的孔洞，在保持連通性的細(xì)化處理過程中將被擴(kuò)大，形成一些圓圈和網(wǎng)狀線．這不僅使細(xì)化圖變得雜亂，而且在特征提取時，將得到許多假特征點(diǎn)，從而影響以后的識別，因此，細(xì)化之前必須采取措施，對指紋二值化圖象進(jìn)行修補(bǔ)．方向中值濾波算法簡介眾所周知，中值濾波算法 由于速度快．便于實(shí)時處理，因而在圖象平滑處理中得到了廣泛應(yīng)用，但直接用中值濾波算法對方向性很強(qiáng)的指紋圖象進(jìn)行濾波．效果并不是很理想，而且指紋紋線容易出現(xiàn)斷線和粘連，因此，處理這類圖象時，有必要引人方向信息，即利用指紋方向圖來指導(dǎo)中值濾波的進(jìn)行，這種方法即所謂的方向中值濾波算法，根據(jù)文獻(xiàn)報(bào)道，方向中值濾波算法主要用于紋理性較強(qiáng)的圖象平滑．它通常采用狹長的窗口，濾波時，窗口方向隨指紋紋線方向的變化而旋轉(zhuǎn)．通常所采用的4種方向?yàn)V波窗口如圖214所示．圖213 方向中值濾波中各方向上的濾波窗口方向中值濾波法存在如下缺點(diǎn)，即若設(shè)濾波窗口長為N．則濾波的結(jié)果就會將那些寬度小于|W／2|的點(diǎn)作為噪聲濾除掉，同時還可將紋線上那些寬度小于N/2的斷裂連接起來；但它具有如下優(yōu)點(diǎn)，由于采用的是有方向的濾波．因此在濾波過程中對紋線損害較少，效果也比較好．方向加權(quán)中值濾波算法前已述，用方向中值濾波算法處理指紋圖象，效果固然好，但這種方法使用的是狹長的窗口，且濾波時