太赫茲成像檢測算法研究

李建蕊 李九生 史葉欣

摘要：由于太赫茲波所處波段位置特殊性以及現階段太赫茲成像系統性能的限制等，太赫茲波成像質量低，無法滿足可視化效果，限制了其發展和應用。結合太赫茲圖像模糊特征和灰度信息提出利用灰度特征對太赫茲圖像進行圖像分割來提高太赫茲圖像質量，抑制太赫茲圖像背景噪聲，保留太赫茲圖像目標重要信息，實現太赫茲成像目標檢測。實驗結果與其他太赫茲圖像處理方法結果對比表明，基于灰度特征的圖像處理算法可以提高圖像清晰度和對比度，實現精確分割，為太赫茲成像在安全檢查和醫學成像等應用中實現快速檢測和提取目標奠定基礎。

關鍵詞：太赫茲波;成像檢測;灰度分析;圖像分割;背景去噪;實驗分析

中圖分類號：TN911.73-34;0436.2

文獻標識碼：A

文章編號：l004-373X（ 2019） 24-0071-04

0 引 言

太赫茲在電磁波中處于特殊位置，但相對于其他波段，太赫茲波的研究起步比較晚，一個重要原因是很難產生存在于遠紅外與微波之間的太赫茲波，并且缺乏相應太赫茲波的檢測手段。隨著科學技術領域的快速發展，太赫茲源和探測器也得到了相應發展。近幾年，太赫茲成像技術在眾多領域均得到重視和應用，通過搭建太赫茲成像系統以實現太赫茲成像和檢查，如醫藥工程領域、工業工程領域以及安全檢查領域[1-5]。太赫茲成像所得太赫茲圖像的質量會受到太赫茲成像系統硬件和外界環境的影響和干擾，使太赫茲圖像出現模糊和對比度低等問題，甚至實驗結果使人眼無法在視覺上辨別目標物體的特征信息。對應的硬件和技術問題，如成像系統中產生太赫茲波的太赫茲源會隨機地出現波動使實驗樣本上的部分灰度特征在讀取過程中出現誤差，使所得太赫茲圖像的灰度分布出現抖動。這些并不是實驗樣本本來的灰度信息，這是由于太赫茲源不穩定造成的。圖像中出現“拖尾和斷點”的現象，是由于太赫茲波探測器所造成的。目前太赫茲成像系統所使用的探測器多為熱釋電型傳感器，雖然熱釋電型傳感器的響應范圍大和性價比高，但響應速度較慢，而響應元本身也存在著噪聲[6]。另外，該類探測器除了對太赫茲波有響應外，也易對環境中的紅外線產生響應，在太赫茲圖像上產生高斯噪聲等一系列隨機噪聲，太赫茲圖像中的條紋干擾噪聲，正是由太赫茲波探測器內部元件之間是相互干擾所致[7-8]。此外，在太赫茲成像實驗過程中，如二維成像平移臺機械移動的抖動和不均勻，數據采集卡的電子噪聲等多種因素導致了太赫茲圖像灰度信息重疊和錯位[9-10]。本文采用自設太赫茲成像系統對隱藏在紙箱中的衣物包裹金屬手槍進行透射式成像。根據實驗所獲得的太赫茲圖像存在的各種噪聲和干擾，假設太赫茲圖像的像素與像素之間的關系是相互獨立的，根據恒定背景圖像的強度分析，太赫茲圖像的噪聲可以建模為高斯噪聲。由于實驗所得太赫茲圖像具有分段平滑的特點，且沒有明顯的紋理，可以選取各向異性擴散濾波、雙邊濾波和低通濾波等方法對太赫茲圖像進行初步的平滑處理以削弱高頻噪聲對后期圖像處理的影響。太赫茲圖像中目標對象與背景噪聲的產生原因不同，在對太赫茲圖像進行灰度直方分析過程中，兩者所屬的像素類別會有所不同，根據太赫茲圖像像素類的變化，實現噪聲與目標的分割，并檢測提取出目標對象的邊界，從而準確檢測出目標對象和有用信息。

1 基本原理

對于存在背景噪聲的太赫茲圖像，先對太赫茲圖像進行一次平滑濾波，以減少高頻噪聲對后期太赫茲圖像檢測提取造成影響。為了保護太赫茲圖像中的局部信息包括邊界和紋理，采用低通濾波器對太赫茲進行首次濾波處理。初步濾波后，對太赫茲圖像進行全局信息和局部信息提取，找到太赫茲圖像的全局峰值和像素類均值，利用斜截差法獲得目標檢測的候選閾值，并通過太赫茲圖像中目標與背景灰度特征的差異來篩選出最終合適的閾值對太赫茲圖像進行分割，利用均衡化增強消除局部信息的非均勻特性。然后，提取圖像邊界，獲得無背景噪聲的太赫茲圖像。基于灰度特征的太赫茲成像檢測算法框圖如圖1所示。

圖像直方圖分布繪制了每個灰度值的像素數，設太赫茲圖像為f（x），歸一化灰度直方為p（x）。

2 實驗結果與分析

本文采用自設太赫茲成像系統如圖2所示，其為基于BWO連續太赫茲成像系統。以返波管振蕩器（BWO）作為太赫茲波輻射源搭建一個連續太赫茲波點掃描透射成像系統，并使用多線程控制技術控制鎖相放大器和二維成像平移臺。由于太赫茲源光斑具有點擴散特性，為了使太赫茲成像系統所得太赫茲圖像具有更完整的目標信息，本實驗采用掃描步長小于成像系統聚焦光斑尺寸的方式對樣品物體進行太赫茲透射成像，掃描步長為0.5 mm。在指定電壓下，返波管輻射出強度最大的連續太赫茲波（頻率達到0.71 THz），經斬波器和衰減片后由兩個聚乙烯透鏡進行校準和匯聚，使太赫茲波聚焦到樣本上，使其攜帶樣本信息。樣本置于可沿x軸和y軸掃描的線性步進二維平移臺上，平移臺載物根據預先設計好的掃描區域、步長和步速進行二維平移運動。該掃描區域必須使成像樣本完全置于光路中，且聚焦后的焦斑與樣本必須在同一個平面上。穿過樣本的光束，經過另一對透鏡的準直并聚焦到探測器上，探測到的光強信號由放大器放大并送到計算機進行保存，后續利用Matlab軟件提取數據以還原太赫茲圖像。為提高信噪比，本實驗在溫度保持298 K，濕度低于1%的條件下，對隱藏在紙箱中由衣物包裹的金屬手槍進行透射式成像，成像系統的分辨率為480 μmx360μm，如圖3所示。圖3a）為隱藏在紙箱中的待測金屬手槍的光學圖片，圖3b）為經過透射成像后所得到手槍的太赫茲圖像，圖3c）為太赫茲圖像的灰度直方圖，圖3d）為歸一化后的灰度直方分布圖。

太赫茲圖像的灰度直方圖可以直觀地反應太赫茲圖像的像素分布和局部與全局信息。通過分析其灰度直方圖，可以更快速地找到分割太赫茲圖像方法。由圖3d）可以看出，太赫茲原圖像中，目標的灰度級主要集中在50-100之間，而背景噪聲則主要集中在190-230之間。因此圖像邊界的灰度級就可能分布在150-190或者100-150之間。為后期快速從候選閾值中篩選出合適閾值創造了條件。圖4為對太赫茲圖像歸一化灰度直方分布進行分析所得到的相關曲線。圖中綠色曲線為太赫茲圖像的歸一化灰度直方分布曲線;藍色曲線為經過初次低通濾波后得到的歸一化灰度直方分布曲線;紅色曲線為全局峰值曲線，其波峰對應為太赫茲圖像的全局峰值，每個波峰代表一個像素類，那么曲線對應的橫坐標即為像素類的均值;黃色曲線為檢測目標所需的候選閾值曲線。本文選擇閾值分別為77，112，182，將太赫茲圖像進行像素分類標記。圖5a）為太赫茲原圖像進行初次濾波處理后的結果，圖5b）為標記后的太赫茲圖像，利用所選閾值對太赫茲圖像進行標記，可以看出太赫茲圖像在標記過程中局部信息的標記是不均勻的。此時利用均衡化，消除不均勻現象，并利用sobel邊界提取算子進行邊界提取，圖5c）所示為本文提出算法處理后的最終結果。

本文采用小波重構分解作為對比，對太赫茲圖像進行處理。圖6a）為太赫茲原圖像，圖6b）為平移不變小波分解后高頻和低頻信息分布，圖6c）為經過維納濾波后，進行小波重構得到的太赫茲圖像。雖然小波分解可以分解出不同頻率的信息，但是由于太赫茲圖像中存在背景噪聲，使得太赫茲圖像的特征沒有明顯顯現出來。從圖5和圖6對比發現，本文所提出的方法可以很好地分割背景和目標物體，并成功識別出測試樣本的輪廓，增強了可視化效果，精確地提取出了目標物體。

3 結語

本文利用自設的太赫茲成像系統對隱藏在紙箱中衣物包裹的實驗樣品進行透射式成像，對所獲得的太赫茲圖像進行灰度特征分析，并與小波分解重構方法獲得的結果進行比較。研究結果表明，本文算法一方面可以有效抑制太赫茲圖像的隨機噪聲和條紋噪聲等復雜噪聲，提高圖像清晰度和對比度;另一方面，根據太赫茲圖像本身的灰度強度分布進行分類，實現可選多閾值的精確分割，根據全局峰值對太赫茲圖像進行標記分類，去除背景噪聲保留太赫茲目標全局信息和細節信息。采用邊界提取算子進一步精準地識別和提取出了太赫茲圖像目標邊界，增強可視化效果，達到了預期目的，為太赫茲成像技術在安全檢測和醫學檢測等應用中實現快速檢測和提取目標奠定了基礎。

注：本文通訊作者為李九生。

參考文獻

[1] DENG F，TANG Q，ZENG G，et al.Effectiveness of digital in—frared thermal imaging in detecting lower extremity deep ve-nous thrombosis [J]. Medical physics， 2015， 42 （5）： 2242 -2248.

[2]王林華，袁明輝，黃慧，等，太赫茲安檢系統人體圖像邊緣物體識別[J]紅外與激光工程，2017，46（ 11）：149-154.

WANG Linhua. YUAN Minghui， HUANG Hui， et al.Recogni-tion of edge object of human bodY in THz security inspectionsystem [J]. Infrared and laser engineering， 2017. 46（11）： 149-154.

[3] ZHANG X C， SHKURINOV A. ZHANG Y.Extreme terahertzscience[J].Nature photonics， 2017. 11（1）： 16-18.

[4] ANJUM A， JAVED S，BHOPAL M F，et al.Terahertz antennaarray as a heat blocking structure [J]. Journal of electromagnet-ic waves and applications， 2015. 29（1）：14-25.

[5]許文忠，鐘凱，梅嘉林，等.太赫茲波在沙塵中的衰減特性[J]紅外與激光工程，2015，44（2）：523-527.

XU Wenzhong， ZHONG Kai， MEI Jialin. et al.THz wave at-tenuation characteristics in sand and dust [J]. Infrared and la-ser engineering， 2015， 44（2）：523-527.

[6] HOSSAIN J，SHARMA S，KISHORE V V N.Multi - peakGaussian fit applicahility to wind speed distribution [J]. Renew-ahle&sustainable energy reviews. 2014. 34（34）： 483-490.

[7]崔珊珊，李琦.基于小波變換的太赫茲數字全息再現像去噪研究[J].紅外與激光T程，2015 ，44（6）：1836-1840.

CUI Shanshan， LI Qi. Denoising research on terahertz digitalholography based on wavelet transform [J]. Infrared and laserengineering， 2015， 44（6）： 1836-1840.

[8] DHILLON S S，VITIELLO M S，LINFIELD E H. et al.The2017 terahertz science and technology roadmap[J].Journal ofphysics D： applied physics， 2017， 50（4）： 1-6.

[9] ZHANG X C， SHKURINOV A. ZHANG Y.Extreme terahertzscience [J]. Nature photonics， 2017. 11（1）： 16-18.

[10] KOWALSKI M， KASTEK M. WALCZAKOWSKI M， et al.Passive imaging of concealed objects in terahertz and long-wavelength infrared [J]. Applied optics， 2015， 54（17）： 3826- 3833.

作者簡介：李建蕊（1984-），女，河北唐山人，碩士，講師，研究方向為太赫茲技術和移動通信。

李九生（1976-），男，廣西桂林人，教授，博士生導師，研究方向為太赫茲技術。

史葉欣（1995-），女，內蒙古人，碩士，研究方向為太赫茲技術。