CN105095615B - 对血管中血流数据进行处理的方法和*** - Google Patents
对血管中血流数据进行处理的方法和*** Download PDFInfo
- Publication number
- CN105095615B CN105095615B CN201410165696.3A CN201410165696A CN105095615B CN 105095615 B CN105095615 B CN 105095615B CN 201410165696 A CN201410165696 A CN 201410165696A CN 105095615 B CN105095615 B CN 105095615B
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- vector
- point
- unit
- data
- 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
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明涉及一种对血管中血流数据进行处理的方法和***。其中,对血管中血流数据进行处理的方法包括:读取血管网格数据;根据血管网格数据生成血管三维模型;将血管血流参数加载到血管三维模型的对应位置。采用本发明的对血管中血流数据进行处理的方法和***,可以直观的体现出血管中的血流数据。
Description
技术领域
本发明涉及一种对血管中血流数据进行处理的技术,特别是对血管中血流数据进行处理的方法和***。
背景技术
冠状动脉性心脏病(简称冠心病)是动脉粥样硬化导致器官病变的最常见类型,也是严重危害人民健康的常见病,其特点是血液中的脂质沉着在原本光滑的动脉内膜上,在动脉内膜一些类似粥样的脂类物质堆积而成白色斑块(称为动脉粥样硬化病变),这些斑块渐渐增多造成动脉狭窄,使血流受阻,导致心脏缺血,产生心绞痛。若狭窄超过60%,在血液供应需求增加时就会出现缺血伴收缩功能减退、心律失常甚至心肌梗死、心脏病突发等。
目前,临床上用来辅助检查冠心病的技术分为侵入式检查和非侵入式检查两大类。侵入式检查包括:IVUS(血管内超声)和冠状动脉造影(Coronary Angiography),血管内超声可显示管壁的组织学特征,精确地测量管腔直径及截面积,发现冠状动脉造影不能显示的血管早期病变;冠状动脉造影是识别有无与粥样硬化性冠状动脉疾病有关的动脉狭窄的“金标准”,能较明确地揭示冠状动脉的解剖畸形及其阻塞性病变的位置、程度与范围。非侵入式检查包括:心电图、心脏损伤标记物检查、运动负荷试验、SPECT(心肌灌注的单光子发射计算机断层成像术)和PET(正电子发射断层成像术),这些检测技术通过测量心脏电位变化、心肌运动或者心肌灌注改变等间接地评估冠状动脉损伤;心脏磁共振显像(CMR)、冠状动脉CT血管造影(CT,电子计算机X射线断层扫描技术)可以得到冠状动脉的各个位置断层成像。
侵入式检查能获得丰富的与病变及狭窄相关的图像信息,但此种方法需要在患者体内置入导管等仪器,增加了患者的医疗费用和身体痛苦,而且导管对血管存在损伤,刺激冠状动脉和粥样硬化斑块也可能导致患者意外死亡。非侵入式检查则不需要在患者体内置入导管等仪器,危险性低,不会给患者身体带来额外的痛苦,但其获得的病情信息较少,不利于准确分析诊断患者病情。
因此,有必要发展一种使用非侵入性手段有效评估冠状动脉损伤状态的方法,这种方法不仅不会给患者身体造成侵入性伤害,还能获得丰富的患者病情信息。将计算流体力学方法和医学影像技术结合,则能实现这一目标,该方法以非侵入检查获取的患者病情信息为依据,经过流体力学计算,能获得冠状动脉血流储备分数(FFR)、血流压力、血液流速、血管壁应力等非常重要的病情诊断参数。
流体力学计算结果具有繁杂、多维的特点,若直接用于临床,作为医生诊断患者病情的依据,无疑会极大地增加医生的工作强度和难度。
发明内容
在下文中给出关于本发明的简要概述,以便提供关于本发明的某些方面的基本理解。应当理解,这个概述并不是关于本发明的穷举性概述。它并不是意图确定本发明的关键或重要部分,也不是意图限定本发明的范围。其目的仅仅是以简化的形式给出某些概念,以此作为稍后论述的更详细描述的前序。
本发明的一个主要目的在于提供一种对血管中血流数据进行处理的方法和***,其能将经过流体力学计算得到的血流参数进行直观展示。
根据本发明的一方面,一种对血管中血流数据进行处理的方法,所述血流数据包括血管网格数据和血管血流参数,包括:
读取血管网格数据;
根据所述血管网格数据生成血管三维模型;
将所述血管血流参数加载到所述血管三维模型的对应位置。
根据本发明的第二方面,一种对血管中血流数据进行处理的***,其中所述血流数据包括血管网格数据和血管血流参数,包括:
读取模块,用于读取血管网格数据;
三维重建模块,用于根据所述血管网格数据生成血管三维模型;
处理模块,用于将所述血管血流参数加载到所述血管三维模型的对应位置。
采用本发明的对血管中血流数据进行处理的方法和***,可以直观的体现出血管中的血流数据。
附图说明
参照下面结合附图对本发明实施例的说明,会更加容易地理解本发明的以上叙述和其它目的、特点和优点。附图中的部件只是为了示出本发明的原理。在附图中,相同的或类似的技术特征或部件将采用相同或类似的附图标记来表示。
图1为本发明对血管中血流数据进行处理的方法的一种实施方式的流程图;
图2为以血流储备分数为例表示对标量参数进行颜色标注法的效果图;
图3为对矢量参数进行箭头标注法的效果图;
图4为本发明对血管中血流数据进行处理的***的一种实施方式的结构图。
具体实施方式
下面参照附图来说明本发明的实施例。在本发明的一个附图或一种实施方式中描述的元素和特征可以与一个或更多个其它附图或实施方式中示出的元素和特征相结合。应当注意,为了清楚的目的,附图和说明中省略了与本发明无关的、本领域普通技术人员已知的部件和处理的表示和描述。
参见图1所示,为本发明对血管中血流数据进行处理的方法的一种实施方式的流程图。
在本实施方式中,对血管中血流数据进行处理的方法,其中的血流数据包括血管网格数据和血管血流参数。该方法包括:
S10:读取血管网格数据;
S20:根据血管网格数据生成血管三维模型;
S30:将血管血流参数加载到血管三维模型的对应位置。
其中,血管网格数据是对流体力学计算中所用血管三维模型的详细描述,包括节点数据和单元数据两部分。节点数据描述了血管三维模型中每一个坐标点的编号和坐标值,决定了模型的几何结构;单元数据描述了血管三维模型中每一个单元的类型、编号和节点组成,明确了模型中各节点之间的拓扑关系(即哪些点连成一条线,或者哪些点围成一个多边形),决定了模型的拓扑结构。以血管网格数据为依据,结合图形图像技术即可生成血管三维模型。
血管血流参数是对流体力学计算结果的记录,可以是标量,也可以是矢量,若为标量则具有编号和数值两个特征值,若为矢量则具有编号、x方向分量、y方向分量、z方向分量共4个特征值。血流参数的编号等于血管网格数据中某一单元的编号,根据这一对应关系,可以将流体力学计算结果作为属性数据加载到血管三维模型上,直观的表达出血管中的血流数据。
作为一种实施方式,血管血流参数可以包括标量参数和矢量参数,步骤S30可以具体包括:
S31:用颜色填充法和/或血流参数数值抽取法对标量参数进行标记;
S32:用箭头标注法和/或流线标注法对矢量参数进行标记。
标量参数例如可以包括血压和/或血流储备分数(FFR)。矢量参数例如可以包括血流速度和/或血管壁应力。
在一种实施方式中,颜色填充法例如可以具体包括:
a)创建一个颜色查询表,设定颜色查询表的色调范围,例如可以是红色到蓝色;
b)设定颜色查询表的颜色数量,例如可以是256种颜色,则颜色索引号范围为0~255;
c)确定血管三维模型中某一标量参数的最大值smax和最小值smin;
d)设定最小值smin所对应颜色为第一颜色(例如红色),最大值smax所对应颜色为第二颜色(例如蓝色);
e)采用公式index=255×(scalar-sminsmax-smin),计算血管三维模型中某一单元的标量数据scalar所对应的颜色索引号index;
f)与索引号index对应的颜色作为当前单元的颜色填充并显示出来。
颜色填充方法用不同的颜色值将不同的血流动力学参数值直观显示在心血管三维模型上,这样可以根据模型上某一处的颜色范围,轻松判断该处血流动力学参数的大致取值。
然而,颜色填充法仅仅能直观的获得血管三维模型中某标量参数的大致取值范围。若要获取模型上某一处某一标量参数的精确取值,则可以采用血流参数数值抽取方法实现。
在一种实施方式中,血流参数数值抽取法可以包括:
a)选取所述血管三维模型中需要进行数值抽取的区域;
b)将所述需要进行数值抽取区域中各单元的某一种血流参数数值求平均,并标注在该区域的中心处;
c)统计所述需要进行数值抽取区域中各单元的某一种血流参数的数值分布,并绘制数值分布直方图。
采用血流参数数值抽取法从庞大的数据量中抽取少量用户感兴趣的数据进行显示,忽略掉大部分冗余的、不感兴趣的数据,减少了需要可视化映射的数据量,同时保证了量的准确性。
在对矢量参数进行标记时,箭头标注法可以具体包括以下步骤:
a)按照每10个点采样一个的采样频率对矢量场(即由血管三维模型的某一类矢量属性数据形成的矢量场)采样;
b)由x方向分量、y方向分量、z方向分量计算各采样点矢量数据的大小和方向;
c)在各采样点上用一个有方向和长短的箭头来表示该点处矢量数据的方向和大小。
在对矢量参数进行标记时,流线标注法可以具体包括:
a)计算整个血管三维模型的范围值[xmin,ymin,zmin,xmax,ymax,zmax];
b)设定流线起点P0(x0,y0,z0),其坐标值须满足条件:xmin<x0<xmax,y0=ymax,zmin<z0<zmax;
c)设定积分步长step为单元长度(即三维空间坐标系中一个单位的长度)的0.05,最大增长量Time为500个时间单位,初始积分时长time为1;
d)采用公式x1=x0+step计算下一个积分点P1(即流线的下一节点)的x坐标值x1;
e)获取P0点所在单元的矢量数据作为速度矢量v0,获取点所在单元的矢量数据作为速度矢量v1y,获取点所在单元的矢量数据作为速度矢量v2y,获取(x1,y0+v2ystep,z0+step)点所在单元的矢量数据作为速度矢量v3y,采用公式计算下一个积分点P1(即流线的下一节点)的y坐标值y1;
f)获取点所在单元的矢量数据作为速度矢量v1z,获取点所在单元的矢量数据作为速度矢量v2z,获取(x1,y0+step,z0+v2zstep)点所在单元的矢量数据作为速度矢量v3z,并采用公式计算下一个积分点P1(即流线的下一节点)的z坐标值z1;
g)绘制点P0到点P1的线段;
h)对比积分时长time与最大增长量Time,若time≠Time则执行time=time+1,并将P1点作为流线起点P0,然后继续执行步骤d、e、f和步骤g,若time=Time则结束循环,流线绘制完毕。
如图2所示,为以血流储备分数(FFR)为例表示对标量参数进行颜色标注法的效果图。通过用不同的色调、饱和度、亮度和透明度来表示不同的FFR值,这样一来,可以直观的看出在三维模型各位置处,FFR的大致取值。
如图3所示,为对矢量参数进行箭头标注法的效果图。通过在三维模型的各位置上对矢量参数采用箭头进行标注,可以直观的体现出矢量参数的方向和大小。
作为一种优选方案,对血管中血流数据进行处理的方法还可以包括:
S40:将某一时间段内的某一种血流参数采用动画显示。
例如,将血管三维模型和特定时刻的某一种血流参数在同一可视化窗口中显示,得到特定时刻血管的状态分布图,直观显示某一时刻血管的不同位置处的属性值(血流储备分数、血流压力、血液流速、血管壁应力等),通过属性值的对比,可以容易的发现异常值。
在连续时间段内,将多个特定时刻的血管三维模型和某一种血流参数在同一可视化窗口中按时间顺序显示,可以得到连续时间段内血管的状态分布图,用来模拟一个或多个心动周期内,心血管的同一位置处的属性值(血流储备分数、血流压力、血液流速、血管壁应力等)的变化情况。由于各个单一计算结果的不同,用户会观看到连续一段时间内的血管三维模型的状态变化情况。
参见图4所示,为本发明一种对血管中血流数据进行处理的***,其中的血流数据包括血管网格数据和血管血流参数。
在本实施方式中,对血管中血流数据进行处理的***包括读取模块10、三维重建模块20和处理模块30。
其中,读取模块10用于读取血管网格数据。三维重建模块20用于根据血管网格数据生成血管三维模型。处理模块30用于将血管血流参数加载到血管三维模型的对应位置。
其中,血管网格数据是对流体力学计算中所用血管三维模型的详细描述,包括节点数据和单元数据两部分。节点数据描述了血管三维模型中每一个坐标点的编号和坐标值,决定了模型的几何结构;单元数据描述了血管三维模型中每一个单元的类型、编号和节点组成,明确了模型中各节点之间的拓扑关系(即哪些点连成一条线,或者哪些点围成一个多边形),决定了模型的拓扑结构。以血管网格数据为依据,结合图形图像技术即可生成血管三维模型。血管血流参数是对流体力学计算结果的记录,可以是标量,也可以是矢量,若为标量则具有编号和数值两个特征值,若为矢量则具有编号、x方向分量、y方向分量、z方向分量共4个特征值。血流参数的编号等于血管网格数据中某一单元的编号,根据这一对应关系,可以将流体力学计算结果作为属性数据加载到血管三维模型上,直观的表达出血管中的血流数据。
在一种实施方式中,血管血流参数可以包括标量参数和矢量参数。
作为一种优选方案,处理模块30可以包括标量参数处理部和矢量参数处理部。
标量参数处理部用于将血管血流参数中的标量参数加载到血管三维模型的对应位置;矢量参数处理部用于将血管血流参数中的矢量参数加载到血管三维模型的对应位置。
在一种实施方式中,标量参数处理部包括颜色填充模块和血流参数数值抽取模块。其中,颜色填充模块用于根据血管三维模型中各单元标量参数值的不同,对各单元上不同的颜色。
血流参数数值抽取模块,用于选取血管三维模型中需要进行数值抽取区域,将所述需要进行数值抽取区域中各单元某一种血流参数数值的平均值标注在该区域的中心处,并统计所需要进行数值抽取区域中各单元的某一种血流参数的数值分布,绘制数值分布直方图。
在一种实施方式中,矢量参数处理部包括箭头标注模块,用于在所述血管三维模型的预定数量的相邻单元中选取其中一个单元,并对该选取出单元的血管血流矢量参数进行大小和方向标注。
矢量参数处理部还可以包括流线标注模块,用于描述在由血管三维模型的矢量场中一个无质量粒子的运动轨迹(即流线),具体包括:
a)计算整个血管三维模型的范围值[xmin,ymin,zmin,xmax,ymax,zmax];
b)设定流线起点P0(x0,y0,z0),其坐标值须满足条件:xmin<x0<xmax,y0=ymax,zmin<z0<zmax;
c)设定积分步长step为单元长度(即三维空间坐标系中一个单位的长度)的0.05,最大增长量Time为500个时间单位,初始积分时长time为1;
d)采用公式x1=x0+step计算下一个积分点P1(即流线的下一节点)的x坐标值x1;
e)获取P0点所在单元的矢量数据作为速度矢量v0,获取点所在单元的矢量数据作为速度矢量v1y,获取点所在单元的矢量数据作为速度矢量v2y,获取(x1,y0+v2ystep,z0+step)点所在单元的矢量数据作为速度矢量v3y,采用公式计算下一个积分点P1(即流线的下一节点)的y坐标值y1;
f)获取点所在单元的矢量数据作为速度矢量v1z,获取点所在单元的矢量数据作为速度矢量v2z,获取(x1,y0+step,z0+v2zstep)点所在单元的矢量数据作为速度矢量v3z,并采用公式计算下一个积分点P1(即流线的下一节点)的z坐标值z1;
g)绘制点P0到点P1的线段;
h)对比积分时长time与最大增长量Time,若time≠Time则执行time=time+1,并将P1点作为流线起点P0,然后继续执行步骤d、e、f和步骤g,若time=Time则结束循环,流线绘制完毕。
作为一种优选方案,对血管中血流数据进行处理的***还可以包括显示模块40。
显示模块40用于将某一时间段内的血管某一种血流参数用动画显示。
例如,在连续时间段内,将多个特定时刻的血管三维模型和某一种血流参数在同一可视化窗口中按时间顺序显示,可以得到连续时间段内血管的状态分布图,用来模拟一个或多个心动周期内,心血管的同一位置处的属性值(血流储备分数、血流压力、血液流速、血管壁应力等)的变化情况。
采用本发明的对血管中血流数据进行处理的方法和***,可以直观的体现出血管中的血流数据。
上面对本发明的一些实施方式进行了详细的描述。如本领域的普通技术人员所能理解的,本发明的方法和装置的全部或者任何步骤或者部件,可以在任何计算设备(包括处理器、存储介质等)或者计算设备的网络中,以硬件、固件、软件或者它们的组合加以实现,这是本领域普通技术人员在了解本发明的内容的情况下运用他们的基本编程技能就能实现的,因此不需在此具体说明。
此外,显而易见的是,在上面的说明中涉及到可能的外部操作的时候,无疑要使用与任何计算设备相连的任何显示设备和任何输入设备、相应的接口和控制程序。总而言之,计算机、计算机***或者计算机网络中的相关硬件、软件和实现本发明的前述方法中的各种操作的硬件、固件、软件或者它们的组合,即构成本发明的设备及其各组成部件。
因此,基于上述理解,本发明的目的还可以通过在任何信息处理设备上运行一个程序或者一组程序来实现。所述信息处理设备可以是公知的通用设备。因此,本发明的目的也可以仅仅通过提供包含实现所述方法或者设备的程序代码的程序产品来实现。也就是说,这样的程序产品也构成本发明,并且存储或者传输这样的程序产品的介质也构成本发明。显然,所述存储或者传输介质可以是本领域技术人员已知的,或者将来所开发出来的任何类型的存储或者传输介质,因此也没有必要在此对各种存储或者传输介质一一列举。
在本发明的设备和方法中,显然,各部件或各步骤是可以分解、组合和/或分解后重新组合的。这些分解和/或重新组合应视为本发明的等效方案。还需要指出的是,执行上述系列处理的步骤可以自然地按照说明的顺序按时间顺序执行,但是并不需要一定按照时间顺序执行。某些步骤可以并行或彼此独立地执行。同时,在上面对本发明具体实施例的描述中,针对一种实施方式描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施方式中使用,与其它实施方式中的特征相组合,或替代其它实施方式中的特征。
应该强调,术语“包括/包含”在本文使用时指特征、要素、步骤或组件的存在,但并不排除一个或更多个其它特征、要素、步骤或组件的存在或附加。
虽然已经详细说明了本发明及其优点,但是应当理解在不超出由所附的权利要求所限定的本发明的精神和范围的情况下可以进行各种改变、替代和变换。而且,本发明的范围不仅限于说明书所描述的过程、设备、手段、方法和步骤的具体实施例。本领域内的普通技术人员从本发明的公开内容将容易理解,根据本发明可以使用执行与在此所述的相应实施例基本相同的功能或者获得与其基本相同的结果的、现有和将来要被开发的过程、设备、手段、方法或者步骤。因此,所附的权利要求旨在在它们的范围内包括这样的过程、设备、手段、方法或者步骤。
Claims (6)
1.一种对血管中血流数据进行处理的方法,所述血流数据包括血管网格数据和血管血流参数,其特征在于,包括:
读取血管网格数据;
根据所述血管网格数据生成血管三维模型;
将所述血管血流参数加载到所述血管三维模型的对应位置,具体包括:
对比血管血流参数编号和血管三维模型的各单元编号,将血管血流参数加载到与之具有相同编号的血管三维模型的单元中;
用颜色填充法和/或血流参数数值抽取法对标量参数进行标记;
用箭头标注法和/或流线标注法对矢量参数进行标记;
所述颜色填充法具体包括:
a)创建一个颜色查询表,设定颜色查询表的色调范围;
b)设定颜色查询表的颜色数量,且每一个颜色对应一个颜色索引号;
c)确定血管三维模型中某一标量参数的最大值smax和最小值smin;
d)设定最小值smin所对应颜色为红色,最大值smax所对应颜色为蓝色;
e)采用公式index=255×(scalar-smin/smax-smin),计算血管三维模型中某一单元的标量数据scalar所对应的颜色索引号index;
f)将与索引号index对应的颜色作为当前单元的颜色填充并显示出来;
所述血流参数数值抽取法具体包括:
a)选取所述血管三维模型中需要进行数值抽取的区域;
b)将所述需要进行数值抽取区域中各单元的某一种血流参数数值求平均,并标注在该区域的中心处;
c)统计所述需要进行数值抽取区域中各单元的某一种血流参数的数值分布,并绘制数值分布直方图。
2.根据权利要求1所述的对血管中血流数据进行处理的方法,其特征在于,
所述箭头标注法具体包括:
a)按照每10个点采样一个的采样频率对矢量场采样,所述矢量场为由血管三维模型的某一类矢量数据形成的矢量场;
b)由x方向分量、y方向分量、z方向分量计算采样点处的矢量数据的大小和方向;
c)在采样点上用一个有方向和长短的箭头来表示该点处矢量数据的方向和大小;
所述流线标注法具体包括:
a)计算整个血管三维模型的范围值[xmin,ymin,zmin,xmax,xmax,xmax];
b)设定流线起点P0(x0,y0,z0),其坐标值须满足条件xmin<x0<xmax,y0=ymax,zmin<z0<zmax;
c)设定积分步长step为单元长度的0.05,最大增长量Time为500个时间单位,初始积分时长time为1,所述单元长度为三维空间坐标系中一个单位的长度;
d)采用公式x1=x0+step计算下一个积分点P1的x坐标值x1;
e)获取P0点所在单元的矢量数据作为速度矢量v0,获取点所在单元的矢量数据作为速度矢量v1y,获取点所在单元的矢量数据作为速度矢量v2y,获取(x1,y0+v2ystep,z0+step)点所在单元的矢量数据作为速度矢量v3y,采用公式计算下一个积分点P1的y坐标值y1;
f)获取点所在单元的矢量数据作为速度矢量v1z,获取点所在单元的矢量数据作为速度矢量v2z,获取(x1,y0+step,z0+v2zstep)点所在单元的矢量数据作为速度矢量v3z,并采用公式计算下一个积分点P1的z坐标值z1;
g)绘制点P0到点P1的线段;
h)对比积分时长time与最大增长量Time,若time≠Time则执行time=time+1,并将P1点作为流线起点P0,然后继续执行步骤d、e、f和步骤g,若time=Time则结束循环,流线绘制完毕。
3.根据权利要求1-2任意一项所述的对血管中血流数据进行处理的方法,其特征在于,还包括:
对某一时间段内的某一种血流参数进行动画显示。
4.一种对血管中血流数据进行处理的***,所述血流数据包括血管网格数据和血管血流参数,其特征在于,包括:
读取模块,读取血管网格数据;
三维重建模块,根据所述血管网格数据生成血管三维模型;
处理模块,将所述血管血流参数加载到所述血管三维模型的对应位置,具体包括:
对比血管血流参数编号和血管三维模型的各单元编号,将血管血流参数加载到与之具有相同编号的血管三维模型的单元中;
用颜色填充法和/或血流参数数值抽取法对标量参数进行标记;
用箭头标注法和/或流线标注法对矢量参数进行标记;
所述处理模块包括标量参数处理部和矢量参数处理部;
所述标量参数处理部用于将所述血管血流参数中的标量参数加载到所述血管三维模型的对应位置;
所述矢量参数处理部用于将所述血管血流参数中的矢量参数加载到所述血管三维模型的对应位置;
所述标量参数处理部包括颜色填充模块,用于根据血管三维模型中各单元标量参数值的不同,对各单元上不同的颜色,具体包括:
a)创建一个颜色查询表,设定颜色查询表的色调范围;
b)设定颜色查询表的颜色数量,且每一个颜色对应一个颜色索引号;
c)确定血管三维模型中某一标量参数的最大值smax和最小值smin;
d)设定最小值smin所对应颜色为红色,最大值smax所对应颜色为蓝色;
e)采用公式index=255×(scalar-smin/smax-smin),计算血管三维模型中某一单元的标量数据scalar所对应的颜色索引号index;
f)与索引号index对应的颜色作为当前单元的颜色填充并显示出来;
所述标量参数处理部还包括血流参数数值抽取模块,用于选取所述血管三维模型中需要进行数值抽取区域,将所述需要进行数值抽取区域中各单元某一种血流参数数值的平均值标注在该区域的中心处,并统计所需要进行数值抽取区域中各单元的某一种血流参数的数值分布,绘制数值分布直方图。
5.根据权利要求4所述的对血管中血流数据进行处理的***,其特征在于,
所述矢量参数处理部包括箭头标注模块,用于在所述血管三维模型的预定数量的相邻单元中选取其中一个单元,并对该选取出单元的血管血流矢量参数进行大小和方向标注,具体包括:
a)按照每10个点采样一个的采样频率对矢量场采样,所述矢量场为由血管三维模型的某一类矢量数据形成的矢量场;
b)由x方向分量、y方向分量、z方向分量计算采样点矢量数据的大小和方向;
c)在采样点上用一个有方向和长短的箭头来表示该点处矢量数据的方向和大小;
所述矢量参数处理部还包括流线标注模块,用于描述在血管三维模型的矢量场中一个无质量粒子的运动轨迹,具体包括:
a)计算整个血管三维模型的范围值[xmin,ymin,zmin,xmax,xmax,xmax];
b)设定流线起点P0(x0,y0,z0),其坐标值须满足条件:xmin<x0<xmax,y0=ymax,zmin<z0<zmax;
c)设定积分步长step为单元长度的0.05,最大增长量Time为500个时间单位,初始积分时长time为1,所述单元长度为三维空间坐标系中一个单位的长度;
d)采用公式x1=x0+step计算下一个积分P1点的x坐标值x1;
e)获取P0点所在单元的矢量数据作为速度矢量v0,获取点所在单元的矢量数据作为速度矢量V1y,获取点所在单元的矢量数据作为速度矢量v2y,获取(x1,y0+v2ystep,z0+step)点所在单元的矢量数据作为速度矢量v3y,采用公式计算下一个积分点P1的y坐标值y1;
f)获取点所在单元的矢量数据作为速度矢量v1z,获取点所在单元的矢量数据作为速度矢量v1z,获取(x1,y0+step,z0+v2zstep)点所在单元的矢量数据作为速度矢量v3z,并采用公式计算下一个积分点P1的z坐标值z1;
g)绘制点P0到点P1的线段;
h)对比积分时长time与最大增长量Time,若time≠Time则执行time=time+1,并将P1点作为流线起点P0,然后继续执行步骤d、e、f和步骤g,若time=Time则结束循环,流线绘制完毕。
6.根据权利要求4-5中任意一项所述的对血管中血流数据进行处理的***,其特征在于,还包括显示模块;
所述显示模块用于对某一时间段内的某一种血流参数进行动画显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410165696.3A CN105095615B (zh) | 2014-04-23 | 2014-04-23 | 对血管中血流数据进行处理的方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410165696.3A CN105095615B (zh) | 2014-04-23 | 2014-04-23 | 对血管中血流数据进行处理的方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105095615A CN105095615A (zh) | 2015-11-25 |
CN105095615B true CN105095615B (zh) | 2019-05-21 |
Family
ID=54576039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410165696.3A Active CN105095615B (zh) | 2014-04-23 | 2014-04-23 | 对血管中血流数据进行处理的方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105095615B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106023202B (zh) * | 2016-05-20 | 2017-06-09 | 苏州润心医疗科技有限公司 | 基于心脏ct图像的冠状动脉血流储备分数计算方法 |
CN106327487B (zh) * | 2016-08-18 | 2018-01-02 | 苏州润迈德医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
US11039804B2 (en) * | 2016-09-16 | 2021-06-22 | Koninklijke Philips N.V. | Apparatus and method for determining a fractional flow reserve |
KR101809149B1 (ko) * | 2016-11-25 | 2017-12-14 | 한국과학기술연구원 | 순환계질환 발생잠재도를 판단하는 장치 및 그 방법 |
CN112116711B (zh) * | 2019-12-05 | 2024-01-23 | 苏州润迈德医疗科技有限公司 | 用于流体力学分析的圆台血管数学模型的合成方法和装置 |
CN113520363B (zh) * | 2021-06-25 | 2023-08-25 | 安影科技(北京)有限公司 | 一种下丘脑及垂体的血流灌注获取方法及*** |
CN114723684B (zh) * | 2022-03-22 | 2023-03-24 | 推想医疗科技股份有限公司 | 模型训练方法及其装置、血管结构生成方法及其装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101283929A (zh) * | 2008-06-05 | 2008-10-15 | 华北电力大学 | 一种血管三维模型的重建方法 |
EP2090864A1 (en) * | 2008-02-13 | 2009-08-19 | Konica Minolta Sensing, Inc. | Three-dimensional processor and method for controlling display ot three-dimensional data in the three-dimensional processor |
CN102411675A (zh) * | 2011-10-18 | 2012-04-11 | 浙江大学 | 一种在3d人体模型上监控人体生理数据的方法和*** |
CA2850189A1 (en) * | 2011-08-26 | 2013-03-07 | EBM Corporation | System for diagnosing blood flow characteristics, method thereof, and computer software program |
CN103247071A (zh) * | 2013-03-29 | 2013-08-14 | 哈尔滨工业大学深圳研究生院 | 一种构建三维血管模型方法及设备 |
CN103270513A (zh) * | 2010-08-12 | 2013-08-28 | 哈特弗罗公司 | 用于患者特异性血流建模的方法和*** |
-
2014
- 2014-04-23 CN CN201410165696.3A patent/CN105095615B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2090864A1 (en) * | 2008-02-13 | 2009-08-19 | Konica Minolta Sensing, Inc. | Three-dimensional processor and method for controlling display ot three-dimensional data in the three-dimensional processor |
CN101283929A (zh) * | 2008-06-05 | 2008-10-15 | 华北电力大学 | 一种血管三维模型的重建方法 |
CN103270513A (zh) * | 2010-08-12 | 2013-08-28 | 哈特弗罗公司 | 用于患者特异性血流建模的方法和*** |
CA2850189A1 (en) * | 2011-08-26 | 2013-03-07 | EBM Corporation | System for diagnosing blood flow characteristics, method thereof, and computer software program |
CN102411675A (zh) * | 2011-10-18 | 2012-04-11 | 浙江大学 | 一种在3d人体模型上监控人体生理数据的方法和*** |
CN103247071A (zh) * | 2013-03-29 | 2013-08-14 | 哈尔滨工业大学深圳研究生院 | 一种构建三维血管模型方法及设备 |
Non-Patent Citations (1)
Title |
---|
"应用于医学三维影像的血管结构自动提取";肖若秀 等;《光学精密工程》;20140215(第2期);第21-24页 |
Also Published As
Publication number | Publication date |
---|---|
CN105095615A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105095615B (zh) | 对血管中血流数据进行处理的方法和*** | |
JP5651600B2 (ja) | 再パラメータ化されたブルズアイ表示 | |
CN104871211B (zh) | 确定通过冠状动脉的血液流量的方法 | |
US10803995B2 (en) | Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in diseased and hypothetical normal anatomical models | |
CN101336844B (zh) | 医用图像处理装置以及医用图像诊断装置 | |
CN101036165B (zh) | 用于树模型显像以检测肺栓塞的***和方法 | |
US8335365B2 (en) | Diagnosis assisting apparatus, diagnosis assisting method, and storage medium having a diagnosis assisting program recorded therein | |
CN103228216B (zh) | 医用图像处理装置以及x射线ct装置 | |
Van Pelt et al. | Exploration of 4D MRI blood flow using stylistic visualization | |
CN118000760A (zh) | 用于可视化有风险的心脏组织的方法和*** | |
CN107582041A (zh) | 管状构造解析装置及管状构造解析方法 | |
Zhang et al. | Direct quantification of coronary artery stenosis through hierarchical attentive multi-view learning | |
CN105518684A (zh) | 用于预测冠状动脉病变的位置、开始、和/或变化的***和方法 | |
CN101887487B (zh) | 针对心脏病的模型产生器 | |
CN104486991B (zh) | 使用光纤形状感测对体素加标签 | |
CN107257655A (zh) | 用于利用从多谱段血液流动和灌注成像获取的信号确定血液动力学状态参数的方法、***和计算机程序产品 | |
JP2007202957A (ja) | 壁運動測定装置及び医用画像診断装置 | |
CN105074777B (zh) | 用于在诊断过程中支持用户的支持装置 | |
CN110494081A (zh) | 基于根据ecg信号对心肌微血管阻力的估计的冠状动脉疾病度量 | |
CN110706770B (zh) | 心脏数据处理设备及处理方法、计算机可读存储介质 | |
CN110731789B (zh) | 一种冠脉微血管阻力指数获取方法、介质及设备 | |
JP2013052245A (ja) | 情報処理装置および情報処理方法 | |
CN115115735B (zh) | 一种基于多期相冠脉ct造影的内皮动应变的快速计算*** | |
CN206557778U (zh) | 基于图形化书写结构化报告的*** | |
JP2015217113A (ja) | 血管解析装置、医用画像診断装置、血管解析方法及び血管解析プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |