CN110260780A - 一种数字全息显微相位误差同时补偿方法及装置 - Google Patents
一种数字全息显微相位误差同时补偿方法及装置 Download PDFInfo
- Publication number
- CN110260780A CN110260780A CN201910531873.8A CN201910531873A CN110260780A CN 110260780 A CN110260780 A CN 110260780A CN 201910531873 A CN201910531873 A CN 201910531873A CN 110260780 A CN110260780 A CN 110260780A
- Authority
- CN
- China
- Prior art keywords
- phase
- compensation
- original
- error
- obtains
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02041—Interferometers characterised by particular imaging or detection techniques
- G01B9/02047—Interferometers characterised by particular imaging or detection techniques using digital holographic imaging, e.g. lensless phase imaging without hologram in the reference path
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/04—Measuring microscopes
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/36—Microscopes arranged for photographic purposes or projection purposes or digital imaging or video purposes including associated control and data processing arrangements
- G02B21/365—Control or image processing arrangements for digital or video microscopes
- G02B21/367—Control or image processing arrangements for digital or video microscopes providing an output produced by processing a plurality of individual source images, e.g. image tiling, montage, composite images, depth sectioning, image comparison
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Optics & Photonics (AREA)
- Holo Graphy (AREA)
Abstract
本发明公开了一种数字全息显微相位误差同时补偿方法,属于光学成像校正技术领域,解决的是现有方法计算量大、不适合动态相位测量的技术问题,所述方法包括步骤如下,S1、采集待测物体的全息图作为原始全息图;S2、对原始全息图进行傅立叶变换处理得到原始展开相位S3、对原始展开相位作旋转变换得到旋转相位将原始展开相位加上旋转相位进行离轴倾斜误差补偿并得到补偿相位S4、对补偿相位作翻转变换得到翻转相位将补偿相位减去翻转相位进行二次相位误差补偿并得到最终相位本发明还公开了一种实现上述方法的数字全息显微相位误差同时补偿装置。本发明能简单且准确去除离轴倾斜误差和二次相位误差。
Description
技术领域
本发明涉及光学成像校正技术领域,更具体地说,它涉及一种数字全息显微相位误差同时补偿方法及装置。
背景技术
在传统离轴数字全息显微系统中,离轴角度造成恢复相位存在离轴倾斜误差,另外显微物镜的使用也会对恢复相位引入二次相位误差,为了得到准确的相位分布,需要同时去除这些相位误差。
补偿二次相位误差可以采用物理补偿法和数值补偿法。然而,物理补偿法不能去除离轴倾斜误差。补偿离轴倾斜误差可以采用频谱中心法和数值补偿法,频谱中心法也不能去除二次相位误差。因此,我们可以采用数值补偿法来同时去除离轴倾斜误差和二次相位误差。数值补偿法的思想大都是通过提前获取物体或系统预先知识、数值拟合算法、深度学习训练等构建相位补偿模板。然而,这些方法存在计算量大、对噪声敏感等问题。二次曝光法可以同时去除这两种误差,但是需要多次采集全息图,不适合动态相位测量。
发明内容
本发明要解决的技术问题是针对现有技术的上述不足,本发明的目的一是提供一种简单且准确地去除离轴倾斜误差和二次相位误差的数字全息显微相位误差同时补偿方法。
本发明的目的二是提供一种简单且准确地去除离轴倾斜误差和二次相位误差的数字全息显微相位误差同时补偿装置。
为了实现上述实目的一,本发明提供一种数字全息显微相位误差同时补偿方法,包括步骤如下,
S1、利用数字全息显微系统采集待测物体的全息图,并作为原始全息图;
S2、对所述原始全息图进行傅立叶变换得到所述原始全息图的频谱,并对所述原始全息图的频谱提取+1级频谱,进行逆傅立叶变换得到包裹相位,通过解包裹算法得到原始展开相位
S3、对所述原始展开相位作旋转变换得到旋转相位将所述原始展开相位加上所述旋转相位进行离轴倾斜误差补偿,并得到补偿相位
S4、对所述补偿相位作翻转变换得到翻转相位将所述补偿相位减去所述翻转相位进行二次相位误差补偿,并得到最终相位
作为进一步地改进,在所述步骤S2中,所述原始全息图+1级像的公式为:
其中T(x,y)为离轴倾斜误差,T(x,y)=kxx+kyy,P(x,y)为二次相位误差,P(x,y)=k(x2+y2)/r,为原始物体相位,所述原始展开相位的公式为:
进一步地,在所述步骤S3中,对所述原始展开相位旋转180°得到旋转相位所述旋转相位的公式为:
所述补偿相位的公式为:
进一步地,在所述步骤S4中,对所述补偿相位作上下翻转或者左右翻转变换得到翻转相位所述翻转相位的公式为:
或
所述最终相位的公式为:
或
为了实现上述实目的二,本发明提供一种数字全息显微相位误差同时补偿装置,包括:
采集模块,用于采集待测物体的全息图并作为原始全息图,所述采集模块包括数字全息显微系统;
傅立叶变换模块,用于对所述原始全息图进行傅立叶变换得到所述原始全息图的频谱,并对所述原始全息图的频谱提取+1级频谱,进行逆傅立叶变换得到包裹相位,通过解包裹算法得到原始展开相位
第一补偿模块,用于对所述原始展开相位作旋转变换得到旋转相位将所述原始展开相位加上所述旋转相位进行离轴倾斜误差补偿,并得到补偿相位
第二补偿模块,用于对所述补偿相位作翻转变换得到翻转相位将所述补偿相位减去所述翻转相位进行二次相位误差补偿,并得到最终相位
有益效果
本发明与现有技术相比,具有的优点为:本发明通过原始全息图经过傅立叶变换处理得到原始展开相位,原始展开相位旋转180°得到旋转相位,原始展开相位加上旋转相位可以去除离轴倾斜误差并得到补偿相位,补偿相位作上下翻转或者左右翻转变换得到翻转相位,补偿相位减去翻转相位可以去除二次相位误差并得到最终相位;本方法可以同时有效去除离轴倾斜误差和二次相位误差,计算简单且快速,避免了物体或系统预先知识的获取或复杂的数值拟合计算问题,只需单幅全息图即可,适用于动态相位测量。
附图说明
图1为本发明中方法的方框图;
图2为本发明中装置的方框图。
具体实施方式
下面结合附图中的具体实施例对本发明做进一步的说明。
参阅图1、2,一种数字全息显微相位误差同时补偿方法,包括步骤如下,
S1、利用数字全息显微系统采集待测物体的全息图,并作为原始全息图;
S2、对原始全息图进行傅立叶变换得到原始全息图的频谱,并对原始全息图的频谱提取+1级频谱,进行逆傅立叶变换得到包裹相位,通过解包裹算法得到原始展开相位
S3、对原始展开相位作旋转变换得到旋转相位将原始展开相位加上旋转相位进行离轴倾斜误差补偿,并得到补偿相位
S4、对补偿相位作翻转变换得到翻转相位将补偿相位减去翻转相位进行二次相位误差补偿,并得到最终相位
在步骤S2中,原始全息图+1级像的公式为:
其中,T(x,y)为离轴倾斜误差,T(x,y)=kxx+kyy,P(x,y)为二次相位误差,P(x,y)=k(x2+y2)/r,为原始物体相位,Iv指的是提取原始全息图频谱+1级像恢复的强度图像,v是virtual image的简称,其中virtual image其实就是+1级像,|O|是物光波光强,|R|是参考光光波强度,i就是一个虚数单位,kx和ky指的是在X方向和Y方向的倾斜因子,描述的是倾斜程度,k为波数k=2π/λ,λ为波长,r是二次相位误差分布的曲率,原始展开相位的公式为:
在步骤S3中,对所述原始展开相位旋转180°得到旋转相位旋转相位的公式为:
补偿相位的公式为:
在步骤S4中,对所述补偿相位作上下翻转或者左右翻转变换得到翻转相位翻转相位的公式为:
或
最终相位的公式为:
或
当然,也可以先对原始展开相位作翻转变换得到翻转相位,原始展开相位减去翻转相位进行二次相位误差补偿并得到补偿相位,再对补偿相位作旋转变换得到旋转相位,补偿相位加上旋转相位进行离轴倾斜误差补偿并得到最终相位。
根据最终相位的公式可知二次相位误差和离轴倾斜误差都补偿了,最终相位除了原始相位分布,还多了旋转原始相位或翻转原始相位。当原始物体位置不超过相机视场1/4范围内,这些相位分布可以独立分开来,不影响原始相位分布。
一种数字全息显微相位误差同时补偿装置,包括:
采集模块,用于采集待测物体的全息图并作为原始全息图,所述采集模块包括数字全息显微系统;
傅立叶变换模块,用于对所述原始全息图进行傅立叶变换得到所述原始全息图的频谱,并对所述原始全息图的频谱提取+1级频谱,进行逆傅立叶变换得到包裹相位,通过解包裹算法得到原始展开相位
第一补偿模块,用于对所述原始展开相位作旋转变换得到旋转相位将所述原始展开相位加上所述旋转相位进行离轴倾斜误差补偿,并得到补偿相位
第二补偿模块,用于对所述补偿相位作翻转变换得到翻转相位将所述补偿相位减去所述翻转相位进行二次相位误差补偿,并得到最终相位
本发明通过原始全息图经过傅立叶变换处理得到原始展开相位,原始展开相位旋转180°得到旋转相位,原始展开相位加上旋转相位可以去除离轴倾斜误差并得到补偿相位,补偿相位作上下翻转或者左右翻转变换得到翻转相位,补偿相位减去翻转相位可以去除二次相位误差并得到最终相位;本方法可以同时有效去除离轴倾斜误差和二次相位误差,计算简单且快速,避免了物体或系统预先知识的获取或复杂的数值拟合计算问题,只需单幅全息图即可,适用于动态相位测量。
以上仅是本发明的优选实施方式,应当指出对于本领域的技术人员来说,在不脱离本发明结构的前提下,还可以作出若干变形和改进,这些都不会影响本发明实施的效果和专利的实用性。
Claims (5)
1.一种数字全息显微相位误差同时补偿方法,其特征在于:包括步骤如下,
S1、利用数字全息显微系统采集待测物体的全息图,并作为原始全息图;
S2、对所述原始全息图进行傅立叶变换得到所述原始全息图的频谱,并对所述原始全息图的频谱提取+1级频谱,进行逆傅立叶变换得到包裹相位,通过解包裹算法得到原始展开相位
S3、对所述原始展开相位作旋转变换得到旋转相位将所述原始展开相位加上所述旋转相位进行离轴倾斜误差补偿,并得到补偿相位
S4、对所述补偿相位作翻转变换得到翻转相位将所述补偿相位减去所述翻转相位进行二次相位误差补偿,并得到最终相位
2.根据权利要求1所述的一种数字全息显微相位误差同时补偿方法,其特征在于:在所述步骤S2中,所述原始全息图+1级像的公式为:
其中T(x,y)为离轴倾斜误差,T(x,y)=kxx+kyy,P(x,y)为二次相位误差,P(x,y)=k(x2+y2)/r,为原始物体相位,所述原始展开相位的公式为:
3.根据权利要求2所述的一种数字全息显微相位误差同时补偿方法,其特征在于:在所述步骤S3中,对所述原始展开相位旋转180°得到旋转相位所述旋转相位的公式为:
所述补偿相位的公式为:
4.根据权利要求3所述的一种数字全息显微相位误差同时补偿方法,其特征在于:在所述步骤S4中,对所述补偿相位作上下翻转或者左右翻转变换得到翻转相位所述翻转相位的公式为:
或
所述最终相位的公式为:
或
5.一种实现权利要求1所述方法的一种数字全息显微相位误差同时补偿装置,其特征在于:包括,
采集模块,用于采集待测物体的全息图并作为原始全息图,所述采集模块包括数字全息显微系统;
傅立叶变换模块,用于对所述原始全息图进行傅立叶变换得到所述原始全息图的频谱,并对所述原始全息图的频谱提取+1级频谱,进行逆傅立叶变换得到包裹相位,通过解包裹算法得到原始展开相位
第一补偿模块,用于对所述原始展开相位作旋转变换得到旋转相位将所述原始展开相位加上所述旋转相位进行离轴倾斜误差补偿,并得到补偿相位
第二补偿模块,用于对所述补偿相位作翻转变换得到翻转相位将所述补偿相位减去所述翻转相位进行二次相位误差补偿,并得到最终相位
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910531873.8A CN110260780B (zh) | 2019-06-19 | 2019-06-19 | 一种数字全息显微相位误差同时补偿方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910531873.8A CN110260780B (zh) | 2019-06-19 | 2019-06-19 | 一种数字全息显微相位误差同时补偿方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110260780A true CN110260780A (zh) | 2019-09-20 |
CN110260780B CN110260780B (zh) | 2020-12-22 |
Family
ID=67919416
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910531873.8A Active CN110260780B (zh) | 2019-06-19 | 2019-06-19 | 一种数字全息显微相位误差同时补偿方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110260780B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114001643A (zh) * | 2021-09-27 | 2022-02-01 | 上海工程技术大学 | 一种数字全息显微相位畸变补偿方法及装置 |
CN114593671A (zh) * | 2022-01-25 | 2022-06-07 | 华南师范大学 | 基于平行平板的植入式同轴与离轴数字全息切换装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101866659A (zh) * | 2009-04-17 | 2010-10-20 | 索尼公司 | 全息设备,倾斜检测方法和倾斜校正方法 |
CN103322941A (zh) * | 2013-07-09 | 2013-09-25 | 河北工程大学 | 一种准确获取三维显微图像的方法 |
CN104331616A (zh) * | 2014-10-28 | 2015-02-04 | 南京理工大学 | 基于求解光强传输方程的数字全息图解调方法 |
US20150077818A1 (en) * | 2013-09-17 | 2015-03-19 | City University Of Hong Kong | Converting complex holograms to phase holograms |
CN105157561A (zh) * | 2015-09-29 | 2015-12-16 | 南京理工大学 | 基于降采样相位畸变主成分分析的快速相位畸变补偿方法 |
CN107036786A (zh) * | 2017-04-13 | 2017-08-11 | 北京航空航天大学 | 一种基于数字全息层析的光纤熔接点三维结构检测方法 |
CN107167095A (zh) * | 2017-07-10 | 2017-09-15 | 天津农学院 | 基于参考透镜法的数字全息显微相位畸变校正方法 |
CN107677201A (zh) * | 2017-08-21 | 2018-02-09 | 深圳大学 | 基于波前旋转的快速倾斜相位误差补偿方法及装置 |
CN109612384A (zh) * | 2018-11-01 | 2019-04-12 | 南京理工大学 | 一种基于频谱亚像素平移的倾斜像差校正补偿方法 |
CN109724511A (zh) * | 2019-01-29 | 2019-05-07 | 嘉应学院 | 一种数字全息显微的二次相位误差补偿方法及其装置 |
-
2019
- 2019-06-19 CN CN201910531873.8A patent/CN110260780B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101866659A (zh) * | 2009-04-17 | 2010-10-20 | 索尼公司 | 全息设备,倾斜检测方法和倾斜校正方法 |
CN103322941A (zh) * | 2013-07-09 | 2013-09-25 | 河北工程大学 | 一种准确获取三维显微图像的方法 |
US20150077818A1 (en) * | 2013-09-17 | 2015-03-19 | City University Of Hong Kong | Converting complex holograms to phase holograms |
CN104331616A (zh) * | 2014-10-28 | 2015-02-04 | 南京理工大学 | 基于求解光强传输方程的数字全息图解调方法 |
CN105157561A (zh) * | 2015-09-29 | 2015-12-16 | 南京理工大学 | 基于降采样相位畸变主成分分析的快速相位畸变补偿方法 |
CN107036786A (zh) * | 2017-04-13 | 2017-08-11 | 北京航空航天大学 | 一种基于数字全息层析的光纤熔接点三维结构检测方法 |
CN107167095A (zh) * | 2017-07-10 | 2017-09-15 | 天津农学院 | 基于参考透镜法的数字全息显微相位畸变校正方法 |
CN107677201A (zh) * | 2017-08-21 | 2018-02-09 | 深圳大学 | 基于波前旋转的快速倾斜相位误差补偿方法及装置 |
CN109612384A (zh) * | 2018-11-01 | 2019-04-12 | 南京理工大学 | 一种基于频谱亚像素平移的倾斜像差校正补偿方法 |
CN109724511A (zh) * | 2019-01-29 | 2019-05-07 | 嘉应学院 | 一种数字全息显微的二次相位误差补偿方法及其装置 |
Non-Patent Citations (5)
Title |
---|
DINGNAN DENG等: "Off-axis tilt compensation in common-path digital holographic microscopy based on hologram rotation", 《OPTICS LETTERS》 * |
DINGNAN DENG等: "Simple and flexible phase compensation for digital holographic microscopy with electrically tunable lens", 《APPLIED OPTICS 》 * |
ETIENNE CUCHE等: "Digital holography for quantitative phase-contrast imaging", 《OPTICS LETTERS》 * |
HUAKUN CUI等: "Phase aberration compensation by spectrum centering in digital holographic microscopy", 《OPTICS COMMUNICATIONS》 * |
周文静等: "数字显微全息中二次项相位误差的补偿", 《光子学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114001643A (zh) * | 2021-09-27 | 2022-02-01 | 上海工程技术大学 | 一种数字全息显微相位畸变补偿方法及装置 |
CN114001643B (zh) * | 2021-09-27 | 2023-07-25 | 上海工程技术大学 | 一种数字全息显微相位畸变补偿方法及装置 |
CN114593671A (zh) * | 2022-01-25 | 2022-06-07 | 华南师范大学 | 基于平行平板的植入式同轴与离轴数字全息切换装置 |
CN114593671B (zh) * | 2022-01-25 | 2023-10-20 | 华南师范大学 | 基于平行平板的植入式同轴与离轴数字全息切换装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110260780B (zh) | 2020-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yan et al. | Fringe pattern denoising based on deep learning | |
CN110260780A (zh) | 一种数字全息显微相位误差同时补偿方法及装置 | |
CN105157561B (zh) | 基于降采样相位畸变主成分分析的快速相位畸变补偿方法 | |
CN107677201A (zh) | 基于波前旋转的快速倾斜相位误差补偿方法及装置 | |
EP3065001A1 (en) | Holographic microscope and data processing method for high-resolution hologram image | |
CN104331616B (zh) | 基于求解光强传输方程的数字全息图解调方法 | |
CN108319009A (zh) | 基于结构光调制的快速超分辨成像方法 | |
Zhong et al. | Fast phase retrieval in slightly off-axis digital holography | |
CN107395933B (zh) | 一种基于lcos空间光调制器的可编程孔径成像系统及超分辨方法 | |
CN111402395A (zh) | 一种基于cnn校正的被动偏振三维重建方法 | |
Wang et al. | Phase singularities in analytic signal of white-light speckle pattern with application to micro-displacement measurement | |
Konda et al. | Multi-aperture Fourier ptychographic microscopy, theory and validation | |
Konda et al. | Parallelized aperture synthesis using multi-aperture Fourier ptychographic microscopy | |
CN112506019B (zh) | 基于克罗内克积插值的离轴数字全息成像重建方法 | |
Li et al. | Deep learning assisted variational Hilbert quantitative phase imaging | |
CN109724511A (zh) | 一种数字全息显微的二次相位误差补偿方法及其装置 | |
CN109116709A (zh) | 一种弱离轴相移数字全息相移值抽取与物光恢复方法 | |
CN113393472A (zh) | 基于数字全息显微成像的划区复合相位解包裹方法 | |
WO2021003380A1 (en) | Calibration-free phase shifting procedure for self-interference holography | |
CN107290058A (zh) | 基于偏振正交全息的琼斯矩阵参量同步测量装置与方法 | |
CN107589542B (zh) | 宽波段相位差图像重建中的中心波长的选择方法 | |
Fan et al. | Accurate dynamic quantitative phase imaging using multi-wavelength multiplexing | |
CN110017776B (zh) | 基于移位和多项式拟合的全息像差绝对校准方法及系统 | |
US20210191099A1 (en) | System for quantitative differential phase contrast microscopy with isotropic transfer function | |
Chen et al. | Time-lapse imaging using dual-color coded quantitative differential phase contrast microscopy |
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 |