转自公众号:photonix
撰稿人 | 薪胆居士
论文题目 | spatially offset optical coherence tomography: leveraging multiple scattering for high-contrast imaging at depth in turbid media
作者 | gavrielle r. untracht,mingzhou chen,philip wijesinghe,josep mas,harold t. yura,dominik marti,
peter e. andersen*,kishan dholakia*
完成单位 | 丹麦工业大学,英国圣安德鲁斯大学
研究背景
光学相干层析成像(optical coherence tomography, oct)是一种利用信号光和参考光之间的干涉获取深度方向的信息,通过信号光在待测目标上的扫描来实现生物组织或特种材料的三维层析成像的技术。在日常生活中,oct一般多见于医院的眼科,在皮肤科、心血管科和牙科等科室也有所应用。眼科医生使用oct对视网膜(尤其是黄斑区)的内部结构进行成像,通过所得层析图像对眼底出血、视网膜裂孔等疾病进行诊断的同时,也作为后续治疗的依据。
oct一般由光源、干涉光路和光电探测器组成。光源发出的光被分为两束,信号光由光学系统聚焦到待测样本内,经内部微观结构的反射和散射以及组成成分的吸收后被成像透镜收集,最终与参考光在光电探测器上发生相干叠加。光电探测器将干涉产生的光信号转化为电信号,经过数据采集进入电脑加以分析。通过调整参考光的光程,可以使其与来自样本内不同深度的信号光发生干涉,实现在深度方向上的扫描。移动信号光的入射位置可以实现对样本的横向扫描。样本内部的光学参数分布,表现为对光的吸收和散射不同,光电探测器所得的电信号就有所差异,最终体现为oct图像上的明暗分布。oct的轴向分辨率由光源的相干长度决定,一般在1至10微米之间;径向分辨率与普通光学显微镜类似,取决于样品内部聚焦光斑的尺寸,一般也在微米量级。
论文导读
基于oct的形态学成像技术在过去几十年中取得了重大进展,在医学、艺术品保护、材料的无损检测等领域有着重要的应用。目前oct技术最主要的瓶颈在于测量深度有限:一般为1毫米左右。造成这一问题的原因在于待测样品(包括生物组织)对光的吸收和散射。吸收是指光在样品中传播时,与其组成成分相互作用,被部分地转化为其他能量形式的物理过程。散射是指光通过不均匀介质时,一部分光偏离原方向传播的现象。这两者都会造成沿着‘原入射方向传播的光’的衰减,导致光只能在一定厚度的介质内传播,因此也就只能对一定厚度的样品成像。
参与oct成像的是进入样品内部经过吸收和散射后原路返回或者以一个香蕉形路线从样品表面再次出射的光(被称为后向散射光)。后向散射光可以根据散射次数的多少分为两大类:弹道散射和多重散射。弹道散射光的散射次数较少,传播方向与入射光路夹角很小,比较容易被成像光路收集,因此大部分弹道散射光都能参与oct成像。多重散射光的传播方向则比较分散,随机性较强,只有很少的一部分参与成像。因此从传统的观点来看,oct图像主要是由弹道光子组成的。但从成像深度的角度来看,散射次数和入射深度成正比:能进入样品深层的光大多发生过多重散射,而弹道散射通常发生在浅层。这也是限制了oct成像深度的原因之一。目前主流的研究方向是提高弹道光子的入射深度,常见的方法有波前整形等。丹麦工业大学的peter e. andersen与英国圣安德鲁斯大学的kishan dholakia等人则另辟蹊径,提出了修改oct成像系统结构、避开弹道散射光、只利用多重散射光成像的思路,提高了oct图像的成像深度及在纵深方向上的能见度。相关成果以spatially offset optical coherence tomography: leveraging multiple scattering for high-contrast imaging at depth in turbid media为题发表在science advances上。
技术突破
该文借鉴了暗场oct、双轴oct等通过修改入射和成像两条光路的相对几何形状来调整反射光、弹道光和多重散射光之间比例的研究成果,尝试仅利用多重散射光生成oct图像。具体实现思路如图1b所示:在入射光路和成像光路之间加入横向偏移将两者分开,从而避免收集到弹道光。作者将这一方法称为空间偏移oct(so-oct),并通过实验分析和理论建模对这一思路进行了系统研究。
图1 传统oct、so-oct和另两种样本臂的几何结构比较。(a) 常规oct的几何结构。(b) 本文所述so-oct的几何结构。s0是入射和成像两条光路之间的横向偏移。(c) 双轴oct,入射和成像光路有一定夹角和重叠。(d) 角度so-oct是so-oct和双轴oct的结合,保留了光路夹角,消除了重叠部分。f,透镜焦距;d,从透镜到样品表面的距离;n1和n2,折射率;α,两条光路之间的半角偏移。
作者实验所用的两套oct系统如图2所示。图2a是一套自行搭建的测试系统,中心波长为800nm,通过透镜l2的平移来实现并控制入射和成像两条光路之间的横向偏移。这一套系统主要用来进行系统研究。图2b是在一套thorlabs公司的频域oct的基础上改造得来的,用以证明作者的方法可以在商用系统上通过一定改进来实现。
图2 so-oct的实验装置图。(a) 自行搭建的so-oct系统的实验装置图,波长为800nm。入射光路和成像光路之间的横向偏移由透镜l2的平移产生和决定,图中双箭头代表了l2的移动方向。右下角插图表明了两条光路之间的偏差,双箭头代表l2移动时偏差变化的方向。(b) 改造商用频域oct实现的so-oct系统示意图,波长为1295nm。两条光路之间的横向偏移通过调整孔径光阑p1的位置来决定。双箭头代表了p1的移动方向。sld,超辐射发光二极管;smf,单模光纤;o1,显微物镜;g1,二维振镜;b1到b4,thorlabs公司的非偏振分束立方;l1到l6,透镜;nd,中性密度滤光片;dc,色散补偿;m1,反射镜;ccd1和ccd2,相机;gr1,透射光栅;cir,环形器;rr1,反射棱镜;ap1,可调孔径光阑。
在so-oct中,作者将后向散射光分为三类:弹道光、多重散射光和漫散射光,其中漫散射光组成图像的背景。作者通过扩展的惠更斯-菲涅尔方程(extended huygens-fresnel ehf)阐述了利用入射与成像光路之间的横向偏置分离多重散射光的理论依据。通过理论分析,oct系统的外差效率因子ψ可表述为:
式中的中括号内第一项对应弹道光、第三项对应多重散射光,第二项则是弹道光和多重散射光的混合项。
图3 通过数值模拟分析常规oct和本文所述so-oct的外差效率因子(ψ)。(a) 常规oct中弹道散射(黄线)和多重散射光(蓝线)在图像成型过程中的相对贡献随成像深度变化的曲线。以图中箭头为分界点,浅层图像中弹道散射贡献比较大,深层成像则以多重散射光为主。(b) 成像过程中多重散射开始占主导地位的深度随着样品散射系数变化的曲线。(c) 当so-oct中两条光路之间的横向偏移取不同值时,ψ随着深度变化的曲线。图中圆圈标出了图像对比度最大的深度位置。(d) so-oct中的归一化对比度,即归一化外差效率因子。图2b所述系统的数值模拟参数如下:波长λ0=1295nm、光强为1/θ处的半径ω0=1.5mm、焦距f=36mm、折射率n=1.4、散射系数μs=10mm、吸收系数μa=0mm、θ=0.3rad。
图3a展示了常规oct中,弹道光和多重散射光对成像过程的贡献随着成像深度变化的曲线。结果表明弹道散射在浅层成像中占优,但在深层成像过程中,多重散射光的贡献更大。如图3b所示,作者经过系统研究发现,多重散射取得优势地位的‘成像深度阈值’主要由样品的散射系数决定。入射与成像光路之间的横向偏移s的变化对ψ的影响如图3c所示。结果表明可以通过调整s来优化不同深度的成像质量,提高对比度噪声比(contrast- to- noise ratio, cnr)。
作者分别用生物密封膜nescofilm和透明胶对so-oct的成像效果进行了研究。图4a到e所用样品为七层nescofilm(单层厚度为120微米),实验装置如图2a所示。实验时,对应不同深度成像的对比度需求,光路偏移s有四种取值:0、25、50和75微米。与之对应的光电探测器曝光时间分别为:30、48、250、1600微秒,以保证多次测量时同一层的oct图像亮度一致。从结果可以看出,随着偏移s的增加,‘oct图像随着深度的衰减’变小,图像质量提高。根据图4a的线性梯度,可以得出so-oct在s=0、25、50和75微米时,衰减系数分别为5.5、5.2、3.5和1.2 mm,最多减小了4.6倍。
图4f到k使用图2b的oct系统测量多层透明胶。因为该系统是商用产品,没有曝光调整功能,所以将所有实验结果分别以第一层透明胶的亮度为基础做归一化来进行横向对比。用这种方法来做处理时,信号和噪声同比缩放,影响了cnr的改善,导致偏移s的取值范围小于图4a对应的实验,而且图4g到k的图像优化效果看着不明显。但从图4f同样可以看出,‘oct图像随着深度的衰减’变小。
图4 用多层膜样品测试so-oct的成像质量。(a) 样品各层的oct图像平均亮度随成像深度变化的曲线。样品由七层生物密封膜nescofilm制成,厚度约1mm。使用图2a所示实验装置。虚线为线性拟合。右上角图例为入射和成像光路之间的横向偏移s(offset)。(b到e) 是横向偏移s取值分别为0、25、50和75微米时所得oct图像。(f) 多层透明胶带实验结果。使用图2b所示实验装置。(g到k) 分别对应偏移s取值0、25、30、35和40微米。比例尺为200微米。
根据光的散射原理,当散射粒子的直径远小于光的波长(d≤λ)时,散射特性符合瑞利散射理论,光会向各个方向散射;当散射粒子的直径与波长有一定可比性(d<λ)时,符合米氏散射,主要是向前方和后方散射;当散射粒子的直径大于波长时(d>λ)时,用几何光学来表述,光会在一个比较宽的范围内反向散射。因此样品的散射特性也可能会影响so-oct的性能。作者在黄油块中200微米深度左右的一个平面内混入大量聚苯乙烯微珠(直径100微米)作为样品,使用图2a所示800nm波长oct展开实验。其中聚苯乙烯微珠直径远大于波长,产生大范围反向散射;黄油蕴含的脂肪分子则符合米氏散射,即主要是前向和后向散射。实验结果如图5所示,a到c为常规oct图像,黄油和微珠混在一起,难以分辨;d到f为so-oct图像,充分体现出了黄油和微珠的散射特性差异,也证实了so-oct确实是以多重散射光为主进行成像的。
图5 so-oct根据散射特性利用多重散射成像提高对比度。(a到f) 在黄油中200微米的深度加入直径100微米的聚苯乙烯微珠,采用800nmoct采集图像。(a和d) s分别为0和80微米时采集到的平面图像。(b和c) a中红线和紫线处的剖面图。(e和f) d中红线和紫线处的剖面图。绿色和黄色圆圈表示与常规oct相比,so-oct中的微珠和背景对比度更高更清楚。
作者使用斑马鱼测试so-oct对柔软的生物组织成像的能力。将斑马鱼用福尔马林固定并置于具有磷酸盐缓冲盐水溶液的培养皿中,用800nmoct系统对背肌和软骨棘成像,结果如图6所示。so-oct图像(图6b)的对比度明显优于常规oct的图像(图6a),背肌和软骨棘的对比度更高,软骨棘的结构更加清晰。
图6 斑马鱼实验证实so-oct可以有效提高组织层之间的cnr。(a) 常规oct图像。(b) so-oct图像。图中方框分别对应背肌和软骨棘,结果表明(a)的cnr为0.45db,(b)的cnr为1.29db。
作者将斑马鱼嵌入加有直径0.8和0.5微米聚苯乙烯微球的琼脂糖凝胶中来对比常规oct (s=0微米)和so-oct (s=50微米)在深层成像方面的性能,结果如图7所示。图7a和e显示,浅层成像时(深度50微米)常规oct较好。随着深度的增加,so-oct的成像效果逐渐反超。根据浅层成像以弹道散射为主、深层成像以多重散射为主的理论基础,再次证明了so-oct主要利用多重散射光成像。从图7a到d可以看到随着深度的增加,传统oct图像的衰减很明显,而so-oct的衰减则要缓和的多。此外,图7d(深度350微米)的传统oct图中能看到浅层结构的‘阴影’,而图7h中则不存在这样的问题。从图8对cnr的分析可以看出,在250微米左右so-oct的cnr提升效果最佳,达到了传统oct的8.8倍。
图7 嵌入在散射体凝胶中的斑马鱼的实验结果。实验装置为800nmoct系统。成像深度分别为50、200、350和500微米。红色和蓝色矩形表示用于计算图8所示cnr的感兴趣区域。绿色比例尺表示长度为500 微米。(a到d) 常规oct在不同深度处的平面图。(e到h) so-oct在不同深度处的平面图。
图8 嵌入在散射体凝胶中的斑马鱼实验中cnr随成像深度变化的曲线。蓝色为常规oct,红色为so-oct。
在oct领域中,因为骨骼和软骨这一类硬质生物材料的强散射性,导致成像深度严重受限,所以相关研究很是少见。但是so-oct的主要潜在优点之一是有机会改善高散射介质oct图像中的cnr。因此,作者使用离体的小鼠股骨来展示so-oct在硬质生物组织成像方面的优势和潜力。该实验使用1295nmoct系统和40微米的偏移,结果如图9a和b所示。首先,so-oct图像中深度方向上的信号衰减小于常规oct图像。骨的上表面和下表面之间的衰减减少超过24 db。其次,硬骨和骨髓之间的边界更加清晰,骨的下表面可见度更高。最后,作者将so-oct系统的偏移分别取0、25、30、35和40微米时,图像在横向上的平均亮度的归一化数值随成像深度的曲线绘制出来,如图9c所示。总体的信号衰减随着偏移量的增加而减小,与图3d所示规律一致。
图9 用oct系统采集的离体小鼠股骨的b扫描图像。中心波长为1295 nm,图2b所示系统。(a) 传统oct。(b) s=40微米。so-oct展示出增强组织层之间的对比度的能力。图中黄色箭头表示硬骨和骨髓之间的边界,红色箭头表示骨的下表面。(c) 光路偏移取不同值时,横向上的平均信号强度随深度变化的曲线。标记为“thorlabs”的线表示通过原商用设备收集的信号。
观点评述
成像深度问题仍然是oct技术的主要瓶颈之一。目前主流的研究思路是通过减小多重散射,提取深层的弹道光子或提高弹道光子的传播深度来增加成像深度。该文作者提出了另一种思路:通过改进光路直接利用样品深层的多重散射光来实现oct成像,从而提高成像深度。
作者不仅仅搭建了so-oct实验平台,对这一新思路进行了验证和系统研究,还改进了商用系统,指明了对已有oct系统进行升级的可行路线。在实验上兼顾了图像强度随深度衰减的问题、散射特性对系统性能的影响问题,甚至还包括了不常见的骨骼和软骨这样的硬质生物组织成像问题。实验结果展示出so-oct具有减小‘图像强度随深度的衰减’,提高深层成像对比度,利用散射特性对特定散射介质成像等优势和应用潜力。
作者在文章最后提出,根据该文的原理可以通过使用功率随深度增加的结构化光束 (如贝塞尔光束),来达到类似的成像效果;还可以与自适应光学或光束整形等提高光入射深度的技术相结合进一步提高成像深度。该文的研究成果为oct的发展指出了一个新的研究方向。
本文出处
发表于:science advances
论文链接:
https://www.science.org/doi/10.1126/sciadv.adh5435
更多原文内容,请点击“阅读原文”
推荐阅读
关于photonix
photonix 属同行评议、开放获取(oa)高影响力国际期刊。是中国光学工程学会会刊,由中国光学工程学会、清华大学、上海理工大学和西湖大学共同主办,由springer nature集团出版。上海理工大学顾敏院士和西湖大学仇旻教授担任期刊主编,庄松林院士担任期刊名誉主编。期刊拥有强大的国际编委和编辑团队。photonix 主要报道国内外光学与光子学技术与信息、能源、材料、生命、精密制造、纳米、光电子器件、微纳米电子等学科交叉融合发展带来的颠覆性科研成果和最新的工程应用进展。以展现具有前沿性、多学科交叉和衍生性特点的技术为核心,成为推动国际前沿“使能技术”的平台。
photonix 已被sci、ei、scopus、doaj、proquest、cnki、inspec、dimensions等10多个数据库收录。2022年6月获得首个影响因子19.818,位列q1区。同时进入《2022年中国科学院文献情报中心期刊分区表》,位列物理与天体物理大类和光学小类双一区,为top期刊。中国科协首次颁布“光学工程和光学领域高质量期刊目录”photonix 位列t1级。
了解photonix最新动态
阅读原文
查看全文
声明:本文所用视频、图片、文字如涉及凯发娱乐的版权问题,请第一时间告知,我们将根据您提供的证明材料确认凯发娱乐的版权并按国家标准支付稿酬或立即删除。邮箱:wanghaiming@csoe.org.cn
公众号投稿请联系:wanghaiming@csoe.org.cn
商务合作/文章转载请联系:022-59013419
加入学会会员
中国光学工程学会
邀您成为会员
中国光学工程学会是中国民政部批准的全国性一级学会。会员涉及光学和光学工程领域的科学家、工程技术人员、教授、学者,以及商界精英。学会秉承“创新、协同、求真、务实”的精神,将通过学术交流、成果展示、成果转化等多种形式,搭建一个产学研大平台,推动光学工程的技术创新、注重光学科技的成果转化、面向应用促进光学前沿的工程化。中国光学工程学会推荐年轻的领军专家作为学会核心力量,培养一批国内外新领军人物。
欢迎分享
↓点赞
↓在看
文章转载自微信公众号:光学工程