散射成像技术的研究进展 下载: 6573次特邀综述
1 引言
生活中,人们观察物体、辨别物体时总是遵循着“光沿直线传播”“所见即所得”的传统光学成像规律,常见的视觉成像系统如人眼、照相机,以及光学成像工具如放大镜、显微镜和望远镜等都以此为物理基础。传统的光学成像主要通过提取弹道光(或抑制散射光)的方式解决透过无散射(或者弱散射)介质的成像问题。当光透过生物组织、烟尘和云雾等强散射介质时,传统的光学成像规律就会失效。其主要原因是光在强散射介质中传播时,介质中大小为波长量级的粒子对入射光波产生的散射作用使原本有序的波前相位发生了严重畸变,出射光场变得随机且紊乱,最终在观测面上只能接收到散斑图案,难以实现对目标的观测或成像。因此,散射效应成为制约透过散射介质成像技术发展的瓶颈和关键问题。
近年来,计算机技术的飞速发展、介观物理研究的深入、计算成像思想的完善和图像处理技术的发展,促进了以物理机制为基础的计算光学成像技术的发展。计算光学成像技术作为新型的成像手段,不仅推动了传统成像技术的发展,而且在解决散射成像方面表现出了得天独厚的优势。经过各国科学家们的不懈努力,计算成像理论及相关实验技术迅速发展,取得了许多突破性的研究成果。在弹道光提取方面,自适应光学成像技术[1]、门选通技术[2-3]、光学相干层析技术[4]、共聚焦显微技术[5]、多光子显微技术[6-7]、光声显微技术[8-9]、复合荧光分子层析技术[10]、多光谱光声层析技术[11-12]等光学成像技术的发展及应用,解决了天文成像、水下探测和生物成像等领域的关键问题。在散射光利用方面,空间光调制器(SLM)、数字微镜器件(DMD)、微机电系统调制器(MEMS)等数字波前整形器件的出现,促进了计算成像技术与散射成像技术的有机结合,涌现出了许多新型的光学成像技术,如光学相位共轭[13]、波前反馈调节[14-17]、光学传输矩阵[18-20]等,极大地促进了散射成像技术在显微成像领域的应用。与此同时,随着对光学记忆效应[21]的深入研究,散斑相关成像技术[22-23]的提出为透过散射介质成像打开了崭新的局面。随着计算光学成像的进一步发展,基于点扩展函数工程[24]和深度学习[25]的方法,也成为了解决透过散射介质成像的重要途径。
纵观散射成像的发展历程,散射成像技术已从早期的基础理论研究发展到了实验室中的模型验证研究,继而又发展到了透过散射介质成像的应用研究。早期的散射成像技术注重于克服散射或抑制散射,通过弹道光与散射光的分离,最终获取有效的目标信息。现阶段的散射成像技术则侧重于散射光的利用,充分挖掘散射光的特性,实现从不可探测到可探测的质的飞跃。值得注意的是,散射成像技术不仅在显微成像和超分辨成像方面有着广泛的应用,而且也将在全息成像、光纤成像和光通信等领域扮演着重要的角色。
目前,透过散射介质的成像方法主要有波前整形和基于光学记忆效应的散射成像技术。波前整形技术包含光学相位共轭、基于反馈优化的波前整形和光学传输矩阵技术三部分,该技术主要研究光波在介质中的传播规律及特性,为散射效应的利用奠定了基础。基于光学记忆效应的散射成像技术包括散斑相关成像技术和点扩展函数工程成像技术两部分,该技术的核心在于对散斑能量的利用,利用散斑的分布特点实现透过散射介质成像及相关工作。本文分别对以上方法在散射成像方面的成果及进展进行归纳和总结,并对其未来的发展趋势进行了展望。
2 基于波前整形的散射成像技术
光在散射介质中传播时会受到散射效应的影响,在像面上形成系列散斑。如何定量或定性地描述散射介质在光传播过程中的影响已成为利用光散射效应的关键问题。虽然光在多重散射介质中的传播具有很高的随机性,但当散射介质处于稳定状态时,光在散射介质中的传播就具有确定性。散射介质的特性可以与多模光纤的特性进行类比,其输出光场可以看作是多种模式的耦合与叠加。为了精确地描述散射介质在光传播过程中的作用,光学传输矩阵的思想被提出,从而可以有效地将入射光场与出射光场联系起来。近年来,随着光散射理论与实验技术的飞速发展,研究人员基于光学相位共轭、反馈优化的波前整形和光学传输矩阵等技术,实现了透过散射介质的聚焦或者成像[13-20]。在未能获得完善的光学传输矩阵的条件下,为实现透过散射介质聚焦或成像,通常采用光学相位共轭技术或者基于反馈优化的波前整形技术;而在测得完备的光学矩阵之后,利用光学传输矩阵技术往往能够有效地实现对出射光场的控制。相关研究表明,波前整形技术在生物医学成像、超分辨成像和光通信等方面有着广阔的应用前景。
2.1 光学相位共轭
光学相位共轭技术是时间反演技术在光学领域的应用,最早的光学相位共轭技术通过在照相板上记录全息图来实现[26-29]。本质上,光学相位共轭技术利用的是光传播的可逆性,获得透过介质后的光场分布后,反向输入透射光场的相位共轭波形,就可以重建原始的入射光场。光学相位共轭的实现步骤可以分为两步:第一步,光场信息的记录;第二步,相位共轭光的生成。按照相位共轭光产生方式的不同,光学相位共轭技术可以分为非线性光学相位共轭技术和数字光学相位共轭技术两类[30]。前者可以使用数字全息或定量相位成像技术来实现,后者则可以通过空间光调制器实现。在实验中,在多重散射材料的情况下即便是只获得散射介质某一侧的部分光场信息,也可以通过光学相位共轭实现初始波前的重建。
根据非线性过程的差异,非线性光学相位共轭技术可以分为三波混频相位共轭技术[31]、四波混频相位共轭技术[32-33]、受激布里渊散射相位共轭技术[34]和光折变晶体相位共轭技术[35]。总体而言,非线性光学相位共轭技术实施起来比较复杂,通常需要非线性晶体、特定波长和强激光光源。虽然实施起来比较复杂,但是非线性相位共轭技术自提出以来就多用于透过复杂介质的聚焦。光折变晶体是光学相位共轭技术中常用的一种材料,虽然其调节速度较慢,但可透过厚生物组织实现聚焦成像。随着材料技术的飞速发展,许多新型的光学共轭材料不断地被研发出来,其调节速度可与SLM相媲美。许多增益介质能够提供快速的光学共轭调节,但受到其物理效应的限制,仅适用于窄谱光源。三波混频相位共轭技术具有速度快、频带宽的特点,但其有效角度较小。与数字光学相位共轭技术相比,非线性光学相位共轭技术在模式耦合效率方面具有较大优势[29],因此,非线性光学相位共轭技术在生物医学方面仍具有很大潜力。
SLM和DMD等元件的出现,使得数字光学相位共轭调制变成可能,数字光学相位共轭技术的工作原理如
数字光学相位共轭技术虽然在调制效率方面具有一定劣势,但是在透过复杂介质成像方面却有着得天独厚的优势。数字光学相位共轭技术利用计算机记录其输出光场分布,利用调制器生成共轭光,可以对无数个输入光场进行重构。相比较而言,非线性光学相位共轭技术不具有这一特点。如何提高数字光学相位共轭技术的效率问题,将决定数字光学相位共轭技术在未来应用中的地位。
2.2 基于反馈优化波前整形的散射介质成像技术
基于反馈优化的波前整形技术利用优化算法(非线性或线性),通过迭代获取目标光场对应的最优波前,从而实现透过散射介质聚焦或成像。本质上,基于反馈优化的波前整形技术将散射介质对光场的调制过程看作“黑箱”,通过迭代算法获取相应的波前,进而实现对输出光场的模式及不同模式之间耦合的控制。
2007年,Mosk课题组[14]利用SLM对入射到随机散射介质中的光波进行波前相位调制,采用反馈控制算法对空间光调制器的SLM像素进行逐个优化,通过不断迭代的方式获得了最优波前,利用所得的最优波前幅值或相位可以适当补偿由介质散射引起的波前畸变,最终得到了亮度高于调制前散斑1000倍的聚焦光斑,远远优于光学透镜的聚焦效果。波前整形示意图和实验结果分别如
图 2. 波前整形示意图[14]。(a)波前整形前;(b)波前整形后
Fig. 2. Schematics of wavefront optimization [14]. (a) Before wavefront optimization;(b) after wavefront optimization
图 3. 波前整形聚焦结果[14]。(a)聚焦前散斑;(b)单点聚焦;(c)多点聚焦;(d)优化的波前相位
Fig. 3. Focusing results after wavefront optimization[14]. (a) Speckle before focusing; (b) focusing on single point; (c) focusing on multiple points; (d) optimized wavefront phase
2010年,Mosk课题组[15]通过在传统透镜光学系统后面放置一块厚度为6 μm的散射介质(如
图 4. 超衍射极限聚焦示意图[15]。(a)传统透镜聚焦光学系统;(b)随机散射介质聚焦光学系统
Fig. 4. Schematics of focusing beyond diffraction-limit [15]. (a) Focusing optical system with conventional lens; (b) focusing optical system with random scattering media
图 5. 超衍射极限聚焦实验结果[15]。(a)传统透镜光学系统聚焦光斑;(b)随机散射介质光学系统调制前的散斑; (c)超衍射极限聚焦光斑;(d)优化的波前相位
Fig. 5. Experimental results of focusing beyond diffraction-limit[15]. (a) Focal spot of focusing optical system with conventional lens; (b) speckle of focusing optical system with random scattering media before optical modulation; (c) focal spot beyond diffraction-limit; (d) optimized wavefront phase
2011年,Mosk课题组[39]在利用随机散射介质实现超衍射极限聚焦工作的基础上,提出了一种基于SLM的随机散射透镜超衍射极限扫描显微成像方法。该课题组利用反馈控制调节聚焦的方法将光源聚焦到显微镜的物平面上,调节SLM并扫描整个物平面,实现了对金纳米粒子的观测,分辨率可达97 nm。该方法不需要使用荧光分子对样本进行标记,可以实现对样本的无损观测,实验结果如
图 6. 利用随机散射透镜实现的超衍射极限成像[39]。(a)传统显微镜成像;(b)超衍射极限成像;(c)图6 (a)左边第一颗粒子的中心切线与图6 (b)左边第一颗粒子的中心切线的对比
Fig. 6. Imaging beyond diffraction-limit using random scattering lens[39]. (a) Imaging using conventional microscope; (b) imaging beyond diffraction-limit; (c) comparison of center tangents of first spots on left of Figs. 6 (a) and (b)
时效性是基于反馈优化的波前整形技术的关键问题。为解决这一问题,Cui[40]提出了一种空间频率域并行波前优化的方法,该方法可以提升原有反馈控制算法的空间聚焦效率;Conkey等[16]采用全局优化的思想代替原有的逐点优化,通过遗传算法实现了最终的聚焦,得到了更好的聚焦效果;随后,Blochet等[17]将基于MEMS的SLM与现场可编程门阵列(FPGA)技术结合,实现了4 kHz左右的调制速度,将该技术应用于波前整形可有效提升成像速度。
值得注意的是,2012年,Katz等[41]基于反馈优化的波前整形技术,采用非相干光源实现了透过散射介质的实时成像,此项工作极大地推动了波前整形技术在实际应用上的进程。同时,基于反馈优化的波前整形技术在透过多模光纤的光学精密控制和成像方面也有着重要意义[42-45]。此外,考虑到实际应用中散射介质的时变特性,如何实现透过动态散射介质的快速聚焦或成像是未来研究中的重要课题。
2.3 光学传输矩阵的散射介质成像技术
在Mosk等利用波前调制技术验证了可见光透过散射介质后依然能够实现聚焦之后,基于波前调制技术的相关研究引起了世界各国科研人员的广泛关注。测量传输矩阵是一种新的波前调制技术,该方法的核心思想是利用一个复杂的矩阵将入射光场与出射光场联系起来,通过测量传输矩阵并结合相位共轭技术,能够在任意位置、任意时刻实现聚焦和成像[46]。
假设光源是线性极化的,则对于光波在任何介质中的传播过程,都可以用格林函数进行描述。离散化的格林函数可以表征介质对入射光波的作用,离散化的格林函数就是上述提到的光学传输矩阵,它表示m个输出单元与n个输入单元的光场信息(振幅和相位)之间的相互关系,可以表示为
式中:
在测得光学传输矩阵K后,假设目标输出光场为
式中:
由以上分析可知,一旦测得散射介质的光学传输矩阵,就可以采用相位共轭的方法实现聚焦。
2010年,法国卡斯特勒-布罗塞尔实验室(Laboratoire Kastler-Brossel)的Popoff等[18]利用四步相移干涉法测得了ZnO散射介质的光学传输矩阵,
在测得ZnO介质的传输矩阵后,Popoff等[18]利用Tikhonov正则化图像恢复算法对散斑进行了重建,实验结果如
基于光学传输矩阵的散射成像方法的优点在于只要测量出成像系统的光学传输矩阵,便可以从任意目标所成的散斑中迅速恢复出待测目标。但是,就现阶段的研究来看,该方法所需系统较复杂,对系统稳定性的要求非常高,任何改变都有可能导致无法重建目标,目前的研究水平还无法做到对实际物体成像。
图 7. 四步相移测量传输矩阵的实验原理图[18]
Fig. 7. Experimental schematic of measuring optical transmission matrix based on four-step phase shift[18]
图 8. 采用相位共轭法得到的聚焦结果[18]。(a)聚焦前散斑;(b)单点聚焦结果;(c)多点聚焦结果
Fig. 8. Experimental results of focusing via phase conjugation[18]. (a) Speckle before focusing; (b) result of focusing on single point; (c) result of focusing on multiple points
虽然光学传输矩阵方法已在实验中得到了验证,但是在实际应用中仍有许多不足之处。随着对光学传输矩阵方法研究的深入,许多光学传输矩阵测量方法被提出,如:2011年Choi等[47]提出的频域传输矩阵测量方法、2015年Drémeau等[48]提出的基于相位恢复算法的测量方法、2015年Yoon等[49]提出的基于波前相位调制的测量方法。不同的测量方法具有不同的特点,如:频域测量方法无法对倏逝波矩阵进行测量,基于相位恢复算法的测量方法的稳定性差,而基于波前调制的测量方法的精确程度受限于调制光的精准程度。2012年,Tripathi等[50]提出了一种测量矢量光学传输矩阵的方法,利用该方法可以获得指定偏振态的聚焦光波;2015年,Andreoli等[51]提出了一种基于测量多光谱传输矩阵的方法,该方法可以实现宽光谱聚焦;2016年,Mounaix等[52]利用多光谱传输矩阵方法实现了控制超短脉冲激光通过散射介质的空时聚焦;在多光谱传输矩阵研究的基础上,2019年Dong等[53]利用多路复用相位反演的光谱方法实现了透过散射介质聚焦。多物理量探测有助于信号的探测和识别,增加光学传输矩阵的测量维度(光谱和偏振)等将在未来的实际应用中起到至关重要的作用。另外,如何保证光学传输矩阵测量的实时性是未来研究的重要方向。
3 基于光学记忆效应的散射成像技术
基于光学记忆效应的散射成像技术是透过散射介质成像的重要组成部分,它可以利用散斑,通过自相关运算,获取目标信息的傅里叶振幅,进而结合有效的相位恢复算法实现目标的重建。与波前整形技术相比,基于光学记忆效应的散射成像技术具有非入侵的特点,且对于光源、介质和系统的要求较低。随着对光学记忆效应研究的深入,利用散射介质的退相关特性,可实现透过散射的光谱成像[54-55]和三维成像[56-57],这将对未来的新型成像系统具有重要意义。
3.1 光学记忆效应
在一定的入射角度范围内,当改变光源入射方向时,经过散射介质在像平面上得到的散斑形状特征保持不变,但整体发生了平移,这一现象称为光学记忆效应。
3.2 基于散斑相关的散射成像技术
图 10. 基于光学记忆效应的非侵入式成像方法示意图[22]
Fig. 10. Schematic of non-invasive imaging based on optical memory effect [22]
2012年,Bertolotti等[22]在Nature杂志上首次提出了利用散斑相关技术实现非侵入式的透过散射介质成像,实验系统如
即散斑可以看成是原目标O(r)与上述成像系统点扩展函数S(r)卷积的结果,r表示空间位置,目标信息已经被“编码”至散斑中,d为散射介质与探测器之间的距离。对散斑作自相关运算,并运用卷积定理可以得到
式中:<·>为均值运算;★为相关运算;*为卷积运算。
依据介观物理理论[58-59],当介质的有效厚度远小于Anderson localization长度,并且散射足够充分时,<S★S>可以表示为类δ函数[22],即目标的自相关信息可近似由散斑自相关得到。由Wiener-Khinchin定理可知,从目标的自相关信息中可以得到目标在频率域的振幅信息,再结合相位恢复算法并采用迭代的方法就可以得到目标在频率域的相位信息,从而实现目标的重建,实验结果如
2014年,Katz等[23]在Bertolotti研究的基础上提出了一种基于单帧散斑自相关的成像方法,成像模型和实验结果分别如
图 11. 基于光学记忆效应的散射成像方法的实验结果[22]。(a)散斑;(b)散斑自相关;(c)原目标;(d)重建目标
Fig. 11. Experimental results of scattering imaging based on optical memory effect[22]. (a) Speckle; (b) autocorrelation of speckle; (c) original object; (d) reconstructed object
图 12. 基于单帧散斑自相关成像的示意图[23]。(a)成像模型示意图;(b)散斑;(c)散斑自相关;(d)重建目标
Fig. 12. Schematic of single-frame imaging based on speckle autocorrelation [23]. (a) Schematic of imaging model; (b) speckle; (c) speckle autocorrelation; (d) reconstructed objects
针对Fienup型相位恢复算法的缺点,Wu等[62]将天文观测中的三阶相关相位恢复技术(又称双谱技术[63])引入到散射成像的相位恢复过程中,用以解决由相位恢复算法所导致的重建结果不确定的问题。采用双谱分析得到的重建结果如
三阶相关相位恢复技术与Fienup型相位恢复算法的散斑相关成像相比具有以下优势:1)能够精确地恢复出目标频率域的相位信息;2)不需要大量的迭代与重复估计;3) 具有更好的抗噪性;4)保留了目标的方向信息。
除了上文所提到的静态散射介质,透过动态散射介质成像是透过散射介质[64-65]成像领域不可或缺的一部分。2016年,美国马里兰大学的Scarcelli课题组[64]提出了一种基于浴帘效应的散射成像方法。该方法通过对多帧散斑作傅里叶逆变换并叠加取平均值,近似求得了原目标场的自相关,然后结合相位恢复算法得到了目标的频率域相位,进而恢复出了目标信息,如
基于散斑相关的散射成像方法具有系统简单且能够实现非侵入式成像的特点,但其缺点也很明显,如:成像的视场受记忆效应范围的限制,与平均传输自由程成反比。鉴于此,许多扩展光学记忆效应的方法和技术应运而生,如:利用参考目标[66]或先验知识[67]实现超记忆效应范围成像的技术,通过光场估计[68]或散斑估计[69]的方法实现透过散射介质成像的景深扩展。除了以上所描述的基于散斑相关的透过散射介质成像技术外,将散斑相关成像技术与光谱编码技术[24]、压缩感知技术[70-72]、双目视觉技术[73-74]等相结合,也可实现透过散射介质的彩色成像、光谱成像和三维成像等。
图 13. 基于单帧散斑自相关成像的实验结果[23]。(a)成像装置;(b)散斑;(c)~(g)第一列为散斑自相关,第二列为重建目标,第三列为原目标
Fig. 13. Experimental results of single-frame imaging based on speckle autocorrelation[23]. (a) Imaging setup; (b) speckle; (c)-(g) first column shows speckle autocorrelation, second column shows reconstructed objects, and third column shows original objects
图 14. 采用双谱分析得到的重建结果[62]。(a)散斑;(b)傅里叶振幅;(c)傅里叶相位;(d)重建目标;(e)原目标
Fig. 14. Reconstructed results obtained by bi-spectral analysis[62]. (a) Speckles; (b) Fourier amplitudes; (c) Fourier phases; (d) reconstructed objects; (e) original objects
3.3 基于点扩展函数工程的透过散射介质成像技术
根据光学效应范围,透过散射介质成像系统的数学模型如(4)式所示,分析其数学模型不难发现,在传统成像方面已成熟应用的去卷积技术有望实现透过散射介质成像。2016年,Edrei等[75]在测量散射系统点扩展函数的基础上,采用Lucy-Richardson去卷积迭代非线性复原方法[76-77],实现了透过散射介质成像,结果如
图 15. 基于浴帘效应的散射成像结果[64]。(a)原目标;(b)薄散射体远离目标;(c)薄散射体紧贴目标;(d)基于浴帘效应的散射成像系统原理图;(e)目标重建过程
Fig. 15. Experimental results of scattering imaging based on shower-curtain effect[64]. (a) Original object; (b) object is far away from thin scatter; (c) object is close to thin scatter; (d) principle of scattering imaging system based on shower-curtain effect; (e) process of object reconstruction
图 16. 基于去卷积的散射成像结果[75]。(a)实验系统示意图;(b)原目标;(c)散斑;(d)系统点扩展函数;(e)重建结果
Fig. 16. Experimental results of scattering imaging based on deconvolution[75]. (a) Schematic of experimental setup; (b) original object; (c) speckle; (d) point spread function of imaging system; (e) reconstructed result
2017年,西安电子科技大学的邵晓鹏课题组[78]借鉴天文成像所应用的相位多样性技术,通过获取多帧散斑,估算得到了散射光学系统的点扩展函数,进而通过去卷积技术实现了透过散射介质成像。基于相位多样性的散斑成像原理、散斑重建及系统点扩展函数估计结果如
除了上文提到的利用去卷积或估计系统点扩展函数实现透过散射介质成像外,利用散射介质点扩展函数的光谱维和空域正交性,在光谱维或空域多次测量散射介质的点扩展函数,将解卷积得到的物体像进行数据重构,也可以实现透过散射介质的光谱成像和三维成像。例如:2017年,Sahoo等[24]分别测量获得了散射系统在不同波长下的点扩展函数,利用散射介质的光谱退相关特性,将去卷积过程作为光谱滤波的手段,实现了单帧透过散射介质多光谱成像;2018年,Antipa等[79]利用空域中点扩展函数的正交特性,分别测量了物空间不同点源目标所产生的点扩展函数,通过3D-TV[71]去卷积方式,实现了单帧透过散射介质的三维成像,重构过程具有稳定性高、抗噪性能较好的特点。
图 18. 基于相位多样性的散射成像实验结果[78]。(a)原目标;(b)多样性散斑;(c)第一列为估计的随机相位,第二列为估计的局部点扩展函数,第三列为重建目标
Fig. 18. Experimental results of scattering imaging based on phase diversity[78]. (a) Original object; (b) speckle with diversity; (c) first column shows estimated random phases, second column shows estimated local point spread functions, and third column shows reconstructed objects
毫无疑问,基于点扩展函数工程的方法与其他透过散射介质成像方法相比具有较大优势:与波前整形技术相比,无需复杂的波前反馈调制及迭代,仅通过单帧散斑便可实现透过散射介质成像;与散斑相关成像技术相比,避免了相位恢复算法带来的不确定性(如位置信息和方向信息的不确定性),在光谱成像和三维成像方面有着巨大的潜力。然而,基于点扩展函数工程的方法也有不足之处:1)需要获取系统的点扩展函数,在实际应用中往往需要对点扩展进行测量,在实际应用中的难度较大;2)目标重建效果受重建算法的影响较大,合适的去卷积算法对重建效果具有重要作用;3)重建效果受限于光学系统的稳定性,任何器件的移动都可能导致重建结果失效。总而言之,基于点扩展函数工程的方法对于透过散射介质的三维成像和光谱成像具有重要意义。
随着散射成像技术的进一步发展,越来越多的研究表明,散射介质也可以作为一种有益信息被加以利用。一方面,在光学加密解密技术研究中,散射介质可以作为一种“编码变换”实现加密过程,而重建的算法则可以被视为“密钥”进行解密[80-81];另一方面,利用散射光的多方向性可以克服光路中的障碍物而成像,实现“散射透镜”的作用。除以上所介绍的散射成像技术,研究人员在透过散射介质实现景深成像、三维成像及全息成像等方面[82-87]同样取得了许多具有代表性的研究成果,这些成像技术极大地拓展了散射成像技术的应用范围,使得散射成像技术与实际应用的联系更加紧密。
4 结束语
本文围绕当前散射成像技术的研究现状展开叙述,对波前整形和基于光学记忆效应的散射成像技术进行介绍,并对散射成像技术在生物医学领域的应用前景进行了展望。波前整形技术(光学相位共轭技术、基于反馈调节的波前整形技术和光学传输矩阵技术)的关键在于研究光在散射介质中的传播特性,利用数学的形式对散射介质的特性进行定量或定性表示。研究结果表明,利用波前整形技术可以实现透过散射介质聚焦和成像,这在生物医学利用上具有重要意义。在基于光学记忆效应的散射成像技术中,散斑相关成像通过自相关计算恢复了目标的傅里叶振幅,利用相位恢复算法恢复了傅里叶相位,从而实现了非侵入式透过散射介质成像;基于点扩展函数工程的透过散射成像技术利用点扩展函数的退相关特性和随机编码特性,实现了光谱成像和三维成像。总体来看,光学记忆效应、传输矩阵特征值的双模分布、时间反演法等物理概念,均是将电子学、声学的研究中的物理概念迁移到光学研究中,并且一一得到了验证,这为科研人员研究光散射问题指明了一个方向,即可以用某些研究电子散射问题的方法来解决光散射问题,这一事实同时也从另一角度证明了光现象和电磁现象的统一性。
综合相关文献及报道来看,散射成像技术有以下重点、难点亟需突破:1)波前整形技术的实时性仍需进一步提高,在实际应用过程中,散射介质往往具有时变特性,难以实现透过动态散射介质成像,可采用一种快速且准确的光学传输矩阵测量方法,实现透过动态散射介质聚焦或成像;2)目前的大部分波前整形技术的能量利用率较低,可设计一种能量利用率高且便于操控的波前整形技术,在更高能量利用率且快速的前提下实现透过散射聚焦;3)目前已有的波前整形技术无法实现透过散射介质的实时观测,需要设计一种无需重建过程的透过散射介质成像方法,如通过光学补偿的方式实现透过散射介质的实时成像;4)目前基于散斑相关成像方法的成像范围受限于光学记忆效应,拓展光学记忆效应的范围是未来研究的趋势,如通过研究散射介质的物理特性,利用介质自身物理特性(开通道和闭通道)实现光学记忆效应的拓展,从而实现透过散射介质的大视场成像;5)所有的基于点扩展函数工程的透过散射介质成像技术都需要提前测量系统的点扩展函数,如何实现无需测量点扩展函数的成像是未来的研究难点。
除了解决上述难点以外,透过散射介质成像未来较有意义的研究方向为:1)如何实现介质内的快速聚焦或成像,这将对生物医学观测或疾病检测具有重要意义;2)目前的部分研究结果表明,大气或云雾等自然界复杂介质具有散射特性,研究大气或云雾等复杂介质形成散斑的条件,实现透过自然界复杂介质成像对远距离探测和信息感知具有重要意义;3)目前大部分透过散射介质成像的方法只能实现单一物理量的探测,如何通过单帧散斑实现多物理量的探测,对目标信号的探测、识别和分辨具有重要意义;4)随着微纳加工技术的进一步发展,通过定制散射介质实现定制成像,将会对新型化、小型化和集成化的探测成像系统的发展具有深远意义;5)随着深度学习和人工智能技术的发展,将人工智能技术引入到散射成像,通过少量先验知识或通过先验知识迁移的方式实现透过散射介质成像。
综上所述,散射成像已在光学显微、光通信和生物医学领域取得了一系列成果,但是许多问题仍待解决。散射成像技术将向着更高分辨率、更远作用距离和更快速的方向发展。可以预见,随着科研人员对光散射问题的深入研究,上述问题都能够被有效解决,并且会有更多的突破性进展,使散射成像技术能够被广泛应用于人们的日常生活及生产中。
[1] 姜文汉. 自适应光学技术[J]. 自然杂志, 2006, 28(1): 7-13.
Jiang W H. Adaptive optical technology[J]. Chinese Journal of Nature, 2006, 28(1): 7-13.
[2] 钟炜, 张晓晖, 管风, 等. 基于高重复频率脉冲激光的水下全选通成像雷达[J]. 中国激光, 2016, 43(11): 1101009.
[3] 黄子恒, 李微, 杨克成, 等. 水下激光距离选通三维成像方法[J]. 激光与红外, 2016, 46(11): 1315-1319.
Huang Z H, Li W, Yang K C, et al. Underwater laser range-gated 3-D imaging method[J]. Laser & Infrared, 2016, 46(11): 1315-1319.
[4] Huang D, Swanson E, Lin C, et al. Optical coherence tomography[J]. Science, 1991, 254(5035): 1178-1181.
[5] Webb R H. Confocal optical microscopy[J]. Reports on Progress in Physics, 1996, 59(3): 427-471.
[6] Denk W, Strickler J, Webb W. Two-photon laser scanning fluorescence microscopy[J]. Science, 1990, 248(4951): 73-76.
[7] Helmchen F, Denk W. Deep tissue two-photon microscopy[J]. Nature Methods, 2005, 2(12): 932-940.
[8] Zhang H F, Maslov K, Stoica G, et al. Functional photoacoustic microscopy for high-resolution and noninvasive in vivo imaging[J]. Nature Biotechnology, 2006, 24(7): 848-851.
[9] Wang L V. Multiscale photoacoustic microscopy and computed tomography[J]. Nature Photonics, 2009, 3(9): 503-509.
[11] Razansky D, Buehler A, Ntziachristos V. Volumetric real-time multispectral optoacoustic tomography of biomarkers[J]. Nature Protocols, 2011, 6(8): 1121-1129.
[12] Ntziachristos V. Going deeper than microscopy: the optical imaging frontier in biology[J]. Nature Methods, 2010, 7(8): 603-614.
[13] Yaqoob Z, Psaltis D, Feld M S, et al. Optical phase conjugation for turbidity suppression in biological samples[J]. Nature Photonics, 2008, 2(2): 110-115.
[14] Vellekoop I M, Mosk A P. Focusing coherent light through opaque strongly scattering media[J]. Optics Letters, 2007, 32(16): 2309-2311.
[15] Vellekoop I M, Lagendijk A, Mosk A P. Exploiting disorder for perfect focusing[J]. Nature Photonics, 2010, 4(5): 320-322.
[17] Blochet B, Bourdieu L, Gigan S. Fast wavefront optimization for focusing through biological tissue (conference presentation)[J]. Proceedings of SPIE, 2017, 10073: 100730T.
[18] Popoff S M, Lerosey G, Carminati R, et al. Measuring the transmission matrix in optics: an approach to the study and control of light propagation in disordered media[J]. Physical Review Letters, 2010, 104(10): 100601.
[19] Popoff M, Lerosey G, Fink M, et al. Controlling light through optical disordered media: transmission matrix approach[J]. New Journal of Physics, 2011, 13(12): 123021.
[20] 杨虹, 黄远辉, 龚昌妹, 等. 散射介质超衍射极限技术研究进展[J]. 中国光学, 2014, 7(1): 1-25.
[21] Freund I, Rosenbluh M, Feng S C. Memory effects in propagation of optical waves through disordered media[J]. Physical Review Letters, 1988, 61(20): 2328-2331.
[24] Sahoo S K, Tang D L, Dang C. Single-shot multispectral imaging with a monochromatic camera[J]. Optica, 2017, 4(10): 1209-1213.
[25] Li Y Z, Xue Y J, Tian L. Deep speckle correlation: a deep learning approach toward scalable imaging through scattering media[J]. Optica, 2018, 5(10): 1181-1190.
[26] Derode A, Roux P, Fink M. Robust acoustic time reversal with high-order multiple scattering[J]. Physical Review Letters, 1995, 75(23): 4206-4209.
[27] Draeger C, Fink M. One-channel time reversal of elastic waves in a chaotic 2D-silicon cavity[J]. Physical Review Letters, 1997, 79(3): 407-410.
[28] Leith E N, Upatnieks J. Holographic imagery through diffusing media[J]. Journal of the Optical Society of America, 1966, 56(4): 523.
[29] Fink M, Prada C. Acoustic time-reversal mirrors[J]. Inverse problems, 2001, 17(1): R1-R38.
[30] Fisher RA. Optical phase conjugation[M]. New York: Academic Press, 2012: 45- 70.
[31] Voronin É S, Ivakhnik V V, Petnikova V M, et al. Compensation for phase distortions by three-frequency parametric interaction[J]. Soviet Journal of Quantum Electronics, 1979, 9(6): 765-768.
[32] Bloom D M, Bjorklund G C. Conjugate wave-front generation and image reconstruction by four-wave mixing[J]. Applied Physics Letters, 1977, 31(9): 592-594.
[33] Yariv A, Pepper D M. Amplified reflection, phase conjugation, and oscillation in degenerate four-wave mixing[J]. Optics Letters, 1977, 1(1): 16-18.
[34] Králiková B, Skála J, Straka P, et al. Image restoration in a highly non-steady-state regime of stimulated Brillouin scattering in a photodissociation iodine laser[J]. Optics Letters, 1997, 22(11): 766-768.
[35] Karaguleff C, Clark G L. Optical aberration correction by real-time holography in liquid crystals[J]. Optics Letters, 1990, 15(14): 820-822.
[36] Paurisse M, Hanna M, Druon F, et al. Phase and amplitude control of a multimode LMA fiber beam by use of digital holography[J]. Optics Express, 2009, 17(15): 13000-13008.
[37] Cui M, Yang C. Implementation of a digital optical phase conjugation system and its application to study the robustness of turbidity suppression by phase conjugation[J]. Optics Express, 2010, 18(4): 3444-3455.
[38] Lhermite J, Suran E, Kermene V, et al. Coherent combining of 49 laser beams from a multiple core optical fiber by a spatial light modulator[J]. Optics Express, 2010, 18(5): 4783-4789.
[39] van Putten E G, Akbulut D, Bertolotti J, et al. Scattering lens resolves sub-100 nm structures with visible light[J]. Physical Review Letters, 2011, 106(19): 193905.
[40] Cui M. Parallel wavefront optimization method for focusing light through random scattering media[J]. Optics Letters, 2011, 36(6): 870-872.
[42] Park J H, Yu Z, Lee K R, et al. Perspective: wavefront shaping techniques for controlling multiple light scattering in biological tissues: toward in vivo applications[J]. APL Photonics, 2018, 3(10): 100901.
[43] Lemoult F. Lerosey G, de Rosny J, et al. Manipulating spatiotemporal degrees of freedom of waves in random media[J]. Physical Review Letters, 2009, 103(17): 173902.
[44] Lerosey G, de Rosny J, Tourin A, et al. Focusing beyond the diffraction limit with far-field time reversal[J]. Science, 2007, 315(5815): 1120-1122.
[45] Kogelnik H, Pennington K S. Holographic imaging through a random medium[J]. Journal of the Optical Society of America, 1968, 58(2): 273-274.
[46] ShengP. Introduction to wave scattering, localization, and mesoscopic phenomena[M]. New York: Academic Press, 1995: 18- 35.
[47] Choi Y, Yang T D, Fang-Yen C, et al. Overcoming the diffraction limit using multiple light scattering in a highly disordered medium[J]. Physical Review Letters, 2011, 107(2): 023902.
[48] Drémeau A, Liutkus A, Martina D, et al. Reference-less measurement of the transmission matrix of a highly scattering material using a DMD and phase retrieval techniques[J]. Optics Express, 2015, 23(9): 11898-11911.
[49] Yoon J, Lee K, Park J, et al. Measuring optical transmission matrices by wavefront shaping[J]. Optics Express, 2015, 23(8): 10158-10167.
[51] Andreoli D, Volpe G, Popoff S, et al. Deterministic control of broadband light through a multiply scattering medium via the multispectral transmission matrix[J]. Scientific Reports, 2015, 5: 10347.
[52] Mounaix M, Andreoli D, Defienne H, et al. Spatiotemporal coherent control of light through a multiple scattering medium with the multispectral transmission matrix[J]. Physical Review Letters, 2016, 116(25): 253901.
[53] DongJ, KrzakalaF, GiganS. Spectral method for multiplexed phase retrieval and application in optical imaging in complex media[C]∥ICASSP 2019-2019 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), May 12-17, 2019, Brighton, United Kingdom. New York: IEEE, 2019: 4963- 4967.
[54] Park H, Crozier K B. Multispectral imaging with vertical silicon nanowires[J]. Scientific Reports, 2013, 3: 2460.
[55] Stewart J W, Akselrod G M, Smith D R, et al. Toward multispectral imaging with colloidal metasurface pixels[J]. Advanced Materials, 2017, 29(6): 1602971.
[57] Pégard N C, Liu H Y, Antipa N, et al. Compressive light-field microscopy for 3D neural activity recording[J]. Optica, 2016, 3(5): 517-524.
[58] AkkermansE, MontambauxG. Mesoscopic physics of electrons and photons[M]. UK: Cambridge University Press, 2007: 35- 45.
[59] Goodman JW. Speckle phenomena in optics: theory and applications[M]. USA: Roberts and Company Publishers, 2007: 15- 70.
[60] Fienup J R. Reconstruction of an object from the modulus of its Fourier transform[J]. Optics Letters, 1978, 3(1): 27-29.
[61] Fienup J R. Phase retrieval algorithms: a comparison[J]. Applied Optics, 1982, 21(15): 2758-2769.
[62] Wu T F, Katz O, Shao X P, et al. Single-shot diffraction-limited imaging through scattering layers via bispectrum analysis[J]. Optics Letters, 2016, 41(21): 5003-5006.
[63] Lohmann A W, Weigelt G, Wirnitzer B. Speckle masking in astronomy: triple correlation theory and applications[J]. Applied Optics, 1983, 22(24): 4028-4037.
[64] Edrei E, Scarcelli G. Optical imaging through dynamic turbid media using the Fourier-domain shower-curtain effect[J]. Optica, 2016, 3(1): 71-74.
[66] Yang W Q, Li G W, Guohai S T. Imaging through scattering media with the auxiliary of a known reference object[J]. Scientific Reports, 2018, 8: 9614.
[67] Wang X Y, Jin X, Li J Q, et al. Prior-information-free single-shot scattering imaging beyond the memory effect[J]. Optics Letters, 2019, 44(6): 1423-1426.
[68] Jin X, Wang Z P, Wang X Y, et al. Depth of field extended scattering imaging by light field estimation[J]. Optics Letters, 2018, 43(20): 4871-4874.
[69] Liao M H, Lu D J, Pedrini G, et al. Extending the depth-of-field of imaging systems with a scattering diffuser[J]. Scientific Reports, 2019, 9: 7165.
[70] Candès E J, Wakin M B. An introduction to compressive sampling[J]. IEEE Signal Processing Magazine, 2008, 25(2): 21-30.
[71] Rudin L I, Osher S, Fatemi E. Nonlinear total variation based noise removal algorithms[J]. Physica D: Nonlinear Phenomena, 1992, 60: 259-268.
[72] Brady D J, Choi K, Marks D L, et al. Compressive holography[J]. Optics Express, 2009, 17(15): 13040-13049.
[73] Mukherjee S, Vijayakumar A, Kumar M, et al. 3D imaging through scatterers with interferenceless optical system[J]. Scientific Reports, 2018, 8: 1134.
[74] Shi Y Y, Liu Y W, Wang J M, et al. Non-invasive depth-resolved imaging through scattering layers via speckle correlations and parallax[J]. Applied Physics Letters, 2017, 110(23): 231101.
[75] Edrei E, Scarcelli G. Memory-effect based deconvolution microscopy for super-resolution imaging through scattering media[J]. Scientific Reports, 2016, 6: 33558.
[76] Richardson W H. Bayesian-based iterative method of image restoration[J]. Journal of the Optical Society of America, 1972, 62(1): 55-59.
[77] Lucy L B. An iterative technique for the rectification of observed distributions[J]. The Astronomical Journal, 1974, 79(6): 745-754.
[78] Wu T F, Dong J, Shao X P, et al. Imaging through a thin scattering layer and jointly retrieving the point-spread-function using phase-diversity[J]. Optics Express, 2017, 25(22): 27182-27194.
[79] Antipa N, Kuo G, Heckel R, et al. DiffuserCam: lensless single-exposure 3D imaging[J]. Optica, 2018, 5(1): 1-9.
[80] Liao M H, He W Q, Lu D J, et al. Ciphertext-only attack on optical cryptosystem with spatially incoherent illumination: from the view of imaging through scattering medium[J]. Scientific Reports, 2017, 7: 41789.
[81] Li G W, Yang W Q, Li D Y, et al. Cyphertext-only attack on the double random-phase encryption: experimental demonstration[J]. Optics Express, 2017, 25(8): 8690-8697.
[82] Wu P F, Liang Z, Zhao X, et al. Lensless wide-field single-shot imaging through turbid media based on object-modulated speckles[J]. Applied Optics, 2017, 56(12): 3335-3341.
[84] Li G, Yang W, Wang H, et al. Image transmission through scattering media using ptychographic iterative engine[J]. Applied Sciences, 2019, 9(5): 849.
[85] Singh A K, Naik D N, Pedrini G, et al. Exploiting scattering media for exploring 3D objects[J]. Light: Science & Applications, 2017, 6(2): e16219.
[86] Lee K, Park Y. Exploiting the speckle-correlation scattering matrix for a compact reference-free holographic image sensor[J]. Nature Communications, 2016, 7: 13359.
[87] Yilmaz H, van Putten E G, Bertolotti J, et al. Speckle correlation resolution enhancement of wide-field fluorescence imaging[J]. Optica, 2015, 2(5): 424-429.
Article Outline
朱磊, 邵晓鹏. 散射成像技术的研究进展[J]. 光学学报, 2020, 40(1): 0111005. Lei Zhu, Xiaopeng Shao. Research Progress on Scattering Imaging Technology[J]. Acta Optica Sinica, 2020, 40(1): 0111005.