CN104457611A - 双波长剪切干涉数字全息显微测量装置及其方法 - Google Patents
双波长剪切干涉数字全息显微测量装置及其方法 Download PDFInfo
- Publication number
- CN104457611A CN104457611A CN201410766554.2A CN201410766554A CN104457611A CN 104457611 A CN104457611 A CN 104457611A CN 201410766554 A CN201410766554 A CN 201410766554A CN 104457611 A CN104457611 A CN 104457611A
- Authority
- CN
- China
- Prior art keywords
- laser
- measuring device
- plane
- wavelength
- ccd camera
- 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.)
- Pending
Links
Abstract
本发明公开了一种双波长剪切干涉数字全息显微测量装置及其方法。两台激光器出射的激光从不同方向照射在分束棱镜上,并从分束棱镜的同一面出射,照射在被测物体上的同一点,然后通过显微物镜放大得到两个重合的光斑,再通过平面反射镜将两束激光照射在平面平晶上并在其前后表面发生反射,反射光垂直照射在CCD相机上,用CCD相机记录单束光的全息图。使用Matlab软件对全息图做傅里叶变换,取出两种波长光正一级的频谱,再对提取出来的信息做逆傅里叶变换得到全息图,使用角谱法再现,提取相位,通过相位相减并做进一步运算求得测物体的三维高度信息。本发明具有光路简捷、系统稳定性好、再现程序简便等特点,适用于生物医学以及材料科学等领域。
Description
技术领域
本发明属于数字全息显微及成像技术领域,具体涉及一种基于剪切干涉与双波长光学相位解包裹的数字全息显微装置及其方法,用于测量分散排列的微小物体三维高度信息。
背景技术
数字全息技术是计算机技术和传统光学全息相结合的产物,采用数字化的记录和再现方法,能够较为方便地获取物体的振幅信息和相位信息,而相位信息是恢复物体三维形貌的重要参数。
传统的显微技术主要是通过光学显微镜来观察微小物体。但这只能得到光强的分布,且光学显微镜的成像焦深短,每次聚焦只能获取物体的二维图像,对被测物体的厚度、位置以及多层细胞的观测非常困难;光学断面显微镜和激光扫描共焦显微镜可以进行三维观测,但由于是采用机械调节光学系统或被测物体来调焦,致使观测系统非常复杂,且对样品厚度和位置的测量麻烦费时;电子显微镜由于价格昂贵、不能观测活细胞、样品制作困难、测量过程比较麻烦等原因,限制了它的应用。数字全息显微技术是数字全息技术与显微技术的融合,成功的解决了这些问题,能同时获得细胞位置、数量、相位以及三维形貌信息。
数字全息显微成像主要分两步进行。首先利用光学显微镜对被测物体预放大,并用CCD记录数字化全息图。然后数字再现物体的三维图像信息。当前数字全息显微成像光路多采用马赫泽德干涉光路,由于复杂的光路和繁琐的相位解包裹算法而导致其难度较大,特别是光路中分束棱镜的运用使物光与参考光的光程与光强不同,这为干涉系统带入了更多的噪声,降低了系统的稳定性。
发明内容
本发明提出了一种双波长剪切干涉数字全息显微测量装置及其方法,利用剪切干涉共光路自干涉的特点,使光路得到简化并增强了系统的稳定性。
本发明的装置采用的技术方案如下:
双波长剪切干涉数字全息显微测量装置,包括第一激光器、第二激光器、分束棱镜、被测物体、显微物镜、平面反射镜、平面平晶和CCD相机,第一激光器出射的第一激光与第二激光器出射的第二激光从不同方向照射在分束棱镜上,第一激光的反射光与第二激光的透射光从分束棱镜的同一面出射,照射在被测物体上的同一点;经过被测物体的两束激光通过显微物镜放大,得到两个重合的光斑,再通过平面反射镜将两束激光照射在平面平晶上,这两束激光在平面平晶的前后表面发生反射后,反射光垂直照射在CCD相机上。
所述被测物体为分布在载玻片上的直径为1.97微米的二氧化硅微球。
进一步地,所述平面平晶的厚度为3毫米。
本发明利用上述装置的测量方法,包括如下步骤:
a)先后挡住两个激光器出射的两束光中的一束光,分别用CCD相机记录单束光的全息图;
b)双波长频谱分离:将CCD相机采集的全息图读入Matlab软件中,分别对两幅全息图做傅里叶变换,然后提取出两种波长光正一级的频谱,再对其做逆傅里叶变换得到全息图;
c)数字再现算法:用角谱再现算法对全息图进行再现得到物体的振幅信息;
d)再现像的相位提取和修正:对再现的物体像进行相位提取,并对提取的包裹相位进行滤波;
e)相位相减:对两波长的相位进行相减,对小于零的相位差加上一个2π,使相位差调整到[0,2π]范围;
f)得到物体三维高度图:根据相位与物体高度之间的关系,计算得到物体三维高度。
本发明的方法以数字全息技术为基础,利用剪切干涉的共光路自干涉特性将同一束光分成物光与参考光,并采用双波长光学相位解包裹方法代替相位解包裹算法,实现微小物体或生物细胞的三维形貌信息的探测。本发明最大的特点就是将平面平晶引入的剪切干涉与双光束相结合,其具有以下优点:
(1)平面平晶的使用有效的简化了光路,减少了光路搭建的时间与难度以及光学元件的使用;
(2)剪切干涉共光路自干涉的特性,使得干涉条纹对比度较好、系统噪声相对较少、稳定性较强;
(3)双光束的使用降低了编程的难度,实现了双波长光学相位解包裹;
(4)本发明具有光路简捷、系统稳定性好、噪声低、再现程序简便等特点,适用于生物医学以及材料科学等领域。
附图说明
图1是本发明具体光路。其中,1-第一激光器,2-第二激光器,3-分束棱镜,4-被测物体,5-显微物镜,6-平面反射镜,7-平面平晶,8-CCD相机。
图2是横向剪切干涉原理,A与B为两个大小相同具有一定错位的波面。
图3是获得物体三维高度的处理流程图。
具体实施方式
如图1所示,本发明的光路由第一激光器1、第二激光器2、分束棱镜3、被测物体4、显微物镜5、平面反射镜6、平面平晶7和CCD相机8组成。分束棱镜3用来改变第一激光器1出射光的方向,使两束光从同一方向出射;被测物体4为稀疏分布在载玻片上的直径为1.97微米的二氧化硅微球,平面反射镜6将第一激光器1与第二激光器2出射的光照射到厚度为3毫米的平面平晶7上,在平面平晶7前表面和后表面上分别发生反射,并照射在CCD相机8上。
剪切干涉是一种重要的光学测量方法,它是用光学装置把一个具有空间相干性的波前分成相似或相同的两个波前,并使这两个波前产生相对错位,错位后的两个波前在空间公共区域产生剪切干涉条纹。其主要特点是两光波共光路,对机械振动、温度扰动以及空气流动都不敏感,便于在实验室为进行现场检查。剪切干涉根据剪切方向可以分成横向剪切、径向剪切、旋转剪切和反向剪切等。本发明采用的是横向剪切,利用平面平晶的前后表面分别发生反射,形成两个大小相同相对错位的波面,在这两个波面的公共区域产生干涉条纹,干涉原理如图2所示。
普通的数字全息显微技术由于光波波长的限制再现后得到的相位信息是包裹的,其值分布在[-π,π],包裹相位中存在阶跃点,而且阶跃点出现在高度等于波长整数倍处,波长不同则阶跃点出现的区域也不同。要得到真实的相位信息就需要采用电脑编程来完成解包裹这一步骤,现如今的相位解包裹算法众多,各种算法适用范围不同,尚无一种广泛适用的相位解包裹算法,而且这些算法相对较为复杂,并不可避免的引入了一些噪声。本发明采用双波长光学相位解包裹技术可以有效的解决这一问题,很大程度上简化了再现过程,不需要对再现像相位图进行解包裹,并可以实时重构物体形貌。其原理如下:
假设λ1,λ2为两相差甚小的波长,对透射式光路该物体的波前相位分别为
将波长λ1,λ2记录所得物体包裹相位Φ1、Φ2相减可得Φ12:
Λ=λ1λ2/|λ1-λ2| (3)
相位Φ12可看作是波长为Λ的相位,Λ称之为综合波长,大于任意一个原始波长。
Φ12分布范围为[-2π,2π],也存在2π阶跃点,但是在Φ12小于零的值加上2π,即可得Φ。
Φ=Φ12+2π(Φ12<0) (4)
由于本发明的剪切干涉参考光与物光的夹角由平面平晶厚度确定,故难以调整不同波长的频谱位置,所以为了简便起见,将两种波长的全息图分别拍摄。具体步骤如下:
第一激光器1出射的第一激光与第二激光器2出射的第二激光从不同方向照射在分束棱镜3上。第一激光的反射光与第二激光透射光从分束棱镜3的同一面出射,照射在被测物体4的同一点上。两束激光通过显微物镜5放大,得到两个重合的光斑。通过平面反射镜将两束光照射在平面平晶上,使之在平面平晶前后表面发生反射,并使反射光垂直照射在CCD相机8上。最后挡住一束光,分别用CCD相机8记录单束光的全息图。
将CCD相机采集的全息图读入Matlab中,分别对其做傅里叶变换,取出两种波长光的正一级频谱,再对提取出来的信息做逆傅里叶变换得到全息图,使用角谱法再现,提取相位,通过相位相减并做进一步运算求的被测物体的三维高度信息。具体流程如图3所示:
1.双波长频谱分离:分别对全息图进行傅里叶变换,然后提取出两种波长光正一级的频谱,再对其做逆傅里叶变换得到全息图;
2.数字再现算法:用角谱再现算法对全息图进行再现得到物体的振幅信息;
3.再现像的相位提取和修正:对再现像进行相位提取,并对提取的包裹相位进行滤波;
4.相位相减:对两波长的相位进行相减,对小于零的相位差加上一个2π,使相位差调整到[0,2π]范围;
5.得到物体三维高度图:根据相位与物体高度之间的关系,计算得到物体三维高度。
Claims (4)
1.双波长剪切干涉数字全息显微测量装置,包括第一激光器(1)、第二激光器(2)、分束棱镜(3)、被测物体(4)、显微物镜(5)、平面反射镜(6)、平面平晶(7)和CCD相机(8),其特征在于,第一激光器(1)出射的第一激光与第二激光器(2)出射的第二激光从不同方向照射在分束棱镜(3)上,第一激光的反射光与第二激光的透射光从分束棱镜(3)的同一面出射,照射在被测物体(4)上的同一点;经过被测物体(4)的两束激光通过显微物镜(5)放大,得到两个重合的光斑,再通过平面反射镜(6)将两束激光照射在平面平晶(7)上,这两束激光在平面平晶(7)的前后表面发生反射后,反射光垂直照射在CCD相机(8)上。
2.根据权利要求1所述的双波长剪切干涉数字全息显微测量装置,其特征在于,所述被测物体(4)为分布在载玻片上的直径为1.97微米的二氧化硅微球。
3.根据权利要求1或2所述的双波长剪切干涉数字全息显微测量装置,其特征在于,所述平面平晶(7)的厚度为3毫米。
4.利用如权1所述双波长剪切干涉数字全息显微测量装置的测量方法,其特征在于,包括如下步骤:
a)先后挡住两个激光器出射的两束光中的一束光,分别用CCD相机(8)记录单束光的全息图;
b)双波长频谱分离:将CCD相机(8)采集的全息图读入Matlab软件中,分别对两幅全息图做傅里叶变换,然后提取出两种波长光正一级的频谱,再对其做逆傅里叶变换得到全息图;
c)数字再现算法:用角谱再现算法对全息图进行再现得到物体的振幅信息;
d)再现像的相位提取和修正:对再现的物体像进行相位提取,并对提取的包裹相位进行滤波;
e)相位相减:对两波长的相位进行相减,对小于零的相位差加上一个2π,使相位差调整到[0,2π]范围;
f)得到物体三维高度图:根据相位与物体高度之间的关系,计算得到物体三维高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410766554.2A CN104457611A (zh) | 2014-12-11 | 2014-12-11 | 双波长剪切干涉数字全息显微测量装置及其方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410766554.2A CN104457611A (zh) | 2014-12-11 | 2014-12-11 | 双波长剪切干涉数字全息显微测量装置及其方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104457611A true CN104457611A (zh) | 2015-03-25 |
Family
ID=52904049
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410766554.2A Pending CN104457611A (zh) | 2014-12-11 | 2014-12-11 | 双波长剪切干涉数字全息显微测量装置及其方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104457611A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107065492A (zh) * | 2017-05-31 | 2017-08-18 | 南京师范大学 | 一种无透镜彩色全息投影方法 |
CN107677201A (zh) * | 2017-08-21 | 2018-02-09 | 深圳大学 | 基于波前旋转的快速倾斜相位误差补偿方法及装置 |
CN108007379A (zh) * | 2017-11-20 | 2018-05-08 | 西北工业大学 | 一种双波长共路数字全息干涉测量方法与系统 |
CN109298411A (zh) * | 2018-11-30 | 2019-02-01 | 武汉理工大学 | 一种全息激光信号声呐滤波方法 |
CN111829453A (zh) * | 2020-06-28 | 2020-10-27 | 西安工业大学 | 一种共光路数字全息显微测量装置及其测量方法 |
CN112833781A (zh) * | 2020-12-31 | 2021-05-25 | 佛山英智莱科技有限公司 | 一种双光源线结构光传感器 |
CN114001646A (zh) * | 2021-10-28 | 2022-02-01 | 山西大学 | 三波长线差分共焦显微探测方法与装置 |
CN114518181A (zh) * | 2022-01-25 | 2022-05-20 | 北京航空航天大学 | 一种基于波长调制剪切干涉的火焰温度分布测量系统与方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1524491A1 (en) * | 2003-10-16 | 2005-04-20 | Universite Libre De Bruxelles | Apparatus coupling an interferometer and a microscope |
CN200970231Y (zh) * | 2006-11-27 | 2007-11-07 | 浙江大学 | 一种扩展光学相干层析成像动态范围的系统 |
CN103148800A (zh) * | 2013-01-28 | 2013-06-12 | 浙江大学 | 一种基于光场传播的非标记三维显微方法和装置 |
CN103335592A (zh) * | 2013-06-21 | 2013-10-02 | 南京师范大学 | 双洛埃镜数字全息显微测量方法 |
-
2014
- 2014-12-11 CN CN201410766554.2A patent/CN104457611A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1524491A1 (en) * | 2003-10-16 | 2005-04-20 | Universite Libre De Bruxelles | Apparatus coupling an interferometer and a microscope |
CN200970231Y (zh) * | 2006-11-27 | 2007-11-07 | 浙江大学 | 一种扩展光学相干层析成像动态范围的系统 |
CN103148800A (zh) * | 2013-01-28 | 2013-06-12 | 浙江大学 | 一种基于光场传播的非标记三维显微方法和装置 |
CN103335592A (zh) * | 2013-06-21 | 2013-10-02 | 南京师范大学 | 双洛埃镜数字全息显微测量方法 |
Non-Patent Citations (2)
Title |
---|
AMARDEEP S.G. SINGH ET.AL: "Lateral shearing digital holographic imaging of small biological specimens", 《OPTICAL EXPRESS》 * |
袁飞等: "双Lloyd镜数字全息显微测量术", 《物理学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107065492A (zh) * | 2017-05-31 | 2017-08-18 | 南京师范大学 | 一种无透镜彩色全息投影方法 |
CN107065492B (zh) * | 2017-05-31 | 2019-04-12 | 南京师范大学 | 一种无透镜彩色全息投影方法 |
CN107677201A (zh) * | 2017-08-21 | 2018-02-09 | 深圳大学 | 基于波前旋转的快速倾斜相位误差补偿方法及装置 |
CN108007379A (zh) * | 2017-11-20 | 2018-05-08 | 西北工业大学 | 一种双波长共路数字全息干涉测量方法与系统 |
CN109298411A (zh) * | 2018-11-30 | 2019-02-01 | 武汉理工大学 | 一种全息激光信号声呐滤波方法 |
CN111829453A (zh) * | 2020-06-28 | 2020-10-27 | 西安工业大学 | 一种共光路数字全息显微测量装置及其测量方法 |
CN112833781A (zh) * | 2020-12-31 | 2021-05-25 | 佛山英智莱科技有限公司 | 一种双光源线结构光传感器 |
CN114001646A (zh) * | 2021-10-28 | 2022-02-01 | 山西大学 | 三波长线差分共焦显微探测方法与装置 |
CN114001646B (zh) * | 2021-10-28 | 2023-12-26 | 山西大学 | 三波长线差分共焦显微探测方法与装置 |
CN114518181A (zh) * | 2022-01-25 | 2022-05-20 | 北京航空航天大学 | 一种基于波长调制剪切干涉的火焰温度分布测量系统与方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104457611A (zh) | 双波长剪切干涉数字全息显微测量装置及其方法 | |
JP7161777B2 (ja) | ホログラフィック撮像装置および同装置に用いるデータ処理方法 | |
EP3065001B1 (en) | Holographic microscope and data processing method for high-resolution hologram image | |
CN109709786B (zh) | 一种超分辨率数字全息成像系统和成像方法 | |
CN105204310B (zh) | 基于光学扫描全息技术的复振幅物体重建装置 | |
Heimbeck et al. | Terahertz digital holographic imaging | |
CN103149827B (zh) | 消除单光束同轴数字全息直流项和共轭像的方法 | |
Zuo et al. | Direct continuous phase demodulation in digital holography with use of the transport-of-intensity equation | |
CN108614405B (zh) | 一种基于环形光差频扫描实现边缘提取的全息方法 | |
JP7352292B2 (ja) | ホログラフィック撮像装置およびホログラフィック撮像方法 | |
CN103034109A (zh) | 双ccd镜像重叠调节及单曝光同轴数字全息记录装置 | |
Di et al. | Improvement of measurement accuracy in digital holographic microscopy by using dual-wavelength technique | |
Abdelsalam et al. | Real-time dual-wavelength digital holographic microscopy based on polarizing separation | |
JPWO2018038064A1 (ja) | エリプソメトリ装置およびエリプソメトリ方法 | |
CN103335592B (zh) | 双洛埃镜数字全息显微测量方法 | |
CN114111625A (zh) | 一种光学成像系统及方法 | |
CN109283821A (zh) | 基于涡旋透镜的相移数字全息单次曝光成像装置及方法 | |
Bordbar et al. | Application of complex field imaging sensor to additive manufacturing | |
JP6230358B2 (ja) | 光学的距離計測装置 | |
CN114324245B (zh) | 基于部分相干结构光照明的定量相位显微装置和方法 | |
JP6629649B2 (ja) | 観察装置および観察方法 | |
Pedrini et al. | Spatial filtering in digital holographic microscopy | |
Huang et al. | Non-interferometric accurate phase imaging via linear-convergence iterative optimization | |
CN114965477B (zh) | 全自动生物信息分析取样系统及方法 | |
CN117705001B (zh) | 基于凹面镜的数字全息物体形貌测量装置及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150325 |
|
RJ01 | Rejection of invention patent application after publication |