CN106974622A - 基于光学相干断层成像的斑块稳定性测量方法及系统 - Google Patents
基于光学相干断层成像的斑块稳定性测量方法及系统 Download PDFInfo
- Publication number
- CN106974622A CN106974622A CN201710219596.8A CN201710219596A CN106974622A CN 106974622 A CN106974622 A CN 106974622A CN 201710219596 A CN201710219596 A CN 201710219596A CN 106974622 A CN106974622 A CN 106974622A
- Authority
- CN
- China
- Prior art keywords
- patch
- image
- module
- coherence tomography
- optical coherence
- 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
Links
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/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
-
- 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
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)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Vascular Medicine (AREA)
- Cardiology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明提供了一种基于光学相干断层成像的斑块稳定性测量方法及系统，该系统包括：影像采集模块、影像接收模块、图像处理模块、结果显示模块，其中：影像采集模块用于冠脉影像采集和图像信号的产生；影像接收模块用于接收影像采集模块产生的图像信号并传输给图像处理模块；图像处理模块用于对接收到的图像信号进行图像处理和分析，包含纤维帽信息获取模块、巨噬细胞信息获取模块和脂质核信息获取几个子模块；结果显示模块用于显示图像处理模块的测量分析结果。本发明突破了传统光学相干断层成像无法获取脂质核大小参数的局限，在保证光学相干断层成像特有高分辨率结果的同时，添加了脂质核大小和斑块应力的信息。
Description
技术领域
本发明应用于医疗器械及信号处理领域，尤其涉及一种基于光学相干断层成像的图像处理和分析系统，以及应用在基于光学相干断层成像影像的对斑块稳定性进行快速高效准确的测量系统及方法。
背景技术
光学相干断层成像(Optical Coherence Tomography，OCT)是一种新兴的高分辨率断层成像技术。和血管内超声(Intravascular Ultrasound，IVUS)原理类似，该技术以近红外线作为光源，利用光波的干涉法则进行成像，将光束扫描组织的光信号转换成电信号，经过计算机处理后显示为灰色图或伪彩色图的二维和三维图像。由于光的波长较声波更短，光学相干断层成像具有更高的分辨率，可达10～20um。和血管内超声相比，光学相干断层成像能够提供更为快速、清晰的图像，并且更精确地显示血管壁的具体形态情况(内部是否堵塞、堵塞面积、大小等)。
在对纤维帽厚度检测方面，OCT的优势非常明显。2006年，日本学者Takashi Kubo等人(KUBO T,IMANISHI T,TAKARADA S,et al.Assessment of culprit lesionmorphology in acute myocardial infarction-Ability of optical coherencetomography compared with intravascular ultrasound and coronary angioscopy[J].Journal of the American College of Cardiology,2007,50(10):933-9.)通过对30例实验分别进行OCT、IVUS和CAS(Coronary Angioscopy)分析。结果表明，在对薄纤维帽斑块的识别上，光学相干断层成像显著优于其他影像技术，具有更高的准确性。
斑块内大量巨噬细胞浸润是评判斑块稳定性的另外一个重要标准。早在2004年，就已经有研究表明，破裂斑块部位的巨噬细胞密度大于非破裂斑块(MACNEILL B D,JANG IK,BOUMA B E,et al.Focal and multi-focal plaque distributions in patients withmacrophage acute and stable presentations of coronary artery disease[J].Journal of the American College of Cardiology,2004,44(5):972-9.)。
同样，在巨噬细胞浸润方面，光学相干断层成像也具有优势。早在2003年，BrettE.Bouma等(Tearney GJ,Yabushita H,Houser SL,Aretz HT,Jang IK,Schlendorf KH,Kauffman CR,Shishkov M,Halpern EF,Bouma BE.Quantification of macrophagecontent in atherosclerotic plaques by optical coherencetomography.Circulation 2003；107:113–9.)报道了光学相干断层成像能够在体定性、定量地判断巨噬细胞的存在，并且提供了基于光学相干断层成像的判断方法。
斑块负荷指的是斑块在血管上所在部位的面积或体积，对于脂质斑块而言，脂质核的大小即为其斑块负荷。在脂质核识别和相关参数获得方面，光学相干断层成像存在一定的技术限制：其组织穿透能力较低，仅限于1～3mm，在某些情况下，并不能获得斑块处脂质核心所在位置，也难以得到脂质核大小数据。因此并不能通过传统光学相干断层成像获得全面的信息。对于近年来较为新兴的近红外光谱成像技术(Near-InfraredSpectroscopy,NIRS)，根据Brugaletta团队在2011年对其的研究结果，虽然NIRS在脂质核心识别方面具有较大优势(BRUGALETTA S,GARCIA-GARCIA H M,SERRUYS P W,et al.NIRSand IVUS for Characterization of Atherosclerosis in Patients UndergoingCoronary Angiography[J].JACC-Cardiovasc Imag,2011,4(6):647-55.)，但是它并不能定量测量脂质核大小。相比较而言，血管内超声虽然成像由于其较优的组织穿透性(10mm左右)，通常被用于进行斑块脂质核大小的检测，但是由于其分辨率超过100um，故无法识别薄纤维帽。
现有的较优的对斑块稳定性的测量及分析方法，是结合光学相干断层成像对纤维帽厚度和巨噬细胞浸润情况的测量分析和血管内超声对脂质核大小的测量分析，从而对斑块稳定性进行全面的分析。这种方法是在2014年提出的，Li和他的研究团队(LI X,LI J W,JING J,et al.Integrated IVUS-OCT Imaging for Atherosclerotic PlaqueCharacterization[J].IEEE J Sel Top Quantum Electron,2014,20(2):8.)开发出了一种结合IVUS和OCT的成像系统，可以同时进行血管内超声以及光学相干断层成像。
但是，上述设备和方法也存在几点局限：1)导丝的硬质前端长度较长，可能导致导丝无法通过，并且容易造成图像中的非均匀旋转失真(Nonuniform RotationalDistortion)；2)上述系统的成像帧率有所局限，这主要是因为现有的超声脉冲发生器以及扫描激光光源重复频率存在限制；3)两种成像系统的结合势必伴随着成本的提升，不利于该种方法的普及。
此外，目前尚无仅通过光学相干断层成像系统获得斑块脂质核大小和斑块应力信息来分析斑块稳定性的方法和系统。
发明内容
有鉴于此，本发明实施例提供一种基于光学相干断层成像的斑块稳定性测量方法及系统，通过调整传统光学相干断层成像步骤，在获得纤维帽厚度和巨噬细胞浸润情况等基本参数的同时，记录斑块所在位置的血管段在心动周期中的变形行为，结合该变形行为与脂质核大小、斑块应力的对应关系，能够另外获得脂质核大小、斑块应力相关数据，突破了传统过程中必须依赖血管内超声成像或近红外光谱成像等其他腔内成像方式以了解脂质核信息的局限。最后，结合上述信息，实现对斑块稳定性进行高效全面且成本较低的测量分析。
具体而言，本发明提供了以下的技术方案：
一方面，本发明提供的一种基于光学相干断层成像的斑块稳定性测量系统，包括：
影像采集模块，用于冠脉影像采集，并产生图像信号；
影像接收模块，用于接收影像采集模块产生的图像信号并传输给图像处理模块；
图像处理模块，对接收到的图像信号进行处理与分析。其中包括对斑块进行可见特征测量分析和潜在特征测量分析，其中斑块可见特征测量分析包括纤维帽厚度获取和巨噬细胞密度获取，斑块潜在特征测量分析通过成像过程中管腔变形数据获取脂质核大小和斑块应力；
结果显示模块，用于显示图像处理模块的测量分析结果。
优选地，图像处理模块通过一种或多种方法进行可见特征测量分析和潜在特征测量分析。其中包括以下子模块：纤维帽信息获取模块，处理斑块处纤维帽所在区域的图像，实现对纤维帽厚度的测量；巨噬细胞信息获取模块，处理斑块处巨噬细胞所在区域的图像，实现对巨噬细胞的识别和密度计算；脂质核信息获取模块，处理斑块处血管管腔所在区域的图像，实现对脂质核大小及斑块应力计算。
优选地，所述系统中的图像处理模块和结果显示模块集成在一起，实现图像自动处理、情况分析和结果显示等功能。
优选地，该斑块稳定性测量系统在影像采集的过程中实时进行数据传输和图像处理及分析，以实现斑块稳定性测量分析的效率最大化。
另一方面，本发明提供的一种基于光学相干断层成像的斑块稳定性测量方法，具体而言，包括：
步骤1、确定感兴趣血管的起始位置、终止位置；
步骤2、确定斑块在感兴趣血管中的所在位置；
步骤3、进行光学相干断层成像过程中，将扫描装置停留在步骤2确定的位置，扫描多个心动周期，以获得斑块所在位置在心动周期中多时刻的断层成像的图像；
步骤4、扫描装置在斑块所在位置停留多个心动周期后，继续进行传统光学相干断层成像过程；
步骤5、获得步骤3时间段内的多时刻的断层成像的图像，在心动周期中，设置多个关键时间点，得到在心动周期各个关键时间点斑块处血管管腔的轮廓形态变化；
步骤6、利用步骤5中获得的血管管腔变形行为数据，获得脂质核大小和斑块应力数据；
步骤7、利用步骤4中获得的斑块图像，通过光学相干断层成像对斑块进行分析，所述分析包括纤维帽厚度和巨噬细胞浸润程度；
步骤8、结合步骤6获得的结果以及步骤7获得的结果，实现对斑块稳定性的全面准确快速测量分析。
优选地，所述步骤5进一步包括：
步骤501、心动周期中的多个关键时间点的选取通过结合心电图各个波段间的分界点进行选择；
步骤502、分析步骤501中选取的各个关键点，获得斑块所在位置血管管壁在该心动周期中的各个关键时间点的变形行为数据，该变形行为数据包括管腔面积和/或周长在心动周期中随时间的变化关系，管腔轮廓的线应变等指标；
步骤503、对步骤502中获得的多组血管管壁变形行为数据进行统计优化，包括但不仅限于去极值、取平均、排除异常数据等；
优选地，所述步骤6进一步包括：将步骤5中获得的心动周期内斑块所在位置血管管腔变形行为数据通过查表法与脂质核情况对应，获得脂质核大小和斑块应力数据。
优选地，所述步骤7进一步包括：
步骤701、精确测量斑块纤维帽厚度，通过斑块所在位置的纤维帽最小厚度和平均厚度来体现纤维帽的厚薄程度；
步骤702、通过发现信号丰富、明显或融合的亮点区域识别巨噬细胞浸润，通过测定光学相干断层成像信号衰减和粒度指数等参数，获得单位面积巨噬细胞数量，通过巨噬细胞密度数值来体现巨噬细胞的浸润程度；
优选地，通过扫描装置停留探测或扫描接收信号，获取感兴趣区域在一段连续时间段内的变形行为数据，并通过变形行为推导出脂质核大小等参数以帮助感兴趣区域的测量分析，上述的连续时间段可以指一个或多个心动周期；
优选地，通过斑块所在位置血管管腔变形行为数据得出斑块脂质核大小和/或斑块应力信息，其中体现血管管腔变形行为的参数包括但不仅限于血管管腔面积、血管壁周长、管腔轮廓的线应变等中的一种或多种组合；
优选地，获取斑块所在位置血管管腔在心动周期中的变形行为数据，以此得出斑块脂质核大小以及斑块应力数据。其中斑块所在位置血管管腔在心动周期中的变形行为可以通过多种途径获得，优选光学相干断层成像、血管内超声成像、近红外光谱成像等类似其他成像方式；
优选地，通过斑块所在位置的纤维帽最小厚度和平均厚度来体现纤维帽厚薄程度、通过斑块所在位置的巨噬细胞密度来体现巨噬细胞浸润程度、通过脂质核面积和/或体积来体现脂质核大小程度，并且结合以上参数以及斑块应力数据，对斑块稳定性进行测量分析；
与现有技术相比，本发明技术方案的有益效果在于：仅通过调整传统光学相干断层成像步骤，通过斑块所在血管管壁在心动周期中的变形行为，就能获得脂质核大小和斑块应力信息，减少了再进行近红外光谱成像或者其他腔内成像以识别脂质核心，和/或再进行血管内超声或者其他腔内成像以获得脂质核大小的必要性。本发明不仅能够实现更全面高效地对斑块稳定性进行测量分析，而且在提高准确性的同时降低了成本，节约了时间，为斑块稳定性的精准化测量分析提供了全新的、高效的方法。
附图说明
图1是本发明所述基于光学相干断层成像的斑块稳定性分析方法流程图；
图2是本发明所述基于光学相干断层成像的斑块稳定性分析系统框架图。
具体实施方式
下面结合附图对本发明实施例中的技术方案进行详细描述。本领域技术人员应当知晓，下述具体实施例或具体实施方式，是本发明为进一步解释具体的发明内容而列举的一系列优化的设置方式，而该些设置方式之间均是可以相互结合或者相互关联使用的，除非在本发明明确提出了其中某些或某一具体实施例或实施方式无法与其他的实施例或实施方式进行关联设置或共同使用。同时，下述的具体实施例或实施方式仅作为最优化的设置方式，而不作为限定本发明的保护范围的理解。
实施例1：
为进一步理解本发明，下面结合附图及实施例对本发明做进一步的阐述和说明。图1为本发明的方法流程图，在一具体实施方式中，该方法可概括为：光学探头送至病变处扫描若干个心动周期，获得同一部位在心动周期多个关键时间点的血管管腔变形情况；优选通过查表法获得斑块部位脂质核大小、斑块应力信息；结合传统光学相干断层成像对斑块纤维帽厚度、巨噬细胞浸润程度等情况的测量分析，对斑块稳定性进行较为完善全面的分析。
表1为上述查表法所使用的表格示例，通过获得的血管管腔面积最大差值，在表格中对应得到脂质核横截面大小；所述表格通过多次重复试验后，获得的统计数据。
表1
| 血管管腔面积最大差值/mm2 | 血管管腔面积最大差值/mm2 |
| 0.10 | 0.24 |
| 0.20 | 0.56 |
| 0.30 | 0.89 |
| 0.40 | 1.21 |
| 0.50 | 1.54 |
| 0.60 | 1.86 |
| 0.70 | 2.19 |
| 0.80 | 2.24 |
| 0.90 | 2.84 |
| 1.00 | 3.16 |
| … | … |
具体而言，在一具体的实施方式中，该方法可以设置为：
步骤1、确定感兴趣血管的起始位置、终止位置，所述感兴趣血管为待测量分析的目标血管；
步骤2、确定斑块在感兴趣血管中的所在位置；
步骤3、进行光学相干断层成像过程中，将扫描装置停留在步骤2确定的位置，扫描多个心动周期，以获得斑块所在位置在心动周期中多时刻的断层成像的图像；
步骤4、扫描装置在斑块所在位置停留多个心动周期后，继续进行传统光学相干断层成像过程；
步骤5、获得步骤3时间段内的多时刻的断层成像的图像，在心动周期中，设置N个关键时间点，得到在心动周期在这N个关键时间点斑块处血管管腔的轮廓形态变化。
其中，N要求大于或等于2，并需要将舒张末期和收缩末期这两个时刻包括在内；优选的，N＝5；
步骤6、将步骤5中获得的心动周期内斑块所在位置血管管腔变形行为数据，获得脂质核大小和斑块应力数据；
步骤7、利用步骤4中获得的斑块图像，通过光学相干断层成像对斑块进行分析，分析重点包括纤维帽厚度和巨噬细胞浸润程度；
步骤8、结合步骤6获得的结果以及步骤7获得的结果，实现对斑块稳定性的全面准确快速测量分析。
在一具体的实施方式中，所述步骤5进一步包括：
步骤501、心动周期中的多个关键时间点的选取通过结合心电图各个波段间的分界点进行选择；
步骤502、分析步骤501中选取的各个关键点，获得斑块所在位置血管管壁在该心动周期中的各个关键时间点的变形行为数据，该变形行为数据包括管腔面积和/或周长在所选择的心动周期各个关键点的数值大小；
步骤503、对步骤502中获得的多组血管管壁变形行为数据进行统计优化，包括但不仅限于去极值、取平均、排除异常数据等；上述去极值、取平均、排除异常数据等方法，可以采用数据处理领域中的常规算法实现，该些具体算法不作为本申请的保护范围限定理解。
在一具体的实施方式中，所述步骤6进一步包括：将该血管管腔变形行为数据通过查表法与脂质核情况对应，获得脂质核大小和斑块应力数据；
在一具体的实施方式中，所述步骤7进一步包括：
步骤701、精确测量斑块纤维帽厚度，通过斑块所在位置的纤维帽最小厚度和平均厚度来体现纤维帽的厚薄程度；
步骤702、通过发现信号丰富、明显或融合的亮点区域识别巨噬细胞浸润，通过测定光学相干断层成像信号衰减和粒度指数等参数，获得单位面积巨噬细胞数量，通过巨噬细胞密度数值来体现巨噬细胞的浸润程度；
在一具体的实施方式中，可以选择血管管腔面积、血管壁周长、管腔轮廓的线应变等参数中的一个或多个，通过计算上述参数在心动周期各个关键点间的最大差值来体现血管管腔变形程度，差值越大，表示血管变形程度越大。由此得出斑块脂质核大小和斑块应力信息；
在一具体的实施方式中，斑块所在位置血管管腔在心动周期中的变形行为可以通过光学相干断层成像、血管内超声成像、近红外光谱成像等多种方式获得；
在一具体的实施方式中，通过斑块所在位置的纤维帽最小厚度和平均厚度来体现纤维帽厚薄程度、通过斑块所在位置的巨噬细胞密度来体现巨噬细胞浸润程度、通过脂质核面积和/或体积来体现脂质核大小程度。利用以上数据，通过公式计算得到斑块稳定性参数r，并参考该参数对斑块稳定性进行分析。
优选的，本发明提供了一种斑块稳定性参数r的计算公式，通过为纤维帽厚度、巨噬细胞密度以及脂质核面积分配不同的权重并求和，得到参数r数值。上述计算公式为：
其中：d表示斑块所在位置纤维帽最小厚度；ρ表示斑块所在位置巨噬细胞密度；S表示脂质核面积；ω1、ω2、ω3分别表示以上参数的加权系数，其中ω1+ω2+ω3＝1。
优选的，ω1＝0.40-0.58、ω2＝0.23-0.47、ω3＝0.15-0.39。
在一具体的实施方式中，获取的所述斑块稳定性测量结果，后续可以结合先前具有类似结果的统计数据，评估该斑块的潜在风险。
实施例2：
在另一个具体的实施例中，利用本发明所述方法对特定血管上的斑块稳定性进行分析。具体而言，该方法可以通过以下步骤来实现：
(1)、选取心脏左冠回旋支为感兴趣血管；
(2)、确定斑块在左冠回旋支中所处的位置；
(3)、进行光学相干断层成像过程中，将扫描装置停留在步骤(2)确定的位置，扫描多个心动周期，以获得斑块所在位置在心动周期中多时刻的断层成像的图像；
(4)、扫描装置在步骤(2)确定的斑块所在位置停留多个心动周期后，继续进行传统光学相干断层成像过程；
(5)、获得步骤(3)时间段内的多时刻的断层成像的图像，在心动周期中，选取N个时刻，得到在这N个时刻斑块处血管管腔的轮廓形态变化。其中，N要求大于或等于2，并需要将舒张末期和收缩末期这次两个时刻包括在内；优选的，N＝5。
(6)、将步骤(5)中获得的心动周期内斑块所在位置血管管腔变形行为数据通过查表法与脂质核情况对应，获得脂质核大小和斑块应力数据；
(7)、利用步骤(4)中获得的心脏左冠回旋支上的斑块图像，通过光学相干断层成像对斑块进行分析，分析重点包括纤维帽厚度和巨噬细胞浸润程度；
(8)、结合步骤(6)获得的结果以及步骤(7)获得的结果，实现对斑块稳定性的全面准确快速测量分析。
实施例3：
在又一具体的实施例中，本发明还提供了一种基于光学相干断层成像的斑块稳定性测量系统，所述系统包括：
影像采集模块，用于冠脉影像采集，并产生图像信号；
影像接收模块，用于接收影像采集模块产生的图像信号并传输给图像处理模块；
图像处理模块，对接收到的图像信号进行处理与分析。其中包括对斑块进行可见特征测量分析和潜在特征测量分析，其中斑块可见特征测量分析包括纤维帽厚度获取和巨噬细胞密度获取，斑块潜在特征测量分析通过成像过程中管腔变形数据获取脂质核大小和斑块应力；
结果显示模块，用于显示图像处理模块的测量分析结果。
优选地，图像处理模块通过一种或多种方法进行可见特征测量分析和潜在特征测量分析。其中包括以下子模块：纤维帽信息获取模块，处理斑块处纤维帽所在区域的图像，实现对纤维帽厚度的测量；巨噬细胞信息获取模块，处理斑块处巨噬细胞所在区域的图像，实现对巨噬细胞的识别和密度计算；脂质核信息获取模块，处理斑块处血管管腔所在区域的图像，实现对脂质核大小及斑块应力计算。
优选的，所述系统还可包括参数输入模块，用于接受影像接收模块的采样频率、像素大小等数据，并将数据传输给图像处理模块，以适应各种不同类型、品牌和型号的影像采集仪器，使斑块稳定性的测量分析更为迅速与便捷。
优选地，所述系统还可包括分析预测模块，利用所述斑块稳定性测量分析结果，预测该部位发生斑块破裂的风险。
优选地，图像处理模块可以在分析方法上作进一步优化，包括但不仅限于纤维帽自动识别和厚度测量、基于光学相干断层成像所得图像的斑块处巨噬细胞自动识别和密度计算、管腔面积/周长参数和脂质核大小和斑块应力自动对应等。
优选地，可以将所述系统中的图像处理模块和结果显示模块结合起来，实现图像自动处理、情况分析和结果显示等功能。
优选地，该斑块稳定性测量系统在影像采集的过程中实时进行数据传输和图像处理及分析，以实现斑块稳定性测量分析的效率最大化。
与现有技术相比，本发明技术方案的有益效果在于：仅通过调整传统光学相干断层成像步骤，通过斑块所在血管管壁在心动周期中的变形行为，就能获得脂质核大小和斑块应力信息，减少了再进行近红外光谱成像或者其他腔内成像以识别脂质核心，和/或再进行血管内超声或者其他腔内成像以获得脂质核大小的必要性。本发明不仅能够实现更全面高效地对斑块稳定性进行测量分析，而且在提高准确性的同时降低了成本，节约了时间，为斑块稳定性的精准化测量分析提供了新方法。
该系统所采用到的各计算方法或方式，可采用上述实施例1至3中的任意方式实现。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random AccessMemory，RAM)等。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。
Claims (10)
1.一种基于光学相干断层成像的斑块稳定性测量系统，其特征在于，所述系统包括：
影像采集模块，用于冠脉影像采集，并产生图像信号；
影像接收模块，用于接收影像采集模块产生的图像信号并传输给图像处理模块；
图像处理模块，对接收到的图像信号进行处理与分析。其中包括对斑块进行可见特征测量分析和潜在特征测量分析，其中斑块可见特征测量分析包括纤维帽厚度获取和巨噬细胞密度获取，斑块潜在特征测量分析通过成像过程中管腔变形数据获取脂质核大小和斑块应力；
结果显示模块，用于图像处理模块的测量分析结果。
2.根据权利要求1所述的系统，其特征在于，所述图像处理模块进一步包括以下子模块：
纤维帽信息获取模块，处理斑块处纤维帽所在区域的图像，实现对纤维帽厚度的测量；
巨噬细胞信息获取模块，处理斑块处巨噬细胞所在区域的图像，实现对巨噬细胞的识别和密度计算；
脂质核信息获取模块，处理斑块处血管管腔所在区域的图像，实现对脂质核大小及斑块应力计算。
3.一种基于光学相干断层成像的斑块稳定性测量方法，其特征在于，该方法包括：
步骤1、确定感兴趣血管的起始位置、终止位置；
步骤2、确定斑块在感兴趣血管中的所在位置；
步骤3、进行光学相干断层成像过程中，将扫描装置停留在步骤2确定的位置，扫描多个心动周期，以获得斑块所在位置在心动周期中多时刻的断层成像的图像；
步骤4、扫描装置在斑块所在位置停留多个心动周期后，继续进行传统光学相干断层成像过程；
步骤5、获得步骤3时间段内的多时刻的断层成像的图像，在心动周期中，设置多个关键时间点，得到在心动周期各个关键时间点斑块处血管管腔的轮廓形态变化数据(即血管管腔变形行为数据)；
步骤6、利用步骤5中获得的血管管腔变形行为数据，获得脂质核大小和斑块应力数据；
步骤7、利用步骤4中获得的斑块图像，通过光学相干断层成像对斑块进行分析，分析重点包括纤维帽厚度和巨噬细胞浸润程度；
步骤8、结合步骤6获得的结果以及步骤7获得的结果，实现对斑块稳定性的全面准确快速测量分析。
4.根据权利要求3所述的方法，其特征在于，所述步骤5进一步包括：
步骤501、心动周期中的多个关键时间点的选取通过结合心电图各个波段间的分界点进行选择；
步骤502、分析步骤501中选取的各个关键点，获得斑块所在位置血管管壁在该心动周期中的各个关键时间点的变形行为数据，该变形行为数据包括管腔面积和/或周长在心动周期中随时间的变化关系、管腔轮廓的线应变指标；
步骤503、对步骤502中获得的多组血管管壁变形行为数据进行统计优化，包括去极值、取平均、排除异常数据。
5.根据权利要求3所述的方法，其特征在于，所述步骤6进一步包括：
将步骤5中获得的心动周期内斑块所在位置血管管腔变形行为数据通过查表法与脂质核情况对应，获得脂质核大小和斑块应力数据。
6.根据权利要求3所述的方法，其特征在于，所述步骤7进一步包括：
步骤701、精确测量斑块纤维帽厚度，通过斑块所在位置的纤维帽最小厚度和平均厚度来体现纤维帽的厚薄程度；
步骤702、通过发现信号丰富、明显或融合的亮点区域识别巨噬细胞浸润，通过测定光学相干断层成像信号衰减和粒度指数参数，获得单位面积巨噬细胞数量，通过巨噬细胞密度数值来体现巨噬细胞的浸润程度。
7.根据权利要求3所述的方法，其特征在于，通过扫描装置停留探测或扫描接收信号，获取感兴趣区域在一段连续时间段内的变形行为数据，并通过变形行为推导出脂质核参数以帮助感兴趣区域的测量分析，所述的连续时间段为一个或多个心动周期。
8.根据权利要求3所述的方法，其特征在于，通过斑块所在位置血管管腔变形行为数据得出斑块脂质核大小和/或斑块应力信息，其中体现血管管腔变形行为的参数包括血管管腔面积、血管壁周长、管腔轮廓的线应变。
9.根据权利要求3所述的方法，其特征在于，获取斑块所在位置血管管腔在心动周期中的变形行为数据，以此得出斑块脂质核大小以及斑块应力数据；其中斑块所在位置血管管腔在心动周期中的变形行为可以通过多种途径获得，优选光学相干断层成像、血管内超声成像、近红外光谱成像等类似其他成像方式。
10.根据权利要求3或5所述的方法，其特征在于，通过斑块所在位置的纤维帽最小厚度和平均厚度来体现纤维帽厚薄程度、通过斑块所在位置的巨噬细胞密度来体现巨噬细胞浸润程度、通过脂质核面积和/或体积来体现脂质核大小程度，并且结合以上参数以及斑块应力数据，对斑块稳定性进行测量分析。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710219596.8A CN106974622B (zh) | 2017-04-06 | 2017-04-06 | 基于光学相干断层成像的斑块稳定性测量系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710219596.8A CN106974622B (zh) | 2017-04-06 | 2017-04-06 | 基于光学相干断层成像的斑块稳定性测量系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106974622A true CN106974622A (zh) | 2017-07-25 |
| CN106974622B CN106974622B (zh) | 2021-08-24 |
Family
ID=59344961
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710219596.8A Active CN106974622B (zh) | 2017-04-06 | 2017-04-06 | 基于光学相干断层成像的斑块稳定性测量系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106974622B (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109064442A (zh) * | 2018-06-21 | 2018-12-21 | 上海遥思企业管理中心 | 血管压力差修正方法、装置及设备 |
| CN109091167A (zh) * | 2018-06-29 | 2018-12-28 | 东南大学 | 冠状动脉粥样硬化斑块增长的预测方法 |
| WO2019109607A1 (zh) | 2017-12-07 | 2019-06-13 | 上海交通大学 | 基于医学影像序列斑块稳定性指标的快速计算方法及系统 |
| CN110827255A (zh) * | 2019-10-31 | 2020-02-21 | 杨本强 | 一种基于冠状动脉ct图像的斑块稳定性预测方法及系统 |
| US10631718B2 (en) | 2015-08-31 | 2020-04-28 | Gentuity, Llc | Imaging system includes imaging probe and delivery devices |
| CN113096115A (zh) * | 2021-04-28 | 2021-07-09 | 博动医学影像科技（上海）有限公司 | 冠状动脉斑块状态评估方法、装置、电子设备 |
| CN113192004A (zh) * | 2021-04-01 | 2021-07-30 | 温州医科大学附属第二医院（温州医科大学附属育英儿童医院） | 一种基于冠脉多模影像裂纹形成和扩展的简化计算方法与斑块破裂风险量化的评估系统 |
| US11278206B2 (en) | 2015-04-16 | 2022-03-22 | Gentuity, Llc | Micro-optic probes for neurology |
| US11684242B2 (en) | 2017-11-28 | 2023-06-27 | Gentuity, Llc | Imaging system |
| CN113096115B (zh) * | 2021-04-28 | 2024-06-04 | 上海博动医疗科技股份有限公司 | 冠状动脉斑块状态评估方法、装置、电子设备 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102046071A (zh) * | 2008-06-02 | 2011-05-04 | 光学实验室成像公司 | 用于从光学相干断层扫描图像获得组织特性的定量方法 |
| CN103190932A (zh) * | 2013-04-22 | 2013-07-10 | 华北电力大学（保定） | 一种冠状动脉血管壁应力和应变的估算方法 |
| CN105074379A (zh) * | 2013-03-14 | 2015-11-18 | 研究发展基金会 | 用于光学相干断层扫描和双光子荧光成像的设备和方法 |
| US20160066862A1 (en) * | 2014-09-09 | 2016-03-10 | Heartflow, Inc. | Method and system for quantifying limitations in coronary artery blood flow during physical activity in patients with coronary artery disease |
| CN106056126A (zh) * | 2015-02-13 | 2016-10-26 | 西门子公司 | 医学成像中的斑块易损性评定 |
| US20170049328A1 (en) * | 2015-08-21 | 2017-02-23 | Robert R. Alfano | Raman and resonant raman detection of vulnerable plaque optical analyzer and imager |
-
2017
- 2017-04-06 CN CN201710219596.8A patent/CN106974622B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102046071A (zh) * | 2008-06-02 | 2011-05-04 | 光学实验室成像公司 | 用于从光学相干断层扫描图像获得组织特性的定量方法 |
| CN105074379A (zh) * | 2013-03-14 | 2015-11-18 | 研究发展基金会 | 用于光学相干断层扫描和双光子荧光成像的设备和方法 |
| CN103190932A (zh) * | 2013-04-22 | 2013-07-10 | 华北电力大学（保定） | 一种冠状动脉血管壁应力和应变的估算方法 |
| US20160066862A1 (en) * | 2014-09-09 | 2016-03-10 | Heartflow, Inc. | Method and system for quantifying limitations in coronary artery blood flow during physical activity in patients with coronary artery disease |
| CN106056126A (zh) * | 2015-02-13 | 2016-10-26 | 西门子公司 | 医学成像中的斑块易损性评定 |
| US20170049328A1 (en) * | 2015-08-21 | 2017-02-23 | Robert R. Alfano | Raman and resonant raman detection of vulnerable plaque optical analyzer and imager |
Non-Patent Citations (3)
| Title |
|---|
| 余寿芬等: "影响颈动脉斑块稳定性的相关因素", 《神经病学与神经康复学杂志》 * |
| 孙前进等: "易损斑块的研究进展", 《医学综述》 * |
| 王旖旎等: "基于OCT的冠状动脉粥样硬化斑块稳定性影响因素的分析", 《临床心血管病杂志》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11278206B2 (en) | 2015-04-16 | 2022-03-22 | Gentuity, Llc | Micro-optic probes for neurology |
| US11583172B2 (en) | 2015-08-31 | 2023-02-21 | Gentuity, Llc | Imaging system includes imaging probe and delivery devices |
| US10631718B2 (en) | 2015-08-31 | 2020-04-28 | Gentuity, Llc | Imaging system includes imaging probe and delivery devices |
| US11064873B2 (en) | 2015-08-31 | 2021-07-20 | Gentuity, Llc | Imaging system includes imaging probe and delivery devices |
| US11937786B2 (en) | 2015-08-31 | 2024-03-26 | Gentuity, Llc | Imaging system includes imaging probe and delivery devices |
| US11684242B2 (en) | 2017-11-28 | 2023-06-27 | Gentuity, Llc | Imaging system |
| WO2019109607A1 (zh) | 2017-12-07 | 2019-06-13 | 上海交通大学 | 基于医学影像序列斑块稳定性指标的快速计算方法及系统 |
| CN109064442B (zh) * | 2018-06-21 | 2022-07-05 | 上海遥思企业管理中心 | 血管压力差修正方法、装置及设备 |
| CN109064442A (zh) * | 2018-06-21 | 2018-12-21 | 上海遥思企业管理中心 | 血管压力差修正方法、装置及设备 |
| CN109091167A (zh) * | 2018-06-29 | 2018-12-28 | 东南大学 | 冠状动脉粥样硬化斑块增长的预测方法 |
| CN110827255A (zh) * | 2019-10-31 | 2020-02-21 | 杨本强 | 一种基于冠状动脉ct图像的斑块稳定性预测方法及系统 |
| CN113192004B (zh) * | 2021-04-01 | 2022-01-18 | 温州医科大学附属第二医院（温州医科大学附属育英儿童医院） | 一种基于冠脉多模影像裂纹形成和扩展的简化计算方法与斑块破裂风险量化的评估系统 |
| CN113192004A (zh) * | 2021-04-01 | 2021-07-30 | 温州医科大学附属第二医院（温州医科大学附属育英儿童医院） | 一种基于冠脉多模影像裂纹形成和扩展的简化计算方法与斑块破裂风险量化的评估系统 |
| CN113096115A (zh) * | 2021-04-28 | 2021-07-09 | 博动医学影像科技（上海）有限公司 | 冠状动脉斑块状态评估方法、装置、电子设备 |
| WO2022228463A1 (zh) * | 2021-04-28 | 2022-11-03 | 上海博动医疗科技股份有限公司 | 冠状动脉斑块状态评估方法、装置、电子设备 |
| CN113096115B (zh) * | 2021-04-28 | 2024-06-04 | 上海博动医疗科技股份有限公司 | 冠状动脉斑块状态评估方法、装置、电子设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106974622B (zh) | 2021-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106974622A (zh) | 基于光学相干断层成像的斑块稳定性测量方法及系统 | |
| US11793462B2 (en) | Intravascular measurement and data collection systems, apparatus and methods | |
| Sathyanarayana et al. | Characterisation of atherosclerotic plaque by spectral similarity of radiofrequency intravascular ultrasound signals | |
| JP7058264B2 (ja) | 血管内光干渉断層撮影法を用いた冠動脈プラークの自動特徴分析およびリスク評価のためのシステムおよび方法 | |
| US9741111B2 (en) | Subject information acquiring device and subject information acquiring method | |
| JP7467466B2 (ja) | マルチチャネル直交畳み込みニューラルネットワーク | |
| US8498462B2 (en) | Image processing apparatus and program | |
| US20030032880A1 (en) | Apparatus and method for ultrasonically identifying vulnerable plaque | |
| JP2011516865A (ja) | 化学分析モダリティーを補正する血管内構造分析のためのシステムおよび方法 | |
| CN108209955A (zh) | 斑块的表征 | |
| Cao et al. | Spectral analysis assisted photoacoustic imaging for lipid composition differentiation | |
| US20210113098A1 (en) | Image processing apparatus, method and storage medium to determine longitudinal orientation | |
| Flueraru et al. | Added soft tissue contrast using signal attenuation and the fractal dimension for optical coherence tomography images of porcine arterial tissue | |
| KR101049915B1 (ko) | 혈관 내 초음파 영상을 이용한 혈관 내 죽상 경화반 성분 측정 장치 및 방법 | |
| US20120230575A1 (en) | Quantification results in multiplane imaging | |
| MacKenzie et al. | B-mode ultrasound measurement of carotid bifurcation stenoses: is it reliable? | |
| Klingensmith et al. | Validation of an automated system for luminal and medial-adventitial border detection in three-dimensional intravascular ultrasound | |
| Courtney et al. | Effects of transducer position on backscattered intensity in coronary arteries | |
| CN108720870A (zh) | 一种基于超声衰减系数的脂肪肝检测系统 | |
| CN113974580A (zh) | 一种基于nirs技术的颈动脉狭窄程度分析方法 | |
| Rogatkin | Basic principles of organization of system software for multifunctional noninvasive spectrophotometric diagnostic devices and systems | |
| Villiger et al. | Future development | |
| Elamragy et al. | Additional diagnostic accuracy of global longitudinal strain at peak dobutamine stress in patients with moderate pretest probability of coronary artery disease | |
| He et al. | Comparison study of photoacoustic and ultrasound spectrum analysis in osteoporosis detection | |
| Kukk et al. | Integrated ultrasound and photoacoustic tomography for 3D imaging of human skin lesions |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230606 Address after: 200235 A6, Building 1, No. 900 Yishan Road, Xuhui District, Shanghai Patentee after: Shanghai Bodong Medical Technology Co.,Ltd. Address before: 200240 No. 800, Dongchuan Road, Shanghai, Minhang District Patentee before: SHANGHAI JIAO TONG University |
|
| TR01 | Transfer of patent right |