CN111436909B - 一种活体组织的光学相干层析成像系统及方法 - Google Patents
一种活体组织的光学相干层析成像系统及方法 Download PDFInfo
- Publication number
- CN111436909B CN111436909B CN202010348181.2A CN202010348181A CN111436909B CN 111436909 B CN111436909 B CN 111436909B CN 202010348181 A CN202010348181 A CN 202010348181A CN 111436909 B CN111436909 B CN 111436909B
- Authority
- CN
- China
- Prior art keywords
- coherence tomography
- optical coherence
- light
- optical
- living tissue
- 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
- 238000012014 optical coherence tomography Methods 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000003384 imaging method Methods 0.000 claims abstract description 52
- 230000017531 blood circulation Effects 0.000 claims abstract description 35
- 108010076504 Protein Sorting Signals Proteins 0.000 claims abstract description 30
- 238000005070 sampling Methods 0.000 claims abstract description 27
- 230000003287 optical effect Effects 0.000 claims description 49
- 238000012545 processing Methods 0.000 claims description 23
- 230000001427 coherent effect Effects 0.000 claims description 21
- 230000010287 polarization Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000002601 radiography Methods 0.000 abstract description 5
- 210000001519 tissue Anatomy 0.000 description 51
- 210000004556 brain Anatomy 0.000 description 11
- 238000002583 angiography Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 10
- 210000004204 blood vessel Anatomy 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 206010064930 age-related macular degeneration Diseases 0.000 description 5
- 210000003710 cerebral cortex Anatomy 0.000 description 5
- 208000002780 macular degeneration Diseases 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000000799 fluorescence microscopy Methods 0.000 description 4
- 230000002792 vascular Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 208000002026 familial multiple nevi flammei Diseases 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 210000001525 retina Anatomy 0.000 description 3
- 210000003491 skin Anatomy 0.000 description 3
- 201000004569 Blindness Diseases 0.000 description 2
- 206010067193 Naevus flammeus Diseases 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000003925 brain function Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000005590 Choroidal Neovascularization Diseases 0.000 description 1
- 206010060823 Choroidal neovascularisation Diseases 0.000 description 1
- 206010061818 Disease progression Diseases 0.000 description 1
- 208000006787 Port-Wine Stain Diseases 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 230000005750 disease progression Effects 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000007850 fluorescent dye Substances 0.000 description 1
- 102000034287 fluorescent proteins Human genes 0.000 description 1
- 108091006047 fluorescent proteins Proteins 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000011503 in vivo imaging Methods 0.000 description 1
- 230000002779 inactivation Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000653 nervous system Anatomy 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000002428 photodynamic therapy Methods 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 210000001210 retinal vessel Anatomy 0.000 description 1
- 210000004761 scalp Anatomy 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 230000004393 visual impairment Effects 0.000 description 1
Images
Classifications
-
- 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/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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/40—Animals
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Cardiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Vascular Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Hematology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供了一种活体组织的光学相干层析成像系统及方法，通过在所述活体组织的同一表面位置进行M次采样，M≥2；获得随波长变化的干涉信号序列，将干涉信号序列进行转换，变为随深度变化的复数信号序列，依据复数信号序列进行简单的运算，即可得到数值在0‑1之间的血流速度参数，无需对血流速度参数进行归一化处理，进而可显示不同位数的血流造影图像。该光学相干层析成像方法可以简单快速的实现对活体组织的血流造影成像。
Description
技术领域
本发明涉及活体组织检测技术领域，更具体地说，涉及一种活体组织的光学相干层析成像(Optical Coherence Tomography，OCT)系统及方法。
背景技术
随着科学技术的不断发展，在医学方面的研究也在不断进步，基于活体组织的检测，例如，对视网膜、脑皮层或皮肤等活体组织进行功能检测时，需要对其血管进行成像。
但是，目前的成像方式比较复杂，那么如何提供一种快速对活体组织进行成像的技术，是本领域技术人员亟待解决的问题。
发明内容
有鉴于此，为解决上述问题，本发明提供了一种活体组织的光学相干层析成像系统及方法，技术方案如下：
一种活体组织的光学相干层析成像方法，所述光学相干层析成像方法包括：
在所述活体组织的同一表面位置进行M次采样，M≥2；
获取M次采样过程中，随波长变化的干涉信号序列；
依据所述干涉信号序列，获得随深度变化的复数信号序列；
依据所述复数信号序列，获取此表面位置不同深度的血流速度参数，所述血流速度参数的值在0-1之间；
扫描不同表面位置后，依据所述血流速度参数，显示不同位数的血流造影图像。
优选的，在上述光学相干层析成像方法中，所述在所述活体组织的同一表面位置进行M次采样，包括：
OCT光束在所述活体组织的预设表面位置形成光斑；
在所述光斑的固定落点处，进行M次采样。
优选的，在上述光学相干层析成像方法中，所述在所述活体组织的预设区域内M个不同表面位置进行1次采样，包括：
OCT光束在所述活体组织的预设区域进行侧向位移；
在M个不同表面位置上均进行1次采样；
所述预设区域的尺寸小于等于所述OCT光束的光斑直径。
优选的，在上述光学相干层析成像方法中，所述深度变化的方向为OCT光束的传播方向，表面位置为OCT光束垂直平面的位置。
一种活体组织的光学相干层析成像系统，所述光学相干层析成像系统包括：光源发射装置、光学耦合器、参考光路、样品光路和信号处理成像装置；
所述光源发射装置用于出射弱相干光；
所述光学耦合器用于将所述弱相干光分束，一部分入射至所述参考光路，另一部分入射至所述样品光路；
所述参考光路用于将光束反射回所述光学耦合器，所述样品光路用于将被样品散射后的光束输送至所述光学耦合器，二者在光学耦合器中发生干涉；
所述信号处理成像装置获取干涉信号，进行成像。
优选的，在上述光学相干层析成像系统中，所述光源发射装置包括：
扫频光源，用于不同时刻输出不同波长的弱相干光；
偏振控制器，用于对所述激光的偏振态进行处理。
优选的，在上述光学相干层析成像系统中，所述信号处理成像装置包括：
光电探测器，用于检测所述干涉信号；
上位机，用于依据所述干涉信号进行成像。
优选的，在上述光学相干层析成像系统中，所述光源发射装置包括：
连续光谱光源，用于在同一时刻输出不同波长的光束；
光学隔离器，用于实现光束的单向通过。
优选的，在上述光学相干层析成像系统中，所述信号处理成像装置包括：
光栅，用于将不同波长的干涉光在空间上分开；
相机，用于检测所述干涉信号；
上位机，用于依据所述干涉信号进行成像。
优选的，在上述光学相干层析成像系统中，所述参考光路包括：第一透镜和反射镜；
所述样品光路包括：第二透镜和样品。
相较于现有技术，本发明实现的有益效果为：
本发明提供的一种活体组织的光学相干层析成像方法，通过在所述活体组织的同一表面位置进行M次采样，M≥2；获得随波长变化的干涉信号序列，将干涉信号序列进行转换，变为随深度变化的复数信号序列，依据复数信号序列进行简单的运算，即可得到数值在0-1之间的血流速度参数，也就是说无需对血流速度参数进行归一化处理，进而可显示不同位数的血流造影图像。
也就是说，该光学相干层析成像方法可以简单快速的实现对活体组织的血流造影成像。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本发明实施例提供的一种活体组织的光学相干层析成像方法的流程示意图；
图2为本发明实施例提供的OCT光束的坐标系示意图；
图3为本发明实施例提供的信号处理原理示意图；
图4为本发明实施例提供的一种在x-y平面内某一点重复采样示意图；
图5为本发明实施例提供的一种在同一y轴，沿x轴往复扫描的示意图；
图6为本发明实施例提供的一种密集扫描的示意图；
图7为本发明实施例提供的利用幅度差分算法的大鼠脑皮层的OCT血流造影投影图；
图8为本发明实施例提供的利用复数差分算法的大鼠脑皮层的OCT血流造影投影图；
图9为本发明实施例提供的利用幅度差分和算法的大鼠脑皮层的OCT血流造影投影图；
图10为本发明实施例提供的利用复数差分和算法的大鼠脑皮层的OCT血流造影投影图；
图11为本发明实施例提供的一种活体组织的光学相干层析成像系统的结构示意图；
图12为本发明实施例提供的另一种活体组织的光学相干层析成像系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
血流造影可以对视网膜、脑皮层和皮肤等活体组织的血管网络进行成像，用于活体组织功能的研究和疾病的诊断。
例如，脑皮层的血流造影可以用于大脑功能研究。大脑是神经系统最重要的功能单元，由相互连接的大量神经元和为神经网络提供血液的血管组成。脑皮层血液供应的变化可造成大脑神经元功能的激活或失活，改变大脑神经网络的结构，并调整信号传输的模式。因此，脑皮层的血流造影成像对于大脑功能的理解和疾病发展的研究具有十分重要的意义。
年龄相关性黄斑变性(Age-related macular degeneration，AMD)是视网膜黄斑区结构的衰老性改变。AMD易发于50岁以上的人群，其患病率随着年龄的增长而升高，可能导致不可逆的视力下降甚至致盲。AMD主要临床观察表现为黄斑区域形成黄色玻璃膜疣或者视网膜黄斑区形成脉络膜新生血管。临床上年龄相关性黄斑变性的诊断方法主要基于眼底血流造影成像。
皮肤的鲜红斑痣是由无数扩张的毛细血管所组成的扁平的斑块，是一种先天性毛细血管畸形。病灶面积随身体生长而相应增大，终生不消退。临床上采用光动力疗法治疗鲜红斑痣。皮肤的血流造影成像可以研究治疗前后的血管网络分布情况，评估治疗的效果。
目前有许多方法用于血流造影，包括磁共振成像(magnetic resonance imaging，MRI)，电子计算机断层扫描(computed tomography，CT)，正电子发射计算机断层扫描(positron emission tomography，PET)。但是传统的MRI、CT和PET的分辨率大于500μm，难以对～100μm直径的血管进行成像。超声成像具有更好的空间分辨率，达到约100μm，并且能够检测大约1cm深度的样品。
当前常规超声成像技术，可以对体内较大血管内的血流进行成像，但是超声分辨率也难以满足直径小于100μm的小血管成像需求。
光学成像技术通常具有更高的空间分辨率，包括光学相干层析成像(opticalcoherence tomography,OCT)、光声成像(photoacoustic imaging，PAT)、荧光显微成像(fluorescence microscopy，FM)、激光散斑成像(laser speckle imaging，LSI)等。
激光散斑成像是一种广场成像技术，散射波动容易受到周围组织的干扰，因此，成像的分辨率和图像信噪比较差。
荧光显微成像提供了高的空间分辨率，通过检测荧光染料、荧光蛋白或自发荧光，可以用于脑血管网络成像，但是生物组织对荧光散射较强，荧光成像的深度很浅，局限在组织表面附近。
光声成像利用激光激发组织，通过超声换能器检测生物组织受激发射的超声信号。这种技术具有光学的对比度、超声的分辨率和穿透深度，但分辨率仍然不能满足微小血管网络的成像需求。
OCT成像是一种非侵入的高分辨三维医学成像技术。与超声不同，OCT利用短相干长度的可见光或近红外光，基于光的干涉原理，对光学散射介质(如生物组织)进行扫描成像，空间分辨率达到约10μm，生物组织的成像深度为2-3mm，适合视网膜血管、脑皮层微血管、皮肤血管网络的在体成像。OCT系统获取的基于组织散射差异的数据，结合不同的数据分析方法，提取血管中运动红细胞的信息，可以获取高分辨率的血管网络分布图。
但是，目前OCT成像对信号的处理过程很复杂，无法实现快速成像。
基于上述问题，本发明提供了一种活体组织的光学相干层析成像系统及方法，可是实现快速对活体组织进行成像的目标。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
参考图1，图1为本发明实施例提供的一种活体组织的光学相干层析成像方法的流程示意图。
所述光学相干层析成像方法包括：
S101：在所述活体组织的同一表面位置进行M次采样，M≥2；
S102：获取M次采样过程中，随波长变化的干涉信号序列；
S103：依据所述干涉信号序列，获得随深度变化的复数信号序列；
S104：依据所述复数信号序列，获取此表面位置不同深度的血流速度参数，所述血流速度参数的值在0-1之间；
S105：扫描不同表面位置后，依据所述血流速度参数，显示不同位数的血流造影图像。
在该实施例中，通过在所述活体组织的同一表面位置进行M次采样，M≥2；获得随波长变化的干涉信号序列，将干涉信号序列进行转换，变为随深度变化的复数信号序列，依据复数信号序列进行简单的运算，即可得到数值在0-1之间的血流速度参数，也就是说无需对血流速度参数进行归一化处理，进而可显示不同位数的血流造影图像。
也就是说，该光学相干层析成像方法可以简单快速的实现对活体组织的血流造影成像。
下面对信号处理的具体方式进行说明：
在OCT信号处理过程中，将采集到的某一表面位置(x，y)在时间t随波长λ变化的干涉信号序列Γx，y，t(λ)，表面位置为OCT光束垂直平面的位置，通过滤波消除噪声，然后进行快速傅立叶变换(Fast Fourier Transform)，可得到随深度z变化的复数信号，依赖成像深度z变化的复数形式的信号可以表示为
包含幅度Ax，y，z，t部分和相位/>
部分。
如图2所示，其中(x，y)表示与OCT光束垂直平面的坐标，z表示OCT光束方向的坐标。
也就是说，如图3所示，OCT结构成像时，每个(x，y)位置只采集一次信号，用Γx，y(λ)表示，快速傅立叶变换后可得到随深度z变化的复数信号Fx，y(z)。
其中A-scan可以得到沿深度方向的一维信息，B-scan可以得到二维截面图像，C-scan可以得到三维图像。
进一步的，结合以下公式，对其参数进行处理：
或，
其中，Iflow表示血流速度参数；
Fx，y，z，t表示位置(x，y，z)第t次采样的OCT复数信号；
Fx，y，z，t+1表示位置(x，y，z)第t+1次采样的OCT复数信号；
|Fx，y，z，t|表示复数Fx，y，z，t的取模运算，得到复数信号的幅度部分。
当考虑到多次采样以减小信号噪声，在本申请实施例中引入平均计算，进而可以变换为以下四种形式：
(1)幅度差分
由于
0＜||Fx，y，z，t|-|Fx，y，z，t+1||＜|Fx，y，z，t|+|Fx，y，z，t+1|
进而可得：
整理可得：
(2)复数差分
由于
0＜|Fx，y，z，t-Fx，y，z，t+1|＜|Fx，y，z，t|+|Fx，y，z，t+1|
进而可得：
整理可得：
(3)幅度差分和
由于
0＜||Fx，y，z，t|-|Fx，y，z，t+1||＜|Fx，y，z，t|+|Fx，y，z，t+1|
进而可得：
整理可得：
(4)复数差分和
由于
0＜|Fx，y，z，t-Fx，y，z，t+1|＜|Fx，y，z，t|+|Fx，y，z，t+1|
进而可得：
整理可得：
上述(1)-(4)中，Fx，y，z，t表示在同一三维空间位置(x，y，z)第t次采样的复数OCT信号。
当在(x，y，z)三维空间位置采样M次，即可获得Fx，y，z，1，Fx，y，z，2，......，Fx，y，z，M复数信号，将M个信号进行上述运算，即可获得Iflow值。
其中，Iflow的取值范围在0-1之间，因此无需后续的归一化处理。
在图像显示时，只需将0-1区间的计算值映射到0-255(8bit图像)、0-4095(12bit图像)、0-65535(16bit图像)或其它位数的图像。
进一步的，基于本发明上述实施例，在本发明另一实施例中，所述在所述活体组织的同一表面位置进行M次采样，包括：
OCT光束在所述活体组织的预设表面位置形成光斑；
在所述光斑的固定落点处，进行M次采样。
或，
所述在所述活体组织的预设区域内M个不同表面位置进行1次采样，包括：
OCT光束在所述活体组织的预设区域进行侧向位移；
在M个不同表面位置上均进行1次采样；
所述预设区域的尺寸小于等于所述OCT光束的光斑直径。
在该实施例中，为了获得同一位置在不同时刻的OCT信号，可以将OCT光束停留在同一坐标(x，y)的位置，并进行M次采样，然后再移动OCT光束，采集下一个位置，如图4所示。
也可以在同一y坐标，沿x轴往复扫描M次，如图5所示，或在同一x坐标，沿y轴往复扫描M次，提取同一坐标(x，y)的M次采样值。
进一步的，如图6所示，还可以采取密集扫描的方式，提取Fx，y，z，1，Fx+1，y，z，1，......，Fx+M-1，y，z，1数据，或Fx，y，z，1，Fx，y+1，z，1，......，Fx，y+M-1，z，1数据进行上述运算。
假设OCT光束的圆形光斑直径为R，扫描时相邻OCT光束的移动距离为D，当D×(M-1)小于或与R接近时，可以认为Fx，y，z，1，Fx+1，y，z，1，......，Fx+M-1，y，z，1是在不同时刻同一位置采集的数据，同样，Fx，y，z，1，Fx，y+1，z，1，......，Fx，y+M-1，z，1是在不同时刻同一位置采集的数据，该数据采用上述的运算方式进行运算。
当沿x轴方向进行密集扫描采样，且D×(M-1)小于或与R接近时，使用下列四种形式公式计算，其中M为平均计算的窗口。
(5)幅度差分
如图7所示，图7为利用幅度差分算法的大鼠脑皮层的OCT血流造影投影图。
在成像时，将大鼠麻醉，去除大鼠头皮，磨薄大鼠的颅骨，然后使用OCT光束分别沿X和y两个方向进行逐点扫描，得到三维空间(x，y，z)分布的OCT复数信号。
在此，采用沿x方向的密集扫描为例说明，每个位置只采样一次，使用Fx，y，z，1表示，M取值为5，因此，在计算中，认为Fx，y，z，1，Fx+1，y，z，1，Fx+2y，z，1，Fx+3，y，z，1，Fx+4，y，z，1这5个顺序采样点对应同一空间位置不同时刻的OCT复数信号。
使用血流造影数据处理流程得到三维的OCT血流参数Iflow，Iflow的数值在0-1之间，直接映射到0-255。
当Iflow等于0时，对应图像灰度为0；当Iflow等于1时，对应图像灰度为255，即可重建出0-255灰度范围的三维血流造影图像。
对三维图像沿z方向进行标准差投影后，得到如图7所示的投影图像。
(6)复数差分
如图8所示，图8为利用复数差分算法的大鼠脑皮层的OCT血流造影投影图。
(7)幅度差分和
如图9所示，图9为利用幅度差分和算法的大鼠脑皮层的OCT血流造影投影图。
(8)复数差分和
如图10所示，图10为利用复数差分和算法的大鼠脑皮层的OCT血流造影投影图。
当沿y轴方向进行密集扫描采样时，D×(M-1)且小于或与R接近时，使用下列四种形式的公式计算，其中M为平均计算的窗口。
(9)幅度差分
(10)复数差分
(11)幅度差分和
(12)复数差分和
通过上述描述可知，本发明提供的一种活体组织的光学相干层析成像方法，通过在所述活体组织的同一表面位置进行M次采样，M≥2；获得随波长变化的干涉信号序列，将干涉信号序列进行转换，变为随深度变化的复数信号序列，依据复数信号序列进行简单的运算，即可得到数值在0-1之间的血流速度参数，也就是说无需对血流速度参数进行归一化处理，进而可显示不同位数的血流造影图像。
也就是说，该光学相干层析成像方法可以简单快速的实现对活体组织的血流造影成像。
基于本发明上述实施例，在本发明另一实施例中还提供了一种活体组织的光学相干层析成像系统。参考图11，图11为本发明实施例提供的一种活体组织的光学相干层析成像系统的结构示意图。
所述光学相干层析成像系统包括：光源发射装置、光学耦合器、参考光路、样品光路和信号处理成像装置；
所述光源发射装置用于出射弱相干光；
所述光学耦合器用于将所述弱相干光分束，一部分入射至所述参考光路，另一部分入射至所述样品光路；
所述参考光路用于将光束反射回所述光学耦合器，所述样品光路用于将被样品散射后的光束输送至所述光学耦合器，二者在光学耦合器中发生干涉；
所述信号处理成像装置获取干涉信号，进行成像。
其中，所述光源发射装置包括：
扫频光源，用于不同时刻输出不同波长的弱相干光；
偏振控制器，用于对所述弱相干光的偏振态进行处理。
所述信号处理成像装置包括：
光电探测器，用于检测所述干涉信号；
上位机，用于依据所述干涉信号进行成像。
在该实施例中，基于扫频光源，扫频光源输出的光束通过光学耦合器分光后，一部分光进入样品光路，一部分光进入参考光路。参考光路的光经过透镜后，由反射镜反射。进入样品光路的光束，经过透镜后，聚焦在生物组织上。样品光路和参考光路返回的散射光和反射光，进入光学耦合器中产生干涉，由光电探测器检测。扫频光源每一时刻输出一个波长的光束，在较短的时间范围内，完成波长的改变，实现不同波长光束的输出。
基于本发明上述实施例，在本发明又一实施例中还提供了另一种活体组织的光学相干层析成像系统。参考图12，图12为本发明实施例提供的另一种活体组织的光学相干层析成像系统的结构示意图。
所述光学相干层析成像系统包括：光源发射装置、光学耦合器、参考光路、样品光路和信号处理成像装置；
所述光源发射装置用于出射弱相干光；
所述光学耦合器用于将所述弱相干光分束，一部分入射至所述参考光路，另一部分入射至所述样品光路；
所述参考光路用于将光束反射回所述光学耦合器，所述样品光路用于将被样品散射后的光束输送至所述光学耦合器，二者在光学耦合器中发生干涉；
所述信号处理成像装置获取干涉信号，进行成像。
其中，所述光源发射装置包括：
连续光谱光源，用于在同一时刻输出不同波长的光束；
光学隔离器，用于实现光束的单向通行。
其中，所述信号处理成像装置包括：
光栅，用于将不同波长的干涉光在空间上分开；
相机，用于检测所述干涉信号；
上位机，用于依据所述干涉信号进行成像。
其中，所述参考光路包括：第一透镜和反射镜；
所述样品光路包括：第二透镜和样品。
在该实施例中，基于连续宽光谱光源，光源输出的光束通过光学隔离器后，由光学耦合器分光，一部分光进入样品光路，一部分光进入参考光路。参考光路的光经过偏振控制器和透镜后，由反射镜反射。进入样品光路的光，经过偏振控制器和透镜后，聚焦在生物组织上。样品光路和参考光路返回的散射光和反射光，在光学耦合器中干涉，干涉信号通过透镜后，由光栅进行分光，将不同波长的干涉光在空间上分开，经过透镜后在相机上进行光电转换。
以上对本发明所提供的一种活体组织的光学相干层析成像系统及方法进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
需要说明的是，本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
还需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备所固有的要素，或者是还包括为这些过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个......”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (9)
1.一种活体组织的光学相干层析成像方法，其特征在于，所述光学相干层析成像方法包括：
在所述活体组织的同一表面位置进行M次采样，M≥2；
获取M次采样过程中，随波长变化的干涉信号序列；
依据所述干涉信号序列，获得随深度变化的复数信号序列，其中，所述深度变化的方向为光学相干层析成像光束的传播方向，表面位置为光学相干层析成像光束垂直平面的位置；
依据所述复数信号序列，获取此表面位置不同深度的血流速度参数，所述血流速度参数的值在0-1之间；
扫描不同表面位置后，依据所述血流速度参数，显示不同位数的血流造影图像；
其中，所述依据所述复数信号序列，获取此表面位置不同深度的血流速度参数，包括：
或，
其中，Iflow表示血流速度参数；
Fx，y，z，t表示位置(x，y，z)第t次采样的光学相干层析成像复数信号；
Fx，y，z，t+1表示位置(x，y，z)第t+1次采样的光学相干层析成像复数信号；
|Fx，y，z，t|表示复数Fx，y，z，t的取模运算，得到复数信号的幅度部分；
(x，y)表示与所述光学相干层析成像光束垂直平面的坐标，z表示所述光学相干层析成像光束方向的坐标。
2.根据权利要求1所述的光学相干层析成像方法，其特征在于，所述在所述活体组织的同一表面位置进行M次采样，包括：
光学相干层析成像光束在所述活体组织的预设表面位置形成光斑；
在所述光斑的固定落点处，进行M次采样。
3.根据权利要求1所述的光学相干层析成像方法，其特征在于，所述在所述活体组织的预设区域内M个不同表面位置进行1次采样，包括：
光学相干层析成像光束在所述活体组织的预设区域进行侧向位移；
在M个不同表面位置上均进行1次采样；
所述预设区域的尺寸小于等于所述光学相干层析成像光束的光斑直径。
4.一种活体组织的光学相干层析成像系统，其特征在于，应用于如权利要求1-3任一所述的光学相干层析成像方法，所述光学相干层析成像系统包括：光源发射装置、光学耦合器、参考光路、样品光路和信号处理成像装置；
所述光源发射装置用于出射弱相干光；
所述光学耦合器用于将所述弱相干光分束，一部分入射至所述参考光路，另一部分入射至所述样品光路；
所述参考光路用于将光束反射回所述光学耦合器，所述样品光路用于将被样品散射后的光束输送至所述光学耦合器，二者在光学耦合器中发生干涉；
所述信号处理成像装置获取干涉信号，进行成像。
5.根据权利要求4所述的光学相干层析成像系统，其特征在于，所述光源发射装置包括：
扫频光源，用于不同时刻输出不同波长的弱相干光；
偏振控制器，用于对所述弱相干光的偏振态进行处理。
6.根据权利要求5所述的光学相干层析成像系统，其特征在于，所述信号处理成像装置包括：
光电探测器，用于检测所述干涉信号；
上位机，用于依据所述干涉信号进行成像。
7.根据权利要求4所述的光学相干层析成像系统，其特征在于，所述光源发射装置包括：
连续光谱光源，用于在同一时刻输出不同波长的光束；
光学隔离器，用于实现光束的单向通过。
8.根据权利要求7所述的光学相干层析成像系统，其特征在于，所述信号处理成像装置包括：
光栅，用于将不同波长的干涉光在空间上分开；
相机，用于检测所述干涉信号；
上位机，用于依据所述干涉信号进行成像。
9.根据权利要求4所述的光学相干层析成像系统，其特征在于，所述参考光路包括：第一透镜和反射镜；
所述样品光路包括：第二透镜和样品。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010348181.2A CN111436909B (zh) | 2020-04-28 | 2020-04-28 | 一种活体组织的光学相干层析成像系统及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010348181.2A CN111436909B (zh) | 2020-04-28 | 2020-04-28 | 一种活体组织的光学相干层析成像系统及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111436909A CN111436909A (zh) | 2020-07-24 |
| CN111436909B true CN111436909B (zh) | 2023-03-28 |
Family
ID=71657676
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010348181.2A Active CN111436909B (zh) | 2020-04-28 | 2020-04-28 | 一种活体组织的光学相干层析成像系统及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111436909B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113686735B (zh) * | 2021-09-16 | 2023-11-07 | 北京信息科技大学 | 测量血液凝固属性的方法及装置 |
| US20230301510A1 (en) * | 2022-03-24 | 2023-09-28 | Worcester Polytechnic Institute | Robotic-assisted optical coherence tomography (oct) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009165710A (ja) * | 2008-01-17 | 2009-07-30 | Univ Of Tsukuba | 眼底血流量の定量測定装置 |
| CN104523233A (zh) * | 2014-12-29 | 2015-04-22 | 浙江大学 | 基于复数互相关的微血管光学造影及抖动补偿方法与系统 |
| CN105310676A (zh) * | 2015-11-18 | 2016-02-10 | 河北大学 | 基于光学相干原理的成像技术中的位相恢复的方法 |
| CN109907731A (zh) * | 2019-01-31 | 2019-06-21 | 浙江大学 | 基于特征空间的光学相干层析的三维血流造影方法及系统 |
| CN110881947A (zh) * | 2019-12-06 | 2020-03-17 | 北京信息科技大学 | 一种光学相干断层成像方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8433393B2 (en) * | 2011-07-07 | 2013-04-30 | Carl Zeiss Meditec, Inc. | Inter-frame complex OCT data analysis techniques |
| US10368734B2 (en) * | 2015-02-19 | 2019-08-06 | Carl Zeiss Meditec, Inc. | Methods and systems for combined morphological and angiographic analyses of retinal features |
| CN105342568B (zh) * | 2015-11-26 | 2018-02-02 | 上海交通大学 | 联合相位和幅值的光学相干造影方法及系统 |
| CN106137134B (zh) * | 2016-08-08 | 2023-05-12 | 浙江大学 | 多角度复合的血流成像方法及系统 |
| WO2019010507A1 (en) * | 2017-07-14 | 2019-01-17 | Wavesense Engineering Gmbh | OPTICAL APPARATUS |
| CN107452029A (zh) * | 2017-07-31 | 2017-12-08 | 中国医学科学院生物医学工程研究所 | 一种光学微血管血流成像方法 |
| CN110292361B (zh) * | 2019-07-30 | 2024-04-30 | 佛山科学技术学院 | 一种超宽范围的皮肤成像设备与方法 |
-
2020
- 2020-04-28 CN CN202010348181.2A patent/CN111436909B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009165710A (ja) * | 2008-01-17 | 2009-07-30 | Univ Of Tsukuba | 眼底血流量の定量測定装置 |
| CN104523233A (zh) * | 2014-12-29 | 2015-04-22 | 浙江大学 | 基于复数互相关的微血管光学造影及抖动补偿方法与系统 |
| CN105310676A (zh) * | 2015-11-18 | 2016-02-10 | 河北大学 | 基于光学相干原理的成像技术中的位相恢复的方法 |
| CN109907731A (zh) * | 2019-01-31 | 2019-06-21 | 浙江大学 | 基于特征空间的光学相干层析的三维血流造影方法及系统 |
| CN110881947A (zh) * | 2019-12-06 | 2020-03-17 | 北京信息科技大学 | 一种光学相干断层成像方法 |
Non-Patent Citations (2)
| Title |
|---|
| 李培.微血管功能网络的三维光学成像技术及脑科学应用研究.《中国优秀硕士学位论文全文数据库 基础科学辑》.2019,全文. * |
| 杨强等.正向切片光学相干断层成像系统的信号提取方法.《光学学报》.2019,第0711001-1至第0711001-17页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111436909A (zh) | 2020-07-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210321874A1 (en) | Transcranial photoacoustic/thermoacoustic tomography brain imaging informed by adjunct image data | |
| Deegan et al. | Microvascular imaging of the skin | |
| Taruttis et al. | Optoacoustic imaging of human vasculature: feasibility by using a handheld probe | |
| US11471057B2 (en) | Single-impulse panoramic photoacoustic computed tomography (SIP-PACT) | |
| Mahmud et al. | Review of speckle and phase variance optical coherence tomography to visualize microvascular networks | |
| JP6148257B2 (ja) | レーザー光音響超音波撮像システム（Ｌｏｕｉｓ）および使用方法 | |
| Hu et al. | Photoacoustic imaging and characterization of the microvasculature | |
| JP6200902B2 (ja) | 生体内の光学的流れイメージング | |
| CN106943124A (zh) | 用光学相干断层扫描血管造影对局部循环进行量化 | |
| CN104323762B (zh) | 一种基于光声显微成像的鲜红斑痣血管定量化检测装置 | |
| JP2013518695A (ja) | 超高感度光学微小血管造影用の方法および装置 | |
| CN107862724B (zh) | 一种改进的微血管血流成像方法 | |
| CN111436909B (zh) | 一种活体组织的光学相干层析成像系统及方法 | |
| Dolganova et al. | Nanoparticle-enabled experimentally trained wavelet-domain denoising method for optical coherence tomography | |
| Ma et al. | Multiscale confocal photoacoustic dermoscopy to evaluate skin health | |
| Sdobnov et al. | Advances in dynamic light scattering imaging of blood flow | |
| CN107137073A (zh) | 一种基于甲壁血流分析的青光眼和糖尿病等微循环系统疾病检测技术 | |
| US10909683B2 (en) | Methods and systems for enhancing optical image quality | |
| CN114209278A (zh) | 一种基于光学相干层析成像的深度学习皮肤病诊断系统 | |
| Choi et al. | Volumetric cutaneous microangiography of human skin in vivo by VCSEL swept-source optical coherence tomography | |
| Jaeger et al. | Real-time clinical clutter reduction in combined epi-optoacoustic and ultrasound imaging | |
| WO2018230409A1 (ja) | 情報処理装置、情報処理方法、及びプログラム | |
| Morales-Vargas et al. | Exposure Time and Depth Effect in Laser Speckle Contrast Images under an Adaptive Processing | |
| Naumovska et al. | Mapping the architecture of the temporal artery with photoacoustic imaging for diagnosing giant cell arteritis | |
| Lal | Advancements in structural and functional imaging using Fourier domain optical coherence tomography |
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 |