Fast axial scanning system based on digital micromirror device
-
摘要: 为解决双光子荧光显微成像系统轴向扫描问题,提出一种基于数字微镜(DMD)的快速轴向扫描系统。该系统采用DMD选择光路,不同光路放置不同焦距的透镜组对光束发散角产生不同的改变,经物镜聚焦后得到不同深度的轴上扫描点。对该系统的轴向扫描距离、扫描点位置及衍射效率进行了理论计算仿真,结果表明:扫描系统采用4个模块以及5个模块时其轴向扫描距离均可达到1 mm,4模块系统中透镜的焦距为297.3 mm,5模块系统中透镜焦距为361.47 mm。轴向扫描点除边缘点外线性分布, 轴向扫描频率达到几十kHz,满足脑神经成像的要求。Abstract: Aiming at solving axial scanning problem of dual photon fluorescence microscopy imaging system, fast axial scanning system based on digital micromirror device (DMD) is proposed in this paper. DMD is used for selecting optical path in this system, and lens groups with different focal lengths are placed in these optical paths to generate different changes of beam divergence angle. Then axial scanning spots with different depth are obtained after beam focused by objective lens. Axial scanning distance, scanning spot and diffraction efficiency are calculated and simulated. Results show that when scanning system adopts four modules or five modules, axial scanning distance can reach 1mm, focal length of lens in four modules and five modules is 297.3 mm and 361.47 mm respectively. Axial scanning spot varies in good linearity except marginal spot. Axial scanning frequency approaches to dozens of kilohertz. This scanning system can meet requirements of cranial nerve imaging.
-
Keywords:
- optical design /
- axial scanning /
- ZEMAX /
- two photon fluorescence imaging system /
- DMD
-
引言
随着技术的不断进步,光电类武器不断进行升级,我们正朝着多波段、宽视角、超远程和高分辨率的方向发展[1],对不同构型的样品进行应力分布测量成为当前的首要解决问题。应力双折射测量的方法主要有干涉法[2]、Searnmont补偿法[3]、调制法[4,5]等。英国夏普应力工程有限公司的Searnmont补偿法偏光应力仪S-66[6,7]采用钠光灯光源,实现的测量精度为0.2 nm;德国生产的自动应力测试系列产品,型号为StrainMatic,测量精度为±0.1 nm,但是样品测量范围只有36×36 mm[8,9];目前国内生产的应力仪主要是基于干涉法和Searnmont补偿法的偏振分析原理,比较有代表性的产品是台偏式应力仪PSV-702,该仪器利用Searnmont补偿法,虽然价格低廉,但在操作过程中需要旋转检偏器,无法满足对于宽波段光学材料应力测试的要求[10,11]。以上所提及的型号产品,无法保证实现对多构型样品的应力双折射分布测量。
而弹光调制器由于拥有调制频谱宽、调制纯度高、波段范围广的优点,已被广泛应用于双折射检测、光谱测量等应用领域中[12,13]。采用弹光调制器的双折射测量技术,精度较高,无需机械转动,仅需一次测量就可以实现对快轴角以及相位延迟量的测量[14]。
针对目前对不同构型样品实现高精度应力分布测量问题,本文在双弹光调制的基础上设计了一种应用于双折射分布测量的多构型扫描系统,在保证高精度的同时设计出不同扫描算法,实现了对样品任一区域进行快速、准确的应力双折射分布测量,并进行了实验验证。
1 测量原理
如图1所示,将待测样品放置于两个弹光调制器中间,激光器的激光经此测量系统后,被探测器接收[15]。
用Muller矩阵能够全面反映目标的偏振光学特性。根据所搭建的光路系统,其Stokes参量与Muller矩阵关系如下:
$$ {S}_{{\rm{out}}}={\boldsymbol{M}}_{{\rm{A}}}{\boldsymbol{M}}_{{\rm{PE}}{{\rm{M}}}{2}}{\boldsymbol{M}}_{{\rm{sample}}}{\boldsymbol{M}}_{{\rm{PE}}{{\rm{M}}}{1}}{\boldsymbol{M}}_{{\rm{p}}}{S}_{{\rm{{\rm{in}}}}} $$ (1) 式中:${\boldsymbol{M}}_{{\rm{p}}}$、${\boldsymbol{M}}_{{\rm{A}}}$为代表起偏与检偏的Muller矩阵;${\boldsymbol{M}}_{{\rm{PE}}{{\rm{M}}}{1}}$、${\boldsymbol{M}}_{{\rm{PE}}{{\rm{M}}}{2}}$为代表弹光调制器1、2的Muller矩阵;${\boldsymbol{M}}_{{\rm{sample}}}$为光学材料的Muller矩阵;${S}_{{\rm{in}}}$为入射光信号。
待测光学样品的Muller矩阵展开式为
$$ {\boldsymbol{M}}_{{\rm{sample}}}=\left[\begin{array}{cccc}1& 0& 1& 0\\ 0& {\mathrm{cos}}^{2}\left(2\rho \right)+{\mathrm{sin}}^{2}\left(2\rho \right)\mathrm{cos}X& \mathrm{s}\mathrm{i}\mathrm{n}\left(4\rho \right){\mathrm{sin}}^{2}\left(\dfrac{X}{2}\right)& -\mathrm{s}\mathrm{i}\mathrm{n}\left(2\rho \right)\mathrm{sin}X\\ 0& \mathrm{s}\mathrm{i}\mathrm{n}\left(4\rho \right){\mathrm{sin}}^{2}\left(\dfrac{X}{2}\right)& {\mathrm{sin}}^{2}\left(2\rho \right)+\mathrm{c}\mathrm{o}\mathrm{s}\left(2\rho \right)\mathrm{cos}X& \mathrm{c}\mathrm{o}\mathrm{s}\left(2\rho \right)\mathrm{sin}X\\ 0& \mathrm{s}\mathrm{i}\mathrm{n}\left(2\rho \right)\mathrm{sin}X& -\mathrm{c}\mathrm{o}\mathrm{s}\left(2\rho \right)\mathrm{sin}X& \mathrm{cos}X\end{array}\right] $$ (2) 式中:$ \rho $表示双折射的快轴方位角;$ X $则表示相位延迟量。
由于探测器得到的光强仅为Stokes矢量的第1个分量,则将各个器件的Muller矩阵带入到式(1)中即可得到探测器探测的光强为
$$ \begin{split} I = &\frac{K{I}_{0}}{2}[1 + \mathrm{sin}{\delta }_{1}\mathrm{sin}{\delta }_{2}\mathrm{cos}X + \mathrm{cos}{\delta }_{1}\mathrm{cos}{\delta }_{2}\mathrm{sin}\left(4\rho \right){\mathrm{sin}}^{2}\left(\frac{X}{2}\right)+\\ &\mathrm{cos}{\delta }_{1}\mathrm{sin}{\delta }_{2}\mathrm{c}\mathrm{o}\mathrm{s}\left(2\rho \right)\mathrm{sin}X+ \mathrm{sin}{\delta }_{1}\mathrm{cos}{\delta }_{2}\mathrm{s}\mathrm{i}\mathrm{n}\left(2\rho \right)\mathrm{sin}X] \end{split} $$ (3) 式中:$ K $为光电常数;$ {I}_{0} $为激光经过起偏器后的光强;$ {\delta }_{1} $和$ {\delta }_{2} $为弹光调制器1、2的相位调制幅度,可进一步表示为
$$ \left\{\begin{array}{l}{\delta }_{1}={\delta }_{10}\mathrm{sin}{w}_{1}t\\ {\delta }_{2}={\delta }_{20}\mathrm{sin}{w}_{2}t\end{array}\right. $$ (4) 式中:$ {w}_{1} $和$ {w}_{2} $分别表示弹光调制器1、2的频率;$ {\delta }_{10} $和$ {\delta }_{20} $分别表示弹光调制器1、2的相位调制幅度。
将式(4)反带入式(3),并用第一类贝塞尔级数展开得:
$$ \begin{split} I=&\frac{K{I}_{0}}{2}[1+{J}_{0}\left({\delta }_{10}\right){J}_{0}\left({\delta }_{20}\right)\mathrm{sin}\left(4\rho \right){\mathrm{sin}}^{2}\left(\frac{X}{2}\right)+K+\\ &2{J}_{1}\left({\delta }_{10}\right){J}_{1}\left({\delta }_{20}\right)\mathrm{c}\mathrm{o}\mathrm{s}\left({w}_{1}t-{w}_{2}t\right)\mathrm{cos}X+\dots +\\ &2{J}_{3}\left({\delta }_{10}\right){J}_{3}\left({\delta }_{20}\right)\mathrm{s}\mathrm{i}\mathrm{n}\left({3w}_{1}t-{2w}_{2}t\right)\mathrm{s}\mathrm{i}\mathrm{n}\left(2\rho \right)\mathrm{cos}X] \end{split} $$ (5) 得到式(5)所示的探测器光强表达式后,通过数字锁相技术,得到不同倍频的频率信号,并据此求出相位延迟量$ \rho $和快轴方位角$ X $。对式(5)中的信号幅值进行处理可以得到:
$$ \left\{\begin{array}{l}{V}_{\left({w}_{2}-{w}_{1}\right)}=K{I}_{0}{J}_{1}\left({\delta }_{10}\right){J}_{1}\left({\delta }_{20}\right)\mathrm{cos}X\\ {V}_{{w}_{1}}=K{I}_{0}{J}_{0}\left({\delta }_{20}\right){J}_{0}\left({\delta }_{20}\right)\sin\left(2\rho \right)\mathrm{sin}X\\ {V}_{{w}_{2}}=K{I}_{0}{J}_{0}\left({\delta }_{20}\right){J}_{1}\left({\delta }_{20}\right)\cos\left(2\rho \right)\mathrm{sin}X\end{array}\right. $$ (6) 定义比值$ {R}_{1} $和$ {R}_{2} $为
$$ \left\{\begin{array}{l}{R}_{1}=\dfrac{{V}_{{w}_{1}}{J}_{1}\left({\delta }_{20}\right)}{{V}_{\left({w}_{2}-{w}_{1}\right)}{J}_{0}\left({\delta }_{20}\right)}=\sin\left(2\rho \right)\mathrm{tan}X\\ {R}_{2}=\dfrac{{V}_{{w}_{2}}{J}_{1}\left({\delta }_{10}\right)}{{V}_{\left({w}_{2}-{w}_{1}\right)}{J}_{0}\left({\delta }_{10}\right)}=\cos\left(2\rho \right)\mathrm{tan}X\end{array}\right. $$ (7) 则快轴方位角$ \rho $和相位延迟量$ X $可表示为
$$ \left\{\begin{array}{l}\rho =\dfrac{1}{2}\mathit{{\rm{arctan}}}\dfrac{{R}_{1}}{{R}_{2}}\\ X={\rm{arctan}}\sqrt{{R}_{1}^{2}+{R}_{2}^{2}}\end{array}\right. $$ (8) 最后,根据求得的快轴方位角$ X $和入射光波长$ \lambda $,可得到双折射的延迟量$ \Delta $为
$$ \Delta =X\frac{\lambda }{2\pi } $$ (9) 由式(9)可以得出,根据弹光调制器的相位调制幅值,通过对不同样品的不同频率幅值进行测量,能够实现对样品单一点同时进行双折射延迟量和快轴方位角的测量。
对于不同构型的待测样品,设计出相应的扫描算法,实现对样品任一区域进行应力双折射分布测量。
2 测量系统整体设计
系统整体设计框图如图2所示,将测量光路搭建完毕后,使用弹光驱动电路使两个弹光调制器处于谐振工作状态。借助二维扫描平台,使用可移动导轨固定样品,采用单一点的双折射测量原理,可实现对不同构型样品的面阵式应力双折射分布测量。
在系统对待测样品进行面阵式测量过程中,二维扫描平台会进行位移,所以位移精度也是测试系统的一个重要参数,该项参数会依据光栅尺返回的位置信息和目标位置信息参考进行对比。在实际测试中,最小的位移单位为mm,所以只需要保证最大位移误差在0.1 mm即5%内,就不会对应力的测量结果产生影响。
2.1 扫描算法设计模块
方形扫描算法的步骤为:首先确定x、y轴的起始与终点坐标信息,判断出扫描距离与寸动次数;x轴开始寸动,在到达终点后进行位置判断,确定终点位置坐标信息后,y轴寸动一次,且x轴随即往反方向运动,并进行终点位置判断;在扫描范围内反复循环,直到x、y轴都运动到终点位置处,完成此次扫描。
圆域扫描算法的步骤为:首先确定待测圆域半径,根据半径推算出扫描起始点以及扫描距离,后续扫描方式与方形扫描方式一致。在该扫描算法运行时,会对圆域的扫描精度进行判断,如果该精度达不到2 mm×2 mm,则舍去该扫描范围对下一区域进行扫描,该判断条件主要针对于圆形待测样品的圆弧部分。扫描模块流程图如图3所示。
2.2 数据处理模块
每次传输时,FPGA模块会通过USB接口向上位机传输一包数据,上位机拿到数据后根据相关公式和算法对其进行处理。根据每个位置的数据所代表的信息,进行相应的计算。其中根据上位机获取平移台的实时位置信息,利用FPGA的数字锁相模块找到各倍频项的正余弦值,然后通过计算得到各倍频项的幅值,根据测量原理利用得到的幅值做除法运算消除无关项,最终计算延迟量。数据处理模块流程图如图4所示。
2.3 显示模块
为了更直观地显示样品应力分布的情况,上位机根据接收到的样品位置信息和计算出的延迟量、快轴方位角信息,在二维显示图中进行直观的显示,其中延迟量用不同的颜色来表示,快轴方位角由每个格子中心斜线的倾斜度来表示。
上位机通过绘制色列图的方式对其进行显示。色列图的3个颜色显示和透明度显示都可以用0~255来表示。结合上述颜色显示的机理,在程序中开辟一个100×4位的数组,其中100代表色列图由100个色块组成,4代表每块色块由A、R、G、B这4个元素组成。
由255×色彩系数就等于该原色的真实数值,A为固定值255,表示完全不透明。按照该色列图填色原理,能得出颜色分布均匀的色列图。当完成含有100个色块的色列图后,只需把计算出的光程差比上探测器的波长即可算出此光程差对应的比例,数值越大,颜色越深,然后根据此比例在色列图的对应位置找到相应颜色即可。
上位机最终设计界面如图5所示,界面左边为括弹光控制和二维扫描平台控模块以及单点扫描模块;界面中间为双折射测量的二维分布显示;界面右边为数据处理模块,用于显示数据分布及占比情况。
3 实验验证
为了验证该系统的可行性,首先对扫描平台进行二维扫描精度测量实验;其次对不同构型样品进行重复性测试,分析其精度。
3.1 二维扫描精度实验
该部分实验主要从连续扫描误差分析及连续寸动扫描误差分析对该系统的测量精度进行验证。
3.1.1 连续扫描误差分析
本次实验在所搭建的系统框架下,对扫描平台的x轴与y轴进行了各自150 mm的连续扫描,将扫描过程中所得到的实际坐标值与标准坐标值进行误差分析,误差分析图分别如图6和图7所示。在连续扫描实验中,x轴与y轴的连续扫描误差均小于0.005 mm,满足连续扫描步进误差要求。
3.1.2 连续寸动扫描误差分析
连续寸动扫描时,x轴与y轴每次步进2 mm,则各自步进150 mm。记录两轴每次步进的位置,与标准坐标位置对比并进行误差分析,误差分析图分别如图8和图9所示。
在连续寸动扫描实验中,x轴与y轴的连续扫描误差均小于0.009 mm,满足连续寸动扫描步进误差要求。
3.2 样品扫描及重复性测试
该实验分为延迟量和快轴方位角重复性测试。在测试过程中,圆形扫描测试样品选择633 nm(1/4)玻片,方形扫描测试样品选择BK7玻璃。首先,对这两种不同构型的样品进行单点扫描,测试其重复性。其中1/4玻片旋转到164°,在该点进行单点扫描,采集2 500组左右的数据进行重复性分析;对BK7玻璃样品施加应力,在其施加应力处进行单点扫描。测试结果如图10和图11所示。
选取1/4玻片的164°及BK7玻璃的89°进行快轴方位角进行重复性测试,测试结果如图12和图13所示。不同样品的重复性测量结果如表1所示。
对重复性试验所测得数据进行分析,633 nm(1/4)玻片的延迟量波动范围为0.12 nm,标准差为0.035 2 nm;快轴方位角波动范围为0.13°,标准差为0.022 5°。BK7玻璃样品的延迟量波动范围为0.25 nm,标准差为0.038 9 nm;快轴方位角波动范围为0.13°,标准差为0.037 0°,均满足重复性测试要求且重复性良好。其中,测量数据产生波动的原因是光源的不稳定会造成数据测试的偏差。
表 1 不同样品的重复性测量结果Table 1. Repeatability measurement results of different samples样品 测量数据类型 波动范围 标准差 633 nm(1/4)玻片 延迟量 0.12 nm 0.035 2 nm 快轴方位角 0.13° 0.022 5° BK7玻璃 延迟量 0.25 nm 0.038 9 nm 快轴方位角 0.13° 0.037 0° 在保证对样品的高精度测量下,对两个不同构型样品进行面阵式扫描,扫描结果如图14和图15所示。
在图14和图15中,x、y轴均为样品所在扫描平台的坐标,扫描平台的x、y轴坐标范围均为0~150 mm,不同颜色代表相位延迟量的不同。该测试玻片的允许误差为$\lambda$%,即测试结果在152 nm~166 nm之间即可满足要求。在图14所示的测试结果中,该玻片的延迟量分布普遍在153 nm左右,最大为164 nm,满足测量精度要求。
4 结论
本文设计了一种应用于双折射分布测量的多构型扫描系统,实现了对不同构型样品进行高精度应力双折射分布测量。首先采用单点扫描验证其测量的高精度性,其次采用面阵式扫描测量其双折射分布。其中,对圆形样品633 nm(1/4)玻片进行测试,测量误差最小为0.79%,最大误差为0.95%,单点测量延迟量波动范围为0.12 nm,标准差为0.035 2 nm;快轴方位角波动范围为0.13°,标准差为0.022 5°。对方形样品BK7玻璃样品施加应力进行测试,在195.94 nm处进行单点测量,延迟量波动范围为0.25 nm,标准差为0.038 9 nm;快轴方位角波动范围为0.13°,标准差为0.037 0°。测量精度方面,连续扫描过程中,误差最大为2.5%,而在连续寸动扫描过程中,最大误差为0.45%,均远高于允许误差范围。验证了该测量系统的高精度性。此系统的主要优势在于无需转动待测样品或检偏器,并且在保证高精度的同时实现对样品的面阵完整性测量。
-
[1] Masters B R, So P T, Gratton E. Multiphoton excitation fluorescence microscopy and spectroscopy of in vivo human skin[J]. Biophysical Journal, 1997, 72(2): 2405-2412. doi: 10.1016-S0006-3495(97)78886-6/
[2] Bewersdorf J, Pick R, Hell S W. Multifocal multiphoton microscopy[J]. Optics Letters, 1998, 23(9): 655-7. doi: 10.1364/OL.23.000655
[3] Denk W, Delaney K R, Gelperin A, et al. Anatomical and functional imaging of neurons using 2-photon laser scanning microscopy[J]. Journal of Neuroscience Methods, 1994, 54(2): 151-62. doi: 10.1016/0165-0270(94)90189-9
[4] Botcherby E J, Juskaitis R, Booth M J, et al. An optical technique for remote focusing in microscopy[J]. Optics Communications, 2008, 281(4): 880-887. doi: 10.1016/j.optcom.2007.10.007
[5] 张运波, 侯文玫, 郑继红.采用变焦液体透镜的共焦检测系统的设计与仿真[J].应用光学, 2012, 33(1):14-18. http://d.old.wanfangdata.com.cn/Periodical/yygx201201004 Zhang Yunbo, Hou Wenmei, Zheng Jihong. Design and simulation of confocal detection system utilizing tunable-focus liquid lens[J]. Journal of Applied Optics, 2012, 33(1):14-18. http://d.old.wanfangdata.com.cn/Periodical/yygx201201004
[6] 李佃蒙.基于纯相位液晶空间光调制器的共焦扫描方法研究[D].哈尔滨: 哈尔滨工业大学, 2013. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=D420960 Li Dianmeng. Research of confocal scanning method based on phase-only liquid crystal spatial light modulator[D]. Harbin: Harbin Institute of Technology, 2013. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=D420960
[7] 董伟辉, 谢永军, 李恩玲.利用镜面形变实现共轴折反射式变焦光学系统设计[J].应用光学, 2010, 31(6):893-897. doi: 10.3969/j.issn.1002-2082.2010.06.005 Dong Weihui, Xie Yongjun, Li Enling. Design of coaxial catadioptric zoom system using deformable mirrors[J]. Journal of Applied Optics, 2010, 21(6):893-897. doi: 10.3969/j.issn.1002-2082.2010.06.005
[8] 胡少磊.基于主动光学的空间变焦成像技术研究[D].西安: 中国科学院西安光学精密机械研究所, 2012. http://cdmd.cnki.com.cn/Article/CDMD-80142-1013010147.htm Hu Shaolei. Based on active optical imaging technology for space zoom[D]. Xi'an: Xi'an Institute of Optics and Precision Mechanics of CAS, 2012. http://cdmd.cnki.com.cn/Article/CDMD-80142-1013010147.htm
[9] Grewe B F, Voigt F F, Hoff M V, et al. Fast two-layer two-photon imaging of neuronal cell populations using an electrically tunable lens[J]. Biomed Opt. Express, 2011, 2(7): 2035-2046. doi: 10.1364/BOE.2.002035
[10] 杜睿.双光子轴向随机扫描方法研究[D].武汉: 华中科技大学, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10487-1011110464.htm Du Rui.The study of axial random-access laser scanning scheme in two photo microscopy[D]. Wuhan: Huazhong University of Science and Technology, 2011. http://cdmd.cnki.com.cn/Article/CDMD-10487-1011110464.htm
[11] Kirkby P A, Srinivas N K M, Silver R A. A compact acousto-optic lens for 2D and 3D femtosecond based 2-photon microscopy[J]. Optics Express, 2010, 18(13):13721-45. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=Open J-Gate000000072359
[12] 吕晓华.随机扫描双光子荧光显微成像系统研究[D].武汉: 华中科技大学, 2006. http://cdmd.cnki.com.cn/Article/CDMD-10487-2009033595.htm Lyu Xiaohua. Research on the random access two-photon fluorescence microscope[D]. Wuhan: Huazhong University of Science and Technology, 2006. http://cdmd.cnki.com.cn/Article/CDMD-10487-2009033595.htm
[13] 姜润华.用于双光子显微成像的宽带声光偏转器研究[D].武汉: 华中科技大学, 2012. http://cdmd.cnki.com.cn/Article/CDMD-10487-1012361713.htm Jiang Runhua. On the wide-band acousto-optic deflector for two-photon microscopy[D]. Wuhan: Huazhong University of Science and Technology, 2012. http://cdmd.cnki.com.cn/Article/CDMD-10487-1012361713.htm
[14] 李金鹏, 陈磊, 方波, 等.动态干涉仪的位相光栅衍射效率研究[J].红外与激光工程, 2015, 44(9):2696-2701. doi: 10.3969/j.issn.1007-2276.2015.09.027 Li Jinpeng, Chen Lei, Fang Bo, et al. Diffraction efficiencies of phase gratings based on dynamic Interferometer[J]. Infrared and Laser Engineering, 2015, 44(9):2696-2701. doi: 10.3969/j.issn.1007-2276.2015.09.027
-
期刊类型引用(5)
1. 杨丰源,李彬华. 一种新的彩色图像幸运区域融合算法. 计算机仿真. 2024(09): 372-377 . 百度学术
2. 李亚伟,张弘,伍凌帆,杨一帆,陈浩. 基于FPGA的实时大气湍流图像复原算法及实现. 应用光学. 2021(06): 1017-1024 . 本站查看
3. 盘亚楠,谌雨章,郭煜玮. 基于像素位移估计的水下湍流退化图像畸形校正算法. 计算机测量与控制. 2018(06): 248-251+256 . 百度学术
4. 朱军,张震,张伦伟. 一种基于图像恢复的水下零件在线测量方法. 机电一体化. 2018(08): 29-33 . 百度学术
5. 朱军,张震,张伦伟. 基于双谱分析优化算法的水下图像复原方法. 机电一体化. 2018(09): 50-54 . 百度学术
其他类型引用(12)