CN116402795A - 一种血流速度的确定方法、装置、电子设备及存储介质 - Google Patents

一种血流速度的确定方法、装置、电子设备及存储介质 Download PDF

Info

Publication number
CN116402795A
CN116402795A CN202310372969.0A CN202310372969A CN116402795A CN 116402795 A CN116402795 A CN 116402795A CN 202310372969 A CN202310372969 A CN 202310372969A CN 116402795 A CN116402795 A CN 116402795A
Authority
CN
China
Prior art keywords
blood vessel
center line
image
centerline
determining
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
Application number
CN202310372969.0A
Other languages
English (en)
Inventor
张瑜
马骏
郑凌霄
兰宏志
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Raysight Intelligent Medical Technology Co Ltd
Original Assignee
Shenzhen Raysight Intelligent Medical Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Raysight Intelligent Medical Technology Co Ltd filed Critical Shenzhen Raysight Intelligent Medical Technology Co Ltd
Priority to CN202310372969.0A priority Critical patent/CN116402795A/zh
Publication of CN116402795A publication Critical patent/CN116402795A/zh
Pending legal-status Critical Current

Links

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/0275—Measuring blood flow using tracers, e.g. dye dilution
    • 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/0285—Measuring or recording phase velocity of blood waves
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • 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/10016—Video; Image sequence
    • 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/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physiology (AREA)
  • Molecular Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Hematology (AREA)
  • Veterinary Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Cardiology (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Mathematical Physics (AREA)
  • Fuzzy Systems (AREA)
  • Geometry (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本申请提供了一种血流速度的确定方法、装置、电子设备及存储介质，该确定方法包括：确定出第一造影图像中的第一血管中心线以及第二造影图像中的第二血管中心线，并将第一血管中心线映射到第二造影图像中，得到中心线映射图像；当第二血管中心线中不存在新的血管分支时，则基于第一血管中心线中的各个中心点以及第二血管中心线中的各个中心点在中心线映射图像中的坐标确定出新增血管段对应的血流速度；当第二血管中心线中存在新的血管分支时，则基于第一血管中心线中的各个中心点以及第二血管中心线中的各个中心点在中心线映射图像中的坐标确定出每个新增血管分段对应的血流速度。根据确定方法和确定装置，提高了血流速度计算的准确性和效率。

Description

一种血流速度的确定方法、装置、电子设备及存储介质
技术领域
本申请涉及医学图像处理技术领域，具体而言，涉及一种血流速度的确定方法、装置、电子设备及存储介质。
背景技术
血管造影图像中，血流速度是评估血管健康程度的一个重要功能学因素，比如医生在出具冠脉造影诊断报告的时候，会给予流速评级；在进行基于造影的FFR计算的时候，血流速度是一个重要的力学参数，准确的流速是精准计算结果的关键。充分利用造影动态序列的特征和不断充盈的冠脉树结构，计算得到全局流速，能够适应复杂的血管冠脉树实际情况。
目前，在现有的血流速度计算方法中，往往只计算单根血管的平均流速。而血管的实际情况复杂，由于多分支的存在，血流速度在血管树内不同位置会发生不同的变化，因此用单根血管的平均流速难以满足一些复杂的实际情况。
发明内容
有鉴于此，本申请的目的在于提供一种血流速度的确定方法、装置、电子设备及存储介质，在血管动态造影图像中，通过对图像的全局动态分析，充分利用造影动态序列的特征和不断充盈的血管树结构，获得动态的中心线长度变化，按照血管无分支和血管有分支两种策略进行局部流速计算，能够适应复杂的血管结构，提高了血流速度计算的准确性和效率。
第一方面，本申请实施例提供了一种血流速度的确定方法，所述确定方法包括：
获取血管造影视频，并将所述血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将所述血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像；
确定出所述第一造影图像中的第一血管中心线以及所述第二造影图像中的第二血管中心线，并将所述第一造影图像中的第一血管中心线映射到所述第二造影图像中，得到中心线映射图像；
当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中不存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度；
当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度。
进一步的，在所述得到中心线映射图像之后，所述确定方法还包括：
将所述第一造影图像中的第一导管区域和所述第二造影图像中的第二导管区域映射到所述中心线映射图像中；
在所述中心线映射图像中对所述第一导管区域进行骨架化，得到第一导管中心线，对所述第二导管区域进行骨架化，得到第二导管中心线；
计算所述第一导管中心线与第二导管中心线之间的平移向量，并在所述中心线映射图像中基于所述平移向量对所述第一导管中心线进行平移，将中心线平移后得到中心线映射图像作为所述中心线映射图像。
进一步的，通过以下步骤判断所述第二血管中心线中是否存在新的血管分支时：
确定所述第一血管中心线对应的第一端点数量，以及所述第二血管中心线对应的第二端点数量；
当所述第二端点数量大于所述第一端点数量时，则判断所述第二血管中心线中存在所述新的血管分支。
进一步的，所述基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度，包括：
确定出所述第二血管中心线中的每个中心点与所述第一血管中心线的血管终点之间的距离，将所述第二血管中心线中与所述第一血管中心线的血管终点之间距离最近的中心点确定为第一中心点，并将从所述第一中心点到所述第二血管中心线的血管终点之间的血管中心线作为所述新增血管段；
从所述第二血管中心线中的各个中心点中确定出所述第二血管中心线的血管起点，并利用所述第二血管中心线中从所述血管起点到所述第一中心点之间的其他中心点在所述中心线映射图像中的坐标计算出从所述血管起点到所述第一中心点的第一长度，并将所述第一长度作为前一帧血管长度；
基于所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标计算所述第二血管中心线的第二长度，并将所述第二长度作为当前帧血管长度；
将所述当前帧血管长度与所述前一帧血管长度之差确定为长度变化值，并将所述长度变化值与造影图像帧率之商作为所述新增血管段对应的血流速度。
进一步的，所述基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度，包括：
确定出所述第二血管中心线中的每个中心点与所述第一血管中心线的血管终点之间的距离，将所述第二血管中心线中与所述第一血管中心线的血管终点之间距离最近的中心点确定为第二中心点，并将从所述第二中心点到所述第二血管中心线的多个血管终点之间的血管中心线作为新增血管分支段；
将所述新增血管分支段划分为多个新增血管分段，并确定出每个新增血管分段对应的血管长度；
针对于每个新增血管分段，基于该新增血管分段对应的血管长度以及该新增血管分段对应的血流时间确定出该新增血管分段对应的血流速度方程表达式；
根据血管分叉前总流量与血管分叉后总流量之间的关系，联立多个所述血流速度方程表达式，得到每个新增血管分段对应的血流速度求解公式，并基于每个新增血管分段对应的血流速度求解公式求解出每个新增血管分段对应的血流速度。
进一步的，所述将所述新增血管分支段划分为多个所述新增血管分段，包括：
针对于所述新增血管分支段中的每个中心点，在所述新增血管分支段中确定与该中心点相邻的其他中心点的中心点数量；
当所述中心点数量大于或等于3时，则将该中心点确定为分支节点；
基于所述分支节点、所述新增血管分支段中的血管起点以及所述新增血管分支段中的多个血管终点，确定出多个所述新增血管分段。
进一步的，在所述确定出所述第二血管中心线中的新增血管段对应的血流速度，或者在所述确定出所述第二血管中心线中的每个新增血管分段对应的血流速度之后，所述确定方法还包括：
将所述第二造影图像作为所述第一造影图像，并返回执行所述将血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像的步骤，直至所述血管造影视频中不存在与所述第一造影图像相邻的下一帧造影图像。
第二方面，本申请实施例还提供了一种血流速度的确定装置，所述确定装置包括：
造影图像获取模块，用于获取血管造影视频，并将所述血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将所述血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像；
中心线映射图像确定模块，用于确定出所述第一造影图像中的第一血管中心线以及所述第二造影图像中的第二血管中心线，并将所述第一造影图像中的第一血管中心线映射到所述第二造影图像中，得到中心线映射图像；
第一血流速度确定模块，用于当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中不存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度；
第二血流速度确定模块，用于当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度。
第三方面，本申请实施例还提供一种电子设备，包括：处理器、存储器和总线，所述存储器存储有所述处理器可执行的机器可读指令，当电子设备运行时，所述处理器与所述存储器之间通过总线通信，所述机器可读指令被所述处理器执行时执行如上述的血流速度的确定方法的步骤。
第四方面，本申请实施例还提供一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行如上述的血流速度的确定方法的步骤。
本申请实施例提供的一种血流速度的确定方法、装置、电子设备及存储介质，首先，获取血管造影视频，并将所述血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将所述血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像；然后，确定出所述第一造影图像中的第一血管中心线以及所述第二造影图像中的第二血管中心线，并将所述第一造影图像中的第一血管中心线映射到所述第二造影图像中，得到中心线映射图像；当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中不存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度；当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度。
与现有技术中的血流速度确定方法相比，本申请在血管动态造影图像中，通过对图像的全局动态分析，充分利用造影动态序列的特征和不断充盈的血管树结构，获得动态的中心线长度变化，按照血管无分支和血管有分支两种策略进行局部流速计算，能够适应复杂的血管结构，提高了血流速度计算的准确性和效率。
为使本申请的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本申请实施例所提供的一种血流速度的确定方法的流程图；
图2为本申请实施例所提供的一种血流速度的确定装置的结构示意图之一；
图3为本申请实施例所提供的一种血流速度的确定装置的结构示意图之二；
图4为本申请实施例所提供的一种电子设备的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围，而是仅仅表示本申请的选定实施例。基于本申请的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的每个其他实施例，都属于本申请保护的范围。
首先，对本申请可适用的应用场景进行介绍。本申请可应用于医学图像处理技术领域。
血管造影图像中，血流速度是评估血管健康程度的一个重要功能学因素，比如医生在出具冠脉造影诊断报告的时候，会给予流速评级；在进行基于造影的FFR计算的时候，血流速度是一个重要的力学参数，准确的流速是精准计算结果的关键。充分利用造影动态序列的特征和不断充盈的冠脉树结构，计算得到全局流速，能够适应复杂的血管冠脉树实际情况。
经研究发现，目前，在现有的血流速度计算方法中，往往只计算单根血管的平均流速。而血管的实际情况复杂，由于多分支的存在，血流速度在血管树内不同位置会发生不同的变化，因此用单根血管的平均流速难以满足一些复杂的实际情况。
基于此，本申请实施例提供了一种血流速度的确定方法，以提高血流速度计算的准确性和效率。
请参阅图1，图1为本申请实施例所提供的一种血流速度的确定方法的流程图。如图1中所示，本申请实施例提供的血流速度的确定方法，包括：
S101，获取血管造影视频，并将所述血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将所述血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像。
需要说明的是，第一造影图像指的是在静脉中注入血管造影剂后，血管造影视频中出现了血管影像时的第一帧造影图像。第二造影图像则是在血管造影视频中与第一造影图像相邻的下一帧造影图像。
针对上述步骤S101，在具体实施时，获取血管造影视频，并将血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像。
S102，确定出所述第一造影图像中的第一血管中心线以及所述第二造影图像中的第二血管中心线，并将所述第一造影图像中的第一血管中心线映射到所述第二造影图像中，得到中心线映射图像。
需要说明的是，血管中心线指的是在血管造影图像中，在血管区域内，沿着血管中心走向的连线，可以表示为一系列的中心点像素。
针对上述步骤S102，在具体实施时，在确定出第一造影图像和第二造影图像后，在第一造影图像中确定出第一血管中心线，在第二造影图像中确定出第二血管中心点，然后将第一造影图像中的第一血管中心线映射到第二造影图像中，这时映射后的第二造影图像中既存在第一血管中心线又存在第二血管中心线，将映射后的第二造影图像确定为中心线映射图像。
针对上述步骤S102中，在确定血管中心线时，需要先在造影图像中获取血管起点。可选的，血管起点的获取方式可有如下多个替代方案：(1)人工选择血管起点；(2)通过有监督学习的方式，学习大量带标签的血管冠脉树起点，然后用神经网络模型来预测起点的位置；(3)通过导管的定位，将导管的尾端点作为血管起点。获得血管造影图像中的血管起点位置后，以血管起点沿着血管的“特征地图”进行区域生长，获取造影图像中的血管中心线。其中“特征地图”的获取可有如下多个替代方案：(1)直接用血管造影图像原图作为“特征地图”；(2)对血管造影图像原图进行frangi滤波增强，再对增强结果进行阈值分割得到二值化结果，作为“特征地图”；(3)通过有监督学习的方式，学习大量带二值分割标签的血管树，然后用神经网络模型来预测血管二值分割区域，作为“特征地图”。若以非二值图作为“特征地图”，则需要设定生长边界条件，比如设置一个地图范围[gmin,gmax]，其中gmin表示符合生长条件的最低地图灰度值，gmax表示符合生长条件的最高地图灰度值，超出该阈值范围则不满足生长条件，停止生长。若以二值图作为“特征地图”，则需要设定生长边界条件较为简单，即地图前景(灰度值＝1)为生长区域，背景区域(灰度值＝0)不满足生长条件，不在该区域生长。在造影图像中获取血管中心线的方法在现有技术中有详细说明，在此不再赘述。
具体的，针对上述步骤S102，在所述得到中心线映射图像之后，所述确定方法还包括：
步骤1021，将所述第一造影图像中的第一导管区域和所述第二造影图像中的第二导管区域映射到所述中心线映射图像中。
步骤1022，在所述中心线映射图像中对所述第一导管区域进行骨架化，得到第一导管中心线，对所述第二导管区域进行骨架化，得到第二导管中心线。
步骤1023，计算所述第一导管中心线与第二导管中心线之间的平移向量，并在所述中心线映射图像中基于所述平移向量对所述第一导管中心线进行平移，将中心线平移后得到中心线映射图像作为所述中心线映射图像。
当进行心脏冠脉造影时，由于心脏的运动，冠脉树区域在造影图像上的位置会随着时间发生有规律的变化，因此导致上一帧图像中的冠脉树血管中心线结构和当前帧图像中的冠脉树血管中心线结构存在一定的偏移和形变。所以需要进行上一帧血管中心线与当前帧血管中心线的匹配。针对上述步骤1021-步骤1023，在具体实施时，首先在第一造影图像中确定出导管对应的第一导管区域，在第二造影图像中确定出导管对应的第二导管区域，然后将第一造影图像中的第一导管区域和第二造影图像中的第二导管区域映射到中心线映射图像中，并在中心线映射图像中对第一导管区域进行骨架化，得到第一导管中心线，对第二导管区域进行骨架化，得到第二导管中心线。这里，进行骨架化的方法在现有技术中有详细说明，在此不再赘述。得到第一导管中心线和第二导管中心线后，计算第一导管中心线与第二导管中心线之间的平移向量。具体的，将第二导管中心线中与第一导管中心线的导管终点距离最近的导管中心点，与第一导管中心线的导管终点之间的距离作为平移向量。然后在中心线映射图像中基于平移向量对第一导管中心线进行平移，将中心线平移后得到中心线映射图像作为中心线映射图像。这样当血管中心线发生偏移时，可以对前后两帧的血管中心线进行匹配，提高后续血流速度计算的准确性。
作为一种可选的实施方式，上一帧血管中心线与当前帧血管中心线的匹配也可有如下两个方式：(1)基于血管的刚体配准，由于相邻两帧血管形变相对较小，可假设为刚体变换，仅存在旋转平移和缩放变换，会带来少量的误差，但计算速度快；(2)基于血管的非刚体配准，血管为软组织结构，这种方式比较符合实际情况，但其计算量大，计算时间长。刚体配准与非刚体配准的技术在现有技术中均有详细介绍，在此不再赘述。
在复杂的血管结构树中，会存在血管分支，为了保证血流速度计算的准确性，本申请按照无分支和有分支两种策略进行局部流速计算。因此根据本申请提供的血流速度的确定方法，在进行血流速度的之前，还需判断在第二帧造影图像中是否存在新的血管分支，也就是判断第二血管中心线中是否存在新的血管分支。
具体的，通过以下步骤判断所述第二血管中心线中是否存在新的血管分支时：
确定所述第一血管中心线对应的第一端点数量，以及所述第二血管中心线对应的第二端点数量；当所述第二端点数量大于所述第一端点数量时，则判断所述第二血管中心线中存在所述新的血管分支。
如果第二血管中心线与第一血管中心线相比存在新的血管分支时，第二血管中心线中会比第一血管中心线多出至少一个端点。针对上述两个步骤，在具体实施时，首先确定第一血管中心线对应的第一端点数量，以及第二血管中心线对应的第二端点数量。然后比较第一端点数量与第二端点数量，当第二端点数量大于第一端点数量时，则判断第二血管中心线中存在新的血管分支，则通过下述步骤S104计算新增血管分支段对应的血流速度。当第一端点数量与第二端点数量相等，则判断第二血管中心线中不存在新的血管分支，则通过下述步骤S103计算新增血管段对应的血流速度。
S103，当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中不存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度。
需要说明的是，新增血管段指的是与第一血管中心线相比，第二血管中心线中新增的血管段。
针对上述步骤S103，在具体实施时，当第二血管中心线与第一血管中心线相比，第二血管中心线中不存在新的血管分支时，则基于第一血管中心线中的各个中心点在中心线映射图像中的坐标以及第二血管中心线中的各个中心点在中心线映射图像中的坐标确定出第二血管中心线中的新增血管段对应的血流速度。
具体的，针对上述步骤S103，所述基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度，包括：
步骤1031，确定出所述第二血管中心线中的每个中心点与所述第一血管中心线的血管终点之间的距离，将所述第二血管中心线中与所述第一血管中心线的血管终点之间距离最近的中心点确定为第一中心点，并将从所述第一中心点到所述第二血管中心线的血管终点之间的血管中心线作为所述新增血管段。
针对上述步骤1031，在具体实施时，首先在中心线映射图像中确定出第一血管中心线的血管终点。一般地，血管造影的导管往往通过造影图像的边缘进入图像，所以血管中心线的起点会距离造影图像的边缘较近，而随着血管的延伸，血管中心线的终点会距离造影图像的边缘较远。因此将第一血管中心线中距离中心线映射图像的边缘较远的端点作为第一血管中心线的血管终点。然后计算出第二血管中心线中的每个中心点与第一血管中心线的血管终点之间的距离，将第二血管中心线中与第一血管中心线的血管终点之间距离最近的中心点确定为第一中心点。第一中心点确定出来后，将从第一中心点到第二血管中心线的血管终点之间的血管中心线作为新增血管段。第二血管中心线的血管终点的确定方法与上述第一血管中心线的血管终点的确定方法相同，在此不再赘述。具体的，通过下述两个公式确定出与第一血管中心线的血管终点之间距离最近的中心点：
Figure BDA0004169324960000131
Figure BDA0004169324960000132
其中，设第一血管中心线为M，第一血管中心线中的中心点的个数为m，M＝{M0,M1,M2,...,Mm}，Mm表示第一血管中心线M的血管终点，Mm＝(xMm,yMm)表示第一血管中心线M的血管终点在中心线映射图像中的坐标。设第二血管中心线为N，第二导管中心线中的中心点的个数为n，N＝{N0,N1,N2,...,Nn}。Nk表示第二血管中心线中的第k个中心点，(xNk,yNk)表示第二血管中心线中的第k个中心点在中心线映射图像中的坐标，Ni表示第一中心点。
步骤1032，从所述第二血管中心线中的各个中心点中确定出所述第二血管中心线的血管起点，并利用所述第二血管中心线中从所述血管起点到所述第一中心点之间的其他中心点在所述中心线映射图像中的坐标计算出从所述血管起点到所述第一中心点的第一长度，并将所述第一长度作为前一帧血管长度。
需要说明的是，前一帧血管长度指的是在第一血管造影图像对应的帧数中已充盈血管段的长度。
针对上述步骤1032，在具体实施时，从第二血管中心线中的各个中心点中确定出第二血管中心线的血管起点，并利用第二血管中心线中从血管起点到第一中心点之间的其他中心点在中心线映射图像中的坐标计算出从血管起点到第一中心点的第一长度，并将该第一长度作为前一帧血管长度。具体的，通过下述两个公式计算第一长度：
Figure BDA0004169324960000141
Figure BDA0004169324960000142
其中，L0表示前一帧血管长度，以上公式为从N0点逐点往后计算与下一个中心点之间的欧氏距离，然后累加，直至到达Ni点。
步骤1033，基于所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标计算所述第二血管中心线的第二长度，并将所述第二长度作为当前帧血管长度。
需要说明的是，当前帧血管长度指的是在第二血管造影图像对应的帧数中已充盈血管段的长度。
针对上述步骤1033，在具体实施时，基于第二血管中心线中的各个中心点在中心线映射图像中的坐标计算第二血管中心线的第二长度，并将第二长度作为当前帧血管长度。具体的，通过下述公式计算第二长度：
Figure BDA0004169324960000143
其中，L1表示当前帧血管长度。
步骤1034，将所述当前帧血管长度与所述前一帧血管长度之差确定为长度变化值，并将所述长度变化值与造影图像帧率之商作为所述新增血管段对应的血流速度。
需要说明的是，造影图像帧率指的是造影图像的帧率，即拍摄一帧图像造影所需要的时间。这里，造影图像帧率一般常见的为66.66毫秒。
针对上述步骤1034，在具体实施时，将当前帧血管长度与前一帧血管长度之差确定为长度变化值，并根据速度公式，将长度变化值与造影图像帧率之商作为新增血管段对应的血流速度。具体的，长度变化值和血流速度通过下述两个公式计算：
ΔL＝L1-Lo
Figure BDA0004169324960000151
其中，ΔL表示长度变化值，va表示新增血管段a段对应的血流速度，Δt表示造影图像帧率。
S104，当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度。
需要说明的是，新增血管分段指的是与第一血管中心线相比，第二血管中心线中新增的分支血管中的各个分支。
针对上述步骤S104，在具体实施时，当第二血管中心线与第一血管中心线相比，第二血管中心线中存在新的血管分支时，则基于第一血管中心线中的各个中心点在中心线映射图像中的坐标以及第二血管中心线中的各个中心点在中心线映射图像中的坐标确定出第二血管中心线中的每个新增血管分段对应的血流速度。
具体的，针对上述步骤S104，所述基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度，包括：
步骤1041，确定出所述第二血管中心线中的每个中心点与所述第一血管中心线的血管终点之间的距离，将所述第二血管中心线中与所述第一血管中心线的血管终点之间距离最近的中心点确定为第二中心点，并将从所述第二中心点到所述第二血管中心线的多个血管终点之间的血管中心线作为新增血管分支段。
针对上述步骤1041，在具体实施时，首先在中心线映射图像中确定出第一血管中心线的血管终点。然后，确定出第二血管中心线中的每个中心点与第一血管中心线的血管终点之间的距离，将第二血管中心线中与第一血管中心线的血管终点之间距离最近的中心点确定为第二中心点，并将从第二中心点到第二血管中心线的多个血管终点之间的血管中心线作为新增血管分支段。这里，对于如何确定第一血管中心线的血管终点以及如何确定与第一血管中心线的血管终点之间距离最近的中心点的描述可以参照上述步骤1031中的描述，并且能达到相同的技术效果，对此不做赘述。
步骤1042，将所述新增血管分支段划分为多个新增血管分段，并确定出每个新增血管分段对应的血管长度。
针对上述步骤1042，在具体实施时，将得到的新增血管分支段划分为多个新增血管分段，并确定出每个新增血管分段对应的血管长度。这里，对于如何确定新增血管分段对应的血管长度的描述可以参照上述步骤1032中计算血管长度的描述，并且能够达到相同的技术效果，对此不再赘述。
具体的，针对上述步骤1042，所述将所述新增血管分支段划分为多个所述新增血管分段，包括：
步骤10421，针对于所述新增血管分支段中的每个中心点，在所述新增血管分支段中确定与该中心点相邻的其他中心点的中心点数量。
步骤10422，当所述中心点数量大于或等于3时，则将该中心点确定为分支节点。
步骤10423基于所述分支节点、所述新增血管分支段中的血管起点以及所述新增血管分支段中的多个血管终点，确定出多个所述新增血管分段。
针对上述步骤10421-步骤10422，在具体实施时，针对于新增血管分支段中的每个中心点，在新增血管分支段中确定与该中心点相邻的其他中心点的中心点数量。当中心点数量大于或等于3时，则将该中心点确定为分支节点。然后，基于分支节点、新增血管分支段中的血管起点以及新增血管分支段中的多个血管终点，确定出多个新增血管分段。具体的，从新增血管分支段中的血管起点到分支节点之间的血管中心线作为一个新增血管分段，为上游段；从分支节点到新增血管分支段中的任意一个血管终点之间的血管中心线作为一个新增血管分段，为下游段。
步骤1043，针对于每个新增血管分段，基于该新增血管分段对应的血管长度以及该新增血管分段对应的血流时间确定出该新增血管分段对应的血流速度方程表达式。
需要说明的是，血流时间指的是在某个新增血管分段中，从起点流到终点所需的时间。
针对上述步骤1043，在具体实施时，针对于每个新增血管分段，基于该新增血管分段对应的血管长度以及该新增血管分段对应的血流时间确定出该新增血管分段对应的血流速度方程表达式。作为示例，新增血管分支段中有三个新增血管分段，分别是c0、c1和c2，其长度分别为c0、c1、c2，其中c0为上游段，c1和c2为下游段。首先假设从c0段的起点流到终点所需的血流时间为t0，那么c0段的速度
Figure BDA0004169324960000171
可表示为：
Figure BDA0004169324960000172
其中，
Figure BDA0004169324960000181
表示c0段对应的血流速度，c0表示c0段对应的血管长度，t0表示c0段对应的血流时间。
由于拍摄一帧图像造影所需要的时间是相同的，造影图像帧率为△t，那么从分叉节点流到c1和c2的时间也是相同的，为△t-t0，c1段和c2段的血流速度可表示为：
Figure BDA0004169324960000182
Figure BDA0004169324960000183
其中，
Figure BDA0004169324960000184
表示c1段对应的血流速度，/>
Figure BDA0004169324960000185
表示c2段对应的血流速度。
步骤1044，根据血管分叉前总流量与血管分叉后总流量之间的关系，联立多个所述血流速度方程表达式，得到每个新增血管分段对应的血流速度求解公式，并基于每个新增血管分段对应的血流速度求解公式求解出每个新增血管分段对应的血流速度。
需要说明的是，血管分叉前总流量与血管分叉后总流量之间的关系指的是血管分叉前的总流量和分叉后的总流量保持不变。
针对上述步骤1044，在具体实施时，根据血管分叉前总流量与血管分叉后总流量之间的关系，联立上述步骤1043中确定出的多个血流速度方程表达式，得到每个新增血管分段对应的血流速度求解公式，并基于每个新增血管分段对应的血流速度求解公式求解出每个新增血管分段对应的血流速度。在流体力学中，血管分叉前的总流量和分叉后的总流量保持不变，延续上述步骤1043中的例子，可以表示为：
QO＝Q1+Q2
Figure BDA0004169324960000186
Figure BDA0004169324960000187
Figure BDA0004169324960000188
其中Q0，Q1，Q2分别表示c0、c1和c2三段的血流量，S0，S1，S2分别表示c0、c1和c2三段的平均面积，而面积可表示为S＝πr2，r为半径，r0，r1，r2分别表示c0、c1和c2三段的平均直径。平均直径的获取可通过该段每个中心点所在的血管直径的平均值，从而除以2得到平均半径，，每个点的血管直径可由如下多个替代方案获得：若以非二值图为“特征地图”，以垂直于血管中心线方向上的两个梯度最大值作为血管的两个边缘，计算两个边缘的欧式距离作为血管的直径；若以二值图为“特征地图”，以垂直于血管中心线方向上的两个前景(像素值为1)边缘，计算两个边缘的欧式距离作为血管的直径。这时存在4个未知数t0，
Figure BDA0004169324960000191
和
Figure BDA0004169324960000192
而联立多个血流速度方程表达式，即可得到下述四个表达式：
Figure BDA0004169324960000193
Figure BDA0004169324960000194
Figure BDA0004169324960000195
Figure BDA0004169324960000196
其中，除了t0，
Figure BDA0004169324960000197
和/>
Figure BDA0004169324960000198
外，其他都是已知参数，可求出每个新增血管分段对应的血流速度求解公式：
Figure BDA0004169324960000199
Figure BDA00041693249600001910
/>
Figure BDA00041693249600001911
基于上述每个新增血管分段对应的血流速度求解公式，即可求解出每个新增血管分段对应的血流速度。
具体的，根据本申请提供的血流速度的确定方法，在所述确定出所述第二血管中心线中的新增血管段对应的血流速度，或者在所述确定出所述第二血管中心线中的每个新增血管分段对应的血流速度之后，所述确定方法还包括：
将所述第二造影图像作为所述第一造影图像，并返回执行所述将血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像的步骤，直至所述血管造影视频中不存在与所述第一造影图像相邻的下一帧造影图像。
血管造影视频中第二造影图像中新增血管段的血流速度计算出后，需继续计算与该第二造影图像相邻的下一帧造影图像中新增血管段的血流速度，这样如此往复向下计算，即可得到血管造影视频中每个帧造影图像中新增血管段的血流速度。针对上述步骤，在具体实施时，将第二造影图像作为第一造影图像，并返回执行步骤S101中所述将血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像的步骤，直至所述血管造影视频中不存在与所述第一造影图像相邻的下一帧造影图像。这样，所有血管造影视频中每个新增血管段的血流速度均计算出后，即可整合到整体血管结构树上，获得全局的血流速度。
本申请实施例提供的血流速度的确定方法，首先，获取血管造影视频，并将所述血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将所述血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像；然后，确定出所述第一造影图像中的第一血管中心线以及所述第二造影图像中的第二血管中心线，并将所述第一造影图像中的第一血管中心线映射到所述第二造影图像中，得到中心线映射图像；当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中不存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度；当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度。
与现有技术中的血流速度确定方法相比，本申请在血管动态造影图像中，通过对图像的全局动态分析，充分利用造影动态序列的特征和不断充盈的血管树结构，获得动态的中心线长度变化，按照血管无分支和血管有分支两种策略进行局部流速计算，能够适应复杂的血管结构，提高了血流速度计算的准确性和效率。
请参阅图2、图3，图2为本申请实施例所提供的一种血流速度的确定装置的结构示意图之一，图3为本申请实施例所提供的一种血流速度的确定装置的结构示意图之二。如图2中所示，所述确定装置200包括：
造影图像获取模块201，用于获取血管造影视频，并将所述血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将所述血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像；
中心线映射图像确定模块202，用于确定出所述第一造影图像中的第一血管中心线以及所述第二造影图像中的第二血管中心线，并将所述第一造影图像中的第一血管中心线映射到所述第二造影图像中，得到中心线映射图像；
第一血流速度确定模块203，用于当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中不存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度；
第二血流速度确定模块204，用于当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度。
进一步的，所述中心线映射图像确定模块202在用于得到中心线映射图像之后，所述中心线映射图像确定模块202还用于：
将所述第一造影图像中的第一导管区域和所述第二造影图像中的第二导管区域映射到所述中心线映射图像中；
在所述中心线映射图像中对所述第一导管区域进行骨架化，得到第一导管中心线，对所述第二导管区域进行骨架化，得到第二导管中心线；
计算所述第一导管中心线与第二导管中心线之间的平移向量，并在所述中心线映射图像中基于所述平移向量对所述第一导管中心线进行平移，将中心线平移后得到中心线映射图像作为所述中心线映射图像。
进一步的，如图3所示，所述确定装置200还包括分支判断模块205，所述分支判断模块205用于通过以下步骤判断所述第二血管中心线中是否存在新的血管分支时：
确定所述第一血管中心线对应的第一端点数量，以及所述第二血管中心线对应的第二端点数量；
当所述第二端点数量大于所述第一端点数量时，则判断所述第二血管中心线中存在所述新的血管分支。
进一步的，所述第一血流速度确定模块203在用于基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度时，所述第一血流速度确定模块203还用于：
确定出所述第二血管中心线中的每个中心点与所述第一血管中心线的血管终点之间的距离，将所述第二血管中心线中与所述第一血管中心线的血管终点之间距离最近的中心点确定为第一中心点，并将从所述第一中心点到所述第二血管中心线的血管终点之间的血管中心线作为所述新增血管段；
从所述第二血管中心线中的各个中心点中确定出所述第二血管中心线的血管起点，并利用所述第二血管中心线中从所述血管起点到所述第一中心点之间的其他中心点在所述中心线映射图像中的坐标计算出从所述血管起点到所述第一中心点的第一长度，并将所述第一长度作为前一帧血管长度；
基于所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标计算所述第二血管中心线的第二长度，并将所述第二长度作为当前帧血管长度；
将所述当前帧血管长度与所述前一帧血管长度之差确定为长度变化值，并将所述长度变化值与造影图像帧率之商作为所述新增血管段对应的血流速度。
进一步的，所述第二血流速度确定模块204在用于基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度时，所述第二血流速度确定模块204还用于：
确定出所述第二血管中心线中的每个中心点与所述第一血管中心线的血管终点之间的距离，将所述第二血管中心线中与所述第一血管中心线的血管终点之间距离最近的中心点确定为第二中心点，并将从所述第二中心点到所述第二血管中心线的多个血管终点之间的血管中心线作为新增血管分支段；
将所述新增血管分支段划分为多个新增血管分段，并确定出每个新增血管分段对应的血管长度；
针对于每个新增血管分段，基于该新增血管分段对应的血管长度以及该新增血管分段对应的血流时间确定出该新增血管分段对应的血流速度方程表达式；
根据血管分叉前总流量与血管分叉后总流量之间的关系，联立多个所述血流速度方程表达式，得到每个新增血管分段对应的血流速度求解公式，并基于每个新增血管分段对应的血流速度求解公式求解出每个新增血管分段对应的血流速度。
进一步的，所述第二血流速度确定模块204在用于将所述新增血管分支段划分为多个所述新增血管分段，所述第二血流速度确定模块204还用于：
针对于所述新增血管分支段中的每个中心点，在所述新增血管分支段中确定与该中心点相邻的其他中心点的中心点数量；
当所述中心点数量大于或等于3时，则将该中心点确定为分支节点；
基于所述分支节点、所述新增血管分支段中的血管起点以及所述新增血管分支段中的多个血管终点，确定出多个所述新增血管分段。
进一步的，如图3所示，在所述确定出所述第二血管中心线中的新增血管段对应的血流速度，或者在所述确定出所述第二血管中心线中的每个新增血管分段对应的血流速度之后，所述造影图像获取模块201还用于：
将所述第二造影图像作为所述第一造影图像，并返回执行所述将血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像的步骤，直至所述血管造影视频中不存在与所述第一造影图像相邻的下一帧造影图像。
请参阅图4，图4为本申请实施例所提供的一种电子设备的结构示意图。如图4中所示，所述电子设备400包括处理器410、存储器420和总线430。
所述存储器420存储有所述处理器410可执行的机器可读指令，当电子设备400运行时，所述处理器410与所述存储器420之间通过总线430通信，所述机器可读指令被所述处理器410执行时，可以执行如上述图1所示方法实施例中的血流速度的确定方法的步骤，具体实现方式可参见方法实施例，在此不再赘述。
本申请实施例还提供一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时可以执行如上述图1所示方法实施例中的血流速度的确定方法的步骤，具体实现方式可参见方法实施例，在此不再赘述。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-OnlyMemory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释，此外，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
最后应说明的是：以上所述实施例，仅为本申请的具体实施方式，用以说明本申请的技术方案，而非对其限制，本申请的保护范围并不局限于此，尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本申请实施例技术方案的精神和范围，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以权利要求的保护范围为准。

Claims (10)

1.一种血流速度的确定方法，其特征在于，所述确定方法包括：
获取血管造影视频，并将所述血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将所述血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像；
确定出所述第一造影图像中的第一血管中心线以及所述第二造影图像中的第二血管中心线，并将所述第一造影图像中的第一血管中心线映射到所述第二造影图像中，得到中心线映射图像；
当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中不存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度；
当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度。
2.根据权利要求1所述的确定方法，其特征在于，在所述得到中心线映射图像之后，所述确定方法还包括：
将所述第一造影图像中的第一导管区域和所述第二造影图像中的第二导管区域映射到所述中心线映射图像中；
在所述中心线映射图像中对所述第一导管区域进行骨架化，得到第一导管中心线，对所述第二导管区域进行骨架化，得到第二导管中心线；
计算所述第一导管中心线与第二导管中心线之间的平移向量，并在所述中心线映射图像中基于所述平移向量对所述第一导管中心线进行平移，将中心线平移后得到中心线映射图像作为所述中心线映射图像。
3.根据权利要求1所述的确定方法，其特征在于，通过以下步骤判断所述第二血管中心线中是否存在新的血管分支时：
确定所述第一血管中心线对应的第一端点数量，以及所述第二血管中心线对应的第二端点数量；
当所述第二端点数量大于所述第一端点数量时，则判断所述第二血管中心线中存在所述新的血管分支。
4.根据权利要求1所述的确定方法，其特征在于，所述基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度，包括：
确定出所述第二血管中心线中的每个中心点与所述第一血管中心线的血管终点之间的距离，将所述第二血管中心线中与所述第一血管中心线的血管终点之间距离最近的中心点确定为第一中心点，并将从所述第一中心点到所述第二血管中心线的血管终点之间的血管中心线作为所述新增血管段；
从所述第二血管中心线中的各个中心点中确定出所述第二血管中心线的血管起点，并利用所述第二血管中心线中从所述血管起点到所述第一中心点之间的其他中心点在所述中心线映射图像中的坐标计算出从所述血管起点到所述第一中心点的第一长度，并将所述第一长度作为前一帧血管长度；
基于所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标计算所述第二血管中心线的第二长度，并将所述第二长度作为当前帧血管长度；
将所述当前帧血管长度与所述前一帧血管长度之差确定为长度变化值，并将所述长度变化值与造影图像帧率之商作为所述新增血管段对应的血流速度。
5.根据权利要求1所述的确定方法，其特征在于，所述基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度，包括：
确定出所述第二血管中心线中的每个中心点与所述第一血管中心线的血管终点之间的距离，将所述第二血管中心线中与所述第一血管中心线的血管终点之间距离最近的中心点确定为第二中心点，并将从所述第二中心点到所述第二血管中心线的多个血管终点之间的血管中心线作为新增血管分支段；
将所述新增血管分支段划分为多个新增血管分段，并确定出每个新增血管分段对应的血管长度；
针对于每个新增血管分段，基于该新增血管分段对应的血管长度以及该新增血管分段对应的血流时间确定出该新增血管分段对应的血流速度方程表达式；
根据血管分叉前总流量与血管分叉后总流量之间的关系，联立多个所述血流速度方程表达式，得到每个新增血管分段对应的血流速度求解公式，并基于每个新增血管分段对应的血流速度求解公式求解出每个新增血管分段对应的血流速度。
6.根据权利要求5所述的确定方法，其特征在于，所述将所述新增血管分支段划分为多个所述新增血管分段，包括：
针对于所述新增血管分支段中的每个中心点，在所述新增血管分支段中确定与该中心点相邻的其他中心点的中心点数量；
当所述中心点数量大于或等于3时，则将该中心点确定为分支节点；
基于所述分支节点、所述新增血管分支段中的血管起点以及所述新增血管分支段中的多个血管终点，确定出多个所述新增血管分段。
7.根据权利要求1所述的确定方法，其特征在于，在所述确定出所述第二血管中心线中的新增血管段对应的血流速度，或者在所述确定出所述第二血管中心线中的每个新增血管分段对应的血流速度之后，所述确定方法还包括：
将所述第二造影图像作为所述第一造影图像，并返回执行所述将血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像的步骤，直至所述血管造影视频中不存在与所述第一造影图像相邻的下一帧造影图像。
8.一种血流速度的确定装置，其特征在于，所述确定装置包括：
造影图像获取模块，用于获取血管造影视频，并将所述血管造影视频中带有血管影像的第一帧造影图像作为第一造影图像，将所述血管造影视频中与所述第一造影图像相邻的下一帧造影图像作为第二造影图像；
中心线映射图像确定模块，用于确定出所述第一造影图像中的第一血管中心线以及所述第二造影图像中的第二血管中心线，并将所述第一造影图像中的第一血管中心线映射到所述第二造影图像中，得到中心线映射图像；
第一血流速度确定模块，用于当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中不存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的新增血管段对应的血流速度；
第二血流速度确定模块，用于当所述第二血管中心线与所述第一血管中心线相比，所述第二血管中心线中存在新的血管分支时，则基于所述第一血管中心线中的各个中心点在所述中心线映射图像中的坐标以及所述第二血管中心线中的各个中心点在所述中心线映射图像中的坐标确定出所述第二血管中心线中的每个新增血管分段对应的血流速度。
9.一种电子设备，其特征在于，包括：处理器、存储器和总线，所述存储器存储有所述处理器可执行的机器可读指令，当电子设备运行时，所述处理器与所述存储器之间通过所述总线进行通信，所述机器可读指令被所述处理器运行时执行如权利要求1至7任一所述的血流速度的确定方法的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器运行时执行如权利要求1至7任一所述的血流速度的确定方法的步骤。
CN202310372969.0A 2023-03-30 2023-03-30 一种血流速度的确定方法、装置、电子设备及存储介质 Pending CN116402795A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310372969.0A CN116402795A (zh) 2023-03-30 2023-03-30 一种血流速度的确定方法、装置、电子设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310372969.0A CN116402795A (zh) 2023-03-30 2023-03-30 一种血流速度的确定方法、装置、电子设备及存储介质

Publications (1)

Publication Number Publication Date
CN116402795A true CN116402795A (zh) 2023-07-07

Family

ID=87015607

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310372969.0A Pending CN116402795A (zh) 2023-03-30 2023-03-30 一种血流速度的确定方法、装置、电子设备及存储介质

Country Status (1)

Country Link
CN (1) CN116402795A (zh)

Similar Documents

Publication Publication Date Title
CN111815599B (zh) 一种图像处理方法、装置、设备及存储介质
Bauer et al. Segmentation of interwoven 3d tubular tree structures utilizing shape priors and graph cuts
Tang et al. Semiautomatic carotid lumen segmentation for quantification of lumen geometry in multispectral MRI
US9974442B2 (en) Method of, and apparatus for, processing volumetric image data
CN112716446B (zh) 高血压性视网膜病变的病变特征的测量方法及测量系统
WO2012174495A2 (en) Physics based image processing and evaluation process of perfusion images from radiology imaging
US20110081057A1 (en) Apparatus for stenosis estimation
KR20210010920A (ko) 심장 관류 자기 공명 이미징을 사용하여 허혈성 심장 질환을 검출하기 위한 완전 정량적 픽셀 단위 심근 혈류 및 심근 관류 예비능 맵을 자동으로 생성하고 분석하는 방법 및 시스템
CN113160189A (zh) 血管中心线提取方法、装置、设备及存储介质
CN113643354B (zh) 基于增强分辨率的眼底图像的血管管径的测量装置
CN117373070B (zh) 血管分段标注的方法及装置、电子设备和存储介质
Jiang et al. Vessel connectivity using Murray’s hypothesis
CN115049807A (zh) 肺部血管模型的建立方法、装置及服务器
US20230113721A1 (en) Functional measurements of vessels using a temporal feature
CN113469963A (zh) 肺动脉图像分割方法及装置
CN112651984A (zh) 血管管腔内膜轮廓提取方法、装置及超声设备和存储介质
CN117036530A (zh) 基于跨模态数据的冠状动脉血流储备分数预测方法及装置
CN116385540A (zh) 一种血管造影非充盈帧的导管定位方法和导管定位装置
KR101154355B1 (ko) 좌심실 자동 분할 방법
CN116402795A (zh) 一种血流速度的确定方法、装置、电子设备及存储介质
CN113421254B (zh) 微循环血管分支长度和直径计算方法及装置和终端设备
CN115841472A (zh) 大脑中动脉高密度征识别方法、装置、设备及存储介质
CN116309264A (zh) 造影图像判定方法及造影图像判定装置
CN111539926B (zh) 一种图像检测方法及装置
KR102000615B1 (ko) 관상동맥의 시작점을 자동으로 추출하기 위한 방법 및 그 장치

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