基于PM方程的指紋掃描儀圖像增強算法研究

劉曉薇

（江西科技師范大學(xué)數(shù)計學(xué)院，江西 南昌330013）

目前的指紋掃描設(shè)備采用的都是基于光學(xué)結(jié)構(gòu)和CMOS（complementary metal oxide semilonductor）圖像傳感器的采集原理。在高質(zhì)量指紋采集掃描儀中CMOS圖像傳感器約占整機總成本的30%～50%，為了減少成本，縮小體積，大多數(shù)指紋掃描儀都會選用CFA 模式的CMOS 圖像傳感器，現(xiàn)在應(yīng)用最多的是采用GRBG 排列布局的Bayer 格式陣列，因為綠色光譜對應(yīng)人眼對可見光光譜敏感度的峰值，采用雙倍的綠色像素數(shù)可以加強感光效果［1］。

即使在CFA（cdor filter array）模式傳感器中，因為感光元的大小差異，分辨率不同，價格也有10～20倍的差異，一些CMOS傳感器在一些特殊指紋的采集上因清晰度及噪聲等達不到要求而造成采集困難，加上圖像光路采集傳輸系統(tǒng)本身的噪聲干擾和非線性失真，對于指紋圖像，容易引起許多偽特征和丟失很多真實特征，嚴(yán)重影響后續(xù)的指紋比對和驗證，因此圖像去噪和指紋細(xì)節(jié)增強是指紋圖像處理中非常重要的環(huán)節(jié)。怎樣用低端CMOS傳感器得到符合FBI（Federal Bureau of Investigation），GA（公安）等標(biāo)準(zhǔn)的圖像，是目前指紋掃描儀中圖像處理算法要著重解決的問題。

1 基于PM方程的圖像增強算法研究

在傳統(tǒng)的指紋圖像噪聲濾波算法，如中值濾波，高斯濾波等［2］，主要將圖像的高頻成分濾除，由于指紋紋理的圖像信號相對背景來說，也分布在高頻區(qū)域，因此在對噪聲進行濾波的同時必定損害到指紋的特征信息。

為了既保留指紋的特征細(xì)節(jié)信息，又能夠提高圖像的對比度及分辨率，這里提出了一種改進的基于PM方程和USM的去噪增強算法，比較好的解決了提高圖像分辨率（或者MTF（modulation transfer funation））和去除噪聲的矛盾，整個算法流程主要包括3個部分，首先采用自適應(yīng)算法消除圖像的馬賽克效應(yīng)，然后采用優(yōu)化的USM（unsharp mask）銳化算法，同提高圖像邊緣的清晰度，結(jié)合改進的PM噪聲濾波算法，幾乎在不損害圖像指紋細(xì)節(jié)的前提下，有效地去除了指紋圖像存在的高斯噪聲和隨機噪聲。

1.1 去馬賽克算法

因為光照和CFA的特性，在采集的圖像中存在通道不均勻的現(xiàn)象，也就是馬賽克效應(yīng)，通常去除馬賽克的算法采用二次或者三次插值的算法［3-5］，即根據(jù)每個點周圍的4個或者16點對該點進行插值，這種算法通過對圖像平滑在解決通道不均勻方面能夠取到一定的效果，但是不可避免的會引起邊緣或者紋理的模糊，導(dǎo)致分辨率的下降，在這里中，采用了基于邊界的自適應(yīng)插值算法，即首先分別計算插值點的X方向和Y方向梯度，決定邊界的方向，然后沿著邊界方向進行插值參數(shù)的提取。參照圖1所示的布列模板，具體算法流程如下。

1）首先計算水平方向的梯度

圖1 CFA制式圖像模板Fig.1 The templet of CFA

2）然后計算垂直方向的梯度

3）進行插值

其中：ΔH為水平方向的梯度；ΔV為垂直方向的梯度；R，G，B變量分別表示不同通道中的分量。

1.2 優(yōu)化的USM圖像銳化

為了抑制噪聲，在采集原始圖像時候，增益不能調(diào)的過高。由于指紋掃描設(shè)備本身光學(xué)系統(tǒng)的畸變，采集的圖像分辨率對一些干性指紋，或者殘缺的指紋不能滿足要求，這里采用了USM 算法對圖像進行銳化，USM的銳化模型可以公式表達為

其中：i、j為像素的坐標(biāo)；f(i,j)為輸入圖像；y(i,j)為輸出圖像；而z(i,j)為校正信號；λ是用于控制圖像邊緣銳化強弱的松弛因子。z(i,j)一般可以通過下式獲取

這種算法的主要缺點是對噪音非常敏感，尤其在圖像變化比較緩慢的地方的噪音。另外對于圖像的邊緣（高對比度區(qū)域）會出現(xiàn)增強過頭的現(xiàn)象。因此這里采用類似于Photoshop的USM算子原理，用二維高斯函數(shù)g(i,j)［6］作模糊濾波器，得到虛化后的圖像image_g(i,j)。最終根據(jù)平滑后的圖像和原始圖像之間的差值大小決定是否對圖像進行銳化，圖像銳化如公式（7）所示，設(shè)，其中σ為高斯平滑半徑，得到

其中：參數(shù)Δv=f(i,j)-image_g(i,j)；amount控制銳化的強弱，過弱銳化效果不明顯，過強容易引起噪聲；Threshold參數(shù)用來判斷是否對該點進行銳化，偏大將忽略弱邊緣的銳化，偏小則容易對一些噪聲進行強化。

隨著σ半徑的增大，圖像的對比度逐漸變強，邊緣越發(fā)明顯，但在σ增大的同時，時間也是幾何增長。為了保證算法的實時性，這里提出了一個與窗口半徑無關(guān)的USM銳化算法，即用局部均值代替高斯平滑函數(shù)，同時因為局部均值可以用積分圖進行計算，在實際應(yīng)用中可以保證算法的實時性。設(shè)

其中：sum為求和運算。

得到局部均值

從公式（9）可看出，計算時間與窗口w大小無關(guān)，同時在實際應(yīng)用中，調(diào)節(jié)后的指紋圖像隨著σ半徑的增大并未出現(xiàn)明顯的噪音增強現(xiàn)象，效果比較理想。

1.3 改進的基于PM方程的圖像噪聲濾波算法

如前所示，傳統(tǒng)的噪聲濾波算法在濾除噪聲的同時，或多或少都可能會損壞圖像中的某些高頻信息，從而引起模糊邊緣和紋理效應(yīng)，特別是對指紋紋理這種非常敏感的信號。綜合考慮多種處理算法，在這里采用了基于PM方程的各向異性擴散偏微分方程濾波算法［7-8］，這個算法由Perona 和Malik提出，根據(jù)圖像不同方向上的梯度，在非線性尺度空間，確定擴散系數(shù)［9］，其方程表達式為

其中：u是圖像對應(yīng)某一個像素點的灰度值；▽u為梯度；div 是散度算子；t表示擴散時間即迭代的步長；c表示擴散系數(shù)，其方程如式（12）所示［10-11］，為圖像在某點梯度的非負(fù)單調(diào)遞減函數(shù)

其中：k表示擴散系數(shù)中的一個恒定參數(shù)，這個參數(shù)控制擴散強度。PM方程的離散化公式可表示成

其中：q表示圖像像素的位置；uq為圖像像素q處對應(yīng)的灰度值；▽uq,p=up-utq；λ表示擴散速率的一個常數(shù)參數(shù)；Nq表示像素q的鄰域，通常取圖像的左、右、上、下四鄰域，如果加上斜對角領(lǐng)域可以得到更好的效果，但影響算法實時性；t為實時迭代次數(shù)。以圖像在不同方向上的梯度的單調(diào)遞減函數(shù)作為各向異性擴散的擴散系數(shù)，使得在同質(zhì)區(qū)域內(nèi)部，它的擴散系數(shù)較大，可以有效濾除該區(qū)域內(nèi)的噪聲，而在邊緣區(qū)域因為灰度值變化劇烈，梯度較大，對應(yīng)的擴散系數(shù)較小，因此能較好地保留指紋圖像紋理的邊緣信息。

PM方程算法濾除噪聲和保留細(xì)節(jié)的性能較傳統(tǒng)算法有了較大的提升，但是隨著迭代的增加，邊緣被平滑的風(fēng)險增大，其邊緣保持的條件是 ||▽u＞k，在實際應(yīng)用中，因為不同區(qū)域噪聲強弱不同，這時對強噪聲區(qū)域，平滑區(qū)域，邊緣區(qū)域的參數(shù)應(yīng)當(dāng)自適應(yīng)調(diào)整，同時為了減少迭代的次數(shù)，在比較各種PM改進算法的基礎(chǔ)上，考慮實時性的要求，這里提出了基于局部歸一化方差的自適應(yīng)擴散系數(shù)

其中：σ20即為局部方差。在這里局部方差也采用積分圖快速算法來進行計算。具體公式這里不推導(dǎo)，只給出計算過程。設(shè)

可推導(dǎo)出局部平方灰度值

得出局部方差

從而得到歸一化的局部方差

從式（18）可推出在平滑區(qū)域，擴散系數(shù)趨于0，在強噪聲區(qū)域，擴散系數(shù)趨于1，滿足擴散系數(shù)的要求，其邊緣自適應(yīng)定位能力和迭代時間皆優(yōu)于前面的算法。

2 實驗結(jié)果分析

為檢驗本算法在實際應(yīng)用中的效果，這里采用FBI對指紋圖像認(rèn)證的Mitre公司的MTF檢測軟件，采用標(biāo)準(zhǔn)檢測模板和CTF評價系數(shù)來檢測本算法的有效性。從圖2可以看出，圖像經(jīng)過1.1提出的去馬賽克算法處理之后，基本沒有存在通道不均衡的現(xiàn)象，圖像的標(biāo)準(zhǔn)偏差大幅度降低。

圖2 去馬賽克算法處理結(jié)果Fig.2 The filtering effects of demosaicking

圖3 圖像經(jīng)USM和PDE處理結(jié)果Fig.3 The processing results of improving USM and PDE

圖3（a）圖像在經(jīng)過優(yōu)化的USM 銳化之后得到圖3（b），雖然邊緣對比度得到了明顯的增強，但是不可避免還是會引起部分噪聲，而在經(jīng)過改進PM濾波之后得到圖3（c），在低頻區(qū)域，噪聲幾乎可以忽略不計，并且迭代次數(shù)較原PM算法縮短了一倍，Mitre的MTF測試軟件得出CTF系數(shù)結(jié)果如表1所示，從結(jié)果可以看出，在經(jīng)過算法處理之后，各個對線頻率（Tar Freq）的指標(biāo)都有了明顯的提升。其中，cy/mm為每毫米的采樣頻率，Peak CTF為對比度函數(shù)的極大值。

表1 MITRE’S MTF Software 測試結(jié)果Tab.1 The testing results of MITRE’S MTF Software

[1] BLANKSBY A J,LOINAZ M J.Performance analysis of a color CMOS photograph image sensor[J].IEEE Trans Electron Devices,2000,47(1):55-64.

[2] 蔣先剛.基于Hessian特征的視網(wǎng)膜血管圖像的增強濾波算法[J].華東交通大學(xué)學(xué)報,2013,30(3):37-43.

[3] RAJEEV RAMANATH,WESLEY E SNYDER.Adaptive demosaicking[J].IEEE Trans Journal of Electronic Imaging, 2003,12(4):633-642.

[4] KAKARALA R, BAHARAV Z.Adaptive demosaicking with the principal vector method[J].IEEE Trans Consumer Electron,2002,48(4):932-937.

[5] 劉曉松,楊新,文俊.一種用于數(shù)字圖像傳感器的彩色插值算法[J].中國圖像圖形學(xué)報,2003,8A(5):516-523.

[6] PERONA P,MALIK J.Scale space and edge detection using ansotropic diffusion[J].IEEE Trans PAM I, 1990,12(7):629-639.

[7] ANDREA POLESEL, GIOVANNI RAMPONI.Image enhancement via adaptive unsharp masking[J].IEEE Trans Image Processing,2000,9(3):505-510.

[8] ILYA POLLAK, ALAN S WILLSKY.Image segmentation and edge enhancement with stabilized inverse diffusion equation[J].IEEE Transaction on Imae Processing,2000,9(2):256-266.

[9] 宋建軍,侯志強.基于人類視覺模型的各向異性擴散濾波方法[J].中國圖向圖形學(xué)報,2012,17(3):321-327.

[10] 陳卿,胡紹湘.基于邊緣擴散的醫(yī)學(xué)圖像非線性去噪算法[J].計算機應(yīng)用研究,2011,28(8):3178-3180.

[11] 蔣先剛,許倫倫.基于三維各向異性擴散的圖像平滑及三維重構(gòu)效果分析[J].華東交通大學(xué)學(xué)報,2010,27(3):78-82.