CN102805613A - 一种基于两次扫描的高分辨率光学扫描全息切片成像方法 - Google Patents
一种基于两次扫描的高分辨率光学扫描全息切片成像方法 Download PDFInfo
- Publication number
- CN102805613A CN102805613A CN2012102860959A CN201210286095A CN102805613A CN 102805613 A CN102805613 A CN 102805613A CN 2012102860959 A CN2012102860959 A CN 2012102860959A CN 201210286095 A CN201210286095 A CN 201210286095A CN 102805613 A CN102805613 A CN 102805613A
- Authority
- CN
- China
- Prior art keywords
- matrix equation
- slice
- matrix
- under test
- object under
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 41
- 230000003287 optical effect Effects 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 54
- 238000000034 method Methods 0.000 claims description 23
- 238000001093 holography Methods 0.000 claims description 13
- 230000010287 polarization Effects 0.000 claims description 10
- 210000001747 pupil Anatomy 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 238000006116 polymerization reaction Methods 0.000 claims description 4
- 230000007547 defect Effects 0.000 abstract 1
- 230000015572 biosynthetic process Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
Images
Landscapes
- Holo Graphy (AREA)
Abstract
本发明公开了一种基于两次扫描的高分辨率光学扫描全息切片成像方法,属于光学扫描领域,主要解决了现有技术中对任意二维切片的图像重构时存在较大离焦噪声的缺陷。本发明在二维扫描镜上对该物体进行第一次二维扫描,得到包含切片信息的第一矩阵方程后将该物体向二维扫描镜的方向移动距离,对该物体进行第二次扫描,得到包含切片信息的第二矩阵方程;然后将第一矩阵方程和第二矩阵方程整合为一个最小化线性方程,将该线性问题的求解转化为最小化问题,通过引入共轭梯度算法,实现切片成像。通过上述方案,本发明实现了高精度的切片成像,并且大大减少了离焦噪声,适用于各个领域。
Description
技术领域
本发明属于光学扫描领域,具体的说,涉及一种基于两次扫描的高分辨率光学扫描全息切片成像方法。
背景技术
光学扫描全息技术,简称OSH,是一种基于菲涅尔波带板扫描的非传统成像技术,即通过二维光学扫描实现对目标的高分辨率三维成像,它在生物医学成像、荧光物体成像、三维全息电视系统以及光学遥感等领域都有广泛的应用前景。
而通过光学扫描全息技术获得的二维全息图,包含了物体完整的三维信息,因此在光学扫描全息技术中对物体全息图的一个重要分析处理步骤就是物体的切片成像,即物体任意二维切面的图像重构。而物体任意二维切面图像重构中的难点在于如何消除来自物体其他层面的噪声,即离焦噪声。切片成像是一个典型的图像处理中的逆问题,同时也是一个不适定问题。
文献《Optical Scanning Holography with MATLAB》提出了一种传统的切片成像方法,即用物体的全息图与待重构切片处的菲涅尔波带板共轭进行卷积运算,从而实现切片成像,但由于无法消除隔离噪声,因此其应用受到极大的制约。
文献《Three-dimensional microscopy and sectional image reconstruction using optical scanning holography》介绍了一种逆成像算法,此迭代算法能够实现轴向分辨率为1毫米左右的切片成像,并能有效抑制离焦噪声,但其无法在更小的轴向尺寸下实现良好成像。
文献《Depth resolution enhancement in optical scanning holography with a dual-wavelength laser source》提出了一种利用双波长激光器提高切片成像轴向分辨率的方法,其利用输出波长分别为632nm和543nm的激光器,获取两组物体全息图,进而将轴向分辨率提高至2.5微米左右,但由于在光学系统中同时工作的两个不同波长引入了较大的噪声,导致其实用性受到极大限制。
发明内容
本发明的目的在于降低切片成像中的离焦噪声,提出一种基于两次扫描的高分辨率光学扫描全息切片成像方法,通过两次扫描,为切片成像这一不适定逆问题引入更多的线性方程组,从而实现高分辨率的切片成像。
本发明采用的技术方案如下:
一种基于两次扫描的高分辨率光学扫描全息切片成像方法,包括以下步骤:
(1)第一偏振分束器将同一光源发出的光分为两束,将该两束光分别经过处理后再通过第二偏振分束器聚光,聚合后的光在待测物体上产生干涉形成菲涅尔波带板;
其中,所述步骤(1)中形成菲涅尔波带板的具体步骤如下:
(1a)通过第一偏振分束器将光分成两束;
(1b)一束光通过第一光瞳形成平面波,另一束光通过第二光瞳形成球面波;
(1c)平面波与球面波通过第二偏振分束器聚合,在待测物体上产生干涉形成时变的菲涅尔波带板,该菲涅尔波带板的值为:
其中x,y,z为该物体的空间坐标,k为光的波数。
为了得到第一矩阵方程,所述步骤(2)中待测物体为两个离散切片的集合,两个切片的轴向位置分别为z 1 和z 2 ,因此得到第一矩阵方程的具体实现方式如下:
(2a)将该待测物体进行第一次二维扫描,并得到二维全息图:
(2b)将菲涅尔波带板在z 1 和z 2 处的值分别转换为矩阵H 1 (z 1 )和H 1 (z 2 );
(2c)将二维全息图与矩阵H 1 (z 1 )和H 1 (z 2 )结合起来得到第一矩阵方程:
其中n 1 为高斯白噪声,该高斯白噪声是长度为N 2 的一维矢量矩阵。
为了得到第二矩阵方程,所述步骤(3)中得到第二矩阵方程的具体实现方法如下:
(3b)按照第一次二维扫描得到第一矩阵方程的方法得到第二矩阵方程:
其中n2为高斯白噪声,该高斯白噪声是长度为N2的一维矢量矩阵。
(4a)将第一矩阵方程和第二矩阵方程整合,得到:
(4b)将所述矩阵方程(7)转化为最小化线性方程:
(8)
其中H + 为H的共轭转置;
与现有技术相比,本发明具有以下有益效果:
(1)本发明使用了光学扫描全息技术,即通过菲涅尔波带板扫描实现对待测物体的高分辨三维成像,菲涅尔波带板应用范围包括生物、机械、光学、电学类等,因此本发明适用于各个领域,应用范围非常广;
(2)本发明中被扫描的物体在轴向位置移动了,并在此位置上进行了第二次扫描,与只扫描一次相比,多获得了一个二维全息图,因此为切片成像这一不适定逆问题引入更多的线性方程组,使得在求解切片信息时其得到的结果进一步逼近切片信息的真实值,这是从根本上提高切片成像轴向分辨率的原因;
(6)本发明不仅实现方式简单、便于操作,同时具有很强的实用性,适合推广使用。
附图说明
图1为本发明基本结构图。
图2为本发明—实施例的待测物体示意图。
图3为本发明—实施例中不同轴向位置处的菲涅尔波带板图。
图4为本发明—实施例中两次扫描分别获得的正弦全息图和余弦全息图。
图5为本发明—实施例中采用不同的方法在z 1 和z 2 处获得的切片成像图。
具体实施方式
下面结合附图与实施例对本发明作进一步说明,本发明的实施方式包括但不限于下列实施例。
实施例
如图2所示,为了简化问题,图2为本发明中待测物体示意图,该待测物体仅包含两个切片信息。在第一次光学扫描中,待测物体两个切片的轴向位置分别为,,且每个切片尺寸为,矩阵尺寸为,其中按照如下所示步骤即可实现扫描过程:
步骤1 对待测物体进行第一次扫描
(1)如图1所示,由同一光源发出的角频率为ω的光被第一偏振分束器BS1分成两束,其中一束通过光瞳P1(x,y)形成平面波;另一束经过声光调制器产生Ω的频移后再通过光瞳P2(x,y)形成球面波;其中使用的单波长光源中心波长为632nm。
(2)该两束被处理的光经第二偏振分束器BS2合在一起,在待测物体上产生干涉形成如图3所示的时变菲涅尔波带板TD-FZP,然后利用该二维扫描镜2D Scanning控制TD-FZP的偏转,从而实现对三维物体的二维扫描。
当物体被扫描时,光电探测器detector接收传播到探测面的光波,并产生外差电流输出,该外差电流为光生电流,所述光生电流经过混频、放大等电学处理,产生解调信息并存储于计算机中。其中储存的借条信息为同轴FZP编码图像,本质上为包含了物体三维信息的全息图。
而该时变菲涅尔波带板可以表示为:
其中x, y, z代表空间坐标,k为光的波数。从上式看出,对于某一轴向位置z,菲涅尔波带板是一个关于x, y的二维函数;而不同轴向位置z上的菲涅尔波带板是不同的。
其中*代表二维卷积,如果将待测物体看作一系列离散切片的集合,即可对轴向坐标z进行离散化处理,表示为z 1, z 2,..., z n.,分别代表不同切片所在的轴向位置。那么(10)式表征的二维全息图可以表示为:
由于只有两个切片,那么上式可以简化为如下形式:
为了将(2)式简化为一组线性方程进行分析,我们可以分别将和分别转换为一维矢量矩阵ψ 1 和ψ 2 。如果待测物体的切片为一个N×N的矩阵,ψ 1 和ψ 2 则为长度为N 2 的一维矢量矩阵。同样,物体的二维全息图g 1 (x,y)也可以转化为一个长度为N 2 一维矢量矩阵G 1 。
上式的卷积运算可以表示为矩阵运算,为此用菲涅尔波带板在z 1 和z 2 处分别为N×N的矩阵h(x,y;z 1 )和N×N的矩阵h(x,y;z 2 ),分别构造如下所示的两个N 2 ×N 2 的二维矩阵H 1 (z 1 )和H 1 (z 2 ):
则通过(12)式和(13)式可以得到第一矩阵方程:
其中n 1 和n 2 代表系统的高斯白噪声,是长度为N 2 的一维矢量矩阵。
所谓切片成像,即要从G 1 中恢复出切片信息的信息,这是一个不适定逆问题,因此我们对物体进行第二次扫描,以获得第二组物体全息图。
步骤2 对物体进行第二次扫描
从(6)式可以看出,通过第二次扫描我们又获得了N 2 个线性方程组。由于第二次扫描是在轴向位置z 1 ’和z 2 ’完成的,因此菲涅尔波带板是不一样的,即H 2 ≠H 1 ,可见两次扫描获得的线性方程是不同的。这样,我们就为该不适定问题的求解添加了更多的有效线性方程组。
步骤3 根据所有的线性方程组进行切片成像
首先将两次二维扫描的矩形方程组合起来,表示为:
该问题的求解可转化为如下的最小化问题,即:
其中(8)式中|| ||表示二阶范数,λ>0为罚系数,C是拉普拉斯算子。该最小化问题的解可表示为:
如图5所示,其中图5(a)—(h)分别展示了待测物体的切片实例以及用三种不同方法进行切片成像的结果。由图5(c)—(d)可以看出,传统切片成像方法不仅不能区分两个切片,而且引入了较大的离焦噪声;而图5(e)—(f)展示的基于单次扫描的逆成像方法虽然能够在一定程度上抑制离焦噪声,但是无法达到较高的即大于100nm的轴向分辨率;而采用本发明的两次扫描方法,结合共轭梯度算法,则可以完全区分出轴向距离为100nm的两个深度切片,其结果如图5(g)—(h)所示,这意味着,本发明方法将光学扫描全息术的轴向分辨率提高到了100nm。
按照上述实施例,便可很好的实现本发明。
Claims (6)
3.根据权利要求1所述的一种基于两次扫描的高分辨率光学扫描全息切片成像方法,其特征在于,所述步骤(2)中待测物体为两个离散切片的集合,两个切片的轴向位置分别为z 1 和z 2 ,因此得到第一矩阵方程的具体实现方式如下:
(2a)将该待测物体进行第一次二维扫描,并得到二维全息图:
(2b)将菲涅尔波带板在z 1 和z 2 处的值分别转换为矩阵H 1 (z 1 )和H 1 (z 2 );
(2c)将二维全息图与矩阵H 1 (z 1 )和H 1 (z 2 )结合起来得到第一矩阵方程:
(3)
其中n 1 为高斯白噪声,该高斯白噪声是长度为N 2 的一维矢量矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210286095.9A CN102805613B (zh) | 2012-08-13 | 2012-08-13 | 一种基于两次扫描的高分辨率光学扫描全息切片成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210286095.9A CN102805613B (zh) | 2012-08-13 | 2012-08-13 | 一种基于两次扫描的高分辨率光学扫描全息切片成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102805613A true CN102805613A (zh) | 2012-12-05 |
CN102805613B CN102805613B (zh) | 2014-05-28 |
Family
ID=47229612
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210286095.9A Expired - Fee Related CN102805613B (zh) | 2012-08-13 | 2012-08-13 | 一种基于两次扫描的高分辨率光学扫描全息切片成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102805613B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104159094A (zh) * | 2014-07-09 | 2014-11-19 | 四川大学 | 一种改善光学扫描全息层析成像效果的方法 |
CN104616267A (zh) * | 2015-02-17 | 2015-05-13 | 海南大学 | 基于二阶总广义变差最小化的欠采样图像重构的方法 |
CN105204311A (zh) * | 2015-07-06 | 2015-12-30 | 电子科技大学 | 一种基于高斯切趾的光学扫描全息边缘检测方法 |
CN107015466A (zh) * | 2017-04-19 | 2017-08-04 | 电子科技大学 | 基于tr‑music算法的光学扫描全息单点定位方法 |
CN107145053A (zh) * | 2017-05-22 | 2017-09-08 | 电子科技大学 | 基于tr‑music算法的光学扫描全息轴向定位方法 |
CN107193002A (zh) * | 2017-05-22 | 2017-09-22 | 电子科技大学 | 一种可抑制宽带相位噪声的一维距离像高分辨成像方法 |
CN107240074A (zh) * | 2017-05-15 | 2017-10-10 | 电子科技大学 | 基于二维最佳直方图熵法与遗传算法的离焦躁声消除方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04212191A (ja) * | 1990-04-06 | 1992-08-03 | Hamamatsu Photonics Kk | 光ヘテロダイン走査型ホログラフィ装置 |
US20060078113A1 (en) * | 2001-03-30 | 2006-04-13 | Bahram Javidi | Information security using digital holography |
CN102499648A (zh) * | 2011-11-16 | 2012-06-20 | 清华大学 | 一种基于菲涅尔光谱仪的频域光学相干层析成像系统 |
-
2012
- 2012-08-13 CN CN201210286095.9A patent/CN102805613B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04212191A (ja) * | 1990-04-06 | 1992-08-03 | Hamamatsu Photonics Kk | 光ヘテロダイン走査型ホログラフィ装置 |
US20060078113A1 (en) * | 2001-03-30 | 2006-04-13 | Bahram Javidi | Information security using digital holography |
CN102499648A (zh) * | 2011-11-16 | 2012-06-20 | 清华大学 | 一种基于菲涅尔光谱仪的频域光学相干层析成像系统 |
Non-Patent Citations (2)
Title |
---|
EDMUND Y. LAM,ECT: "Three-dimensional microscopy and sectional image reconstruction using optical scanning holography", 《APPLIED OPTICS》 * |
平子良,等: "用非线性晶体的相位共轭特性记录扫描全息图", 《中国激光》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104159094A (zh) * | 2014-07-09 | 2014-11-19 | 四川大学 | 一种改善光学扫描全息层析成像效果的方法 |
CN104616267A (zh) * | 2015-02-17 | 2015-05-13 | 海南大学 | 基于二阶总广义变差最小化的欠采样图像重构的方法 |
CN105204311A (zh) * | 2015-07-06 | 2015-12-30 | 电子科技大学 | 一种基于高斯切趾的光学扫描全息边缘检测方法 |
CN105204311B (zh) * | 2015-07-06 | 2018-05-18 | 电子科技大学 | 一种基于高斯切趾的光学扫描全息边缘检测方法 |
CN107015466A (zh) * | 2017-04-19 | 2017-08-04 | 电子科技大学 | 基于tr‑music算法的光学扫描全息单点定位方法 |
CN107015466B (zh) * | 2017-04-19 | 2019-07-19 | 电子科技大学 | 基于tr-music算法的光学扫描全息单点定位方法 |
CN107240074A (zh) * | 2017-05-15 | 2017-10-10 | 电子科技大学 | 基于二维最佳直方图熵法与遗传算法的离焦躁声消除方法 |
CN107145053A (zh) * | 2017-05-22 | 2017-09-08 | 电子科技大学 | 基于tr‑music算法的光学扫描全息轴向定位方法 |
CN107193002A (zh) * | 2017-05-22 | 2017-09-22 | 电子科技大学 | 一种可抑制宽带相位噪声的一维距离像高分辨成像方法 |
CN107193002B (zh) * | 2017-05-22 | 2019-04-26 | 电子科技大学 | 一种可抑制宽带相位噪声的一维距离像高分辨成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102805613B (zh) | 2014-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102920438B (zh) | 一种基于可变光瞳的高分辨率光学扫描全息切片成像方法 | |
CN102805613B (zh) | 一种基于两次扫描的高分辨率光学扫描全息切片成像方法 | |
US10082662B2 (en) | 3D refractive index tomography and structured illumination microscopy system using wavefront shaper and method thereof | |
US9892812B2 (en) | Fourier ptychographic x-ray imaging systems, devices, and methods | |
CN104459971B (zh) | 一种基于集成光子芯片的结构光照明显微成像系统 | |
CN106885796B (zh) | 一种超分辨率荧光数字全息断层显微成像系统与方法 | |
CN105184295B (zh) | 一种基于小波变换与连通域的全息扫描空间距离提取方法 | |
CN108089425A (zh) | 一种基于深度学习的消除光学扫描全息离焦噪声的方法 | |
CN110160440A (zh) | 一种基于频域oct技术的三维彩色动态成像装置及方法 | |
KR102082747B1 (ko) | 초점거리 조절이 가능한 led 어레이 기반 3차원 이미징 장치 및 방법 | |
CN104503070A (zh) | 一种基于压缩感知的宽视场显微成像系统及方法 | |
CN108982455A (zh) | 一种多焦点光切片荧光显微成像方法和装置 | |
CN107967703B (zh) | 一种基于连通域的光学扫描全息的自聚焦和重建方法 | |
US9696255B2 (en) | Image processing method of two-photon structured illumination point scanning microscopy | |
Liu et al. | Multiview three-dimensional imaging using a Risley-prism-based spatially adaptive virtual camera field | |
Wei et al. | Efficient shape reconstruction of microlens using optical microscopy | |
CN107835074B (zh) | 一种消除随机加密光学扫描全息离焦噪声的方法 | |
CN106768342B (zh) | 基于偏振复用实现不等间隔多平面成像的装置及方法 | |
CN113009801B (zh) | 高速多方向线共焦数字全息三维显微成像方法及装置 | |
CN108153132B (zh) | 一种基于均值梯度函数的光学扫描全息自聚焦方法 | |
CN110383044A (zh) | 细胞观察装置 | |
Chu et al. | High-resolution measurement based on the combination of multi-vision system and synthetic aperture imaging | |
Li et al. | Single shot 3D imaging by active illumination and point-spread-function engineering | |
CN108594413A (zh) | 一种基于双光束频率调制的三维高分辨率成像方法与装置 | |
Kostencka et al. | Off-axis illumination in object-rotation diffraction tomography for enhanced alignment and resolution |
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 |
Granted publication date: 20140528 Termination date: 20150813 |
|
EXPY | Termination of patent right or utility model |