大视场线扫描共聚焦全息显微成像

大视场线扫描共聚焦全息显微成像 Large field of view line-scanning confocal holographic microscopy 1 2 first-author 纵浩天（1996-），男，安徽宿州人，硕士研究生，2014年于安徽大学获得学士学位，主要从事数字全息成像方面的研究。Email：zht12138@mail.ustc.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444723&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444727&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444725&type= 22.01330566 32.00408936 纵浩天（1996-），男，安徽宿州人，硕士研究生，2014年于安徽大学获得学士学位，主要从事数字全息成像方面的研究。Email： zht12138@mail.ustc.edu.cn zht12138@mail.ustc.edu.cn 2 corresp E-mail： zhangyh@sibet.ac.cn E-mail： zhangyh@sibet.ac.cn 张运海（1975-），男，湖北襄阳人，教授，研究员，博士生导师， 1998年、2006年于南京航空航天大学分别获得学士、博士学位，现为江苏省医用光学重点实验室副主任，主要从事激光扫描共聚焦成像、超分辨显微光学方面的研究。Email：zhangyh@sibet.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444729&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444735&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444732&type= 22.01330566 32.00454712 张运海（1975 - ），男，湖北襄阳人，教授，研究员，博士生导师， 1998年、2006年于南京航空航天大学分别获得学士、博士学位，现为江苏省医用光学重点实验室副主任，主要从事激光扫描共聚焦成像、超分辨显微光学方面的研究。Email：zhangyh@sibet.ac.cn zhangyh@sibet.ac.cn 1 2 1 2 中国科学技术大学 生物医学工程学院，安徽 合肥 230026 中国科学技术大学 生物医学工程学院，安徽 合肥 230026 School of Biomedical Engineering， University of Science and Technology of China， Hefei 230026， China School of Biomedical Engineering， University of Science and Technology of China， Hefei 230026， China 中国科学院 苏州生物医学工程技术研究所 江苏省医学光学重点实验室， 江苏 苏州 215163 中国科学院 苏州生物医学工程技术研究所 江苏省医学光学重点实验室， 江苏 苏州 215163 Jiangsu Key Lab of Medical Optics， Suzhou Institute of Biomedical Engineering and Technology， Chinese Academy of Sciences， Suzhou 215163， China Jiangsu Key Lab of Medical Optics， Suzhou Institute of Biomedical Engineering and Technology， Chinese Academy of Sciences， Suzhou 215163， China 传统显微成像一般记录样本的强度信息，对于半透明或相位组织成像对比度较差。为实现相位组织非荧光标记成像，采用线扫描共聚焦全息成像方法，在线扫描共聚焦成像的基础上增加一路参考光，在共聚焦狭缝处形成离轴像面数字全息，通过控制样本的移动实现对样本的扫描，将获得的干涉线合成为二维全息图，通过频域滤波的方式获得振幅与相位分布，采用相邻剖面相似的特性校正环境振动引起的相位横纹，并且通过多区域扫描拼接实现大视场全息成像。对USAF1951分辨率板进行线扫描共聚焦全息成像，采用抖动校正算法，使本实验重建相位图中的抖动横纹降低了84.7%，获取3个子区域图，通过图像拼接达到1 160 μm×1 043 μm的成像视场，扫描更多的子区域可以获取更大的视场，并且对洋葱表皮细胞实现共聚焦相位成像。实验结果表明了该线扫描共聚焦全息成像方法可以实现对半透明样本的大视场相位成像，为相关仪器研制提供了指导与依据。 Traditional microscopic imaging generally records the intensity information of a sample， but the contrast is poor for translucent or phase tissue. To realize phase tissue non-fluorescence labeling imaging， a line-scanning confocal holographic imaging method was adopted. In line-scanning confocal imaging， a reference light was added to form off-axis image plane digital holography at the confocal slit. The sample was scanned by controlling the movement of the sample. The interference lines obtained were combined to form a two-dimensional hologram， and the distribution of amplitude and phase was obtained by filtering in the frequency domain. The phase striation caused by environmental vibration was rectified using the characteristics of adjacent sections， and large field-of-view holographic imaging was realized by multi-area scanning and splicing. A USAF1951 resolution plate was imaged by line scanning confocal holography， and a jitter correction algorithm was used to reduce the jitter stripes in the reconstructed phase map by 84.7% to obtain three sub-area maps. An imaging field of 1 160 μm × 1 043 μm was achieved by splicing； scanning more sub-regions can help in obtaining a larger field of view and in realizing confocal phase imaging of onion epidermal cells. The experimental results show that the line-scanning confocal holographic imaging method can realize large-field phase imaging of translucent samples and provide guidance and a basis for the development of related instruments. 数字全息 光学显微成像 线扫描共聚焦 扫描拼接 大视场 digital holography optical microscopy line-scanning confocal scanning splicing large field of view 1 引　言 1 显微镜是人类观察与探索微观世界的重要工具，它在生物医学和工业检测等领域都发挥着重要的作用。然而，大部分细胞或生物组织的光学吸收系数很小，利用传统显微镜观察时需要对样本进行染色或荧光标记处理，但这种处理往往会带来光毒性和光污染，从而对细胞的活性产生不利的影响。针对这一问题，Zernike相衬（Phase Contrast ，PC）显微镜和微分干涉（Differential Interference Contrast， DIC）显微镜应运而生，为活细胞非标记成像研究提供了有效的手段 ［ 1 1 - 2 2 1-2 ］ 。它们基于光的衍射和干涉原理将样本的相位差信息转化为振幅信息，从而实现对无染色样本的观察。随着光学显微技术的发展，相位成像逐渐成为研究热点，定量相位成像技术得到了长足的发展 ［ 3 3 - 11 11 3-11 ］ 。 20世纪50年代，Marvin Minsky发明了共聚焦显微镜，由于它能够获得高分辨率、高对比度的图像和真正的光学切片，目前已经在工业检查、生物和医学领域得到广泛应用 ［ 12 12 - 14 14 12-14 ］ 。然而，传统的共焦显微镜采用基于强度的成像方式，这意味着丢失了样品的相位信息。为了恢复相位信息，有研究者在共焦显微镜中采用了差分干涉法和电光相位调制法 ［ 15 15 ］ 。近年来，数字全息技术（Digital Holography， DH）作为一种新兴的定量相位成像方式，在光学计量、生物学和医学中得到了广泛的研究和应用 ［ 16 16 - 18 18 16-18 ］ 。DH作为一种宽场成像方式，虽然可以实现数字重聚焦，但仍不能抑制离焦光来实现真正的光学切片。 共焦全息显微镜（Confocal Holographic Microscopy，CHM）的概念是由Herring在1997年 ［ 19 19 - 20 20 19-20 ］ 提出的，用于研究微重力下的流体。CHM在成像中需要一个非常复杂的光学系统，这使得它很难应用在其他领域。2012年，Goy等 ［ 21 21 ］ 将DH合并到点扫描共聚焦显微镜（Point-scanning Confocal Microscope，PCM）中，提出了一种数字点扫描共聚焦显微镜（Digital Point-scanning Confocal Microscope，DPCM），并于2013年 ［ 22 22 ］ 采集到生物样品的相位图。2014年，Schnell等 ［ 23 23 ］ 将合成光学全息术（Synthetic Optical Holography，SOH）与点扫描共聚焦成像系统结合，通过移动参考光束实现了同轴全息，完成了合成光学全息定量共焦相位成像。2016年，Liu初步验证了线扫描共聚焦与DH结合的可行性 ［ 24 24 ］ ，在实验中用分辨率板做测试，但并未使用生物样本，实际应用中还存在问题，比如受光路系统和CCD或CMOS等光电传感器件的阵列尺寸和像素总数的限制成像视场过小，环境振动会引起相位抖动横纹等，从而限制了这种成像技术在大视场测量中的应用与推广。 本文将线扫描共聚焦与离轴DH相结合，利用位移台扫描获取多条干涉线并合成二维共聚焦全息图，采用频域滤波法提取样本的振幅与相位分布，实现线扫描共聚焦全息成像。采用多区域拼接法，在不影响分辨率的情况下扩大了线扫描共聚焦全息成像的视场；且利用相邻剖面相似的特点校正了环境、机械振动等外部因素所造成的相位抖动横纹问题，为相关仪器的研制提供了依据与指导。 2 方　法 1 线扫描共焦全息成像以DH为基础，以高相干光源激光器作为光源（目前多为He-Ne激光器），迈克尔逊干涉仪或马赫-曾德尔（Mach-Zehnder， M-Z）干涉仪作为基本光路框架，将高相干的激光分为两路，分别为探测光和参考光，利用线探测光和线参考光发生干涉从而产生一维干涉条纹。干涉条纹包含了物体的振幅与相位信息。通过狭缝滤除非焦面杂散光，用CCD作为探测器记录干涉条纹进而得到全息图。 2.1　数字全息 2 在离轴数字全息的过程中，物光波与参考光波以一定的夹角在全息面上发生干涉并被CCD记录下来。物光波的复振幅含有振幅和相位信息，但是CCD只能记录强度信息，因此使用参考光与物光波发生干涉，使物光波的相位信息转化为强度信息，即全息图，于是相位信息便被记录在干涉条纹中。后续的重建过程即从全息图的条纹信息中将振幅与相位信息提取出来并处理。 定义到达CCD记录面时物光波 <math id="M1"><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444750&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444746&type= 12.61533356 3.81000018 和参考光波 <math id="M2"><mi>R</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444754&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444752&type= 12.19199944 3.81000018 的复振幅分别为： <math display="block" id="M3"><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">[</mo><mo>-</mo><mi mathvariant="normal">j</mi><msub><mrow><mi>ϕ</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444758&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444756&type= 59.18199921 4.74133348 ， <math display="block" id="M4"><mi>R</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">[</mo><mo>-</mo><mi mathvariant="normal">j</mi><msub><mrow><mi>ϕ</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444762&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444760&type= 58.50466537 4.74133348 ， 式中： <math id="M5"><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444766&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444764&type= 13.88533401 4.74133348 和 <math id="M6"><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444770&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444768&type= 13.88533401 4.74133348 分别是物光和参考光的振幅分布， <math id="M7"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444775&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444773&type= 13.20799923 4.74133348 和 <math id="M8"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444780&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444778&type= 12.86933327 4.74133348 分别是物光和参考光的相位分布。CCD记录面上的光强可以表示为： <math display="block" id="M9"><mi>I</mi><mo>=</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>O</mi><mo>+</mo><mi>R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>=</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>O</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mi>O</mi><msup><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo>+</mo><msup><mrow><mi>O</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444782&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444781&type= 68.49533081 4.74133348 . 将 式（1） 式（1） 和 式（2） 式（2） 中两束光波的复振幅代入 式（3） 式（3） 并化简可得： <math display="block" id="M10"><mi>I</mi><mo>=</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>O</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mn mathvariant="normal">2</mn><mfenced open="|" close="|" separators="|"><mrow><mi>O</mi></mrow></mfenced><mfenced open="|" close="|" separators="|"><mrow><mi>R</mi></mrow></mfenced><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><mo>-</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444784&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444783&type= 63.33066559 4.74133348 . 式（4） 式（4） 为全息图的分布，前两项表征的是物光波和参考光波的振幅；而第三项则是物光与参考光的干涉项，在全息图中表现为干涉条纹的分布，包含物光的振幅和相位信息，并且分别受到参考光振幅和相位的调制。由于物光与参考光存在一定的夹角，引入了参考光的载频，因此干涉得到的全息图经傅里叶变换后，所得的空间频谱图的各干涉像相互分离，在全息图衍射场中有0级和±1级干涉项。在滤波重建的过程中，需要消去+1级项中参考光的载频，消去 <math id="M11"><msup><mrow><mi>R</mi></mrow><mrow><mi>*</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444791&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444789&type= 3.81000018 3.04799986 的方法可以表示为： <math display="block" id="M12"><mi>I</mi><mi>R</mi><mo>=</mo><mo stretchy="false">(</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>O</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mi>R</mi><mo>+</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>O</mi><mo>+</mo><mi>R</mi><mi>R</mi><msup><mrow><mi>O</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444795&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444793&type= 62.22999954 4.74133348 ， 式中，第一项为直透光或者零级像； 第二项为+1级项或者原始像，包含物光的原始信息； 第三项是产生畸变的共轭像，也称为 - 1级项。在频谱中将+1级项 <math id="M13"><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>R</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>O</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444801&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444798&type= 9.56733322 4.74133348 单独滤出，由于参考光一般为平面波，所以其振幅常作为常数忽略。一般情况下需要通过数值衍射的方式重建物光场，但线扫描共焦全息显微成像利用的是像面全息，逆变换到空间域后无需进行衍射重建。强度与相位分布分别为： <math display="block" id="M14"><msub><mrow><mi>A</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444809&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444807&type= 35.39066696 5.58799982 ， <math display="block" id="M15"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfenced open="[" close="]" separators="|"><mrow><mfrac><mrow><mi mathvariant="normal">I</mi><mi mathvariant="normal">m</mi><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi><mo stretchy="false">(</mo><mi>O</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444813&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444811&type= 51.64666748 9.73666668 . 2.2　线扫描共聚焦全息 2 2 . 2 . 1　光路结构 3 线扫描共焦全息显微镜是在线扫描共聚焦系统的基础上，增加一路参考光，形成迈克尔逊干涉仪框架，如 图2 图2 所示。激光经扩束准直后，经柱面镜CL单方向会聚，即一个方向为平行光，另一个方向会聚，目的是为了在最终样本面形成一条焦线，后续需要利用这条焦线扫描成像。由分光棱镜BS分光后，反射的一路为常规线扫描共聚焦光路，经物镜在样本面会聚成一条焦线，由位移台控制样本移动实现扫描（见 图4 图4 ）。透射的参考光路中采用和探测光路中相同的结构，目的是为使两路的光程差维持在光源的相干长度内，并且使参考光在CCD靶面上会聚成一条焦线，与物光的焦线发生干涉，生成共聚焦干涉线并由CCD记录。为了实现一维离轴DH，在 yoz 平面内稍微倾斜分光棱镜BS来引入一个离轴角。 2 . 2 . 2　处理过程 3 由线扫描共聚焦成像的结构可知，在探测端需要有狭缝滤除离焦光。本文采用的是虚拟狭缝，利用CCD的若干行像素宽度作为狭缝。优势为后期数据处理时可以根据具体情况自由调节狭缝的宽度。具体公式如下： <math display="block" id="M16"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">n</mi><mi mathvariant="normal">f</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>n</mi><mo stretchy="false">)</mo><mo>=</mo><munder><mo movablelimits="false">∑</mo><mrow><mi>I</mi><mo>∈</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">l</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">t</mi></mrow></munder><mrow><msub><mrow><mi>I</mi></mrow><mrow><mi>n</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444837&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444835&type= 37.50733566 6.77333355 . 式（8） 式（8） 的含义是将虚拟狭缝内的若干行像素对应相加作为当前第 n 行的像素值。虚拟狭缝的宽度为： <math display="block" id="M17"><mi>S</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1.22</mn><mi>λ</mi><mo>×</mo><mi>M</mi></mrow><mrow><mi>P</mi><mo>×</mo><mi>N</mi><mi>A</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444842&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444840&type= 24.63800049 7.62000036 ， 式中： λ 为光源波长， M 为系统放大倍数， NA 为物镜的数值孔径， P 为CCD的像素尺寸。 线探测光与线参考光干涉生成一维干涉线并由CCD记录，随着位移台的移动扫描，CCD再记录下一个位置的一维干涉线，以此类推。相对于针对每一条线都进行全息重建，分别提取振幅与相位然后生成二维图 ［ 24 24 ］ ，本文采用 式（8） 式（8） 和 式（9） 式（9） 处理每条一维干涉线，然后合成二维全息图，后续只需进行一次全息重建便可分别获得振幅与相位图，处理过程更加简洁高效，处理流程如 图3 图3 所示。 2 . 2 . 3　抖动校正 3 在位移台扫描的过程中，由于干涉条纹很敏感，环境振动、机械振动可能会引起线上干涉条纹的整体抖动，导致不同线的条纹分布会出现错位。而相位信息正包含在条纹分布中，这种错位会引起相位的偏移。由于是线扫描，所以对于每一行的像素，其相位偏移值是相同的，利用相邻行剖面相似的特性，采用下述步骤可以很大程度上校正抖动： （1）用第2行像素的相位减去第1行对应像素的相位，然后对得到的结果取平均，得到第2行相对于第1行的相位偏移量 <math id="M18"><mo>∆</mo><msub><mrow><mi>h</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444846&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444844&type= 5.92666674 4.06400013 ； （2）第2行所有像素相位值都减去相位偏移量 <math id="M19"><mo>∆</mo><msub><mrow><mi>h</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444853&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444850&type= 5.92666674 4.06400013 ； （3）用第 n +1行像素的相位减去第 n 行对应像素的相位，然后对得到的结果取平均，得到第 n +1行相对于第 n 行的相位偏移量 <math id="M20"><mo>∆</mo><msub><mrow><mi>h</mi></mrow><mrow><mi>n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444862&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444855&type= 5.92666674 4.06400013 ； （4）第 n +1行减去相位偏移量 <math id="M21"><mo>∆</mo><msub><mrow><mi>h</mi></mrow><mrow><mi>n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444868&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=9444865&type= 5.92666674 4.06400013 ；以此类推，直至最后一行。 这是一种较为近似的算法，能很大程度上去除抖动所产生的横纹，对于部分阶跃型样本可能存在较小的误差，但在绝大多数情况下可以达到良好的效果。 2.3　大视场扫描拼接方法 2 受到线扫描共焦成像系统与CCD靶面尺寸的限制，无法同时满足分辨率与视野的实际应用需求。为了解决这一问题，本文采用子图拼接的方式，在保证CCD采样分辨率的情况下获取更大的成像视场。采用位移台移动扫描的方法实现多幅子图的采集，对采集的子图拼接，达到大视场成像的目的，并且结构相对简单。 在线扫描共焦全息显微系统中，光源发出的光束通过系统后在样本面形成一条焦线，利用位移台控制样本移动实现焦线对样本的扫描，如 图4 图4 所示。根据重叠区域要求和单次成像区域大小确定位移台沿 X 方向的位移量，位移台每移动一个位置，控制位移台沿 Y 方向移动扫描并采集一幅图像。相比于其他拼接方式 ［ 16 16 ］ ，本文为了充分利用线扫描成像的优势，使后续拼接过程更加高效与便捷，利用位移台沿 Y 方向连续扫描出一矩形区域，使用矩形区域拼接的方式，过程如下： （1）控制位移台移动到某一待成像位置后，然后沿 Y 方向扫描一矩形区域，命名为 A 1 ，位移台复位； （2）控制位移台沿 X 方向平移一固定步进，以步骤1中相同位移量控制位移台扫描一矩形区域，命名为 A 2 ，位移台复位； 以此类推，至最后一个子区域 A n ，后续对 A 1 ~ A n 做拼接处理。为了避免信息丢失，子区域之间需要有一定的重叠区域，重叠区域的大小可以根据沿 X 轴的视场，通过控制位移台的步进调节。 常规拼接算法通常采用基于特征点的图像配准方法去确定图像之间的重叠区域，然后再进行拼接处理。但本文中重叠区域的大小可根据位移台横向步进值和子区域横向视场的大小计算得到，所以直接进行拼接效率更高，拼接算法流程如 图5 图5 所示。 根据图像重叠部分的大小，将具有含有重叠区域的子区域分为非重叠区域和重叠区域，例如A与B、C与D。其中B与C为重叠区域，两者叠加并平均可得重叠部分E图，随后将三部分A，E，D拼合，便可完成子区域拼接。多幅子区域图像之间的拼接与上述相同。重叠区域的大小是根据位移台横向步进值和子区域横向视场计算得到的，因此拼接精度取决于位移台的位移精度。 3 实验结果 1 线扫描共焦全息显微镜系统原理如 图2 图2 所示。光源选用He-Ne激光器，波长为632.8 nm；显微物镜放大倍数为10；探测器采用面阵CMOS（滨松flash 4.0 V2），像素数为2 048×2 048，像素尺寸为6.5 μm×6.5 μm，靶面尺寸约为13 mm，具有自定义工作区域的功能，本实验采用的虚拟狭缝就是利用这一功能，工作在条状区域下的采集速度可达7 894 frame/s。柱面镜CL的焦距 f =100 mm，透镜L4的焦距 f =100 mm，筒镜为OLYMPUS X71显微镜内置，焦距为180 mm。L3的焦距 f =300 mm，可以算得系统的总放大倍数为30。 局部狭缝内线干涉线如 图6 图6 所示。根据实验系统的相关参数，由 式（9） 式（9） 算得狭缝宽度约为8.9行像素，实验中取整数8行，将面阵CMOS工作区域设置为2 048×8像素。首先对分辨率板进行成像，由以上过程合成整幅线扫描二维共聚焦全息图，如 图7 图7 所示。 以 图7 图7 （a）中的子区域 A 1 为例，利用离轴全息重建算法，首先利用傅里叶变换得到其频谱，可以看到 - 1级项、0级项、+1级项依次排布。根据全息重建原理，需要将+1级项单独滤出并频移到频谱中心，利用傅里叶逆变换到空间域便可得到物光波复振幅，可以分别根据 式（6） 式（6） 、 式（7） 式（7） 提取物光的振幅与相位，如 图8 图8 所示。 实验中采用减去系统空采相位图的方法去除系统的相位畸变。用离轴全息重建算法分别对子区域 A 1 ， A 2 ， A 3 进行处理，提取出振幅和相位信息。其中相位图中的横纹是由于在位移台扫描过程中，环境、机械振动引起的条纹抖动，导致最终重建出的线扫描共聚焦全息图中的条纹产生了横向上的微小错位，使最终重建的相位图中产生横纹，符合前期的预测。用第二节中的方法来校正抖动所带来的横纹，效果如 图8 图8 （d）所示。 图9 图9 为 图8 图8 （c）中沿白线和8（d）中对应位置的剖面对比，校正前后对比效果明显，计算方差可知本次实验中抖动程度降低了约84.7%。 对子区域重建后的图像进行拼接处理，拼得大视场图像如 图10 图10 所示。拼接所得的视场达1 160 μm×1 043 μm，本文采用线扫描共聚焦成像方式，所以根据图像线扫描方向可分为共聚焦方向和非共聚焦方向，理论上共聚焦方向和非共聚焦方向有着不同的成像分辨率，共聚焦方向的成像分辨率要略高于非共聚焦方向 ［ 25 25 ］ 。根据 图10 图10 （d）可以看出，分辨率板第8组第6条可以分辨清楚，共聚焦方向分辨率高于2.19 μm，非共聚焦方向可分辨第8组第4条，非共聚焦方向分辨率高于2.73 μm。由此可见，线扫描共聚焦全息成像的分辨率符合传统线扫描共聚焦显微镜的分辨率的规律。通过在物镜焦面附近沿 Z 轴方向按相同步进值移动反射镜，测量其成像总强度，然后绘制可得其归一化强度曲线，如 图10 图10 （e）所示。测两者的半峰全宽（Full Width Half Maximum，FWHM）可知，线扫描共聚焦全息显微镜的轴向分辨率与传统线扫描共聚焦显微镜相当。 图11 图11 为扫描洋葱表皮细胞表面的结果。常见的动植物细胞直径约为5~100 μm，理论上本成像系统的分辨率足以对一般的生物细胞成像。这里以洋葱表皮细胞为例，图中细胞近似长方形，宽度约为60~100 μm，长度约为150~300 μm。 图11 图11 （a）为线扫描共焦扫描并拼接所得的洋葱细胞成像图； 图11 图11 （b）为线扫描共焦全息重建并拼接所得的洋葱细胞成像图； 图11 图11 （c）为全息重建算法提取出的相位信息，利用最小二乘解包裹算法 ［ 26 26 ］ 解包裹处理后拼接所得的相位成像图。 4 结　论 1 本文采用线扫描共聚焦与离轴DH相结合的方法，将获得的共聚焦干涉线合成二维共聚焦全息图，采用频域滤波法成功提取出分辨率板和洋葱表皮细胞的相位信息，实现了基于共聚焦的定量相位成像。采用相邻剖面相似的特性校正了由于环境或机械振动所造成的相位横纹，使本实验中重建所得相位图的抖动程度降低了84.7%。并且提出了一种针对线扫描共聚焦全息成像的大视场成像系统，采用控制位移台移动样本扫描的方法达到对样本的大范围线扫描的目的，实现了对大视场子区域的采集与拼接，达到了大视场成像的目的。该系统通过控制位移台扫描3个矩形并且有一定重叠区域的子区域进行拼接，在不影响成像分辨率的前提下获得了1 160 μm×1 043 μm的大视场，并且采集更多子区域可以获得更大的视场。它具有结构简单、易于应用的优势，可为相关仪器的研制提供参考。 参考文献 ZERNIKE F . How I discovered phase contrast ［J］. Science ， 1955 ， 121 （ 3141 ）： 345 - 349 . How I discovered phase contrast BURCH C R ， STOCK J P . Phase-contrast microscopy ［J］. Journal of Scientific Instruments ， 1942 ， 19 （ 5 ）： 71 - 75 . Phase-contrast microscopy 李建素 ， 党长营 . 预放大数字全息系统记录距离的设计 ［J］. 光学 精密工程 ， 2019 ， 27 （ 4 ）： 784 - 794 . 预放大数字全息系统记录距离的设计 LI J S ， DANG CH Y . Design of recording distance in pre-magnification digital holography ［J］. Optics and Precision Engineering ， 2019 ， 27 （ 4 ）： 784 - 794 . （in Chinese） Design of recording distance in pre-magnification digital holography LUO Z P ， MA J S ， SU P ， et al . . Digital holographic phase imaging based on phase iteratively enhanced compressive sensing ［J］. Optics Letters ， 2019 ， 44 （ 6 ）： 1395 - 1398 . . Digital holographic phase imaging based on phase iteratively enhanced compressive sensing 张璐 ， 汤其剑 ， 邓定南 ， 等 . 基于光强传输方程的双相机动态相位成像的视场校正 ［J］. 中国激光 ， 2019 ， 46 （ 8 ）： 164 - 169 . 基于光强传输方程的双相机动态相位成像的视场校正 ZHANG L ， TANG Q J ， DENG D N ， et al . . Field-of-view correction for dual-camera dynamic phase imaging based on transport of intensity equation ［J］. Chinese Journal of Lasers ， 2019 ， 46 （ 8 ）： 164 - 169 . （in Chinese） . Field-of-view correction for dual-camera dynamic phase imaging based on transport of intensity equation KIM M K . Wide area quantitative phase microscopy by spatial phase scanning digital holography ［J］. Optics Letters ， 2020 ， 45 （ 3 ）： 784 - 786 . Wide area quantitative phase microscopy by spatial phase scanning digital holography CHENG H ， DENG H L ， SHEN CH ， et al . . Phase retrieval based on transport of intensity equation and image interpolation ［J］. Infrared and Laser Engineering ， 2018 ， 47 （ 10 ）： 1026003 . . Phase retrieval based on transport of intensity equation and image interpolation 范瑶 ， 陈钱 ， 孙佳嵩 ， 等 . 差分相衬显微成像技术发展综述 ［J］. 红外与激光工程 ， 2019 ， 48 （ 6 ）： 0603014 . 差分相衬显微成像技术发展综述 FAN Y ， CHEN Q ， SUN J S ， et al . . Review of the development of differential phase contrast microscopy ［J］. Infrared and Laser Engineering ， 2019 ， 48 （ 6 ）： 0603014 . （in Chinese） . Review of the development of differential phase contrast microscopy 李通 ， 赵巨峰 ， 毛海锋 ， 等 . 基于角度照明优化的傅里叶叠层显微成像方法 ［J］. 激光与光电子学进展 ， 2020 ， 57 （ 8 ）： 276 - 285 . 基于角度照明优化的傅里叶叠层显微成像方法 LI T ， ZHAO J F ， MAO H F ， et al . . An efficient Fourier ptychographic microscopy imaging method based on angle illumination optimization ［J］. Laser & Optoelectronics Progress ， 2020 ， 57 （ 8 ）： 276 - 285 . （in Chinese） . An efficient Fourier ptychographic microscopy imaging method based on angle illumination optimization SUN M L ， CHEN X ， ZHU Y Q ， et al . . Neural network model combined with pupil recovery for Fourier ptychographic microscopy ［J］. Optics Express ， 2019 ， 27 （ 17 ）： 24161 - 24174 . . Neural network model combined with pupil recovery for Fourier ptychographic microscopy CLAVEAU R ， MANESCU P ， ELMI M ， et al . . Digital refocusing and extended depth of field reconstruction in Fourier ptychographic microscopy ［J］. Biomedical Optics Express ， 2020 ， 11 （ 1 ）： 215 - 226 . . Digital refocusing and extended depth of field reconstruction in Fourier ptychographic microscopy 张运海 ， 杨皓旻 ， 孔晨晖 . 激光扫描共聚焦光谱成像系统 ［J］. 光学 精密工程 ， 2014 ， 22 （ 6 ）： 1446 - 1453 . 激光扫描共聚焦光谱成像系统 ZHANG Y H ， YANG H M ， KONG CH H . Spectral imaging system on laser scanning confocal microscopy ［J］. Optics and Precision Engineering ， 2014 ， 22 （ 6 ）： 1446 - 1453 . （in Chinese） Spectral imaging system on laser scanning confocal microscopy 缪新 ， 张运海 ， 黄维 . 皮肤反射式共聚焦显微成像自适应图像亮度调节 ［J］. 光学 精密工程 ， 2019 ， 27 （ 6 ）： 1270 - 1276 . 皮肤反射式共聚焦显微成像自适应图像亮度调节 MIAO X ， ZHANG Y H ， HUANG W . Image brightness adaptive adjustment during skin imaging by reflectance confocal microscopy ［J］. Optics and Precision Engineering ， 2019 ， 27 （ 6 ）： 1270 - 1276 . （in Chinese） Image brightness adaptive adjustment during skin imaging by reflectance confocal microscopy 叶寒 ， 翁祖昕 ， 张运海 ， 等 . 修正边界区域误差的共聚焦显微成像粗糙度测量 ［J］. 激光与光电子学进展 ， 2020 ， 57 （ 21 ）： 211203 . 修正边界区域误差的共聚焦显微成像粗糙度测量 YE H ， WENG Z X ， ZHANG Y H ， et al . . Surface roughness measurement using laser confocal microscope with boundary area correction ［J］. Laser & Optoelectronics Progress ， 2020 ， 57 （ 21 ）： 211203 . （in Chinese） . Surface roughness measurement using laser confocal microscope with boundary area correction CORLE T R ， KINO G S . Introduction ［M］// Confocal Scanning Optical Microscopy and Related Imaging Systems . Amsterdam ： Elsevier ， 1996 ： 1 - 66 . Introduction 冯方 ， 田爱玲 ， 刘丙才 ， 等 . 基于数字全息扫描成像的划痕缺陷全场三维测试 ［J］. 中国激光 ， 2020 ， 47 （ 4 ）： 261 - 268 . 基于数字全息扫描成像的划痕缺陷全场三维测试 FENG F ， TIAN A L ， LIU B C ， et al . . Full-field three-dimensional test for scratch defects using digital holographic scanning imaging system ［J］. Chinese Journal of Lasers ， 2020 ， 47 （ 4 ）： 261 - 268 . （in Chinese） . Full-field three-dimensional test for scratch defects using digital holographic scanning imaging system WITTKOPP J M ， KHOO T C ， CARNEY S ， et al . . Comparative phase imaging of live cells by digital holographic microscopy and transport of intensity equation methods ［J］. Optics Express ， 2020 ， 28 （ 5 ）： 6123 - 6133 . . Comparative phase imaging of live cells by digital holographic microscopy and transport of intensity equation methods 肖文 ， 李解 ， 潘锋 ， 等 . 基于USENet实现数字全息细胞再现相位像超分辨重构 ［J］. 光子学报 ， 2020 ， 49 （ 6 ）： 179 - 190 . 基于USENet实现数字全息细胞再现相位像超分辨重构 XIAO W ， LI J ， PAN F ， et al . . Super-resolution in digital holographic phase cell image based on USENet ［J］. Acta Photonica Sinica ， 2020 ， 49 （ 6 ）： 179 - 190 . （in Chinese） . Super-resolution in digital holographic phase cell image based on USENet HERRING R . A.Confocal scanning LASER holography and associated microscope ［J］. Optik ， 1997 ， 105 ： 65 - 68 . A.Confocal scanning LASER holography and associated microscope HERRING R A ， JACQUEMIN P ， SAWICKA B D ， et al . . Developing a confocal acoustic holography microscope for non-invasive 3D temperature and composition measurements ［J］. Ultramicroscopy ， 2009 ， 109 （ 7 ）： 830 - 836 . . Developing a confocal acoustic holography microscope for non-invasive 3D temperature and composition measurements GOY A S ， PSALTIS D . Digital confocal microscope ［J］. Optics Express ， 2012 ， 20 （ 20 ）： 22720 . Digital confocal microscope GOY A S ， UNSER M ， PSALTIS D . Multiple contrast metrics from the measurements of a digital confocal microscope ［J］. Biomedical Optics Express ， 2013 ， 4 （ 7 ）： 1091 - 1103 . Multiple contrast metrics from the measurements of a digital confocal microscope SCHNELL M ， PEREZ-ROLDAN M J ， CARNEY P S ， et al . . Quantitative confocal phase imaging by synthetic optical holography ［J］. Optics Express ， 2014 ， 22 （ 12 ）： 15267 - 15276 . . Quantitative confocal phase imaging by synthetic optical holography LIU C G ， KNITTER S ， CONG Z L ， et al . . High-speed line-field confocal holographic microscope for quantitative phase imaging ［J］. Optics Express ， 2016 ， 24 （ 9 ）： 9251 - 9265 . . High-speed line-field confocal holographic microscope for quantitative phase imaging 赵家旺 ， 张运海 ， 王发民 ， 等 . 线扫描虚拟结构调制共聚焦显微成像 ［J］. 中国光学 ， 2020 ， 13 （ 6 ）： 1 - 9 . 线扫描虚拟结构调制共聚焦显微成像 ZHAO J W ， ZHANG Y H ， WANG F M ， et al . Line-scanning confocal microscopy through virtual structured modulation ［J］. Chinese Optics ， 2020 ， 13 （ 6 ）： 1 - 9 . （in Chinese） Line-scanning confocal microscopy through virtual structured modulation 钱晓凡 ， 饶帆 ， 李兴华 ， 等 . 精确最小二乘相位解包裹算法 ［J］. 中国激光 ， 2012 ， 39 （ 2 ）： 0209001 . 精确最小二乘相位解包裹算法 QIAN X F ， RAO F ， LI X H ， et al . . Accurate least-squares phase unwrapping algorithm ［J］. Chinese Journal of Lasers ， 2012 ， 39 （ 2 ）： 0209001 . （in Chinese） . Accurate least-squares phase unwrapping algorithm 10.37188/OPE.20212901.0001 O426.1 O438.1 A 1004-924X（2021）01-0001-09 纵浩天,张运海,王发民等.大视场线扫描共聚焦全息显微成像[J].光学精密工程,2021,29(01):1-9. 纵浩天,张运海,王发民等.大视场线扫描共聚焦全息显微成像[J].光学精密工程,2021,29(01):1-9. ZONG Hao-tian,ZHANG Yun-hai,WANG Fa-min,et al.Large field of view line-scanning confocal holographic microscopy[J].Optics and Precision Engineering,2021,29(01):1-9. ZONG Hao-tian,ZHANG Yun-hai,WANG Fa-min,et al.Large field of view line-scanning confocal holographic microscopy[J].Optics and Precision Engineering,2021,29(01):1-9. 江苏省国际合作项目（No.BZ2020004）；国家重点研发计划资助项目（No.2017YFC0110303）；济南市“高校20条”资助项目（No.2018GXRC018）；山东省自然科学基金资助项目（No.ZR2019BF012） 2020-08-28 2020-10-21 2021-01-15 2021-03-25 光学精密工程 1 29