简介这是一篇发表于《重庆邮电大学学报自然科学版》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 μm0.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 GHz0.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/(n1)²)·(1/|T(ω)|)]其中 T(ω) 是复透射函数n 是该频率下的折射率。这里 4n/(n1)² 是两界面透射系数的乘积修正。没有这项修正算出来的吸收系数会明显偏高尤其是在低含水量的干态切片上误差可达 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 就停下来检修”写进了实验手册的第一行。做医疗方向的太赫兹研究最难的不是脉冲走得多远而是每一步都有可核查的基准希望这段流程能帮到你。本文还有配套的精品资源点击获取
