CN113180631A - 基于血管内成像的血流速度、血流储备分数的分析方法 - Google Patents
基于血管内成像的血流速度、血流储备分数的分析方法 Download PDFInfo
- Publication number
- CN113180631A CN113180631A CN202110474542.2A CN202110474542A CN113180631A CN 113180631 A CN113180631 A CN 113180631A CN 202110474542 A CN202110474542 A CN 202110474542A CN 113180631 A CN113180631 A CN 113180631A
- Authority
- CN
- China
- Prior art keywords
- blood
- blood flow
- image
- blood vessel
- flow velocity
- 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
- 230000017531 blood circulation Effects 0.000 title claims abstract description 110
- 238000003384 imaging method Methods 0.000 title claims abstract description 54
- 238000004458 analytical method Methods 0.000 title claims abstract description 44
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 83
- 238000012014 optical coherence tomography Methods 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 33
- 230000004089 microcirculation Effects 0.000 claims abstract description 16
- 238000002608 intravascular ultrasound Methods 0.000 claims abstract description 13
- 230000000747 cardiac effect Effects 0.000 claims description 22
- 210000004369 blood Anatomy 0.000 claims description 19
- 239000008280 blood Substances 0.000 claims description 19
- 238000011010 flushing procedure Methods 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000003795 chemical substances by application Substances 0.000 description 11
- 239000002872 contrast media Substances 0.000 description 10
- 238000001514 detection method Methods 0.000 description 8
- 210000004351 coronary vessel Anatomy 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 6
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 239000011780 sodium chloride Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000013146 percutaneous coronary intervention Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000013076 target substance Substances 0.000 description 3
- 208000031481 Pathologic Constriction Diseases 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000036770 blood supply Effects 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000002565 electrocardiography Methods 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000002504 physiological saline solution Substances 0.000 description 2
- 238000002601 radiography Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 208000037804 stenosis Diseases 0.000 description 2
- 230000036262 stenosis Effects 0.000 description 2
- 201000001320 Atherosclerosis Diseases 0.000 description 1
- 201000000057 Coronary Stenosis Diseases 0.000 description 1
- 244000208734 Pisonia aculeata Species 0.000 description 1
- 206010047141 Vasodilatation Diseases 0.000 description 1
- 230000003143 atherosclerotic effect Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000035487 diastolic blood pressure Effects 0.000 description 1
- 238000003113 dilution method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 238000010230 functional analysis Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 208000031225 myocardial ischemia Diseases 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035488 systolic blood pressure Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 231100000216 vascular lesion Toxicity 0.000 description 1
- 230000024883 vasodilation 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/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/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
- A61B5/0035—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 adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- 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
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
Abstract
本发明提供一种基于血管内成像的血流速度的分析方法、血流储备分数的分析方法、以及微循环阻力分析方法。其中基于血管内成像的血流速度的分析方法包括以下步骤：步骤S1，获取待检测血管的第一腔内影像，所述第一腔内影像为光学相干断层影像或血管内超声影像；步骤S2，基于所述第一腔内影像，获取血流信号，并基于所述血流信号计算所述待检测血管内的血流速度。根据本发明实施例的基于血管内成像的血流速度的分析方法，能够准确地测得血管内的动态血流值也即血流速度，能够获得更多的血管内血流信息，且可以用于优化计算FFR、微循环阻力的计算等。
Description
技术领域
本发明涉及待检测血管内介入影像医疗器械领域，具体涉及一种基于血管内成像的血流速度、血流储备分数、微循环阻力的分析方法。
背景技术
冠状动脉疾病已经成为全球第一大致死疾病。目前，经皮冠状动脉介入治疗(percutaneous coronary intervention，简称PCI)是冠状动脉疾病有效的治疗方法之一。
光学相干断层成像(OCT)、血管内超声(IVUS)等腔内影像技术，对于动脉粥样硬化病变斑块的识别，支架的放置与评估有着重大的临床意义，已经用于指导PCI治疗。其中因OCT具有高分辨率(纵向：10μm,横向：20-40μm)，在识别冠脉的结构形态特征及易损斑块的检测上具有绝对的优势，但仅局限于在形态与结构上对动脉粥样硬化及冠脉狭窄进行评估，并不能准确反映待检测血管形态变化引起的冠脉血流与供血等功能性的变化。
另一方面，血流储备分数(fractional flow reserve,简称FFR)通过测量冠状动脉狭窄段远端与近端的压力差值，可有效反映狭窄病变对待检测血管供血功能的影响，评估是否导致冠脉灌注心肌的缺血。
目前，FFR已成为临床上诊断，指导并评估PCI治疗的金标准。然而，传统的FFR需要通过压力导丝来测定血液压力，但是这种检查操作复杂、耗时较长，所需手术耗材(FFR导丝)价格昂贵。另外，注射待检测血管扩张药物产生的副作用会使患者出现不适反应，导丝在介入过程中容易造成病人待检测血管损伤。上述原因限制了FFR的推广。
此外，现有的通过影像学计算获得FFR的方法通常是基于使用固定的平均血流压力作为计算流体力学血流模型的边界条件。但是，实际的血流压力并非恒定不变的，被检测者个体化的差异以及下游血管的病变程度都影响着实际的血流情况。
因此，迫切需要一种能够准确地测得血管内的动态血流值的方法，以便进一步优化FFR计算，同时还可以用于了解心动周期中不同状态之间的压差、准确计算微循环阻力等。
发明内容
有鉴于此，本发明的目的在于提供一种简单、快速且能够准确计算基于血管内成像的血流速度的分析方法。
本发明还提供一种快速、准确的血流储备分数的分析方法。
此外，本发明还提供一种快速、准确的微循环阻力的分析方法。
为了实现上述目的，本发明采用以下技术方案：
根据本发明第一方面实施例的基于血管内成像的血流速度的分析方法，包括以下步骤：
步骤S1，获取待检测血管的第一腔内影像，所述第一腔内影像为光学相干断层影像或血管内超声影像；
步骤S2，基于所述第一腔内影像，获取血流信号，并基于所述血流信号计算所述待检测血管内的血流速度。
根据本发明实施例的血流速度的分析方法，能够准确地测得血管内的动态血流值也即血流速度能够获得更多的血管内血流信息。
其中，作为第一腔内影像，可以为光学相干断层影像(OCT影像)或血管内超声影像(IVUS影像)。其中OCT具有高分辨率(纵向：10μm,横向：20-40μm)，故而优选。
进一步地，所述步骤S2中，所述血流速度通过多普勒频移获得。也就是说，通过测量测量两次图像中相同位置处的所述血流信号的相位差来计算血流速度。通过多普勒频移能够准确地计算得到动态的血流速度。
进一步地，所述步骤S1包括：
步骤S11，将成像导管伸入所述待检测血管的预定位置；
步骤S12，使得所述成像导管在所述预定位置转动的同时在预定时间内以预定时间间隔扫描，获得所述第一腔内影像，其中，所述预定时间为一个或多个心动周期，
步骤S2中，基于相邻两次扫描得到的所述血流信号获取多普勒频移，基于所述多普勒频移确定所述血流速度。
也就是说，首先将成像导管伸入待检测血管的预定位置处，此后在预定时间内以预定时间间隔进行扫描，通过确定相邻两次扫描得到的所述血流信号即可获取多普勒频移，由此来确定血流速度。其中，由于在一个心动周期内的不同时期例如舒张压或收缩压状态下，其血流速度是不同的，因此通过监测一个或多个心动周期，能够获得更多的血流信息。
其中，所述心动周期例如通过x射线造影或心电图进行标定。也就是说，结合x射线造影技术或心电图技术来监测心动周期，并基于监测得到的心动周期的信息，确定所测得的血流速度所对应的心动周期的状态，能够获取更多的血管及血流信息。
进一步地，所述多普勒频移通过相同位置相邻时刻两次扫描的图像中的所述血流信号获得，或者，所述多普勒频移通过不同帧图像中相同位置处的所述血流信号获得。
也就是说，测量多普勒频移可以通过相同位置处相近时刻两次扫描的血流信号获得，例如成像导管位于所检测的待检测血管远端时趋于慢速旋转的状态，以保证相邻的A扫描之间有着足够的重叠，从而能够近似获得相同位置处，相近时刻信号的相位的变化。或者，多普勒频移通过不同帧图像中相同位置处的血流信号获得，因为不同帧图像之间也会存在着成像区域的重叠，所以也可获得不同帧图像相对应的同一位置的信号。
进一步地，所述第一腔内影像为光学相干断层影像，其中，所述步骤S12中，在扫描之前，向所述待检测血管注入第一预定量的冲洗剂，使得所述冲洗剂与血液混合。
当测量远端的血流速度时，首先将操作成像导管伸入至待检测血管的远端，因为血液对OCT光的强散射与吸收，因此需要向待检测血管中注入少量的冲洗剂与血液混合以便冲走部分血液从而提高光信号的穿透深度，作为冲洗剂例如可以为X射线造影剂或生理盐水。当冲洗剂与待检测血管中的血液中混合后，成像导管进行旋转的同时获取相同位置处的多个血流信号，基于多个血流信号获取多普勒频移，并基于多普勒频移确定最大血流速度。
现有血流速度测量方法已知有温度稀释法，将已知温度的冷盐水按一定速率通过导管注射至冠状动脉，通过插入在冠脉内的温度传感元件测量血流温度下降幅度和时间，来估算血流速度。如果使用该方法来进行血流速度检测，则血流速度检测、腔内影像获取需要分步进行。根据本实施例的分析方法，采用OCT影像来获取多普勒频移，在扫描之前，向待检测血管中注入少量的冲洗剂例如X射线造影剂、生理盐水来稀释待检测血管中的血液，从而能够获取其中的血流信号的多普勒频移，继而基于该多普勒频移获取动态的血流值。
具体地，通过在血液中加入适量的标的物(即X射线造影剂、生理盐水)，通过心跳带动该标的物进入冠脉到达测量位置即OCT镜头位置，由于先期到达的标的物与血液混合，这部分混合液体返回的OCT信号可由系统检测到，不同时刻获取的信号之间由于血液的流动产生多普勒频移，通过测定多普勒频移则可通过下述公式计算出混合液体的流速(也即血流速度)：
其中Vz为沿待检测血管方向流动的血流速度，λc为OCT光源的中心波长，
为两次检测之间的相位变化，ΔT为两次检测之间的时间间隔，n为被测液体的折射率，θ为导管出射光束与血流方向的夹角。
由计算的公式可以看出，由于已知中心波长、折射率、角度，只要确定两次检测的时间间隔，以及该两次检测所对应的多普勒频移(也即相位变化)，即可得到对应的血流速度。
当系统对于相位分辨的能力一定时，所能测得的血流速度的最大值以及精度将由时间间隔ΔT所决定。也就是说，通过调整时间间隔，可以获得更高的精度。另外，ΔT实际对应相邻线间扫描通过的时间间隔，如果ΔT设置过小，则受分辨率限制不能准确识别，此时可以通过适当增加时间间隔来减小采样密度。
根据本发明第二方面实施例的血流储备分数的分析方法，包括以下步骤：
步骤S10，获取待检测血管中的一个或多个心动周期内的血流速度的平均流速，所述血流速度根据上述任一项所述的血流速度的分析方法分析得到；
步骤S20，获取所述待检测血管的第二腔内影像，基于所述第二腔内影像以及所述平均流速进行流体力学分析计算，获得所述血流储备分数，
其中，所述第二腔内影像为光学相干断层影像或血管内超声影像。
基于腔内影像直接进行计算流体动力学仿真模拟时，通常将血流速度设定为一个常数。而实际血流速度因被检测者个体差异，病变状况，及心动周期等情况发生变化，实时的血流测量结果可以更准确地测定和分析待检测血管内的血流分析，进一步优化流体力学模拟和血流储备分数的计算。也就是说，根据本发明的血流储备分数的分析方法，通过第一腔内影像确定整个心动周期的血流速度的平均流速，并且结合该待检测血管的第二腔内影像进行流体动力学分析计算，获得血流储备分数(FFR)。相比于现有的FFR使用固定的平均血流压力作为计算流体力学血流模型的边界条件而言，基于待检测血管的腔内影像来计算其实时血流情况，结合腔内影像获得的待检测血管的具体的结构，可以更准确地计算血流储备分数。
进一步地，所述步骤S10中，所述血流速度包括所述待检测血管远端和/或近端的平均速度，所述步骤S20中，在获得所述第二腔内影像后，基于所述第二腔内影像进行血管分割，三维重建出血管的结构形态，并基于所述平均流速进行流体动力学分析计算，获得所述血流储备分数。
也就是说，既可以采用远端的平均速度，也可以采用近端的平均速度，还可以同时采用远端以及近端的平均速度来与由第二腔内影像得到的血管的结构心态进行结合，通过流体动力学分析计算获得FFR值。其中，远端位置处的血管较细，成像导管可能对其平均流速的影响较大，但是从程序上而言，远端测速、回撤获得第二腔内影像，可以一步完成，程序更加简化。因此可以结合具体条件选用远端平均流速、近端平均流速、远端与近端的综合平均流速。
进一步地，所述步骤S20包括：
使成像导管到达所述待检测血管的远端；
使所述成像导管旋转的同时回撤至所述待检测血管的近端，并获得所述待检测血管的所述第二腔内影像。
也就是说，作为获取第二腔内影像的具体流程，先将成像导管置于所述待检测血管的远端，此后使成像导管旋转的同时回撤至所述待检测血管的近端，回撤途中进行扫描获得反应了整个待检测血管的结构信息的第二腔内影像。
结合测量血流流速而言，既可以如下：先将成像导管到达血管远端，此后使得成像导管旋转获取该位置处的第一腔内影像，此后使得成像导管旋转的同时回撤，并进行扫描，获得第二腔内影像，基于第一腔内影像解析远端的血流速度的平均流速，基于第二腔内影像获得血管的结构形态，基于平均流速及结构形态，通过流体动力学解析，例如CFD仿真模拟，即可得到FFR值。
由此，通过一次回撤即可同时获取血流速度以及腔内影像，减少了测量步骤和时间，并大大提高了对血管病变功能学分析的准确率和效率。
进一步地，所述第二腔内影像为光学相干断层影像，其中，在使所述成像导管回撤之前，向所述待检测血管内注入第二预定量的冲洗剂以冲走所述血液。
也就是说，在采用OCT影像作为第二腔内影像时，在回撤之前，需要向血管内注入大量的冲洗剂(例如X光造影剂、生理盐水等)来冲走血液，也更准确地反应血管内的结构形态。
根据本发明第三方面实施例的微循环阻力的分析方法，包括如下步骤：
获取待检测血管近端的压力；
基于上述任一项所述的血流储备分数的分析方法分析血流储备分数；
基于近端的压力与所述血流储备分数，获得远端压力；
基于所述远端压力，以及基于上述任一项所述的血流速度的分析方法获得待检测血管的远端的一个或多个心动周期相应的血流速度的平均流速，确定所述待检测血管的所述微循环阻力。
也就是说，通过血流速度的平均流速来解析得到FFR值之后，结合血管近端的压力来确定远端压力，通过远端压力与平均流速即可确定微循环阻力。
附图说明
图1为本发明实施例1的基于血管内成像的血流速度的分析方法的流程示意图；
图2为使用高速OCT系统进行血流储备分数分析的示意图；
图3为根据实施例2的血流储备分数分析方法的流程示意图；
图4为根据实施例3的血流储备分数分析方法的流程示意图；
图5为根据实施例4的微循环阻力的分析方法的流程示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。
除非另作定义，本发明中使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本发明中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。同样，“一个”或者“一”等类似词语也不表示数量限制，而是表示存在至少一个。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电性的连接，不管是直接的还是间接的。“上”、“下”、“左”、“右”等仅用于表示相对位置关系，当被描述对象的绝对位置改变后，则该相对位置关系也相应地改变。
下面，通过具体实施例进一步详细说明本发明的血流速度的分析方法、血流储备分数的分析方法、以及微循环阻力的分析方法。
实施例1通过OCT影像的多普勒频移分析血流速度
如图1所示，本实施例的血流速度的分析包括如下步骤：
1a)将成像导管伸入带检测血管的预定位置处；
1b)向血管内注入少量的X光造影剂；
1c)在X光造影剂与血液的混合液到达该预定位置处时，使得成像导管原地转动，并以预定时间间隔例如1s进行扫描，
同时，通过心电图监测心动周期，确定所述成像导管获取1个或多个心动周期的第一腔内影像；
1d)对于相邻两次扫描得到的第一腔内影像，通过下述公式计算血流速度：
其中Vz为沿待检测血管方向流动的血流速度，λc为OCT光源的中心波长，
为两次检测之间的相位变化，ΔT为两次检测之间的时间间隔，n为被测液体的折射率，θ为导管出射光束与血流方向的夹角。
得到一整个心动周期中分别对应不同状态的血流速度。
实施例2通过OCT影像分析血流储备分数
如图2所示，通过高速OCT系统来分析血流储备分数。本实施例中采用高速OCT系统通过一次检测，获取血流速度的平均流速以及用于解析结构形态的腔内影像(即OCT回撤影像)。
其中，如图2所示，高速OCT系统包括OCT成像系统101、扫描装置102、电动控制机构103、以及导管104构成。其中，OCT成像系统101包括OCT光源、图像采集系统。
在获取OCT回撤影像的操作中，电动控制机构103带动导管104(也即带动成像镜头)运动到待检测血管105的远端，此后电动控制机构103控制扫描装置102进行360度旋转扫描，并通过回撤实现螺旋形3-D扫描，获得OCT图像。考虑到血液对OCT光的强散射和吸收，在导管104到达指定位置即远端后，向血管内注入大量冲洗剂例如x射线造影剂或生理盐水对血液进行冲洗，当冲洗剂到达远端(也就是，当OCT成像系统101能够清楚地获取OCT信号时)，通过电动控制机构103控制扫描装置102进行360度旋转扫描并回撤。
另外，为了获取血流信号的多普勒频移，具体操作如下：电动控制机构103带动导管104(也即带动成像镜头)运动到待检测血管105的预定位置例如远端，此后向血管内注入少量冲洗剂例如x射线造影剂或生理盐水以使得该冲洗剂与血液混合，当混合的血液到达远端，电动控制机构103控制扫描装置102在预定时间内(例如相当于1个或多个心动周期的时间)以预定的时间间隔进行360度旋转扫描，获取不同时刻的血流信号，从不同时刻的血流信号获取多普勒频移，基于多普勒频移来确定血流速度。
具体地，图3示出了本实施例的血流储备分数分析方法的流程示意图。
如图3所示，本实施例的血流储备分数分析方法包括如下步骤：
2A)获取待检测血管的远端的血流速度的整个心动周期内的平均流速；
具体地，采用跟实施例1同样的方法进行分析计算获得，
2B)向血管内注入大量X光造影剂以冲走血液；
2C)当X光造影剂到达远端，使得成像导管旋转的同时回撤，直至近端，以获得OCT回撤影像；
2D)对于OCT回撤影像进行血管分割，三维重建出血管的结构形态；
2E)基于上述得到的远端的平均流速、以及血管的结构形态，进行流体动力学仿真模拟，获得FFR值。
实施例3通过血管内超声影像(IVUS影像)分析血流储备分数
图4示出了本实施例的血流储备分数的分析方法的流程示意图。
如图4所示，根据本实施例的血流储备分数的分析方法，包括如下步骤：
3A)获取待检测血管的远端的血流速度的整个心动周期内的平均流速。
该步骤中，与实施例1不同的是，不需要在血管内注入冲洗剂。也就是说，无需执行1b)步骤，在成像导管到达远端后，可以立刻使得该成像导管旋转获取血管内超声影像，此后基于该血管内超声影像获取多普勒频移，并基于多普勒频移确定血流速度。具体的其他细节可以参考实施例1进行，在此不再赘述。
3B)使得成像导管旋转的同时回撤，直至近端，获得IVUS回撤影像；
3C)对于IVUS回撤影像进行血管分割，三维重建出血管的结构形态；
3D)基于上述得到的远端的平均流速、以及血管的结构形态，进行流体动力学仿真模拟，获得FFR值。
实施例4微循环阻力分析
图5示出了本实施例的微循环阻力分析的流程示意图。
如图5所示，根据本实施例的微循环阻力分析方法，包括如下步骤：
4A)确定平均流速V以及FFR值；
具体地，可以参考实施例1的方法获得平均速度V，并参考实施例2的方法获得FFR值，
4B)获取近端压力Pa，并计算远端压力Pd，其中Pd＝Pa*FFR；
其中，近端压力Pa可以通过压力导丝测量得到，也可以根据经验设定等。
4C)基于平均流速以及远端压力，计算微循环阻力IMR，其中，
IMR＝Pd/V。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种基于血管内成像的血流速度的分析方法，其特征在于，包括以下步骤：
步骤S1，获取待检测血管的第一腔内影像，所述第一腔内影像为光学相干断层影像或血管内超声影像；
步骤S2，基于所述第一腔内影像，获取血流信号，并基于所述血流信号计算所述待检测血管内的血流速度。
2.根据权利要求1所述的分析方法，其特征在于，所述步骤S2中，所述血流速度通过多普勒频移获得。
3.根据权利要求2所述的分析方法，其特征在于，所述步骤S1包括：
步骤S11，将成像导管伸入所述待检测血管的预定位置；
步骤S12，使得所述成像导管在所述预定位置转动的同时在预定时间内以预定时间间隔扫描，获得所述第一腔内影像，其中，所述预定时间为一个或多个心动周期；
步骤S2中，基于相邻两次扫描得到的所述血流信号获取多普勒频移，基于所述多普勒频移确定所述血流速度。
4.根据权利要求3所述的分析方法，其特征在于，所述心动周期通过x射线造影或心电图进行标定，
或者，所述多普勒频移通过相同位置相邻时刻两次扫描的图像中的所述血流信号获得，
或者，所述多普勒频移通过不同帧图像中相同位置处的所述血流信号获得。
5.根据权利要求3所述的分析方法，其特征在于，所述第一腔内影像为光学相干断层影像，其中，所述步骤S12中，在扫描之前，向所述待检测血管注入第一预定量的冲洗剂，使得所述冲洗剂与血液混合。
6.一种血流储备分数的分析方法，其特征在于，包括如下步骤：
步骤S10，获取待检测血管中的一个或多个心动周期内的血流速度的平均流速，所述血流速度根据权利要求1至5任一项所述的基于血管内成像的血流速度的分析方法分析得到；
步骤S20，获取所述待检测血管的第二腔内影像，基于所述第二腔内影像以及所述平均流速进行流体力学分析计算，获得所述血流储备分数，
其中，所述第二腔内影像为光学相干断层影像或血管内超声影像。
7.根据权利要求6所述的血流储备分数的分析方法，其特征在于，
所述步骤S10中，所述血流速度包括所述待检测血管远端和/或近端的平均速度，
所述步骤S20中，在获得所述第二腔内影像后，基于所述第二腔内影像进行血管分割，三维重建出血管的结构形态，并基于所述平均流速进行流体动力学分析计算，获得所述血流储备分数。
8.根据权利要求7所述的血流储备分数的分析方法，其特征在于，所述步骤S20包括：
使成像导管到达所述待检测血管的远端；
使所述成像导管旋转的同时回撤至所述待检测血管的近端，获得所述待检测血管的所述第二腔内影像。
9.根据权利要求8所述的血流储备分数的分析方法，其特征在于，所述第二腔内影像为光学相干断层影像，其中，在使所述成像导管回撤之前，向所述待检测血管内注入第二预定量的冲洗剂以冲走所述血液。
10.一种微循环阻力的分析方法，其特征在于，包括如下步骤：
获取待检测血管近端的压力；
基于权利要求6至9任一项所述的血流储备分数的分析方法分析血流储备分数；
基于近端的压力与所述血流储备分数，获得远端压力；
基于所述远端压力，以及权利要求1至4任一项所述的血流速度的分析方法获得待检测血管的远端的一个或多个心动周期相应的血流速度的平均流速，确定所述待检测血管的所述微循环阻力。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110474542.2A CN113180631A (zh) | 2021-04-29 | 2021-04-29 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
| JP2023535452A JP2023553556A (ja) | 2021-04-29 | 2022-04-27 | 血管内イメージングに基づく血流速度、血流予備量比の分析方法 |
| PCT/CN2022/089546 WO2022228464A1 (zh) | 2021-04-29 | 2022-04-27 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110474542.2A CN113180631A (zh) | 2021-04-29 | 2021-04-29 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113180631A true CN113180631A (zh) | 2021-07-30 |
Family
ID=76980914
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110474542.2A Pending CN113180631A (zh) | 2021-04-29 | 2021-04-29 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
Country Status (3)
| Country | Link |
|---|---|
| JP (1) | JP2023553556A (zh) |
| CN (1) | CN113180631A (zh) |
| WO (1) | WO2022228464A1 (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113616160A (zh) * | 2021-09-14 | 2021-11-09 | 昆山戎影医疗科技有限公司 | 基于多模态医学影像的ffr确定方法、装置、设备及介质 |
| WO2022228464A1 (zh) * | 2021-04-29 | 2022-11-03 | 上海博动医疗科技股份有限公司 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105283119A (zh) * | 2013-04-12 | 2016-01-27 | 皇家飞利浦有限公司 | 用于血流储备分数模拟的形状感测超声探头 |
| CN105559810A (zh) * | 2015-12-10 | 2016-05-11 | 上海交通大学 | 血管单位时间血流量与血流速度的计算方法 |
| CN106175700A (zh) * | 2016-09-27 | 2016-12-07 | 广东永士达医疗科技有限公司 | 一种应用于人体开放管腔道的oct探头 |
| CN107550519A (zh) * | 2017-08-22 | 2018-01-09 | 深圳先进技术研究院 | 一种多功能血管内超声成像装置 |
| CN107978371A (zh) * | 2017-11-30 | 2018-05-01 | 博动医学影像科技（上海）有限公司 | 快速计算微循环阻力的方法及系统 |
| CN111134651A (zh) * | 2019-12-09 | 2020-05-12 | 杭州脉流科技有限公司 | 基于腔内影像计算血流储备分数的方法、装置、系统以及计算机存储介质 |
| CN111134614A (zh) * | 2019-12-20 | 2020-05-12 | 苏州大学 | 基于oct测量眼球血管中血流的绝对速度的方法及系统 |
| CN111462117A (zh) * | 2020-05-19 | 2020-07-28 | 深圳北芯生命科技有限公司 | 基于血管图像的数据处理系统及其数据处理方法 |
| CN111754506A (zh) * | 2020-07-01 | 2020-10-09 | 杭州脉流科技有限公司 | 基于腔内影像的冠脉狭窄率计算方法、装置、系统和计算机存储介质 |
| CN212365043U (zh) * | 2020-05-19 | 2021-01-15 | 深圳北芯生命科技有限公司 | 基于血管图像的数据处理系统 |
| CN112450882A (zh) * | 2020-11-20 | 2021-03-09 | 深圳先进技术研究院 | 超声探头、内窥镜、内窥成像系统以及内窥成像方法 |
| CN112535466A (zh) * | 2020-12-16 | 2021-03-23 | 成都全景恒升科技有限公司 | 一种基于血管影像的血流储备分数计算方法 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10307056B2 (en) * | 2013-12-05 | 2019-06-04 | Bioptigen, Inc. | Systems and methods for quantitative doppler optical coherence tomography |
| US9687147B2 (en) * | 2014-12-02 | 2017-06-27 | Nidek Co., Ltd. | Optical coherence tomography device and control program |
| CN105286779B (zh) * | 2015-10-30 | 2017-06-30 | 温州医科大学 | 一种在体视网膜血流动力学的成像与绝对流速测量方法 |
| CN105326486B (zh) * | 2015-12-08 | 2017-08-25 | 博动医学影像科技（上海）有限公司 | 血管压力差与血流储备分数的计算方法及系统 |
| CN107730540B (zh) * | 2017-10-09 | 2020-11-17 | 全景恒升（北京）科学技术有限公司 | 基于高精度匹配模型的冠脉参数的计算方法 |
| CN109065170B (zh) * | 2018-06-20 | 2021-11-19 | 博动医学影像科技（上海）有限公司 | 获取血管压力差的方法及装置 |
| WO2020077381A1 (en) * | 2018-10-17 | 2020-04-23 | The University Of Melbourne | Coronary evaluation system and process |
| CN112384138B (zh) * | 2018-10-25 | 2023-07-25 | 苏州润迈德医疗科技有限公司 | 获取心表大动脉的血流量的方法、装置、系统及存储介质 |
| CN112155580B (zh) * | 2019-11-20 | 2021-07-13 | 苏州润迈德医疗科技有限公司 | 基于造影图像修正血流速度和微循环参数的方法及装置 |
| CN113180631A (zh) * | 2021-04-29 | 2021-07-30 | 博动医学影像科技（上海）有限公司 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
-
2021
- 2021-04-29 CN CN202110474542.2A patent/CN113180631A/zh active Pending
-
2022
- 2022-04-27 JP JP2023535452A patent/JP2023553556A/ja active Pending
- 2022-04-27 WO PCT/CN2022/089546 patent/WO2022228464A1/zh active Application Filing
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105283119A (zh) * | 2013-04-12 | 2016-01-27 | 皇家飞利浦有限公司 | 用于血流储备分数模拟的形状感测超声探头 |
| CN105559810A (zh) * | 2015-12-10 | 2016-05-11 | 上海交通大学 | 血管单位时间血流量与血流速度的计算方法 |
| CN106175700A (zh) * | 2016-09-27 | 2016-12-07 | 广东永士达医疗科技有限公司 | 一种应用于人体开放管腔道的oct探头 |
| CN107550519A (zh) * | 2017-08-22 | 2018-01-09 | 深圳先进技术研究院 | 一种多功能血管内超声成像装置 |
| CN107978371A (zh) * | 2017-11-30 | 2018-05-01 | 博动医学影像科技（上海）有限公司 | 快速计算微循环阻力的方法及系统 |
| CN111134651A (zh) * | 2019-12-09 | 2020-05-12 | 杭州脉流科技有限公司 | 基于腔内影像计算血流储备分数的方法、装置、系统以及计算机存储介质 |
| CN111134614A (zh) * | 2019-12-20 | 2020-05-12 | 苏州大学 | 基于oct测量眼球血管中血流的绝对速度的方法及系统 |
| CN111462117A (zh) * | 2020-05-19 | 2020-07-28 | 深圳北芯生命科技有限公司 | 基于血管图像的数据处理系统及其数据处理方法 |
| CN212365043U (zh) * | 2020-05-19 | 2021-01-15 | 深圳北芯生命科技有限公司 | 基于血管图像的数据处理系统 |
| CN111754506A (zh) * | 2020-07-01 | 2020-10-09 | 杭州脉流科技有限公司 | 基于腔内影像的冠脉狭窄率计算方法、装置、系统和计算机存储介质 |
| CN112450882A (zh) * | 2020-11-20 | 2021-03-09 | 深圳先进技术研究院 | 超声探头、内窥镜、内窥成像系统以及内窥成像方法 |
| CN112535466A (zh) * | 2020-12-16 | 2021-03-23 | 成都全景恒升科技有限公司 | 一种基于血管影像的血流储备分数计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 涂圣贤，等: "血管内光学成像与血流的融合：冠心病介入诊疗评估新技术研制", 《中国医学物理学杂志》, vol. 33, no. 12, pages 1210 - 1211 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022228464A1 (zh) * | 2021-04-29 | 2022-11-03 | 上海博动医疗科技股份有限公司 | 基于血管内成像的血流速度、血流储备分数的分析方法 |
| CN113616160A (zh) * | 2021-09-14 | 2021-11-09 | 昆山戎影医疗科技有限公司 | 基于多模态医学影像的ffr确定方法、装置、设备及介质 |
| WO2023040560A1 (zh) * | 2021-09-14 | 2023-03-23 | 苏州博动戎影医疗科技有限公司 | 基于多模态医学影像的ffr确定方法、装置、设备及介质 |
| CN113616160B (zh) * | 2021-09-14 | 2024-02-06 | 苏州博动戎影医疗科技有限公司 | 基于多模态医学影像的ffr确定方法、装置、设备及介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022228464A9 (zh) | 2023-12-07 |
| JP2023553556A (ja) | 2023-12-22 |
| WO2022228464A1 (zh) | 2022-11-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107730540B (zh) | 基于高精度匹配模型的冠脉参数的计算方法 | |
| JP7069236B2 (ja) | イメージングシステムの動作を制御する方法及びイメージを取得するシステム | |
| EP2453791B1 (en) | Apparatus for measuring flow and pressure within a vessel | |
| EP2854649B1 (en) | Fluid flow measurement systems and methods | |
| JP5997232B2 (ja) | 光干渉断層撮影プローブの作動方法、コンピュータシステム、及び光干渉断層撮影データ収集システム | |
| US8029446B2 (en) | Ultrasonic imaging system and imaging method | |
| US20140276059A1 (en) | Externally imaging a body structure within a patient | |
| WO2022228464A1 (zh) | 基于血管内成像的血流速度、血流储备分数的分析方法 | |
| JP5886581B2 (ja) | 血管プラーク画像診断装置 | |
| CN111462117A (zh) | 基于血管图像的数据处理系统及其数据处理方法 | |
| US11944778B2 (en) | Methods and systems for automatic pullback trigger | |
| CN212365043U (zh) | 基于血管图像的数据处理系统 | |
| Mackle et al. | Fibre optic intravascular measurements of blood flow: A review | |
| Pandian et al. | Intravascular ultrasound estimation of arterial stenosis | |
| EP4159112A1 (en) | Fluorescence calibration based on manual lumen detection | |
| CN217938192U (zh) | 一种多模态腔内成像系统 | |
| AU2010298337B2 (en) | Apparatus, systems, and methods of in-vivo blood clearing in a lumen | |
| RIZZON | 19. Transesophageal echocardiographic assessment of coronary flow reserve and coronary artery stenosis using echo enhancement CARLO CAIATI, SABINO ILICETO, NAVIN C. NANDA & |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 200233 16th floor, building 82, No. 1198, Qinzhou North Road, Xuhui District, Shanghai Applicant after: Shanghai Bodong Medical Technology Co.,Ltd. Address before: 200233 16th floor, building 82, No. 1198, Qinzhou North Road, Xuhui District, Shanghai Applicant before: PULSE MEDICAL IMAGING TECHNOLOGY (SHANGHAI) Co.,Ltd. |