CN109297595B - 一种光学相干层析相位解卷绕的方法及装置 - Google Patents

一种光学相干层析相位解卷绕的方法及装置 Download PDF

Info

Publication number
CN109297595B
CN109297595B CN201811024739.0A CN201811024739A CN109297595B CN 109297595 B CN109297595 B CN 109297595B CN 201811024739 A CN201811024739 A CN 201811024739A CN 109297595 B CN109297595 B CN 109297595B
Authority
CN
China
Prior art keywords
phase
light
spectrum
lens
theta
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.)
Active
Application number
CN201811024739.0A
Other languages
English (en)
Other versions
CN109297595A (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.)
Northeastern University Qinhuangdao Branch
Original Assignee
Northeastern University Qinhuangdao Branch
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 Northeastern University Qinhuangdao Branch filed Critical Northeastern University Qinhuangdao Branch
Priority to CN201811024739.0A priority Critical patent/CN109297595B/zh
Publication of CN109297595A publication Critical patent/CN109297595A/zh
Application granted granted Critical
Publication of CN109297595B publication Critical patent/CN109297595B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/12Generating the spectrum; Monochromators
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0062Arrangements for scanning
    • A61B5/0066Optical coherence imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/12Generating the spectrum; Monochromators
    • G01J3/26Generating the spectrum; Monochromators using multiple reflection, e.g. Fabry-Perot interferometer, variable interference filters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/12Generating the spectrum; Monochromators
    • G01J2003/1204Grating and filter

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种光学相干层析相位解卷绕的装置,包括,宽带光源,用于发出光线;耦合器,用于将光线进行耦合;第一透镜,用于将发散光线变为准直光线;第二透镜,用于将发散光线变为准直光线;第三透镜,用于将准直光线进行汇集;光学玻璃,用于实现部分光线的透射和部分光线的反射;样品台,用于盛放样品;透射光栅,用于产生衍射光谱;第四透镜,用于将光线进行汇集;CCD线阵相机,用于采集光谱信息,并将光谱信息传输到处理器;处理器,用于对光谱信息进行分析。本发明能够改进现有技术的不足,解决SDOCT中相位卷绕的问题。

Description

一种光学相干层析相位解卷绕的方法及装置
技术领域
本发明涉及三维层析成像技术领域,特别是一种光学相干层析相位解卷绕的方法及装置。
背景技术
在光学相干层析(optical coherence tomography,OCT)中,相位信息是其中重要的组成部分,利用OCT的相位信息可以实现纳米级材料表面轮廓测量、细胞相位成像、血液流速检测等。近年来,发展了多种相位测量技术,如衍射相位显微镜、傅里叶相位显微镜、希尔伯特相位显微镜以及数字全息技术等,能够实现纳米级别的高精度相位成像。然而相位卷绕是限制上述相位成像技术应用的主要问题,相位卷绕是干涉检测的固有问题,其描述的是当相位变化超过一个周期2π之后又会重复变化,从而产生相位多值性。长期以来,有许多数值去卷绕方法被提出,其基本原理都是根据相位图的连续性对原卷绕位图进行补偿。但是当相位噪声较大,或者相位急剧变化造成过度卷绕时,这些数值去卷绕方法就会受到影响而失效,尤其是当相位急剧变化超过π时,理论上数值去卷绕方法已无法恢复真实相位。
近来,一种基于低相干光干涉的谱域光学相干层析(Spectral Domain OpticalCoherence Tomography,SDOCT)提供了简便的相位计算方法,该技术在提高相位成像的稳定性同时,实现了亚纳米级的精度。但和其他相位成像方法类似,SDOCT也存在相位卷绕,例如相位成像时,相位物体的尖锐变化区域和多普勒OCT测血流时流速湍急处,都会由于相位的过度卷绕而影响结果。目前已有数种方法被提出以解决SD-OCT中的相位卷绕问题。王毅等(基于谱域相位分辨光学相干层析的纳米级表面形貌成像,物理学报,2017)利用相邻两点复共轭相减,将相邻两点相位差绝对值的限制条件由目前的π扩大到2π;Hansford等(Synthetic wavelength based phase unwrapping in spectral domain opticalcoherence tomography,OPTICS EXPRESS,2009)利用合成波长方法拓展了测量范围;Jun等(High-dynamic-range quantitative phase imaging with spectral domain phasemicroscopy,OPTICS LETTERS,2009)由干涉光谱相位函数的斜率计算绝对相位。但是,以上方法在拓展了测量范围的同时也成倍放大了噪声,故得到的相位只能作为参照,通过2π分段对卷绕相位进行补偿,如Jun报道的“High-dynamic-range quantitative phaseimaging with spectral domain phase microscopy” (OPTICS LETTERS,2009),Yangzhi等报道的“High-sensitive and broad-dynamic-range quantitative phase imagingwith spectral domain phase microscopy”(OPTICS EXPRESS,2013)。这种分段补偿方法的缺陷在于在每一段分段的边界处,放大的噪声会导致卷绕周期的误判,产生2π的补偿误差,从而使得最终噪声被放大到2π,该边界错误即使在较小的噪声条件下依然无法避免。
发明内容
本发明要解决的技术问题是提供一种光学相干层析相位解卷绕的方法及装置,解决SDOCT中相位卷绕的问题。
为解决上述技术问题,本发明所采取的技术方案如下。
一种光学相干层析相位解卷绕的装置,包括,
宽带光源,用于发出光线;
耦合器,用于分光;
第一透镜,用于将发散光线变为准直光线;
第二透镜,用于将发散光线变为准直光线;
第三透镜,用于将准直光线进行汇集;
光学玻璃,用于实现部分光线的透射和部分光线的反射;由其下表面作为参考面,其下表面反射的光为参考光,其透射光为探测光;
样品台,用于盛放样品;
透射光栅,用于产生衍射光谱;
第四透镜,用于将经过透射光栅衍射的光成像于CCD线阵相机;
CCD线阵相机,用于采集光谱信息,并将光谱信息传输到计算机;
计算机,用于对光谱信息进行分析。
一种上述的光学相干层析相位解卷绕的装置的相位解卷绕的方法,包括以下步骤,
A、宽带光源发出的光经过耦合器后射出,由第一透镜准直,由第三透镜汇聚后经过光学玻璃,聚焦在搭载在样品台的样品上,经样品反射的样品光沿原路返回,经过耦合器、光栅,进入CCD线阵相机;
B、由光学玻璃的下表面反射的光,沿原路返回,经过耦合器、光栅,进入CCD线阵相机;
C、进入CCD线阵相机的参考光和样品光进行干涉,形成干涉光谱, CCD线阵相机将采集到的干涉光谱传输到计算机进行分析。
作为优选,步骤C中,对光谱信息进行分析包括以下步骤,
C1、获取谱域光学相干层析的干涉光谱P0(k),k为波数,经过滤波,消除直流分量,经傅里叶变换,得到干涉光谱的幅度谱和相位谱,设定幅度谱的极大值点坐标为(x0,a0),x0和a0分辨表示横坐标和纵坐标的值,设定在相位谱中,横坐标x0对应的相位为θ0;由傅里叶变换得到的θ0为卷绕相位,设定实际相位为θ=θ0+2n0π,n0为待确定的卷绕次数;
C2、生成仿真光谱Pn(k),
Figure BDA0001788229020000031
上式中k、θ′、kc为系统参数,k表示光谱仪的光谱轴,用波数表示,θ′为初始相位,kc为k的中心值;
C3、将光谱轴k分成相等的两部分,k1和k2,k1和kz的中心值分别为 k1c和k2c。将测量到的干涉光谱P0(k)和仿真光谱Pn(k)按照k1和k2分成两部分,分别为P0(k1)和P0(k2)、Pn(k1)和Pn(k2);分别如步骤C1 所示方法计算P0(k1)、P0(k2)、Pn(k1)和Pn(k2)的相位,分别为θ01、θ02、θn1和θn2
C4、在公式(1)中n取不同的值,计算Pn(k),经步骤C2,得到不同的θn1和θn2,令
Figure BDA0001788229020000041
为:
Figure BDA0001788229020000042
当n等于nmin
Figure BDA0001788229020000043
为最小值时,此时的nmin即为待确定的卷绕次数n0
得到实际相位为(θ0+2πnmin)。
采用上述技术方案所带来的有益效果在于:本发明提供的方法在进行相位补偿时,无需参考高噪声的非卷绕相位,不存在边界误判问题;也不需要参考临近点相位做卷绕判断,可以在单点采集时直接得到非卷绕的相位信息,同时保留了原相位信息的低噪声水平,不存在噪声放大问题。
附图说明
图1是本发明一个具体实施方式的结构图。
具体实施方式
参照图1,本发明一个具体实施方式包括,
宽带光源1,用于发出光线;
耦合器2,用于分光;
第一透镜3,用于将发散光线变为准直光线;
第二透镜4,用于将发散光线变为准直光线;
第三透镜5,用于将准直光线进行汇集;光学玻璃6,用于实现部分光线的透射和部分光线的反射,由其下表面作为参考面,其下表面反射的光为参考光,其透射光为探测光;
样品台8,用于盛放样品7;
透射光栅9,用于产生衍射光谱;
第四透镜10,用于将经过透射光栅9衍射的光成像于CCD线阵相机11;
CCD线阵相机11,用于采集光谱信息,并将光谱信息传输到计算机12;
计算机12,用于对光谱信息进行分析。
一种上述的光学相干层析相位解卷绕的装置的相位解卷绕的方法,包括以下步骤,
A、宽带光源1发出的光经过耦合器2后射出,由第一透镜3准直,由第三透镜5汇聚后经过光学玻璃6,聚焦在搭载在样品台8的样品7上,经样品反射的样品光沿原路返回,经过耦合器2、光栅9,进入CCD线阵相机11;
B、由光学玻璃6的下表面反射的参考光,沿原路返回,经过耦合器 2、光栅9,进入CCD线阵相机11;
C、进入CCD线阵相机11的参考光和样品光进行干涉,形成干涉光谱,CCD线阵相机11将采集到的干涉光谱传输到计算机12进行分析。
步骤C中,对光谱信息进行分析包括以下步骤,
C1、获取谱域光学相干层析的干涉光谱P0(k),k为波数,经过滤波,消除直流分量,经傅里叶变换,得到干涉光谱的幅度谱和相位谱,设定幅度谱的极大值点坐标为(x0,a0),x0和a0分辨表示横坐标和纵坐标的值,设定在相位谱中,横坐标x0对应的相位为θ0;由傅里叶变换得到的θ0为卷绕相位,设定实际相位为
Figure BDA0001788229020000062
n0为待确定的卷绕次数;
C2、生成仿真光谱Pn(k),
Figure BDA0001788229020000061
上式中k、θ′、kc为系统参数,k表示光谱仪的光谱轴,用波数表示,θ′为初始相位,kc为k的中心值;
C3、将光谱轴k分成相等的两部分,k1和k2,k1和k2的中心值分别为 k1c和k2c;将测量到的干涉光谱P0(k)和仿真光谱Pn(k)按照k1和k2分成两部分,分别为P0(k1)和P0(k2)、Pn(k1)和Pn(k2);分别如步骤C1 所示方法计算P0(k1)、P0(k2)、Pn(k1)和Pn(k2)的相位,分别为θ01、θ02、θn1和θn2
C4、在公式(1)中n取不同的值,计算Pn(k),经步骤C2,得到不同的θn1和θn2,令
Figure BDA0001788229020000071
为:
Figure BDA0001788229020000072
当n等于nmin
Figure BDA0001788229020000073
为最小值时,此时的nmin即为待确定的卷绕次数n0
得到实际相位为(θ0+2πnmin)。
与现有的相位去卷绕方法相比,本实施例具有以下优点:(1)本发明在提取相位信息时,无需做相邻点判断和相位参照,可以在单点采集时直接得到非卷绕相位;(2)不存在边界误判问题,同时保留了原相位信息的低噪声水平,不存在噪声放大问题。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (1)

1.一种应用于光学相干层析相位解卷绕的装置的相位解卷绕的方法,其特征在于:所述光学相干层析相位解卷绕的装置包括,
宽带光源(1),用于发出光线;
耦合器(2),用于分光;
第一透镜(3),用于将发散光线变为准直光线;
第二透镜(4),用于将发散光线变为准直光线;
第三透镜(5),用于将准直光线进行汇集;光学玻璃(6),用于实现部分光线的透射和部分光线的反射,由其下表面作为参考面,其下表面反射的光为参考光,其透射光为探测光;
样品台(8),用于盛放样品(7);
透射光栅(9),用于产生衍射光谱;
第四透镜(10),用于将经过透射光栅(9)衍射的光成像于CCD线阵相机(11);
CCD线阵相机(11),用于采集光谱信息,并将光谱信息传输到计算机(12);
计算机(12),用于对光谱信息进行分析;
该方法包括以下步骤,
A、宽带光源(1)发出的光经过耦合器(2)后射出,由第一透镜(3)准直,由第三透镜(5)汇聚后经过光学玻璃(6),聚焦在搭载在样品台(8)的样品(7)上,经样品反射的样品光沿原路返回,经过耦合器(2)、光栅(9),进入CCD线阵相机(11);
B、由光学玻璃(6)的下表面反射的参考光,沿原路返回,经过耦合器(2)、光栅(9),进入CCD线阵相机(11);
C、进入CCD线阵相机(11)的参考光和样品光进行干涉,形成干涉光谱,CCD线阵相机(11)将采集到的干涉光谱传输到计算机(12)进行分析;
步骤C中,对光谱信息进行分析包括以下步骤,
C1、获取谱域光学相干层析的干涉光谱P0(k),k为波数,经过滤波,消除直流分量,经傅里叶变换,得到干涉光谱的幅度谱和相位谱,设定幅度谱的极大值点坐标为(x0,a0),x0和a0分辨表示横坐标和纵坐标的值,设定在相位谱中,横坐标x0对应的相位为θ0;由傅里叶变换得到的θ0为卷绕相位,设定实际相位为θ=θ0+2n0π,n0为待确定的卷绕次数;
C2、生成仿真光谱Pn(k),
Figure FDA0002677375240000021
上式中k、θ′、kc为系统参数,k表示光谱仪的光谱轴,用波数表示,θ′为初始相位,kc为k的中心值;
C3、将光谱轴k分成相等的两部分,k1和k2,k1和k2的中心值分别为k1c和k2c;将测量到的干涉光谱P0(k)和仿真光谱Pn(k)按照k1和k2分成两部分,分别为P0(k1)和P0(k2)、Pn(k1)和Pn(k2);分别如步骤C1所示方法计算P0(k1)、P0(k2)、Pn(k1)和Pn(k2)的相位,分别为θ01、θ02、θn1和θn2
C4、在公式(1)中n取不同的值,计算Pn(k),经步骤C3,得到不同的θn1和θn2,令
Figure FDA0002677375240000022
为:
Figure FDA0002677375240000023
当n等于nmin,
Figure FDA0002677375240000024
为最小值时,此时的nmin即为待确定的卷绕次数n0
得到实际相位为(θ0+2πnmin)。
CN201811024739.0A 2018-09-04 2018-09-04 一种光学相干层析相位解卷绕的方法及装置 Active CN109297595B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811024739.0A CN109297595B (zh) 2018-09-04 2018-09-04 一种光学相干层析相位解卷绕的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811024739.0A CN109297595B (zh) 2018-09-04 2018-09-04 一种光学相干层析相位解卷绕的方法及装置

Publications (2)

Publication Number Publication Date
CN109297595A CN109297595A (zh) 2019-02-01
CN109297595B true CN109297595B (zh) 2020-11-03

Family

ID=65166197

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811024739.0A Active CN109297595B (zh) 2018-09-04 2018-09-04 一种光学相干层析相位解卷绕的方法及装置

Country Status (1)

Country Link
CN (1) CN109297595B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110261309B (zh) * 2019-07-02 2021-09-28 东北大学秦皇岛分校 一种磁动式光学相干层析成像系统及其磁场调制方法
CN112748089B (zh) * 2019-10-31 2022-02-01 北京理工大学 多普勒光学相干层析成像中的相位解缠绕方法及装置
CN113251945B (zh) * 2021-05-17 2022-07-12 东北大学秦皇岛分校 一种线轮廓成像装置的解调方法及成像装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100150467A1 (en) * 2008-07-21 2010-06-17 Mingtao Zhao Methods, systems, and computer readable media for synthetic wavelength-based phase unwrapping in optical coherence tomography and spectral domain phase microscopy
CN106932363A (zh) * 2017-03-10 2017-07-07 南京市计量监督检测院 一种钻石纯度检测系统及方法
CN107661089A (zh) * 2017-09-19 2018-02-06 北京工业大学 一种频域光学相干层析连续色散补偿成像方法和系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100150467A1 (en) * 2008-07-21 2010-06-17 Mingtao Zhao Methods, systems, and computer readable media for synthetic wavelength-based phase unwrapping in optical coherence tomography and spectral domain phase microscopy
CN106932363A (zh) * 2017-03-10 2017-07-07 南京市计量监督检测院 一种钻石纯度检测系统及方法
CN107661089A (zh) * 2017-09-19 2018-02-06 北京工业大学 一种频域光学相干层析连续色散补偿成像方法和系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《基于谱域相位分辨光学相干层析的纳米级表面形貌成像》;王毅 等;《物理学报》;20171231;第66卷(第15期);第154202-1至154202-10页 *

Also Published As

Publication number Publication date
CN109297595A (zh) 2019-02-01

Similar Documents

Publication Publication Date Title
Vo et al. Surface recovery algorithm in white light interferometry based on combined white light phase shifting and fast Fourier transform algorithms
US8565499B2 (en) Methods, systems, and computer readable media for synthetic wavelength-based phase unwrapping in optical coherence tomography and spectral domain phase microscopy
CN109297595B (zh) 一种光学相干层析相位解卷绕的方法及装置
CN101832817B (zh) 并行复频域光学相干层析成像方法与系统
CN110779464B (zh) 一种时域频域联合分析宽光谱相干测量方法及系统
CN102818786B (zh) 正弦相位调制的并行复频域光学相干层析成像系统和方法
CN110361099B (zh) 一种谱域低相干光干涉光程差解调方法
JP7138734B2 (ja) デュアル干渉法試料厚さ計
CN107228632B (zh) 一种基于加窗傅里叶变换的位移场层析测量装置及方法
CN115046469B (zh) 一种面向光纤白光干涉的干涉条纹包络提取方法
Hendargo et al. Synthetic wavelength based phase unwrapping in spectral domain optical coherence tomography
Gianto et al. Comparison of envelope detection techniques in coherence scanning interferometry
CA2899651A1 (en) N-wavelength interrogation system and method for multiple wavelength interferometers
JP3934490B2 (ja) 低コヒーレント干渉縞解析方法
Kim et al. Fast and reliable measurement of thin film thickness profile based on wavelet transform in spectrally resolved white-light interferometry
CN114894308A (zh) 一种基于低相干干涉的光谱仪标定方法与系统
CN104280139A (zh) 一种动态相位测量装置及方法
Hovell et al. From light to displacement: A design framework for Optimising spectral-domain low-coherence interferometric sensors for in situ measurement
JP2009031229A (ja) 測定信号のノイズ処理方法
JP6274555B2 (ja) 群遅延演算を用いたofdr方式光ファイバ計測方法及びそれを実施する装置
Xu et al. A zero-crossing detection method applied to Doppler OCT
CN109557035B (zh) 一种正弦相位差分解调的多普勒优化成像方法
US7864327B1 (en) In-phase/in-quadrature demodulator for spectral information of interference signal
Liang et al. Pseudo Wigner-Ville distribution for 3D white light scanning interferometric measurement
CN102680958A (zh) 一种适用于多普勒测速仪的测速信号处理方法

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