CN105559756B - 基于全空间调制谱分割角度复合的微血管造影方法与系统 - Google Patents
基于全空间调制谱分割角度复合的微血管造影方法与系统 Download PDFInfo
- Publication number
- CN105559756B CN105559756B CN201610080511.8A CN201610080511A CN105559756B CN 105559756 B CN105559756 B CN 105559756B CN 201610080511 A CN201610080511 A CN 201610080511A CN 105559756 B CN105559756 B CN 105559756B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- modulation spectrum
- light
- arm
- sample
- 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
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 92
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000011218 segmentation Effects 0.000 title claims abstract description 29
- 150000001875 compounds Chemical class 0.000 title claims abstract description 19
- 230000003287 optical effect Effects 0.000 claims abstract description 27
- 230000017531 blood circulation Effects 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 15
- 238000005516 engineering process Methods 0.000 claims abstract description 15
- 238000002583 angiography Methods 0.000 claims abstract description 11
- 230000003595 spectral effect Effects 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 230000021615 conjugation Effects 0.000 claims abstract description 6
- 239000000835 fiber Substances 0.000 claims description 18
- 230000010287 polarization Effects 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 6
- 238000009738 saturating Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 239000000571 coke Substances 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 239000013307 optical fiber Substances 0.000 claims 1
- 238000012014 optical coherence tomography Methods 0.000 abstract description 32
- 230000003014 reinforcing effect Effects 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 50
- 241000208199 Buxus sempervirens Species 0.000 description 11
- 238000003384 imaging method Methods 0.000 description 10
- 238000002601 radiography Methods 0.000 description 10
- 230000007423 decrease Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 206010052128 Glare Diseases 0.000 description 2
- 238000000149 argon plasma sintering Methods 0.000 description 2
- 239000012472 biological sample Substances 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004313 glare Effects 0.000 description 2
- 230000000004 hemodynamic effect Effects 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 208000018152 Cerebral disease Diseases 0.000 description 1
- 208000010412 Glaucoma Diseases 0.000 description 1
- 206010064930 age-related macular degeneration Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008081 blood perfusion Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 210000004220 fundus oculi Anatomy 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 208000002780 macular degeneration Diseases 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000002969 morbid Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0275—Measuring blood flow using tracers, e.g. dye dilution
- A61B5/02755—Radioactive tracers
-
- 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/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
-
- 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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0077—Devices for viewing the surface of the body, e.g. camera, magnifying lens
- A61B5/0079—Devices for viewing the surface of the body, e.g. camera, magnifying lens using mirrors, i.e. for self-examination
-
- 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/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- 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
-
- 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/7271—Specific aspects of physiological measurement analysis
- A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
-
- 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/02044—Imaging in the frequency domain, e.g. by using a spectrometer
-
- 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/02083—Interferometers characterised by particular signal processing and presentation
-
- 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/0209—Low-coherence interferometers
- G01B9/02091—Tomographic interferometers, e.g. based on optical coherence
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Physiology (AREA)
- General Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Hematology (AREA)
- Cardiology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于全空间调制谱分割角度复合的微血管造影方法与系统。结合光学相干层析技术的三维空间分辨能力和动态散射技术的运动识别能力,实现无标记的三维光学微血管造影。利用空间频率域的横向扫描调制谱编码不同入射角度的探测光,通过分割调制谱获得入射角度分辨的相互独立的血管造影子图,实现多个空间角度复合的微血管造影图。通过在深度(z)域去除共轭镜像,重构复数值的OCT干涉光谱,在空间频率域得到全空间的调制谱,避免调制谱共轭镜像的干扰。本发明中的空间角度复合技术可以增强血管造影的对比度;空间频率域的全空间技术提供最大化的调制谱带宽和对比度增强效果;多个角度分辨的探测技术可以实现血流绝对流速的测量。
Description
技术领域
本发明涉及光学相干层析成像技术(Optical Coherence Tomography,OCT)以及在此基础之上的无标记三维光学微血管造影技术(OCT Angiography,Angio-OCT),尤其涉及基于全空间调制谱分割角度复合的微血管造影方法与系统。
技术背景
OCT技术是自从上世纪90年代以后逐步发展起来的在生物医学领域具有重大应用价值的一种新型成像技术。凭借着无标记物、非侵入性、非接触性、高成像分辨率和高探测灵敏度等优点,它吸引着越来越多的国内外科研工作者更加深入地对其进行研究,目前在医学临床中得到了广泛应用。OCT技术主要依据所探测到的由于生物样品的光学不均匀性所导致的样品背向散射光光强变化,获取样品内部的反射率信息,进而重构出样品的断层结构图像。然而通常在疾病的早期阶段,由于正常的和病变的生物组织间光散射特性的差异性较小,运用传统的结构型OCT较难加以检测并区分,故而在临床应用上其存在诸多的局限性。为了能够获得除生物样品组织形态结构之外更多的生理信息,功能型OCT技术得到了开发并广泛应用。
诸多的疾病(如青光眼、老年黄斑病变等眼底疾病,脑中风等脑疾病)与血流灌注的病态变化存在密切的相关性。如若能够实时监测病人血管的变化,并提供血管的三维影像,对该类疾病的早期诊断和控制具有重要的意义。Angio-OCT技术是一项具有巨大前景的可有助于血管疾病诊断的工具,其能够较好的区分静态的组织和动态的血流信号,实现了OCT在获取血管内血流信息方面的功能性拓展。相较于传统的造影技术,其优势体现在无需注入显影剂和利用X射线,并结合前面所述OCT技术的特点,能够实现微血管三维深度分辨的高对比度的成像。通常,血流对比度的模型建立是运用数学的手段分析光散射信号的时间统计特性,并用阈值分割动态的血流信号和静态的组织背景信号来实现的。然而由于动态和静态信号统计曲线之间存在重叠,导致动态和静态信号分割错误,血流对比度受到抑制。因此,采取有效的方法提高造影图的血流对比,有助于更加清晰地解释说明图像的特征。
基于Angio-OCT信号的时间统计特性研究得出:多个独立的血管造影子图的平均可以提高造影的对比度。参照一些消除散斑的方法,可以通过波长多样性、角度多样性和偏振多样性等方法来获取独立的造影子图。Jia等人提出了一种类似于波长多样性的光谱分割Angio-OCT方法。该方法将OCT干涉信号的全波长光谱分割成不同的子光谱,每个子光谱可产生独立的造影子图,通过复合形成新的血管造影图。然而相对于原来的全光谱,每个子光谱带宽变窄,导致图像的轴向分辨率下降。
发明内容
本发明针对现有技术的不足,提出了基于全空间调制谱分割角度复合的微血管造影方法与系统。
基于全空间调制谱分割角度复合的微血管造影方法,具体包括以下步骤:
1)通过在OCT深度(z)域去除共轭镜像,重构复数值的干涉光谱,在空间频率域得到全空间的横向扫描调制谱;
2)利用空间频率域的横向扫描调制谱编码不同入射角度的探测光,通过分割调制谱获得入射角度分辨的相互独立的OCT干涉子光谱;
3)结合2)中的调制谱分割技术与光学微血管造影技术,得到基于横向扫描调制谱分割的多个角度分辨的独立的血管造影子图;
4)结合1)中的全空间调制谱和2)中的调制谱分割技术和与光学微血管造影技术,获得基于全空间横向扫描调制谱分割的多个角度分辨的独立的血管造影子图;
5)对3)或4)中得到的血管造影子图进行平均,得到多个空间角度复合的血管造影图。
根据2)中的调制谱分割技术,实现基于横向扫描调制谱分割的多个角度分辨的血流绝对流速测量;
根据1)中的全空间调制谱和2)中的调制谱分割技术,实现基于全空间横向扫描调制谱分割的多个角度分辨的血流绝对流速测量。
本发明的目的是通过如下技术方案实现的:
基于全空间调制谱分割角度复合的微血管造影方法:
1)在OCT成像中,对应于一次横向扫描得到的干涉光谱实数信号,可表示为S(k,x),其中k代表波数,x代表横向的快扫方向。由于直接对光谱实数信号沿横向x方向作傅里叶变换,所得到的空间频率分布会产生正负频率的混叠,因此需要重构出复数值的OCT干涉光谱。为了重构出复数谱,对实数谱S(k,x)沿着k方向进行傅里叶变换,得到深度域(z)的空间结构信息。通常在谱域OCT中,一半的z空间对应约3mm的成像范围,能够满足绝大多数的应用要求,于是将待测样品完全置于零光程差的一边,以分辨复共轭镜像。通过去除z域中的共轭项,沿k方向进行傅里叶逆变换,最终得到复数值的光谱信号对该复数谱沿x方向进行傅立叶变换,在空间频率域获得全空间的横向扫描调制谱。
2)通常在OCT系统的样品臂中,为获得较高的横向成像分辨率,要求照射到物镜上的光束直径尺寸尽量大。因此,当一定宽度的准直光束照射到样品臂扫描振镜的转轴中心时,光束距离转轴中心存在一定的偏移量δ。而偏移量δ会引入样品臂光程的调制,导致在空间频率域的横向扫描调制,其调制频率fm与偏移量δ成线性函数关系:
其中,k表示光源的中心波数,ω表示扫描振镜的角速度。因此,在全空间的横向扫描调制谱中,不同的调制频率对应着不同入射角度的探测光。对该调制谱进行分割,可获得入射角度分辨的相互独立的OCT干涉子光谱。分割得到的子光谱个数及子光谱间的重叠视最终的图像效果而定:分割的个数越多会导致横向分辨率减小得越多,此外还增加了处理算法的复杂度;子光谱间相互重叠的部分越大,会导致子光谱相互之间的独立性下降,影响编码不同入射角度的探测光。
3)对获得的入射角度分辨的相互独立的子光谱,利用目前的微血管造影技术,如幅值差分法、复数差分法、去互相关法和散斑方差法等,分别产生相对应的独立的造影子图。对这些空间角度分辨的造影子图进行复合平均,得到新的微血管造影图。
传统的多普勒频移测流速法只能测量与探测光束方向平行的速度分量。要测量血流的绝对速度,还需测得探测光束方向与血流方向的夹角(即多普勒角)。基于全空间调制谱分割,能够分辨不同入射角度的探测光束。对分割得到的角度分辨的调制子光谱,利用多普勒血流速度计算方法,可以获得不同入射角度探测方向的流速分量,进而根据几何关系可确定绝对的血流速度值。
基于全空间调制谱分割角度复合的微血管造影系统,包括低相干宽带光源、光环行器、光纤耦合器、参考臂、样品臂、光谱仪和信号处理模块;其中参考臂包括参考臂偏振控制器件、参考臂准直透镜、参考臂聚焦透镜和平面反射镜;样品臂包括样品臂偏振控制器、样品臂准直透镜、正交扫描振镜和样品臂聚焦物镜;光谱仪包括光谱仪准直透镜、闪耀光栅、傅里叶透镜和线阵CMOS相机;
低相干宽带光源经过光环行器与光纤耦合器一侧的输入端相连接;光纤耦合器另一侧的其中一个输出端口经过参考臂偏振控制器件与参考臂准直透镜相连接;参考臂聚焦透镜的光轴与参考臂准直透镜的光轴重合,平面反射镜置于参考臂聚焦透镜的焦平面处;光纤耦合器另一侧的另外一个输出端口,经过样品臂偏振控制器与样品臂准直透镜相连接;正交扫描振镜的第一扫描镜转轴中心位于样品臂准直透镜的光轴上,第二扫描镜转轴中心位于样品臂聚焦物镜的光轴上,待测样品置于样品臂聚焦物镜的焦平面处;光环行器的另一个输出端口与光谱仪准直透镜相连接,闪耀光栅根据分光原理放置在光谱仪准直透镜的出射光路上,傅里叶透镜置于闪耀光栅的出射光路上,线阵CMOS相机的采光面与傅里叶透镜的后焦平面重合;线阵CMOS相机后面接信号处理模块;
低相干宽带光源发出的光经过光环行器后入射到光纤耦合器,出射光分成两部分,一部分进入参考臂,经准直聚焦后照射到平面反射镜;另一部分进入样品臂,经过准直聚焦后照射到待测样品上;样品臂的正交扫描振镜实现样品臂光束对待测样品的三维扫描;参考臂反射镜反射回的光与待测样品背向散射回的光在光纤耦合器处发生干涉,出射的干涉光经过光谱仪后由信号处理模块采集处理。
与背景技术相比,本发明具有的有益效果是:
1)利用本发明获得的多个空间角度复合的微血管造影图,具有增强的对比度和血管连通性
2)相比于背景技术中谈及的类似于波长多样性的光谱分割Angio-OCT方法,其对光谱分割会导致轴向分辨率的下降。本发明对轴向分辨率不产生影响。
3)本发明谈及的全空间的横向调制谱分割技术,可以提供最大化的空间频率域的横向调制谱,消除由于调制谱分割导致的横向分辨率下降。
4)利用本发明谈及的基于全空间调制谱分割的多个角度分辨的探测技术,可以在单光束单次测量的情况下实现血流绝对流速的测量。
附图说明
图1为样品臂扫描光束示意图;
图2为偏移量δ引入的调制频率fm图;
图3为横向扫描调制谱正负频率混叠图;
图4是本发明的方法实现流程图;
图5是本发明的成像系统示意图。
图1中:①、探测光束;②、扫描振镜;③、振镜转轴中心;④、聚焦物镜。
图2中:偏移量δ引入的调制频率fm。
图3中:正负频率混叠的横向扫描调制谱。
图4中:A表示OCT干涉光谱实数信号;B表示深度(z)域的空间结构信息;C表示去除共轭项后的空间结构信息;D表示复数值的干涉光谱;E表示空间频率域的全空间横向扫描调制谱;F1至Fn分别表示分割得到的调制子光谱;G1至Gn分别表示z域空间结构信号;H1至Hn分别表示光学造影子图;I表示光学造影图;
图5中:1、低相干宽带光源;2、光环行器;3、光纤耦合器;4、参考臂偏振控制器件;5、参考臂准直透镜;6、参考臂聚焦透镜;7、平面反射镜;8、样品臂偏振控制器件;9、样品臂准直透镜;10、正交扫描振镜;11、样品臂聚焦物镜;12、待测样品;13、准直透镜;14、闪耀光栅;15、傅里叶透镜;16、线阵CMOS相机;17、信号处理模块。
具体实施方式
下面结合附图和实施例子对本发明作进一步的说明。
图1、图2和图3所示为本发明的横向调制原理示意图。下面对图1、图2和图3中的示意图作详细说明。
1)通常在OCT系统的样品臂中,为获得较高的横向成像分辨率,要求照射到物镜上的光束直径尺寸尽量大。因此,在图1中,当一定宽度的准直光束照射到样品臂扫描振镜的转轴中心时,光束距离转轴中心存在一定的偏移量δ。而偏移量δ会引入样品臂光程的调制,导致在空间频率域的横向扫描调制(如图2),其调制频率fm与偏移量δ成线性函数关系:
其中,k表示光源的中心波数,ω表示扫描振镜的角速度。因此,在全空间的横向扫描调制谱中,不同的调制频率对应着不同入射角度的探测光。
2)在OCT成像中,对应于一次横向扫描得到的干涉光谱为实数值信号。由于直接对光谱实数信号沿横向快扫方向作傅里叶变换,所得到的空间频率分布会产生正负频率的混叠(如图3),因此需要重构出复数值的OCT干涉光谱。
图4所示为本发明的方法的具体实现流程图。下面对图4中的过程进行详细说明:
1)在OCT成像中,对应于一次横向扫描得到的是干涉光谱实数信号(即方框A)。为了重构出复数值的光谱信号,对实数谱沿着波数k方向进行傅里叶变换,得到深度z域的空间结构信息(即方框B)。通常在谱域OCT中,一半的z空间对应约3mm的成像范围,能够满足绝大多数的应用要求,于是将待测样品完全置于零光程差的一边,以分辨复共轭镜像。通过去除z域中的共轭项,得到新的z域空间结构信息(即方框C)。而后沿着k方向进行傅里叶逆变换,最终得到复数值的光谱信号(即方框D)。对该复数谱沿x方向进行傅立叶变换,在空间频率域获得全空间的横向扫描调制谱(即方框E)。
2)对得到的横向调制谱进行分割,可获得入射角度分辨的相互独立的OCT干涉子光谱(即方框F1至Fn)。分割得到的子光谱个数n及子光谱间的重叠视最终的图像效果而定:分割的个数越多会导致横向分辨率减小得越多,此外还增加了处理算法的复杂度;子光谱间相互重叠的部分越大,会导致子光谱相互之间的独立性下降,影响编码不同入射角度的探测光。
3)对获得的入射角度分辨的相互独立的子光谱(即方框F1至Fn),分别沿着x方向进行傅里叶逆变换,再沿着k方向进行傅里叶变换,得到z域的空间结构信息(即方框G1至Gn)。利用目前的光学微血管造影技术,如幅值差分法、互相关法等,分别产生相对应的独立的造影子图(即方框H1和Hn)。对这些空间角度分辨的造影子图进行复合平均,得到新的微血管造影图(即方框I),其血流对比度得到了提高。
传统的多普勒频移测流速法只能测量与探测光束方向平行的速度分量。要测量血流的绝对速度,还需测得探测光束方向与血流方向的夹角(即多普勒角)。基于全空间调制谱分割,能够分辨不同入射角度的探测光束。对分割得到的调制子光谱(即方框F1至Fn),利用多普勒血流速度计算方法,可以获得不同入射角度探测方向的流速分量,根据几何关系可确定绝对的血流速度值。
图5所示为本发明的成像系统示意图。下面对图5进行详细说明。
基于全空间调制谱分割角度复合的微血管造影系统,包括低相干宽带光源1、光环行器2、光纤耦合器3、参考臂、样品臂、光谱仪和信号处理模块17;其中参考臂包括参考臂偏振控制器件4、参考臂准直透镜5、参考臂聚焦透镜6和平面反射镜7;样品臂包括样品臂偏振控制器8、样品臂准直透镜9、正交扫描振镜10和样品臂聚焦物镜11;光谱仪包括光谱仪准直透镜13、闪耀光栅14、傅里叶透镜15和线阵CMOS相机;
其特征在于:低相干宽带光源1经过光环行器2与光纤耦合器3一侧的输入端相连接;光纤耦合器3另一侧的其中一个输出端口经过参考臂偏振控制器件4与参考臂准直透镜5相连接;参考臂聚焦透镜6的光轴与参考臂准直透镜5的光轴重合,平面反射镜7置于参考臂聚焦透镜6的焦平面处;光纤耦合器3另一侧的另外一个输出端口,经过样品臂偏振控制器8与样品臂准直透镜9相连接;正交扫描振镜10的第一扫描镜转轴中心位于样品臂准直透镜9的光轴上,第二扫描镜转轴中心位于样品臂聚焦物镜11的光轴上,待测样品12置于样品臂聚焦物镜11的焦平面处;光环行器2的另一个输出端口与光谱仪准直透镜13相连接,闪耀光栅14根据分光原理放置在光谱仪准直透镜13的出射光路上,傅里叶透镜15置于闪耀光栅14的出射光路上,线阵CMOS相机16的采光面与傅里叶透镜15的后焦平面重合;线阵CMOS相机16后面接信号处理模块17;
低相干宽带光源1发出的光经过光环行器2后入射到光纤耦合器3,出射光分成两部分,一部分进入参考臂,经准直聚焦后照射到平面反射镜7;另一部分进入样品臂,经过准直聚焦后照射到待测样品12上;样品臂的正交扫描振镜10实现样品臂光束对待测样品12的三维扫描;参考臂反射镜7反射回的光与待测样品12背向散射回的光在光纤耦合器3处发生干涉,出射的干涉光经过光谱仪后由信号处理模块17采集处理。
Claims (4)
1.基于全空间调制谱分割角度复合的微血管造影方法,其特征在于,该方法具体包括以下步骤:
1)、在OCT成像中,对应于一次横向扫描得到的干涉光谱实数信号S(k,x),其中k代表波数,x代表横向的快扫方向,对实数谱S(k,x)沿着k方向进行傅里叶变换,得到深度域(z)的空间结构信息,去除z域中的共轭项,沿k方向进行傅里叶逆变换,最终得到复数值的光谱信号对该沿x方向进行傅立叶变换,在空间频率域获得全空间的横向扫描调制谱;
2)、利用空间频率域的横向扫描调制谱编码不同入射角度的探测光,通过分割调制谱获得入射角度分辨的相互独立的OCT干涉子光谱;
3)、结合2)中的调制谱分割技术与光学微血管造影技术,得到基于横向扫描调制谱分割的多个角度分辨的独立的血管造影子图;
4)、结合1)中的全空间调制谱和2)中的调制谱分割技术和光学微血管造影技术,获得基于全空间横向扫描调制谱分割的多个角度分辨的独立的血管造影子图;
5)、对3)或4)中得到的血管造影子图进行平均,得到多个空间角度复合的血管造影图。
2.根据权利要求1所述的基于全空间调制谱分割角度复合的微血管造影方法,其特征在于:根据2)中的调制谱分割技术,实现基于横向扫描调制谱分割的多个角度分辨的血流绝对流速测量。
3.根据权利要求1所述的基于全空间调制谱分割角度复合的微血管造影方法,其特征在于:根据1)中的全空间调制谱和2)中的调制谱分割技术,实现基于全空间横向扫描调制谱分割的多个角度分辨的血流绝对流速测量。
4.根据权利要求1所述的基于全空间调制谱分割角度复合的微血管造影方法,其所基于的系统包括低相干宽带光源、光环行器、光纤耦合器、参考臂、样品臂、光谱仪和信号处理模块;其中参考臂包括参考臂偏振控制器件、参考臂准直透镜、参考臂聚焦透镜和平面反射镜;样品臂包括样品臂偏振控制器、样品臂准直透镜、正交扫描振镜和样品臂聚焦物镜;光谱仪包括光谱仪准直透镜、闪耀光栅、傅里叶透镜和线阵CMOS相机;
其特征在于:低相干宽带光源经过光环行器与光纤耦合器一侧的输入端相连接;光纤耦合器另一侧的其中一个输出端口经过参考臂偏振控制器件与参考臂准直透镜相连接;参考臂聚焦透镜的光轴与参考臂准直透镜的光轴重合,平面反射镜置于参考臂聚焦透镜的焦平面处;光纤耦合器另一侧的另外一个输出端口,经过样品臂偏振控制器与样品臂准直透镜相连接;正交扫描振镜的第一扫描镜转轴中心位于样品臂准直透镜的光轴上,第二扫描镜转轴中心位于样品臂聚焦物镜的光轴上,待测样品置于样品臂聚焦物镜的焦平面处;光环行器的另一个输出端口与光谱仪准直透镜相连接,闪耀光栅根据分光原理放置在光谱仪准直透镜的出射光路上,傅里叶透镜置于闪耀光栅的出射光路上,线阵CMOS相机的采光面与傅里叶透镜的后焦平面重合;线阵CMOS相机后面接信号处理模块;
低相干宽带光源发出的光经过光环行器后入射到光纤耦合器,出射光分成两部分,一部分进入参考臂,经准直聚焦后照射到平面反射镜;另一部分进入样品臂,经过准直聚焦后照射到待测样品上;样品臂的正交扫描振镜实现样品臂光束对待测样品的三维扫描;参考臂反射镜反射回的光与待测样品背向散射回的光在光纤耦合器处发生干涉,出射的干涉光经过光谱仪后由信号处理模块采集处理。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610080511.8A CN105559756B (zh) | 2016-02-05 | 2016-02-05 | 基于全空间调制谱分割角度复合的微血管造影方法与系统 |
US15/547,922 US10667707B2 (en) | 2016-02-05 | 2016-04-20 | Microangiography method and system based on full-space modulation spectrum splitting and angle compounding |
PCT/CN2016/079693 WO2017133083A1 (zh) | 2016-02-05 | 2016-04-20 | 基于全空间调制谱分割角度复合的微血管造影方法与系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610080511.8A CN105559756B (zh) | 2016-02-05 | 2016-02-05 | 基于全空间调制谱分割角度复合的微血管造影方法与系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105559756A CN105559756A (zh) | 2016-05-11 |
CN105559756B true CN105559756B (zh) | 2019-11-15 |
Family
ID=55870810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610080511.8A Active CN105559756B (zh) | 2016-02-05 | 2016-02-05 | 基于全空间调制谱分割角度复合的微血管造影方法与系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10667707B2 (zh) |
CN (1) | CN105559756B (zh) |
WO (1) | WO2017133083A1 (zh) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105852816A (zh) * | 2016-05-16 | 2016-08-17 | 浙江大学 | 一种全通道调制编码的角度复合散斑降噪方法与系统 |
CN106236205A (zh) * | 2016-07-27 | 2016-12-21 | 深圳市中科微光医疗器械技术有限公司 | 一种基于近红外相干断层成像技术的血管导航系统及方法 |
CN106137134B (zh) * | 2016-08-08 | 2023-05-12 | 浙江大学 | 多角度复合的血流成像方法及系统 |
CN107822600A (zh) * | 2017-12-14 | 2018-03-23 | 天津海仁医疗技术有限公司 | 应用于血管性皮肤病检测、定位的装置和系统及工作方法 |
CN108852285B (zh) * | 2018-04-16 | 2020-10-23 | 清华大学深圳研究生院 | 一种基于频域oct的微血管造影方法 |
CN108956533A (zh) * | 2018-07-13 | 2018-12-07 | 福州大学 | 一种可消除多色误差的oct共轭镜像去除装置及方法 |
CN109008982B (zh) * | 2018-08-09 | 2021-01-29 | 中南大学 | 一种多普勒光学相干层析的绝对流速测量装置与方法 |
CN110873956B (zh) * | 2018-08-29 | 2021-07-23 | 暨南大学 | 一种超高速正交偏振成像装置及方法 |
CN109297599B (zh) * | 2018-11-09 | 2023-08-18 | 福州大学 | 可消除oct共轭镜像的定相差线阵谱域oct装置及方法 |
CN109691971A (zh) * | 2018-12-29 | 2019-04-30 | 佛山科学技术学院 | 眼前节组织动态参数快速成像的光学相干层析成像系统 |
CN109620134B (zh) * | 2019-01-21 | 2020-05-22 | 浙江大学 | 基于光纤阵列多通道并行探测的微血管造影方法和系统 |
CN109965839A (zh) * | 2019-04-08 | 2019-07-05 | 广东唯仁医疗科技有限公司 | 眼前节oct成像装置 |
CN110013212B (zh) * | 2019-04-23 | 2024-04-16 | 广东唯仁医疗科技有限公司 | 基于光学相干断层扫描的多参数、多功能眼睛测量仪 |
CN110151133B (zh) * | 2019-05-24 | 2021-10-01 | 哈尔滨工业大学 | 基于图像分割与时频信息融合的乳腺光学成像装置及方法 |
CN111134614A (zh) * | 2019-12-20 | 2020-05-12 | 苏州大学 | 基于oct测量眼球血管中血流的绝对速度的方法及系统 |
CN111568373A (zh) * | 2020-04-20 | 2020-08-25 | 南京师范大学 | 一种重复扫描的octa毛细血管网成像方法 |
CN112075925A (zh) * | 2020-09-21 | 2020-12-15 | 北京脑科学与类脑研究中心 | 基于散斑原理的荧光成像照明装置、成像系统及成像方法 |
CN112535465B (zh) * | 2020-11-03 | 2023-05-30 | 佛山科学技术学院 | 一种基于片层光的三维血流速度成像方法及装置 |
CN112396622B (zh) * | 2020-11-24 | 2023-10-31 | 浙江大学 | 基于多维特征空间的微血流图像分割量化方法和系统 |
CN112842286B (zh) * | 2020-12-31 | 2022-08-05 | 温州医科大学 | 一种基于功能oct的无创血流粘滞度测量方法 |
CN113116303A (zh) * | 2021-04-16 | 2021-07-16 | 佛山科学技术学院 | 一种甲襞微循环血管三维扫描分析系统及方法 |
CN113267142B (zh) * | 2021-05-17 | 2022-08-19 | 东北大学秦皇岛分校 | 一种表面轮廓成像装置及成像方法 |
CN113665150B (zh) * | 2021-08-31 | 2022-07-12 | 广州永士达医疗科技有限责任公司 | 一种软管制作方法及探头软管 |
WO2023211382A1 (en) * | 2022-04-25 | 2023-11-02 | Nanyang Technological University | Phase-restoring translational shifting of multidimensional images |
CN115089118A (zh) * | 2022-06-20 | 2022-09-23 | 暨南大学附属第一医院(广州华侨医院) | 一种双模态融合的血管成像系统 |
CN115639198B (zh) * | 2022-11-15 | 2024-05-07 | 南京理工大学 | 一种全场光学时空相干编码动态体成像装置及方法 |
CN117297555B (zh) * | 2023-11-29 | 2024-02-09 | 北京理工大学 | 一种基于并行光计算的光谱处理系统及在oct中的应用 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101803908A (zh) * | 2010-03-01 | 2010-08-18 | 浙江大学 | 基于色散调制的无镜像光学频域成像系统及方法 |
CN103267732A (zh) * | 2013-05-21 | 2013-08-28 | 浙江大学 | 基于移动光栅空间载频谱域oct全量程成像方法及系统 |
CN104411233A (zh) * | 2012-02-03 | 2015-03-11 | 俄勒冈健康科学大学 | 体内光学流成像 |
CN205458608U (zh) * | 2016-02-05 | 2016-08-17 | 浙江大学 | 基于全空间调制谱分割角度复合的微血管造影系统 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9086264B2 (en) * | 2007-06-06 | 2015-07-21 | Oregon Health & Science University | Polarization sensitive spectral domain OCT using an interference signal modulated at a constant frequency and a two-path reference arm with a single reference mirror |
JP5675268B2 (ja) * | 2010-10-21 | 2015-02-25 | キヤノン株式会社 | 光干渉断層撮像装置、光干渉断層撮像方法、補償方法およびプログラム |
WO2013081902A1 (en) * | 2011-11-30 | 2013-06-06 | University Of Miami | System and method for improving image quality in vivo oct imaging |
US8896841B2 (en) * | 2012-03-13 | 2014-11-25 | Kabushiki Kaisha Topcon | Optical imaging method and optical imaging apparatus using optical coherence tomography |
US9016862B2 (en) * | 2012-05-10 | 2015-04-28 | Sonomed Ip Holdings, Inc. | Multimodality correlation of optical coherence tomography using secondary reference images |
CN103565405B (zh) * | 2013-11-15 | 2015-12-09 | 浙江大学 | 基于分段光谱光程编码的谱域oct探测方法 |
CN104523239B (zh) * | 2015-01-12 | 2017-02-22 | 南京理工大学 | 全深度谱域光学相干层析成像装置及方法 |
-
2016
- 2016-02-05 CN CN201610080511.8A patent/CN105559756B/zh active Active
- 2016-04-20 US US15/547,922 patent/US10667707B2/en active Active
- 2016-04-20 WO PCT/CN2016/079693 patent/WO2017133083A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101803908A (zh) * | 2010-03-01 | 2010-08-18 | 浙江大学 | 基于色散调制的无镜像光学频域成像系统及方法 |
CN104411233A (zh) * | 2012-02-03 | 2015-03-11 | 俄勒冈健康科学大学 | 体内光学流成像 |
CN103267732A (zh) * | 2013-05-21 | 2013-08-28 | 浙江大学 | 基于移动光栅空间载频谱域oct全量程成像方法及系统 |
CN205458608U (zh) * | 2016-02-05 | 2016-08-17 | 浙江大学 | 基于全空间调制谱分割角度复合的微血管造影系统 |
Non-Patent Citations (1)
Title |
---|
全量程谱域光学相干层析成像方法研究;洪威;《中国优秀硕士学位论文全文数据库基础科学辑》;20130815(第8期);正文第13-14页,图2.2 * |
Also Published As
Publication number | Publication date |
---|---|
US20180070842A1 (en) | 2018-03-15 |
WO2017133083A1 (zh) | 2017-08-10 |
US10667707B2 (en) | 2020-06-02 |
CN105559756A (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105559756B (zh) | 基于全空间调制谱分割角度复合的微血管造影方法与系统 | |
CN107615005B (zh) | 高分辨率3-d谱域光学成像设备和方法 | |
JP5887006B2 (ja) | 角度分解型のフーリエドメイン光干渉断層撮影法を遂行する方法及びシステム | |
CN205458608U (zh) | 基于全空间调制谱分割角度复合的微血管造影系统 | |
JP4344829B2 (ja) | 偏光感受光画像計測装置 | |
CN107595250B (zh) | 基于运动与图形混合对比度的血流成像方法与系统 | |
JP5787255B2 (ja) | Ps−octの計測データを補正するプログラム及び該プログラムを搭載したps−octシステム | |
CN105147241B (zh) | 基于双空间载频技术拓展oct成像深度的方法与系统 | |
US7995814B2 (en) | Dynamic motion contrast and transverse flow estimation using optical coherence tomography | |
CN106137134B (zh) | 多角度复合的血流成像方法及系统 | |
US9226655B2 (en) | Image processing apparatus and image processing method | |
CN105342568B (zh) | 联合相位和幅值的光学相干造影方法及系统 | |
CN107788950B (zh) | 基于自适应阈值分割的血流成像方法与系统 | |
CN208837916U (zh) | 一种血流成像系统 | |
CN109115804A (zh) | 一种定量检测玻璃亚表面缺陷的装置及方法 | |
CN105996999B (zh) | 基于oct测量样品深度分辨衰减系数的方法和系统 | |
JPWO2010143601A1 (ja) | 2ビーム型光コヒーレンストモグラフィー装置 | |
CN108245130B (zh) | 一种光学相干断层血管造影装置及方法 | |
EP3111833A1 (en) | Method and apparatus for performing multidimensional velocity measurements based on amplitude and phase signals in the field of optical interferometry | |
CN109620134A (zh) | 基于光纤阵列多通道并行探测的微血管造影方法和系统 | |
CN102657519A (zh) | 基于oct的大动态范围流速的高灵敏度测量系统及方法 | |
CN106491078B (zh) | 去除血流图像中组织抖动噪声的方法及装置 | |
JP2012502262A (ja) | 半透明物質の画像化の方法および装置 | |
CN104545872B (zh) | 基于线性相关系数来重构三维微血流分布的方法及装置 | |
JP2010151684A (ja) | 局所的な複屈折情報を抽出可能な偏光感受光画像計測装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |