光学相干断层扫描成像中散斑噪声的抑制系统和方法
技术领域
本发明涉及低相干光学干涉成像技术领域,特别是一种光学相干断层扫描成像中散斑噪声的抑制系统和方法。
背景技术
光学相干断层扫描成像(optical coherence tomography,OCT)是一种低相干光学干涉成像技术。它可非接触、非侵入性地对人体等生物组织的散射结构进行实时的三维成像。该技术在检测眼底视网膜、消化道系统、呼吸道系统、心血管、牙齿和皮肤等方面的疾病上有很大的潜力。目前OCT技术已经被广泛应用到视网膜疾病诊断、随访观察及治疗效果评价方面。由于需要对更小的病灶区域进行细微观察分析以及对感兴趣病变区域进行定量化评价,因此对OCT成像质量要求变得越来越高。然而OCT作为一种利用相干光的空间相干性和时间相干性通过干涉来实现对生物组织进行成像的技术,必然存在散斑现象。OCT中的散斑作为有用的信号存在时,其携带着部分生物组织的信息,但其同时也是一种噪声存在,作为噪声存在时,这种较强的强度的随机涨落,降低了图像对比度,干扰了对图像细节的分辨,严重影响了成像质量和后期定量分析的准确性。
现今OCT已经发展到第三代,基于扫频光源的OCT(swept-source OCT,SS-OCT),第一、二代分别是时域OCT和谱域OCT。谱域OCT和扫频光源OCT又共称为傅里叶域OCT,扫频光源OCT相比于谱域OCT,它采用扫频光源提供的k时钟作为采样时钟,从而使得测量得到的干涉光谱在波数空间线性均匀分布,因而可以直接对测得的干涉光谱进行离散傅里叶变化来获取深度反射轮廓,大大简化了后期的信号处理,而谱域OCT则需要在计算离散傅里叶变换之前对测量得到的光谱进行重采样来使得光谱在波数空间线性均匀分布。在SS-OCT系统中,散斑噪声主要来源于散射粒子多次的随机散射引起的相位的叠加。
目前,主要有两大类的方法来有效抑制散斑噪声。一种方法是需要对OCT的硬件系统的光路进行相关的改进或扫描模式的改变,从而获得在相同的位置具有低相关散斑模式的多幅断层图像,最后将这些具有低相关散斑模式的断层图像进行平均便能有效的减少散斑对比度。其中低相关的散斑模式可以来源于不同的方法:时间变化、空间变化、频率变化、角度变化和极化状态变化等(例如,见Michael Pircher等人的“Speckle reduction inoptical coherence tomography by frequency compounding(用频域复合方法减少光学相干断层成像中的散斑)”,Journal of Biomedical Optics(生物医学光学期刊),2003;8(3):565-569.)另外一种方法是基于计算机软件处理的方法,主要是通过对得到的断层图像进行相关的图像处理来有效减少散斑对比度。主要包括:小波变换、曲波变换、图像配准和滤波算法等(例如,见Aydogan Ozcan等人的“Speckle reduction in opticalcoherence tomography images using digital filtering(用数字滤波方法来减少光学断层图像中的散斑)”,Optical Society of America(美国光学学会),2007;24(7):1901-10.)
发明内容
本发明要解决的技术问题为:针对基于扫频光源的OCT即SS-OCT,提出一种散斑噪声抑制方法,在不牺牲成像速度和横向分辨率的同时,有效地降低散斑噪声,从而提高OCT断层图像质量。
本发明采取的技术方案具体为:一种光学相干断层扫描成像中散斑噪声的抑制系统,包括扫频光源、样品臂光路、参考臂光路、平衡探测器、数据采集器、第一光分路器和第二光分路器;
扫频光源通过第一光分路器分别连接样品臂光路和参考臂光路;第一光分路器包括样品背向散射光输出端,所述样品背向散射光输出端与参考臂光路的输出端,经第二光分路器连接平衡探测器的输入端,平衡探测器的输出端连接数据采集器;
所述样品臂光路上设有第一准直器、分光镜组、第一反射镜组、X-Y振镜扫描仪;第一准直器连接第一光分路器;由准直器传递出的光束经分光镜组分光,和第一反射镜组反射后形成至少2个相互平行的光束,所述至少2个相互平行的光束经X-Y振镜扫描仪传递至样品上。
本发明的数据采集器可利用现有图像采集技术。在应用时,样品臂光路中入射光和背向散射光可以通过两平行光束传递路径的两个光路进行传递,这样便形成了三种不同的光程:1)最短光程:入射光和背向散射返回光都通过两光路中较短路径的一个光路;2)中间光程:入射光和背向散射光分别通过两光路中的一个光路;3)最长光程:入射光和背向散射光都通过两光路中较长路径的一个光路。三个不同光路对于样品产生的图像是不同角度,因而三个图像具有不相关的散斑模式;同时三个不同光路具有不等的光程,因而在一次采集的二维断层图像(即一个B-scan)中即表现为三个不同深度位置的图像,通过对这样三个图像进行处理即可得到改善散斑对比度的图像。
进一步的,本发明样品臂光路上还设有透镜组,透镜组设置于X-Y振镜扫描仪与样品之间的光路上。透镜组可采用两个焦距不同的透镜,以调整两束相互平行的光束之间的距离和光束直径的大小,便于两光束同时入射到样品上。
更进一步的,本发明的参考臂光路上依次设有第二准直器、第二反射镜组、第三准直器和光纤耦合器;光纤耦合器的输出端即参考臂光路的输出端。所述第二反射镜组包含2个法线相互垂直的反射镜,参考臂光束经第二准直器聚焦后输出为平行光束,然后经第二反射镜组反射后形成与前述平行光束相互平行的光束,再经第三准直器传输至光纤耦合器与样品臂光路的背向散射光发生干涉,干涉信号由平衡探测器测量获得。
优选的,本发明的样品臂光路中,分光镜组包括第一分光镜和第二分光镜,第一反射镜组包含两个反射镜;第一分光镜将第一准直器传递来的光束分成两个光束,两反射镜分别设置于上述两个光束的传递路径上,并分别将相应光束反射至第二分光镜;两反射镜反射出的两个光束经第二分光镜后,至少有两个光束相互平行射出;X-Y振镜扫描仪设置于上述两个相互平行光束的传递路径上,以将两个相互平行的光束传递至样本上。两平行光束之间的距离可通过平移反射镜来实现,以保证两束光都能够通过成像物镜聚焦到样品上。在此基础上,平移两个反射镜以改变两平行光束之间的距离,可调整最终获得的多幅样品图像之间的深度方向分离度。
样品臂光路中的第一分光镜和第二分光镜的分光比皆为50/50;两分光镜的分光面以及两反射镜的反射面四者之间相互平行,且与准直器输出的光束之间呈45度夹角。
优选的,第一光分路器和第二光分路器的分光比皆为50/50。有利于提高光纤的耦合效率,更便于根据所获取的图像数据进行后续的处理。
本发明还公开基于上述系统的散斑噪声抑制方法,数据采集器至少获得样品的三幅基于不同散斑模式的图像,对所述三幅图像的处理包括以下步骤:
步骤一,对三幅图像进行分割,分成三个相同尺寸的子图;
步骤二,选择三幅子图中的其中一幅作为参考子图,将另外两幅子图分别向参考子图进行配准,得到配准后的三幅图像;
步骤三,对配准后的三幅图像进行两两平均或三者平均,得到最终散斑已改善的样品图像。
本发明的有益效果为:通过对样品臂光路进行改造,利用两个分光镜和一对反射镜得到两个平行的光路,使得入射光和返回光的光路能够具有不同的选择(如两平行光路为光路1和光路2,入射光和背向散射返回光即可选择:光路1入射原路返回,光路2入射原路返回,光路1入射光路2返回或光路2入射光路1返回),从而形成三个不同光程差的子图,这三个子图具有不同的散斑模式,对这三个子图进行平均处理后即可有效降低散斑噪声。本发明在硬件方面对传统硬件系统的改动很小,并且在一次采集的二维断层图像(即一个B-scan)中可以获得三个子图,提高了成像速度,这三个子图由于不同的角度关系具有不相关的散斑模式,三个子图取平均后能够有效地降低散斑噪声。即本发明在不牺牲成像速度和横向分辨率的同时,能够有效地降低散斑噪声,从而提高OCT断层图像质量。
附图说明
图1所示为本发明系统结构示意图;
图2所示为样品臂光路结构示意图;
图3所示为实施例所得样品图像处理过程示意图;
图4所示为双光路情况下,三个子图的相关性系数随光束分开距离的变化示意图;
图5所示为双光路情况下,三个子图平均处理后相对于原图散斑对比度的提高倍数示意图;
图6所示为单光路情况下,三个子图的相关性系数随光束分开距离的变化趋势示意图。
具体实施方式
以下结合附图和具体实施例进一步描述。
结合图1和图2,本发明的光学相干断层扫描成像中散斑噪声的抑制系统,包括扫频光源SS、样品臂光路Sample ARM、参考臂光路、平衡探测器BD、数据采集器、第一光分路器FC1和第二光分路器FC2;
如图1,扫频光源通过第一光分路器分别连接样品臂光路和参考臂光路;第一光分路器包括样品背向散射光输出端,所述样品背向散射光输出端与参考臂光路的输出端,经第二光分路器连接平衡探测器的输入端,平衡探测器的输出端连接数据采集器;
如图2,所述样品臂光路上设有第一准直器CL3、分光镜组、第一反射镜组、X-Y振镜扫描仪GS;第一准直器CL3连接第一光分路器FC1;由准直器CL3传递出的光束经分光镜组分光,和第一反射镜组反射后形成至少2个相互平行的光束,所述至少2个相互平行的光束经X-Y振镜扫描仪传递至样品上。
本发明的数据采集器可利用现有图像采集技术。在应用时,样品臂光路中入射光和入射样品后返回的背向散射光可以通过两平行光束传递路径的两个光路进行传递,这样便形成了三种不同的光程:1)最短光程:入射光和背向散射返回光都通过两光路中较短路径的一个光路;2)中间光程:入射光和背向散射光分别通过两光路中的一个光路;3)最长光程:入射光和背向散射光都通过两光路中较长路径的一个光路。三个不同光路的光程对于样品产生的图像是不同角度图像的复合,由于不同的角度关系即具有不相关的散斑模式,在一次采集的二维断层图像(即一个B-scan)中即表现为三个不同深度方向的图像,通过对这样三个图像进行处理即可得到改善散斑对比度的图像。
实施例
参考图1,本实施例以眼球模型作为样品。
参考臂光路上依次设有第二准直器CL1、第二反射镜组、第三准直器和光纤耦合器FC3;光纤耦合器的输出端即参考臂光路的输出端。所述第二反射镜组包含2个法线相互垂直的反射镜M1和M2,参考臂光束经第二准直器CL1聚焦后输出为平行光束,然后经第二反射镜组反射后形成与前述平行光束相互平行的光束,再经第三准直器CL2传输至光纤耦合器与样品臂光路的背向散射光发生干涉,干涉信号由平衡探测器测量获得。
样品臂光路中,还设有透镜组,透镜组采用两个焦距不同的透镜L,以调整两束相互平行的光束之间的距离和光束直径的大小,使得两光束同时入射到样品眼球模型上。透镜组设置于X-Y振镜扫描仪GS与样品之间的光路上。分光镜组包括第一分光镜BS1和第二分光镜BS2,第一反射镜组包含两个反射镜M3和M4;第一分光镜将第一准直器CL3传递来的光束分成两个光束,两反射镜分别设置于上述两个光束的传递路径上,并分别将相应光束反射至第二分光镜;两反射镜反射出的两个光束经第二分光镜后,至少有两个光束相互平行射出;X-Y振镜扫描仪设置于上述两个相互平行光束的传递路径上,以将两个相互平行的光束传递至透镜组,进而传递至样品上。
样品臂光路中的第一分光镜和第二分光镜的分光比皆为50/50;两分光镜的分光面以及两反射镜的反射面四者之间相互平行,且与准直器输出的光束之间呈45度夹角。
第一光分路器和第二光分路器的分光比皆为50/50,有利于提高光纤的耦合效率,更便于根据所获取的图像数据进行后续的处理。
扫频光源的光源中心波长为1060纳米,波长范围从985纳米到1095纳米。激光器平均输出功率为20毫瓦,扫频速度是100千赫兹。从扫频激光源中出来的光由FC1分出两束光,分别进入到样品臂(50%)和参考臂(50%),然后经由两臂后的光通过50/50的FC2发生干涉,干涉信号由平衡探测器BD探测得到,并将干涉光信号转变为电信号传递给数据采集器。数据采集器可采用14位的数据采集卡,并由扫频光源提供外部k时钟驱动,通过k时钟把光谱采样到线性波数空间(K空间)并做傅立叶变换计算。数据采集程序可建立在LabVIEW平台上,用来采集数据和控制X-Y振镜扫描仪,并且用LabVIEW进行数据处理可以实时展现干涉光谱图像。上述图像采集及显示处理为现有技术。
样品臂光路中,FC1分出的其中一束光首先经过准直器,然后通过50/50分光镜分成相等的两部分光(反射部分:实线表示;透射部分:虚线表示)。反射和透射的光分别通过反射镜传递给下一个50/50分光镜,然后形成两束平行的光束传递给X-Y振镜扫描仪,最后通过一对聚焦长度75毫米和35毫米的透镜缩小光束后入射到样品上。
参考图2,样品臂光路形成的两平行光束可通过平移反射镜来改变两光束之间的分开距离,从而使得三幅图像深度方向分开更大。同时入射光和入射样品后返回的背向散射光可以通过以上两个光路(实线和虚线),便形成了三种不同的光程:1)最短光程:入射光和背向散射返回光都通过虚线光路;2)中间光程:入射光通过实线光路背向散射返回光通过虚线光路或者入射光通过虚线光路背向散射返回光通过实线光路;3)最长光程:入射光和背向散射返回光都通过实线光路。三个不同的光程对于样品产生的图像是不用角度图像的复合,在一次采集的二维断层图像(即一个B-scan)中表现为三个不同深度方向的图像,如图3(a)所示。
参考图3,本发明基于上述系统的散斑噪声抑制方法,数据采集器至少获得样品的三幅基于不同散斑模式的图像,对所述三幅图像的处理包括以下步骤:
步骤一,对三幅图像进行分割,分成三个相同尺寸的子图;
步骤二,选择三幅子图中的其中一幅作为参考子图,将另外两幅子图分别向参考子图进行配准,得到配准后的三幅图像;如图3(a);
步骤三,对配准后的三幅图像进行两两平均或三者平均,得到最终散斑已改善的样品图像,如图3(b)(c)(d)(e)。图3(a)中的方框区域为可自行定义的感兴趣区域,用来计算散斑对比度。
实验验证
通过选用眼球模型作为实验样品,所采集的B-scan的大小为512*2000像素,它包含3个不同散斑模式的眼球模型的OCT图像,如图3(a)所示。通过将不同散斑模式的图像两两平均或三个平均后得到的图像,可以看出散斑噪声被有效抑制。
进一步的,为了了解光束分开距离的大小对于我们方法的影响。我们首先分开两束光直到三幅图像填满整个B-scan,然后以25微米的步长同时移动两个反射镜来不断减小光束分开的距离。评估每一个光束分开距离的情况对于散斑的抑制效果。三个子图两两之间的相关性系数随这光束分开距离的变化趋势,如图4所示,光束分开距离在700-1150微米时,彼此的相关性系数表现出一个轻微的波动。同时我们计算了1+2平均,1+3平均,2+3平均和1+2+3平均的散斑对比度(C),散斑对比度等于感兴趣区域的均值(<Ι>)与感兴趣区域标准差(σI)的比值,即C=<I>/σI。我们选择如图3(a)中的方框区域作为感兴趣区域来计算散斑对比度,它避开了底部最黑的镜面反射的干扰。
为了对比散斑的抑制效果,我们计算了平均后图像的散斑对比度对于原始图像散斑对比度提高的倍数如图5所示,子图两两平均后的散斑对比度提高倍数在1.28倍左右轻微波动。这正是因为它们之间的相关性系数的轻微波动的关系。也可以发现散斑对比度提高倍数与相关性系数成反比例关系,即相关性系数越低平均后散斑对比度提高倍数越高,抑制散斑效果越好。
除此之外,我们使用单个光路(另一个光路被挡住)来研究更大范围的光束分离对于彼此相关性系数的影响。如图6所示,当光束分离小于600微米时,相关性系数呈现出明显的下降趋势;当光束分离大于600微米时,相关性系数几乎保持不变,这与双光路情况下得到的结论一致。
用我们的方法在三个子图平均后能够使得散斑对比度提高1.46倍,如图5所示,能够在提高成像速度的同时有效降低散斑噪声。
综上所述,本发明通过对样品臂光路进行改造,得到样品基于不同散斑模式的三个子图,再对这三个子图进行处理,从而有效降低散斑噪声。在硬件方面对传统硬件系统的改动很小,在不牺牲成像速度和横向分辨率的同时,提高了OCT断层图像质量。