Design of time domain full field OCT imaging system based on LED illumination

光学系统的相干层析成像技术(optical coherence tomography，OCT)在迈克尔逊干涉的理论上发展起来，通过采集相干光与物体的散射光汇聚产生干涉后的特定信息，实现对所需物体特定位置的层析成像，具有非接触、无损伤、分辨率高、探测率高^[1]等优点，可用于视网膜成像^[2]、胃肠内窥镜^[3]、泌尿科等^[4]，还能进行光学薄膜长度标定^[5]、电路板刻制等^[6]，在生物医学、物理学、材料学等方面得到了大范围的使用^[7]。

基于步进电机控制的时域OCT系统由于其方便、快速、精确而被用在各个领域^[8]，但如今OCT的应用也遭遇了一些困难。其中一点是扫描的速率亟待继续加快，以适应同步成像的需要；另一点是整个系统实现的成本还亟待继续降低，但在此基础上要满足适当效率的扫描显示功能，使其能用于欠发达地区^[9]。目前，OCT一般选宽带光源，如白光以及超辐射发光二极管(super-luminescent diode，SLD)等，白光频带较宽，精度高，但是稳定性较差。SLD的探测深度大，但相干长度不够短，导致系统的分辨率降低。华中科技大学的WU等人曾使用波长为1.3 μm、功率为20 μW的通信用单模发光二极管(light-emitting diode，LED)作为OCT系统的光源，并使用该系统成功对人体皮肤进行了扫描测量，但其轴向分辨率只有20 μm^[10]，与目前主流OCT系统的10 μm左右的轴向分辨率还有一定差距。2020年，NANDAKUMAR等人使用波长为670 nm、带宽为20 nm的LED作为光源搭建了时域全场OCT系统，并使用程序代码进行图像的后续处理，且避免使用昂贵的器件搭建系统，实现了以较低成本获取高分辨率的断层扫描图^[11]。但该系统对干涉信号使用直接采集的方式，实际扫描过程中很容易受到零级像与共轭像的干扰，需要后续再进行图像的降噪处理等操作(而本实验中使用待反馈的闭环四步移相方法采集信号，可以避免此干扰)。目前，更多高性能、低成本的光源仍有待研究，并应用于光学层析成像系统中^[12]。

本文中提出了一种基于LED照明的时域全场OCT成像系统，使用LED作为系统的入射光。与常见的SLD相比，LED有着成本低廉、损耗小、耐用性长等特性，且LED有着优秀的时间相干性，所以也被允许使用在OCT上。除此之外，LED不是激光，而且它的光谱较宽，扫描时不会受到散斑噪声的制约，所以在OCT系统中利用LED作为入射光可以改善图像的品质，还能减小系统成本。本文中的信号采集方面使用带反馈的闭环四步移相法。最早由PATIL等人提出了四步甚至多步移相的广义定义，但四步移相法是在干涉测量中使用最多的一种，移动的次数较少，同时又保证了系统的稳定性以及抗噪声的能力。

相比传统的点扫描OCT，时域全场OCT可以省略横向扫描过程，减少扫描时间^[13]。时域全场OCT采用宽场照明、平行探测方式，适于高散射或多层样品的成像^[14]。此外，本文中使用带反馈的闭环四步移相法进行信号采集，比普通的开环、不带反馈的移相法更加精确，另外使用平面镜、洋葱内表皮进行了扫描实验，得到了较为理想的2维层析图，证明了该扫描显像技术的合理性和准确性。与已存在的方式比较，该方式不仅进一步减小了成本，还兼具了更快的扫描速率，并确保了较理想的图像清晰度，具有很大的实用意义。

图 1是本文中提出的基于LED照明的时域全场OCT系统。图中，a和b是两个非球面透镜，c和d是两个聚焦物镜。透镜b之后放置一个孔径光阑，可以通过调节不同直径的小孔来调整光束的大小，可以进一步准直光束，并提高光场的空间相干性。由于参考臂上的平面镜与样品臂上的生物样本的反射率不同，所以在物镜c与d之前分别放置一个可调衰减器，用以平衡两臂的光强到同一水平。

Figure 1.  Time domain full field OCT system based on LED illumination

本系统基于Michelson干涉仪结构，主体是由宽带光源、分束镜、样品臂、参考臂和探测器组成。系统入射光用的是索雷博公司型号为M810L3的LED，它的中心波长为810 nm，其带宽为25 nm，相干长度为23 μm，LED的输出功率为32 mW。物镜采用的是索雷博公司型号为RMS10X的Olympus平场消色差物镜，可以提供10倍的放大倍率，此外，它有着平阔的视场和对近红外光谱中波长的像差修正，所以它有着优质的球面像差和色差修正以及十分完美的视场平阔度。除此之外，它具有高数值孔径(0.25)，可以获得高空间分辨率。

在光束路径上，LED光的强度一般会持续降低且会发散^[15]，若安装未作处理的LED作为入射光，会影响系统的扫描质量，所以有必要对LED出射的光采取准直操作。因此，本系统中添加了两个非球面透镜以及一个毛玻璃，用于对LED的出射光实施整形，如图 2所示。LED光进入第1个透镜之后在毛玻璃上聚焦为一点，然后经过第2个透镜实现进一步整形。在透镜后的感光片上能够发现准直之后的LED光束，而且由于添加了毛玻璃，光束分布更加均匀，这就完成了对入射光的准直。

Figure 2.  Structure diagram of a LED beam collimation

四步移相法是采用光学干涉术与数字相移技术融合的一种相位干涉测量方法，采用零差探测模式，使用步进移相的方法进行解调处理，信号处理简单，信息存储量小，因此非常适合应用于OCT系统中进行信号采集^[16]。其原理是使用步进电机在一定范围内移动参考镜，轻微改变两臂的光程差，让产生的干涉条纹同时进行对应的偏移，使信号的光强按照余弦规律变化，在探测器端就可以接收到多幅具有不同相位的干涉图像，从而能够避免零级像和共轭像的阻碍^[17]。四步移相法还有个好处是仅对光探测到的一层样本成像，而另外层的样本仅是背景像，等同于避免了非相干层样本的信息^[18]，可以更高效率地完成层析成像。

探测器所记录下的干涉场强度分布为：

式中，I₁(x，y)、I₂(x，y)分别为两个相干光的光强，φ(x，y)为待测样品层的主相位分布函数，Δφ为参考镜的变化相位位移。

假设背景光强为I₀(x，y)，γ₀(x，y)为调制度，则(1)式中干涉场中某点的光强可表示为：

针对样品的某一被测状态，每隔π/2采集一次干涉信号，在0~2π内一共采集到4幅干涉图，分别为0，π/2，π，3π/2，对应得到4个干涉条纹的光强分布见下^[19]：

可得产生的干涉信号中的样品层的主相位φ₀(x，y)为：

振幅为：

每次平移台改变λ/8的长度，相位就会更换λ/4，其中λ为光源的中心波长。当更换了3次相位，就实现了一次完整的四步移相过程，在探测器中保存得到的干涉条纹图，使用上述算法计算，就能够获得样本的相位和振幅。

本系统采用的LED的波长为810 nm，带宽为25 nm，则根据公式:

式中，λ为LED的中心波长，Δλ为LED的带宽，可求得系统的轴向分辨率约为11.8 μm。

OCT的横、纵向分辨率是互不相关的，纵向分辨率主要由光源的波长以及带宽决定，而横向分辨率除了决定于光源的波长之外，还由物镜的焦距及其上的光束尺寸共同决定，即：

式中，f为物镜的焦距，d为光斑尺寸。由(10)式可知，系统的横向分辨率约为19.8 μm。

本实验中，首先使用平面镜作为扫描对象，固定在样品臂的平移台上，使用程序以及驱动完成平移台位置的改变，通过调整参考臂的距离，使从样品臂和参考臂反射回来的光在分束镜的位置出现干涉，使用四步移相法得到的4幅干涉图如图 3所示。实验中使用的互补金属氧化物半导体(complementary metal-oxide-semiconductor，CMOS)相机的曝光时间为2.15 ms，即采集1幅样品图的时间为2.15 ms。

Figure 3.  Plane mirror interferogram

图 3a~图 3d分别为由CMOS相机获取的4幅干涉条纹图。每幅干涉图的相位间隔202.5 nm波长。改变参考臂的位置，依次改变101.25 nm，实现完整的四步移相。通过(7)式和(8)式计算得到4幅干涉图，其扫描平面镜的相位图和振幅图如图 4所示。

Figure 4.  a—phase diagram of plane mirror  b—amplitude diagram of plane mirror

由于平面反射镜结构简单、反射系数高，所以干涉条纹比较均匀。从图 4a可知，相位条纹分明，表明四步相移完成的效果较理想。

然后，将成像样品换成一层洋葱内表皮，继续使用四步移相法测量。洋葱鳞叶表皮细胞构造直观，层次多样，十分方便OCT成像。将洋葱划开，割裂一块鳞叶，用镊子取下一小块表皮，然后放在载玻片上进行成像，得到的4幅干涉图如图 5所示。其中图 5a~图 5d分别为相位间隔202.5 nm波长的干涉图。

Figure 5.  Interferogram of the inner epidermis of a layer of onion

由于洋葱组织结构厚度不均匀，因此样品每个位置的光程长度都不相等，从而产生不规则的干涉条纹。根据计算得出的相位和振幅图如图 6所示。

Figure 6.  a—phase diagram of the inner epidermis of a layer of onion  b—amplitude diagram of the inner epidermis of a layer of onion

从图 6a中可以看到，计算出的相位图比较清晰；从图 6b中可以较为清晰地看到洋葱的组织结构，其中A处为细胞，B处为细胞壁，细胞壁的边缘清晰可见，说明对一层洋葱内表皮的成像效果较为理想。

对相位图进行解包运算，得到如图 7所示的解包裹相位图。图中最高点纵坐标为85 rad，最低点纵坐标为-33 rad。根据(9)式计算可得轴向分辨率为11.8 μm。

Figure 7.  The unwrapped phase diagram of inner epidermis of a layer of onion

进一步地，在对一层洋葱表皮成像的基础上，覆盖第2层洋葱内表皮对其成像。在OCT系统中，确保样品臂和参考臂的光程差尽可能小是至关重要的，然而，当使用两层洋葱内表皮作样品成像时，样品的厚度增加，在样品臂中，空气中的光路逐渐被洋葱组织中具有更高折射率的光路所取代。因此，对于这种折射率的变化，需要再次调整样品臂的距离，确保两臂的光程差在光源的一个相干长度内，才能继续得到干涉条纹信息。获取到的两层洋葱内表皮的干涉条纹如图 8所示。其中图 8a~图 8d分别为相位相差202.5 nm波长距离的干涉图。

Figure 8.  Interferogram of the inner epidermis of two layers of onion

计算得出的相位和振幅图如图 9所示。

Figure 9.  a—phase diagram of the inner epidermis of two layers of onion  b—amplitude diagram of the inner epidermis of two layers of onion

图 9b中，A处为细胞，B处为细胞壁。从图 9中可以看出，对两层洋葱内表皮进行成像，依然可以得到清晰的相位图和组织结构图，说明对两层洋葱内表皮的成像效果较为理想。对相位图进行解包运算得到的解包裹相位图，如图 10所示。图中最高点纵坐标为43 rad，最低点纵坐标为-75 rad。

Figure 10.  The unwrapped phase diagram of inner epidermis of two layers of onion

关于灵敏度的测定，以洋葱表皮细胞图像为例，从图 6和图 9中可以看出洋葱表皮中的细胞和细胞壁等构造，可知本系统对组织内部的构造特征有着较为有效的灵敏度，对组织的微观构造能有所表现。

本文中提出了一种基于LED照明的时域全场OCT成像系统，阐述了其扫描的实现原理，开展了系统结构研究、方法推导和实验研究，分别对平面镜和一层、两层洋葱表皮进行了扫描成像，获取了2维层析图，表明了本方式的合理性和准确性。与以往的OCT扫描方式相比，由于使用LED作为入射光，所以减小了实现成本，但对光束采取了均匀准直处理，也能够让LED有着较理想的光束质量，进而保证了扫描图像的质量；同时，采用宽场照明和四步移相的方法，使得扫描的速率和图像的清晰度得到了大幅度增强。总的来说，本文中研究的OCT系统不仅减小了成本，还改善了扫描速率和质量，为今后研究速率更快、成本更小的OCT系统添加了一种实现思路。