基于Gabor濾波和形態(tài)學變換的手臂靜脈線提取

彭曉光 王彪 唐超穎 蘇菡 陳曉騰

摘 要：提出了一種針對近紅外圖像的手臂靜脈提取算法。首先利用限制對比度自適應(yīng)直方圖均衡化（CLAHE）對近紅外圖像進行預(yù)處理，增強對比度。然后利用Gabor濾波器對增強后的圖像進行濾波，由Gabor濾波后的方向圖和能量圖獲得手臂靜脈位置。接著，對能量圖進行高低帽變換，然后將整幅圖像進行二值化，采用開閉運算處理二值化圖像，利用形態(tài)學骨架提取的方法提取手臂靜脈線，最后對骨架化后的圖像進行毛刺修剪，得到比較完整有效的靜脈線。實驗結(jié)果表明，本文提出的算法可以有效地保證靜脈線的連續(xù)性和完整性，有效地提高手臂靜脈的匹配率。

關(guān)鍵詞：靜脈提取;圖像增強;Gabor濾波;靜脈匹配;

中圖分類號：TP391

文獻標志碼：A

引 言

隨著計算機視覺的快速發(fā)展，生物特征識別的應(yīng)用也越來越廣泛，生物特征識別技術(shù)是利用人體固有的生理特征和行為特征進行身份識別。由于傳統(tǒng)身份識別技術(shù)不可靠、不安全等缺點，生物特征識別技術(shù)已經(jīng)成為國內(nèi)外較為熱門的研究方向。靜脈識別是通過近紅外相機獲得對應(yīng)靜脈部位近紅外圖像，利用靜脈的結(jié)構(gòu)特征完成對人的身份鑒別認證。與其他生物特征識別相比，靜脈屬于人體皮下組織的固有特征，基本不會發(fā)生改變，并且必須在活體下進行采集，具有難以偽造等優(yōu)點，因此越來越受到研究人員的重視。

靜脈識別首先要從近紅外圖像中獲得較為完整的靜脈結(jié)構(gòu)，然后提取靜脈完成身份識別，因此，從近紅外圖像中獲得較為完好的靜脈線結(jié)構(gòu)成為至關(guān)重要的工作，靜脈提取算法的優(yōu)劣直接關(guān)系到匹配識別效果的好壞，甚至對整個靜脈識別系統(tǒng)的性能產(chǎn)生重大影響。近年來對手部靜脈提取與識別算法的研究越來越多。文獻[1]構(gòu)造兩個代價函數(shù)：一個是基于靜脈測度和強度信息的頻率，另一個是基于區(qū)域信息，針對兩個不同的代價函數(shù)，利用最小代價路徑的方法提取靜脈線，運用在冠狀動脈造影圖像（ CTAT）;文獻[2]提出了一種基于梯度的靜脈模式分割算法，最終得到的是二值化模式的手背靜脈;文獻[3]使用基于局部曲率閾值的方法提取手背靜脈，利用形態(tài)學骨架化的提取靜脈線;文獻[4]首先使用CLAHE增強手背靜脈近紅外圖像對比度，然后利用形態(tài)學骨架化的方法提取手背靜脈;文獻[5]同樣采用CLAHE增強手背靜脈圖像，接著使用PDSVP（Palm Dorsa Subcutane-ous Vein Pattern）和形態(tài)學處理提取手背靜脈;文獻[6]采用對比度自適應(yīng)直方圖均衡化（CLAHE）對手背靜脈圖像增強后二值化，然后利用局部靜態(tài)閾值與Niblcak相結(jié)合的方法提取手背靜脈;文獻[7]結(jié)合手指靜脈紋路特征設(shè)計方向圖和方向濾波器進行圖像增強，然后對增強后的圖像提取手指靜脈紋路;文獻[8]利用靜脈灰度圖像靜脈所在位置灰度值與其紋路垂直方向上的灰度值成谷形的特征提取手指靜脈線。

目前大部分對靜脈識別的研究都是基于手指、手掌和手背區(qū)域，這些區(qū)域的皮膚薄，靜脈距離皮膚較近，相對比較容易提取，但是靜脈數(shù)量相對較少。本文的研究對象是手臂部位的靜脈特征，相對手指、手背和手掌區(qū)域，皮膚較厚，脂肪較多，提取相對較為困難，但是手臂部位具有更加豐富的靜脈信息，手臂靜脈的匹配更具有魯棒性。本文首先使用CLAHE方法對手臂近紅外圖像進行增強處理，然后使用Gabor進行濾波，利用濾波后的方向圖和能量圖獲得靜脈結(jié)構(gòu)，使用OTSU方法對靜脈結(jié)構(gòu)圖像進行二值化，最后采用形態(tài)學的處理方法，獲得較為完整的靜脈線。

1 手臂靜脈提取

1.1 手臂近紅外靜脈圖像預(yù)處理

本文采用的是新加坡南洋理工大學的近紅外手臂圖像數(shù)據(jù)庫，采用JAI-AD080CL相機拍攝，包含來自不同國家、不同人種的樣本，原始數(shù)據(jù)庫的手臂近紅外如圖1所示，圖中a和b膚色有一定差異，拍照時光照強度也有差別。靜脈隱藏于皮膚下層，從圖1中可以看出，靜脈線與皮膚的對比度不高，直接提取靜脈紋理，通過二值化的方法提取靜脈線，將會導(dǎo)致一部分靜脈信息的丟失，也有可能將背景信息誤識為靜脈結(jié)構(gòu)，這對后續(xù)的特征提取和識別造成巨大的影響，因此需要首先對靜脈圖像進行增強預(yù)處理。

圖像增強一般分為空間域和頻域增強，空域增強方法中直方圖均衡（HE）是實用有效的工具，HE以概率論為基礎(chǔ)，通過每一點灰度值概率的運算，把目標圖像中直方圖比較集中的分布改變?yōu)榫鶆蛑狈綀D分布，同時使輸出的圖像像素灰度的概率密度均勻分布。直方圖均衡分為全局法和自適應(yīng)法（ Adaptive HE，AHE），全局法在實際應(yīng)用中效果不佳。文獻[9]中提出了對比度受限的局部直方圖均衡法（Contrast-Limited Adaptive HE，CLA-HE），結(jié)合自適應(yīng)直方圖和對比度受限兩種方法，從整幅圖像的視覺效果出發(fā)，改善圖片質(zhì)量。這種方法考慮窗口內(nèi)像素直方圖，使圖像增強效果適應(yīng)性更好，其表達式為：

（1）對圖像中任一點，根據(jù)窗口大小W確定窗口區(qū)域，使用式（1）計算窗口的直方圖。

（2）對上述窗口內(nèi)的直方圖進行均衡化，實現(xiàn)對窗口中心像素的處理。

（3）移動窗口到下一個相鄰像素，重復(fù)上述過程，直至遍歷整幅圖像。

如圖2（a）所示，（從圖1經(jīng)過手臂掩膜處理，主要針對前臂所在區(qū)域進行處理，本文之后圖像均對主要的手臂區(qū)域進行處理）原始手臂近紅外圖像中靜脈雖然已經(jīng)顯現(xiàn)出來了，但與皮膚對比度并不強烈。經(jīng)過自適應(yīng)直方圖均衡化處理，處理結(jié)果如圖2（b）所示，靜脈與手臂皮膚已經(jīng)呈現(xiàn)明顯的差異，有利于下一步的靜脈提取。

1.2 Gabor濾波和增強

Gabor變換與人類視覺系統(tǒng)中的視覺刺激響應(yīng)非常相似，它在提取目標的局部空間和頻率域信息方面具有良好的特性，被廣泛應(yīng)用于生物醫(yī)學圖像處理，比如指紋、人臉、虹膜等，均得到比較理想的結(jié)果[10]。本文利用Gabor濾波器提取預(yù)處理圖像中的靜脈信息，由于預(yù)處理近紅外圖像中的靜脈為較為深色的線條，故采用含有16個不同方向和尺度的Gabor濾波器的實部進行處理[11]。下式為Gabor濾波器的空域數(shù)學表達式：

1.3 靜脈線提取

高帽（ top-hat）變換是原始圖像與圖像的開運算的結(jié)果之差，低帽（ bottom-hat）變換是圖像的閉運算結(jié)果與原始圖像之差，高低帽變換的定義[13]為：

高低帽變換重要用途是校正不均勻光照對圖像的影響，高帽變換適用于暗背景中的亮物體，而低帽變換則適用于相反的情況。因此高低帽變換可以用于對Gabor濾波增強后的圖像進行處理，高帽變換具有高通濾波特性，強調(diào)靜脈圖像的灰度峰值，增強圖像中的靜脈血管部分;而低帽變換可以求出圖像的灰度谷值，突出靜脈相互連接的目標間的界限，所以高低帽變換結(jié)合使用，可以使圖像前景和背景進一步拉伸，盡可能突顯出靜脈的所有部分，最后達到靜脈增強的作用。Gabor濾波增強后的靜脈加上高帽變換的結(jié)果，再減去低帽變換的結(jié)果，可以有效提高圖像的對比度，如式（9）所示。

整個處理過程及結(jié)果如圖4所示，對Gabor濾波增強后的靜脈圖像（a）進行高帽變換（b）和低帽變換（b），高帽變換使主要的靜脈血管亮度增加，次要的靜脈血管變暗，通過原圖加上高帽變換結(jié)果圖，再減去低帽變換結(jié)果圖，最終的處理結(jié)果（d）明顯可以看出整體的靜脈血管部分變亮，既突出了重要血管，又去除了大部分噪聲。

使用Otsu方法[14]對上述結(jié)果圖進行二值化處理。由于得到的靜脈脈絡(luò)上有小的空洞，所以對二值圖像進行形態(tài)學閉運算處理，可以平滑輪廓，有利于后續(xù)的形態(tài)學骨架提取，而且可以融合狹窄間斷和細長溝壑，并消除小的空洞。接著，對閉運算后的圖像進行開運算，消除圖像中小的孤島和尖刺。

為了提取靜脈線特征，采用形態(tài)學的骨架提取獲得單像素靜脈線。集合A的骨架[14]可以用腐蝕和開運算操作來表達，即骨架可以表示為：

2 實驗結(jié)果及分析

2.1 靜脈提取算法對比

首先將本文算法與直接利用Otsu算法進行靜脈提取的結(jié)果進行比較。圖5（a）是一幅近紅外手臂圖像，圖5（b）為經(jīng)過CLAHE處理后的圖像。圖5（c）是使用傳統(tǒng)的Otsu算法進行提取的結(jié)果?？梢姡唵蔚拈撝捣椒ㄖ荒馨咽直蹍^(qū)域分割出來，基本上無法分出靜脈區(qū)域，無法進行后續(xù)的靜脈提取。

將本文算法與另外三種常用的靜脈提取算法進行比較，結(jié)果如圖6所示，第一行四幅手臂近紅外圖像，第二行為采用局部自適應(yīng)閾值法對原始的NIR手臂圖像進行處理，然后提取的靜脈線，雖然提取出了較為完整的靜脈線，但是引入了過多的噪聲和毛刺，這對后續(xù)的匹配識別都有較大的不良影響。第三行為采用局部梯度極值法提取的靜脈線，只提取出了主要的靜脈線，漏掉了許多重要的靜脈信息，另外靜脈的連續(xù)性較差。第四行為采用Ni-black算法對處理后的圖像進行分割，靜脈的連續(xù)性也比較差，毛刺也比較多，但總體比前兩種效果要好。最后一行為本文提出算法的實驗結(jié)果，與前邊三組實驗對比，可以明顯看出，提起出的靜脈線明顯比較連續(xù)。由于經(jīng)過形態(tài)學的處理，融合了二值圖像上一些細長的溝壑，平滑了輪廓，消除了一部分靜脈線上的空洞，在骨架提取時這些毛刺也相對較少。

2.2 靜脈線修復(fù)

雖然經(jīng)過本文提取算法已經(jīng)提取出較為完整的靜脈線，但是仍然存在少量的毛刺，這對后續(xù)的匹配會產(chǎn)生不良的影響，因此需要刪除靜脈線狀圖中的毛刺[15]。二值靜脈線的端點定義為某點8鄰域像素值的和為1的點，分叉點定義為某點8鄰域像素值的和大于或者等于3的點，端點到端點或者分叉點的線長度小于閾值T時，這段曲線定義為毛刺，需要對這些毛刺進行修剪。如圖7（b）所示，黑色圓圈標記的為端點，紅色圓圈標記的分叉點。本文經(jīng)過多次實驗統(tǒng)計，把毛刺閾值T設(shè)為15個像素，可以達到較好的效果，既刪除了多余的毛刺，又保留了重要的靜脈線部分。經(jīng)過毛刺修剪后的靜脈圖像如圖7（c）所示，與圖7（a）相比可以明顯看出，靜脈線上的許多毛刺已經(jīng)被刪除，這對后續(xù)的特征提取非常有利。

2.3 匹配實驗

針對100人的手臂近紅外圖像數(shù)據(jù)庫進行了提取和匹配實驗，測試人員包括女性8人，男性92人，年齡范圍為15 - 70周歲，來自8個不同的國家，每一個人采集兩次圖像，兩次拍攝的光照與拍攝角度均有變化，平均時間間隔為兩周，圖像采集距離約為2米，圖像尺寸為1024* 768。匹配實驗中把第一周期采集的100幅手臂近紅外圖像作為參考數(shù)據(jù)庫，把第二周期采集的圖像作為待匹配圖像。采用文獻[16]的匹配算法，分別對自適應(yīng)閾值、局部梯度閾值和本文算法提取的手臂靜脈圖像進行匹配實驗，最終的匹配結(jié)果如圖8所示，紅色曲線為本文提取算法的識別率，可以明顯看出，優(yōu)于其他兩種算法。

結(jié)束語

針對手臂靜脈這一較為新穎的生物特征，采用CLAHE對手臂近紅外圖像進行預(yù)處理，并進行Gabor濾波和增強，然后使用一系列形態(tài)的高低帽變換、開閉運算和骨架提取操作，提取手臂靜脈曲線。運用本文算法進行了靜脈的提取，通過匹配實驗對比驗證，獲得了比其他提取方法相對較高的識別率，證明本文提出算法的有效性。

參考文獻

[1]METZ C T，SCHAAP M，WEUSTINK A C，et aL Coronary cen-terline extraction from CT coronary angiography images using aminimum cost path approach[J]. Medical physics，2009， 36（ 12）：5568- 5579.

[2]WANG Y，ZHANG K，SHARK L K. Personal identificationbased on multiple keypoint sets of dorsal hand vein images[J]. IET Biometrics，2014，3（4）： 234-245.

[3]LAJEVARDI S M，ARAKALA A，DAVIS S，et al. Handvein authentication using biometric graph matching [J]. IETBiometrics，2014，3（4）： 302-313.

[4] JOARDAR S，CHATTERJEE A，RAKSHIT A.A real-timepalm dorsa subcutaneous vein pattern recognition systemusing collaborative representation-based classification [J].IEEE Transactions on Instrumentation and Measurement，2015，64（4）：959-966.

[5] JOARDAR S，CHATTERJEE A，RAKSHIT A. Real-timeNIR imaging of Palm Dorsa subcutaneous vein pattern basedbiometrics： An SRC based approach [J]. IEEE Instrumenta-tion& Measurement Magazine，2016，19（2）： 13-19.

[6] 喻恒，鄭均輝，基于CLAHE和改進的NIBLACK算法的手背靜脈提取方法[J].計算機與現(xiàn)代化，2015，（4）：74-77.

[7]王科俊，馬慧，李雪峰，等，使用方向濾波技術(shù)的手指靜脈紋路提取方法[J].中國圖象圖形學報，2011， 16（7）：1206 -1212.

[8]余成波，秦華鋒，手指靜脈圖像特征提取算法的研究[J].計算機工程與應(yīng)用，2008，（24） ：175-177.

[9] CROMARTIE R，PIZER S M. Structurrsensitive adaptivecontrast enhancement methods and their evaluation [J].Image and Vision Computing，1993，11（8）： 460-467.

[10] DAUGMAN J G.Complete discrete 2-D Gabor transformsby neural networks for image analysis and compression [J].IEEE Transactions on Acoustics， Speech， and Signal Pro-cessing，1988，36（7）： 1169-1179.

[11]唐超穎，一種基于近紅外一可見光映射的快捷靜脈提取算法[J].東南大學學報：自然科學版，2016，46（1）：35-41.

[12] ZHANG H，TANG C，KONG A W K，et al. Matching veinpatterns from color images for forensic investigation[C]//Biometrics： Theory， Applications and Systems（ BTAS），2012 IEEE Fifth International Conference on. IEEE，2012：77 - 84.

[13]索愛黎，王小鵬，形態(tài)學圖像分析原理與應(yīng)用[M].北京：清華大學出版社，2008.

[14] 阮秋琦.數(shù)字圖像處理學[M].北京：電子工業(yè)出版社，2011.

[15] 趙珊，手臂靜脈紅外圖像特征檢測與匹配算法研究[D].南京航空航天大學，2016.

[16] TANG C，ZHANG H，KONG A W，et al. Using multiplemodels to uncover blood vessel patterns in color images forforensic analysis [J]. Information Fusion，2016， 32（ 32）： 26 -39.