全文HTML
--> --> -->根据样品与太赫兹波相互作用方式的不同, 太赫兹成像模式有三种: 透射成像、反射成像和衰减全反射成像. 透射成像是发展最早的成像模式, 通过太赫兹波与样品相互作用后的透射信号成像. 对于在检测波长具有高吸收特性的样品, 透射成像需要对样品进行处理(如: 生物组织需要进行切片处理), 这使得样品制作复杂, 也限制了透射成像的应用范围. 与透射成像相比, 反射成像和衰减全反射成像不需要复杂的样品制作过程. 反射成像是发展较早和应用最广泛的成像模式[9,10]. 反射成像是通过太赫兹波与样品相互作用后的反射信号进行成像, 获得的样品信号通常较弱. 因此, 反射成像需要输出功率较高的太赫兹辐射源[11,12]. 基于样品不同深度位置的折射率和吸收系数不同, 反射成像还可实现样品一定深度的太赫兹成像[13,14]. 对于获取样品反射率, 反射成像需要检测对太赫兹波有高反射特性的样品(如金镜)获得参考信号[15], 这使得实验步骤复杂. 反射成像受样品表面不均匀性引起的漫反射现象影响严重, 使得反射信号较弱[16]. 另外, 连续太赫兹反射成像受反射窗口厚度引起的干涉现象影响严重, 导致获得样品的准确反射率较困难. 衰减全反射成像是基于倏逝波与样品的相互作用而成像, 其可通过检测有无样品的信号比值获得样品反射率[11]. 衰减全反射成像通常采用具有高折射率材料的全反射棱镜进行成像, 故其有望实现高分辨率成像[17]. 衰减全反射成像的穿透深度较低, 仅适用于样品表面检测. 因穿透深度低, 衰减全反射成像对太赫兹波的衰减较小, 故对辐射源的输出功率要求较低[18]. 由以上分析可知, 反射和衰减全反射成像模式具有不同的优缺点. 为了结合这两种成像模式的优点, 研究共光路反射和衰减全反射成像是有必要的.
连续太赫兹反射成像(TR)发展较早. 2003年, 张希成等[19]为了检测哥伦比亚号航天飞机失事的原因, 搭建了第一套小型化、便携式连续太赫兹反射成像系统, 并证明了太赫兹成像技术在无损检测方面的优势. 此后, 连续太赫兹反射成像得到了广泛关注并进行深入研究. 为了实现连续太赫兹反射快速成像, Lee等[20]采用太赫兹气体激光器和焦平面阵列探测器, 实现了连续太赫兹实时反射成像. 然而, 此系统成像分辨率较低, 无法实现样品边界的有效识别. 进一步, Watts等[21]将压缩感知成像与近场成像相结合, 实现分辨率达9 μm的高速太赫兹反射成像. 尽管近场成像可实现高分辨率, 然而压缩感知的成像能力因受算法限制, 将其应用于所有样品较困难. 与之相比, 基于点扫描的连续太赫兹反射成像更易于样品的高精确度检测. 为了实现连续太赫兹反射高对比度成像, Doradla等[22]研究了连续太赫兹反射偏振成像, 实验证明采用正交偏振有助于提高太赫兹反射成像的图像对比度. 然而, 此方法会造成极大的能量损失, 导致成像信噪比较低. 为了提高连续太赫兹反射成像分辨率, Chernomyrdin等[23]基于固体浸没透镜技术, 将连续太赫兹反射成像的空间分辨率提高到了0.15λ—0.3λ之间. 与连续太赫兹反射成像相比, 连续太赫兹衰减全反射成像(ATR)发展较晚. 2016年, 本课题组开发了基于垂直点扫描的连续太赫兹衰减全反射成像系统[11]. 为了避免二次反射对衰减全反射成像的影响, 我们又研究了基于水平点扫描的连续太赫兹衰减全反射成像系统[24]. 目前, 关于反射成像和衰减全反射成像的研究均为相互独立, 相互独立的成像系统具有系统复杂、操作繁琐和不易于样品快速检测的局限性. 为了更好地比较样品在反射和衰减全反射模式下的成像效果, 研究共光路连续太赫兹反射和衰减全反射双模式成像系统更有助于实际应用.
本文理论研究了成像角度对反射和衰减全反射成像的影响, 获得适用于连续太赫兹反射和衰减全反射成像的太赫兹波成像角度. 通过反射窗口和全反射棱镜快速切换的方式, 搭建了基于太赫兹光泵气体激光器为辐射源的共光路连续太赫兹双模式成像系统. 以蒸馏水和猪肉组织为样品, 实验分析了双模式成像的成像分辨率、获得反射信号稳定性、图像对比度和获得样品反射率的准确性. 结果表明, 衰减全反射模式更有助于实现样品的高灵敏度成像.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2021/11/PIC/11-20210182-1_mini.jpg)
Figure1. Schematic diagram of terahertz wave propagation in the imaging window (take reflection as an example).
为了使得更多的太赫兹波能量进入反射窗口, 反射窗口前表面需要对太赫兹波有高透射特性. 为了获得高反射信号, 反射窗口后表面需对太赫兹波有高反射特性. 太赫兹波在反射窗口界面的反射系数、透射系数, 可由菲涅耳公式获得, 其公式为[25]:
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2021/11/PIC/11-20210182-2_mini.jpg)
Figure2. (a) The relation between the reflection coefficient and transmission coefficient of S and P polarized waves at the interface and the incident angle θ1; (b) the relation between the reflection coefficient r and incident angle θ1 of terahertz wave after the front surface transmission, back surface reflection and front surface transmission through reflection window.
衰减全反射成像中, 太赫兹波在成像窗口与样品界面的入射角θ3需大于临界角θC, 实验中通常采用高于样品折射率的材料作为全反射成像窗口, 如高阻硅材料(silicon, n = 3.42 @ 2.52 THz). 衰减全反射模式是基于倏逝波的衰减而成像, 倏逝波的振幅随穿透深度的增加急骤衰减. 穿透深度x0公式如下[17]:
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2021/11/PIC/11-20210182-3_mini.jpg)
Figure3. The relationship between the penetration depth and the incident angle of evanescent wave in different samples, when using the high resistance silicon total reflection prism.
由图3可知, 当θ3 ≥ 43°时穿透深度值变化较小. 因此, 基于此角度进行了全反射成像窗口设计. 通常生物组织的折射率小于水的折射率(n = 2.05), 水为待测样品时的临界角为34.8°. 为了实现共光路反射和衰减全反射成像研究, 考虑到适用于反射成像的入射角θ1, 设计了底角为49°的等腰三角形状的全反射棱镜. 此时太赫兹波在成像面的入射角大于临界角(43.5° > 34.8°), 适用于生物组织的衰减全反射成像研究. 进一步, P偏振下连续太赫兹衰减全反射成像的样品反射率RP的公式可表示为[27]:
3.1.实验装置
图4为基于点扫描的共光路连续太赫兹衰减全反射和反射成像系统示意图, 通过将全反射棱镜替换为反射窗口即可实现双模式成像. 此研究中, 采用CO2激光泵浦的连续THz波气体激光器(FIRL100, Edinburgh Instruments Ltd, UK)作为THz源, 其具有输出功率高、连续波运转、可调谐、光束质量好及稳定性好的特点. 在频率为2.52 THz时, 输出功率最高可达150 mW. 因此, 我们将2.52 THz作为工作频率. 采用可在室温工作且稳定性好的高莱探测器(GC-1P, Tydex Ltd.)接收太赫兹波. 考虑到高莱探测器的慢响应特性, 采用频率为50 Hz的斩波器对连续太赫兹波进行调制. 为了降低激光器的功率波动对成像结果的影响, 采用太赫兹波线栅偏振器(Microtech Instruments, Inc.)将太赫兹波分成两束: 参考光和信号光. 信号光通过平面镜和离轴抛物面镜的准直和聚焦后, 以30°角斜入射至成像窗口表面, 太赫兹波焦点位于成像窗口下表面, 即待测样品与成像窗口的界面. 反射成像采用材料为石英的矩形成像窗口, 其厚度为0.5 mm. 对于全反射成像, 采用等腰三角形状的全反射棱镜作为成像窗口. 其材料为高阻硅, 大小为34.8 mm × 34.8 mm × 20.0 mm且底角为49°. 携带样品信息的太赫兹波分别经过两个离轴抛物面镜的准直和聚焦后, 最终被太赫兹波探测器接收. 离轴抛物面镜1和2的焦距为2 in(1 in = 25.4 mm), F值为1. 离轴抛物面镜3的焦距为4 in, F值为2. 离轴抛物面镜1和2的离轴角度为30°, 离轴抛物面镜3的离轴角度为90°. 待测样品被置于载物台上, 待测样品上表面与成像窗口下表面紧密接触, 载物台、样品和成像窗口均被固定于x-y线性扫描移动平台(SIGMA KOKI CO., LTD.)上进行逐点扫描成像, 扫描速率约为10 pixel/s. 实验中, 扫描步长设置为200 μm, 温度控制在23 ℃.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2021/11/PIC/11-20210182-4_mini.jpg)
Figure4. Schematic diagram of common path continuous terahertz reflection and attenuation total reflection imaging system.
2
3.2.数据分析
为了降低激光器的功率波动对成像结果的影响, 采用样品信号与参考信号的比值作为样品二维成像的像素值. 考虑到成像窗口的不均匀性, 通过有无样品太赫兹成像图中对应的像素比值(即相对反射率)获得样品太赫兹图像. 将有无样品的太赫兹波成像图分别定义为Fsample和Fsubstrate, 样品的相对反射率为Fsample/Fsubstrate. 便于比较, 反射和衰减全反射模式中样品的相对反射率分别表示为Rr和Ra. 为了比较不同模式的图像对比度, 采用文献[22]中的图像对比度公式:![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20210182_Z-20210603041006-1.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20210182_Z-20210603041129-1.png)
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2021/11/PIC/11-20210182-5_mini.jpg)
Figure5. The resolution of continuous terahertz reflection and attenuated total reflection imaging systems, when the quartz and the high resistance silicon materials are used as imaging windows.
为了对不同成像模式进行对比, 基于高阻硅材料的全反射棱镜获得衰减全反射成像系统分辨率, 如图5所示. 由图可知, 衰减全反射成像的分辨率为400 μm. 与反射成像相比, 衰减全反射模式的成像分辨率较大. 这源于在成像光路和聚焦装置一定的情况下, 太赫兹波在全反射成像中的入射角小于反射成像中的入射角, 使得太赫兹波在全反射成像中的空间角较大. 故选用较高折射率的成像窗口或采用较大的太赫兹波空间角, 均可获得较高的成像分辨率.
进一步, 分析了不同成像模式对待测样品反射率的影响. 根据(3)式和(6)式, 样品的反射率主要与成像窗口折射率、入射角度和偏振有关. 为了避免成像窗口折射率对反射率的影响, 反射和衰减全反射成像窗口均采用高阻硅材料. 反射窗口的厚度均为0.5 mm, 此厚度在检测波长下满足干涉相减. 图6为实验获得连续太赫兹反射和衰减全反射成像模式下, 蒸馏水的相对反射率与入射角θ3的关系, 同时计算了蒸馏水在此双模式下的理论反射率与入射角θ3的关系. 由图可知, 衰减全反射模式下的实验相对反射率Ra与理论反射率RP一致, 如图中蓝三角和绿色线所示. 这是因为衰减全反射是基于样品对太赫兹波的衰减而成像, 当无样品时仅存在全反射棱镜介质的散射和吸收损耗, 对样品来说此过程入射的太赫兹波无衰减; 有无样品信号的比值去除了全反射棱镜衰减对样品检测结果的影响. 然而反射模式下, 实验相对反射率Rr比理论反射率RP大, 如图中红三角和黑线所示. 考虑到反射系统中通常采用金镜获得参考信号, 本实验将金镜置于样品成像处进行成像获得参考信号, 并采用样品信号与参考信号相除的方式获得样品反射率, 如图中紫色三角所示. 实验获得的反射率仍然高于理论反射率, 这是因为连续太赫兹反射成像中获得的样品信号是太赫兹波在反射窗口前表面和后表面反射信号矢量叠加后的强度信号, 反射窗口厚度引起的干涉现象使得准确获得样品反射率较困难. 通过将反射和衰减全反射模式合并在同一个系统中, 可弥补反射成像无法准确获得样品反射率的缺点.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2021/11/PIC/11-20210182-6_mini.jpg)
Figure6. The reflectivity and relative reflectivity of distilled water obtained theoretically and experimentally under the reflection and attenuation total reflection modes.
图7为连续太赫兹反射和衰减全反射模式下水滴的太赫兹成像. 其中, 反射成像采用石英和高阻硅材料的反射窗口且其厚度均为0.5 mm, 衰减全反射模式采用材料为高阻硅的全反射棱镜. 图7(a)—(c)为水滴的可见光图(visual image), 图7(a)—(c)中的水滴面积分别为4 mm × 4 mm, 6 mm × 6 mm和7 mm × 7 mm. 图7(d)—(f)分别为与图7(a)—(c)相对应的太赫兹波成像图. 右边的颜色条表示相对反射率, 相对反射率越低表示样品对太赫兹波的吸收越大. 图7(d)—(f)中成像区域对应的像素值分别为40 × 40 = 1600, 40 × 40 = 1600和60 × 60 = 3600, 检测时间分别约为2.6, 2.6, 6.0 min. 双模式成像过程中, 需要进行反射窗口和全反射棱镜的快速切换及样品重新放置, 此操作过程大约1 min. 总之, 通过此共光路连续太赫兹双模式成像系统获得图7(e)和图7(f)共需要约10 min. 太赫兹波高吸收区域为水滴区域, 不同成像模式均可实现水滴与成像窗口的明显区分. 图7(g)—(i)分别为图7(d)—(f)中白色虚线处对应的相对反射率值. 采用石英和高阻硅材料为反射窗口时, 水滴区域的相对反射率分别为47% ± 3%和65% ± 5%. 根据(3)式, 采用石英反射窗口时水滴的理论反射率高于采用高阻硅反射窗口时的理论反射率, 这与实验结果一致. 另外采用高阻硅反射窗口时水滴相对反射率值的动态范围较大, 这源于太赫兹波在石英反射窗口前表面的反射率(9%)小于采用高阻硅反射窗口的反射率(23%), 石英反射窗口能较好的避免镜面反射. 衰减全反射模式下, 水滴区域的相对反射率为24% ± 1%. 与反射模式相比, 衰减全反射具有较好的信号稳定性. 这可能源于反射模式受干涉现象和表面漫反射影响严重.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2021/11/PIC/11-20210182-7_mini.jpg)
Figure7. (a)–(c) Visible image and (d)–(f) terahertz image of water droplets in different imaging modes, (g)–(i) the relative reflectivity map corresponding to the white dotted line in Fig.7(d)–(f), respectively.
进一步, 基于共光路连续太赫兹双模式成像系统进行生物组织成像研究. 图8(a)为猪肉组织置于石英反射窗口后表面的实物图. 由图可知, 猪肉组织的部分区域并未与反射窗口紧密贴合. 图8(b)和图8(c)为分别采用高阻硅和石英材料为反射窗口时, 猪肉组织的连续太赫兹反射成像. 由图可知, 不同材料反射窗口均可实现肌肉和脂肪组织区域的明显区分, 然而石英反射窗口下的太赫兹反射成像能够更多的反映样品细节信息. 结果表明, 尽管高阻硅反射窗口可获得较高的成像分辨率, 然而石英反射窗口更有助于抑制镜面反射来获得样品信息. 为了对比连续太赫兹反射和衰减全反射成像模式在生物组织成像时的效果, 采用内部具有较小脂肪区域(图8(d)中①处虚线间的距离约为1 mm)的猪肉组织进行成像, 如图8(d)所示. 因高阻硅对可见光的透过率较低, 故仅给出待测猪肉组织的实物图. 图8(e)和图8(f)为分别采用石英反射窗口和高阻硅全反射棱镜, 获得的猪肉组织太赫兹成像图. 与反射成像相比, 太赫兹衰减全反射成像能更好的显示组织中的细微结构信息. 这证明了连续太赫兹衰减全反射成像具有较高的成像分辨率和检测灵敏度. 根据(7)式, 图8(b),(c),(e),(f)的图像对比度分别为0.07, 0.17, 0.17, 0.64. 这证明了反射成像模式采用石英反射窗口有助于增强图像对比度, 衰减全反射成像模式的图像对比度高于反射成像模式. 在衰减全反射成像中, 因成像面不同位置对应的太赫兹波在全反射棱镜中的传输路径不同, 导致衰减全反射成像具有分辨率退化的缺点[28]; 而反射成像中成像面不同位置的成像分辨率相同, 可实现样品大面积成像. 因此, 对于大尺寸样本可首先采用反射成像进行样品初步鉴定, 然后采用衰减全反射成像进行细节准确分析.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2021/11/PIC/11-20210182-8_mini.jpg)
Figure8. (a) Visible image of pork tissue in close contact with the reflection window; (b), (c) the terahertz images of pork tissue in reflective and attenuated total reflection modes using high resistance silicon and quartz materials as reflection windows, (d) visible image of pork tissue when the imaging window is not covered; (e), (f) the terahertz images of pork tissue in reflective and attenuated total reflection modes respectively.