鄢圣玥 陳森森

摘? 要：針對臨床中大多數(shù)手部紅外靜脈圖像增強(qiáng)算法存在成像質(zhì)量差、成像過程中容易斷點而無法連續(xù)高質(zhì)量成像的問題，文章在已有AGC和Frangi算法的基礎(chǔ)上，設(shè)計一種基于雙波長的手部紅外靜脈圖像增強(qiáng)算法，該算法能夠有效解決單一的AGC和Frangi算法在圖像成像上不連續(xù)、圖像成像質(zhì)量差的問題。在對不同波長進(jìn)行測試后，采用波長為850 nm+940 nm的紅外光照射獲得原始靜脈圖像，可提高紅外靜脈圖像增強(qiáng)算法的性能，對提高紅外靜脈圖像產(chǎn)品的臨床應(yīng)用效果具有重要價值。

關(guān)鍵詞：雙波長；手部靜脈；圖像增強(qiáng)；紅外成像

中圖分類號：TP391.4? 文獻(xiàn)標(biāo)識碼：A? 文章編號：2096-4706（2023）05-0082-04

Research on Hand Infrared Vein Image Enhancement Algorithm Based on Dual-wavelength

YAN Shengyue， CHEN Sensen

（Wuhan University of Technology， Wuhan? 430070， China）

Abstract： In view of the problems of poor imaging quality， it is easy to break points in the imaging process and cannot continuously high-quality imaging in most of the hand infrared vein image enhancement algorithms in clinical practice， on the basis of the existing AGC and Frangi algorithms， this paper designs a hand infrared vein image enhancement algorithm based on dual-wavelength. This algorithm can effectively solve the problems of image's imaging discontinuity and poor imaging quality of a single AGC and Frangi algorithm. After testing different wavelengths， the original vein image is obtained by infrared irradiation with wavelength of 850 nm+940 nm， which can improve the performance of infrared vein image enhancement algorithm and has important value for improving the clinical application effect of infrared vein image products.

Keywords： dual-wavelength; hand vein; image enhancement; infrared imaging

0? 引? 言

靜脈紅外圖像處理技術(shù)在臨床醫(yī)療中應(yīng)用廣泛，可以用于觀察患者的靜脈血管的狀態(tài)，輔助穿刺和進(jìn)行相關(guān)疾病的診斷及排查，同時也是基于靜脈AI診斷相關(guān)技術(shù)的基礎(chǔ)。

現(xiàn)有的靜脈紅外圖像增強(qiáng)方法主要是從提升圖像對比度，增強(qiáng)圖像紋理特征，基于圖像融合的增強(qiáng)三個方面來進(jìn)行設(shè)計。其中提升圖像對比度主要為各種改進(jìn)的直方圖均衡法；增強(qiáng)圖像紋理特征主要基于一系列濾波技術(shù)；基于圖像融合的增強(qiáng)方法主要為將不同方法處理的多張圖像進(jìn)行融合，由于不同方法處理的圖像在某些特征上效果明顯，通過圖像融合可以把圖像的增強(qiáng)效果明顯區(qū)域疊加，同時減弱單張圖像的噪聲部分，進(jìn)而得到效果更佳的圖像。以上靜脈圖像增強(qiáng)方法，對不同部位的靜脈圖像處理效果差別較大，局限于身體某一部位的血管，導(dǎo)致普適性較差，這給相關(guān)產(chǎn)品的實際的臨床應(yīng)用帶來了巨大的挑戰(zhàn)，如何進(jìn)一步提高算法的穩(wěn)定性和成像質(zhì)量，成了當(dāng)前該領(lǐng)域需要解決的主要問題。

在上述背景下，本文提出了一種融合AGC和Frangi的算法，能夠有效地改善當(dāng)前手部紅外靜脈圖像增強(qiáng)算法存在的斷點和噪聲問題，同時使用雙波長紅外光照射，有效地提高手部紅外靜脈圖像成像質(zhì)量，具有非常重要的意義和價值。

1? 算法設(shè)計思路

通過光源的對比，可以發(fā)現(xiàn)850 nm+940 nm的雙波長組合能獲得最佳的靜脈圖像，通過前文對AGC和Frangi兩種算法的分析，其中AGC算法能夠有效地處理照度不均勻的紅外圖像，達(dá)到降噪的目的，F(xiàn)rangi能夠?qū)D像中的特征數(shù)據(jù)進(jìn)行識別從而檢測靜脈血管，融合這兩種算法，將提高手部紅外靜脈圖像的成像效果，解決成像后的圖像斷點及噪聲過大的問題?；谠撍悸?，本文構(gòu)建了融合了AGC和Frangi算法的雙波長的手部紅外靜脈圖像增強(qiáng)算法，其實現(xiàn)流程如圖1所示。

基于上述流程，在對手部紅外靜脈圖像進(jìn)行處理的過程中，首先采集雙波長紅外光850 nm+940 nm照射下的手部靜脈圖像，對初始圖像進(jìn)行引導(dǎo)濾波，起到銳化圖像邊緣的作用，再通過AGC（自適應(yīng)伽馬校正）拉伸圖像的對比度，接著對處理后的圖像進(jìn)行多尺度高斯濾波降噪處理，然后對處理后圖像運(yùn)用Frangi算法進(jìn)行血管識別，提取出管狀的血管輪廓，最后對圖像依次進(jìn)行局部二值化分割和開運(yùn)算降噪，得到視覺效果明顯的分割后的血管形狀。在整個圖像處理中，將AGC算法與Frangi算法創(chuàng)新性的結(jié)合，二者結(jié)合可以有效減少Frangi處理后產(chǎn)生的斷點。在進(jìn)行Frangi算法處理后，圖像存在很多斷點，因此在進(jìn)行Frangi算法處理前，需要對圖像進(jìn)行預(yù)處理，除了本發(fā)明使用的引導(dǎo)濾波、多尺度高斯濾波外，輔以AGC算法，能夠更加明顯地減少Frangi處理后產(chǎn)生的斷點。

2? 算法詳細(xì)設(shè)計實現(xiàn)

對紅外靜脈圖像增強(qiáng)算法的實現(xiàn)過程進(jìn)行分析，其核心步驟詳細(xì)設(shè)計如下。

2.1? 圖像獲取及ROI提取

使用850 nm+940 nm組合雙波長近紅外光照射，獲得手部靜脈圖像，對圖像進(jìn)行ROI感興趣區(qū)域提取，獲取大小為320×240的圖像區(qū)域，采集的圖像中不僅含有手部靜脈的感興趣區(qū)域，還含有背景等冗余信息，會增加不必要的計算量，進(jìn)而降低處理速度。對圖像進(jìn)行ROI區(qū)域提取可以簡單有效的解決這一問題。

2.2? 引導(dǎo)濾波處理

在完成圖像預(yù)處理及部分區(qū)域提取后，還未去除噪聲，本文通過對ROI區(qū)域進(jìn)行引導(dǎo)濾波，以解決其識別血管存在的噪聲的問題，使用輸入圖像作為引導(dǎo)圖像，達(dá)到保持邊緣的效果。具體原理為：

其中? 為經(jīng)濾波后的圖像，Gi（x， y）為引導(dǎo)圖像，k為引導(dǎo)圖像某一點的像素，ωk為以點k為中心、半徑為r的局部方形窗口，aK（x， y）和bK（x， y） 為引導(dǎo)。

圖像在ωk內(nèi)的線性系數(shù)。提取得到的圖像存在光照不均勻、對比度低等問題，因此需要進(jìn)一步處理，運(yùn)用引導(dǎo)濾波增強(qiáng)圖像邊。引導(dǎo)濾波的過程中保留了圖像的邊緣信息，同時抑制了圖像的噪聲，對圖像有增強(qiáng)效果，處理后的圖像邊緣部分加深，靜脈血管更加明顯，圖像對比度有提升，然而部分噪聲的邊緣也一并放大，所以需要更進(jìn)一步處理。

2.3? 圖像色彩空間轉(zhuǎn)換及AGC算法

在ROI區(qū)域引導(dǎo)濾波處理后的圖像，運(yùn)用自適應(yīng)伽馬校正（AGC）進(jìn)行對比度拉伸。將圖像色彩空間轉(zhuǎn)換至HSV，因為其H分量最為清晰，故提取H分量，對H分量進(jìn)行自適應(yīng)伽馬校正，之后進(jìn)行通道融合，其中伽馬值根據(jù)圖像特點確定?？紤]到圖像受到的光照不均勻，傳統(tǒng)的灰度拉伸方法會過曝或是過暗。因而利用自適應(yīng)的方法來處理圖像。經(jīng)過算法處理后的血管部分與背景有著更加明顯的對比度，便于后續(xù)處理。

2.4? 多尺度高斯濾波處理

對色彩空間轉(zhuǎn)換后的圖像進(jìn)行多尺度高斯濾波，經(jīng)測試驗證得到合適的高斯核參數(shù)以及相應(yīng)權(quán)重。本文選擇三個高斯核對圖像進(jìn)行加權(quán)處理。三個高斯核的σk值分別選取15、80、200，對應(yīng)權(quán)重為1/2、1/4、1/4。具體公式為：

在上述模型中，G（x， y， σ）為高斯函數(shù)，x、y、σ為高斯函數(shù)的參數(shù)。

其中P（x， y）為輸出圖像P在像素點（x， y）處的灰度值，G（x， y， 15）為σ值分為15的高斯函數(shù)，G（x， y， 80）為σ值為80的高斯函數(shù)，G（x， y， 200）為σ值為200的高斯函數(shù)，s（x， y）為輸入源圖像s在像素點（x， y）處的灰度值。

基于上述模型，圖像色彩空間轉(zhuǎn)換后得到對比度較高三通道的圖像，隨后采用多尺度高斯濾波的算法，分別對三個通道進(jìn)行加權(quán)的多尺度高斯濾波，不斷調(diào)試找到最好的去噪效果，從而得到噪點較少的圖像，將分別處理后的三個通道的圖像進(jìn)行融合。

2.5? 利用Frangi二維濾波器對血管特征進(jìn)行提取

在完成圖像融合后，利用Frangi二維濾波器，對血管特征進(jìn)行提取，從而區(qū)分血管和噪點，增強(qiáng)圖像中的靜脈信息。具體操作是通過對圖像中的像素點求解其Hessian矩陣特征值，即λ1和λ2。Hessian矩陣特征值和對應(yīng)的特征向量分別代表該點處沿某一方向上圖形曲率大小和方向。圖像中像素點（x， y）的Hessian矩陣H為：

其中?2f（x，y）/（?x2）為圖像關(guān)于x的二階偏導(dǎo)數(shù)，?2f（x，y）/（?y2）為圖像關(guān)于y的二階偏導(dǎo)數(shù)，?2f（x，y）/（?x?y）為圖像關(guān)于x、y的二階混合偏導(dǎo)數(shù)。

可以將Hessian矩陣中的元素求解，函數(shù)f （x， y）的導(dǎo)數(shù)等于函數(shù)f （x， y）與高斯函數(shù)的導(dǎo)數(shù)卷積。即，

其中f （x， y）為圖像f在點（x， y）處的灰度值，其中G（x， y）為高斯函數(shù)，?2f（x，y）/（?x2）為高斯函數(shù)關(guān)于x的二階偏導(dǎo)數(shù)，?2f（x，y）/（?y2）為高斯函數(shù)關(guān)于y的二階偏導(dǎo)數(shù)，?2f（x，y）/（?x?y）為高斯函數(shù)關(guān)于x、y的二階混合偏導(dǎo)數(shù)。經(jīng)上述模型處理后得到圖像的海森矩陣，進(jìn)而可以求得海森矩陣的特征值λ1和λ2。

根據(jù)上述模型和方法，根據(jù)求得的特征值λ1和λ2，本文進(jìn)一步將圖像形態(tài)大致分為三類，具體為：

（1）背景，灰度分布較均勻。任意方向上特征值λ1和λ2都較小。

（2）孤立的點，它在任意方向上的特征值λ1和λ2都很大。

（3）血管，特征值λ1、λ2一大一小。

經(jīng)過該算法處理后的圖像能將血管部分與背景區(qū)分，形成管狀結(jié)構(gòu)，進(jìn)一步濾去了非血管形態(tài)的噪聲。處理后的圖像還存在背景與血管結(jié)構(gòu)區(qū)分度不明顯的問題，一些噪聲被轉(zhuǎn)換為較小的管狀結(jié)構(gòu)，這些都需要進(jìn)一步處理來去除。進(jìn)行局部閾值的二值化分割（NiBlack）處理，對像素點選取r×r的鄰域，計算r×r鄰域內(nèi)的均值和方差，并結(jié)合修正值得到圖像分割的閾值，計算公式為：

T（x， y）=m（x， y）+k*s（x， y）

其中，T（x， y）為所求局部閾值，m（x， y）為局部均值，s（x， y）為局部方差。根據(jù)測試參數(shù)選取r=20，k=-0.2。分割后的圖像，血管與背景黑白分明，比較好辨認(rèn)，較原圖像有著很好的提高效果。由于同時當(dāng)進(jìn)行圖像遍歷時，如果域r×r范圍內(nèi)都是背景，經(jīng)NiBlack計算后必有一部分被確定為目標(biāo)，產(chǎn)生偽噪聲。對圖像進(jìn)行腐蝕膨脹等操作以去除偽噪聲，達(dá)到了對圖像進(jìn)行優(yōu)化的目的。最后將處理后的圖像輸出完成。

以上就是整個算法的詳細(xì)設(shè)計實現(xiàn)，基于本文的算法實現(xiàn)了AGC和Frangi算法的融合設(shè)計，能夠有效地提高算法的圖像處理效果，提取ROI區(qū)域降低圖像實時處理的速度，提高輸出圖像的平穩(wěn)性和清晰度。

3? 算法實驗仿真分析

基于本文設(shè)計的算法，利用OpenCv視覺處理庫對其進(jìn)行了實驗仿真分析，如圖2所示，是采集到的在三種不同波長紅外光照射下的原始圖像。經(jīng)過分析可知，其中850 nm+940 nm組合光照下獲得的圖像質(zhì)量最佳，利于后續(xù)處理。

圖3是利用本文算法實現(xiàn)的濾波處理后的結(jié)果，通過比較可以發(fā)現(xiàn)，730 nm+850 nm雙波長組合下圖像的細(xì)節(jié)不夠清晰，730 nm+940 nm雙波長組合下的圖像會存在照度不均勻所帶來的陰影區(qū)域，綜合比較后得出850 nm+940 nm雙波長組合更適合本文算法處理。同時經(jīng)過濾波處理后的圖片血管與周圍組織的對比更加清晰明顯。

進(jìn)一步對手部靜脈的圖像進(jìn)行增強(qiáng)處理，在提取的ROI區(qū)域，進(jìn)一步進(jìn)行靜脈圖像的分割以及去噪，其中圖

4（a）是本文選取的850 nm+940 nm雙波長照射下的圖像，圖4（b）中是通過本文算法生成的靜脈增強(qiáng)圖像的效果。從上述的結(jié)果來看，基于本文算法能夠很好地實現(xiàn)對紅外靜脈圖像進(jìn)行增強(qiáng)，將原本肉眼難以觀察到的靜脈變得清晰。有效地去除掉周邊毛細(xì)血管和圖像噪聲，形成清晰的靜脈圖像。

4? 結(jié)? 論

手部紅外靜脈圖像增強(qiáng)算法直接關(guān)系到相關(guān)靜脈穿刺輔助醫(yī)療產(chǎn)品的準(zhǔn)確性和精度，在臨床上使用相關(guān)的靜脈穿刺輔助設(shè)備，容易出現(xiàn)較大概率的錯誤穿刺現(xiàn)象。在這一背景下，本文從紅外靜脈圖像增強(qiáng)算法入手，通過對已有的AGC和Frangi算法進(jìn)行深入研究，設(shè)計了一種基于雙波長的手部紅外靜脈圖像增強(qiáng)算法，通過該算法能夠有效地解決單一的AGC和Frangi算法在圖像成像上不連續(xù)，圖像成像質(zhì)量差的問題，能夠有效地提取靜脈圖像，對提高靜脈穿刺輔助設(shè)備的性能具有非常重要的意義。

參考文獻(xiàn)：

[1] 張大鵬，盧光明，郭振華.基于手掌靜脈與掌紋的近紅外成像設(shè)備與身份識別方法：CN101251889 [P].2008-08-27.

[2] 孟琳，劉靜，曹慧，等.基于Frangi濾波器和Otsu視網(wǎng)膜血管分割 [J].激光與光電子學(xué)進(jìn)展，2019，56（18）：127-133.

[3] 王依人，鄧國慶，夏營威，等.基于方向可調(diào)濾波器的血管圖像增強(qiáng)算法 [J].系統(tǒng)仿真學(xué)報，2018，30（6）：2095-2101.

[4] 黃琨，陳書旺.基于紅外成像的輸液用血管定位系統(tǒng)研究 [J].計算機(jī)與數(shù)字工程，2018，46（5）：1025-1028.

[5] 彭曉光，王彪，唐超穎，等.基于Gabor濾波和形態(tài)學(xué)變換的手臂靜脈線提取 [J].計算技術(shù)與自動化，2018，37（1）：121-125+151.

[6] 岳林，萬新軍，張晨皓，等.復(fù)雜背景下紅外靜脈圖像的分割與增強(qiáng) [J].電子科技，2017，30（3）：118-120.

[7] RAO H，ZHANG P L，SUN C S. Contrast enhancement for the infrared vein image of leg based on the optical angular spectrum theory [J].Signal Image and Video Processing，2017，11（3）：423-429.

[8] 曹偉，王華彬，石軍，等.基于邊緣檢測加權(quán)引導(dǎo)濾波的指靜脈圖像增強(qiáng)算法 [J].激光與光電子學(xué)進(jìn)展，2017，54（2）：172-180.

[9] 嚴(yán)嬌嬌，種蘭祥，李婷.一種人手背靜脈特征識別方法 [J].計算機(jī)科學(xué)，2018，45（S1）：206-209.

作者簡介：鄢圣玥（2002.05—），女，漢族，湖北荊州人，本科在讀，研究方向：電子科學(xué)與技術(shù)。

收稿日期：2022-11-01

基金項目：國家大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計劃（S202210497296）