CN105310676A - 基于光学相干原理的成像技术中的位相恢复的方法 - Google Patents

基于光学相干原理的成像技术中的位相恢复的方法 Download PDF

Info

Publication number
CN105310676A
CN105310676A CN201510795750.7A CN201510795750A CN105310676A CN 105310676 A CN105310676 A CN 105310676A CN 201510795750 A CN201510795750 A CN 201510795750A CN 105310676 A CN105310676 A CN 105310676A
Authority
CN
China
Prior art keywords
phase
optical coherence
vertical direction
doppler
imaging
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
Application number
CN201510795750.7A
Other languages
English (en)
Other versions
CN105310676B (zh
Inventor
王益民
姚晓天
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hebei University
Original Assignee
Hebei University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hebei University filed Critical Hebei University
Priority to CN201510795750.7A priority Critical patent/CN105310676B/zh
Publication of CN105310676A publication Critical patent/CN105310676A/zh
Application granted granted Critical
Publication of CN105310676B publication Critical patent/CN105310676B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于光学相干原理的成像技术中的位相恢复的方法，所述方法对位相折叠引起的信号误差进行消除，其特征在于：所述方法包括获取样品内部的结构和位相；利用得到的振幅和相位值，将每个图像点所对应的复数值的实部和虚部计算出来；在竖直方向或水平方向上，设定边界处的相位值在竖直方向或水平方向上，利用获得的复数值，将相邻两点之间的相位差计算出来；在竖直方向或水平方向上，利用设定的初始位相以及计算出来的位相差利用迭代计算的方法，将实际的位相信号恢复出来。

Description

基于光学相干原理的成像技术中的位相恢复的方法
技术领域
本发明涉及一种基于光学相干原理的成像技术中的位相进行恢复的方法。
背景技术
光学相干层析成像技术(OCT)是一种非接触，无损伤的医用光学成像技术。通过对样品内部散射光的检测，OCT可以对生物样品内部的结构进行高分辨的成像，分辨率在微米量级。除了能够获得生物组织的结构信息，傅里叶域的多普勒OCT(FD-OCT)技术可以测量组织中的血液流动[1,2],这项技术对于人体内部以及生物组织中的血液循环系统的研究具有重要意义。
在多普勒FD-OCT中[3,4],被运动的血流散射出来的探测光会产生一个多普勒频率移动Δν,它正比于血流的流速V.这个频率移动会对探测到的FD-OCT光谱信号产生一个相移。经过快速傅里叶变换，变换结果是一个复函数F(z)，它包含振幅A(z)和位相Φ(z).利用公式(1)可以计算相邻两条竖直采样线之间的位相差,
从而得到多普勒频率移动其中T是相邻两条竖直采样线之间的时间间隔，N是用于平均的竖直采样线的条数。很遗憾，基于位相的多普勒FD-OCT成像技术有一个缺陷。测量到的位相值在数学上被限定在-π到+π的范围，这是反正切函数的主值。因此，多普勒FD-OCT的位相包含折叠和非连续部分。在多普勒FD-OCT中,位相折叠问题将最大可探测多普勒频率移动限定在Δνmax＝1/(2T)。
在光学干涉技术中，人们发展了不同的方法对折叠的位相进行恢复。大部分位相恢复方法是通过复杂算法在折叠的位相信息中寻找位相跳变，然后对相关的数据点加上一个2π整数倍的数值[6-8].这些方法很复杂，而且计算量也比较大。在FD-OCT中，通过将探测到的光谱进行分割，人们提出了一种合成波长法对折叠位相进行恢复[9]。但是，如果合成波长不够长，这个方法仍然会产生位相折叠。另外，分割光谱会降低FD-OCT的成像分辨率，这是人们不原意看到的。
在多普勒FD-OCT中,有效的位相恢复对于自动数据处理和实时显示至关重要。然而，到现在为止，在多普勒OCT领域，还没有关于位相折叠问题的深入研究报导。
发明内容
针对上述问题，本发明提供一种基于光学相干原理的成像的位相进行恢复方法。
为达到上述目的，本发明基于光学相干原理的成像的位相进行恢复方法，所述方法对位相折叠引起的信号误差进行消除，所述方法包括
获取样品内部的结构和位相；
利用得到的振幅和相位值，将每个图像点所对应的复数值的实部和虚部计算出来；
在竖直方向或水平方向上，设定边界处的相位值
在竖直方向或水平方向上，利用获得的复数值，将相邻两点之间的相位差计算出来；
在竖直方向或水平方向上，利用设定的初始位相以及计算出来的位相差利用迭代计算的方法，将实际的位相信号恢复出来。
进一步地，利用多普勒OCT的复数值进行，包括实部和虚部。
进一步地，设定边界条件处的初始相位值。
进一步地，相邻两点之间的相位差是利用多普勒信号的复数值的实部和虚部，利用反正切函数运算得到。
进一步地，相位恢复是在多普勒OCT图像的竖直方向进行的。
进一步地，在多普勒OCT图像的竖直方向上，相邻两点的物理距离应该足够小，以使得它们之间的相位差满足：
进一步地，光学相干原理的成像包括多普勒OCT成像、偏振光敏感的OCT成像、相位敏感的光学显微成像。
本发明基于光学相干原理的成像的位相进行恢复方法与现有技术具备如下有益效果：
本发明利用多普勒FD-OCT信号的复数值，通过计算竖直方向上相邻两点之间的位相差，真实的位相信息可以通过迭代计算恢复出来。这个算法也可应用于其他基于位相的光学干涉检测技术中。
附图说明
图1是多普勒FD-OCT位相示意图；
图2是在图1中位置P1Pn处各点的相位值示意图；
图3是在复数坐标系SrSj中显示的多普勒FD-OCT的复数值分布，数据点为图1中所示血管上P1和Pc之间的点；
图4是在多普勒FD-OCT位相图中的某条竖线上的相邻两点Pm和Pm-1；
图5是在复数坐标系中位相(对于点Pm)和(对于点Pm-1)的关系。
具体实施方式
下面结合附图对本发明做进一步的描述。
实施例1
本实施例基于光学相干原理的成像的位相进行恢复方法，以用于对多普勒光学相干层析成像(OCT)中的位相进行恢复，对位相折叠引起的信号误差进行消除为例，所述方法包括如下步骤：
(1)用多普勒OCT得到样品内部的结构和位相，其中的位相值与样品内部的机械运动有关，例如血液流动；
(2)利用(1)中得到的多普勒OCT的振幅和相位值，将每个OCT图像点所对应的复数值的实部和虚部计算出来；
(3)在竖直方向上，设定边界处的相位值
(4)在竖直方向上，利用(2)中获得的复数值，将相邻两点之间的相位差计算出来。
(5)在竖直方向上，利用(3)中设定的初始位相以及(4)中计算出来的位相差利用迭代计算的方法，将实际的多普勒OCT位相信号恢复出来。
进一步地，所述的位相恢复的方法需要利用多普勒OCT的复数值进行，包括实部和虚部。
进一步地，需要设定边界条件处的初始相位值。
进一步地，相邻两点之间的相位差是利用多普勒信号的复数值的实部和虚部，利用反正切函数运算得到。
进一步地，相位恢复是在多普勒OCT图像的竖直方向进行的。
进一步地，在多普勒OCT图像的竖直方向上，相邻两点的物理距离应该足够小，以使得它们之间的相位差满足：
本实施例相位恢复方法也可以在多普勒OCT图像的水平方向进行。
实施例
如图1所示，图1显示了一幅多普勒FD-OCT位相示意图，样品组织为人体或动物组织，例如人眼的视网膜。对于在样品组织内部的血管Vr,由虚线P1Pn所标记位置的各点的模拟位相值被显示在图2中，其中，P1是在组织上边缘的某一点，Pn则是在血管Vr下面的某一点。根据图2,血管Vr的位相值为负，当位相接近–π时，它跳变到+π区域并发生位相折叠。
在多普勒FD-OCT中，利用多普勒OCT的振幅A(z)和位相对于每一个多普勒图像点，我们可以利用公式(2)计算出它的复数值f(z)＝(Sr(z),Sj(z))，
其中Sr(z)和Sj(z)分别是复数f(z)的实部和虚部。我们模拟计算了从点P1到血管中心Pc之间各点的复数值，并把它们显示在图3所示的复数坐标系中。P1和Pc分别在第IV和第II象限中，由白色圆圈标记出来。
对于一个多普勒位相点，例如图3中的点Pa，当它的复数值落在第III或第IV象限时,它的相位值φa在0到–π的区间中,利用公式(1)，可以把它正确计算出来。然而，在血管中心区的点，由于血流速度较高，其复数值可能会转动到第二象限，例如图3中的Pb点。在这种情况下，如果直接使用公式(1),我们将会得到位相值φ’,而不是φb。很明显，位相值φ’在物理上是不正确的，位相折叠因此产生。尽管在多普勒位相曲线中有2π的位相跳变，但在图3所示的复数坐标系中，多普勒复数值的旋转并没有显示出非连续性，或者大的跳变。这表明真实的位相信息被隐含在多普勒OCT的复数值中。通过对复数值的分析，有可能将折叠后的位相无损失的恢复出来。
在图4中,对于样品组织，考虑血管Vr处，在一个竖直方向上相邻的两个位相点Pm和Pm-1,它们的复数值被用黑点显示在图5中。从图5中可以看出，点Pm的位相值可以写为其中是点Pm-1的位相，是相邻两点Pm和Pm-1之间的位相差。位相差可以利用两点的复数值进行计算。对于血管，血流速度沿着血管的横截面连续分布，流速的一维分布可以被拟合成抛物线函数。[10]所以，在一个血管的横截面中，相邻两点的多普勒位相差会是一个小量，不足以引起位相折叠。因此，在多普勒位相图的竖直方向，设定组织上边缘P1点的位相值为对于深度方向第mth点的位相值可以利用以下公式进行计算：
其中fm和fm-1分别是多普勒OCT图像点Pm和Pm-1处的复数值。公式(3)显示了一种迭代算法可对多普勒FD-OCT的位相折叠进行恢复。通过设定组织边缘的相位值或者血管外部任何相关区域的边缘处的位相值，通过计算相邻两点间的位相差组织内部深度方向的任意一点的位相值可以通过迭代方法恢复出来。
为了能够成功的使用本专利所述的迭代算法对多普勒FD-OCT的位相进行恢复，两个相邻多普勒数值点之间的物理间隔必须足够小，以使得它们之间的位相差满足条件：在竖直方向上，如果OCT系统的分辨率足够高，这个条件就能获得满足。例如，对于一个半径为50微米的血管，它的中心最大血流速度为20毫米/秒。假设多普勒角度为80度，那么在OCT探测光束的方向上的最大速度分量为3473微米/秒。考虑抛物线形流速分布[10],测量到的血流速度的函数分布为V＝-1.39r2+3473,其中r是以血管中心为原点沿血管半径方向的位值坐标。流速分布的一阶导数为δV/δr＝-2.78r.所以，在竖直方向上，相邻两个像素之间血流速度的最大变化量为δV/δr|max＝139(微米/秒)/微米。这个最大值发生在血管壁处，那里r＝50微米。对于一个线采样时间间隔为36.7微秒的FD-OCT系统，在位相π的限制下，最大可探测流速为4.2毫米/秒。用δV/δr|max去除4.2毫米/秒，我们得到间距值30.2微米。因此，在竖直方向，如果多普勒FD-OCT图像上相邻两点之间的物理间距小于30.2微米，迭代算法(3)就可以安全使用。在现有的商用FD-OCT系统中，系统在竖直方向的分辨率约为5微米，远小于条件值30.2微米，因此，迭代算法的使用条件完全满足。在以上讨论中，迭代计算是在竖直方向进行的，这个方法也可以在水平方向进行。在水平方向使用时，FD-OCT在水平方向的采样步长应该足够小，除了多普勒FD-OCT领域，本专利描述的方法也可以应用于其它光学干涉仪技术中，例如偏振光敏感的OCT，相位敏感的显微成像技术等。在这种情况下，样品中待测的物理特性参数应具有连续分布。
对本发明应当理解的是，以上所述的实施例，对本发明的目的、技术方案和有益效果进行了进一步详细的说明，以上仅为本发明的实施例而已，并不用于限定本发明，凡是在本发明的精神原则之内，所作出的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内，本发明的保护范围应该以权利要求所界定的保护范围为准。

Claims (7)

1.基于光学相干原理的成像技术中的位相恢复的方法，所述方法对位相折叠引起的信号误差进行消除，其特征在于：所述方法包括
获取样品内部的结构和位相；
利用得到的振幅和相位值，将每个图像点所对应的复数值的实部和虚部计算出来；
在竖直方向或水平方向上，设定边界处的相位值
在竖直方向或水平方向上，利用获得的复数值，将相邻两点之间的相位差计算出来；
在竖直方向或水平方向上，利用设定的初始位相以及计算出来的位相差利用迭代计算的方法，将实际的位相信号恢复出来。
2.根据权利要求1所述的基于光学相干原理的成像技术中的位相恢复的方法，其特征在于：利用多普勒OCT的复数值进行，包括实部和虚部。
3.根据权利要求1所述的基于光学相干原理的成像技术中的位相恢复的方法，其特征在于：设定边界条件处的初始相位值。
4.根据权利要求2所述的基于光学相干原理的成像技术中的位相恢复的方法，其特征在于：相邻两点之间的相位差是利用多普勒信号的复数值的实部和虚部，利用反正切函数运算得到。
5.根据权利要求1所述的基于光学相干原理的成像技术中的位相恢复的方法，其特征在于：相位恢复是在多普勒OCT图像的竖直方向进行的。
6.根据权利要求1所述的基于光学相干原理的成像技术中的位相恢复的方法，其特征在于：在多普勒OCT图像的竖直方向上，相邻两点的物理距离应该足够小，以使得它们之间的相位差满足：
7.根据权利要求1所述的基于光学相干原理的成像技术中的位相恢复的方法，其特征在于：光学相干原理的成像包括多普勒OCT成像、偏振光敏感的OCT成像、相位敏感的光学显微成像。
CN201510795750.7A 2015-11-18 2015-11-18 基于光学相干原理的成像技术中的位相恢复的方法 Active CN105310676B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510795750.7A CN105310676B (zh) 2015-11-18 2015-11-18 基于光学相干原理的成像技术中的位相恢复的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510795750.7A CN105310676B (zh) 2015-11-18 2015-11-18 基于光学相干原理的成像技术中的位相恢复的方法

Publications (2)

Publication Number Publication Date
CN105310676A true CN105310676A (zh) 2016-02-10
CN105310676B CN105310676B (zh) 2018-03-16

Family

ID=55239692

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510795750.7A Active CN105310676B (zh) 2015-11-18 2015-11-18 基于光学相干原理的成像技术中的位相恢复的方法

Country Status (1)

Country Link
CN (1) CN105310676B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109223044A (zh) * 2018-07-12 2019-01-18 郑州轻工业学院 一种光学相干层析成像系统的多普勒流速检测方法
CN111436909A (zh) * 2020-04-28 2020-07-24 北京信息科技大学 一种活体组织的光学相干层析成像系统及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6075408A (en) * 1998-12-30 2000-06-13 International Business Machines Corp. OQPSK phase and timing detection
CN101181153A (zh) * 2007-12-12 2008-05-21 中国科学院上海光学精密机械研究所 高分辨率光学相干层析成像方法
CN203016921U (zh) * 2013-01-18 2013-06-26 杭州电子科技大学 一种基于级联马赫曾德干涉仪的光谱相位校准系统
US20140211838A1 (en) * 2013-01-29 2014-07-31 Zte (Usa) Inc. Method and apparatus for algorithm on flexible square-qam coherent detection
CN104935385A (zh) * 2014-03-19 2015-09-23 中兴通讯股份有限公司 用于高阶正交幅度调制信号的盲均衡的技术
CN105030201A (zh) * 2015-07-27 2015-11-11 浙江大学 基于静态区域信息的扫频oct数字相位矫正方法与系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6075408A (en) * 1998-12-30 2000-06-13 International Business Machines Corp. OQPSK phase and timing detection
CN101181153A (zh) * 2007-12-12 2008-05-21 中国科学院上海光学精密机械研究所 高分辨率光学相干层析成像方法
CN203016921U (zh) * 2013-01-18 2013-06-26 杭州电子科技大学 一种基于级联马赫曾德干涉仪的光谱相位校准系统
US20140211838A1 (en) * 2013-01-29 2014-07-31 Zte (Usa) Inc. Method and apparatus for algorithm on flexible square-qam coherent detection
CN104935385A (zh) * 2014-03-19 2015-09-23 中兴通讯股份有限公司 用于高阶正交幅度调制信号的盲均衡的技术
CN105030201A (zh) * 2015-07-27 2015-11-11 浙江大学 基于静态区域信息的扫频oct数字相位矫正方法与系统

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109223044A (zh) * 2018-07-12 2019-01-18 郑州轻工业学院 一种光学相干层析成像系统的多普勒流速检测方法
CN109223044B (zh) * 2018-07-12 2021-08-24 郑州轻工业学院 一种光学相干层析成像系统的多普勒流速检测方法
CN111436909A (zh) * 2020-04-28 2020-07-24 北京信息科技大学 一种活体组织的光学相干层析成像系统及方法
CN111436909B (zh) * 2020-04-28 2023-03-28 北京信息科技大学 一种活体组织的光学相干层析成像系统及方法

Also Published As

Publication number Publication date
CN105310676B (zh) 2018-03-16

Similar Documents

Publication Publication Date Title
CN105342568B (zh) 联合相位和幅值的光学相干造影方法及系统
US8107694B2 (en) Method of tracking speckle displacement between two images
Buracas et al. Contribution of area MT to perception of three-dimensional shape: a computational study
CN106226764A (zh) 一种基于D‑InSAR的煤矿开采地沉陷区域的测定方法
CN107066961B (zh) 指纹配准方法及装置
CN104634460A (zh) 一种分布式fbg传感网络多峰自适应精确寻峰方法
Liu et al. Automated extraction of tidal creeks from airborne laser altimetry data
CN105974430A (zh) 基于随机介质表面散射光的运动目标跟踪系统及跟踪方法
CN106097317A (zh) 一种基于离散余弦相位信息的多光斑检测和定位方法
CN107689051A (zh) 一种基于变化因子的多时相sar影像变化检测方法
CN105310676A (zh) 基于光学相干原理的成像技术中的位相恢复的方法
CN103630885A (zh) 合成孔径雷达的目标识别方法和系统
CN104836990B (zh) 一种桥墩防撞图像监控系统及监控方法
Miao et al. A novel marker for robust and accurate phase-based 2D motion estimation from noisy image data
Warner et al. Pine Island Glacier (Antarctica) velocities from Landsat7 images between 2001 and 2011: FFT-based image correlation for images with data gaps
Miyachi et al. Automated detection of arterial wall boundaries based on correlation between adjacent receive scan lines for elasticity imaging
CN104463863B (zh) 基于时间外差投影的运动干涉场的标定方法与系统
Marron et al. Extended-range digital holographic imaging
CN111175724A (zh) 一种基于质心计算的目标微多普勒分量提取方法
Lawton Outputs of paired Gabor filters summed across the background frame of reference predict the direction of movement (vision)
CN102288117B (zh) 基于光学联合变换相关器的像移探测成像系统
Demi New approach to automatic contour detection from image sequences: An application to ventriculographic images
CN202133375U (zh) 基于光学联合变换相关器的像移探测成像系统
Krotkov Visual hyperacuity: Representation and computation of high precision position information
Cordero-Grande et al. Improving harmonic phase imaging by the windowed Fourier transform

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wang Yimin

Inventor after: Yao Xiaotian

Inventor after: Su Ya

Inventor before: Wang Yimin

Inventor before: Yao Xiaotian

GR01 Patent grant
GR01 Patent grant