重构定心法提高微球轴向位移测量信噪比.pdf

Optoelectronics 光电子光电子,2023,13(2),43-51 Published Online June 2023 in Hans.https:/www.hanspub.org/journal/oe https:/doi.org/10.12677/oe.2023.132006 文章引用文章引用:原熙英,曾雅楠,赵瑞,邵佳琪,郭昊明,郭启航.重构定心法提高微球轴向位移测量信噪比J.光电子,2023,13(2):43-51.DOI:10.12677/oe.2023.132006 重构定心法提高微球轴向位移测量信噪比重构定心法提高微球轴向位移测量信噪比 原熙英原熙英1，曾雅楠，曾雅楠1,2*，赵赵 瑞瑞1，邵佳琪，邵佳琪1，郭昊明，郭昊明1，郭启航，郭启航1 1天津农学院工程技术学院，天津 2天津大学精密测试技术及仪器国家重点实验室，天津 收稿日期：2023年4月21日；录用日期：2023年6月13日；发布日期：2023年6月28日 摘摘 要要 微球轴向位移的测量容易受到环境噪声干扰，从机理上提高测量的精准性和信噪比，是解决噪声问题的微球轴向位移的测量容易受到环境噪声干扰，从机理上提高测量的精准性和信噪比，是解决噪声问题的重要途经。离轴数字全息显微技术能够测量微球轴向位移重要途经。离轴数字全息显微技术能够测量微球轴向位移1，测量分辨力达到纳米级别，而且对微球的，测量分辨力达到纳米级别，而且对微球的位置没有特殊要求。全息定位标志的精准度是测量精准的重要保障。然而，利用传统的图像定心法定位位置没有特殊要求。全息定位标志的精准度是测量精准的重要保障。然而，利用传统的图像定心法定位微球中心，为全息定位标志的计算引入误差。本文提出重构定心法，用重构的相位中心定位微球中心，微球中心，为全息定位标志的计算引入误差。本文提出重构定心法，用重构的相位中心定位微球中心，校正全息定位标志，提高测量信噪比。校正全息定位标志，提高测量信噪比。通过对通过对5 m直径直径微球测量轴向位移，测量信噪比提升微球测量轴向位移，测量信噪比提升20%，验证，验证重构定心法的有效性。重构定心法的有效性。关键词关键词 微球，数字全息，位移微球，数字全息，位移 The Reconstruction Centering Method to Improve the Signal-to-Noise Ratio of the Axial Displacement Measurement of Microspheres Xiying Yuan1,Yanan Zeng1,2*,Rui Zhao1,Jiaqi Shao1,Haoming Guo1,Qihang Guo1 1College of Engineering and Technology,Tianjin Agricultural College,Tianjin 2State Key Laboratory of Precision Measuring Technology and Instruments,Tianjin University,Tianjin Received:Apr.21st,2023;accepted:Jun.13th,2023;published:Jun.28th,2023 Abstract The measurement of the axial displacement of microspheres is easily disturbed by environmental *通讯作者。原熙英 等 DOI:10.12677/oe.2023.132006 44 光电子 noise,and improving the accuracy and signal-to-noise ratio of the measurement mechanism is an important way to solve the noise problem.Off-axis digital holographic microscopy can measure the axial displacement of microspheres,and the measurement resolution can reach the nanometer level,and there are no special requirements for the position of the microspheres.The accuracy of holo-graphic positioning marks is an important guarantee for accurate measurement.However,using the traditional picture centering method to locate the center of the microspheres can introduce errors in the calculation of holographic positioning marks.This article proposes a reconstruction centering method,which uses the reconstructed phase center to locate the center of the microsphere,corrects the holographic positioning mark,and improves the measurement signal-to-noise ratio.By measur-ing the axial displacement of 5 m diameter microspheres,the signal-to-noise ratio was improved by 20%to verify the effectiveness of the reconstruction centering method.Keywords Microsphere,Digital Holography,Displacement Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 微球作为一种主要粒子，是表征分子扭矩 1 2，流体动力学 3 4，单分子力谱 5，空气中颗粒物研究 6 7等的重要标志。它也是显微镜中用于超分辨率成像的有利工具 8 9。跟踪微球并定量测量其三维运动对于揭示微球实验中涉及的物理或生物学原理至关重要。由于微球的特征尺寸在微米或亚微米级别，因此，光学显微技术是测量移动微球三维位移的主要方法。微球位移的高分辨测量是一项具有挑战性但意义重大的工作。数字全息显微技术(Digital Holographic Microscopy,DHM)具有纳米分辨率 10，而且 DHM 并不需要扫描方式 11。因此，DHM 是一种适用于动态测量微球位移的方法。实际上，基础测量精度的提高为应用提供了更准确的数据基础。关于测量机理的讨论已经发展到研究者更加关注实际测量过程引起的实际因素与理想测量模型之间的差异的地步，例如倾斜面的重建 12和孪生图像的消除 13。离轴数字全息显微技术与传统的重聚焦方式测量微球轴向位移不同，是对透射光线进行追迹，从而定位标记的变化，测量轴向位移。离轴全息提供了在微球位于特殊位置(例如微球显微聚焦像附近时)轴向位移测量的解决方法，克服了同轴全息测量微粒的难点问题。微球轴向位移的测量容易受到环境噪声干扰，从机理上提高测量的精准性和信噪比，是解决噪声问题的重要途经。本文重点关注微球光路搭建中，非平行光及离轴光照射微球造成的微球成像中心与光程长最大值点不一致，因而造成实际光程长差分分析结果与理论结果有差别的问题，从而通过建立重构定心法，重新计算全息定位标志，通过提高标志定位的精准度，从而提高轴向位移测量的测量信噪比。2.重重构定心法测量微球轴向位移构定心法测量微球轴向位移 对实际取得的微球全息图进行重构(5 m 直径二氧化硅微球，浸没在水中，放大倍数 54.77 倍)，微球重构相位最大值与微球成像几何中心并不是同一点，微球相位最大值像素坐标(448,371)，微球成像几何中心(453,367)，相差 x 方向 5 像素，y 方向 4 像素，如图 1 所示。入射微球中心的中心光线具备最大的光程长值，因此中心光线与成像面的交点应为微球重构像为最大值所在的像素点。微球重构像几何中Open AccessOpen Access原熙英 等 DOI:10.12677/oe.2023.132006 45 光电子 心应为微球投影在成像面的圆心。而复合式测量微球位移的方法需要分析光程长最大值与周边像素光程长差值，以此作为轴向位移的移动标志物。Figure 1.The hologram(the interference fringes are removed)and there constructed phase of microsphere 图图 1.微球全息图(已去除干涉条纹)和微球的重构相位 因此传统方法中，中心像素点并不是光程长最大值的点，需重新分析具体情况，评估具体情况与理想情况的偏差，建立新的测量方法，提高测量准确度和信噪比。分析上述现象可知，入射微球的中心光线与成像面不垂直，且入射光束并非平行光束。分析原因：1)机械安装造成的成像光轴与理想情况存在偏差；2)前面光路的平行扩束与理想情况存在偏差。倾斜入射的非平行光束等效于轴外点光源发出的发散光束，具备孔径角和视场。由于轴外物点宽光束经系统成像后存在失对称的情况，即对称与主光线的一对光线经光学系统后，失去对主光线的对称性，交点不再位于主光线上。在光学设计中，这种像差称为彗差。如图 2 所示。Figure 2.Schematic diagram of actual incident microsphere rays 图图 2.实际入射微球光线示意图 原熙英 等 DOI:10.12677/oe.2023.132006 46 光电子 子午面是 AB 起点构成的面，弧矢面是 AB起点构成的面。ABABOO均为边缘光线及主光线与成像面的交点。以环带像素光程长差值的二阶导数尖峰处作为测量微球轴向位移的定位标志。光程长：()()()()2222222222,cos2cossin2sinsinsin cot2sincosddyOPL z L u ynLruruLuLruuzLrLruun rn Lruu=+(1)其中()()()()21212121!sinsin2 1212!21mmmmmnunuuuLrLrnrnrmm+=+(2)环带像素光程长值为()(),1111,&,22221111,&,2222111,&22npRPOPLzOPLx y zxnxnxynynyxnxnxynynyxnxnxyn=+1,221111,&,.2222ynyxnxnxynyny+(3)由于光程长与相位呈线性关系，光程长最大值点即为重构相位中心，同时为重构强度最大值。追踪重构像中心的强度值，多项式拟合计算微球轴向位移。3.系系统与实验统与实验 3.1.测量系测量系统与测量模型统与测量模型 Figure 3.Digital holography microscope experimental system 图图 3.数字全息显微实验系统 原熙英 等 DOI:10.12677/oe.2023.132006 47 光电子 如上图 3 所示，Laser 为激光器件，其发出的光源波长为 692 nm，该光源可被中灰镜 NF 调整光强大小，后被分光棱镜 BS1 分成样本光路和参考光路，样本光通过反射镜 M1 折射后照到被测样本，后被MO1 显微物镜放大。参考光被 MO2 显微物镜放大，参考光路图与样本光路几乎相同。使用分光棱镜 BS2将参考光和样本光进行合光，后照射在 CCD 上，发生干涉后，得到数字全息图(干涉图)。像面数字全息图强度为()22*,HIx yROR OO R=+(4)其中22RO+为零级项，O 为频谱滤波时的虚像，O*为频谱滤波时的实像，R 为当物光路没有样本时的参考光。参考光的入射角 从不同的方向进入光路，使得零级项会出现在频谱的中央位置。所以全息图的实像与虚像，分别在频谱中以零级项为中心在频谱中互不交叉。对于波前重构成像的获取信息。需要用分离出来的后两项*R OO R+通过圆形窗进行逆傅里叶变换，模拟参考光并引入误差，求得含有被测生物样本相位信息的物光虚像。再现像面上的光场复振幅的公式为：()()()()()*,u x yx y R x yRx y O x y=(5)(),x y用于校正畸变因子、重构面上的光波(),I x y为强度、相位为(),x y Figure 4.Schematic of transmitted digital holographic mi-croscopy 图图 4.透射数字全息显微的原理图 因为相位和光程差呈线性关系(见图 4)，则()()()()(),seeeesex yn d x yndd x yn dnnd x y=+=(6)公式(6)中量样本折射率为 ns，样本所处环境折射率为 ne。样本厚度(),d x y，样品池的厚度 de。则利用相位与光程差的关系得被测物体像面相位与厚度的关系为()(),sex yd x ynn=。3.2.相位中心校正法提高相位型样本测量准确度相位中心校正法提高相位型样本测量准确度 透明微球是典型的相位型样本，在生物学、生物医学领域，微球作为观测生物大分子行为的标志物，其实时位移是研究关注的重点。利用像素差分光程长法能够突破传统全息同轴技术不能在聚焦像面附近测量位移的限制。像素差分光程长法追迹两像素点折射光线，计算两像素点光程差，分析光程差二阶导数变化，随着重构距离的改变，光程长差二阶导数出现尖峰，且尖峰位置随微球的轴向移动发生等距位移。如图 5 为直径为 5 m 的 SiO2微球的全息图、重构相位图和重构强度图。依据相位图可以计算两像素点差分光程长，如图 6 为微球中心像元和距离中心 7 个像素的像元的差分光程长位移标志位，即二阶导数值。随着微球的位移，二阶导数尖峰值也发生改变。原熙英 等 DOI:10.12677/oe.2023.132006 48 光电子 (a)(b)(c)Figure 5.Microspheres(a)hologram,(b)phasediagram,and(c)intensity diagram 图图 5.微球(a)全息图、(b)相位图和(c)强度图 Figure 6.Second derivative of differential optical path length 图图 6.差分光程长二阶导数 图 6 中微球中心为利用霍夫变换求取的微粒全息图的圆心。这样定义微球中心是基于光线能够垂直入射 CCD，使 CCD 中心及光轴在同一直线上。然而在实际光路中，由于搭建光路存在角度误差，光轴不能垂直于 CCD 平面，因此不能正入射微球。倾斜角度的入射带来的后果是微球相位中心与微球成像面圆心并不是同一点。而差分光程长方法利用中心像元与其他像元做差分所获取的二阶导数值，其数值对比最为显著，以此作为标志位能够使差分光程长计算微球位移的测量获得最优准确度。但基于以上分析，中心像元的定义应保证光程长最大化，即相位中心点应为中心像元。以实验验证以上推论。原熙英 等 DOI:10.12677/oe.2023.132006 49 光电子 Figure 7.Center of microsphere image and reconstructed phase center 图图 7.微球像面圆心与重构相位中心 图 7 为重构的相位图，当重构距离为 d=3.5 mm 时，在微球区域内，相位最大值像素点的坐标为(407,379)，该点相对光程长值为329.81 nm，而霍夫变换计算的像面圆心坐标为(407,379)，对应相对光程长值为329.81 nm，相对误差达到 13.1%。利用重构相位的中心作为中心像素点，通过校正中心像元的方式提升测量准确度，抑制由环境引起的测量噪声。图 8 重构相位中心与微粒圆心差分光程长二阶导数值。由图 8 可知，校正中心像元为重构相位中心后，二阶导数标志位的尖峰相对于校正前更易区分，尖峰值更高，且相对帧数与校正前并不同，两种方法重构轴向位置相距约 0.478 nm。Figure 8.Reconstruction of the second derivative of the difference optical path length between the phase center and the particle center 图图 8.重构相位中心与微粒圆心差分光程长二阶导数值 原熙英 等 DOI:10.12677/oe.2023.132006 50 光电子 利用校正中心像元的方法，测量微粒位移，位移测量的结果图如图 9 所示。而以微粒圆心为中心像元的光程长差分方法得到的测量分辨力为 5 nm，但校正中心像元后，轴向位移测量分辨力能够达到 4 nm，提升 20%，测量信噪比相对于校正前亦有提升。Figure 9.Diagram of axial displacement of the particle after center pixel correction 图图 9.中心像元校正后微球轴向位移图 4.结论结论 通过对直径为 5 m 的 SiO2微球进行实验，验证了微粒像面圆心与重构相位中心并不在同一点，这是由光轴与 CCD 面并非垂直、光路装置安装存在偏差造成的。同时通过计算微粒像面圆心与重构相位中心对应的光程长值，证实重构相位中心的光程长值高于微粒像面圆心的对应值。因此，校正中心像元为重构相位中心，能够提高差分像素光程长的相对差值。经过实验证实，中心像元校正后，微球位移测量分辨力能够提高约 20%，信噪比亦有所提升，能够将校正中心像元法用于对浑浊介质中相位型样本的测量，提升传统光程长差分方法的成像质量及测量性能。基金项目基金项目 本文由天津市技术创新引导专项基金(22YDTPJC00980)，天津市中央引导地方科技发展资金项目(22ZYCGSN00370)，天津市研究生科研创新项目(2021YJSS141)，天津农学院研究生教育教学研究与改革项目(2021-YC-54)，国家自然科学基金项目(51775381)资助。参考文献参考文献 1 Chen,H.,et al.(2011)Improved High-Force Magnetic Tweezers for Stretching and Refolding of Proteins and Short DNA.Biophysical Journal,100,517-523.https:/doi.org/10.1016/j.bpj.2010.12.3700 原熙英 等 DOI:10.12677/oe.2023.132006 51 光电子 2 Chen,D.F.,et al.(2022)Resolution and Contrast Enhancement for Lensless Digital Holographic Microscopy and Its Application in Biomedicine.Photonics,9,Article 358.https:/doi.org/10.3390/photonics9050358 3 Scharnowski,S.and Khler,C.J.(2020)Particle Image Velocimetry-Classical Operating Rules from Todays Perspec-tive.Optics and Lasers in Engineering,135,106185.https:/doi.org/10.1016/j.optlaseng.2020.106185 4 Beresh,S.J.(2021)Time-Resolved Particle Image Velocimetry.Measurement Science and Technology,32,Article 102003.https:/doi.org/10.1088/1361-6501/ac08c5 5 Flewellen,J.L.,Minoughan,S.,Garcia,I.L.and Tolar,P.(2022)Digital Holography-Based 3D Particle Localization for Single-Molecule Tweezer Techniques.Biophysical Journal,121,2538-2549.https:/doi.org/10.1016/j.bpj.2022.06.001 6 Go,T.,Kim,J.and Sang,J.L.(2020)Three-Dimensional Volumetric Monitoring of Settling Particulate Matters on a Leaf Using Digital in-Line Holographic Microscopy.Journal of Hazardous Materials,404,124116.https:/doi.org/10.1016/j.jhazmat.2020.124116 7 Kim,J.Go,T.and Sang,J.L.(2021)Volumetric Monitoring of Airborne Particulate Matter Concentration Using Smart Phone-Based Digital Holographic Microscopy and Deep Learning.Journal of Hazardous Materials,418,126351.https:/doi.org/10.1016/j.jhazmat.2021.126351 8 Guo,S.,Wang,Y.X.,et al.(2016)Resolution Enhancement Phase-Contrastimaging by Microsphere Digital Hologra-phy.Optics Communications,36,81-87.https:/doi.org/10.1016/j.optcom.2015.12.031 9 Wang,Z.B.,Guo,W.,Lin,L.,Boris,L.,Ashfaq,K.,Liu,Z.,Chen,Z.and Hong,M.(2011)Optical Virtual Imaging at 50 nm Lateral Resolution with a White-Light Nanoscope,Nature Communications,2,Article No:218.https:/doi.org/10.1038/ncomms1211 10 Zeng,Y.N.,Lei,H.,Hu,X.T.,et al.(2015)Three-Dimensional Particle Tracking by Pixel Difference Method of Opti-cal Path Length Basedon Digital Holographic Microscopy,Journal of Vacuum Science&Technology B,33,051808.https:/doi.org/10.1116/1.4929690 11 Panahi,M.,Jamali,R.,et al.(2021)3D Monitoring of the Surface Slippage Effect on Micro-Particle Sedimentation by Digital Holographic Microscopy.Scientific Reports,11,129161.12 Yuan,H.,et al.(2021)Accurate Reconstruction for the Measurement of Tilt Surfaces with Digital Holography.Optics Communications,496,127135.https:/doi.org/10.1016/j.optcom.2021.127135 13 Shangraw,M.,and Ling,H.(2021)Separating Twin Images in Digital Holographic Microscopy Using Weak Scatterers.Applied Optics,60,626-634.https:/doi.org/10.1364/AO.410167