简介:这是一篇发表于《重庆邮电大学学报(自然科学版)》2021年第33卷第2期的太赫兹成像与光谱技术医疗应用综述,面向生物医学、中医药学与智慧医疗方向的研究人员、临床科研工作者及高年级本科生。资源为1个PDF文档,压缩包约2.43MB,内容系统介绍了太赫兹时域光谱与成像技术的基本原理,并重点阐释其在人体皮肤组织、肿瘤组织、骨密度测量、癌症早期诊断以及中草药真伪与质量鉴别等场景中的应用现状;同时分析了设备成本高、成像速度较慢、生物样本适应性不足等现实挑战,并对未来结合人工智能和大数据分析实现精准医疗、手术导航与实时监测的趋势作出展望。目前已有134人学习下载,适合作为医疗光子学、生物医学检测及智慧医疗方向的参考文献或专业指导,可帮助读者快速把握太赫兹技术在医学领域的研究脉络、关键方法与待解决的核心问题。
1. 太赫兹医疗应用:为什么怕水的波段反而能形成诊断对比度
太赫兹最怕水。一个太赫兹脉冲穿过几百微米厚的含水组织,信号就衰减到几乎读不出来;可医学上拿它做病灶成像,靠的恰恰是水。这个结论听起来像玄学,但太赫兹成像和光谱在医疗领域的应用综述里,几乎所有带临床对照数据的实验都指向同一件事:良恶性组织含水量的微小差异,会在太赫兹波段的折射率和吸收系数上被放大成肉眼可辨的对比度。皮肤基底细胞癌、乳腺癌、结直肠黏膜病变,肿瘤区域往往比正常组织高出百分之几到十几的含水量,这正好落在太赫兹最灵敏的探测范围里。这篇博文写给两类人:一类在做医疗仪器立项,想评估太赫兹适不适合作为术中切缘判断或体表病灶筛查手段;另一类是生物光子学方向的研究者,手头有太赫兹时域光谱设备,想拿它测离体活检组织。下文从物理基础、成像与光谱两条路线选型,讲到样本制备、数据处理、常见翻车点,最后给一套可复现的设备验证流程。
2. 成像还是光谱:两条路线的原理、边界和医疗场景匹配
2.1 水是“敌人”也是“对比剂”:太赫兹看病灶的物理基础
太赫兹波段处于微波与红外之间,频率覆盖 0.1 THz 到 10 THz。以 1 THz 为例,单个光子能量只有约 4 meV,远低于分子电离能,所以不损伤 DNA,这是它在医学影像里被反复讨论的前提。水里氢键网络的集体振动模式在 0.5 THz 到 2.0 THz 之间有很强的吸收和显著的介电色散,也就是说,含水量稍有变化,样品的复折射率就跟着变。肿瘤组织的典型特征是血管密度增加、细胞密度增大、细胞外间隙改变,这些都会推高局部水的体积分数。太赫兹脉冲照上去,组织区域比正常区域反射更强、透射更弱,飞行时间也拖得更长。
这里有一个容易误判的点:很多人以为太赫兹医疗检测是在看“分子指纹”或“蛋白质构象”,实际上在体测量和新鲜离体组织里,水信号把分子振动特征压得差不多平了。真正能稳定复现的对比度来自水的介电响应,而不是某个特异性的化学键吸收峰。何时才能看到分子指纹?把组织脱水、压成粉末或者放进低温环境,水的贡献降下来,蛋白质骨架振动模式才显出来。带着这个概念去读综述里的实验曲线,就很容易分辨哪些数据来自新鲜含水样本,哪些来自固定脱水样本。
2.2 成像路线的选型逻辑:透射、反射和太赫兹 OCT 各看什么
太赫兹成像最基本的形式是逐点扫描:把飞秒激光驱动的光电导天线产生的脉冲聚焦到组织上,二维移动平移台,每个像素保存一条完整的时域波形,事后提取特征成像。这个思路的优点是对比度来源可调,你可以画时间延迟图、峰幅度图、0.8 THz 频率点透射图;缺点是采集慢,一张 2 cm x 2 cm 的切片,步进 200 微米,逐点测下来要几十分钟,对离体组织来说这是真实的时间压力。
透射构型适合离体薄片,信号路径短,物理模型是清晰的比尔-朗伯吸收关系,数据反演最稳定。反射构型适合在体皮肤、术中创面和厚组织块,太赫兹波从体表射入,在组织界面和内部折射率不连续处产生回波,但相位受表面平整度影响极大,体表的一个微小水膜都会导致波形畸变。太赫兹 OCT 是另一个分支,利用低相干度脉冲做深度剖面,可以分辨表皮和真皮的层状结构,目前更多用在烧伤深度评估里。选哪条路线,核心判断只有两个:目标组织能不能拿到薄切片,以及你需要的对比度是来自透射损耗还是界面反射。
2.3 光谱路线的选型逻辑:THz-TDS、ATR 和单频成像的关系
太赫兹时域光谱(THz-TDS)提供的是宽频复杂折射率数据,一次测量覆盖 0.1 THz 到 3 THz 的连续频谱,适合做组织成分的定量趋势分析。衰减全反射(ATR)是医疗场景里经常被低估的构型:红外波段的 ATR 设计可以直接搬到太赫兹波段,用高阻硅或金刚石晶体做内反射元件,典型的探测深度在 10 到 50 微米,只接触组织表面。这个特性让 ATR 特别适合两类样本:一类是含水量极高的样本,透射式会因为强吸收信号归零,ATR 只取表面一层倏逝波,把吸收拉到可测范围;另一类是粉末或干态活检组织,想比较不同病理状态之间的光谱细微差异,ATR 的接触压力可控,重复性比透射式气压夹持更好。
单频成像则是在确定某个频率点有区分度之后,把问题简化成“只测这个频点的强度分布”。比如综述里常见的是拿 0.6 THz 到 1.0 THz 之间的某个频率做全图映射,因为低频段水的吸收相对弱,穿透好一些;高频段对比度更高,但噪声也上来。我的建议是先用 THz-TDS 把样品的频谱数据完整扫一遍,确定区分度最好的频段,再考虑要不要换成单频快速成像系统——临床场景里时间有限,没有人愿意在手术室里等逐点扫描二十分钟。
2.4 与现有影像学对比:太赫兹不能替代谁,只能补谁的位
X 光看电子密度差异,MRI 看氢质子弛豫,超声看声阻抗界面,红外光谱看化学键振动。太赫兹的位置在它们中间:它比 X 光安全,又比 MRI 对水含量敏感得多;比红外光子能量低,适合含水组织。但所有综述都会诚实指出穿透深度的问题。太赫兹在含水组织的穿透深度在几百微米量级,超过一毫米信号就几乎没有有效信息了,所以它天生不是做脑部和腹部成像的料,真正发挥价值的是表皮、黏膜、口腔、宫颈、消化道内镜可及区域,以及离体活检组织的术中快速病理。
这决定了评估太赫兹医疗方案时的第一原则:不要拿它跟 CT 比分辨率,也不要跟 MRI 比穿透深度,只看“离体薄切片”和“体表浅层组织”这两个目标换热区。想通这一点再往下做样本制备和数据处理,基本不会走偏。
3. 从综述到实验台:样本制备与三种测量构型的落地选择
3.1 样本制备决定成败:厚度、含水量和病理金标准的对齐
太赫兹透射测量对厚度极其挑剔。太薄,组织与正常组织的光学差异不够,区分度被噪声吃掉;太厚,水吸收饱和,透过率掉到探测器底噪以下,连波形都解不出来。常见做法的经验区间是:新鲜或冷冻离体组织切片取 200 到 500 微米,皮肤和黏膜组织取 200 到 300 微米比较顺手,含有较多胶原的乳腺组织可以放宽到 300 到 500 微米。切片方式上,振动切片优于传统的石蜡切片,因为不需要脱水步骤,能保证样本维持含水状态。
在体测量和离体测量的样本逻辑完全不同。体表反射测的是表皮层,角质层含水量低、微观结构不平整,个体差异很大,综述里的在体数据通常要按部位、年龄、环境湿度做分组,否则基线漂移比病灶信号还大。离体活检组织的测量最好在样本离体后两小时内完成,期间用生理盐水浸润的纱布包裹,避免表面风干。如果做不到新鲜测量,冷冻保存后复温测也还能用,但复温过程要控制速度,避免表面凝结水珠破坏光学接触。
样本的病理金标准对齐是整个实验设计里最容易被轻视的一步。太赫兹扫描得到的是一张光学对比度图,它本身不告诉你哪里是癌,只有把同一张切片做好 HE 染色,由病理医生划出肿瘤区域、间质区域和正常区域,再和太赫兹特征图做空间配准,才能算有效数据。我见过不止一个项目,光学数据漂亮得像艺术品,结果没有对应的病理标注,最后论文里只能写“初步观察”,审稿人一句“how was ground truth defined”就把工作打回原形。
3.2 三种测量构型的参数设置与使用边界
| 构型 | 适用样本 | 厚度/状态 | 可测频段 | 关键参数 | 最大的坑 |
|---|---|---|---|---|---|
| 透射式 | 新鲜离体切片 | 200-500 μm,夹持于石英/TPX窗口 | 0.1-2.5 THz | 时间步长 0.05 ps,时间窗 20-60 ps | 水吸收饱和、窗口间多重反射 |
| 反射式 | 在体皮肤、组织块表面 | 无厚度要求,表面尽量平整 | 0.2-2.5 THz | 焦点位置对准表面,参考面用金属镜 | 表面水膜、焦点偏移引入相位误差 |
| ATR | 粉末、干态切片、皮肤表面 | 接触深度 10-50 μm | 0.5-3.5 THz | 压紧力恒定,晶面清洁 | 样本与晶体间微小空气隙毁掉全反射 |
透射式的窗口材料优先选高阻硅或者石英。高阻硅折射率约 3.4,阻抗匹配不如石英,但机械强度好,反复夹持不易碎;石英损耗略高,透射谱在 2 THz 以上会有明显起伏。夹持时必须在窗口边缘用薄垫片控制压紧程度,压太紧样本被挤出测量区域,压太松形成空气间隙,波形上会多出一个负向小峰。ATR 构型里,高阻硅晶体和金刚石晶体之间的选型标准是:研究含水量高的新鲜组织用金刚石,因为折射率更高、机械硬度大;研究干燥粉末用高阻硅,因为价格差一个量级,损耗更低。压紧力的重复性是 ATR 测组织时最影响结果的一环,压力变了,组织与晶面的接触面积就变,吸收谱基线跟着漂。
3.3 最小可行实验:透射式 THz-TDS 测量一份离体组织切片
搭一套能复现综述里基本趋势的最小实验,不需要直接上临床原型机,一套实验室 THz-TDS 系统就可以启动。测量协议我会按下面几步走,每一步都对应一个具体的检查点。
第一步,测量空气参考。把空载的样本窗口放进光路,保持室内湿度稳定,以 0.05 ps 的时间步长扫描,时间窗总长不少于 40 ps。扫描完成后检查参考波形峰值幅度,如果峰幅度低于系统指标的一半,先检查光电导天线有没有被激光对准,再检查湿度是不是太高——太赫兹波对空气水分很敏感,实验室湿度超过 60% 时参考信号会明显缩水。
第二步,测量窗口背景。将两片干净的石英窗片叠放,测量透过信号。这个信号和空气参考之间的差异,反映窗口材料的吸收和反射损耗,后续样本测量要用这一步做校正,而不是直接用空气参考。很多新手跳过这一步,结果提取吸收系数时出现系统性偏高。
第三步,固定组织切片。把切片从缓冲液中捞出,在滤纸上轻轻蘸一下去掉表面明水,立刻夹入窗口。整个操作控制在 15 秒内。夹好后先目视检查有没有明显气泡,有气泡就重新夹一次。
第四步,采集样本时域波形。同样的时间步长和窗口参数,做 32 到 128 次平均。对含水切片,信号衰减很厉害,128 次平均比 32 次平均的波形会更平滑,代价是测量时间翻倍。
第五步,每张切片至少测 2 到 3 个不同位置,并拍照记录。太赫兹光斑聚焦后直径约几百微米,组织本身不均匀,单点测量结果波动很大,多点平均才能代表整片组织。
第六步,重复测量同一位置,间隔 10 分钟,用于评估脱水带来的漂移。如果两次测量的主峰幅度差超过 10%,说明脱水速度太快,后来的数据必须做时间校正,或者换用保湿环境重新测。这套流程下来,一份样本从开始到结束大约占用 15 到 25 分钟,正好压在离体组织可接受的时间窗口内。
4. 时域波形到诊断特征:光谱数据处理的关键步骤与参数
4.1 预处理:去直流、零点漂移和参考归一化
从 THz-TDS 系统里拿到的原始数据是两组随时间变化的电场幅度波形:参考波形和样本波形。第一步是去直流,把每个波形的平均值减掉,消除探测器的电子偏置。第二步是时间原点对齐,因为样本比空气慢,样本波形的主峰会比参考波形的主峰落后一段距离,这个时间差本身就是重要的物理量,但不能把波形直接硬对齐而抹掉它。正确做法是以参考波形主峰的时刻为时间零点,样本波形保持原始时间戳,后续所有计算都从零点开始取窗。
第三步是判断样本信号是否有效。把样本波形主峰幅度和参考波形主峰幅度做比值,如果比值低于 1% 到 3%,说明组织太厚或者含水量太高,透射信号已经掉到噪声门槛里,继续提取吸收系数没有意义。这种情况下可以改为反射构型,或者直接把样本切薄。
参考归一化的操作是频域除法:分别对参考波形和样本波形做快速傅里叶变换,得到复频谱后相除,得到样品的复透射函数。要注意频域除法对噪声非常敏感,在频谱的两端和水的强吸收凹陷处,分母接近零,相除会放大噪声,所以做除法之前必须确认有效频段的范围,超出有效范围的数据点直接截断。
4.2 FFT 参数怎么设:窗函数、时间窗长和补零
FFT 的参数直接决定光谱质量,这里最容易出现的错误是把分辨率调得太高,看起来曲线细腻,实际全是噪声。频率分辨率等于扫秒时长的倒数。扫描时间窗 20 ps,分辨率为 50 GHz,你要看 0.5 THz 处的吸收,只有 10 个点,够用但不富余;时间窗拉到 60 ps,分辨率提高到 16.7 GHz,0.1 THz 以下的低频段才真正可读。代价是扫描时间变长、组织脱水更严重,通常要在 30 到 50 ps 之间折中。
窗函数选择遵循一个原则:关注幅度谱用汉宁窗,关注时间形状和相位用矩形窗。矩形窗旁瓣泄漏会让吸收谱出现周期性的伪纹,增大频谱的“被子”问题;汉宁窗把主峰展宽了一点,但频谱干净得多。对组织样本,我一般先用汉宁窗做全频段扫描,找到有区分度的频点之后,再改用矩形窗复核该频点的时间响应,确认不是旁瓣泄漏造成的假特征。
补零这个操作在论文里经常被过度使用,本质只是插值,让曲线变得平滑,不提高真实分辨率。很多综述给的频谱图纵轴间隔很细,那是补零后的视觉效果,不是采集时达到了那个分辨率。判断系统真实分辨率的方式只有一个,看原始时间窗的长度,不看 FFT 点数。这一点在对比不同文献数据时尤其重要,两台系统如果时间窗长度不同,低频细节就没有直接可比性。
4.3 折射率和吸收系数提取:公式逻辑与合理性校验
拿到频域复透射函数之后,提取光学常数是按下述逻辑做:折射率主要由相位延迟决定。设样本厚度为 d,真空光速为 c,样本相对空气的相位差为 Δφ(ω),折射率 n(ω) 满足 n(ω) = 1 + Δφ(ω)·c / (ω·d)。时间上更直观的验证方式是看时域的主峰延迟量 Δt,它与折射率的关系是 Δt = d·(n - 1) / c。如果测一个 300 微米的含水组织切片,主峰延迟大约在 0.3 到 0.6 皮秒之间,对应 n 在 1.3 到 2.0 之间;如果算出来的 n 大于 2.5,大概率是厚度测量不准或者相位解包跳了周期。
吸收系数 α(ω) 从幅度信息提取,要先把界面反射损耗扣除。公式的核心形式是 α(ω) = (2/d)·ln[(4n/(n+1)²)·(1/|T(ω)|)],其中 T(ω) 是复透射函数,n 是该频率下的折射率。这里 4n/(n+1)² 是两界面透射系数的乘积修正。没有这项修正,算出来的吸收系数会明显偏高,尤其是在低含水量的干态切片上,误差可达 30% 以上。
相位解包是所有步骤里最玄的部分。时域上组织样本带来的延迟不可能超过一个周期时,相位谱会出现 2π 的跳跃,表现为曲线呈锯齿状。解包算法可以修正,但前提是频点间隔足够密,让相邻频点间的相位差小于 π。如果频点间距大于这个阈值,解包会失败,表现为折射率曲线突然跳变。遇到这种情况,先把时间窗拉长,提高频域采样密度,再重新解包,不要强行平滑曲线。
4.4 成像特征图怎么出:像素特征抽取和病理图对齐
逐点扫描成像的操作要点是每个像素保存完整波形再离线重建。重建时选择的特征量直接决定图像能呈现什么对比度。常用特征有主峰幅度、主峰时间延迟、指定频率下的幅度或相位、指定频段内吸收系数的积分值。主峰幅度图对应组织对太赫兹的总损耗,时间延迟图对应含水量的折射率累计效应,这两个特征在综述里出现频率最高,因为它们对水的响应敏感,实现简单。
重建单频特征图之前,先固定一个频率点,比如 0.8 THz,然后把每个像素频谱在该频率的幅值取出来,按空间位置排布。关键要对准每个像素的取值范围,同一个样本内正常区和肿瘤区在该频点的幅值可能只差 10% 到 20%,伪彩色标尺的上下限要设置得当,否则图像会被一两个异常像素拉满对比度,病灶区域反而看不清。
与病理切片对齐时,一个常见做法是先在透射光下扫描切片和标尺一起拍一张白光照,再做太赫兹逐点扫描,得到太赫兹图之后,用高通滤波或边缘检测找到切片轮廓,和照片里的轮廓做刚性配准。误差控制在 1 到 2 个像素以内就够了,毕竟病理切片的 HE 染色和太赫兹扫描不是同一时刻完成的,切片的细微形变无法完全消除,定性对比为主,逐像素定量在方法论上站不住脚。
5. 太赫兹医疗实验的五个坑:现象、原因与排查办法
5.1 测量中信号幅度越来越弱:样本在窗口里脱水
现象:同一张切片连续测三个位置,第三个位置的波形主峰幅度比第一个低了 15% 以上,折射率曲线也在缓慢向上漂。原因:新鲜组织在开放光路里不断失水,含水量下降,太赫兹吸收减弱,折射率改变。解决:把单样本测量顺序固定为从远到近扫描,把关键位置的测量排在前面;两片石英窗口边缘用少量凡士林密封,延缓水分蒸发;如果样本量大,给载物台加一个小型保湿舱,通入氮气泡石加湿,湿度保持在 80% 以上。每测三个位置,补测一次原始参考位置,用重复测量差做漂移校正。
5.2 频谱上出现周期性抖动条纹:法布里-珀罗干涉
现象:吸收谱曲线叠加了一串等间距的周期波动,间隔大约几十 GHz,看起来像梳子齿。原因:太赫兹脉冲在平行的样本窗口之间来回反射,和主脉冲相干叠加,FFT 后在频谱上形成周期性调制。解决:时间窗口只截取到主脉冲后 5 到 8 皮秒,把后续的多次反射回波切掉;或者把窗口材料换成楔形片,让反射光束偏出光路。如果必须保留长时间窗,就在数据处理里用窗函数把主脉冲之后的尾部逐渐衰减掉。
5.3 时域波形平顶或削顶:探测器饱和
现象:参考波形主峰平坦得像方波,幅值明显低于预期的线性响应。原因:泵浦光功率或者偏置电压太高,光电导天线或探测晶体的响应饱和。解决:先降低泵浦光功率,观察波形峰值是否按比例下降;在线性范围确认之前,不要把增益调满。判断依据是连续调低功率,波形峰值同步下降,说明系统处于线性区。记录数据时,保证最强参考波形峰值不超过探测器饱和阈值的 70%。
5.4 同一张片子重新定位结果差异大:光斑尺寸和组织不均匀
现象:对同一张组织切片测了五次,每次移动后再回到第一个位置,谱线有起伏,尤其是高频段的吸收系数变化明显。原因:组织本身的纤维结构、微血管分布不均匀,聚焦光斑只有几百微米,定位稍偏就落在另一簇纤维结构上;另外,重新夹持样本会改变组织与窗口的接触压力。解决:每个采样点做 3 次重复测量并取平均;用相机标定样本位置,固定样本在窗口上的摆放方位;高频段的吸收系数本身波动就大,看趋势不要逐点抠数。
5.5 折射率曲线在低频端呈锯齿状:相位解包和时间窗起点错误
现象:提取出的折射率在 0.2 THz 以下剧烈抖动,甚至出现负值;时域上样本主峰没有完全落在时间窗内,靠近窗口边缘。原因:时间轴零点设得比参考主峰实际位置早太多,低频相位不能正确展开;或样本信号太弱,在噪声中定位主峰不准。解决:先用峰值搜索算法精确找到参考和样本波形的主峰位置,重新设置时间窗,确保样本主峰位于时间窗起始后的前三分之一;解包前检查相位谱是否单调递增,不单调就先调整时间偏移,不要硬解。这一步做对了,折射率曲线在低频段会稳定很多,材料参数提取的可信度同时也上来了。
6. 守住信噪比和重复性:让太赫兹实验从“能出图”到“能诊断”
6.1 用标准样本给系统做“体检”:硅片、水和 TPX
从综述读到“系统动态范围 60 dB”这种话时,你无法确认它对自己意味着什么,所以每套设备都自己标定。找一片抛光高阻硅,厚度 0.5 毫米左右,用卡尺量准到微米级,测它的透射谱,算出折射率,和公认值 3.4 比较,偏差应该在 1% 以内。再测一次去离子水在室温下的吸收谱,水在 1 THz 附近的吸收系数应该在 100 到 150 cm⁻¹ 量级,峰形平缓。这两个测试分别校验系统的相位准确度和幅度准确度,通过了再谈组织样本。
6.2 从定性区分到定量诊断:ROI 标注、双盲和 AUC
单点光谱或者整图对比度都只是定性趋势,要拿到能支撑论文和立项的结论,办法是把问题变成分类任务:让病理医生在 HE 切片上划出病灶区域,映射到太赫兹图上做感兴趣区(ROI)标注,提取每个 ROI 的均值、方差和特征频点幅值,做双盲分类。用受试者工作特征曲线下面积(AUC)评估区分能力。AUC 低于 0.7 的特征只能算初步观察,0.8 以上才值得继续做。我个人的习惯是先看单频点幅值图有没有肉眼可见的边界反差,再决定要不要上主成分分析做多频点联合分类。多频点模型容易过拟合,样本量少于 30 例时,特征数量超过三个就要提高警惕。
6.3 从单点特征到全像素映射:一个可复用的验证路径
具体技巧是先在一组已知病理的样本上,确定一个或两个区分度最优的频率点和阈值,然后在全像素图像上做逻辑回归映射,输出每个像素的病理概率伪彩图。概率图的边界用病理切片的轮廓做限制,只统计组织区域内的像素。最后一步是拿独立的验证样本集测一遍,不参与建模的数据才有说服力。我刚开始做太赫兹组织测量时,最贵的教训是连续两周的数据因为探测晶片老化出现基线漂移,前一周的数据怎么算都和病理不一致。后来我把“每天开机先测一片硅片,折射率超出 3.38 到 3.43 就停下来检修”写进了实验手册的第一行。做医疗方向的太赫兹研究,最难的不是脉冲走得多远,而是每一步都有可核查的基准,希望这段流程能帮到你。
本文还有配套的精品资源,点击获取