拓十年匠心定制 · 商业建站与技术教学双线并行 咨询热线:400-886-1026 service@lmnt.cn
ARTICLE DETAIL

资讯详情

深耕网站建设与运营推广的一线实战洞察。

偏振敏感OCT与相位成像:让光学相干断层扫描从形态走向功能

偏振敏感OCT与相位成像:让光学相干断层扫描从形态走向功能

光学相干断层扫描(OCT)这些年已经成了很多实验室,特别是眼科和心内科手里的标配工具。可真正把它用明白的人不算多,原因倒不是扫描不会做,而是很多人把OCT当成一台超快显微镜,导出个结构图就结束了。我自己做OCT系统集成和信号处理的时候也走过这个弯路:盯着高分辨率的B-scan看了半天,层次分明,可问题是有些病变在早期阶段,折射率和散射强度几乎没有变化,结构图上完全看不出正常和异常的差别。到了这一步才意识到,我们一直在用的是OCT一半的能力。

所谓增强型成像,就是把另一半能力捡起来。PI 这个词在各类技术文档里通常指向两条路线:一条是 Polarization Imaging,也就是偏振成像;另一条是 Phase Imaging,也就是相位成像。偏振路线通过探测后向散射光偏振态的变化,反推组织内部的纤维结构和退偏特性;相位路线则通过提取干涉信号的相位变化,测量血流、微振动这一类动态信息。两条路线都能让OCT从“只看形态”进化成“同时看结构和功能”,这也是今天这篇东西的核心内容。下面我从物理原理、系统搭建、数据处理到实际应用,把整个链条拆开来讲。

1. 先搞清楚:普通OCT为什么看不清“增强信息”

1.1 结构成像的边界在哪

OCT的成像机理并不复杂:从低相干光源发出的光被分束后,一路进入样品,一路进入参考臂,两者反射回来在探测器上干涉。因为低相干光只会在两臂光程差小于相干长度的地方形成干涉条纹,所以通过对干涉光谱做傅里叶变换,就能重建出样品内部的轴向反射轮廓,也就是A-scan。把光束在横向扫起来,就有了B-scan和三维体数据。

这个过程拿到的本质上是“后向散射强度的深度分布”。强度图上的亮暗,反映的是样品内部的折射率差异、散射体尺寸和浓度,但偏偏不反映分子排列、纤维走向这种更偏组织学的信息。好比拍照只能拍出物体的轮廓和颜色,拍不出这块布料是纯棉还是化纤。病变早期的纤维结构损伤、胶原降解、轴突退行,很多时候并不伴随明显的强度改变,普通OCT自然就被“打回原形”。

1.2 偏振维度里藏着组织学特征

偏振信息和结构信息最大的不同,在于它探测的是“有序性”。人体里大量重要组织具有双折射特性:角膜和巩膜的胶原板层、肌腱和韧带的致密胶原束、视网膜神经纤维层的微管与轴突膜、动脉中膜里的平滑肌与弹性纤维,全都是高度取向、规则排列的结构。当一束光穿过这类各向异性介质时,正交偏振分量之间会积累一个相位差,这个相位差就是双折射产生的延迟量,通常叫retardance。

另外还有一类信息叫退偏(depolarization)。双折射对应的介质还能保持光的偏振有序性,只是改变了偏振态;而很多颗粒状、囊泡状的微结构,比如皮肤里的黑色素细胞、某些钙化组织、脂质沉积,会把光散射得乱七八糟,导致原本偏振态明确的光变成部分偏振光甚至是完全退偏的光。这两个特性都能从OCT的偏振通道里定量提取出来,正好补上强度图缺失的那块拼图。

所以,“增强型成像”的底层逻辑并不玄乎:结构灰度图告诉你哪里有散射差异,偏振延迟图告诉你哪里有纤维排列和形态异常,退偏图告诉你哪里是颗粒性或钙化性结构。三者叠在一起,才是真正意义上的多参数OCT。

1.3 相位维度是另一条“增强”路

相位信息属于另一个方向。干涉信号本身是一个复数,除了幅值,还带相位。只要散射体在两次相邻扫描之间发生了微小的轴向位移,干涉相位就会跟着变化,这个变化量正比于位移量。利用这个机制,OCT可以测出血流速度、组织微振动,甚至做成不依赖造影剂的微血管成像。这一类PI——相位成像,和偏振成像并不冲突,两者共用同一套OCT引擎,只是后端解算的物理量不同。

我在后面专门用一章讲相位路线,因为很多人在做偏振增强时会忽略它,实际上把两者结合起来才是最强的形态。就好比一台车既有方向盘又有油门,一个管方向,一个管动力。

2. 偏振敏感OCT(PS-OCT)的系统架构与器件选型

2.1 描述偏振的三个数学工具:Jones、Stokes、Mueller

想把偏振测量做扎实,数学语言先得统一。常用的三套工具是Jones矢量、Stokes参数和Mueller矩阵。

Jones矢量用两个复数分量描述完全偏振光的电场,操作最方便,但只能描述完全偏振光,遇到退偏就无能为力。Stokes参数则用四个实数(I, Q, U, V)描述光的总强度、水平/垂直偏振差、45度方向偏振差以及圆偏振分量,它不但能描述完全偏振光,还能描述部分偏振光。Mueller矩阵是4x4的变换矩阵,把入射光的Stokes矢量映射到出射光的Stokes矢量,样品的各种偏振特性——双折射、退偏、旋光——都能在Mueller矩阵元素里找到对应关系。

在PS-OCT里,我自己的体会是:Stokes参数是日常主力。因为OCT信号天然带散斑噪声,散射过程又复杂,样品往往同时表现出双折射和退偏特性,Stokes的四分量平均处理比Jones的纯复数模型更鲁棒。Mueller矩阵则适合做完整表征,代价是至少需要多个不同偏振态的入射光,采集和处理都要重不少。

2.2 三种主流架构怎么选

市面上的PS-OCT系统五花八门,但按入射偏振态的管理方式,基本可以归成三类。

第一种是单入射偏振态加偏振分集接收,也叫偏振分集OCT(PD-OCT)。只用一个偏振态(通常是圆偏振或线偏振)照射样品,返回光在接收端由偏振分束器拆成两个正交通道分别探测。优点是结构简单,只需要在标准OCT的样品臂里加偏振控制器、在接收端加偏振分束器,速度也不打折。缺点是只能测相对双折射,对样品慢轴的绝对方向不敏感,而且如果入射偏振态和样品主轴刚好重合,延迟值会掉进盲区。

第二种是多入射偏振态方案。用偏振调制器或电动波片轮流把入射状态切换到两个或多个不同的偏振态,分别采集后再做矩阵重建。它能够同时给出双折射量、光轴方向和退偏度,信息量最完整。缺点是采集时间成倍增加,系统成本和标定难度也上去了。

第三种是扫频源OCT里的波长或时分复用方案。扫频光源本身扫得极快,可以在相邻波长窗口里编码不同入射偏振态,再用同步分时办法拆包。这类系统能做很高的速度,但时序控制、偏振开关带宽和通道一致性都是不小的工程难点。

选型逻辑其实很简单:如果你主要看眼科视神经纤维层,单入射偏振态加偏振分集就够了,省下来的成本可以花在提升扫频速度上;如果你要做斑块分类、组织病理分级这种对光轴方向敏感的研究,直接上多入射态方案,别在第一种架构上反复调参。

2.3 关键元件选型与偏振校准

偏振敏感通路上的核心器件不外乎偏振控制器、偏振分束器、波片和起偏器。偏振控制器非常关键,因为普通单模光纤在弯曲、温度变化时会产生随机的偏振态漂移,光纤写错了、绕线盘弯曲半径太小,都能让最终到达样品的偏振态偏到完全想象不到的地方。做系统集成时,我建议所有光纤尽量用保偏光纤,至少样品臂和参考臂的关键段要保偏;实在要用单模光纤,偏振控制器必须靠近样品端放置。

光路上的波片选型也有讲究。四分之一波片负责把线偏振转换成圆偏振,但波长带宽必须匹配你的光源。如果光源中心波长和波片标称波长差得远,圆偏振度就不纯,后续延迟计算结果会系统性地偏高或偏低。

校准这一步不能省。我会用一块已知延迟量的双折射晶体(比如云母片或石英波片)作为标准样品,先测它的理论延迟曲线,再看系统实测值偏多少。偏差超过5%,优先检查入射偏振态是否纯净,其次看两个探测通道的增益是否一致,再不行就要重新做色散补偿。

3. 从干涉光谱到偏振增强图:完整数据管线怎么搭

3.1 偏振通道的信号处理链路

不管用哪种硬件架构,PS-OCT的底层数据处理都遵循同一套流程:先对两个偏振通道的干涉光谱做波数域均匀化重采样,再做色散补偿和傅里叶变换,得到对应每个像素的复信号。假设接收端是H和V两个正交通道,那么每个像素点都能得到一对复数幅值,记作A_H和A_V。

接下来就是Stokes参数的计算。I等于两个通道强度之和,Q等于强度差,U和V则是两个复信号之间的交叉相关项。这一步要特别小心符号约定,不同的论文对角频率正负和通道方向的定义存在差异,如果你照搬代码不核对定义,出来的相位延迟图可能整个镜像翻转。我自己通常会在代码里先跑一个线偏振标准样,确保计算出的光轴方向跟实际摆放方向一致,再放心处理正式数据。

3.2 延迟、光轴方向、退偏度各自怎么读

当入射光是圆偏振时,Stokes参数和双折射参数之间有非常直接的对应关系:相位延迟量由Q、U、V三个分量的组合经过反正切运算得到,反映的是双折射强度;光轴方向则从U和V的相对相位里解出来,反映的是纤维结构在横向平面里的朝向;退偏度则通过对局部窗口内Stokes矢量长度做平均来估算,值越接近1说明该点的偏振保持越好,值越低说明退偏越强。

图像判读上有一个常见的误区:延迟图里的“亮”不能直接等价于“纤维多”。延迟亮区只代表双折射量大,也就是各向异性强的有序结构;但如果胶原排列混乱、退偏严重,延迟会反而降下来。所以延迟图和退偏图必须对着看,一个太亮一个太暗,往往说明组织已经发生了定向结构的破坏或钙化。

3.3 融合显示和后处理避坑

多参数OCT的显示一般不会只出单张灰图,而是像彩超一样做叠加。主流做法是:结构强度图做底图、用灰度显示,延迟量映射到色调(hue),光轴方向映射到饱和度或方向彩条,退偏度则单独出一张彩色图。这样医生看一张图就能快速判断:哪个区域结构轮廓清晰,哪个区域纤维双折射异常,哪个区域是退偏主导的颗粒性病变。

后处理里最容易被忽视的是散斑噪声。OCT的结构图经过中值滤波、各向异性滤波都还有一定噪声,偏振参数对噪声更敏感。低信噪比的像素点,噪声本身就近似随机偏振,计算出来的退偏度会严重低估,延迟也会飘。我的做法是先做一个强度掩膜,把信噪比低于阈值的像素直接置为无效,只在有效区域内计算偏振参数,否则后期统计全部被拖偏。

3.4 实时性与算力怎么平衡

PS-OCT的数据量比普通OCT翻倍还不止,两个偏振通道意味着同等扫描条件下要处理两套干涉光谱。对于100kHz A-scan速率的扫频源系统,实时计算每个像素的Stokes参数和延迟量,CPU基本跟不上,GPU加速是现在的主流选择。如果你只是离线处理离线数据,用Python的NumPy向量化操作也能接受;但如果系统要在眼科门诊环境里对病人实时出图,我建议直接把Stokes和延迟解算放进CUDA或者OpenCL管线里,同时把强度掩膜和空间滤波也一并搬上GPU,避免主机和设备之间来回传数据。

4. 增强型成像的落地场景:从眼科到工业探伤

4.1 眼科:神经纤维层的双折射与青光眼

眼科大概是PS-OCT应用最成熟的战场。视网膜神经纤维层(RNFL)由大量平行排列的轴突组成,轴突内的微管结构和细胞膜使得RNFL具有明显的双折射特性。多个研究团队发现,青光眼患者早期RNFL的双折射延迟会先于厚度发生变化,轴突损伤导致微管崩解,延迟值下降,但结构图上厚度可能还在正常范围。

这意味着偏振增强可以成为比普通厚度测量更灵敏的早期预警指标。扫描时让入射圆偏振光穿过RNFL,从延迟图上直接算出纤维层的双折射强度,再结合视盘周围扫描的剖面数据,就能做出延迟距离曲线。这种分析不是替代视网膜厚度测量,而是两者互补:厚度看“还剩多少组织”,延迟看“组织还有没有正常结构”。

4.2 心血管和皮肤组织:胶原与斑块

血管内OCT现在越来越普及,用来评估动脉粥样硬化斑块的纤维帽厚度。问题是脂质池和纤维帽在强度图上的区分并不总是可靠,脂质区域的信号衰减容易和纤维组织重叠。偏振增强在这里的优势很明显:纤维帽富含胶原,表现出规律的双折射;脂质池则以各向同性的液态环境为主,双折射极弱;钙化结节则表现出强退偏。三种特征在延迟图和退偏图上几乎是三块截然不同的区域,判读起来远比只看灰度清晰。

皮肤应用同理。皮肤真皮层胶原网络的退偏特性和方向特性,在老年性皮肤松弛、瘢痕增生、肿瘤侵犯深度评估里都有研究价值;黑色素颗粒的强退偏特性也能帮助定位色素性病变。对于烧伤深度评估,PS-OCT可以在不接触创面的情况下区分表层坏死和深层存活组织,这个场景普通OCT很难做到。

4.3 工业无损检测:应力双折射与内部缺陷

别以为偏振OCT只能看生物组织。很多透明或半透明的聚合物材料存在加工残余应力,应力会让原本各向同性的材料产生受力方向的折射率差异,也就是应力双折射。利用PS-OCT扫描注塑件、薄膜涂层和纤维增强复合材料,延迟图直接显示应力分布区域,光轴方向还能看出主应力方向。这比传统的偏光显微镜有很大的优势:不需要把样品切片,可以直接看三维内部。

当然工业场景也有边界。OCT的成像深度受散射限制,透明聚合物能看到几毫米,高散射材料可能连几百微米都穿透不了;再加上需要接触或贴合窗口,检测效率不如某些快速线扫方案。但在无创、高分辨率、内部层析这三个要求同时出现的场景里,偏振增强OCT基本没有对手。

4.4 什么时候别用偏振增强

不是所有样品都需要上偏振。各向同性的均质材料、低散射强度样本、主要靠吸收对比的病灶,偏振信息往往不丰富,硬做只会增加系统复杂度和处理时间。我见过不少项目一上来就追求多参数成像,结果样品根本没有明显双折射,白折腾一圈。

做出这个判断并不难:先在结构图上找高亮、有序、方向性明确的结构;再用一块起偏板配合普通OCT粗略看一下回波随偏振角的变化;变化不明显就果断放弃偏振模块。增强型成像的“增强”二字,永远是建立在结构图像足够稳固、样本物理特性足够丰富的前提之下的。

5. 相位型PI:多普勒OCT和相位方差造影怎么跟偏振互补

5.1 多普勒OCT的测速原理

光学相干断层扫描里的相位PI,最经典的就是多普勒OCT。当散射体沿探束方向运动时,相邻两次A-scan之间的干涉相位会变化,这个相位差和散射体轴向速度成正比:速度越快,单位时间内的相位漂移越大。利用这个公式,就能把静止组织的高反射信号滤掉,只显示运动的血流信号,而且能给出血流方向和相对速度。

一个绕不开的坑是相位缠绕。相位是周期函数,相位差超过π就会被折叠回-π到π之间,测出来的速度会突然从高速“翻转”成反方向高速。解决思路要么把A-scan速率提上去,让每两次扫描之间的真实相移小于π;要么用相位解缠算法把折叠恢复回来。后者在血管曲折复杂、血流方向杂乱时特别容易失败,所以我更倾向于把系统速率冗余留足,让硬件避坑。

5.2 相位方差造影:不做多普勒也能看血管

多普勒OCT有个天然限制:只有当血流方向与探束方向存在投影分量时,相位差才明显;对于垂直于探束方向的血管,相位差接近零,信号就会漏掉。这让微血管网成像变得非常吃力。

相位方差法绕开了这个限制。它不依赖相位差的平均值,而是统计同一位置在多次重复扫描中相位和幅值的波动方差。流动血液里的红细胞不断改变散射信号,方差明显增高;静止组织则方差很低。这个特性对方向不敏感,哪怕是横跨扫描平面的毛细血管也能被识别。现在很多商用OCT系统的血流量模式,底层就是这个原理,它等于把“动态增强”做成了无需注射造影剂的三维血管成像。

5.3 偏振与相位协同:一鱼两吃

偏振PI和相位PI能在同一个OCT系统里同时实现,因为两者共用干涉信号,只是分析维度不同。相位通道拿来做血流和血管造影,偏振通道拿来做纤维双折射和退偏评估。在视网膜疾病场景里,这种组合的价值特别直接:一边看视网膜微血管密度,一边看神经纤维层延迟量,血管和神经的状态同时拿到,对缺血性视神经病变、糖尿病视网膜病变这一类既有血供异常又伴随纤维损伤的疾病,信息维度比单个参数完整得多。

实现协同的额外成本主要是通道数量和数据管线。接收端需要同时支持偏振分集和复数相位提取,后端需要把血管造影算法和Stokes解算并行处理。对现代扫频源OCT来说,算力不再是瓶颈,真正的瓶颈反而是实验设计:怎么保证病人眼球运动不会把相位信号污染,怎么确保扫描轨迹在偏振通道和相位通道之间完全同步。这些细节,做实验之前就得在方案里写清楚。

6. 自己搭系统时踩过的坑:排查链路与实用建议

6.1 偏振延迟测量不准时,按什么顺序查

我在调试PS-OCT时遇到最多的现象是:标准样品的延迟值偏大或偏小,重复性还差。这时候千万别先怀疑算法,而是按链路从前往后查。

第一步查入射偏振态。拿一个线偏振片放在样品臂末端,旋转检偏器,看探测器输出是不是标准的正弦变化;如果不是,说明光纤里的偏振态已经被扭曲了,先调偏振控制器。第二步查通道一致性。两个偏振探测通道的响应如果不匹配,Q值就会带系统误差,延迟必然不准;在样品位置放一个反射镜,记录H和V通道对同一线性偏振光的响应,做增益归一化。第三步查色散补偿。偏振通道之间的残余色散会让Stokes参数的交叉项出现额外的相位项,导致延迟计算中出现周期的震荡条纹;这种条纹往往在图像上表现为规则排列的明暗交替,很好认,也很好修,重新拟合色散系数就行。最后再回到标准样品验证。

如果以上都正常但图像仍有区域性噪声异常,就要检查是不是样品表面的强反射造成了偏振通道之间的串扰。处理办法要么降低参考臂功率避免探测器饱和,要么在算法里把浅层强反射区域做掩膜,不要让它污染后端的深度统计。

6.2 一套实用的每日开机校准流程

偏振敏感系统不像普通OCT那样“开机就能扫”。参考光路里的光纤会随温度和使用时间发生变化,我就是吃过这个亏的人,前一天晚上校准好的状态,第二天早上起来延迟图整体漂移了几个百分点。所以养成每日校准的习惯非常重要。

我的流程大致如下:先做暗电流和背景光谱采集,然后让系统空扫空气路径,确认两个通道的噪声底大致一致;接着放标准波片样品,连续采集100帧,计算延迟均值和标准差,均值偏离理论值超过5%就微调偏振控制器;最后采集一个活体或仿体样品的快速扫描,确认结构图、延迟图、退偏图三张图都正常。整套流程十分钟内完事,却能帮你躲掉大半天返工。

6.3 入门路线建议:不一定要从零造系统

给刚开始做PI增强OCT的同行一个实在建议:别急着从零搭一套PS-OCT硬件。商用的普通OCT系统大多预留了光纤接口,你可以在样品臂前面串一个四分之一波片加偏振分束器,在接收端把正交通道分离出来,先做line-of-sight的简易偏振测量。这个外挂方案虽然牺牲了一部分成像质量和解析能力,但足以帮你建立对偏振数据的直觉:延迟图怎么读、退偏图为什么噪点大、光轴方向为什么在血管区域会乱跳。

等技术直觉有了,再去啃全系统集成或者直接采购偏振OCT平台,就顺理成章得多。处理上建议先拿仿体验算法,用已知双折射的透明聚合物薄膜做标准数据,把Stokes计算流程跑通,再碰生物组织。生物组织的退偏和多重散射会把波动放大,跳过仿体阶段直接上组织数据,几乎必然在调参泥潭里耗掉大量时间。

6.4 最后的经验总结

剥开那些精密的硬件和花哨的图像,偏振增强和相位增强的PS-OCT并没有想象中那么高不可攀。它比普通OCT多出的,其实是一套更完整的物理视角:强度看散射,偏振看秩序,相位看运动。这三页拼起来,才是一个更接近组织真实状态的三维画像。

我个人的体会是,这个领域的项目成败往往不取决于某一个高精尖器件,而取决于你能不能把偏振态控制得稳、把通道校准做得勤、把掩膜和噪声处理想得早。硬件商会把系统说得再成熟,实际调试时每一个光纤弯曲、每一组增益失配都会跳到脸上。从最简单的偏振分集通道开始,一步一步把每个环节的误差摸清,等这些坑都填平了,多参数OCT给你的回报会比想象中大得多。

返回列表