CN102402172A - 一种三维实时超分辨数字全息记录系统 - Google Patents
一种三维实时超分辨数字全息记录系统 Download PDFInfo
- Publication number
- CN102402172A CN102402172A CN201110389124XA CN201110389124A CN102402172A CN 102402172 A CN102402172 A CN 102402172A CN 201110389124X A CN201110389124X A CN 201110389124XA CN 201110389124 A CN201110389124 A CN 201110389124A CN 102402172 A CN102402172 A CN 102402172A
- Authority
- CN
- China
- Prior art keywords
- light
- wavelength
- beam splitter
- pulse laser
- recording system
- 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.)
- Granted
Links
- 238000001093 holography Methods 0.000 title abstract description 8
- 230000010287 polarization Effects 0.000 claims abstract description 65
- 230000003287 optical effect Effects 0.000 claims abstract description 33
- 239000013078 crystal Substances 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 7
- 230000001678 irradiating effect Effects 0.000 claims description 13
- 238000001228 spectrum Methods 0.000 claims description 10
- 238000002310 reflectometry Methods 0.000 claims description 8
- 230000003595 spectral effect Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 abstract description 10
- 238000003384 imaging method Methods 0.000 abstract description 10
- 238000000034 method Methods 0.000 description 20
- 230000008569 process Effects 0.000 description 11
- 238000009647 digital holographic microscopy Methods 0.000 description 7
- 238000005286 illumination Methods 0.000 description 7
- 230000001052 transient effect Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 3
- 230000002035 prolonged effect Effects 0.000 description 3
- 230000001681 protective effect Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 229910052709 silver Inorganic materials 0.000 description 3
- 239000004332 silver Substances 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000009022 nonlinear effect Effects 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241001123862 Mico Species 0.000 description 1
- 241000656145 Thyrsites atun Species 0.000 description 1
- QBLDFAIABQKINO-UHFFFAOYSA-N barium borate Chemical compound [Ba+2].[O-]B=O.[O-]B=O QBLDFAIABQKINO-UHFFFAOYSA-N 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- VCZFPTGOQQOZGI-UHFFFAOYSA-N lithium bis(oxoboranyloxy)borinate Chemical compound [Li+].[O-]B(OB=O)OB=O VCZFPTGOQQOZGI-UHFFFAOYSA-N 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000000399 optical microscopy Methods 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- WYOHGPUPVHHUGO-UHFFFAOYSA-K potassium;oxygen(2-);titanium(4+);phosphate Chemical compound [O-2].[K+].[Ti+4].[O-]P([O-])([O-])=O WYOHGPUPVHHUGO-UHFFFAOYSA-K 0.000 description 1
Images
Landscapes
- Holo Graphy (AREA)
Abstract
本发明涉及一种三维实时超分辨数字全息记录系统，属三维显微成像技术领域。由入射光产生系统和记录系统组成，入射光产生系统包括飞秒脉冲激光器、倍频晶体、二向色镜、光程延迟系统和非偏振分束镜，记录系统包括显微物镜和系列偏振分束镜、非偏振分束镜、反射镜、数字相机和计算机全息图像处理再现系统。利用波分、角分、时分和偏振复用技术，采用光学倍频晶体，将激光器发出的激光分为两个不同波长的激光先后进入记录系统，经系列分束、高分辨率、大视场和景深，在CCD的一帧画面中实现三维实时超分辨成像的数字全息记录。具有系统简洁、记录信息多、可实现空间和时间分辨率共同提高等优点，可广泛用于显微物体的高分辨率实时三维成像。
Description
技术领域
本发明涉及一种三维实时超分辨数字全息记录系统，具体地说，是利用数字全息显微技术实现三维实时超高空间分辨率的一种成像系统，属于三维显微成像技术领域。
背景技术
数字全息显微术是近些年发展起来的一种新的显微成像技术，现已在生物细胞成像、MEMS器件和微光学器件的检测、微小物体的变形和振动测量中得到应用。它利用CCD或者CMOS等数字相机记录物体与参考光干涉形成的全息图，并利用计算机模拟衍射过程再现被记录显微物体的波前信息（相位和振幅），利用二者的数据信息能够直接获得物体的三维分布。数字全息显微术具有如下优点：（1）不需要对物体进行扫描，利用复振幅信息直接可以得到三维数据信息，因而三维分布信息的获取速度较快；（2）通过数字调焦技术，可获得待测样品各个层面的分布，不需要对物体进行切片；（3）对测量环境没有特殊要求，而且可以直接对活体生物进行成像，不需要在物体表面镀膜或者对物体染色；（4）测量范围较宽可以从毫米到亚微米。
通过十几年的发展，数字全息显微术已走出实验室逐步向产品化发展。瑞士Lyncee Tec公司开发了世界上第一部数字全息显微镜，借助它可以完成微小物体的实时记录，每秒可连续记录变形过程中15幅全息图；2009年美国专利局授权了一种可以实时记录显微物体全息图的系统，该系统采用三部激光器发出的三个不同波长的光作为光源，实时记录MEMS器件的变形过程；2008年日本研究人员在论文《Real-time digital holographic microscopy using the graphic processing unit》（Opt. Express (2008)16, 11776-11781）中，提出采用GPU实时再现动态全息图的方法，利用GPU系统显著提高了再现速度，每秒可以再现24幅全息图，但全息图的采集速度等同于CCD的采集速度；2010年印度研究人员在论文《Real-Time digital holographic microscopy for phase contrast 3D imaging of dynamic phenomena》（Journal of display technology (2010)6,500-505）中，提出的三维实时记录系统，它的记录速度也受限于CCD的采集速度。上述几种典型的数字全息实时全息记录系统，虽然都可以对显微物体进行动态记录，但它们所能达到的时间分辨率也仅为CCD的帧频，空间分辨率决定于显微物镜的数值孔径。2006年Xiaolei Wang等在论文《Pulsed digital holography system recording ultrafast process of the femtosecond order》（Opt. Lett. (2006)31, 1636-1638）以及2010年Linwei Zhu等在论文《Femtosecond off-axis digital holography for monitoring dynamic surface deformation》（Appl. Opt. (2010) 49, 2510-2518）提出的飞秒数字全息记录系统虽然可以突破CCD的帧频限制，达到飞秒量级的时间分辨率，但是这些系统的空间分辨率没有突破显微镜的分辨极限，而且不是实时采集，只能获取几个时间点的信息。
数字全息显微术的分辨率和传统光学显微系统一样，取决于显微物镜的数值孔径（NA）以及光源的波长。当光源的波长固定，它的分辨率只与数值孔径有关，使用大数值孔径的显微物镜可以获得高分辨率的再现像，但大数值孔径的物镜会使物方视场缩小、工作距离和焦深变短，这就限制了该技术在生物、医学和材料科学方面的应用。因为在这些测量领域中，尤其是记录一些微小物体的变形或者是微生物体的运动轨迹，不仅要求实时和高分辨率，而且要求大视场和大景深。
为了克服空间分辨率与物方视场、工作距离和焦深之间的矛盾，在数字全息显微系统中引入合成孔径技术，它虽然可以克服这个矛盾，获得超过衍射极限的分辨率，但是系统不具有实时记录的能力。在文献《Single-exposure two-dimensional superresolution in digital holography using a vertical cavity surface-emitting laser source array》（Opt. Lett. (2011) 7, 1149-1151）中，V. Mico等提出了利用VCSEL光源提供的多角度照射光可以替代移动CCD，从而达到改进传统数字全息显微合成孔径系统的目的，但当被测物体频谱分布比较复杂时，很容易造成频谱混叠，因而缺乏普适性。另外，在数字全息显微记录光路中插入光栅的方法，也被文献《Super-resolution in digital holography by a two-dimensional dynamic phase grating》（Opt. Express (2008) 16, 17107-17118）报道；作为一种可以提高分辨率的方法，它不仅结构简单，而且具有较高普适性，但是为了避免频谱混叠，需要对物体的尺寸以及频谱范围进行限制，进而物方视场也受限。
总之，尽管现有的数字全息显微系统可以独立拥有超高的空间分辨率，或者较高的时间分辨率，但都没能在空间分辨率和时间分辨率的提高方面取得统一。
发明内容
本发明的目的是针对已有数字全息显微系统不能同时具有高时间、高空间分辨率、大视场和大景深的缺点，提供一种三维实时超分辨数字全息记录系统，采用波分、角分、时分和偏振复用混合技术，在数字相机的同一帧画面中同时记录包含两个瞬态时刻、包含每个瞬态时刻高、低频信息的四幅全息图，使数字全息显微系统同时具有高时间、空间分辨率、大视场、长工作距离和大景深的特点。
本发明描述的三维实时超分辨数字全息记录系统，包括入射照射光源和记录系统两个部分。入射光产生部分包括飞秒级脉冲激光器1，可按一定比例产生另一波长激光脉冲的倍频晶体2，对两波长脉冲激光中某一波长脉冲激光进行反射、同时允许另一波长脉冲激光直接透过的二向色镜3，用于延长两波长脉冲激光中某一波长脉冲激光光程的光程延迟系统（反射镜组合4和5），将长、短两波长脉冲激光汇合成一个方向传播的非偏振分束镜6；记录系统包括可将进入该系统的毎束入射脉冲激光分成两束偏正方向沿X和Y方向、相互垂直的线性偏振光的偏振分束镜7，位于偏振分束镜7与反射镜9和12之间、将来自偏振分束镜7的沿X方向水平偏正的光束分成参考光A和照射光B的非偏振分束镜8，位于非偏振分束镜8与被测物体13之间、将照射光B反射后沿一定方向照射被测物体13的反射镜12，位于偏振分束镜7与反射镜16和17之间、将来自偏振分束镜7的沿Y方向垂直偏正的光分成照射光C和参考光D的非偏振分束镜15，位于被测物体13与非偏振分束镜15之间、将照射光C反射后沿另一方向照射被测物体13的反射镜17，将参考光A和D反射后射向非偏振分束镜10的反射镜9和16；位于被测物体13和非偏振分束镜10之间、将物光放大的显微物镜14,用于调整来自显微物镜14且分别携带物体高、低频信息的物光和来自反射镜9、16的参考光A和D之间夹角的非偏振分束镜10，用于记录来自非偏振分束镜10的物光和参考光所形成全息图的数字相机CCD11，以及计算机三维数字全息图像处理再现系统。
所述的飞秒级脉冲激光器1是可输出中心波长为775nm～1030nm、单脉冲宽度在飞秒量级的脉冲激光器。可采用符合上述条件的普通市售激光器。
所述的倍频晶体2为可在透过原长波波长脉冲激光的同时产生另一波长为1/2原波长的短波脉冲激光、倍频效率为30%～70%的普通非线性光学晶体2，如：BBO，KTP等。通过倍频晶体2，可将激光器发出的光分割成长波长和短波长两束不同波长的脉冲激光，长波长脉冲的中心波长为短波长脉冲中心波长的2倍。
所述光程延迟系统包括可使两波长脉冲激光中某一波长脉冲激光到达记录系统的时间滞后或超前于另一波长脉冲激光的反射镜4和5，反射镜4和5分别与二向色镜3和非偏振分束镜6相对应。光程延迟系统使两波长脉冲激光中某一波长脉冲激光到达记录系统滞后或超前于另一波长脉冲激光的时间差为50fs～33ps。具体通过调整反射镜4、5之间的间距，以及反射镜4、5与二向色镜3和非偏振分束镜6之间的距离，从而延长或缩短光程，使两束脉冲激光产生时间差。具体时间差科根据实际需要，作适应性调整。
所述两块反射镜4和5可用一块屋脊棱镜代替，根据短波长脉冲滞后或超前长波长脉冲的实际需要，反射镜4、5或屋脊棱镜可选择性地置于入射光产生部分的长波长脉冲激光光程或短波长脉冲激光光程中。两波长脉冲激光中某一波长脉冲激光在二向色镜3反射镜4、5以及非偏振分束镜6之间的反射过程中经历较长的光程，而另一波长脉冲激光直接透过二向色镜3和非偏振分束镜6，经历的光程较短。两不同波长脉冲经历的光程差决定了输出时间差，其中某一波长脉冲激光滞后于另一波长激光脉冲到达记录系统，它们的时间间隔在飞秒到皮秒量级之间，可使其后所形成的四幅全息图记录在数字相机CCD的一帧画面中，并且不发生混叠。
所述分别从不同方向照射被测物体的两束照射光B和C中，照射光B沿物体平面的法线方向照射被测物体13、照射光C与物平面法线呈sin-1（NA）角度照射被测物体13，其中，NA是显微物镜（14）的数值孔径。照射光C的具体照射角度根据实际需要，由选定的显微物镜14的数值孔径决定。
所述非偏振分束镜6、7、8、10、15均为直透光与反射光分光比为1:1的普通宽光谱分光棱镜，反射镜4、5、9、12、16和17均为普通宽光谱、高反射率反射镜，显微物镜14为可在记录前对被记录物体进行预先放大的普通小数值孔径、大视场、长距离和焦深的物镜，数字相机CCD11为普通黑白面阵型数字相机。非偏振分束镜6、7、8、10、15和反射镜9、12、16和17的光谱范围与激光光源的光谱范围相适应可在不改变光原有偏振状态的情况下，对通过它们的中心波长为387nm～1030nm激光脉冲进行分束或高效率反射。
本系统使用时，飞秒激光器1输出波长较长的激光脉冲，该光透过倍频晶体2（如：BBO，KTP等）后，由于倍频晶体的非线性效应，在直接透过波长较长的激光脉冲同时，还将产生中心波长为原波长的一半的另一激光脉冲；不同波长的两束激光脉冲同时射向二向色镜3，其中，波长较长的脉冲将直接透过二向色镜3和非偏振分束镜6，进入记录系统部分；而短波长脉冲将被二向色镜3反射，依次被第一块反射镜4、第二块反射镜5和非偏振分束镜6反射后，形成飞秒到皮秒量级的输出时间差，两束脉冲将先后进入记录系统部分。
因波长较长的激光脉冲经历的光程短，它将先于波长较短的激光脉冲进入记录系统，它将被偏振分束镜7分解成两束线偏振光，分别沿水平X和垂直Y方向偏振。沿水平X方向偏振的激光脉冲被非偏振分束镜8分束后形成A和B两束光，A光束被反射镜9和非偏振分束镜10反射后射向数字相机CCD11，B光束被反射镜12反射后依次水平透过被测物体13、显微物镜14和非偏振反射镜10，射向CCD11，与A光束发生干涉形成第一幅全息图。沿垂直Y方向偏振的激光脉冲被非偏振分束镜15分束后形成C和D两束光，D光束被反射镜16和非偏振分束镜10反射后射向CCD11，C光束被反射镜17反射后倾斜照射被测物体13、再透过显微物镜14和非偏振反射镜10，射向CCD11，与D光束发生干涉形成第二幅全息图。波长较短的激光脉冲经飞秒到皮秒的时间延迟后进入记录系统部分，并沿着与波长较长脉冲光相同的路径传播，最终形成第三幅和第四幅全息图。两个瞬态时刻同一动态过程的四幅全息图，均由数字相机CCD11复合记录在一帧画面中，并传输到计算机中进行数字再现，通过对复合全息图频域滤波处理，将高频和低频信息相互分离，合成复振幅从而获得三维超高时间和空间分辨率再现图像。
本发明通过控制不同波长两束激光脉冲的延迟时间，保证了四幅全息图记录在数字相机CCD的一帧画面中。同时，由于采用偏振复用技术，使得在一个瞬态时刻下记录到的包含高频和低频信息的全息图相互之间不发生干扰；采用时分复用技术，保证两个瞬态过程的两组全息图(每组包括一幅高频信息全息图和一幅低频信息全息图)相互之间不发生干扰；采用角分和波分复用技术，保证了记录在一帧画面中的四幅全息图的再现像不发生混叠。
本发明实现了超高空间分辨率的实时三维成像，与现有技术相比具有如下优点：
（1）由于本发明只使用一部激光器，较使用多部激光器的方法和系统来说，更加简洁；
（2）用同轴光和离轴光同时照射物体，可以保证CCD能同时记录到物体高频和低频信息，保证了在显微物镜的数值孔径不大的情况下，也能获得在系统截止频率以外的物体信息，保证了系统具有很高的分辨率而且具有大视场和景深；较只能记录到物体的低频或者高频信息的传统方法相比，可以记录到更多的物体信息，并且空间分辨率能突破成像系统的极限分辨率，可以实现空间和时间分辨率的共同提高，完成三维超高空间分辨率的再现；
（3）利用偏振方向不同而且到达CCD时间不同的四组参物光组，可以保证四幅全息图以非相干的形式记录在CCD上而不相互干扰；
（4）通过控制波长较长和较短两束激光脉冲的光程（时间）差，可以实时记录微小物体的变化过程，实时记录速度高于CCD的采集速度，而且可以捕捉在飞秒到皮秒量级瞬态变化的信息。
附图说明
图1为本发明系统示意图；
图2为本发明入射光产生系统示意图。
图中：1-飞秒激光器，2-倍频晶体，3-二向色镜，4-第一反射镜，5-第二反射镜，6-非偏振分束镜，7-偏振分束镜，8-非偏振分束镜，9-反射镜，10-非偏振分束镜，11-数字相机CCD，12-反射镜，13-被测物体，14-显微物镜，15-非偏振分束镜，16-反射镜，17-反射镜。
具体实施方式
下面结合附图和实施例对本发明作进一步阐述，但本发明的保护范围不限于所述内容。
实施例1：如图1、2所示，该三维实时超分辨数字全息记录系统，包括入射光产生系统和记录系统两个部分。入射光产生部分包括飞秒脉冲激光器1，改变入射激光脉冲波长的倍频晶体2，能将来自倍频晶体2的两不同波长激光脉冲分开并具有输出时间差的二向色镜3，反射镜4、5以及非偏振分束镜6；记录系统包括将进入该系统（来自于非偏振分束镜6）的毎束入射的激光脉冲分解成水平X方向和垂直Y方向线偏振光束的偏振分束镜7，位于偏振分束镜7与反射镜 9之间、将来自偏振分束镜7沿X方向偏振光分解成参考光A和照射光B的非偏振分束镜8，位于非偏振分束镜8与被测物体13之间、将照射光B反射后垂直于物体平面（沿物体平面的法线方向）照射被测物体13的反射镜12，位于偏振分束镜7与反射镜 16和17之间、将来自偏振分束镜7的沿Y方向偏振的光分解成照射光C和参考光D的非偏振分束镜15，位于被测物体13与非偏振分束镜15之间、将照射光C反射后倾斜sin-1（NA）=7o照射被测物体13的反射镜17，将参考光A和D反射后射向非偏振分束镜10的反射镜9和16，位于被测物体13和非偏振分束镜10之间、将物光放大后射向非偏振分束镜10的显微物镜14,将来自显微物镜14的低频信息（垂直物光产生）和高频信息（倾斜物光产生）与来自反射镜9、16的参考光A和D汇合的非偏振分束镜10，用于记录来自非偏振分束镜10的参考光（A和D）和物光（携带物体低频和高频信息）干涉图像的数字相机CCD11，以及普通计算机三维数字全息图像处理再现系统。
该系统中，激光器1为输出中心波长为800nm、脉冲间隔为35fs的Tsunami 3941-35型飞秒激光振荡器，倍频晶体2为可使原输入脉冲激光频率提高两倍，也就是使通过它的激光波长变为原有波长的1/2倍的BBO偏硼酸钡晶体（尺寸为：0.01 mm×7mm×4mm，晶体响应波长范围为189-3500nm、倍频效率约为40%）。所产生的激光脉冲光学延迟由一块二向色镜3、反射面对称放置的两块反射镜4、5以及一块非偏振分束镜实现，二向色镜是具有两向色性的分光板，其对波长为400nm的光具有高反射率、对800nm的光具有较高的透过率，通过反射过程延长中心波长400nm的激光脉冲经历的光程，使其滞后于波长为800nm的激光脉冲50fs输出，保证其后所形成的四幅全息图记录在数字相机CCD的一帧画面中，时间间隔可以通过改变光程差长短调整。非偏振分束镜6、7、8、10、15均为分光比为1:1的普通直透光与反射光宽宽带偏振分光棱镜，它能够将中心波长为800nm和400nm的入射光分别分成两束光，但是不改变其偏振方向，保持其原有偏振方向；反射镜9、12、16和17为能对近红外800nm的长波长激光和400nm的较短波长的光具有大于90%反射率的保护银反射镜；显微物镜14为可在记录前对被记录物体进行预先放大的普通小数值孔径、大视场、长工作距离和深焦物镜（NA=0.12，3×），因此，倾斜照射物体的光束与物平面法线的夹角sin-1（NA）=7o；数字相机CCD11为用于记录全息图的普通黑白面阵型CCD，可以选用PIKE F-505B型CCD，其感光面尺寸为8.5 mm×7.1mm，像素为2452×2054，像元尺寸为3.45μm×3.45μm，帧频为每秒15幅，电子快门70μs-67s，按照每帧画面可以记录两个时间点信息的特点，每秒可以记录30幅全息图，可实现高速实时采集。
在该系统中，由于飞秒激光器1输出的中心波长为800nm的激光，直接射向BBO晶体2，由于晶体的非线性效应，飞秒激光透过BBO晶体后，将产生400nm的倍频激光脉冲，但由于BBO的倍频效率（40%）不能达到100%，因此，800nm的激光脉冲也以60%的比例透过BBO晶体，两波长的脉冲激光将同时射向二向色镜3。400nm的激光脉冲将被二向色镜3反射，经过一定的光路延迟后，依次被第一块反射镜4、第二块反射镜5以及被非偏振分束镜6反射后，将进入记录系统部分，通过调整两块反射镜的位置，控制波长为800nm和400nm激光脉冲的时间间隔为50fs；另外一束800nm的脉冲光将直接透过二向色镜3和透过非偏振分束镜6，将直接进入记录系统部分。
由于中心波长为800nm的脉冲经历的光程短，因此该束光将在t时刻进入记录系统部分，该束脉冲将被偏振分束镜7分成两束偏振方向相互垂直的线性偏振光（即，一束沿x方向偏振，另一束沿y方向偏振）。
沿x方向偏振的脉冲光，被非偏振分束镜8分成两束光（A和B），A光束被反射镜9和非偏振分束镜10反射后射向CCD11。在CCD记录平面（x-y面）上，A光束的分布可以表示成：
其中：
，
是参考光束A与x轴的夹角，l1=800nm。B光束被反射镜12反射后，透过物体13和显微物镜14，在显微物镜后得到一个物体放大的实像，这个光场经反射镜10反射后，也将射向CCD11，当放大像和CCD之间的距离满足菲涅耳近似条件，B光束在CCD记录平面上的分布可以表示成：
其中，t(x, y)是放大后的物体所成实像的复振幅分布，h(x, y)是系统点扩散函数，
是二维卷积。
沿y方向偏振、波长为800nm的脉冲光，被非偏振分束镜15分束后形成两束光（C和D），其中C光束被反射镜16和非偏振分束镜10反射后，射向记录介质CCD11，C光束在CCD记录面上的分布可以表示为：
其中，
，
是参考光束C与y轴的夹角。
D光束被反射镜17反射后，将倾斜照射被测物体13、而后透过显微物镜14和非偏振反射镜10，射向CCD11，D光束在CCD记录面上的分布可以表示为：
其中，为倾斜照射物体的光波的复振幅分布。从式可以看出，倾斜照射物体会使物光场发生平移，依据角谱理论可知，CCD所能记录到的物体频谱范围也将发生变化，原来记录不到的频谱信息，会因为使用了倾斜照射，进入到CCD记录范围之内，通过控制倾斜照射的角度，传统照射方式记录不到的高频信息可以被记录。通过合成高频和低频信息，相当于系统的数值孔径由NA扩大到2NA，因此，系统的空间分辨率会提高一倍，甚至超出了系统的极限分辨率。
中心波长400nm的高频光经历一定时间延迟Dt=50fs后（Dt=DL/c, DL 是光程差，c是光速，DL=30μm），即t+Dt时刻，进入记录系统部分，这束脉冲将被偏振分束镜7分成沿x方向偏振和沿y方向偏振的两束线性偏振光。
沿x方向偏振的脉冲光，被非偏振分束镜8分成两束光（A和B），A光束被反射镜9和非偏振分束镜10反射后射向CCD11。在CCD记录平面（x-y面）上，A光束的分布可以表示成：
其中 l2=400nm。B光束被反射镜12反射后，透过物体13和显微物镜14，在显微物镜后得到一个物体放大的实像，这个光场经反射镜10反射后，也将射向CCD11，当放大像和CCD之间的距离满足菲涅耳近似条件，B光束在CCD记录平面上的分布可以表示成：
沿y方向偏振，波长为400 nm的脉冲光，被非偏振分束镜15分束后形成两束光（C和D），其中C光束被反射镜16和非偏振分束镜10反射后，射向数字相机CCD11，C光束在CCD记录面上的分布可以表示为：
（7）
D光束被反射镜17反射后，将倾斜照射被测物体13、而后透过显微物镜14和非偏振反射镜10，射向CCD11，D光束在CCD记录面上的分布可以表示为：
其中，
为倾斜照射物体的光波的复振幅分布。
若400nm和800nm之间的时间间隔在CCD一帧画面的积分时间内，八束光将会将会先后到达CCD，但由于偏振方向相互垂直的两束光，以及达到时间不同的两束光之间也不能发生干涉，因此，八束光之间，只能形成四幅全息图。
其中，I1、I2、I3和I4代表四幅全息图的强度分布。
再现中，对复合全息图进行傅里叶变换，I1和I2参考光的载波频率分别在x和y方向， I1和I2的正负一级再现像的中心分别位于频域的横坐标和纵坐标轴上，因此，它们的再现像在频域不会相互干扰。I3和I4的正负一级像也是同样的情况，它们在频域里也不会相互干扰。
由于中心波长为800nm和400nm激光脉冲的波长相差一倍，则两个波长形成的干涉条纹频率也会相差一倍，因此，I1和I2以及I3和I4的频谱分布在频域中不会发生频谱重叠。
通过频域滤波，可以将t时刻记录在I1和I2全息图中的正一级频谱信息滤出来，再通过傅里叶逆变换和相位校正，将两复振幅相加，获得了该时刻超高分辨率的再现像的复振幅分布。同样，通过频域滤波，可以将t+Dt时刻记录在I3和I4全息图中的正一级频谱信息滤出来，再通过傅里叶逆变换和相位校正，将两复振幅相加，获得第二个时刻的超高分辨率再现像的复振幅分布。由于Dt=100fs在飞秒到皮秒范围内，所以该方法与系统可以捕捉瞬态变化信息。通过对复振幅取模平方可以获得再现像的强度分布，通过取复振幅的幅角可以获得其相位信息，综合相位和振幅信息可以得到不同时刻的被记录物体的三维信息。这样，CCD一帧画面上可以记录四幅全息图，也就是两个瞬态时刻的超高分辨率的三维信息，通过CCD连续采集，将会得到微小物体变化过程高分辨率、大视场，实时三维分布图像。
实施例2：如图1、2所示，该三维实时超分辨数字全息记录系统与实施例1相同。所用激光器1为输出激光脉冲中心波长为1030nm、脉冲时间间隔小于100fs的ORANGE型掺镱光纤振荡器，倍频晶体2为可使原输入脉冲激光频率提高两倍，也就是使通过它的激光波长变为原有波长的1/2倍的KTP磷酸钛氧钾（尺寸为：2 mm×5 mm×5 mm，晶体响应波长范围为350-4500nm、倍频效率约为70%）。所产生的激光脉冲光学延迟系统由一块二向色镜3、反射面对称放置的两块反射镜4、5以及一块非偏振分束镜6实现，二向色镜是具有两向色性的分光板，其对波长为515nm的光具有高反射率、对1030nm的光具有较高的透过率，通过反射过程延长中心波长515nm的激光脉冲经历的光程，将滞后于波长为1030nm的激光脉冲800fs输出。非偏振分束镜和反射镜的布置使每两束照射光中的一束沿垂直于物体平面的方向照射被测物体，另一束倾斜12o照射被测物体。非偏振分束镜6、7、8、10和15均为分光比为1:1的普通直透光与反射光宽宽带偏振分光棱镜，它能够将波长较长的脉冲激光1030nm和波长较短的脉冲激光515nm的入射光分成两束光，但是不改变其偏振方向，保持其原有偏振方向；反射镜4、5、9、12、16和17为能对近红外1030nm的光和515nm的光具有大于95%反射率的保护银反射镜；显微物镜14为可在记录前对被记录物体进行预先放大的普通小数值孔径、大视场、长工作距离和深焦物镜（5×NA：0.2），倾斜照射物体的光束与物平面法线的夹角sin-1（NA）=12o；数字相机CCD11为用于记录全息图的普通黑白面阵型CCD，可选用Prosilica EC1600，其感光面尺寸为7.2mm×5.4mm，像素为1620×1220，像元尺寸为4.4μm×4.4μm，帧频为每秒15幅，电子快门10μs～10s，按照每帧画面可以记录两个时间信息的特点，每秒可以记录30幅全息图。
实施例3：如图1、2所示，该三维实时超分辨数字全息记录系统与实施例1相同。所用激光器1为输出激光脉冲中心波长为775nm、脉冲时间间隔小于150fs的Clarke-MXR CPA2010型飞秒激光振荡器，倍频晶体2为可在透过775nm高频激光脉冲的同时、产生另一波长频率为387nm的低频激光脉冲的LBO三硼酸锂晶体（尺寸为：5mm×3mm×8mm，晶体响应波长范围为160～2600nm、倍频效率约为30%）；所产生的激光脉冲光学延迟由一块二向色镜3、一块屋脊棱镜以及一块非偏振分束镜6实现，二向色镜是具有两向色性的分光板，其对波长为387nm的光具有高反射率、对775nm的光具有较高的透过率，通过反射过程延长中心波长387nm的激光脉冲经历的光程，使其滞后于波长为775nm的激光脉冲33ps输出。非偏振分束镜和反射镜的布置使每两束照射光中的一束沿垂直于物体平面的方向照射被测物体，另一束倾斜16o照射被测物体。非偏振分束镜6、7、8、10和15均为分光比为1:1的普通直透光与反射光宽宽带偏振分光棱镜，它能够将波长较长的脉冲激光775nm和波长较短的脉冲激光387nm的入射光分成两束光，但是不改变其偏振方向，保持其原有偏振方向；反射镜9、12、16和17为能对近红外775nm的光和387nm的光具有大于90%反射率的保护银反射镜；显微物镜14为可在记录前对被记录物体进行预先放大的普通小数值孔径、大视场、长工作距离和深焦物镜（NA：0.28,10×），倾斜照射物体的光束与物平面法线的夹角sin-1（NA）=16o；数字相机CCD11为用于记录全息图的普通黑白面阵型CCD，型号为：Basler Scout scA1390-17 单色相机，其感光面尺寸为6.4mm×4.8mm，像素为1392×1040，像元尺寸为4.65μm×4.65μm，帧频为每秒17幅，按照每帧画面可以记录两个时间信息的特点，每秒可以记录34幅全息图。
Claims (8)
1.一种三维实时超分辨数字全息记录系统，其特征在于：装置包括入射光产生系统和记录系统两个部分；入射光产生系统包括飞秒脉冲激光器（1），可按一定比例产生另一波长为原波长一半的脉冲激光的倍频晶体（2），对两波长脉冲激光中某一波长脉冲激光进行反射、同时允许另一波长的脉冲激光按一定比例直接透过的二向色镜（3），用于延长两波长脉冲激光中某一波长脉冲激光光程的光程延迟系统，将长、短两波长脉冲激光汇合成一个方向传播的非偏振分束镜（6）；记录系统包括可将进入该系统的毎束入射脉冲激光分成两束偏正方向分别沿X和Y方向、相互垂直的线性偏振光的偏振分束镜（7），位于偏振分束镜（7）与反射镜（9）和（12）之间、将来自偏振分束镜(7)的沿X方向水平偏正的光束分成参考光A和照射光B的非偏振分束镜(8)，位于非偏振分束镜(8)与被测物体(13)之间、将照射光B反射后沿物体平面法线方向照射被测物体(13)的反射镜(12)，位于偏振分束镜(7)与反射镜 (16)和(17)之间、将来自偏振分束镜(7)的沿Y方向垂直偏正的光分成照射光C和参考光D的非偏振分束镜(15)，位于被测物体(13)与非偏振分束镜(15)之间、将照射光C反射后沿一倾斜方向照射被测物体(13)的反射镜(17)，将参考光A和D反射后射向非偏振分束镜(10)的反射镜(9)和(16)，位于被测物体(13)和非偏振分束镜（10）之间、将物体放大的显微物镜（14）,用于调整来自显微物镜（14）且分别携带物体高、低频信息的物光和来自反射镜（9）、（16）的参考光A和D之间夹角的非偏振分束镜（10），用于记录来自非偏振分束镜（10）的物光和参考光所形成全息图的数字相机CCD（11），以及计算机三维数字全息图像处理再现系统。
2.根据权利要求1所述的三维实时超分辨数字全息记录系统，其特征是：飞秒脉冲激光器（1）为可输出中心波长为775nm～1030nm、单脉冲宽度在飞秒量级的脉冲激光器。
3.根据权利要求1或2所述的三维实时超分辨数字全息记录系统，其特征是：倍频晶体（2）是倍频效率为30%～70%的非线性光学晶体。
4.根据权利要求1所述的三维实时超分辨数字全息记录系统，其特征在于：光程延迟系统包括可使两波长脉冲激光中某一波长脉冲激光到达记录系统的时间滞后或超前于另一波长脉冲激光的反射镜组合（4）和（5），反射镜（4）和（5）分别与二向色镜（3）和非偏振分束镜（6）相对应。
5.根据权利要求4所述的三维实时超分辨数字全息记录系统，其特征在于：光程延迟系统使两波长脉冲激光中某一波长脉冲激光到达记录系统滞后或超前于另一波长脉冲激光的时间差为50fs～33ps。
6.根据权利要求4或5所述的三维实时超分辨数字全息记录系统，其特征是：光程延迟系统中，两块反射镜（4）和（5）可用一块屋脊棱镜代替，反射镜（4）、（5）或屋脊棱镜置于入射照射光源部分的长波长脉冲激光光程或短波长脉冲激光光程中。
7.根据权利要求1所述的三维实时超分辨数字全息记录系统，其特征在于：分别从不同方向照射被测物体的两束照射光B和C中，照射光B沿物体平面的法线方向照射被测物体（13）、照射光C与物平面法线呈sin-1（NA）角度照射被测物体（13），其中，NA是显微物镜（14）的数值孔径。
8.根据权利要求1所述的三维实时超分辨数字全息记录系统，其特征是：非偏振分束镜（6）、（7）、（8）、（10）、（15）均为直透光与反射光分光比为1:1的宽光谱分光棱镜，反射镜（4）、（5）、（9）、（12）、（16）和（17）均为宽光谱、高反射率反射镜，显微物镜（14）为可在记录前对被记录物体进行预先放大的小数值孔径、大视场、长距离和焦深的物镜，数字相机CCD（11）为黑白面阵型数字相机，各非偏振分束镜和反射镜的光谱范围与激光光源的光谱范围相适应。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110389124.XA CN102402172B (zh) | 2011-11-30 | 2011-11-30 | 一种三维实时超分辨数字全息记录系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110389124.XA CN102402172B (zh) | 2011-11-30 | 2011-11-30 | 一种三维实时超分辨数字全息记录系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102402172A true CN102402172A (zh) | 2012-04-04 |
| CN102402172B CN102402172B (zh) | 2015-03-11 |
Family
ID=45884491
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110389124.XA Expired - Fee Related CN102402172B (zh) | 2011-11-30 | 2011-11-30 | 一种三维实时超分辨数字全息记录系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102402172B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102411298A (zh) * | 2011-11-30 | 2012-04-11 | 昆明理工大学 | 一种三维实时超分辨数字全息记录方法 |
| CN103353675A (zh) * | 2013-07-02 | 2013-10-16 | 浙江大学 | 基于微结构的移频超分辨显微成像方法和装置 |
| CN104538827A (zh) * | 2014-12-15 | 2015-04-22 | 中国航空工业集团公司北京航空制造工程研究所 | 一种基于数字全息技术的光束净化方法以及固体激光器 |
| CN106872408A (zh) * | 2017-04-26 | 2017-06-20 | 赣南师范大学 | 一种浮游生物成像检测装置 |
| CN108195761A (zh) * | 2018-03-06 | 2018-06-22 | 南京信息工程大学 | 一种多维可调的分子准直实验系统 |
| CN111819836A (zh) * | 2018-03-14 | 2020-10-23 | 株式会社富士 | 拍摄单元及元件安装机 |
| CN112525104A (zh) * | 2020-12-18 | 2021-03-19 | 昆明理工大学 | 一种数字全息三维形貌测量装置及方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19744368A1 (de) * | 1997-10-08 | 1999-05-20 | Lzh Laserzentrum Hannover Ev | Verfahren und Vorrichtung zur Mikrobearbeitung von Werkstücken mittels Laserstrahlung, insbesondere zum Bilden von im wesentlichen rotationssymmetrischen Ausnehmungen in Werkstücken |
| US6037579A (en) * | 1997-11-13 | 2000-03-14 | Biophotonics Information Laboratories, Ltd. | Optical interferometer employing multiple detectors to detect spatially distorted wavefront in imaging of scattering media |
| EP1162519A1 (en) * | 1999-12-17 | 2001-12-12 | Japan Science and Technology Corporation | Production method and device for hologram |
| CN1542568A (zh) * | 2003-11-07 | 2004-11-03 | 中国科学院上海光学精密机械研究所 | 飞秒脉冲全息装置 |
| CN101788273A (zh) * | 2010-02-05 | 2010-07-28 | 北京航空航天大学 | 一种基于多偏振态合成的数字全息三维显微观测装置 |
| CN101806733A (zh) * | 2010-03-11 | 2010-08-18 | 中国科学院上海光学精密机械研究所 | 飞秒数字全息动态观察测量装置 |
| CN101916570A (zh) * | 2010-07-02 | 2010-12-15 | 中国科学院上海光学精密机械研究所 | 合成孔径飞秒脉冲数字全息记录装置 |
| CN202583708U (zh) * | 2011-11-30 | 2012-12-05 | 昆明理工大学 | 一种三维实时超分辨数字全息记录系统 |
-
2011
- 2011-11-30 CN CN201110389124.XA patent/CN102402172B/zh not_active Expired - Fee Related
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19744368A1 (de) * | 1997-10-08 | 1999-05-20 | Lzh Laserzentrum Hannover Ev | Verfahren und Vorrichtung zur Mikrobearbeitung von Werkstücken mittels Laserstrahlung, insbesondere zum Bilden von im wesentlichen rotationssymmetrischen Ausnehmungen in Werkstücken |
| US6037579A (en) * | 1997-11-13 | 2000-03-14 | Biophotonics Information Laboratories, Ltd. | Optical interferometer employing multiple detectors to detect spatially distorted wavefront in imaging of scattering media |
| EP1162519A1 (en) * | 1999-12-17 | 2001-12-12 | Japan Science and Technology Corporation | Production method and device for hologram |
| CN1542568A (zh) * | 2003-11-07 | 2004-11-03 | 中国科学院上海光学精密机械研究所 | 飞秒脉冲全息装置 |
| CN101788273A (zh) * | 2010-02-05 | 2010-07-28 | 北京航空航天大学 | 一种基于多偏振态合成的数字全息三维显微观测装置 |
| CN101806733A (zh) * | 2010-03-11 | 2010-08-18 | 中国科学院上海光学精密机械研究所 | 飞秒数字全息动态观察测量装置 |
| CN101916570A (zh) * | 2010-07-02 | 2010-12-15 | 中国科学院上海光学精密机械研究所 | 合成孔径飞秒脉冲数字全息记录装置 |
| CN202583708U (zh) * | 2011-11-30 | 2012-12-05 | 昆明理工大学 | 一种三维实时超分辨数字全息记录系统 |
Non-Patent Citations (3)
| Title |
|---|
| 王晓雷: "应用波分复用技术实现飞秒级超快动态过程的脉冲数字显微全息记录", 《光学学报》, vol. 26, no. 11, 30 November 2006 (2006-11-30), pages 1632 - 1635 * |
| 胡浩丰等: "强飞秒激光烧蚀石英玻璃的超快时间分辨光学诊断", 《物理学报》, vol. 60, no. 1, 15 January 2011 (2011-01-15) * |
| 袁操今等: "利用相位模板实现数字全息超分辨成像", 《光子学报》, vol. 39, no. 5, 31 May 2010 (2010-05-31), pages 893 - 896 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102411298B (zh) * | 2011-11-30 | 2014-11-05 | 昆明理工大学 | 一种三维实时超分辨数字全息记录方法 |
| CN102411298A (zh) * | 2011-11-30 | 2012-04-11 | 昆明理工大学 | 一种三维实时超分辨数字全息记录方法 |
| CN103353675B (zh) * | 2013-07-02 | 2015-11-11 | 浙江大学 | 基于微结构的移频超分辨显微成像方法和装置 |
| CN103353675A (zh) * | 2013-07-02 | 2013-10-16 | 浙江大学 | 基于微结构的移频超分辨显微成像方法和装置 |
| CN104538827B (zh) * | 2014-12-15 | 2017-08-04 | 中国航空工业集团公司北京航空制造工程研究所 | 一种基于数字全息技术的光束净化方法以及固体激光器 |
| CN104538827A (zh) * | 2014-12-15 | 2015-04-22 | 中国航空工业集团公司北京航空制造工程研究所 | 一种基于数字全息技术的光束净化方法以及固体激光器 |
| CN106872408A (zh) * | 2017-04-26 | 2017-06-20 | 赣南师范大学 | 一种浮游生物成像检测装置 |
| CN106872408B (zh) * | 2017-04-26 | 2023-05-30 | 赣南师范大学 | 一种浮游生物成像检测装置 |
| CN108195761A (zh) * | 2018-03-06 | 2018-06-22 | 南京信息工程大学 | 一种多维可调的分子准直实验系统 |
| CN108195761B (zh) * | 2018-03-06 | 2023-08-11 | 南京信息工程大学 | 一种多维可调的分子准直实验系统 |
| CN111819836A (zh) * | 2018-03-14 | 2020-10-23 | 株式会社富士 | 拍摄单元及元件安装机 |
| CN111819836B (zh) * | 2018-03-14 | 2021-06-18 | 株式会社富士 | 拍摄单元及元件安装机 |
| CN112525104A (zh) * | 2020-12-18 | 2021-03-19 | 昆明理工大学 | 一种数字全息三维形貌测量装置及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102402172B (zh) | 2015-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102402172B (zh) | 一种三维实时超分辨数字全息记录系统 | |
| CN105241857B (zh) | 一种超分辨成像系统 | |
| Memmolo et al. | Breakthroughs in photonics 2013: holographic imaging | |
| CN105807412B (zh) | 一种基于自由曲面整形的全内反射显微方法与装置 | |
| CN106164784B (zh) | 数字全息设备 | |
| CN101634548B (zh) | 一种光纤-空间光合成孔径的数字全息三维显微观测装置 | |
| CN103926225B (zh) | 一种基于倏逝波照明的荧光发射微分显微方法及装置 | |
| CN101514892B (zh) | 一种基于数字全息的长工作距离原位三维显微观测装置 | |
| CN101788273B (zh) | 一种基于多偏振态合成的数字全息三维显微观测装置 | |
| CN107966110A (zh) | 一种双模态数字全息显微装置 | |
| CN102411298A (zh) | 一种三维实时超分辨数字全息记录方法 | |
| CN107463081B (zh) | 一种利用数字全息技术的三维实时显微成像检测系统 | |
| JP2013531816A5 (zh) | ||
| CN110108201A (zh) | 透反射双模式的高精度离轴数字全息显微装置及成像方法 | |
| CN111250873B (zh) | 一种基于gb-sted的深层超分辨激光直写系统及其实现方法 | |
| CN104457611A (zh) | 双波长剪切干涉数字全息显微测量装置及其方法 | |
| CN106404723B (zh) | 一种二次谐波高分辨率成像方法及系统 | |
| US11762180B2 (en) | Three-photon light sheet imaging | |
| US20230356324A1 (en) | Laser processing system having optical diffraction tomography function | |
| CN201407997Y (zh) | 一种三维显微观测装置 | |
| CN202583708U (zh) | 一种三维实时超分辨数字全息记录系统 | |
| CN107356194B (zh) | 基于二维周期光栅和点衍射的四视场数字全息检测装置与方法 | |
| CN103217209A (zh) | 一种基于数字全息的微振动实时面探测装置 | |
| CN117760559A (zh) | 一种单次曝光定量微分干涉成像装置和方法 | |
| KR102052414B1 (ko) | 반사 위상 현미경 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150311 Termination date: 20201130 |