連續(xù)太赫茲波數字全息相襯成像

史曉玉，王大勇,2*，戎 路,2，趙 潔，王云新,2

連續(xù)太赫茲波數字全息相襯成像

史曉玉1，王大勇1,2*，戎 路1,2，趙 潔1，王云新1,2

1北京工業(yè)大學應用數理學院，北京 100124；2北京工業(yè)大學應用數理學院，北京市精密測控技術與儀器工程技術研究中心，北京 100124

太赫茲波具有獨特的低能性、高穿透性、懼水性等成像特性，將其應用于相襯成像能夠反映物體的內部結構和更加豐富全面的生物信息，在生物醫(yī)學檢測等領域具有重要的應用。其中，太赫茲波數字全息成像是一種可以給出定量的振幅和相位信息的非接觸、全場相襯成像方法，是太赫茲成像技術領域的重要研究方向之一。本文基于連續(xù)太赫茲源，從離軸式和同軸式數字全息成像的相襯成像原理、光路系統(tǒng)和再現算法多個方面，介紹了相關技術的研究現狀，分析了太赫茲源、再現算法等因素對成像分辨率的影響，并對太赫茲數字全息的發(fā)展趨勢進行了展望。

太赫茲成像；數字全息；分辨率；相襯成像

1 引 言

太赫茲波(Terahertz，THz)是指頻率為0.1 THz~10 THz(波長在30 μm~3 mm)范圍內的電磁波[1-2]，位于毫米波與紅外波之間，是電磁波譜中唯一尚未完全開發(fā)的頻段。太赫茲波具有低能性、穿透性、懼水性、指紋吸收譜等重要的特征，在生物醫(yī)學、物理科學、電子信息、國防航天等諸多領域具有廣闊的應用前景，引起了世界各國科研人員的廣泛關注。太赫茲技術被評為“改變未來世界的十大技術之首”，也是國際上重要的前沿交叉領域。其中，太赫茲成像技術被認為是太赫茲波科學與技術中最具有應用前景的研究領域之一[3-4]。

太赫茲波相襯成像不僅可以利用材料吸收獲得空間密度分布，還可以通過相位測量得到折射率的空間分布，能夠反映更準確的內部結構和更加豐富的材料信息。傳統(tǒng)的相襯成像主要是基于太赫茲波時域脈沖成像系統(tǒng)(terahertz time-domain spectroscopy，TDS)實現的，這種成像方法通過相干延遲掃描探測的方式可以同時獲得物體的振幅和相位分布[5-6]。此外，該方法還可以通過光譜成像探測物體的光譜信息。然而，太赫茲波時域成像技術通常需要時間延遲采樣，難以快速全場獲得樣品的幅值和相位信息。為滿足生物醫(yī)學和生命科學研究中的實際需求，需要進一步探索其他有效的相襯成像方法。

可見光波段的數字全息成像是一種較為成熟的相襯成像方法，采用CCD和CMOS等二維面陣探測器替代傳統(tǒng)記錄材料獲取全息圖，并通過計算機實現衍射傳播和數值再現，可以給出定量的振幅和相位信息，具有全場、無機械掃描、實時成像的特點，在生物醫(yī)學成像、無損檢測等領域取得了豐富的研究成果[7-11]。在太赫茲波段引入數字全息成像技術能充分利用后者的優(yōu)勢，從所獲取的太赫茲干涉強度圖像中提取出樣品的定量振幅與相位信息，可為太赫茲波相襯成像提供一種新的技術途徑。該方法既能穿透可見光波段非透明的非極性材料，又能實現全視場定量相襯成像，有望在生物醫(yī)學領域得到廣泛應用。

近年來，太赫茲器件得到了迅速發(fā)展，出現了可靠、穩(wěn)定、高功率的連續(xù)(continuous-wave, CW)太赫茲發(fā)射源，也出現了能在室溫下工作、像元尺寸更小、像素個數更多的面陣式探測器，為太赫茲波與數字全息的結合提供了充分的可能。1) 在連續(xù)太赫茲激光源方面：太赫茲氣體激光器(far-infrared laser，FIRL)是被二氧化碳激光器泵浦的大功率可調諧連續(xù)太赫茲波輻射源，主要應用輻射頻率在2.52 THz (118.8 μm)、輸出功率可以達到500 mW，是在室溫下運行的高功率連續(xù)太赫茲輻射源[12]。量子級聯激光器(quantum cascade laser，QCL)是基于多層半導體結構電子躍遷的電泵浦激光器[13-15]。這類激光器大部分是在低溫下實現的，在低溫下運行的輸出功率可以達到幾瓦，是太赫茲波相干成像的理想光源。耿氏二極管(Gunn diode)振蕩器是一類基于半導體技術的緊湊型連續(xù)太赫茲源。基于砷化鎵和氮化鎵的耿氏二極管輸出頻率最高分別可達到200 GHz和3 THz，在室溫下可實現0.1 THz~1 THz范圍內的單色連續(xù)波輻射[16-17]。倍頻鏈(multiplier chain)通常與耿氏二極管等電子源結合，產生更高頻率的輻射波。倍頻器是基于肖特基二極管或晶體管的電子電路，對電磁波具有非線性響應[17]。入射波通過倍頻器產生高次諧波分量，并利用濾波器或波導進行收集。毫米波矢量網絡分析儀(vector network analyzer，VNA)的發(fā)射和接收部分常通過集成倍頻模塊來進行擴頻?；夭ㄕ袷幤?backward-wave oscillator，BWO)可以通過改變加速電壓產生覆蓋微波到太赫茲波頻率范圍的可調諧激光輸出[18-19]，電子槍附近的輸出功率在1 mW (1 THz)到50 mW (0.2 THz)之間。這類太赫茲源的輸出功率高、譜線窄、頻率可調以及其產生的波前質量好，常被用作相干太赫茲成像的光源[20]。2) 在太赫茲面陣探測器方面：熱釋電探測器是基于LiTaO3熱釋電晶體實現的面陣式探測器，能夠探測到從紅外到太赫茲波段的輻射光。商用的熱釋探測器包括Ophir-Spiricon公司Pyrocam III(像素個數124×124，像素間隔100 μm)和Pyrocam IV(像素個數320×320，像素間隔80 μm)兩種型號[21]。微測熱輻射計是集成在芯片上的測熱輻射計探測器，借助對溫度敏感的電阻探測器來間接測量電磁輻射功率?；谖y熱輻射計熱探測器成像得到的圖像中噪聲主要來源于熱輻射噪聲。目前常用的微測熱輻射計包括Xenics Gobi-640(像素個數640 × 480，像素間距17 μm)[22]，Devitech IR-032(像素個數640 × 480，像素間距25 μm)[23-24]，Miricle 307(像素個數640 × 480，像素間距25 μm)[25]，IRay Technology(像素個數640×512，像素間距17 μm)[26]，NEC IRV-T0831(像素個數320 × 240，像素間距23.5 μm)[27]。

本文主要介紹目前連續(xù)太赫茲波數字全息相襯成像技術的研究進展，從離軸式全息和同軸式全息兩方面展開介紹，包括連續(xù)太赫茲波數字全息成像記錄與再現的機理，實驗裝置設計、數值再現算法和成像結果等。

2 連續(xù)太赫茲波離軸數字全息成像

連續(xù)太赫茲波數字全息成像包括記錄和再現兩個過程。太赫茲波離軸數字全息成像的全息圖記錄過程如圖1所示[28]，太赫茲波通過分束后，一束光照射在物平面(0,0)上產生受樣品調制的太赫茲波，即物光波；另一束光以一定角度傳播到距離為的記錄平面(,)，即參考光波。攜帶物體信息的物光與參考光干涉，干涉光場的強度分布即為全息圖，通過探測器記錄下來，記錄面的太赫茲波全息圖強度分布可表示為

連續(xù)太赫茲波數字全息的再現過程是利用計算機完成的。現階段太赫茲波段的面陣式探測器像素個數較少、靶面尺寸較小，直接記錄的圖像信噪比較低，需要在進行數值再現前對記錄的全息圖進行預處理，提高全息圖中條紋的對比度，抑制噪聲。然后通過頻譜濾波等后期處理把實像和孿生像分離出來，通過基于傅里葉變換運算的衍射傳播再現算法模擬太赫茲波的傳播過程，基于瑞利?索末菲衍射傳播條件將記錄面復振幅回傳至物平面。常用的再現算法主要包括基于菲涅耳衍射理論的菲涅耳變換再現方法，基于瑞利-索末菲衍射積分公式的卷積法和基于標量衍射的亥姆霍茲方程的角譜法[29]。再現過程中可借助自聚焦算法獲得準確的再現距離和物體信息，得到的物體強度分布(0,0)和包裹相位(0,0)分布可分別表示為

對于光學厚度在一個波長以上的樣品，還需要通過相位解包裹算法來獲取實際的相位值。

圖2 (a) 基于耿氏型二極管振蕩器的數字全息成像裝置；(b) 鏤空金屬字母的重建圖像[31]

2.1 透射式連續(xù)太赫茲波離軸數字全息成像

早期的連續(xù)太赫茲波數字全息成像系統(tǒng)是通過單點探測器的逐點掃描來完成數字全息圖。愛爾蘭國立梅努斯大學的Mahon等人在2004年利用100 GHz的耿氏(Gunn)二極管振蕩器，搭建了離軸無透鏡菲涅耳太赫茲數字全息成像系統(tǒng)[30]。太赫茲波照射聚四氟乙烯樣品，利用肖特基二極管逐點掃描來記錄全息圖，光路如圖2(a)所示[31]。重建得到鏤空金屬字母M的振幅圖像如圖2(b)所示，成像分辨率約為9 mm (3)。

美國阿拉巴馬州立大學Heimbeck等[32]基于調諧頻率范圍為0.66 THz~0.76 THz的50 μW相干連續(xù)太赫茲源與肖特基二極管探測器，搭建了馬赫-曾德型離軸太赫茲數字全息裝置(如圖3(a)所示)。實驗中，分別在0.68 THz和0.725 THz頻率下記錄相應的太赫茲數字全息圖，利用雙波長解包裹方法重建得到聚甲基戊烯材料樣品的相襯圖像，實驗結果如圖3(b)所示，其軸向相位精度約為1/40、橫向分辨率約為0.64 mm，解決了2π相位解包裹問題。然而，基于肖特基二極管探測的成像系統(tǒng)需要逐點掃描來記錄全息圖，較難實現快速全場成像。

圖3 (a) 馬赫-曾德型連續(xù)太赫茲波離軸數字全息裝置；(b) 解包裹方法重建得到聚甲基戊烯平凸透鏡的相位圖像[32]

圖4 (a) 基于面陣式熱釋電探測器的太赫茲波離軸成像系統(tǒng)；(b) 間隔為0.4 mm的金屬橫靶條強度再現結果；(c) 間隔為0.4 mm的金屬豎靶條強度再現結果[33]

2011年，哈爾濱工業(yè)大學利用商業(yè)二氧化碳泵浦連續(xù)太赫茲源(2.52 THz)和面陣式熱釋電探測器(像素個數124×124，像素間隔100 μm，像素尺寸85 μm × 85 μm)搭建了太赫茲離軸數字全息成像系統(tǒng)，其橫向分辨率達到0.4 mm (3.4)，如圖4所示[33]。這種太赫茲源和探測器在連續(xù)太赫茲全息成像中得到了廣泛的應用。此后，也可以使用面陣非掃描記錄太赫茲強度圖像，使得太赫茲數字全息成像速率和分辨率都有了較大的提高。隨后，該課題組還驗證了離軸連續(xù)太赫茲波數字全息成像中菲涅耳角譜算法、瑞利–索末菲卷積算法和角譜積分等數字全息再現算法的有效性[34]。考慮到太赫茲波的波長相對較長，研究者在實驗中通常選擇較短的記錄距離來提高分辨率。研究表明，角譜積分法比其他兩種再現算法更適合在這種實驗條件下來重建圖像。該課題組還將記錄距離從39 mm[33]減小到27 mm[35]，相應的成像分辨率進一步提升到0.245 mm (2.1)，高于同一時期太赫茲遠場焦平面成像的分辨率。

隨著太赫茲波段新型光源的發(fā)展，具有小型化、高功率和高頻率特點的太赫茲量子級聯激光器逐漸商業(yè)化，并被用于太赫茲數字全息全場成像系統(tǒng)。從2014年起，瑞士聯邦材料科學與技術實驗室的Hack等人基于太赫茲量子級聯激光器(3 THz)和非制冷微測熱輻射計(像素個數640×480，像素尺寸25 μm×25 μm)搭建了如圖5(a)所示的勞埃德鏡干涉光路，并開展了一系列太赫茲離軸數字全息成像研究。實驗中分別通過模擬和實驗的方式獲得金屬西門子和聚丙烯薄板的再現圖像，結果如圖5(b)所示，其相應的橫向分辨率為0.280 mm (2.8)，相襯成像精度約為0.5 rad[23]。日本德島大學的Yamagiwa等人在2015年基于相同的太赫茲源以及非制冷微測熱輻射計(像素個數320×240，像素尺寸23.5 μm×23.5 μm)開展了太赫茲數字全息成像研究[36]，對塑料和硅片的復雜結構進行成像，系統(tǒng)的縱向精度達到1.7 μm (0.017)。

2.2 反射式連續(xù)太赫茲波離軸數字全息成像

作為另一種離軸全息成像方法，反射式連續(xù)太赫茲波離軸數字全息也得到了廣泛的研究。2005年，俄羅斯新西伯利亞國立大學的Cherkassky等人基于準連續(xù)的高功率太赫茲源(波長范圍為120 μm~180 μm)提出了反射式離軸太赫茲波數字全息成像方法[37]。通過增強型CCD和InAs紅外熱錄像儀記錄純幅值物體的太赫茲全息圖。探測過程中的熱轉換弛豫時間降低了時間分辨率。因此，該系統(tǒng)的曝光時間和動態(tài)范圍嚴重影響了系統(tǒng)的成像空間分辨率和時間分辨率。2008年，芬蘭赫爾辛基理工大學的Tamminen等[38]提出了一種長距離探測的反射式太赫茲波數字全息成像方法。毫米波擴展的矢量網絡分析儀通過兩個波紋喇叭天線生成照明光束和參考光束(310 GHz)，兩束光的強度比為9:1。接收端的矢量網絡分析儀通過開放式波導天線記錄全息圖(400 mm×400 mm)，距離接收器1.5 m處的橫向分辨率理論上可以達到2 mm (~2)。

另一方面，太赫茲波單點探測器可以直接獲取連續(xù)太赫茲波強度信息，然而這種探測器需要對數字全息圖進行逐點掃描，難以直接全場記錄。2010年，哈爾濱工業(yè)大學基于0.1 THz和0.12 THz的連續(xù)太赫茲源實現了雙波長解包裹[39]，通過掃描參考光束來逐點記錄樣品的干涉條紋。2011年，清華大學研究組利用光泵連續(xù)太赫茲源(2.52 THz)和高萊點探測器(Golay cell)設計了如圖6(a)所示的四步相移邁克爾遜像面干涉裝置[40]，實驗結果如圖6(b)所示，獲得了塑料薄片的相襯信息，橫向分辨能力達到0.2 mm，軸向測量范圍達到40 μm。

微測熱輻射計等太赫茲波面陣探測器的興起，為直接記錄全場全息圖在硬件上提供了支持。2015年，瑞士Zolliker等人采用光泵連續(xù)太赫茲源(FIRL100型，2.52 THz)與非制冷氧化釩微測熱輻射計(像素個數640×480，像素尺寸25 μm×25 μm)搭建了離軸太赫茲波數字全息成像系統(tǒng)(如圖7所示)[24]。2015年，意大利國家光學研究所的Locatelli等人搭建了反射式和透射式兩種離軸太赫茲數字全息裝置[25]，利用商業(yè)化的太赫茲量子級聯激光器(2.8 THz)和高分辨率的微測熱輻射計(像素個數640×480，像素尺寸25 μm×25 μm)，對黑色聚丙烯掩模板覆蓋的金屬薄片進行實時成像，得到了再現振幅像。

圖5 (a) 基于太赫茲量子級聯激光器的太赫茲波離軸成像系統(tǒng)；(b) 實驗結果：金屬刻蝕西門子星照片(左)，再現強度像(中)，再現相位像(右)[23]

圖6 (a) 實驗裝置圖；(b) 塑料薄片的三維深度形貌[40]

圖7 基于光泵太赫茲源和微測熱輻射計的太赫茲數字全息成像光路[24]

此外，2013年，中國科學院電子所基于低頻太赫茲源(0.3 THz)和點探測方式也對反射式離軸太赫茲全息成像開展了數值仿真研究，通過菲涅耳衍射積分的卷積算法再現全息圖，分析了束腰與分辨率的關系，理論成像分辨率可以達到2.36 mm[41]。2018年，英國格拉斯哥大學的Humphreys等人利用GPU計算提高反射式太赫茲全息成像的再現速度，以視頻速度(50 Hz)采集樣品信息細節(jié)，分辨能力達到280 μm (2.4)[42]。2019年，本課題組利用亞像素圖像配準[43]和圖像融合算法擴展反射式太赫茲數字全息的視場范圍[44]，對聚四氟乙烯(PTFE)板覆蓋的鍍金書簽進行成像，并與采用三維表面輪廓儀的測量結果進行比較，覆蓋部分的相對測量誤差為5.98%，未覆蓋部分的相對測量誤差為2.93%。

3 連續(xù)太赫茲波同軸數字全息成像

與離軸全息成像結構相比，連續(xù)太赫茲波同軸數字全息成像的參考光與物光共路，因此具有所用元器件較少、光路結構緊湊、對太赫茲源的相干性要求較低、可充分利用探測器有限的空間帶寬積、成像分辨率較高等優(yōu)點[45]。在連續(xù)太赫茲波同軸全息圖的記錄過程中(如圖8)[28]，太赫茲源輻射出的太赫茲波直接照在樣品上，物體的衍射光作為被物體調制的物光波(,)，未衍射的直透太赫茲波作為參考光波(,)，物光波和參考光波共路傳播。由太赫茲探測器采集樣品的衍射光和透射光干涉產生的全息圖。

同軸全息同樣通過數值再現獲得定量的復振幅分布，主要缺點是存在孿生像的串擾。為了再現出物體的相位信息，需要利用相位復原算法重建全息圖的相位分布。相位復原方法的一般流程為迭代過程，大致可分為GS算法、ER算法和混合輸入輸出HIO算法[46]。

連續(xù)太赫茲波同軸數字全息成像方法最早是由哈爾濱工業(yè)大學提出的[47]。2012年，該研究組基于光泵太赫茲激光器(2.52 THz)和熱釋電探測器(像素個數124×124，像素間隔100 μm，像素尺寸85 μm×85 μm像素)對鍍在特氟龍基底的“HIT”圖標進行強度成像，如圖9(a)所示。與太赫茲離軸全息方法相比，同軸數字全息成像系統(tǒng)的記錄距離更短，其再現圖像的分辨率從0.4 mm (3.4)提高到0.2 mm (1.7)，如圖9(b)所示。隨后，他們將太赫茲波同軸數字全息用于遮擋隱蔽目標成像[48]，獲取了多種隱藏物體的幅值再現像，孿生像重疊問題還有待解決，且沒有獲得相位分布。

圖8 連續(xù)太赫茲波同軸數字全息原理示意圖[28]

圖9 (a) 連續(xù)太赫茲波同軸數字全息成像系統(tǒng)；(b) 字母“TIH”的照片；(c) 圖(b)中物體的再現圖像[47]

針對同軸數字全息孿生像混疊問題，本課題組在2014年提出了適用于太赫茲數字全息的相位復原方法，在記錄面與物體平面之間往復迭代去除孿生像[49]，所用相同的太赫茲源(2.52 THz)和熱釋電探測器(像素個數124×124，像素尺寸85 μm×85 μm)搭建了與哈工大類似的同軸數字全息成像系統(tǒng)，實現了蜻蜓后翅標本的連續(xù)太赫茲波數字全息成像?？紤]到熱釋電探測器的靈敏度較低，提出多幅累加和歸一化方法對全息圖進行預處理，抑制了系統(tǒng)太赫茲光斑的不均勻性，提高了全息圖的條紋對比度。2015年，哈爾濱工業(yè)大學研究了再現像的均方差與記錄距離、多幅全息圖之間位置間隔的關系，結果表明只有當記錄距離小于或等于2 mm時，記錄距離的選取才會對再現像產生影響[50]。

為了改進太赫茲波同軸數字全息成像的再現質量和分辨率，本課題組提出了同軸全息圖外推迭代算法，擴大了探測器的數值孔徑，獲得了高分辨率的翅脈結構信息(如圖10(a)所示)，揭示了連續(xù)太赫茲波同軸數字全息技術在生物醫(yī)學成像領域應用的可行性。隨后，又將外推法用于人肝癌組織切片成像[51]，采用外推法和亞像素微位移法，將理論分辨率從240 μm (~2)提升為158 μm (1.3)，從再現的相位再現像中可以觀察到人體原發(fā)性肝癌組織切片中的高密度纖維化組織(如圖10(b)所示)。由于肝纖維化是肝癌和其他疾病的早期征象，因此連續(xù)太赫茲波數字全息技術有望應用于癌癥診斷醫(yī)學成像領域。

圖10 (a) 高分辨率的翅脈再現相位圖；(b) 人體肝癌切片的再現結果[49-51]

后續(xù)的研究主要聚焦在改進物平面以及記錄面的約束條件，從而提高太赫茲數字全息成像的分辨率以及成像質量。哈爾濱工業(yè)大學研究了物平面的幅值、相位和支持域約束條件及其對迭代算法收斂性的影響，此外還分析了補零擴展、邊界復制擴展和切趾操作對再現圖像的影響[52]。2016年，該研究組在物平面上使用支持域約束條件對太赫茲同軸數字全息圖進行迭代再現[53]。隨后，他們采用雙曝光法進行了仿真和實驗研究[54]。與單次曝光的振幅約束相位復原算法相比，雙曝光的相位約束相位復原算法重建的圖像對比度增加了0.146，雙曝光的振幅約束相位復原算法重建的圖像對比度增加了0.225。2019年，中國工程物理研究院基于自制的光泵5.24 THz激光器以及高分辨微測熱輻射計(像素個數640×512，像素尺寸17 μm×17 μm)構建了片上連續(xù)太赫茲波同軸數字全息系統(tǒng)[26]，記錄距離為3 mm。采用物面約束條件和L1稀疏約束條件，在連續(xù)太赫茲波數字全息成像中實現了亞波長分辨率(40 μm，0.7)，如圖11所示。

4 提高連續(xù)太赫茲波數字全息成像性能及其他研究

單幅全息圖的外推改進再現圖像分辨率的程度是有限的，為此，研究者們提出通過移動探測器合成孔徑的方法來提高分辨率，該方法在太赫茲離軸數字全息[24]和太赫茲同軸數字全息中[55-56]都得到了驗證。2015年，瑞士Zolliker等[24]提出了合成孔徑太赫茲數字全息成像方法。為了擴大參考光的照明范圍，在實驗中移動參考光的反射鏡來完成相移成像，并通過菱形路徑掃描移動平移臺來移動微測熱輻射計記錄全息圖，拼接所有子全息圖完成合成孔徑成像。與基于單幅全息圖的再現結果相比，基于合成孔徑方法大幅提升了系統(tǒng)成像分辨率(如圖12)，系統(tǒng)的橫向分辨率達到200 μm (1.68)，相位精度達到0.4 rad的，對應的軸向分辨率為6 μm。2016年，北京工業(yè)大學與中國工程物理研究院合作，基于單色量子級聯激光器搭建了太赫茲波同軸數字全息成像系統(tǒng)[55]，激光器的波長在97 μm，97.6 μm和98.9 μm三處對應功率之比為4:1:2。利用微測熱輻射計(像素個數320×240，像素尺寸23.5 μm×23.5 μm)將9幅子全息圖(7.52 mm×5.64 mm)圖像合成為一幅全息圖(15.5 mm×13.6 mm)，獲得了125 μm (1.28)的分辨率，成功檢測到微結構硅臂中的金屬污染物。2017年，中國工程物理研究院研發(fā)出4.3 THz的量子級聯激光器，并將其運用于連續(xù)太赫茲波同軸全息成像系統(tǒng)中，利用合成孔徑方法將分辨率提高到70 μm (1)[27]。

圖11 (a), (d) 分辨率為50 μm和40 μm物體的光學顯微圖像；(b), (c) 分辨率為50 μm物體的再現復振幅分布；(e), (f) 分辨率為40 μm物體的再現復振幅分布[26]

圖12 西門子星樣品的(a)強度再現圖像和(b)相位再現圖像[26]

為了進一步提高連續(xù)太赫茲波同軸數字全息的成像質量，研究者對太赫茲波數字全息圖像的去噪方法開展了研究。哈爾濱工業(yè)大學提出通過馬爾可夫鏈蒙特卡爾采樣理論對太赫茲數字全息成像圖像進行去噪處理[57]。2015年，他們基于小波變換對太赫茲同軸數字全息再現像的去噪進行了研究。結果證明采用“bior2.2”小波基的同態(tài)濾波對太赫茲同軸數字全息圖的去噪效果較好[58]。

另一方面，在連續(xù)太赫茲波同軸數字全息成像中，再現距離是一個重要的參數。再現距離存在誤差將無法準確得到物平面上的復振幅分布，嚴重的離焦情況還會導致觀察不到樣品的真實形貌。針對這一問題，本課題組研究了連續(xù)太赫茲波同軸數字全息自動聚焦算法[59]。2017年，該課題組基于離軸式和同軸式兩種光路，提出了一種基于自聚焦算法的多平面太赫茲波數字全息成像方法[60]。利用自動聚焦判據曲線處理圖像，得到了樣品不同位置處(15 mm和41 mm)的聚焦復振幅分布，并合成為三維空間分布，有效地擴展了太赫茲波同軸全息成像的景深。2016年，中國工程物理研究院采用雙平面記錄法，在不同的探測距離處記錄了兩組同軸全息圖[61]。他們利用亞像素微位移法對全息圖進行精細采樣。記錄平面的數據冗余度越大，迭代收斂速度越快，再現圖像的質量越好，再現圖像的分辨率達到150 μm (1.3)。

為了解決2π相位解包裹問題，埃及的Ibrahim等人基于3 THz的量子級聯激光器和非制冷微測熱輻射計搭建了反射式太赫茲數字全息光路，采用直接傅里葉以及埃爾米特多項式算法解決了這一問題，并測量了楔形(從80 μm到110 μm)聚合物薄膜的厚度[62]。

除此之外，基于脈沖太赫茲源也可以實現數字全息成像，首都師范大學于2008年基于脈沖式太赫茲源搭建了數字全息成像系統(tǒng)[63]。該系統(tǒng)中飛秒激光器(800 nm, 100 fs)產生的脈沖光被分為泵浦光和探測光。泵浦光擴束后照在碲化鋅晶體上，通過整流效應產生的太赫茲脈沖光束照射樣品，攜帶樣品波前信息的太赫茲透射光波照在另一片碲化鋅上，相對應的探測光通過偏振調節(jié)后也照射到晶體上。太赫茲電場的波前復振幅信息通過偏振態(tài)的改變加載到探測光上。攜帶樣品信息的探測光被晶體的左側表面反射并投影到CCD上。通過連續(xù)調節(jié)太赫茲光與探測光的光程差，測量不同時間點的太赫茲圖像，對太赫茲時域信號進行傅里葉變換，獲得不同頻率成分的太赫茲幅值和相位信息。由于脈沖太赫茲成像可直接獲取樣品的復振幅信息，因此全息圖的記錄過程中無需使用參考光。此外，由于太赫茲脈沖是可見光通過晶體的相互作用激發(fā)的，該方法還可以利用可見光的光場調控方法間接調控太赫茲脈沖的波前信息，從而豐富太赫茲數字全息成像的應用。國內首都師范大學[64-70]和俄羅斯圣光機大學的Petrov等[71-73]陸續(xù)在脈沖太赫茲波數字全息成像及其應用方面做了大量的研究工作，取得了重要的進展。

5 結 論

通過世界各個研究組的不斷探索研究，連續(xù)太赫茲波數字全息成像技術快速發(fā)展，取得了顯著的成果和關鍵性的進展。本文總結了離軸透射式和反射式、以及同軸式連續(xù)太赫茲波數字全息成像的研究進展，包括數字全息成像的成像原理、光路系統(tǒng)和再現算法。對于離軸數字全息，借助新型的面陣式太赫茲波探測器，全息圖記錄由逐點掃描發(fā)展為全場成像，提高了成像速率，并對提高分辨率、擴大視場范圍、遮擋下的形貌檢測開展研究。在連續(xù)太赫茲波同軸數字全息成像研究中，主要去除孿生像的相位復原算法計算效率和再現精度開展研究。此外，對于太赫茲波數字全息成像分辨率的提高、成像質量的改進以及成像速度的提升也積累了一定的研究成果。這些研究成果為進一步推動太赫茲成像在生物醫(yī)學、智能制造、無損檢測等領域的應用提供了技術儲備。

連續(xù)太赫茲波數字全息技術的發(fā)展還面臨著一系列科學技術問題有待解決：1) 針對太赫茲波的傳輸特性優(yōu)化衍射傳播和再現算法，深入研究太赫茲波在物體邊界和內部的傳輸機理，以及與材料的相互作用機理。2) 提高太赫茲數字全息圖的條紋對比度和系統(tǒng)穩(wěn)定性，目前的太赫茲光泵激光器和量子級聯激光器等連續(xù)太赫茲源存在光斑質量較差和功率不穩(wěn)定性等問題，以熱釋電探測器為代表面陣探測器的靈敏度、響應速率和動態(tài)范圍較差，迫切需要開發(fā)和提升太赫茲波段元器件性能，提高連續(xù)太赫茲數字全息的成像質量。3) 現有的太赫茲數字全息成像分辨率還不足以完全滿足實際成像應用的需求，一方面需要發(fā)展大數值孔徑的太赫茲波段顯微物鏡，另一方面可以將在其它波段已經充分驗證的遠場無透鏡超分辨率成像方法，如表面等離激元、照明光場調控等引入到太赫茲波段，提高太赫茲數字全息成像分辨率。4) 現有的連續(xù)太赫茲波數字全息只能獲得物光波的定量復振幅分布，需要解耦算法才能得到折射率分布，同時可以與衍射層析等三維成像方法相結合，獲得樣品內部吸收系數和折射率三維分布，獲得物體內部更豐富的信息。5) 連續(xù)太赫茲數字全息還處于應用探索階段，需要進一步建立光學參數與樣品內部特性的聯系，開發(fā)其在生物醫(yī)學、材料表征、無損檢測等領域的應用。

[1] Zhang X C, Shkurinov A, Zhang Y. Extreme terahertz science[J]., 2017, 11(1): 16–18.

[2] Tonouchi M. Cutting-edge terahertz technology[J]., 2007, 1(2): 97–105.

[3] De Cumis U S, Xu J H, Masini L,. Terahertz confocal microscopy with a quantum cascade laser source[J]., 2012, 20(20): 21924–21931.

[4] Suga M, Sasaki Y, Sasahara T,. THz phase-contrast computed tomography based on Mach-Zehnder interferometer using continuous wave source: proof of the concept[J].2013, 21(21): 25389–25402.

[5] Fischer B M, Hoffmann M, Helm H,. Terahertz time-domain spectroscopy and imaging of artificial RNA[J]., 2005, 13(14): 5205–5215.

[6] Zhong H, Redo-Sanchez A, Zhang X C. Identification and classification of chemicals using terahertz reflective spectroscopic focal-plane imaging system[J]., 2006, 14(20): 9130–9141.

[7] Bianco V, Memmolo P, Leo M,. Strategies for reducing speckle noise in digital holography[J]., 2018, 7(1): 48.

[8] Nelson J W, Knefelkamp G R, Brolo A G,. Digital plasmonic holography[J]., 2018, 7(1): 52.

[9] Tikan A, Bielawski S, Szwaj C,. Single-shot measurement of phase and amplitude by using a heterodyne time-lens system and ultrafast digital time-holography[J].2018, 12(4): 228–234.

[10] Schnars U, Falldorf C, Watson J,.[M]. New York: Springer, 2015.

[11] Asundi A.[M]. Chichester/Hoboken: Wiley, 2011.

[12] Minamide H, Ito H. Frequency-agile terahertz-wave generation and detection using a nonlinear optical conversion, and their applications for imaging[J]., 2010, 11(7–8): 457–471.

[13] Faist J, Capasso F, Sivco D L,. Quantum cascade laser[J]., 1994, 264(5158): 553–556.

[14] Rochat M, Ajili L, Willenberg H,. Low-threshold terahertz quantum-cascade lasers[J].2002, 81(8): 1381–1383.

[15] Kumar S. Recent progress in terahertz quantum cascade lasers[J]., 2011, 17(1): 38–47.

[16] Gribnikov Z S, Bashirov R R, Mitin V V. Negative effective mass mechanism of negative differential drift velocity and terahertz generation[J]., 2001, 7(4): 630–640.

[17] Maestrini A, Thomas B, Wang H,. Schottky diode-based terahertz frequency multipliers and mixers[J]., 2010, 11(7–8): 480–495.

[18] Dobroiu A, Yamashita M, Ohshima Y N,. Terahertz imaging system based on a backward-wave oscillator[J]., 2004, 43(30): 5637–5646.

[19] Martens S, Gompf B, Dressel M. Characterization of continuous-wave terahertz sources: laser mixing versus backward-wave oscillators[J].2009, 48(29): 5490–5496.

[20] Baik C W, Son Y M, Kim S I,. Microfabricated coupled-cavity backward-wave oscillator for terahertz imaging[C]//, 2008: 398–399.

[21] Wan M, Muniraj I, Malallah R,. Sparsity based terahertz reflective off-axis digital holography[J]., 2017, 10233: 102330T.

[22] Valzania L, Feurer T, Zolliker P,. Terahertz ptychography[J]., 2018, 43(3): 543–546.

[23] Hack E, Zolliker P. Terahertz holography for imaging amplitude and phase objects[J]., 2014, 22(13): 16079–16086.

[24] Zolliker P, Hack E. THz holography in reflection using a high resolution microbolometer array[J]., 2015, 23(9): 10957–10967.

[25] Locatelli M, Ravaro M, Bartalini S,. Real-time terahertz digital holography with a quantum cascade laser[J].2015, 5: 13566.

[26] Li Z Y, Yan Q, Qin Y,. Sparsity-based continuous wave terahertz lens-free on-chip holography with sub-wavelength resolution[J]., 2019, 27(2): 702–713.

[27] Deng Q H, Li W H, Wang X M,. High-resolution terahertz inline digital holography based on quantum cascade laser[J].2017, 56(11): 113102.

[28] Huang H C. The research on continuous-wave terahertz in-line digital holographic imaging method[D]. Beijing: Beijing University of Technology, 2017.

黃昊翀. 連續(xù)太赫茲波同軸數字全息成像方法的研究[D]. 北京: 北京工業(yè)大學, 2017.

[29] Goodman J W.[M]. 3rd ed. Greenwoood Village: Roberts & Company Publishers, 2005.

[30] Mahon R, Murphy A, Lanigan W. Terahertz holographic image reconstruction and analysis[C]//, 2004: 749–750.

[31] Mahon R J, Murphy J A, Lanigan W. Digital holography at millimetre wavelengths[J].2006, 260(2): 469–473.

[32] Heimbeck M S, Kim M K, Gregory D A,. Terahertz digital holography using angular spectrum and dual wavelength reconstruction methods[J]., 2011, 19(10): 9192–9200.

[33] Ding S H, Li Q, Li Y D,. Continuous-wave terahertz digital holography by use of a pyroelectric array camera[J]., 2011, 36(11): 1993–1995.

[34] Li Q, Li Y D, Ding S H,. Terahertz computed tomography using a continuous-wave gas laser[J].2012, 33(5): 548–558.

[35] Li Q, Ding S H, Li Y D,. Experimental research on resolution improvement in CW THz digital holography[J]., 2012, 107: 103–110.

[36] Yamagiwa M, Ogawa T, Minamikawa T,. Real-time amplitude and phase imaging of optically opaque objects by combining full-field off-axis terahertz digital holography with angular spectrum reconstruction[J]., 2018, 39(6): 561–572.

[37] Cherkassky V S, Knyazev B A, Kubarev V V,. Imaging techniques for a high-power THz free electron laser[J]., 2005, 543(1): 102–109.

[38] Tamminen A, Ala-Laurinaho J, Raisanen A V. Indirect holographic imaging at 310 GHz[C]//, 2008: 168–171.

[39] Wang X K, Hou L, Zhang Y. Continuous-wave terahertz interferometry with multiwavelength phase unwrapping[J].2010, 49(27): 5095–5102.

[40] Wang Y X, Zhao Z R, Chen Z Q,. Continuous-wave terahertz phase imaging using a far-infrared laser interferometer[J].2011, 50(35): 6452–6460.

[41] Gao X, Li C, Fang G Y. Study of image reconstruction for terahertz indirect holography with quasi-optics receiver[J].2013, 30(6): 1291–1296.

[42] Humphreys M, Grant J P, Escorcia-Carranza I,. Video-rate terahertz digital holographic imaging system[J]., 2018, 26(20): 25805–25813.

[43] Guizar-Sicairos M, Thurman S T, Fienup J R. Efficient subpixel image registration algorithms[J].2008, 33(2): 156–158.

[44] Wang D Y, Zhao Y L, Rong L,. Expanding the field-of-view and profile measurement of covered objects in continuous-wave terahertz reflective digital holography[J]., 2019, 58(2): 023111.

[45] Wang D Y, Huang H C, Zhou X,. Phase-contrast imaging by the continuous-wave terahertz line digital holography[J]., 2014, 41(8): 232–237.

王大勇, 黃昊翀, 周遜, 等. 連續(xù)太赫茲波同軸數字全息相襯成像[J]. 中國激光, 2014, 41(8): 232–237.

[46] Fienup J R. Reconstruction of an object from the modulus of its Fourier transform[J]., 1978, 3(1): 27–29.

[47] Xue K, Li Q, Li Y D,. Continuous-wave terahertz in-line digital holography[J].2012, 37(15): 3228–3230.

[48] Li Q, Xue K, Li Y D,. Experimental research on terahertz Gabor inline digital holography of concealed objects[J].2012, 51(29): 7052–7058.

[49] Rong L, Latychevskaia T, Wang D Y,. Terahertz in-line digital holography of dragonfly hindwing: amplitude and phase reconstruction at enhanced resolution by extrapolation[J]., 2014, 22(14): 17236–17245.

[50] Hu J Q, Li Q, Yang Y F. Simulation research on continuous Terahertz inline digital holography imaging based on phase retrieval algorithm[J]., 2015, 52(1): 92–98.

胡佳琦, 李琦, 楊永發(fā). 基于相位恢復法的連續(xù)太赫茲同軸數字全息成像仿真研究[J]. 激光與光電子學進展, 2015, 52(1): 92–98.

[51] Rong L, Latychevskaia T, Chen C H,. Terahertz in-line digital holography of human hepatocellular carcinoma tissue[J].2015, 5: 8445.

[52] Hu J Q, Li Q, Cui S S. Research on object-plane constraints and hologram expansion in phase retrieval algorithms for continuous-wave terahertz inline digital holography reconstruction[J]., 2014, 53(30): 7112–7119.

[53] Hu J Q, Li Q, Zhou Y. Support-domain constrained phase retrieval algorithms in terahertz in-line digital holography reconstruction of a nonisolated amplitude object[J].2016, 55(2): 379–386.

[54] Hu J Q, Li Q, Chen G H. Reconstruction of double-exposed terahertz hologram of non-isolated object[J]., 2016, 37(4): 328–339.

[55] Huang H C, Rong L, Wang D Y,. Synthetic aperture in terahertz in-line digital holography for resolution enhancement[J].2016, 55(3): A43–A48.

[56] Wan M, Li W H, Wang D Y,. Digital holography of continuous terahertz wave synthetic aperture[J]., 2017, 15(3): 358–363.

萬敏, 黎維華, 王大勇, 等. 連續(xù)太赫茲波合成孔徑數字全息成像方法[J]. 太赫茲科學與電子信息學報, 2017, 15(3): 358–363.

[57] Chen G H, Li Q. Markov chain Monte Carlo sampling based terahertz holography image denoising[J].2015, 54(14): 4345–4351.

[58] Cui S S, Li Q. De-noising research on terahertz digital holography based on wavelet transform[J]., 2015, 44(6): 1836–1840.

崔珊珊, 李琦. 基于小波變換的太赫茲數字全息再現像去噪研究[J]. 紅外與激光工程, 2015, 44(6): 1836–1840.

[59] Huang H C, Wang D Y, Rong L,. Application of autofocusing methods in continuous-wave terahertz in-line digital holography[J].2015, 346: 93–98.

[60] Huang H C, Wang D Y, Li W H,. Continuous-wave terahertz multi-plane in-line digital holography[J].2017, 94: 76–81.

[61] Li Z Y, Li L, Qin Y,. Resolution and quality enhancement in terahertz in-line holography by sub-pixel sampling with double-distance reconstruction[J]., 2016, 24(18): 21134–21146.

[62] Ibrahim D G A. Steep large film thickness measurement with off-axis terahertz digital holography reconstructed by a direct Fourier and Hermite polynomial[J].2018, 57(10): 2533–2538.

[63] Huang H C, Wang D Y, Rong L,. Continuous-wave off-axis and in-line terahertz digital holography with phase unwrapping and phase autofocusing[J].2018, 426: 612–622.

[64] Zhang Y, Zhou W, Wang X,. Terahertz digital holography[J].2008, 44(5): 380–385.

[65] Wang X K, Cui Y, Hu D,. Terahertz quasi-near-field real-time imaging[J].2009, 282(24): 4683–4687.

[66] Wang X K, Cui Y, Sun W F,. Terahertz real-time imaging with balanced electro-optic detection[J].2010, 283(23): 4626–4632.

[67] Wang X K, Cui Y, Sun W F,. Terahertz polarization real-time imaging based on balanced electro-optic detection[J].2010, 27(11): 2387–2393.

[68] Guo L H, Wand X K, Zhang Y. Terahertz digital holographic imaging of biological tissues[J]., 2017, 25(3): 611–615.

郭力菡, 王新柯, 張巖. 生物組織的太赫茲數字全息成像[J]. 光學精密工程, 2017, 25(3): 611–615.

[69] Zheng X H, Wang X K, Sun W F,. Developments and applications of the Terahertz digital holography[J]., 2014, 41(2): 24–34.

鄭顯華, 王新柯, 孫文峰, 等. 太赫茲數字全息術的研發(fā)與應用[J]. 中國激光, 2014, 41(2): 24–34.

[70] Shi J, Wang X K, Zheng X H,. Recent advances in terahertz digital holography[J]., 2017, 10(1): 131–147.

石敬, 王新柯, 鄭顯華, 等. 太赫茲數字全息術的研究進展[J]. 中國光學, 2017, 10(1): 131–147.

[71] Petrov N V, Kulya M S, Tsypkin A N,. Application of terahertz pulse time-domain holography for phase imaging[J]., 2016, 6(3): 464–472.

[72] Balbekin N S, Kulya M S, Petrov N V. Terahertz pulse time-domain holography in dispersive media[J]., 2017, 41(3): 348–355.

[73] Kulya M, Petrov N V, Tsypkin A,. Hyperspectral data denoising for terahertz pulse time-domain holography[J]., 2019, 27(13): 18456–18476.

Phase contrast imaging based on continuous-wave terahertz digital holography

Shi Xiaoyu1, Wang Dayong1,2*, Rong Lu1,2, Zhao Jie1, Wang Yunxin1,2

1College of Applied Science, Beijing University of Technology, Beijing 100124, China;2Institute of Information Photonics Technology, Beijing University of Technology, Beijing 100124, China

圖1 The experimental set-up of the off-axis terahertz digital holography

Overview:Terahertz (THz) radiation is characterized with low-energy, high-penetration, and water-absorption, which could provide internal structure of objects and comprehensive biological information by THz phase contrast imaging. Due to this unique feature, THz radiation has been applied in biomedical imaging, non-destructive testing, and other fields. As an important part of THz imaging technology, continuous-wave THz digital holography (TDH) by recording the complex amplitude in the hologram and numerically retrieving the corresponding phase-shift properties of object, is qualified as a non-invasive and whole-field phase contrast imaging method. With the development of continuous-wave THz sources, detectors, and imaging components, the continuous-wave TDH imaging technology is well developed.

The development and status of off-axis and in-line TDH are reviewed, including the recording and reconstruction theory, experimental setup, and reconstruction algorithms. For the off-axis TDH, the reflective and transmitted TDH has all been introduced. The firstly off-axis TDH configuration is attempted using a 100 GHz Gunn diode oscillator and Schottky-barrier diode. Complete “full-field” phase imaging with higher lateral resolution is achieved using an optically pumped FIR laser and a pyroelectrial array detector, which is then widely applied in continuous-wave TDH configuration later. And the effectiveness of Rayleigh-Sommerfeld convolution algorithm, Fresnel angular spectrum algorithm, and angular spectrum integral are evaluated for off-axis TDH. Miniaturized THz quantum cascade laser with high power and high frequency are also used in TDH full-field imaging system to improve the imaging resolution. For the in-line TDH, the scattered beam by sample interferes with the unscattered part to form the in-line hologram, which is quite suitable for isolated objects imaging. Compared with off-axis TDH, the recording geometry of in-line TDH is more compact, and reconstruction resolution is higher. The twin image is one of the most important problems for in-line TDH, which is solved by iterations with proper constraints in these two planes and phase retrieval algorithm. Constraints on the object plane and extrapolation algorithms for iterative reconstruction of in-line TDH are studied to achieve higher resolution. In addition, different methods are presented to improve the imaging quality of continuous-wave TDH, such as synthetic aperture, denoising method, auto-focusing algorithms, and multi-plane imaging.

In conclusion, the paper summarized the progress of continuous-wave THz digital holography, including the influence of existing THz imaging components and the reconstruction algorithms on resolution and fidelity of imaging are analyzed. And the future development trend of continuous-wave TDH is discussed in the end of the paper.

Citation: Shi X Y, Wang D Y, Rong L,. Phase contrast imaging based on continuous-wave terahertz digital holography[J]., 2020,47(5): 190543

Phase contrast imaging based on continuous-wave terahertz digital holography

Shi Xiaoyu1, Wang Dayong1,2*, Rong Lu1,2, Zhao Jie1, Wang Yunxin1,2

1College of Applied Science, Beijing University of Technology, Beijing 100124, China;2Institute of Information Photonics Technology, Beijing University of Technology, Beijing 100124, China

Terahertz (THz) radiation, due to its unique propagation characters of low-energy, high-penetration, water-absorption, provides internal structure of objects and comprehensive biological information in phase contrast imaging. It has been applied in biomedical imaging, non-destructive testing, and other fields. As an important part of THz imaging technology, continuous-wave (CW) THz digital holography (TDH) is qualified as a non-invasive and whole-field phase contrast imaging method. In this paper, we review the development and status of off-axis and in-line TDH, including the recording and reconstruction theory, experimental setup, and reconstruction algorithms. The influence of existing THz sources and the reconstruction algorithms on resolution and fidelity of imaging are analyzed. And the development trend of TDH is prospected in the end.

terahertz imaging; digital holography; resolution; phase contrast imaging

National Natural Science Foundation of China (61675010) and Beijing Nova Program (2018072)

* E-mail: wdyong@bjut.edu.cn

O438.1；TB877.1

A

史曉玉，王大勇，戎路，等. 連續(xù)太赫茲波數字全息相襯成像[J]. 光電工程，2020，47(5): 190543

10.12086/oee.2020.190543

: Shi X Y, Wang D Y, Rong L,Phase contrast imaging based on continuous-wave terahertz digital holography[J]., 2020, 47(5): 190543

2019-09-12；

2019-12-24

國家自然科學基金面上項目(61675010)；北京市科技新星項目(2018072)

史曉玉(1989-)，女，博士研究生，助理研究員，主要從事隨機激光，微納光學和太赫茲成像的研究。E-mail：xyshi@bjut.edu.cn

王大勇(1968-)，男，博士研究生，教授，主要從事實時光電混合信息處理、光學全息、數字全息、太赫茲成像、微波光子學及空間光調制器等的研究。E-mail：wdyong@bjut.edu.cn

版權所有?2020中國科學院光電技術研究所