CN113009801B - 高速多方向线共焦数字全息三维显微成像方法及装置 - Google Patents
高速多方向线共焦数字全息三维显微成像方法及装置 Download PDFInfo
- Publication number
- CN113009801B CN113009801B CN202110218940.8A CN202110218940A CN113009801B CN 113009801 B CN113009801 B CN 113009801B CN 202110218940 A CN202110218940 A CN 202110218940A CN 113009801 B CN113009801 B CN 113009801B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- lens
- light signal
- speed
- digital holographic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 69
- 238000001228 spectrum Methods 0.000 claims abstract description 59
- 238000001514 detection method Methods 0.000 claims abstract description 54
- 230000010287 polarization Effects 0.000 claims abstract description 35
- 230000003287 optical effect Effects 0.000 claims abstract description 19
- 238000006243 chemical reaction Methods 0.000 claims description 23
- 230000005540 biological transmission Effects 0.000 claims description 20
- 238000001914 filtration Methods 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 15
- 238000012545 processing Methods 0.000 claims description 15
- 230000003595 spectral effect Effects 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 7
- 238000013519 translation Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 5
- 101100173586 Schizosaccharomyces pombe (strain 972 / ATCC 24843) fft2 gene Proteins 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 3
- 230000007123 defense Effects 0.000 claims description 2
- 239000000523 sample Substances 0.000 description 37
- 238000005516 engineering process Methods 0.000 description 13
- 238000001093 holography Methods 0.000 description 5
- 238000005286 illumination Methods 0.000 description 5
- 238000010226 confocal imaging Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004624 confocal microscopy Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/0005—Adaptation of holography to specific applications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0808—Methods of numerical synthesis, e.g. coherent ray tracing [CRT], diffraction specific
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/0005—Adaptation of holography to specific applications
- G03H2001/0033—Adaptation of holography to specific applications in hologrammetry for measuring or analysing
- G03H2001/0038—Adaptation of holography to specific applications in hologrammetry for measuring or analysing analogue or digital holobjects
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/0005—Adaptation of holography to specific applications
- G03H2001/005—Adaptation of holography to specific applications in microscopy, e.g. digital holographic microscope [DHM]
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/0443—Digital holography, i.e. recording holograms with digital recording means
- G03H2001/0445—Off-axis recording arrangement
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Microscoopes, Condenser (AREA)
Abstract
本发明公开了一种高速多方向线共焦数字全息三维显微成像方法及装置，光信号经起偏、扩束和聚焦后分为参考光信号和探测光信号；将参考光信号中继至全反射镜并接收参考返回光信号；将探测光信号中继至样品平台并接收探测扫描光信号；参考返回光信号和探测扫描光信号合束后滤除杂散光，再中继到高速线阵相机的光敏面，得到多帧一维线共焦干涉图；多帧一维线共焦干涉图经过采集和叠加生成多个扫描方向的重建干涉场，再进行处理实现样品平台的物光频谱拓展，然后对拓展后的物光频谱做处理，重建物光波场的三维相位分布。本发明可以有效地实现高速多方向线共焦数字全息显微成像，具有分辨率高的优点。
Description
技术领域
本发明涉及三维显微成像技术领域，特别涉及高速多方向线共焦数字全息三维显微成像方法及装置。
背景技术
数字全息三维成像技术沿袭了光学全息记录理论，利用光电探测器代替普通全息记录介质记录干涉图，通过计算机数值计算取代光学衍射来再现所记录的物光波场，同时获得物光波场的振幅信息和相位信息，整个记录和再现过程都数字化，是一种无接触、无破坏且具有较高分辨率的测量介质中相位变化的3D成像技术，有望成为3D微纳成像领域研究的有力工具。
空间分辨率是评价数字全息三维显微成像***的重要指标，由成像***的数值孔径(NA)和记录波长共同确定。在数值孔径一定时，采用深紫外波长记录光也可以有效提高记录***的分辨率，但是对波长要求较高，同时记录光可能会对样品产生影响；在波长一定时，使用大数值孔径的成像物镜可以有效提高记录***的分辨率，但是大数值孔径的成像物镜却有着工作距离和焦深短以及视场小的缺点，而且成像物镜的像差对再现质量造成的影响也难以消除。另外还受到记录器件(如CCD,CMOS)有限的光敏面尺寸及像元尺寸的限制，这在很大程度上影响了数字全息三维显微技术的进一步发展及应用。
为了突破数字全息中成像***空间带宽积的限制，研究人员提出了光栅调制数字全息、合成孔径数字全息、结构光照明数字全息、倾斜照明数字全息等方法。
从空间频域角度出发，再现像的横向分辨率由数字全息显微***所能记录下的物体衍射的最高空间频率决定。光栅调制数字全息三维成像技术通过在物体和线阵相机光敏面之间放置衍射光栅，记录过程中利用光栅的衍射将物体衍射的高频光波再次衍射到线阵相机光敏面的靶面上，压缩物体的空间频率，在数值再现过程中加入“数字光栅”，实现横向超分辨的再现成像。但是尚不能应用于高速的数字全息成像。
合成孔径数字全息三维成像技术是一种不需要在记录***中使用显微放大***即可获得高分辨率再现像的技术，能突破由数值孔径决定的***分辨率极限，但是由于传统方法需要移动成像器件或者被测物体，在记录过程中会耗费大量的时间，而且对***的稳定性要求也会很高。
结构光照明数字全息三维成像技术，采用类似于莫尔条纹原理，通过对宽场显微镜照明的预先余弦调制，拓展接收***频域带宽，提高截止频率，能够使***横向分辨率提高到原来的2倍。但是在记录过程中，需要在2π范围内准确移动调制光的相位获得多幅图像，将获取的各幅干涉图像进行快速傅里叶变换。同时为了实现分辨率的各向同性，需要在多个方向上对称地采用结构光照明。在提高成像分辨率的同时却增加了物体记录与再现方面的复杂性。
上述数字全息研究方法由于采用宽场成像，纵向分辨力低，并且噪声影响严重。共焦显微技术由于具备了良好的纵向分辨力而受到了广泛应用，并且横向分辨力也可以得到提高。在先技术(A.S.Goy and D.Psaltis,Digital confocal microscope,Opt.Express2012,20(20):22720-22727)将点扫描共聚焦显微技术与数字全息三维成像技术相结合，对扫描点进行实时数字全息记录，获得目标的所有相位和强度信息，可以重建目标的三维图像。为了进一步提高成像速度，在先技术(Changgeng Liu,Sebastian Knitter,ZhilongCong,Ikbal Sencan,Hui Cao,and Michael A.Choma,High-speed line-field confocalholographic microscope for quantitative phase imaging,Opt.Express 2016,24(9):9251-9265)将一维线共焦显微技术与数字全息三维成像技术相结合，集成两种技术的优点，通过逐线记录一维离轴干涉图，实现快速高分辨低噪声的数字全息三维成像。但是在先技术均只能实现一维方向线共焦，无法实现二维横向分辨率的各向同性。
发明内容
本发明的目的在于，提供高速多方向线共焦数字全息三维显微成像方法及装置。本发明可以有效地实现高速多方向线共焦数字全息显微成像；此外，本发明提高了三维成像的横向和纵向空间分辨率，实现高分辨率的数字全息三维成像，同时还实现了横向分辨率各向同性，具有信噪比高、抗干扰能力强的优点。
本发明的技术方案：高速多方向线共焦数字全息三维显微成像方法，其特征在于：激光信号首先经过起偏、扩束和聚焦后，由偏振分束棱镜分为参考光信号和探测光信号；将参考光信号通过第一透镜组中继至全反射镜并接收参考返回光信号；将探测光信号通过第二透镜组中继至样品平台并接收探测扫描光信号，其中探测光信号在中继过程中经过一维快速扫描器和线扫描方向快速转换单元，利用一维快速扫描器实现对样品平台的线扫描，利用线扫描方向快速转换单元实现多个扫描方向的快速切换，完成对样品的多方向共焦线扫描；参考返回光信号和探测扫描光信号在偏振分束棱镜合束进入接收光路；在接收光路中，合束光信号先经过狭缝滤除非聚焦平面的杂散光，再通过检偏形成离轴干涉场，然后通过接收透镜组中继到高速线阵相机，并得到多帧一维线共焦干涉图，经过高速数字采集后，在主控计算机中进行叠加生成多方向重建干涉场，并通过频谱叠加和角谱后向传输处理，实现高速多方向线共焦数字全息显微成像。
上述的高速多方向线共焦数字全息三维显微成像方法，在主控计算机中进行叠加生成多方向重建干涉场后，先进行二维傅里叶变换，再进行频谱叠加、滤波和平均处理，实现样品平台的物光频谱拓展，然后对拓展后的物光频谱做角谱后向传输处理，最后采用解包裹运算重建物光波场的三维相位分布，完成三维成像，同时通过样品平台轴向平移获得样品不同深度层次的三维相位信息。
前述的高速多方向线共焦数字全息三维显微成像方法，所述第一透镜组包括沿参考光信号传输方向依次设置的第一透镜、第一四分之一波片和第二透镜；所述第一四分之一波片用于调整参考返回光信号的光强。
前述的高速多方向线共焦数字全息三维显微成像方法，所述第二透镜组包括沿探测光信号传输方向依次设置的第三透镜、一维快速扫描器、第四透镜、线扫描方向快速转换单元、第五透镜和显微物镜；所述一维快速扫描器处于第三透镜和第四透镜之间；所述线扫描方向快速转换单元处于第四透镜和第五透镜之间。
前述的高速多方向线共焦数字全息三维显微成像方法，所述线扫描方向快速转换单元由第一二维快速扫描器、中空正多边形内反射镜和第二二维快速扫描器组成，入射的探测光信号首先经过第一二维快速扫描器偏转，然后发射到中空正多边形内反射镜的单个反射镜面，再反射到第二二维快速扫描器，由第二二维快速扫描器完成角度补偿后出射，出射的探测光信号相对于入射的探测光信号实现二维扫描视场内一个固定角度的旋转。
前述的高速多方向线共焦数字全息三维显微成像方法，所述接收透镜组包括沿参考返回光信号和探测扫描光信号传输方向设置的第六透镜、第七透镜和检偏器；所述的狭缝位于第六透镜的前焦面，用于滤除非聚焦平面的杂散光。
前述的高速多方向线共焦数字全息三维显微成像方法，所述的重建干涉场的强度分布表示为：
iH(kΔxH，lΔyH)
＝a(kΔxH，lΔyH)+c(kΔxH，lΔyH)exp[-j2π(ξxΔxH+ξyΔyH)]
+c*(kΔxH，lΔyH)exp[j2π(ξxΔxH+ξyΔyH)]；
式中：
其中O0(kΔxH，lΔyH)为探测扫描光信号的离散振幅分布，
为探测扫描光信号的离散相位分布，r0(kΔxH，lΔyH)为参考返回光信号的离散振幅分布，ξx和ξy分别是x和y方向的空间频率，由参考返回光信号和探测扫描光信号的离轴夹角θ决定；k和l都是整数，-M/2≤k≤M/2，-N/2≤l≤N/2，M和N分别为两个方向的像素数；ΔxH和ΔyH为重建干涉图平面的采样间隔，ΔxH＝ΔyH＝Lx/M＝Ly/N，Lx×Ly是重建干涉场的面积，*表示卷积运算。
前述的高速多方向线共焦数字全息三维显微成像方法，将M个扫描方向的重建干涉场在主控计算机中通过二维傅里叶变换分别求出频谱，将M个频谱叠加、滤波和加权平均处理，实现样品平台的物光频谱拓展，其中二维傅里叶变换计算如下：
IH(pΔξ，qΔη)＝FFT2{iH(kΔxH，lΔyH)}p，q
＝A(pΔξ，qΔη)+C+1(pΔξ-ξx，qΔη-ξy)+C-1(pΔξ+ξx，qΔη+ξy)；
式中：IH表示iH的二维傅里叶变换，p和q都是整数，-M/2≤p≤M/2，-N/2≤q≤N/2；Δξ和Δη分别为频谱域内两个轴的取样间隔，Δξ＝1/(MΔxH)＝1/Lx，Δη＝1/(NΔyH)＝1/Ly；A为干涉频谱的零级峰，C+1为干涉频谱的+1级峰，C-1为干涉频谱-1级峰；
通过参考返回光信号和探测扫描光信号的离轴夹角θ调整ξx和ξy，使其A、C+1和C-1在频谱域内足够分开，然后将频谱域内的A和C-1滤除，再将C+1平移至频谱域的中心，再进行角谱后向传输处理，即利用角谱传递函数将C+1向后传播以得到样品平台所在平面的物光频谱；最后采用解包裹运算，即对物光频谱作二维逆傅里叶变换得到重建后的物光波场的复振幅分布：
式中：n和n都是整数，-M/2≤m≤M/2；-N/2≤n≤N/2；λ是波长，H(pΔξ，qΔη)是频谱域内的角谱传递函数；IFFT2是二维逆傅里叶变换；
和
为成像面的采样间隔，大小与高速线阵相机光敏面像素尺寸相同；
从而得到物光波场的相位分布：
实现如前述的高速多方向线共焦数字全息三维显微成像方法的装置，包括单频单模连续激光器，单频单模连续激光器经起偏器连接有扩束器，扩束器的输出端依次连接有半波片、矩形光阑、柱面透镜和偏振分束棱镜；所述偏振分束棱镜连接有第一透镜组、第二透镜组和接收透镜组；所述第一透镜组包括与偏振分束棱镜连接的第一透镜，第一透镜经第一四分之一波片连接有第二透镜；所述第二透镜连接有全反射镜；所述第二透镜组包括与偏振分束棱镜连接的第二四分之一波片，第二四分之一波片连接有第三透镜，第三透镜经一维快速扫描器连接有第四透镜，第四透镜经线扫描方向快速转换单元连接有第五透镜；所述第五透镜经显微物镜连接有样品平台；所述接收透镜组的前侧设有狭缝，所述接收透镜组包括与位于狭缝后侧的第六透镜，第六透镜经第七透镜连接有检偏器；所述检偏器连接有高速线阵相机，高速线阵相机经高速数字采集卡连接有主机计算机；所述主控计算机分别还与一维快速扫描器和线扫描防线快速转换单元连接。
前述的高速多方向线共焦数字全息三维显微成像方法的装置，所述线扫描方向快速转换单元由依次连接的第一二维快速扫描器、中空正多边形内反射镜和第二二维快速扫描器组成。
与现有技术相比，本发明中激光信号首先经过起偏、扩束和聚焦后，由偏振分束棱镜分为参考光信号和探测光信号；将参考光信号通过第一透镜组中继至全反射镜并接收参考返回光信号；将探测光信号通过第二透镜组中继至样品平台并接收探测扫描光信号，其中探测光信号在中继过程中经过一维快速扫描器和线扫描方向快速转换单元，利用一维快速扫描器实现对样品平台的线扫描，利用线扫描方向快速转换单元实现多个扫描方向的快速切换，完成对样品的多方向共焦线扫描；参考返回光信号和探测扫描光信号在偏振分束棱镜合束进入接收光路；在接收光路中，合束光信号先经过狭缝滤除非聚焦平面的杂散光，再通过检偏形成离轴干涉场，然后通过接收透镜组中继到高速线阵相机，并得到多帧一维线共焦干涉图，经过高速数字采集后，在主控计算机中进行叠加生成多方向重建干涉场，并通过频谱叠加和角谱后向传输处理，由此本发明可以实现高速多方向线共焦数字全息显微成像。此外，本发明在主控计算机中进行叠加生成多方向重建干涉场后，先进行二维傅里叶变换，再进行频谱叠加、滤波和平均处理，从而实现样品平台的物光频谱拓展，然后对拓展后的物光频谱做角谱后向传输算法处理，最后采用解包裹运算重建物光波场的三维相位分布，完成三维成像，同时通过样品平台轴向平移可以获得样品不同深度层次的三维相位信息。因此，本发明通过将离轴数字全息三维成像技术与高速多方向线共焦成像技术相结合，利用多方向线共焦频谱叠加、滤波和平均处理实现横向成像频谱范围拓展，从而提高了横向空间分辨率，通过线共焦成像提高了纵向空间分辨率，从而实现了高分辨率的数字全息三维成像。同时仅通过一维扫描便实现了横向分辨率各向同性，具有高速成像、信噪比高、抗干扰能力强的优点，在高速高分辨高精度三维检测领域具有重要的发展前景。
附图说明
图1是本发明的框架结构示意图；
图2是扫描方向快速转换单元的结构示意图；
图3是中空正三角形内反射镜原理示意图；
图4是高速线阵相机的光敏面、探测扫描光信号和参考返回光信号之间的几何关系图；
图5是主控计算机中的算法流程图。
附图中的标记为：
1、单频单模连续激光器；2、起偏器；3、扩束器；4、半波片；5、矩形光阑；6、柱面透镜；7、偏振分束棱镜；8、第一透镜；9、第一四分之一波片；10、第二透镜；11、全反射镜；12、第三透镜；13、一维快速扫描器；14、第四透镜；15、扫描方向快速转换单元；16、第五透镜；17、显微物镜；18、样品平台；19、第六透镜；20、检偏器；21、第七透镜；22、高速线阵相机；23、高速数字采集卡；24、主机计算机；25、第一二维快速扫描器；26、中空正多边形内反射镜；27、第二二维快速扫描器；28、第二四分之一波片；29、狭缝。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明，但并不作为对本发明限制的依据。
实施例1：单频单模连续激光器产生的光信号首先经过起偏、扩束和柱面透镜聚焦后，由偏振分束棱镜分为参考光信号和探测光信号；将参考光信号通过第一透镜组中继至全反射镜并接收参考返回光信号；将探测光信号通过第二透镜组中继至样品平台并接收探测扫描光信号，其中第二透镜组内设置有一维快速扫描器和线扫描方向快速转换单元，利用一维快速扫描器实现对样品平台的某个方向的线扫描，利用线扫描方向快速转换单元实现多个扫描方向的快速切换；参考返回光信号和探测扫描光信号在偏振分束棱镜合束进入接收光路；在接收光路中，合束光信号先经过狭缝滤除非聚焦平面的杂散光，再通过接收透镜组中继到高速线阵相机的光敏面，中继时经过检偏器的检偏形成离轴干涉场，得到多帧一维线共焦干涉图；多帧一维线共焦干涉图经过高速数字采集，在主控计算机中进行叠加生成多方向重建干涉场，再进行二维傅里叶变换、频谱叠加、滤波和平均处理，实现样品平台的物光频谱拓展，然后对拓展后的物光频谱做角谱后向传输算法处理，最后采用解包裹运算重建物光波场的三维相位分布，完成三维成像，同时通过样品平台轴向平移获得样品不同深度层次的三维相位信息。通过样品平台轴向平移获得样品不同深度层次的三维相位信息，最终实现高速多方向线共焦数字全息高分辨率三维显微成像。
实现上述方法的装置，如图1所示，包括单频单模连续激光器1，单频单模连续激光器1经起偏器2连接有扩束器3，扩束器3的输出端依次连接有半波片4、矩形光阑5、柱面透镜6和偏振分束棱镜7；所述偏振分束棱镜7连接有第一透镜组、第二透镜组和接收透镜组；所述第一透镜组包括与偏振分束棱镜7连接的第一透镜8，第一透镜8经第一四分之一波片9连接有第二透镜10；所述第二透镜10连接有全反射镜11；所述第二透镜组包括与偏振分束棱镜7连接的第二四分之一波片28，第二四分之一波片28连接有第三透镜12，第三透镜12经一维快速扫描器13连接有第四透镜14，第四透镜14经线扫描方向快速转换单元15连接有第五透镜16；如图2所示，所述线扫描方向快速转换单元5由依次连接的第一二维快速扫描器25、中空正多边形内反射镜26和第二二维快速扫描器27组成；所述第五透镜16经显微物镜17连接有样品平台18；所述接收透镜组的前侧设有狭缝29，所述接收透镜组包括与位于狭缝29后侧的第六透镜19，第六透镜19经第七透镜21连接有检偏器20；所述检偏器20连接有高速线阵相机22，高速线阵相机22经高速数字采集卡23连接有主机计算机24；所述主控计算机24分别还与一维快速扫描器13和线扫描防线快速转换单元15连接。
实施例2：在实施例1的基础上，具体的，采用633nm单频单模连续激光器，激光器线宽10kHz，输出功率为20mW，单频单模连续激光器产生的光信号首先经过起偏，确保偏振消光比大于25dB，经过10×扩束器进行扩束，然后再经过半波片和矩形光阑，矩形光阑位于柱面透镜的前焦面，柱面透镜焦距50mm，矩形光阑长度5mm；光信号经过柱面透镜后，在后焦面产生一条聚焦线光束入射至1×2偏振分束棱镜，由1×2偏振分束棱镜分为参考光信号和探测光信号，两路信号分束比由半波片进行控制，大部分能量作为探测光信号，小部分能量作为参考光信号。
将参考光信号通过第一透镜组中继至全反射镜并接收参考返回光信号；所述第一透镜组包括沿参考光信号传输方向依次设置的第一透镜、第一四分之一波片和第二透镜；其中柱面透镜、第一透镜和第二透镜共同构成了4f***，第一四分之一波片用于调整参考返回光信号的光强；第一透镜的焦距为150mm，第二透镜的焦距为20mm；
将探测光信号通过第二透镜组中继至样品平台并接收探测扫描光信号；具体是，所述第二透镜组包括沿探测光信号传输方向依次设置的第三透镜、一维快速扫描器、第四透镜、线扫描方向快速转换单元、第五透镜和显微物镜；其中柱面透镜、第三透镜、第四透镜、第五透镜和显微物镜共同构成了4f***，一维快速扫描器处于第三透镜和第四透镜之间，即第三透镜的后焦面和第四透镜前焦面(两焦平面重合)，一维快速扫描器采用THORLABS的一维振镜***GVS001；线扫描方向快速转换单元处于第四透镜和第五透镜之间，即第四透镜后焦面和第五透镜的前焦面。利用一维快速扫描器实现对样品平台的某个固定角度的线扫描；利用线扫描方向快速转换单元实现多个扫描方向的快速切换，线扫描方向快速转换单元由第一二维快速扫描器、中空正多边形内反射镜和第二二维快速扫描器组成，入射的探测光信号首先经过第一二维快速扫描器偏转，然后发射到中空正多边形内反射镜的某个反射镜面，假设该反射镜面的倾斜角度为θ，再反射到第二二维快速扫描器，由第二二维快速扫描器完成角度补偿后出射，出射探测光信号相对于入射探测光信号实现二维扫描视场2θ角度的旋转；如图3所示，中空正三角形内反射镜可以实现3个角度(0度、120度和240度)的角度旋转，本实施例中的中空正多边形内反射镜为中空正十六边形内反射镜，可以实现4个角度(0度、45度、90度和135度)的旋转，出射的探测光信号通过第五透镜和显微物镜中继到样品平台，然后接受样品平台反射回的探测扫描光信号；
参考返回光信号和探测扫描光信号在偏振分束棱镜合束进入接收光路；在接收光路中，信号先经过狭缝滤除非聚焦平面的杂散光，再通过接收透镜组中继到高速线阵相机的光敏面，中继时经过检偏器的检偏形成离轴干涉场，得到多帧一维线共焦干涉图；所述接收透镜组包括沿参考返回光信号和探测扫描光信号传输方向设置的第六透镜、第七透镜和检偏器；其中柱面透镜、第六透镜和第七透镜共同构成了4f***，狭缝位于第六透镜的前焦面，用于滤除非聚焦平面的杂散光；所述第六透镜焦距50mm，第七透镜焦距300mm；高速线阵相机选用AVIIVA EM2，像素数目512个，像素尺寸14μm×14μm，扫描行频126kHz，扫描线数512条，成像***放大率为45倍，像素分辨率为0.31μm。
多帧一维线共焦干涉图经过高速数字采集(NI的PCIe-1433图像采集卡)，在主控计算机中进行叠加生成某个扫描方向的重建干涉场；
重建干涉场的强度分布表示为：
iHxkΔxH，lΔyH)
＝a(kΔxH，lΔyH)+c(kΔxH，lΔyH)exp[-j2π(ξxΔxH+ξyΔyH)]
+c*(kΔxH，lΔyH)exp[j2π(ξxΔxH+ξyΔyH)]；
式中：
其中O0(kΔxH，lΔyH)为探测扫描光信号的离散振幅分布，
为探测扫描光信号的离散相位分布，r0(kΔxH，lΔyH)为参考返回光信号的离散振幅分布，ξx和ξy分别是x和y方向的空间频率，如图4所示，由参考返回光信号和探测扫描光信号的离轴夹角θ决定；k和l都是整数，-M/2≤k≤M/2，-N/2≤l≤N/2，M和N分别为两个方向的像素数；ΔxH和ΔyH为重建干涉图平面的采样间隔，ΔxH＝ΔyH＝Lx/M＝Ly/N，Lx×Ly是重建干涉场的面积，*表示卷积运算。
在主控计算机中进行叠加生成多个扫描方向的重建干涉场，再进行二维傅里叶变换、频谱叠加、滤波和平均处理，实现样品平台物光频谱的拓展，然后对拓展后的物光频谱做角谱后向传输算法处理，最后采用解包裹运算重建物光波场的三维相位分布，通过样品平台轴向平移获得样品不同深度层次的三维相位信息。如图5所示，具体是：将M个扫描方向的重建干涉场在主控计算机中通过二维傅里叶变换分别求出频谱，将M个频谱叠加、滤波和加权平均，实现样品平台的物光频谱拓展，提高横向空间分辨率，其中二维傅里叶变换计算如下：
IH(pΔξ，qΔη)＝FFT2{iH(kΔxH，lΔyH)}p，q
＝A(pΔξ，qΔη)+C+1(pΔξ-ξx，qΔη-ξy)+C-1(pΔξ+ξx，qΔη+ξy)；
式中：IH表示iH的二维傅里叶变换，p和q都是整数，-M/2≤p≤M/2，-N/2≤q≤N/2；Δξ和Δη分别为二维频谱域内的两个轴的取样间隔，Δξ＝1/(MΔxH)＝1/Lx，Δη＝1/(NΔyH)＝1/Ly；A为干涉频谱的零级峰，C+1为干涉频谱的+1级峰，C-1为干涉频谱-1级峰。通过参考返回光信号和探测扫描光信号的离轴夹角θ调整ξx和ξy，使其A、C+1和C-1在频谱域内足够分开，然后将频谱域内的A和C-1滤除，再将C+1平移至频谱域的中心，再进行角谱后向传输处理，即利用角谱传递函数将C+1向后传播以得到样品平台所在平面的物光频谱；最后采用解包裹运算，即对物光频谱作二维逆傅里叶变换得到重建后的物光波场的复振幅分布：
式中：m和n都是整数，-M/2≤m≤M/2；-N/2≤n≤N/2；λ是波长，H(pΔξ，qΔη)是频谱域内的角谱传递函数；IFFT2是二维逆傅里叶变换；
和
为成像面的采样间隔，大小与线阵相机光敏面像素尺寸相同；
从而得到物光波场的相位分布：
最后进行解包裹运算，得到物光波场的三维相位分布，完成三维成像。通过样品平台轴向平移获得样品不同深度层次的三维相位信息。
本实施例最终的扫描成像速约为67帧/秒。
综上所述，本发明通过将离轴数字全息三维成像技术与高速多方向线共焦成像技术相结合，利用多方向线共焦频谱叠加、滤波和平均处理实现横向成像频谱范围拓展，提高了横向空间分辨率；通过线共焦成像提高了纵向空间分辨率。从而实现了高分辨率的数字全息三维成像。同时还实现了横向分辨率各向同性，具有信噪比高、抗干扰能力强的优点，在高速高分辨高精度三维检测领域具有重要的发展前景。
Claims (10)
1.高速多方向线共焦数字全息三维显微成像方法，其特征在于：激光信号首先经过起偏、扩束和聚焦后，由偏振分束棱镜分为参考光信号和探测光信号；将参考光信号通过第一透镜组中继至全反射镜并接收参考返回光信号；将探测光信号通过第二透镜组中继至样品平台并接收探测扫描光信号，其中探测光信号在中继过程中经过一维快速扫描器和线扫描方向快速转换单元，利用一维快速扫描器实现对样品平台的线扫描，利用线扫描方向快速转换单元实现多个扫描方向的快速切换，完成对样品的多方向共焦线扫描；参考返回光信号和探测扫描光信号在偏振分束棱镜合束进入接收光路；在接收光路中，合束光信号先经过狭缝滤除非聚焦平面的杂散光，再通过检偏形成离轴干涉场，然后通过接收透镜组中继到高速线阵相机，并得到多帧一维线共焦干涉图，经过高速数字采集后，在主控计算机中进行叠加生成多方向重建干涉场，并通过频谱叠加和角谱后向传输处理，实现高速多方向线共焦数字全息显微成像。
2.根据权利要求1所述的高速多方向线共焦数字全息三维显微成像方法，其特征在于：在主控计算机中进行叠加生成多方向重建干涉场后，先进行二维傅里叶变换，再进行频谱叠加、滤波和平均处理，实现样品平台的物光频谱拓展，然后对拓展后的物光频谱做角谱后向传输处理，最后采用解包裹运算重建物光波场的三维相位分布，完成三维成像，同时通过样品平台轴向平移获得样品不同深度层次的三维相位信息。
3.根据权利要求1所述的高速多方向线共焦数字全息三维显微成像方法，其特征在于：所述第一透镜组包括沿参考光信号传输方向依次设置的第一透镜、第一四分之一波片和第二透镜；所述第一四分之一波片用于调整参考返回光信号的光强。
4.根据权利要求1所述的高速多方向线共焦数字全息三维显微成像方法，其特征在于：所述第二透镜组包括沿探测光信号传输方向依次设置的第三透镜、一维快速扫描器、第四透镜、线扫描方向快速转换单元、第五透镜和显微物镜；所述一维快速扫描器处于第三透镜和第四透镜之间；所述线扫描方向快速转换单元处于第四透镜和第五透镜之间。
5.根据权利要求1所述的高速多方向线共焦数字全息三维显微成像方法，其特征在于：所述线扫描方向快速转换单元由第一二维快速扫描器、中空正多边形内反射镜和第二二维快速扫描器组成，入射的探测光信号首先经过第一二维快速扫描器偏转，然后发射到中空正多边形内反射镜的单个反射镜面，再反射到第二二维快速扫描器，由第二二维快速扫描器完成角度补偿后出射，出射的探测光信号相对于入射的探测光信号实现二维扫描视场内一个固定角度的旋转。
6.根据权利要求1所述的高速多方向线共焦数字全息三维显微成像方法，其特征在于：所述接收透镜组包括沿参考返回光信号和探测扫描光信号传输方向设置的第六透镜、第七透镜和检偏器；所述的狭缝位于第六透镜的前焦面，用于滤除非聚焦平面的杂散光。
7.根据权利要求1所述的高速多方向线共焦数字全息三维显微成像方法，其特征在于：所述的重建干涉场的强度分布表示为：
iH(kΔxH，lΔyH)＝a(kΔxH，lΔyH)+c(kΔxH，lΔyH)exp[-j2π(ξxΔxH+ξyΔyH)]+c*(kΔxH，lΔyH)exp[j2π(ξxΔxH+ξyΔyH)]；
式中：
其中O0(kΔxH，lΔyH)为探测扫描光信号的离散振幅分布，
为探测扫描光信号的离散相位分布，r0(kΔxH，lΔyH)为参考返回光信号的离散振幅分布，ξx和ξy分别是x和y方向的空间频率，由参考返回光信号和探测扫描光信号的离轴夹角θ决定；k和l都是整数，-M/2≤k≤M/2，-N/2≤l≤N/2，M和N分别为两个方向的像素数；ΔxH和ΔyH为重建干涉图平面的采样间隔，ΔxH＝ΔyH＝Lx/M＝Ly/N，Lx×Ly是重建干涉场的面积，*表示卷积运算。
8.根据权利要求2所述的高速多方向线共焦数字全息三维显微成像方法，其特征在于：将M个扫描方向的重建干涉场在主控计算机中通过二维傅里叶变换分别求出频谱，将M个频谱叠加、滤波和加权平均处理，实现样品平台的物光频谱拓展，其中二维傅里叶变换计算如下：
IH(pΔξ，qΔη)＝FFT2{iH(kΔxH，lΔyH)}p，q＝A(pΔξ，qΔη)+C+1(pΔξ-ξx，qΔη-ξy)+C-1(pΔξ+ξx，qΔη+ξy)；
式中：IH表示iH的二维傅里叶变换，p和q都是整数，-M/2≤p≤M/2，-N/2≤q≤N/2；Δξ和Δη分别为频谱域内两个轴的取样间隔，Δξ＝1/(MΔxH)＝1/Lx，Δη＝1/(NΔyH)＝1/Ly；A为干涉频谱的零级峰，C+1为干涉频谱的+1级峰，C-1为干涉频谱-1级峰；
通过参考返回光信号和探测扫描光信号的离轴夹角θ调整ξx和ξy，使其A、C+1和C-1在频谱域内足够分开，然后将频谱域内的A和C-1滤除，再将C+1平移至频谱域的中心，再进行角谱后向传输处理，即利用角谱传递函数将C+1向后传播以得到样品平台所在平面的物光频谱；最后采用解包裹运算，即对物光频谱作二维逆傅里叶变换得到重建后的物光波场的复振幅分布：
式中：m和n都是整数，-M/2≤m≤M/2；-N/2≤n≤N/2；λ是波长，H(pΔξ，qΔη)是频谱域内的角谱传递函数；IFFT2是二维逆傅里叶变换；
和
为成像面的采样间隔，大小与高速线阵相机光敏面像素尺寸相同；
从而得到物光波场的相位分布：
9.根据权利要求1-8任一项所述的高速多方向线共焦数字全息三维显微成像方法的装置，其特征在于：包括单频单模连续激光器(1)，单频单模连续激光器(1)经起偏器(2)连接有扩束器(3)，扩束器(3)的输出端依次连接有半波片(4)、矩形光阑(5)、柱面透镜(6)和偏振分束棱镜(7)；所述偏振分束棱镜(7)连接有第一透镜组、第二透镜组和接收透镜组；所述第一透镜组包括与偏振分束棱镜(7)连接的第一透镜(8)，第一透镜(8)经第一四分之一波片(9)连接有第二透镜(10)；所述第二透镜(10)连接有全反射镜(11)；所述第二透镜组包括与偏振分束棱镜(7)连接的第二四分之一波片(28)，第二四分之一波片(28)连接有第三透镜(12)，第三透镜(12)经一维快速扫描器(13)连接有第四透镜(14)，第四透镜(14)经线扫描方向快速转换单元(15)连接有第五透镜(16)；所述第五透镜(16)经显微物镜(17)连接有样品平台(18)；所述接收透镜组的前侧设有狭缝(29)，所述接收透镜组包括与位于狭缝(29)后侧的第六透镜(19)，第六透镜(19)经第七透镜(21)连接有检偏器(20)；所述检偏器(20)连接有高速线阵相机(22)，高速线阵相机(22)经高速数字采集卡(23)连接有主机计算机(24)；所述主控计算机(24)分别还与一维快速扫描器(13)和线扫描防线快速转换单元(15)连接。
10.根据权利要求9所述的高速多方向线共焦数字全息三维显微成像方法的装置，其特征在于：所述线扫描方向快速转换单元(5)由依次连接的第一二维快速扫描器(25)、中空正多边形内反射镜(26)和第二二维快速扫描器(27)组成。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110218940.8A CN113009801B (zh) | 2021-02-26 | 2021-02-26 | 高速多方向线共焦数字全息三维显微成像方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110218940.8A CN113009801B (zh) | 2021-02-26 | 2021-02-26 | 高速多方向线共焦数字全息三维显微成像方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113009801A CN113009801A (zh) | 2021-06-22 |
| CN113009801B true CN113009801B (zh) | 2022-03-11 |
Family
ID=76386436
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110218940.8A Active CN113009801B (zh) | 2021-02-26 | 2021-02-26 | 高速多方向线共焦数字全息三维显微成像方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113009801B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115481361A (zh) * | 2021-05-31 | 2022-12-16 | 华为技术有限公司 | 一种光计算装置和方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1475878A (zh) * | 2003-06-27 | 2004-02-18 | 中国科学院上海光学精密机械研究所 | 重构x射线全息图的共焦扫描装置 |
| US8896840B2 (en) * | 2012-04-25 | 2014-11-25 | Canon Kabushiki Kaisha | Interferometric method and digital holographic microscope |
| CN102944989B (zh) * | 2012-10-23 | 2015-07-01 | 深圳大学 | 一种相移数字全息高速成像方法及*** |
| CN107024848B (zh) * | 2017-04-26 | 2020-08-18 | 中国工程物理研究院激光聚变研究中心 | 基于逐点扫描的太赫兹透射式全息成像***及成像方法 |
-
2021
- 2021-02-26 CN CN202110218940.8A patent/CN113009801B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN113009801A (zh) | 2021-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11644791B2 (en) | Holographic imaging device and data processing method therefor | |
| CN109709786B (zh) | 一种超分辨率数字全息成像***和成像方法 | |
| US8405059B2 (en) | Method and apparatus for improving the resolution and/or sectioning ability of an imaging system | |
| CN101983313A (zh) | 测量物体表面形貌的设备和方法 | |
| CN111257227B (zh) | 基于偏振自相关的暗场共焦显微测量装置和方法 | |
| CN111121675B (zh) | 一种用于微球表面显微干涉测量的视场扩展方法 | |
| JP7352292B2 (ja) | ホログラフィック撮像装置およびホログラフィック撮像方法 | |
| CN107885070B (zh) | 一种基于slm的非相干数字全息单次曝光成像方法与*** | |
| CN111650203A (zh) | 一种微球内表面缺陷测量方法 | |
| CN116183568B (zh) | 一种三维结构光照明超分辨显微成像的高保真重构的方法和装置 | |
| CN113009801B (zh) | 高速多方向线共焦数字全息三维显微成像方法及装置 | |
| CN116448728A (zh) | 一种基于振镜扫描的双光束场干涉光片显微成像方法和装置 | |
| CN109283821B (zh) | 基于涡旋透镜的相移数字全息单次曝光成像装置及方法 | |
| CN101122774B (zh) | 高分辨率数字全息像的获取装置 | |
| CN113031422B (zh) | 一种全息成像装置 | |
| CN108957999B (zh) | 基于相位型涡旋透镜的相移全息装置及成像方法 | |
| CN114624981B (zh) | 一种基于时间展宽的超快全息显微成像方法及其*** | |
| Phan et al. | Super-resolution digital holographic microscopy for three dimensional sample using multipoint light source illumination | |
| CN113075871B (zh) | 高速双光束正交线共焦数字全息三维显微成像方法及装置 | |
| CN113218517B (zh) | 一种计算成像同轴全息光纤模式分解方法及装置 | |
| CN111122510B (zh) | 基于f-p干涉仪透射式正交偏振相位显微成像装置 | |
| CN111459004B (zh) | 双光路合成孔径全息图拼接方法 | |
| TW201638880A (zh) | 影像處理方法 | |
| CN111025876A (zh) | 基于压电陶瓷的透射式相位显微成像测量*** | |
| CN116223374B (zh) | 一种光学相干断层多点分布式成像方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |