CN113984713A - 基于单线扫的光学相干层析成像色散失配探测及补偿方法 - Google Patents
基于单线扫的光学相干层析成像色散失配探测及补偿方法 Download PDFInfo
- Publication number
- CN113984713A CN113984713A CN202111226587.4A CN202111226587A CN113984713A CN 113984713 A CN113984713 A CN 113984713A CN 202111226587 A CN202111226587 A CN 202111226587A CN 113984713 A CN113984713 A CN 113984713A
- Authority
- CN
- China
- Prior art keywords
- line
- dispersion
- compensation
- dispersion mismatch
- scan
- 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
- 239000006185 dispersion Substances 0.000 title claims abstract description 64
- 238000012014 optical coherence tomography Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000001514 detection method Methods 0.000 title claims abstract description 20
- 238000001228 spectrum Methods 0.000 claims description 22
- 238000003384 imaging method Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 230000000007 visual effect Effects 0.000 abstract description 2
- 239000011521 glass Substances 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 4
- 239000005357 flat glass Substances 0.000 description 4
- 239000000835 fiber Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000000106 sweat gland Anatomy 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000002612 dispersion medium Substances 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0064—Body surface scanning
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/44—Detecting, measuring or recording for evaluating the integumentary system, e.g. skin, hair or nails
- A61B5/441—Skin evaluation, e.g. for skin disorder diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Psychiatry (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Dermatology (AREA)
- Mathematical Physics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Chemical & Material Sciences (AREA)
- Signal Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明属于光学相干层析成像(Optical coherence tomography,OCT)领域,具体公开了一种基于单线扫的光学相干层析成像色散失配探测及补偿方法,本发明提出的指数通过计算单线扫(Single A‑line)频谱图中亮线的线斜率,实现色散失配探测及补偿。通过使用新的色散失配指数,所提出的单线扫色散失配探测及补偿技术仅需使用单线扫,可以实现快速地、可视化地探测傅里叶域光学相干层析成像(Fourier‑domain optical coherence tomography,FD‑OCT)设置的色散失配状态并执行基于单线扫的、对样本结构无依赖(Sample‑structure‑nondependent,SSNd)的色散失配补偿方法。
Description
技术领域
本发明属于光学相干层析成像技术领域,具体涉及一种基于单线扫的光学相干层析成像色散失配探测及补偿方法。
背景技术
OCT可以实现强大的无创、高分辨率的生物组织和材料三维成像,已广泛应用于生物医学等领域。在新一代的FD-OCT中使用超宽带光源进行高分辨率成像存在的一个困难是色散失配。当样品臂(Sample arm)和参考臂(Reference arm)含有不同长度的光色散介质时,会发生色散失配,进而恶化FD-OCT系统轴向点扩散函数的半峰全宽(Full width athalf maxima,FWHM),导致成像分辨率降低。因此,色散补偿对于实现OCT的高分辨率起着至关重要的作用,尤其是对于高分辨率OCT来说。在过去的几十年里,人们提出并研究了各种色散补偿方法。有基于硬件的方法,如在一个臂中增加弥散补偿材料、使用双光纤拉伸器和使用声光调制器;有基于软件的方法,如分数傅里叶变换法、相项法、MSI法、色散校准法和网格搜索法等。
然而,基于硬件的方法操作不便,适用性差,需要较高的硬件成本;基于软件的方法经常受到散斑和样品结构的影响。同时,这些传统的基于硬件和基于软件的方法无法可视化地指示色散失配状态。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于单线扫的光学相干层析成像色散失配探测及补偿方法。
本发明的目的是通过以下技术方案来实现的:
基于单线扫的光学相干层析成像色散失配探测及补偿方法,包括以下步骤:
S1:探测色散失配并进行可视化:对单线扫原始数据进行谱空间分析,利用加窗傅里叶变换得到频谱图;然后将频谱图进行二值化和骨架化处理;检测出频谱图中的线并对其进行数值拟合;计算出线的斜率,所述斜率表示色散失配指数;
S2:进行色散补偿:向单线扫原始数据中引入了补偿相位,所述补偿相位公式如下:
Φ=-a2(ω-ω0)2-a3(ω-ω0)3;
S3:将进行色散补偿后的频谱重复步骤S1的步骤,当频谱中斜率不为0时,重复步骤S2,当频谱中斜率为0时,步骤终止,通过重复步骤进行色散补偿实现了快速色散失配补偿。
进一步的,所述步骤S1还包括对单线扫原始数据进行谱空间分析的方法:线扫中的一个亮点对应单线扫频谱图中的一条线,所述线的线斜率值表示所述亮点成像深度下的色散失配状态,所以线扫不同点对应的频谱图中不同线用于可视化地显示不同成像深度的色散失配。
本发明的有益效果是:
1)本发明无需获取和处理OCT(Optical coherence tomography,OCT)帧扫描数据,仅需使用单线扫数据。本发明对样本结构无依赖,通过使用新提出的色散失配指数即可快速地、可视化地实现色散失配探测和补偿。
附图说明
图1为本发明单线扫的色散失配探测的步骤流程图;
图2为本发明色散失配补偿的步骤流程图;
图3为本发明实施例实验装置结构图;
图4为本发明实施例未使用本发明和使用本发明的镜面玻璃的单线扫频谱对比图和OCT对比图;
图5为本发明实施例未使用本发明和使用本发明的光学窗口玻璃的单线扫频谱对比图和OCT对比图;
图6为本发明实施例未使用本发明和使用本发明的人体皮肤的单线扫频谱对比图和OCT对比图。
具体实施方式
下面将结合实施例,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
参阅图1-5,本发明提供一种技术方案:基于单线扫的光学相干层析成像色散失配探测及补偿方法,包括以下步骤:
S1:探测色散失配并进行可视化:对单线扫原始数据进行谱空间分析,利用加窗傅里叶变换得到频谱图;然后将频谱图进行二值化和骨架化处理;检测出频谱图中的线并对其进行数值拟合;计算出线的斜率,所述斜率表示色散失配指数;
S2:进行色散补偿:向单线扫原始数据中引入了补偿相位,所述补偿相位公式如下:
Φ=-a2(ω-ω0)2-a3(ω-ω0)3;
S3:将进行色散补偿后的频谱重复步骤S1的步骤,当频谱中斜率不为0时,重复步骤S2,当频谱中斜率为0时,步骤终止,通过重复步骤进行色散补偿实现了快速色散失配补偿。
进一步的,所述步骤S1还包括对单线扫原始数据进行谱空间分析的方法:线扫中的一个亮点对应单线扫频谱图中的一条线,所述线的线斜率值表示所述亮点成像深度下的色散失配状态,所以线扫不同点对应的频谱图中不同线用于用于可视化地显示不同成像深度的色散失配。
为更清楚的对本发明进行进一步的说明,本实施例构建了一个频谱域OCT(Spectral-domain OCT,SD-OCT)来进行说明,如图3所示;光源采用中心波长为850nm,半峰全宽为160nm的宽带光源(Superlum Broadlighters T-850-HP),输出光束被一个50:50的光纤耦合器(TW850R5A2,Thorlabs)分束到样本臂和参考臂。干涉光进入光谱仪(Cobra-S800,Wasatch Photonics),光谱仪探测到的信号被传输至计算机中。在实验中,光谱仪中具有2048个像素的相机和二维扫描振镜(Galvanometer,GVSM002-EC,Thorlabs)都通过计算机生成的触发信号进行同步。曝光时间为50μs,灵敏度为102dB。该装置在空气中的轴向分辨率为~2.5μm,横向分辨率为~3μm。
通过对不同样本进行实验,以证明该技术的性能。首先进行了镜面玻璃和0.5mm厚度的光学窗口玻璃的成像实验。在对镜面玻璃成像时,在样品臂添加了一个中性密度滤光片,以避免相机曝光饱和。图4和图5分别为镜面玻璃和光学窗口玻璃的频谱图及对应的OCT图像。实验获得了800条线扫的帧扫描数据,然后在帧扫描数据中选择中心线扫进行色散失配探测和补偿。但值得注意的是,任何来自帧扫描的线扫都可以用于该方法,不局限于中心线扫。为了得到线扫的频谱图,使用加窗傅里叶变换对线扫数据进行处理。当使用加窗傅里叶变换进行光谱空间分析时,如果使用0.25nm的步长移动高斯窗口,会导致频谱图的轴向分辨率降低。为了平衡光谱分辨率和空间分辨率,采用了步长为0.43nm的窗口。频谱图的轴向分辨率约为6.4μm@850nm。
图4中分别为未采用本发明进行色散补偿和采用本发明进行色散补偿的镜面玻璃的线扫频谱图和对应的OCT图像。从图4可以看出,当发生色散失配时,频谱图中的直线是倾斜的,而当良好补偿了色散失配时,频谱图中的直线斜率为0。由此可见,该技术仅用单线扫数据就能探测色散失配,并能实现探测结果的可视化,更进一步地,仅使用单线扫数据就能很好地进行色散补偿。
图5分别为厚度为0.5mm的光学窗口玻璃在未采用本发明进行色散补偿和采用本发明进行色散补偿的中心线扫频谱图和OCT图像。可以看到图5中显示有两条直线,两条直线分别表示了光学窗口的表面和底面的色散失配状态。从图5可以看出,本发明仅用单线扫数据就可以实现深度解析的色散失配探测,并实现探测结果的可视化。
为了证明本发明能够很好地对复杂结构样本进行色散失配探测和补偿,本实例提供了本发明在人体皮肤上实验的实施案例,如图6所示。图6中分别为未采用本发明进行色散补偿和采用本发明进行色散补偿的人体皮肤的中心线扫频谱图和对应的OCT图像。在图6中,我们可以看到在对人体皮肤进行实验时,出现了相同现象:色散失配时频谱图中的直线是倾斜的,而良好补偿色散失配时频谱图中的直线是倾斜的。人体皮肤中的汗腺在使用本发明进行色散补偿后清晰可见,而在人体皮肤中的汗腺在使用本发明进行色散补偿前却是模糊不清的,这说明本发明仅使用单线扫数据就可以很好地实现复杂结构样本的可视化快速色散失配探测和补偿。
为了验证本发明能够快速进行色散失配探测和补偿,本实例进一步比较了使用帧扫描的光谱分析和本发明提出的单线扫的色散失配及补偿技术的处理时间。使用镜面玻璃实验中的数据进行处理,通过在计算机工作站(Intel(R)Xeon(R)Silver 4114CPU@2.20GHz,64.0GB RAM)上运行MATLAB代码,使用帧扫描的方法耗时约为506.78秒,而本发明仅需1.11s即可得到补偿结果。
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (2)
1.基于单线扫的光学相干层析成像色散失配探测及补偿方法,其特征在于:包括以下步骤:
S1:探测色散失配并进行可视化:对单线扫原始数据进行谱空间分析,利用加窗傅里叶变换得到频谱图;然后将频谱图进行二值化和骨架化处理;检测出频谱图中的线并对其进行数值拟合;计算出线的斜率,所述斜率表示色散失配指数;
S2:进行色散补偿:向单线扫原始数据中引入了补偿相位,所述补偿相位公式如下:
Φ=-a2(ω-ω0)2-a3(ω-ω0)3;
S3:将进行色散补偿后的频谱重复步骤S1的步骤,当频谱中斜率不为0时,重复步骤S2,当频谱中斜率为0时,步骤终止。
2.根据权利要求1所述的基于单线扫的光学相干层析成像色散失配探测及补偿方法,其特征在于:所述步骤S1还包括对单线扫原始数据进行谱空间分析的方法:线扫中的一个亮点对应单线扫频谱图中的一条线,所述线的线斜率值表示所述亮点成像深度下的色散失配状态,线扫不同点对应的频谱图中不同线用于可视化的显示不同成像深度的色散失配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111226587.4A CN113984713A (zh) | 2021-10-21 | 2021-10-21 | 基于单线扫的光学相干层析成像色散失配探测及补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111226587.4A CN113984713A (zh) | 2021-10-21 | 2021-10-21 | 基于单线扫的光学相干层析成像色散失配探测及补偿方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113984713A true CN113984713A (zh) | 2022-01-28 |
Family
ID=79739933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111226587.4A Pending CN113984713A (zh) | 2021-10-21 | 2021-10-21 | 基于单线扫的光学相干层析成像色散失配探测及补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113984713A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130169971A1 (en) * | 2011-11-23 | 2013-07-04 | Wasatch Photonics, Inc. | Oct imaging system for curved samples |
CN104771144A (zh) * | 2015-04-13 | 2015-07-15 | 中国科学院上海光学精密机械研究所 | 频域光学相干层析成像系统的色散补偿方法 |
CN107661089A (zh) * | 2017-09-19 | 2018-02-06 | 北京工业大学 | 一种频域光学相干层析连续色散补偿成像方法和系统 |
CN111879730A (zh) * | 2020-08-14 | 2020-11-03 | 中国科学院上海光学精密机械研究所 | 基于矩形窗函数优化的光学相干层析成像信号处理方法 |
CN112597947A (zh) * | 2020-12-29 | 2021-04-02 | 电子科技大学 | 一种基于傅里叶域光学相干层析成像技术的色散补偿方法 |
CN112683848A (zh) * | 2020-12-21 | 2021-04-20 | 中国科学院上海光学精密机械研究所 | 光学相干层析成像系统色散补偿方法 |
-
2021
- 2021-10-21 CN CN202111226587.4A patent/CN113984713A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130169971A1 (en) * | 2011-11-23 | 2013-07-04 | Wasatch Photonics, Inc. | Oct imaging system for curved samples |
CN104771144A (zh) * | 2015-04-13 | 2015-07-15 | 中国科学院上海光学精密机械研究所 | 频域光学相干层析成像系统的色散补偿方法 |
CN107661089A (zh) * | 2017-09-19 | 2018-02-06 | 北京工业大学 | 一种频域光学相干层析连续色散补偿成像方法和系统 |
CN111879730A (zh) * | 2020-08-14 | 2020-11-03 | 中国科学院上海光学精密机械研究所 | 基于矩形窗函数优化的光学相干层析成像信号处理方法 |
CN112683848A (zh) * | 2020-12-21 | 2021-04-20 | 中国科学院上海光学精密机械研究所 | 光学相干层析成像系统色散补偿方法 |
CN112597947A (zh) * | 2020-12-29 | 2021-04-02 | 电子科技大学 | 一种基于傅里叶域光学相干层析成像技术的色散补偿方法 |
Non-Patent Citations (3)
Title |
---|
NORMAN LIPPOK ET AL.: "Single-shot speckle reduction and dispersion compensation in optical coherence tomography by compounding fractional Fourier domains", 《OPTICS LETTERS》 * |
姜盼秋 等: "谱域光学相干层析系统的色散补偿技术研究", 《光电工程》 * |
张仙玲: "光学相干层析成像系统与像质改善方法研究", 《中国优秀博士学位论文全文数据库 信息科技辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7859679B2 (en) | System, method and arrangement which can use spectral encoding heterodyne interferometry techniques for imaging | |
Li et al. | Robust reconstruction of local optic axis orientation with fiber-based polarization-sensitive optical coherence tomography | |
Zhang et al. | Graphics processing unit accelerated non-uniform fast Fourier transform for ultrahigh-speed, real-time Fourier-domain OCT | |
Desjardins et al. | Speckle reduction in OCT using massively-parallel detection and frequency-domain ranging | |
Yun et al. | High-speed spectral-domain optical coherence tomography at 1.3 µm wavelength | |
US7982879B2 (en) | Methods and systems for performing angle-resolved fourier-domain optical coherence tomography | |
Choi et al. | Spectral domain optical coherence tomography of multi-MHz A-scan rates at 1310 nm range and real-time 4D-display up to 41 volumes/second | |
Kuo et al. | Balanced detection for spectral domain optical coherence tomography | |
JP6840520B2 (ja) | 画像処理装置、撮像装置、画像処理方法及びプログラム | |
Pan et al. | Depth-dependent dispersion compensation for full-depth OCT image | |
US8797551B2 (en) | Compact, affordable optical test, measurement or imaging device | |
Shen et al. | Improving lateral resolution and image quality of optical coherence tomography by the multi-frame superresolution technique for 3D tissue imaging | |
Dubois et al. | Full-field optical coherence tomography | |
EP3055644B1 (en) | Wavelength-encoded tomography | |
Watanabe et al. | Three-dimensional imaging by ultrahigh-speed axial-lateral parallel time domain optical coherence tomography | |
Ho et al. | Mirau-type full-field optical coherence tomography with switchable partially spatially coherent illumination modes | |
Jelly et al. | Deep imaging with 1.3 µm dual-axis optical coherence tomography and an enhanced depth of focus | |
Ni et al. | Single A-line method for fast sample-structure-nondependent dispersion compensation of FD-OCT | |
CN113984713A (zh) | 基于单线扫的光学相干层析成像色散失配探测及补偿方法 | |
Liu et al. | Dispersion correction for optical coherence tomography by parameter estimation in fractional fourier domain | |
Steiner et al. | Automatic estimation of noise parameters in Fourier-domain optical coherence tomography cross sectional images using statistical information | |
Iyer et al. | Dispersion mismatch correction for evident chromatic anomaly in low coherence interferometry | |
US20230414095A1 (en) | Processing system for oct imaging, oct imaging system and method for oct imaging | |
Yurdakul et al. | Line-field coherent sensing with LED illumination | |
Yu et al. | Tri-band spectroscopic optical coherence tomography based on optical parametric amplification for lipid and vessel visualization |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220128 |
|
RJ01 | Rejection of invention patent application after publication |