CN107123112B - 血流状态分析系统及方法 - Google Patents
血流状态分析系统及方法 Download PDFInfo
- Publication number
- CN107123112B CN107123112B CN201710257724.8A CN201710257724A CN107123112B CN 107123112 B CN107123112 B CN 107123112B CN 201710257724 A CN201710257724 A CN 201710257724A CN 107123112 B CN107123112 B CN 107123112B
- Authority
- CN
- China
- Prior art keywords
- vessel
- model
- blood vessel
- blood
- region
- 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
- 238000000034 method Methods 0.000 title claims abstract description 142
- 230000017531 blood circulation Effects 0.000 title claims abstract description 84
- 238000004458 analytical method Methods 0.000 title claims abstract description 55
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 136
- 238000012545 processing Methods 0.000 claims description 87
- 210000004351 coronary vessel Anatomy 0.000 claims description 45
- 239000012530 fluid Substances 0.000 claims description 29
- 230000002792 vascular Effects 0.000 claims description 21
- 230000036772 blood pressure Effects 0.000 claims description 19
- 208000031481 Pathologic Constriction Diseases 0.000 claims description 16
- 210000001367 artery Anatomy 0.000 claims description 15
- 208000037804 stenosis Diseases 0.000 claims description 15
- 230000036262 stenosis Effects 0.000 claims description 15
- 230000003187 abdominal effect Effects 0.000 claims description 9
- 210000003141 lower extremity Anatomy 0.000 claims description 9
- 210000001627 cerebral artery Anatomy 0.000 claims description 7
- 230000008569 process Effects 0.000 abstract description 34
- 238000004364 calculation method Methods 0.000 description 37
- 238000003860 storage Methods 0.000 description 37
- 230000000875 corresponding effect Effects 0.000 description 34
- 230000000004 hemodynamic effect Effects 0.000 description 31
- 230000003993 interaction Effects 0.000 description 30
- 210000004369 blood Anatomy 0.000 description 20
- 239000008280 blood Substances 0.000 description 20
- 210000001519 tissue Anatomy 0.000 description 19
- 210000000056 organ Anatomy 0.000 description 18
- 230000011218 segmentation Effects 0.000 description 18
- 230000002159 abnormal effect Effects 0.000 description 17
- 238000004422 calculation algorithm Methods 0.000 description 16
- 238000003709 image segmentation Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 15
- 210000002216 heart Anatomy 0.000 description 15
- 238000004891 communication Methods 0.000 description 10
- 230000000747 cardiac effect Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 8
- 238000012986 modification Methods 0.000 description 8
- 230000004048 modification Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 230000009467 reduction Effects 0.000 description 7
- 230000008878 coupling Effects 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 6
- 238000005859 coupling reaction Methods 0.000 description 6
- 238000013480 data collection Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 238000002591 computed tomography Methods 0.000 description 5
- 238000009499 grossing Methods 0.000 description 5
- 238000012805 post-processing Methods 0.000 description 5
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000007405 data analysis Methods 0.000 description 4
- 230000010339 dilation Effects 0.000 description 4
- 230000002452 interceptive effect Effects 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 230000033764 rhythmic process Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000002604 ultrasonography Methods 0.000 description 4
- 230000005856 abnormality Effects 0.000 description 3
- 238000010968 computed tomography angiography Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000003902 lesion Effects 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 210000004165 myocardium Anatomy 0.000 description 3
- 238000002600 positron emission tomography Methods 0.000 description 3
- 206010002329 Aneurysm Diseases 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 2
- 208000007536 Thrombosis Diseases 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 238000002583 angiography Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 230000002490 cerebral effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000002526 effect on cardiovascular system Effects 0.000 description 2
- 201000011066 hemangioma Diseases 0.000 description 2
- 238000003706 image smoothing Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 238000012856 packing Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000002603 single-photon emission computed tomography Methods 0.000 description 2
- 238000012731 temporal analysis Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 208000002251 Dissecting Aneurysm Diseases 0.000 description 1
- 240000004282 Grewia occidentalis Species 0.000 description 1
- 208000018262 Peripheral vascular disease Diseases 0.000 description 1
- 206010057469 Vascular stenosis Diseases 0.000 description 1
- 206010047330 Vertebral artery stenosis Diseases 0.000 description 1
- 201000006559 Vertebrobasilar insufficiency Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 210000005242 cardiac chamber Anatomy 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 210000000038 chest Anatomy 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 208000019553 vascular disease Diseases 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
-
- 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
-
- 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
-
- 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
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- 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/021—Measuring pressure in heart or blood vessels
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- 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
-
- 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
Abstract
本发明提供一种血流状态分析方法和系统。在确定血流状态的过程中,获取多个时相的图像,建立与多个时相分别对应的多个血管模型。关联所述多个血管模型,根据关联结果分别设定所述多个血管模型的边界条件,确定所述多个血管模型的血管状态。
Description
【技术领域】
本发明涉及一种血流状态分析系统和方法,尤其是一种利用计算流体力学获取多时相血流动力学参数的系统及方法。
【背景技术】
CTA、MRA成像技术广泛应用在外周血管疾病诊断中,尤其在血管狭窄(椎动脉狭窄)、动脉瘤、夹层动脉瘤、肿瘤供血动脉等血管疾病等方面。血管分析软件应用提供了血管分析的工具,用于细小四级血管的精确提取、完整去骨、快速自动测量等。医学图像分析软件系统在针对血管的分析中一般采用图像分割技术,以及图像显示技术对相应对象的血管进行形态学上的三维模拟重建。医生通过针对血管形态学上的指标(血管狭窄程度、血管瘤膨胀程度等)对病变程度进行分析处理,可利用信息单一。
计算流体力学(Computational Fluid Dynamics/CFD)是一个20世纪50年代以来,随着计算机的发展而产生的一个介于数学、流体力学和计算机之间的交叉学科。其主要研究内容是通过计算机和数值方法来求解流体力学的控制方程,对流体力学问题进行模拟和分析。血管模型或血流模型是一个利用计算流体力学的新兴应用。利用计算流体力学对单时相数据进行分析难以全面的体现分析区域的实际情况和变化规律,也会由于选择时相不准确而使结果出现偏差。
【发明内容】
本发明解决的问题是提供一种血管流体分析方法,执行于至少一个包含处理器和存储器的设备,所述方法包括:获取一个第一时相的图像和一个第二时相的图像;在所述第一时相的图像中选择一个第一血管区域,其中,所述第一血管区域包含一个血管;在所述第二时相的图像中选择一个第二血管区域,其中,所述第二血管区域包含至少一部分所述血管;建立一个第一血管模型,其中,所述第一血管模型与所述第一血管区域相对应;建立一个第二血管模型,其中,所述第二血管模型与所述第二血管区域相对应;设定所述第一血管模型的边界条件和所述第二血管模型的边界条件;根据所述第一血管模型的边界条件,确定所述第一血管模型中所述血管在所述第一时相的状态;基于所述血管在所述第一时相的状态,关联所述第一血管模型和所述第二血管模型;以及所述关联结果以及所述第二血管模型的边界条件,确定所述第二血管模型中所述血管在所述第二时相的状态。
可选的,所述第一血管区域和所述第二血管区域包括冠状动脉、腹部动脉、大脑动脉或下肢动脉。
可选的,所述关联所述第一血管模型和所述第二血管模型包括关联所述第一血管模型和所述第二血管模型的所述血管的入口段、所述血管的分叉段、所述血管的狭窄段、或所述血管的出口段。
可选的,所述方法进一步包括对所述第一血管模型或所述第二血管模型进行网格化处理。
可选的,所述关联所述第一血管模型和所述第二血管模型包括将第一血管模型的网格和所述第二血管模型的网格进行匹配。
可选的,所述血管状态包括血流速、血压力、血管壁应力、血管壁切应力或血流储备系数(FFR)。
可选的,所述确定所述第一血管模型中所述血管在所述第一时相的状态或所述确定所述第二血管模型中所述血管在所述第二时相的状态包括使用计算流体力学(CFD)方法。
可选的,一个血流状态的分析系统,包括至少一个处理器和存储设备,所述系统包括:一个接收模块,被配置为:获取一个第一时相的图像和一个第二时相的图像;以及一个多时相特征生成模块,被配置为:在所述第一时相的图像中选择一个第一血管区域,其中,所述第一血管区域包含一个血管;在所述第二时相的图像中选择一个第二血管区域,其中,所述第二血管区域包含至少一部分所述血管;建立一个第一血管模型,其中,所述第一血管模型与所述第一血管区域相对应;建立一个第二血管模型,其中,所述第二血管模型与所述第二血管区域相对应;设定所述第一血管模型的边界条件和所述第二血管模型的边界条件;根据所述第一血管模型的边界条件,确定所述第一血管模型中所述血管在所述第一时相的状态;基于所述血管在所述第一时相的状态,关联所述第一血管模型和所述第二血管模型;以及根据所述关联结果以及所述第二血管模型的边界条件,确定所述第二血管模型中所述血管在所述第二时相的状态。
可选的,一种血流状态的分析方法,所述分析方法包括:获取多个时相的血管图像,包括一个第一时相的血管图像和一个第二时相的血管图像,所述多个时相的血管图像分别对应于同一个血管或其一部分;建立多个血管模型,所述多个血管模型分别对应于所述多个时相的血管图像;根据所述多个血管模型,获得所述血管或其一部分的多个状态,包括一个第一血管状态和一个第二血管状态,其中,所述第一血管状态对应于所述第一时相的血管图像,所述第二血管状态对应于所述第二时相的血管图像;根据所述血管或其一部分的所述多个状态,获得所述血管或其一部分的状态随时间变化的关系;以及根据所述关系,获得所述血管或其一部分的第三血管状态。
可选的,所述分析方法进一步包括:关联所述多个血管模型;以及根据所述关联结果,采用计算流体力学(CFD)分析所述多个血管状态,所述关联所述多个血管模型包括关联所述多个血管模型中至少两个血管模型的入口段、分叉段、狭窄段或出口段。
【附图说明】
图1A和图1B是根据本申请的一些实施例所示的包括血管状态分析系统;
图2是根据本申请的一些实施例所示的一个计算设备的结构,该计算设备可以实施本申请中披露的特定系统;
图3是根据本申请的一些实施例所示的一个移动设备的结构示意图,该移动设备可以实施本申请中披露的特定系统;
图4A是根据本申请的一些实施例所示的处理设备的示例性模块示意图;
图4B是根据本申请的一些实施例所示的多时相特征处理的示例性流程图;
图5是根据本申请的一些实施例所示的多时相特征生成模块的示例性模块示意图;
图6是根据本申请的一些实施例所示的获取多时相特征的示例性流程图;
图7是根据本申请的一些实施例所示的获取多时相特征的示意图;
图8是根据本申请的一些实施例所示的设定边界条件的示例性流程图;
图9是根据本申请的一些实施例所示的血流模型的示意图;
图10是根据本申请的一些实施例所示的网格化处理的示例性流程图;
图11是根据本申请的一些实施例所示的边界区域网格化处理的示意图;
图12是根据本申请的一些实施例所示的网格划分的示例性流程图;
图13是根据本申请的一些实施例所示的获得任意位点参数的示例性流程图;以及
图14是根据本申请的一些实施例所示的获得位点血流动力学参数的示意图。
【具体实施方式】
在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施的限制。
在图像数据处理过程中,“图像分割”、“图像提取”、“图像分类”可以相互转化,均表达从大范围区域内选取符合某条件的图像。在一些实施例中,成像系统可以包括一种或多种形态。所述形态包括但不限于,数字减影血管造影(DSA)、磁共振成像(MRI)、磁共振血管造影(MRA)、计算机断层扫描(CT)、计算机断层扫描血管造影(CTA)、超声波扫描(US)、正电子发射断层扫描术(PET)、单光子发射计算机断层扫描(SPECT)、SPECT-MR、CT-PET、CE-SPECT、DSA-MR、PET-MR、PET-US、SPECT-US、TMS-MR、US-CT、US-MR、X射线-CT、X射线-PET、X射线-US、视频-CT、视频-US和/或类似的一种或多种的组合。在一些实施例中,成像扫描的目标可以是器官、机体、物体、损伤部位、肿瘤等一种或多种的组合。在一些实施例中,成像扫描的目标可以是头部、胸腔、腹部、器官、骨骼、血管等一种或多种的组合。在一些实施例中,扫描的目标可以为一个或多个部位的血管组织。在一些实施例中,图像可以是二维图像和/或三维图像。在二维图像中,最细微的可分辨元素可以为像素点(pixel)。在三维图像中,最细微的可分辨元素可以为体素点(voxel)。在三维图像中,图像可由一系列的二维切片或二维图层构成。
图像分割过程可以基于图像的像素点(或体素点)的相应特征进行。在一些实施例中,所述像素点(或体素点)的相应特征可以包括纹理结构、灰度、平均灰度、信号强度、颜色饱和度、对比度、亮度等一种或多种的组合。在一些实施例中,所述像素点(或体素点)的空间位置特征也可以用于图像分割过程。
本申请涉及一种获取血流状态的方法和系统。在确定血流状态的过程中,获取多个时相的图像,建立与多个时相分别对应的多个血管模型。关联所述多个血管模型,根据关联结果分别设定所述多个血管模型的边界条件,确定所述多个血管模型的血管状态。所述关联所述多个血管模型包括配准多个血管模型。
图1A是根据本申请的一些实施例所示的包括血流状态分析系统100。该血流状态分析系统100可以包括数据采集设备110、处理设备120、存储设备130和显示设备140。数据采集设备110、处理设备120、存储设备130和交互设备140相互之间可以通过网络180进行通信。
数据采集设备110可以是一个采集数据的设备。所述数据可以包括图像数据、对象特征数据等。在一些实施例中,所述数据采集设备110可以包括一个成像设备。所述成像设备可以采集所述图像数据。所述成像设备可以是磁共振成像仪(magnetic resonanceimaging,MRI)、电子计算机断层扫描仪(computed tomography,CT)、正电子发射型计算机断层显像仪(positron emission computed tomography,PET)、B超仪(b-scanultrasonography)、超声诊断仪(diasonography)热断层扫描仪(Thermal texture maps,TTM)、医用电子内窥镜(medical electronic endoscope,MEE)等中的一种或多种的组合。所述图像数据可以是包括对象的血管、组织或器官的图片或数据。在一些实施例中,所述数据采集设备可以包括一个对象特征采集设备。所述对象特征采集设备可以采集对象的心率、心律、血压、血流速率、血液粘稠度、心输出量、心肌质量、血管流阻,以及/或其他与血管、组织或器官相关的对象特征数据。在一些实施例中,所述对象特征采集设备可以获取对象年龄、身高、体重、性别等其他对象特征数据。在一些实施例中,所述图像数据和对象特征数据可以是多时相数据。例如,所述多时相数据可以是在不同的时间点或时相获得的对象身上相同或近似位置的数据。在一些实施例中,所述对象特征采集设备可以集成在所述成像设备中,从而同时采集图像数据和对象特征数据。在一些实施例中,所述数据采集设备110可以通过网络180将其所采集的数据发送至处理设备120、存储设备130和/或交互设备140。
处理设备120可以对数据进行处理。所述数据可以是通过数据采集设备110采集到的数据,从存储设备130中读取的数据,从交互设备140中获得的反馈数据,如用户的输入数据,或通过网络180从云端或者外接设备中获得的数据等。在一些实施例中,所述数据可以包括图像数据、对象特征数据、用户输入数据等。所述处理可以包括在图像数据中选择感兴趣的区域。所述感兴趣的区域可以由处理设备120自行选择或根据用户输入数据选择。在一些实施例中,选择的感兴趣区域可以是血管、组织或者器官等。例如,所述感兴趣区域可以是动脉血管,如冠状动脉、腹部动脉、大脑动脉、下肢动脉等。处理设备120可以进一步对所述图像中对感兴趣的区域进行分割。图像分割的方法可以包括基于边缘的图像分割方法,如Perwitt算子法、Sobel算子法、梯度算子法、Kirch算子法等,基于区域的图像分割方法,如区域生长法、阈值法、聚类法等以及其他分割方法,如基于模糊集、神经网络的方法等。
处理设备120可以对所述感兴趣区域进行模型重建。模型的选择可以基于对象特征数据、感兴趣区域的特征等。例如,如果选定了感兴趣的区域为冠状动脉,处理设备120可以对包含冠状动脉的图像进行分割从而提取出冠状动脉的图像。然后,处理设备120可以根据对象特征、冠状动脉一般特征、冠状动脉图像特征等进行模型的重建。重建的模型可以与冠状动脉血管的形状相对应,也可以与冠状动脉中血液流动的形态相对应。在建立感兴趣区域的模型后,处理设备120可以根据模型进行分析和计算。所述分析和计算的方法可以包括计算流体力学(Computed fluid dynamics)等。
在一些实施例中,处理设备120可以获得多时相的数据,如对象在5个不同时间点上冠状动脉区域的图像。在这种情况下,处理设备120可以对不同时相的感兴趣区域(例如,整个冠状动脉,冠状动脉上的分支,或者冠状动脉的血液入口截面等)的图像分别构建模型,再对模型依次进行分析和计算。在一些实施例中,处理设备120可以对所述不同时相的模型进行网格化处理,并对网格化处理后的模型进行相互关联,从而降低计算量、提高计算准确度。关于网格化处理和模型相互关联的说明可以参见本申请其他地方的描述,例如,图6,图11及其描述。在一些实施例中,所述分析和计算的结果可以包括血管、组织、或器官的物理状态和相关系数或参数。例如,对冠状动脉模型进行分析和计算的结果可以包括冠状动脉的血流动力学参数,如血流速率、血液压力、血管壁应力、血管壁切应力、血流储备系数或血流储备分数(Fractional Flow Reserve,FFR)、冠状动脉血流储备(Coronary FlowReserve,CFR)等中的一种或多种的组合。在一些实施例中,处理设备120可以根据不同时相的分析和计算结果生成所述物理状态和/或相关系数或参数与时相或时间的关系(例如,血液动力学参数随时间的变化)。该关系可以用曲线或者对照表的方式体现。基于所述曲线或对照表,处理设备120可以获得任意时相的感兴趣区域的物理状态和/或相关系数或参数。
在一些实施例中,处理设备120可以对其获得的数据或处理结果进行降噪或平滑处理。在一些实施例中,处理设备120可以将其获得的数据或处理结果发送至存储设备130进行存储,或者发送至交互设备140进行显示。所述处理结果可以是处理过程中产生的中间结果,如感兴趣区域的模型,也可以是处理的最终结果,如分析和计算得出的血流动力学参数等。在一些实施例中,处理设备120可以是一个或多个处理元件或设备,如中央处理器(central processing unit,CPU)、图形处理器(graphics processing unit,GPU)、数字信号处理器(digital signal processor,DSP)、系统芯片(system on a chip,SoC)、微控制器(microcontroller unit,MCU)等。在一些实施例中,处理设备120也可以是特殊设计的具备特殊功能的处理元件或设备。处理设备120可以是本地的,或相对于数据采集设备110是远程的。
存储设备130可以储存数据或信息。所述数据或信息可以包括数据采集设备110获取的数据、处理设备120产生的处理结果或控制指令、以及交互设备140所接收到的用户输入数据等。存储设备130可以是一种或多种可以读取或写入的存储媒介,包括但不限于静态随机存储器(static random access memory,SRAM),随机存储器(random-access memory,RAM)、只读存储器(read-only memory,ROM)、硬盘、闪存等。在一些实施例中,存储设备130也可以是远程的存储器,如云盘等。
交互设备140可以接收、发送,以及/或显示数据或信息。所述接收的数据或信息可以包括数据采集设备110获取的数据、处理设备120产生的处理结果、存储设备130存储的数据等。例如,交互设备140显示的数据或信息可以包括数据采集设备110获得的心血管的实际图像150、处理设备120根据实际图像150所建立的心血管模型160,以及处理设备120从心血管模型160中提取出的冠状动脉模型170等。显示的形式可以包括但不限于二维或三维的医学图像、几何模型及其网格分析、矢量图(如速度矢量线)、等值线图、填充型的等值线图(云图)、XY散点图、粒子轨迹图、模拟流动效果等一种或多种组合。又例如,交互设备140发送的数据或信息可以包括用户的输入信息。交互设备140可以接收用户输入的处理设备120的一个或多个运行参数,并发送到处理设备120。
在一些实施例中,交互设备140可以包括一个用户交互界面。用户可以通过特定的交互装置,如鼠标、键盘、触摸板、麦克风等向交互设备140输入一个用户输入数据。例如,用户可以点击交互设备140所显示的模型并选择模型中感兴趣的区域。又例如,用户可以选择交互设备140所显示的血管模型中任意的位置,交互设备140可以从处理设备120获取并显示该位置的血流速、血压、血流量等。
在一些实施例中,交互设备140可以是显示屏等具有显示功能的设备。在一些实施例中,交互设备140可以具有处理设备120部分或全部的功能。例如,交互设备140可以对处理设备120生成的结果进行平滑、降噪、变色等操作。举例说明,变色操作可以将一个灰度图变成彩图,或将一个彩图变成一个灰度图。在一些实施例中,交互设备140与处理设备120可以是一个集成的设备。所述集成的设备可以同时实现处理设备120和交互设备140的功能。在一些实施例中,交互设备140可以包括台式电脑、服务器、移动设备等。移动设备可以包括笔记本电脑、平板电脑、ipad、交通工具(例如,机动车、船、飞机等)的内置设备、可穿戴设备等。在一些实施例中,交互设备140可以包括或连接到显示装置、打印机、传真等。
网络180可以用于血流状态分析系统100内部的通信,接收系统外部的信息,向系统外部发送信息等。在一些实施例中,数据采集设备110、处理设备120和交互设备140之间可以通过有线连接、无线连接、或其结合的方式接入网络180。网络180可以是单一网络,也可以是多种网络的组合。在一些实施例中,网络180可以包括但不限于局域网、广域网、公用网络、专用网络、无线局域网、虚拟网络、都市城域网、公用开关电话网络等中的一种或几种的组合。在一些实施例中,网络180可以包括多种网络接入点,例如有线或无线接入点、基站或网络交换点,通过以上接入点使数据源连接网络180并通过网络发送信息。
图1B所示的是一个血流状态分析系统100的另一个示意图。图1B与图1A类似。图1B中,处理设备120可以与数据采集设备110直接相连,而数据采集设备110不与网络180直接相连。
以上的描述仅仅是本发明的具体实施例,不应被视为是唯一的实施例。显然,对于本领域的专业人员来说,在了解本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变。例如,数据采集设备110、处理设备120、交互设备140之间可以不通过网络180而直接进行数据或信息的交换。又例如,这些设备也可以通过可移动存储器或其他中间媒介的方式进行数据或信息的交换。
图2是根据本申请的一些实施例所示的一个计算设备200的结构。该计算设备200可以实施本申请中披露的特定系统。本实施例中的特定系统利用功能框图解释了一个包含用户界面的硬件平台。计算设备200可以实施当前描述血流状态分析系统100中的一个或多个组件、模块、单元、子单元(例如,处理设备120,交互设备140等)。另外,血流状态分析系统100中的一个或多个组件、模块、单元、子单元(例如,处理设备120,交互设备140等)能够被计算设备200通过其硬件设备、软件程序、固件以及它们的组合所实现。这种计算机可以是一个通用目的的计算机,也可以是一个有特定目的的计算机。两种计算机都可以被用于实现本实施例中的特定系统。为了方便起见,图2中只绘制了一台计算设备,但是本实施例所描述的进行信息处理并推送信息的相关计算机功能是可以以分布的方式、由一组相似的平台所实施的,分散系统的处理负荷。
如图2所示,计算设备200可以包括内部通信总线210,处理器(processor)220,只读存储器(ROM)230,随机存取存储器(RAM)240,通信端口250,输入/输出组件260,硬盘270,用户界面280。内部通信总线210可以实现计算设备200组件间的数据通信。处理器220可以执行程序指令完成在此披露书中所描述的血流状态分析系统100的一个或多个功能、组件、模块、单元、子单元。处理器220由一个或多个处理器组成。通信端口250可以配置实现计算设备200与血流状态分析系统100其他部件(比如数据采集设备110)之间数据通信(比如通过网络180)。计算设备200还可以包括不同形式的程序储存单元以及数据储存单元,例如硬盘270,只读存储器(ROM)230,随机存取存储器(RAM)240,能够用于计算机处理和/或通信使用的各种数据文件,以及处理器220所执行的可能的程序指令。输入/输出组件260支持计算设备200与其他组件(如用户界面280),和/或与血流状态分析系统100其他组件(如数据库140)之间的输入/输出数据流。计算设备200也可以通过通信端口250从网络180发送和接收信息及数据。
图3描述了一种移动设备的结构,该移动设备能够用于实现实施本申请中披露的特定系统。在本例中,用于显示和交互位置相关信息的用户设备是一个移动设备300。移动设备300可以包括智能手机、平板电脑、音乐播放器、便携游戏机、全球定位系统(GPS)接收器、可穿戴计算设备(如眼镜、手表等),或者其他形式。本例中的移动设备300包括一个或多个中央处理器(CPUs)340,一个或多个图形处理器(graphical processing units(GPUs))330,一个显示320,一个内存360,一个天线310,例如一个无线通信单元,存储单元390,以及一个或多个输入/输出(input output(I/O))设备350。任何其他合适的组件,包括但不限于系统总线或控制器(图上未显示),也可能被包括在移动设备300中。如图3所示,一个移动操作系统370,如iOS、Android、Windows Phone等,以及一个或多个应用380可以从存储单元390加载进内存360中,并被中央处理器340所执行。应用380可能包括一个浏览器或其他适合在移动设备300上接收并处理图像或血液状态分析相关信息的移动应用。用户与血流状态分析系统100中一个或多个组件关于图像或血液状态分析相关信息的交互可以通过输入/输出系统设备350获得并提供给处理设备120,以及/或血流状态分析系统100中的其他组件,例如:通过网络180。
根据本申请的一些实施例,图4A是处理设备的示例性模块示意图。处理设备120可以包括接收模块410、控制模块420、多时相特征生成模块430、多时相特征处理模块440和输出模块450。
接收模块410可以从数据采集设备110和/或存储设备130获取图像数据、对象特征数据等。所述图像数据可以是包括对象的血管、组织或器官的图片或数据。所述对象特征数据可以包括对象的心率、心律、血压、血流速率、血液粘稠度、心输出量、心肌质量、血管流阻以及其他与血管、组织或器官相关的对象特征数据以及对象年龄、身高、体重、性别等其他对象特征数据。在一些实施例中,所述图像数据和对象特征数据可以是多时相数据。例如,所述多时相数据可以是在不同的时间点或时相获得的对象身上相同或近似位置的数据。
控制模块420可以发出控制指令。控制指令可以控制其他模块进行输入、输出、存储、处理等操作。例如,所述控制指令可以控制接收模块410接收所需要的数据。又例如,所述控制指令可以控制多时相特征生成模块430生成多时相的特征等。
多时相特征生成模块430可以生成多时相特征。所述多时相特征可以包括多时相模型、多时相参数、多时相边界条件、多时相的分析结果等。更具体的,多时相特征生成模块430可以在多时相图像数据中分别选择感兴趣的区域。所述感兴趣的区域可以由多时相特征生成模块430自行选择或根据用户输入数据选择。在一些实施例中,选择的感兴趣区域可以是血管、组织或者器官等。例如,所述感兴趣区域可以是动脉血管,如冠状动脉、腹部动脉、大脑动脉、下肢动脉等。所述多时相图像中选择的感兴趣区域可以是相对应的。例如,可以包含至少部分相同的血管、组织、或器官等。多时相特征生成模块430可以进一步对所述多时相图像中的感兴趣的区域进行分割。图像分割的方法可以包括基于边缘的图像分割方法(如Perwitt算子法、Sobel算子法、梯度算子法、Kirch算子法等),基于区域的图像分割方法(如区域生长法、阈值法、聚类法等),以及其他分割方法,如基于模糊集、神经网络的方法等。在一些实施例中,多时相特征生成模块430可以同时对多时相图像中感兴趣的区域进行分割。在一些实施例中,多时相特征生成模块430可以依次对多时相图像中感兴趣的区域进行分割。
多时相特征生成模块430可以对所述感兴趣区域进行模型重建,从而生成多时相模型。模型的选择可以基于对象特征数据、感兴趣区域的特征等。例如,如果选定了感兴趣的区域为冠状动脉,多时相特征生成模块430可以对包含冠状动脉的图像进行分割从而提取出冠状动脉的图像。然后,多时相特征生成模块430可以根据对象特征、冠状动脉一般特征、冠状动脉图像特征等进行模型的重建。重建的模型可以与冠状动脉血管的形状相对应,或与冠状动脉中血液流动的形态相对应。在建立感兴趣区域的模型后,多时相特征生成模块430可以设置参数和边界条件并根据模型进行分析和计算。具体的参数和边界条件设置的方法以及分析方法可以参见本申请其他部分的描述。
多时相特征处理模块440可以对生成的多时相计算结果进行处理(也被称为后处理)。所述处理可以包括利用拟合、插值等方法生成模型的计算结果与时相的关系曲线或对照表。根据所述关系曲线或对照表,多时相特征处理模块440可以进一步生成任意时相的分析结果的估计值。具体的后处理的流程和结果可以参考图13及其描述。在一些实施例中,所述多时相特征处理模块440可以将生成的多时相计算结果(例如,血管状态)与一个参考结果进行比对并生成一个比对结论。所述参考结果可以是存储在存储设备130中的数据,可以是存储在网络180中的数据,也可以是用户自行输入的数据。在一些实施例中,所述参考结果及相关比对结论可以存储在一个表中。例如,所述计算结果为血流速率时,所述参考结果可以为一个血流速率范围与危险程度的对应关系。所述危险程度可以分为正常、预警、危险、极度危险等。在一些实施例中,用户可以根据临床经验手动输入所述对应关系。在一些实施例中,所述比对可以是同一对象不同时期的血流速率的计算结果的比对。
输出模块450可以将生成的多时相计算结果或数据进行输出。例如,输出模块450可以将多时相计算结果或特征发送至存储设备130进行存储,或者发送至交互设备140进行显示。在一些实施例中,多时相特征处理模块440或输出模块450可以在输出前对所述多时相特征或计算结果进行降噪或平滑处理。所述多时相计算结果可以是生成的中间结果,如感兴趣区域的模型,或生成的最终结果,如分析和计算得出的血流动力学参数或计算结果与时相的关系曲线或对照表等。
根据本申请的一些实施例,图4B是多时相特征处理的示例性流程图。在一些实施例中,流程可以通过处理设备120实现。
在462中,可以产生一个或多个控制指令。在一些实施例中,462可以通过控制模块420实现。所述控制指令可以控制流程中其他步骤的进行。
在464中,可以接收多时相数据。在一些实施例中,464可以通过接收模块410实现。所述多时相数据可以包括多时相图像数据和多时相对象特征数据。在一些实施例中,所述多时相对象特征数据可以是时间上连续的对象特征数据或特征曲线。
在466中,可以生成多时相特征。在一些实施例中,466可以通过多时相特征生成模块430实现。所述多时相特征可以包括多时相模型、多时相参数、多时相边界条件、多时相的分析结果等。
在468中,可以对生成的多时相特征进行处理。在一些实施例中,468可以通过多时相特征处理模块440实现。所述处理可以包括利用拟合、插值等方法生成多时相特征与时相的关系曲线或对照表。
在470中,可以输出多时相特征或处理结果。在一些实施例中,470可以通过输出模块450实现。在一些实施例中,可以跳过468,直接将生成的多时相特征输出。
根据本申请的一些实施例,图5是多时相特征生成模块的示例性模块示意图。多时相特征生成模块430可以包括数据获取单元510、参数设置单元520、计算单元530、网格生成单元540、匹配单元550、区域选择单元560、输出单元570和判断单元580。
数据获取单元510可以从多时相特征生成模块430中其他单元、血流状态分析系统100中其他设备或模块或外界设备或模块中获取数据。所述数据可以包括图像数据、对象特征数据、用户输入数据等。所述图像数据可以是包括对象的血管、组织或器官的图片或数据。所述对象特征数据可以包括对象的心率、心律、血压、血流速率、血液粘稠度、心输出量、心肌质量、血管流阻以及其他与血管、组织或器官相关的数据。在一些实施例中,所述图像数据和对象特征数据可以是多时相数据。在一些实施例中,数据获取单元510可以从存储设备130中获取经过处理后的数据,如重建好的血管模型等。在一些实施例中,数据获取单元510可以对获取的图像数据进行预处理。所述预处理可以包括图像增强、图像降噪、图像平滑等。
参数设置单元520可以选择模型以及设置参数和边界条件。所述模型的选择可以包括根据具体分析的病变部位(感兴趣区域)以及对象特征数据(如血液粘稠度等)选择适合的血液粘性模型和流速边界模型。所述血液粘性模型可以包括牛顿流体模型、非牛顿流体模型以及用户自定义的其他流体模型。所述牛顿流体模型可以用来模拟对象体内血液粘稠度较为恒定的区域,而非牛顿流体模型可以用来模拟对象体内血液粘稠度不恒定的区域。所述流速边界模型可以包括抛物线模型、双曲线模型、椭圆模型、平均流模型、Womersley分布模型、Reynolds模型、混合模型等。在一些实施例中,所述参数的设置可以包括选定的模型中参数的设置,例如牛顿流体模型中的血液粘性系数,牛顿流体模型中的血液密度,模拟计算中的时间步数,模拟计算中的时间步长等。
所述边界条件的设置可以包括对边界区域进行初始条件或限制条件的设置。所述边界区域指的是感兴趣区域的边缘区域。例如,如果选定的感兴趣区域是一个血管区域或与血管对应的血流区域,边界区域可以是血管的出口、入口、血管壁等。设置的边界条件可以包括边界区域的血压、血流速率、流阻、压强、应力等。在一些实施例中,存储设备130或血流状态分析系统100内部或外部的存储设备可能包括一个边界条件库。用户或参数设置单元520可以根据对象特征数据自行设置一个边界条件或从所述边界条件库中选择一个已有的边界条件。在一些实施例中,用户或参数设置单元520也可以根据感兴趣的区域将一个低阶耦合模型作为边界条件。所述低阶耦合模型可以利用与感兴趣区域相耦合的区域或组织的经验性模型作为边界条件。所述低阶耦合模型可以是二阶模型,一阶模型,零阶模型(集中参数模型)或这些低阶模型的组合形式。
计算单元530可以对多时相特征生成模块430中其他单元中产生的数据或信息进行计算。在一些实施例中,计算单元530可以根据图像数据生成一个与之相对应的模型。所述模型的生成可以基于参数设置单元520所选择的模型类型和设置的参数。在一些实施例中,计算单元530可以在建立感兴趣区域的模型后,对模型进行分析和计算。所述分析和计算的方法可以包括计算流体力学(Computed fluid dynamics)等。在一些实施例中,所述分析和计算的结果可以包括血管、组织、或器官的物理状态和相关系数或参数。例如,对冠状动脉模型进行分析和计算的结果可以包括冠状动脉的血流动力学参数,如血流速率、血液压力、血管壁应力、血管壁切应力、血流储备系数(FFR)、冠状动脉血流储备(CFR)等中的一种或多种的组合。
在一些实施例中,计算单元530所计算的信息和数据可以是多时相的。计算单元530可以分别对所述多时相的信息和数据进行分析和计算。在一些实施例中,计算单元530可以根据不同时相的分析和计算结果生成感兴趣区域的物理状态和/或相关系数或参数与时相或时间的关系。在一些实施例中,该关系可以用曲线或者对照表的方式体现。基于所述曲线或对照表,可以获得任意时相的感兴趣区域的物理状态和/或相关系数或参数。在一些实施例中,所述曲线、对照表或者任意时相的感兴趣区域的物理状态和相关参数系数可以通过输出单元570发送到血流状态分析系统100其他模块或单元或血流状态分析系统100外界的模块或单元中。
网格生成单元540可以在生成模型的网格。在一些实施例中,网格生成单元540可以在模型上生成二维或三维的网格。例如,网格生成单元540可以在模型的边界区域(入口、出口等)生成二维的网格,而在模型的其他区域生成三维的网格。所述三维网格可以是基于所述二维网格而建立的。具体的关于网格形成的方法和流程可以参见,例如,图10、图12及其描述。
匹配单元550可以对多时相的数据进行匹配。在一些实施例中,匹配单元550可以将不同时相的模型进行相互关联。所述不同时相的模型可以是经过网格化处理后的模型。在一些实施例中,所述不同时相的模型相互关联的过程可以包括先识别出不同时相的模型中的特征区域。然后,将不同时相相对应的特征区域进行关联。例如,如果所述不同时相的模型是血流模型(即感兴趣的血管中血流覆盖区域的模型),所述特征区域可以包括血流入口区域、血流分叉区域、血流出口区域、血流狭窄区域、血流扩张区域等。然后,匹配单元550可以将不同时相的相应特征区域进行关联。在一些实施例中,一个特征区域在不同时相可能对应不同个数的网格。在这种情况下,所述特征区域在不同时相中的网格可以依靠某种算法或方法进行关联。例如,当第一时相中多个网格对应于第二时相中一个网格或少数网格时,匹配单元550可以对第一时相中多个网格的数值进行平均化处理后再与所述第二时相中一个或少数网格的数值进行对应。在一些实施例中,在初始时相的计算中,可以将内部网格(网格模型中不包括边界区域的网格)的初始值(例如,压强初始值,速率初始值等)设置为0。在后续时相的计算中,可以利用网格匹配将上一时相的内部网格计算结果映射或匹配至当前时相的内部网格相对应的网格中,并作为当前时相网格的初始值。在一些实施例中,匹配单元550在完成匹配后可以提示用户确认匹配是否准确。如果用户确认匹配准确则进行后续流程。如果用户认为匹配不准确,则用户可以对匹配结果进行修正或调整。用户也可以选择在其参与的情况下重新进行网格匹配。
区域选择单元560可以在图像数据中选择感兴趣的区域。所述感兴趣的区域可以由区域选择单元560自行选择或根据用户输入的信息选择。在一些实施例中,选择的感兴趣区域可以是血管、组织或者器官等。区域选择单元560可以进一步对所述图像中对感兴趣的区域进行分割。图像分割的方法可以包括基于边缘的图像分割方法,如Perwitt算子法、Sobel算子法、梯度算子法、Kirch算子法等,基于区域的图像分割方法,如区域生长法、阈值法、聚类法等以及其他分割方法,如基于模糊集、神经网络的方法等。区域选择单元560可以进全自动分割或半自动分割。例如,如果选择的感兴趣区域是冠状动脉、腹部动脉、大脑动脉、下肢动脉等,区域选择单元560可以自动进行分割。如果选择的感兴趣区域是其他机器较难准确分割的血管或部位,则可以进行半自动分割,由用户在分割过程中进行修正。在一些实施例中,区域选择单元560可以对根据图像数据重建的三维模型进行区域选择和分割。
输出单元570可以将多时相特征生成模块430中一个或多个单元所产生的信息、数据或处理结果发送至血流状态分析系统100中的其他模块或单元中。例如,输出单元570可以将计算单元530所生成的模型发送至交互设备140进行显示。又例如,输出单元570可以将网格生成单元540进行网格化处理后的模型发送至存储设备130中进行存储。
判断单元580可以进行逻辑判断。例如,血流状态分析系统100中其它模块或单元可以发送一个判断请求至判断单元580。判断单元580可以根据所述判断请求对相应内容进行判断。当判断出特定情况或生成判断结果后,判断单元580可以将判断结果或相应的操作指令发送至相应的模块或单元(例如,发送判断请求的模块或单元)。例如,判断单元580可以判断区域选择单元560所要分析的血管是否存在异常(如血管狭窄,动脉瘤等)。如果判断出该血管存在异常,判断单元560可以突出显示(如用不同的颜色表示)该异常血管,并同时提示用户确认该异常血管是否满足用户的要求。如果满足,则进行后续操作;如果不满足,则用户可以手动选择异常血管,然后进行后续操作。例如,区域选择单元560可以将用户选择的感兴趣区域和其自身生成的区域发送给判断单元580,判断单元580可以判断所述感兴趣区域和区域选择单元560生成的区域是否相同。如果判断单元580判断出用户和区域选择单元560生成的区域相同,则可以发送指令给区域选择单元560使其进行进一步地分割处理。否则,可以通过显示设备140再次确认用户的选择。
以上的描述仅仅是本发明的具体实施例,不应被视为是唯一的实施例。显然,对于本领域的专业人员来说,在了解本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变。例如,上述各单元以单时相为例进行了说明,但是可以理解的是各单元所接收、处理或输出的数据可以是多时相的。对于多时相的数据,上述各单元可以对不同时相的数据分别进行相应操作从而产生多时相的特征。例如,网格生成单元530可以对多时相的模型分别进行对应的网格化处理,从而生成多时相的网格化处理后的模型。又例如,参数设置单元520可以对多时相的模型或数据分别设置相对应的参数或边界条件。
根据本申请的一些实施例,图6是获取多时相特征的示例性流程图。在一些实施例中,获取多时相特征的流程可以通过多时相特征生成模块430实现。
在602中,可以接收多时相数据。所述多时相数据可以包括多时相图像数据、多时相对象特征数据等。多时相图像数据可以是包括对象的血管、组织或器官的多个不同时间点的图片或数据。所述对象特征数据可以包括对象的心率、心律、血压、血流速率、血液粘稠度、心输出量、心肌质量、血管流阻以及其他与血管、组织或器官相关的数据。如图7所示,702包含3个时相的心脏的图像。704是对象在一个心动周期内的血压曲线。在一些实施例中,所述多时相图像可以包含至少部分相同的血管、组织或器官。在一些实施例中,可以对获取的多时相图像进行预处理。所述预处理可以包括图像增强、图像降噪、图像平滑等。
在604中,可以选择多时相图像中感兴趣的血管区域。所述感兴趣的区域可以由区域选择单元560自行选择或根据用户输入数据选择。对于不同时相的图像,可以依次进行感兴趣区域的选择。不同时相的图像中选择的感兴趣区域可能是相同的。在一些实施例中,可以进一步对选择出的感兴趣的区域进行分割。图像分割的方法可以包括基于边缘的图像分割方法(如Perwitt算子法、Sobel算子法、梯度算子法、Kirch算子法等),基于区域的图像分割方法(如区域生长法、阈值法、聚类法等),以及其他分割方法,如基于模糊集、神经网络的方法等。所述分割的方法可以是全自动分割或半自动分割。例如,如果选择的感兴趣区域是冠状动脉、腹部动脉、大脑动脉、下肢动脉等,可以自动进行分割。如果选择的感兴趣区域是其他机器较难准确分割的血管或部位,则可以进行半自动分割,由用户在分割过程中进行修正。在一些实施例中,可以依次对不同时相的图像进行分割,也可以同时对不同时相的图像进行分割。
在606中,可以建立血管区域的多时相模型。所述多时相模型可以是血管模型,或血流模型。所述血管区域可以是冠状动脉血管区域、腹部动脉血管区域、大脑动脉血管区域、下肢动脉血管区域等。在一些实施例中,所述血管区域可以是上述血管的全部区域或者一部分区域。例如,所述血管区域可以是整个冠状动脉模型、左冠状动脉模型或右冠状动脉模型或者冠状动脉的分支的模型,如左前降支、左回旋支、对角支等。如图7所示,708中从左到右分别是图像、模型和分割后的模型。在一些实施例中,可以根据具体分析的病变部位(例如,感兴趣区域)以及对象特征数据(例如,血液粘稠度等)选择适合的血液粘性模型和流速边界模型。所述血液粘性模型可以包括牛顿流体模型、非牛顿流体模型以及用户自定义的其他流体模型。所述流速边界模型可以包括但不限于抛物线模型、双曲线模型、椭圆模型、平均流模型、Womersley分布模型、Reynolds模型、混合模型等。在一些实施例中,不同时相的图像可能会分别建立与其相对应的模型。
在608中,可以对建立的多时相模型进行网格化处理。在一些实施例中,可以在模型的边界区域(例如,血管的入口、出口等)生成二维的网格,而在模型中其他区域形成三维的网格。所述三维网格可以是基于所述二维网格而建立的。如图7所示,710是模型708进行网格化处理后的结果。具体的关于网格形成的方法和流程可以参见,例如,图10、图12及其描述。
在610中,可以设定多时相的参数和边界条件。在一些实施例中,所述参数的设置可以包括选定的模型中参数的设置,如流速u、密度ρ、血压P、截面积S等。所述边界条件的设置可以包括对边界区域进行初始条件或限制条件的设置。所述边界区域指的是感兴趣区域的边缘区域。例如,边界区域可以是血管的出口、入口、血管壁等。设置的边界条件可以包括边界区域的血压、血流量、血流速率、血管流阻、压强、应力等。在一些实施例中,可以根据感兴趣的区域将一个低阶耦合模型作为边界条件。所述低阶耦合模型可以利用与感兴趣区域相耦合的区域或组织的经验性模型作为边界条件。所述低阶耦合模型可以是二阶模型,一阶模型,零阶模型(集中参数模型)或这些低阶模型的组合形式。。如图7所示,712可以是模型的选择和参数的设置的一个实施例。714可以是边界条件的设定的一个实施例。
在612中,可以选择一个当前时相(在第一次选择的时候也可称作初始时相)。在一些实施例中,初始时相可以依据一些特定的规则进行选择出来。例如,对于冠状动脉模型,可以选择一个模型变化较慢或较平缓的时相作为初始时相,例如,离心脏收缩前期或者舒张末期最近的时相。初始时相可以由机器(例如,多时相特征生成模块430)自行判断选择,或由用户选择。如果机器和用户均没有或无法选择一个初始时相,可以任意选择一个时相或者选择多时相特征生成模块430第一个接收到的时相作为初始时相。
在614中,可以对当前时相(在第一次执行时为初始时相)进行分析,例如,计算流体力学(CFD)分析。根据设定的模型、边界条件以及参数可以对三维血管模型的血流动力学进行求解。所述求解的控制方程可以基于Euler方程组、纳维尔-斯托克斯方程组(Navier-Stokes Equations)或格子玻尔兹曼方程(Lattice Boltzmann Method)。求解计算的离散化方法可以包括有限差分法、有限体积法、有限元法、边界元方法、谱方法、格子玻尔兹曼法、无网格法等一种或多种组合。求解计算流场的流体可以是无粘的或有粘性的,可以是可压缩流体也可以是不可压缩流体,可以是层流也可以是湍流,可以是定常流动也可以是非定常流动。可以根据所模拟流体的物理特性,相应地选择相应的控制方程以及模拟方法。例如,针对无粘流体的流场计算可以选用Euler方程组或者格子玻尔兹曼方程,针对有粘性流体的流场计算可以选用N-S方程组或者玻尔兹曼方程。例如,对于冠状动脉进行的计算流体力学(CFD)计算可以使用纳维尔-斯托克斯方程(Navier-Stokes Equations):
和
其中,ρ代表血液密度,u代表血流的速度,t代表时间,以及σ代表血流应力项(由血压p和血液粘性力项决定)。在一些实施例中,在初始时相的计算中,可以将所述模型中的流速的初始值设置为0。在后续时相的计算中,所述模型中流速的初始值不再设置为0,而是匹配前后时相的网格,并根据上一时相的计算结果对当前时相的对应网格进行赋值,作为当前时相流速的初始值。
分析结果可以包括当前时相该模型任意一个区域或者位点的物理状态和相关系数或参数。例如,对冠状动脉模型进行分析的结果可以包括冠状动脉任意一个区域或者位点的血流动力学参数,如血流速率、血液压力、血管壁应力、血管壁切应力、血流储备系数(FFR)、冠状动脉血流储备(CFR)等中的一种或多种的组合。如图7所示,718显示一种血流动力学的分析计算。
在616中,可以判断是否已遍历所有时相。在判断出已遍历所有时相后,可以执行618。在判断出未遍历所有时相后,可以执行620。
在618中,可以将分析结果进行输出。例如,将分析结果发送至血流状态分析系统100其他模块或设备中。在一些实施例中,可以对分析结果进行后处理。所述后处理可以包括生成分析结果与时相的关系曲线或对照表。所述后处理也可以包括根据所述关系曲线或对照表进一步输出任意时相的分析结果的预估值。如图7所示,716显示一种后处理的结果。后处理的流程和结果的示例可以参考图13及其描述。在一些实施例中,618可以进一步包括将生成的多时相计算结果(例如,血管状态)与一个参考结果进行比对并生成一个比对结论。所述参考结果可以是存储在存储设备130中的数据,可以是存储在网络180中的数据,也可以是用户自行输入的数据。在一些实施例中,所述参考结果及相关比对结论可以存储在一个表中。例如,所述计算结果为血流速率时,所述参考结果可以为一个血流速率范围与危险程度的对应关系。所述危险程度可以分为正常、预警、危险、极度危险等。在一些实施例中,用户可以根据临床经验手动输入所述对应关系。在一些实施例中,所述比对可以是同一对象不同时期的血流速率的计算结果的比对。
在620中,可以将下一时相的模型与当前时相的模型进行匹配。在一些实施例中,所述不同时相的模型相互匹配的过程可以包括先识别出不同时相的模型中的特征区域。然后,将不同时相相对应的特征区域进行关联。在一些实施例中,所述特征区域可以包括血流入口区域、血流分叉区域、血流出口区域、血流狭窄区域、血流扩张区域等。操作620可以包括将不同时相的相应特征区域进行关联。如图7所示,706可以是模型以及其特征区域相互关联的一个实施例。在一些实施例中,一个特征区域在不同时相可能对应不同个数的网格。在这种情况下,所述特征区域在不同时相中的网格可以依靠某种算法或方法进行关联。例如,当第一时相中多个网格对应于第二时相中一个网格或少数网格时,匹配单元550可以对第一时相中多个网格的数值进行平均化处理后再与所述第二时相中一个或少数网格的数值进行对应。所述网格数值的对应包括将第一时相的网格的数值作为第二时相对应网格的输入。
在622中,可以将下一时相设为当前时相并进一步执行614。
以上的描述仅仅是本发明的具体实施例,不应被视为是唯一的实施例。显然,对于本领域的专业人员来说,在了解本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变。例如,606可以在604之前执行,可以先基于图像数据建立整个区域的模型,再在模型中选择和分割出感兴趣的部分。
根据本申请的一些实施例,图8是设定边界条件的示例性流程图。在一些实施例中,设定边界条件的流程与610相对应。在一些实施例中,设定边界条件的流程可以通过参数设置单元520实现。
在802中,可以获取血管区域的一个模型。在一些实施例中,所述模型可以通过606获得。所述模型可以是血管模型,也可以是血流模型。所述血管区域可以是冠状动脉血管区域、腹部动脉血管区域、大脑动脉血管区域、下肢动脉血管区域等。在一些实施例中,所述血管区域可以是上述血管的全部区域或者一部分区域。例如,所述血管区域可以是整个冠状动脉模型、左冠状动脉模型或右冠状动脉模型或者冠状动脉的分支的模型,如左前降支、左回旋支、对角支等。在一些实施例中,所述模型可以是用掩膜(mask)表示的血管区域。在一些实施例中,所述模型可以是用网格表示的血管模型。
在804中,可以判断所述血管区域的模型是否存在异常情况。所述异常情况可以包括血管狭窄、血栓、血管扩张、血管瘤等。如图9所示,模型910是一个存在狭窄的冠状动脉模型,该模型包括930、940两处狭窄区域。模型920是一个正常的冠状动脉模型,该模型与930、940相对应的区域(935、945)没有明显狭窄。在一些实施例中,804可以包括提取所述血管区域的模型的中心线。在一些实施例中,血管中心线可以指位于血管内部的沿着血管走向的一条假想的线。血管中心线可以包括血管中一个或多个像素点(或体素点)的集合。在一些实施例中,血管中心线可以包括血管中心或靠近血管中心的像素点(或体素点)的集合或组成的一条线。在一些实施例中,血管中心线可以包括一个或多个血管端点。所述血管中心线可以为所述端点之间的一条路径。在一些实施例中,关于血管中心线提取的示例性方法可以参考2016年8月30日提交的申请号为PCT/CN2016/097294的国际申请,其内容以引用的方式被包含于此。然后,可以沿所述中心线设置多个特征点。进一步地,可以计算模型在所述特征点的横截面积。根据特征点的横截面积,可以判断出模型是否存在异常情况。例如,如果计算出特征点的横截面积出现不正常的降低(例如,两个横截面积正常的特征点中间包括一个横截面积较低的特征点),可以判断模型可能存在狭窄或血栓。在一些实施例中,所述特征点的个数应该至少达到可以使得横截面积的变化足够判断是否存在异常情况的数量(例如,选择的相邻特征点的距离小于狭窄区域的长度)。在判断出模型不存在异常情况之后,可以执行822。在判断出模型存在异常情况之后,可以执行806。
在806中,可以确定存在异常情况的区域。例如,可以根据存在异常情况的模型(也称作异常模型)横截面积的变化情况,判断并标记模型中可能的狭窄或膨大区域。血管狭窄或膨大区域可能是某段血管的局部面积最小值或最大值,或血管截面变化较大的区域。在一些实施例中,可以将确定的异常区域发送给用户。如果发现异常区域不准确,用户可以对所述异常区域进行修正。例如,用户可以手动选择异常模型中的一个或多个点或者一段区间作为异常区域。
在808中,可以获取与异常模型相关的数据。所述数据可以包括血管入口血流流速、入口血流流量、入口血压、入口流阻、分流个数、出口个数、血液粘稠度、血液密度等。在一些实施例中,血管进口流量可以通过与之相连的组织或器官的相关参数或特征计算获得。例如,冠状动脉的进口流量可以用心输出量来估测,而心输出量可以通过分析心脏腔室在一个心动周期的容积变化获得。一些经验性的生理规律也可以用于估测这些物理量。例如,冠脉血流量与心肌质量成正比,Q∝QoMα,其中Q表示冠脉血流量,Qo为一常数,M表示心肌质量,指数α为一预定义的变化因子。在一些实施例中,心肌质量M可以通过无创的方法获得。例如,通过心肌体积与心肌密度相乘获得。在一些实施例中,冠状动脉的入口血压可以通过血压计等测量获得。
在810中,可以基于异常模型重建一个正常模型。在一些实施例中,重建的区域仅限于确定的异常区域,而模型中其它区域保持不变。在一些实施例中,重建的方法可以是基于血管异常区域前后的管径和其中心线进行放样或拉伸操作从而形成正常区域。在一些实施例中,重建的方法也可以是对异常区域进行膨胀或缩小操作。所述膨胀操作后的血管截面不大于狭窄段前后血管的截面。所述缩小操作后的血管截面不小于扩张段前后血管的截面。在一些实施例中,可以对重建后的异常区域周围进行平滑处理从而避免其出现明显的突变。在一些实施例中,可以将重建的正常模型发送给用户。如果发现重建不准确,用户可以对部分或全部重建的正常模型进行修正。例如,用户可以对重建后的正常模型进行局部膨胀、局部腐蚀,局部平滑等。
在812中,可以获取与正常模型相对应的数据。所述数据可以包括正常模型的边界条件、各出入口流阻等。边界条件可以包括入口和出口的血压、流速、流量等。在一些实施例中,入口的血压、流速可以在808中获得。出口的流量可以通过计算获得。例如,可以假设分支血管的流量分配与分支管径正相关,即Q∝dk,其中d为分支血管近端(靠近分支点)的平均管径,k为放大系数。然后,可以从进口开始在每个血管分叉口按照上述正相关的关系进行流量分配,直到将所述进口的流量分配到每个出口为止。基于正常模型以及所述边界条件,可以进行计算流体力学(CFD)模拟计算,从而获得正常模型每个出口的流阻(出口压强与出口流量的比值)。
在816中,可以确定模型的总流阻。总流阻可以通过以下公式计算:
其中R表示模型的总流阻,Pinlet表示入口血压压强,以及Q为入口血流量。Pinlet和Q可以在808中获得,在此不再赘述。
在818中,可以通过对血管中心线、血管截面以及血管异常区域的分析来确定对象血管在正常情况下各级各分支血管近端(靠近分支点)的管径大小。
在820中,可以根据正常模型管径的大小对出口的流阻进行分配。流阻可以依照以下公式进行分配:
其中d表示管径,i表示本级分叉血管的编号,j表示本次流阻分配所处级别,以及k表示流阻分配指数(比如,冠状动脉中k可以被设置为2.7)。
在822中,可以生成与实际血管模型相对应的边界流阻。在一些实施例中,所述与实际血管模型相对应的边界流阻可以根据正常模型相对应的边界流阻获得。例如,实际血管模型相对应的边界流阻可以与正常模型相对应的边界流阻相同。在一些实施例中,血管模型的边界流阻可以根据步骤820所描述的流阻分配方式获得。
根据本申请的一些实施例,图10是网格化处理的示例性流程图;在一些实施例中,网格化处理的流程与608相对应。所述网格化处理的流程可以通过网格生成单元540实现。
在1002中,可以获取一个模型。所述模型可以是本申请其他实施例中说明的模型,例如,对象血管或血流、组织或器官或其他感兴趣区域重建后的模型。如图11所示,1110可以是一个冠脉血流模型,即模型1110可以表示冠脉血管中血流所覆盖的区域。在不考虑血管壁厚度以及血管堵塞等情况下,模型1110也可以近似表示一个冠脉血管模型。
在1004中,可以确定所述模型的边界区域。如果模型是一个血管或与血管对应的血流区域,边界区域可以是血管的出口、入口、血管壁等。如图11所示,模型1110的入口1120可以在1004中被确定为模型1110的边界区域。
在1006中,可以对确定的边界区域进行面网格划分(也被称为二维网格划分)。所述面网格划分可以是将边界区域对应的平面用网格进行划分。网格的划分算法包括三角网格划分,四角网格划分,六角网格划分,或者类似的,或者一个或多个的组合。示例性的网格划分算法包括Loop算法,蝶型细分算法,Catmull-Clark算法,Doo-Sabin算法,Delaunay三角划分算法等。网格划分的方法示例可以参见图12及其描述。如图11所示,1130是模型1110的入口截面图,以及1140是1130网格划分的示例性结果。
在1008中,可以对模型侧壁进行面网格划分。在一些实施例中,侧壁和边界区域可以使用不同的网格划分方法进行划分。例如,侧壁可以使用曲面网格剖分算法进行网格划分。所述曲面网格剖分的算法可以包括映射法和自动网格生成法等。所述映射法可以包括将侧壁映射至平面,利用二维网格划分方法对平面进行划分后,再将划分后的网格映射回侧壁。所述自动网格生成法可以根据侧壁中不同区域的曲度将侧壁分成若干个近似平面再分别进行二维网格划分。平面网格划分可以参照本申请中其他地方的描述,例如,图12以及描述。
在1010中,可以根据边界区域和侧壁的面网格划分结果对模型进行体网格划分(也被称为三维网格划分)。所述体网格划分可以是将模型用三维的网格进行划分。所述三维网格可以包括四面体网格、六面体网格、棱柱体网格(边界层网格)、四面体与六面体混合网格、笛卡尔网格、球填充法网格等。在一些实施例中,1004至1008可以被跳过,即可以直接对模型进行体网格划分。
根据本申请的一些实施例,图12是网格划分的示例性流程图。在一些实施例中,流程1200可以通过多时相特征生成模块430实现。在一些实施例中,图6中的608,图10中的1006等可以根据流程1200实施。
在1202中,可以获取一个二维图像。在一些实施例中,所述二维图像可以通过数据获取单元510获得。在一些实施例中,二维图像可以是一个二维医学图像,或者其中用户感兴趣的部分(例如,冠脉血管所在的区域,大脑区域等)。仅仅作为示例,二维图像可以是一个CT图像、MRI图像、PET图像等。二维图像可以以灰度或者彩色的方式呈现出来。在一些实施例中,二维图像可以是一个时相模型的二维显示。例如,二维图像可以是一个与1006中的时相模型的边界区域有关的图像。更具体地,二维图像可以显示一个血流模型的入口/出口区域(例如,如图9所示)。二维图像可以是一个由图像处理设备(例如,处理设备120)重建的图像。二维图像可以来自一个本地存储装置或者外界的存储装置(例如,存储设备130)。
在1204中,网格生成单元540可以提取二维图像中感兴趣区域的轮廓点。在一些实施例中,提取二维图像中感兴趣区域的轮廓点包括先分割感兴趣区域,再对分割后的感兴趣区域提取轮廓点。所述对感兴趣区域的分割方法见本申请中其他地方的描述。在一些实施例中,感兴趣区域的轮廓点可以包括一个或多个位于感兴趣区域边界的像素点(也称为“轮廓像素点”)。例如,冠脉入口截面处的轮廓点可以包括一个或多个位于冠脉壁的轮廓像素点。在一些实施例中,感兴趣区域的轮廓像素点可以是连续的,部分连续的,或者不连续的。这里所说的连续是指一个轮廓像素点与至少一个或多个其他轮廓像素点相临。提取出的轮廓点的信息可以保存在一个或者多个存储装置中(例如,存储设备130,存储模块260等)。轮廓点的信息可以在后续流程中被网格生成单元540或者其他可以进行数据分析的单元/模块使用。示例性的轮廓点的信息可以包括轮廓点的位置,轮廓点的数量,或者类似的,或者一个或多个的组合。
在1206中,可以根据轮廓点确定一个或多个区域。在一些实施例中,所述一个或多个区域的确定可以通过网格生成单元540实现。所述一个或多个区域可以通过依次连接感兴趣区域的轮廓点形成。仅仅作为示例,确定一个或多个区域可以包括确定一个感兴趣区域的初始轮廓像素点(例如,初始轮廓像素点可以选为轮廓像素点中x/y坐标最小的点)。按照顺时针或者逆时针方向,将感兴趣区域的轮廓像素点进行排序。从初始轮廓像素点开始,用线条连接前一个轮廓像素点与后一个轮廓像素点以形成短边。当最后一个轮廓像素点与初始轮廓像素点连接形成短边后,可以形成一个封闭的轮廓曲线。在一些实施例中,感兴趣区域可以位于一条封闭的轮廓曲线的内部。例如,图11中模型入口1140的感兴趣区域(即,进行网格划分的区域)位于轮廓曲线的内部。在一些实施例中,感兴趣区域可以是位于两条封闭的轮廓曲线之间的区域。例如,感兴趣区域可以是一个二维环状结构,或者与二维环状结构拓扑等价的结构。一个或多个区域的信息(例如,区域对应的网格曲线)可以保存在一个或者多个存储装置中(例如,存储设备130,存储模块260等)。一个或多个区域的信息可以在后续流程中被网格生成单元540或者其他可以进行数据分析的单元/模块使用。
在1208中,可以对一个区域是否需要进行网格划分进行判断。在一些实施例中,所述判断可以由判断单元580实现。在判定出该区域不需要进行网格划分之后,流程1200进入1210。在判定出该区域需要进行网格划分之后,流程1200进入1212。在一些实施例中,判断单元580进行判断的条件可以是该区域是否是感兴趣区域。当该区域是感兴趣区域时,则被判定为需要进行网格划分。如本披露书中在其他地方描述的,感兴趣区域可以包括需要进行血液状态分析的区域,例如,特定血管中血液流动的区域。
在1210中,可以对该区域进行标记。在一些实施例中,对区域的标记可以由网格生成单元540实施。标记可以采用计算机可读取的代码或者可执行的指令的形式。被标记的区域信息可以被保存在一个或多个存储装置(例如,存储设备130,存储模块260等)。被标记的区域信息可以在后续流程中被网格生成单元540或者其他可以进行数据分析的单元/模块读取,从而在进行网格划分时剔除该被标记的区域。
在1212中,可以对该区域进行网格划分。在一些实施例中,网格划分可以由网格生成单元540实施。在一些实施例中,网格的划分可以基于该区域的轮廓点进行。网格的划分算法包括三角网格划分,四角网格划分,六角网格划分,或者类似的,或者一个或多个的组合。示例性的网格划分算法包括Loop算法,蝶型细分算法,Catmull-Clark算法,Doo-Sabin算法,Delaunay三角划分算法等。仅仅作为举例,网格生成单元540可以利用Delaunay三角划分(Delaunay triangulation)算法对该区域的所有轮廓点进行网格划分。又例如,网格生成单元540可以先将该区域的轮廓点分成不同子集,并对各个子集轮廓点进行网格划分。然后,网格生成单元540可以将各个子集的网格划分合并形成该区域的网格划分。具体地,可以将该区域所有的轮廓点按照x/y坐标进行排序(例如,先按照x坐标进行非递减排序,对于x坐标相同的点,按照y坐标非递减排序)。将排序后的轮廓点按照数量分成子集A和子集B,并分别完成Delaunay三角划分。再将子集A和子集B的Delaunay三角划分合并成所有轮廓点的Delaunay三角划分。在一些实施例中,网格划分还可以包括将该区域的轮廓曲线与通过划分算法划分的网格进行叠加,从而划分的网格中保留该区域的轮廓曲线(例如,1206中提到的由轮廓像素点连接形成的一个或多个短边)。
在一些实施例中,对一个区域的网格的划分可以采用并行技术的网格生成方法。例如,应用区域分裂或类似算法将该区域分为若干个子区域,在每个子区域内独立进行网格划分,再修复相邻子区域的边界网格,进而获得该区域的完整网格。
在1214中,可以对该区域设置网格划分控制条件。在一些实施例中,网格划分控制条件的设置可以由网格生成单元540实现。网格划分控制条件可以控制网格的数量,大小,分布,形状,或者类似的,或者一个或多个的组合。在一些实施例中,网格生成单元540可以设置一个网格单元的面积约束条件,使得任意网格单元的面积满足该面积约束条件。例如,网格生成单元540可以设置一个面积约束值,使得任意网格单元的面积都不大于该面积约束值。在一些实施例中,网格生成单元540可以设置一个网格单元的内角约束条件,使得任意网格单元的内角都满足该内角约束条件。例如,网格生成单元540可以设置一个三角网格单元的内角约束值,使得任意三角网格单元的最小内角都不小于该内角约束值。在一些实施例中,网格划分控制条件可以由用户通过,例如,交互设备140,输入后获得。网格划分控制条件也可以由网格生成单元540或者其他具有数据分析功能的单元/模块根据特定条件分析获得。特定条件可以包括生成网格需要的时间,生成的网格的数量,根据生成的网格进行模型计算的时间,根据生成的网格计算获得的结果的精确程度等。
在1216中,可以判断划分的网格是否满足控制条件。在一些实施例中,对网格划分的判断可以由网格生成单元540完成。如果划分的网格不满足控制条件,则流程1200进入1218。
在1218中,可以对网格进行处理。在一些实施例中,对网格的处理可以由网格生成单元540完成。网格处理可以包括调整网格数量,改变网格大小等一个或多个操作。调整网格数量可以包括提高网格密度,降低网格密度等。改变网格大小可以包括分割网格,合并网格,对网格进行重组等。
在一些实施例中,如果一个三角网格单元不满足面积约束条件(例如,该三角网格单元的面积大于面积约束值),则可以在该三角网格单元中插入一个或多个辅助点。辅助点的插入可以是随机的,也可以按照原三角网格的特征点的位置插入。网格生成单元540可以根据辅助点生成新的网格。例如,可以在三角网格单元的内部,例如,重心处插入一个辅助点。连接辅助点和原三角网格的顶点可以生成三个新的三角网格单元。又例如,可以在三角网格单元内部随机性地或者与非随机性地插入多个辅助点。根据多个辅助点利用Delaunay三角划分(Delaunay triangulation)算法划分出Delaunay三角网格。在一些实施例中,如果一个三角网格单元不满足内角约束条件,则可以采用特定的算法对三角网格单元进行处理。例如,可以采用flip算法更新三角网格单元。具体地,flip算法可以包括通过选择一个包含两个相邻的三角网格单元的四边形(四边形的一条对角线为两个三角网格单元相邻的边),选择该四边形中另一条对角线作为两个新三角形的边),获得两个新的三角网格单元。内角约束条件可以包括三角网格单元的最小内角不小于一个内角约束值。内角约束值可以是5°、10°、15°、20°、25°,或者其他数值。
网格经过处理后重新返回1216,由网格生成单元540判断处理后的网格是否满足控制条件。直到网格满足控制条件,则流程1200进入1220。
在1220,网格生成单元540可以判断是否已经遍历所有区域。如果否,则流程1200返回到1208,对没有处理到的区域进行是否划分网格的判断。如果已经遍历所有区域,则在1222中,由网格生成单元540生成感兴趣区域的网格。在一些实施例中,划分不同区域的网格类型的算法可以相同,或不同。例如,所有区域可以都采用Delaunay三角划分算法划分网格。又例如,一部分区域可以采用Delaunay三角划分算法划分网格,一部分区域可以采用四角网格划分算法或者六角网格划分算法进行网格划分。在一些实施例中,不同区域的网格划分控制条件可以相同,或不同。例如,所有区域的网格划分控制条件可以包括面积约束条件和/或内角约束条件。不同区域的面积约束条件和/或内角约束条件可以相同,或不同。具体地,所有区域的内角约束条件可以是任意三角网格的最小内角不小于一个内角约束值(例如,20°)。又例如,脑部图像的面积约束条件可以包括最大三角网格单元的面积不大于A,血管图像的面积约束条件可以包括最大三角网格单元的面积不大于B,其中A小于B。
以上的描述仅仅是本发明的具体实施例,不应被视为是唯一的实施例。显然,对于本领域的专业人员来说,在了解本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变。在一些实施例中,1210可以被跳过。例如,如果一个区域被判断为不需要进行网格划分,则可以直接删除该区域。在一些实施例中,1214可以移到1208之前,即网格生成单元540可以对所有需要划分网格的区域设置相同的网格划分控制条件。在一些实施例中,流程1200可以对三维图像进行网格划分。例如,对三维区域的网格划分可以采用快速Delaunay的球充填法(ball-packing)。网格生成单元540可以在三维几何区域内用球充填法布点,依据几何模型的几何特征和空间关系自适应地分布疏密合适的节点,然后采用快速Delaunay插入技术生成三维网格。
根据本申请的一些实施例,图13是描述获得一个位点对应的血流动力学参数的示例性流程图。在一些实施例中,获得一个位点对应的血流动力学参数的流程可以通过多时相特征处理模块440实现。在一些实施例中,图4B中的468可以根据获得一个位点对应的血流动力学参数的流程实施。
在1302中,多时相特征处理模块440可以获得多时相的血流动力学参数值。在一些实施例中,多时相的血流动力学参数值可以与所述获取多时相特征的流程中614-618有关。所述血流动力学参数值可以表示冠状动脉血管区域、腹部动脉血管区域、大脑动脉血管区域、下肢动脉血管区域等血管区域的血流状态。血流动力学参数可以包括血液流速、血液压力、血管壁应力、血管壁切应力、血流储备系数(FFR)、冠状动脉血流储备(CFR)等中的一种或多种的组合。在一些实施例中,多时相的血液动力学参数值可以对应一个特定时间段内的血流状态。例如,可以获取一个心脏跳动周期内不同时相的血流动力学参数值,从而获得在该心脏跳动周期内的血流状态。获取的时相数量可以是3个,或者5个,或者8个,或者10个,或者15个等。
在1304中,多时相特征处理模块440可以确定一个血管中的位点。所述位点可以是血管入口/出口的平面上的任意一点,或血管壁或者血管内部空间上的任意一点。在一些实施例中,所述位点的确定可以由用户通过,例如,交互设备140,输入后实现。
在1306中,多时相特征处理模块440可以拟合位点的血流动力学参数曲线。在一些实施例中,所述血流动力学参数曲线可以表示一个心脏跳动周期内的血流状态。所述拟合可以包括将多时相的血流动力学参数按照一个函数进行拟合。用于拟合的函数可以是线性的,或非线性的。适用的非线性函数可以包括多项式函数,对数函数,指数函数,或者类似的,或者其中多个的组合。例如,根据冠脉血管入口平面上一个位点多时相的FFR值,可以拟合出在一定时间范围内所述位点的FFR曲线。进一步的,在获得所述位点在一个心脏跳动周期内的FFR拟合曲线后,可以根据心脏跳动的周期性,生成所述为点在任意时刻的FFR曲线。
在1308中,多时相特征处理模块440可以根据参数曲线获得位点一个感兴趣时相的血流动力学参数值(例如,FFR值)。所述感兴趣时相可以不同于1302中获得的多时相中任意一个。在一些实施例中,所述感兴趣时相的选择可以是由用户通过,例如,交互设备140,输入后实现。在一些实施例中,多时相特征处理模块440可以根据血流动力学参数曲线对所述位点的血流动力学参数值进行处理。例如,可以根据一段时间的血流动力学参数值(例如,FFR值),获取平均血流动力学参数值(例如,平均FFR值)。
[1]以上的描述仅仅是本发明的具体实施例,不应被视为是唯一的实施例。显然,对于本领域的专业人员来说,在了解本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修正和改变。在一些实施例中,在进行拟合位点的血流动力学参数曲线之前,多时相特征处理模块440可以获得为点额外的血流动力学参数。所述额外的血流动力学参数可以是通过插值方法处理后获得,也可以是用户通过,交互设备140,输入后获得。
根据本申请的一些实施例,图14是描述获得位点对应的血流动力学参数的示意图。图片1402显示包含心脏区域和腹部区域的多时相图像数据。图片1406显示一个时相的图像中对应的冠脉以及建立的冠脉模型。图像1408显示所述冠脉在不同时相的FFR分布。所述不同时相的血液流动状态可以根据获取多时相特征的流程获得。图片1404显示对象的特异性临床数据,包括对象的主动脉压随时间的变化曲线,相位冠脉血流量随时间的变化曲线。图片1410显示对象的FFR随时间变化的曲线。
上文已对基本概念做了描述,显然,对于本领域技术人员来说,上述发明披露仅仅作为示例,而并不构成对本申请的限定。虽然此处并没有明确说明,本领域技术人员可能会对本申请进行各种修改、改进和修正。该类修改、改进和修正在本申请中被建议,所以该类修改、改进、修正仍属于本申请示范实施例的精神和范围。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
Claims (10)
1.一种血流状态分析方法,其特征在于:所述方法包括:
获取一个第一时相的图像和一个第二时相的图像;
在所述第一时相的图像中选择一个第一血管区域,其中,所述第一血管区域包含一个血管;
在所述第二时相的图像中选择一个第二血管区域,其中,所述第二血管区域包含至少一部分所述血管;
建立一个第一血管模型,其中,所述第一血管模型与所述第一血管区域相对应;
建立一个第二血管模型,其中,所述第二血管模型与所述第二血管区域相对应;
设定所述第一血管模型的边界条件和所述第二血管模型的边界条件;
根据所述第一血管模型的边界条件,确定所述第一血管模型中所述血管在所述第一时相的状态;
基于所述血管在所述第一时相的状态,关联所述第一血管模型和所述第二血管模型;以及
根据所述关联结果以及所述第二血管模型的边界条件,确定所述第二血管模型中所述血管在所述第二时相的状态。
2.如权利要求1所述的方法,其特征在于:所述第一血管区域和所述第二血管区域包括冠状动脉、腹部动脉、大脑动脉或下肢动脉。
3.如权利要求1所述的方法,其特征在于:所述关联所述第一血管模型和所述第二血管模型包括关联所述第一血管模型和所述第二血管模型的所述血管的入口段、所述血管的分叉段、所述血管的狭窄段、或所述血管的出口段。
4.如权利要求1所述的方法,其特征在于:所述方法进一步包括对所述第一血管模型得到第一血管模型的网格和对所述第二血管模型进行网格化处理后得到第二血管模型的网格,所述关联所述第一血管模型和所述第二血管模型包括将第一血管模型的网格和所述第二血管模型的网格进行匹配。
5.如权利要求1所述的方法,其特征在于:所述血管在所述第一时相的状态或在第二时相的状态分别包括血流速、血压力、血管壁应力、血管壁切应力或血流储备系数(FFR)。
6.如权利要求1所述的方法,其特征在于:所述确定所述第一血管模型中所述血管在所述第一时相的状态或所述确定所述第二血管模型中所述血管在所述第二时相的状态包括使用计算流体力学(CFD)方法。
7.一个血流状态分析系统,包括至少一个处理器和存储设备,其特征在于:所述系统包括:
一个接收模块,被配置为:
获取一个第一时相的图像和一个第二时相的图像;以及
一个多时相特征生成模块,被配置为:
在所述第一时相的图像中选择一个第一血管区域,其中,所述第一血管区域包含一个血管;
在所述第二时相的图像中选择一个第二血管区域,其中,所述第二血管区域包含至少一部分所述血管;
建立一个第一血管模型,其中,所述第一血管模型与所述第一血管区域相对应;
建立一个第二血管模型,其中,所述第二血管模型与所述第二血管区域相对应;
设定所述第一血管模型的边界条件和所述第二血管模型的边界条件;
根据所述第一血管模型的边界条件,确定所述第一血管模型中所述血管在所述第一时相的状态;
基于所述血管在所述第一时相的状态,关联所述第一血管模型和所述第二血管模型;以及
根据所述关联结果以及所述第二血管模型的边界条件,确定所述第二血管模型中所述血管在所述第二时相的状态。
8.一种血流状态分析方法,所述分析方法包括:
获取多个时相的血管图像,包括一个第一时相的血管图像和一个第二时相的血管图像,所述多个时相的血管图像分别对应于同一个血管或其一部分;
建立多个血管模型,所述多个血管模型分别对应于所述多个时相的血管图像;
根据所述多个血管模型,获得所述血管或其一部分的多个状态,包括一个第一血管状态和一个第二血管状态,其中,所述第一血管状态对应于所述第一时相的血管图像,所述第二血管状态对应于所述第二时相的血管图像;
根据所述血管或其一部分的所述多个状态,获得所述血管或其一部分的状态随时间变化的关系;以及
根据所述关系,获得所述血管或其一部分的第三血管状态。
9.如权利要求8所述的方法,其特征在于:所述方法进一步包括:
关联所述多个血管模型;以及
根据所述关联结果,采用计算流体力学(CFD)分析所述多个血管状态,所述关联所述多个血管模型包括关联所述多个血管模型中至少两个血管模型的入口段、分叉段、狭窄段或出口段。
10.如权利要求8所述的方法,其特征在于:所述方法进一步包括对所述多个血管模型分别进行网格化处理得到多个血管模型的网格,以及将所述多个血管模型的网格进行匹配。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNPCT/CN2017/072256 | 2017-01-23 | ||
PCT/CN2017/072256 WO2018133118A1 (zh) | 2017-01-23 | 2017-01-23 | 血流状态分析系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107123112A CN107123112A (zh) | 2017-09-01 |
CN107123112B true CN107123112B (zh) | 2020-08-04 |
Family
ID=59724862
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710257724.8A Active CN107123112B (zh) | 2017-01-23 | 2017-04-19 | 血流状态分析系统及方法 |
CN201710257262.XA Active CN107115111B (zh) | 2017-01-23 | 2017-04-19 | 血流状态分析系统及方法 |
CN201710258086.1A Pending CN107123159A (zh) | 2017-01-23 | 2017-04-19 | 血流状态分析系统及方法 |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710257262.XA Active CN107115111B (zh) | 2017-01-23 | 2017-04-19 | 血流状态分析系统及方法 |
CN201710258086.1A Pending CN107123159A (zh) | 2017-01-23 | 2017-04-19 | 血流状态分析系统及方法 |
Country Status (4)
Country | Link |
---|---|
US (4) | US10360682B2 (zh) |
EP (1) | EP3375364A4 (zh) |
CN (3) | CN107123112B (zh) |
WO (1) | WO2018133118A1 (zh) |
Families Citing this family (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10210956B2 (en) * | 2012-10-24 | 2019-02-19 | Cathworks Ltd. | Diagnostically useful results in real time |
US11071501B2 (en) | 2015-08-14 | 2021-07-27 | Elucid Bioiwaging Inc. | Quantitative imaging for determining time to adverse event (TTE) |
US11094058B2 (en) | 2015-08-14 | 2021-08-17 | Elucid Bioimaging Inc. | Systems and method for computer-aided phenotyping (CAP) using radiologic images |
US11676359B2 (en) | 2015-08-14 | 2023-06-13 | Elucid Bioimaging Inc. | Non-invasive quantitative imaging biomarkers of atherosclerotic plaque biology |
US11113812B2 (en) | 2015-08-14 | 2021-09-07 | Elucid Bioimaging Inc. | Quantitative imaging for detecting vulnerable plaque |
US10176408B2 (en) | 2015-08-14 | 2019-01-08 | Elucid Bioimaging Inc. | Systems and methods for analyzing pathologies utilizing quantitative imaging |
US11087459B2 (en) | 2015-08-14 | 2021-08-10 | Elucid Bioimaging Inc. | Quantitative imaging for fractional flow reserve (FFR) |
JP6653673B2 (ja) * | 2017-02-28 | 2020-02-26 | 富士フイルム株式会社 | 血流解析装置および方法並びにプログラム |
CN107610773A (zh) * | 2017-09-12 | 2018-01-19 | 北京即刻叁维数据科技股份有限公司 | 一种基于主动脉医学影像的血管夹层辅助诊断方法 |
CN108186038B (zh) * | 2018-02-11 | 2020-11-17 | 杭州脉流科技有限公司 | 基于动脉造影影像计算冠脉血流储备分数的系统 |
CN108257214A (zh) * | 2018-02-12 | 2018-07-06 | 成都信息工程大学 | 心脏左心室流场域数值模拟方法、计算机、计算机程序 |
CN109065170B (zh) * | 2018-06-20 | 2021-11-19 | 博动医学影像科技(上海)有限公司 | 获取血管压力差的方法及装置 |
CN109102893A (zh) * | 2018-07-04 | 2018-12-28 | 中山大学 | 一种基于Cross模型的多粒子混合修正的血栓模拟方法 |
CN110599444B (zh) * | 2018-08-23 | 2022-04-19 | 深圳科亚医疗科技有限公司 | 预测血管树的血流储备分数的设备、系统以及非暂时性可读存储介质 |
CN112584738B (zh) * | 2018-08-30 | 2024-04-23 | 奥林巴斯株式会社 | 记录装置、图像观察装置、观察系统、观察系统的控制方法及存储介质 |
CN110428420B (zh) * | 2018-09-05 | 2022-05-17 | 深圳科亚医疗科技有限公司 | 基于患者的冠状动脉ct血管造影图像来确定冠状动脉的流动信息的方法、装置和介质 |
CN110384493A (zh) * | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | 测量微循环阻力指数的系统以及冠脉分析系统 |
CN112384138B (zh) * | 2018-10-25 | 2023-07-25 | 苏州润迈德医疗科技有限公司 | 获取心表大动脉的血流量的方法、装置、系统及存储介质 |
CN109616200A (zh) * | 2018-11-06 | 2019-04-12 | 北京三普威盛科技有限公司 | 用于冠脉狭窄评估的方法,装置,存储介质及电子设备 |
CN112022137B (zh) | 2018-11-30 | 2021-07-13 | 博动医学影像科技(上海)有限公司 | 建立血管截面函数和血管应力的方法及装置 |
US11819279B2 (en) * | 2018-11-30 | 2023-11-21 | Koninklijke Philips N.V. | Patient lumen system monitoring |
CN113520353B (zh) * | 2018-12-31 | 2023-06-16 | 深圳北芯生命科技股份有限公司 | 导管模拟器的使用方法 |
CN109948622B (zh) * | 2019-03-26 | 2020-12-11 | 数坤(北京)网络科技有限公司 | 一种头颈体动脉瘤检测方法、装置及计算机可读存储介质 |
US10861157B2 (en) | 2019-04-04 | 2020-12-08 | Medtronic Vascular, Inc. | System and methods for determining modified fractional flow reserve values |
CN110400638B (zh) * | 2019-05-30 | 2022-09-30 | 合肥工业大学 | 用于研究颅内动脉瘤生长过程以及探究生长机制的方法 |
JP2022543330A (ja) | 2019-08-05 | 2022-10-12 | エルシド バイオイメージング インコーポレイテッド | 形態学的および血管周囲疾患の複合評価 |
CN110910398B (zh) * | 2019-10-28 | 2021-07-20 | 衢州学院 | 一种基于决策层融合的视频复杂场景区域分割方法和装置 |
US11200669B2 (en) * | 2019-11-19 | 2021-12-14 | Uih America, Inc. | Systems and methods for determining plasma input function used in positron emission tomography imaging |
CN111161342B (zh) * | 2019-12-09 | 2023-08-29 | 杭州脉流科技有限公司 | 基于冠脉造影图像获取血流储备分数的方法、装置、设备、系统及可读存储介质 |
CN113487588B (zh) * | 2020-03-02 | 2024-03-22 | 上海博动医疗科技股份有限公司 | 冠状动脉的血流速度的计算方法、装置及电子设备 |
WO2022001026A1 (en) * | 2020-06-30 | 2022-01-06 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining blood vessel parameters |
CN112001893B (zh) * | 2020-07-31 | 2024-02-20 | 上海联影医疗科技股份有限公司 | 一种血管参数的计算方法、装置、设备及存储介质 |
US11439308B2 (en) * | 2020-07-13 | 2022-09-13 | GE Precision Healthcare LLC | Methods and systems for thermal monitoring of tissue with an ultrasound imaging system |
US20220044408A1 (en) * | 2020-08-07 | 2022-02-10 | Canon Medical Systems Corporation | Medical image processing apparatus, system, and method |
CN112561871B (zh) * | 2020-12-08 | 2021-09-03 | 中国医学科学院北京协和医院 | 一种基于平扫ct图像的主动脉夹层分割方法和装置 |
CN112434676B (zh) * | 2021-01-26 | 2021-04-20 | 北京圣点云信息技术有限公司 | 一种基于汇流累计量信息的体内静脉识别方法及装置 |
US20230079079A1 (en) * | 2021-09-14 | 2023-03-16 | United Therapeutics Corporation | Systems and methods for tissue evaluation |
CN113963001B (zh) * | 2021-10-25 | 2023-05-23 | 中国人民解放军总医院第一医学中心 | 一种基于多时相cta的侧支循环血管定量方法 |
CN113706539A (zh) * | 2021-10-29 | 2021-11-26 | 南京裕隆生物医学发展有限公司 | 一种用于识别肿瘤的人工智能辅助系统 |
WO2023152688A1 (en) * | 2022-02-10 | 2023-08-17 | Cathworks Ltd. | System and method for machine-learning based sensor analysis and vascular tree segmentation |
EP4311496A1 (en) * | 2022-07-26 | 2024-01-31 | Koninklijke Philips N.V. | Blood flow parameters |
CN116636826B (zh) * | 2023-07-27 | 2023-09-22 | 深圳市爱保护科技有限公司 | 血压估算模型的训练方法、装置、电子设备及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101347343A (zh) * | 2007-07-17 | 2009-01-21 | 株式会社东芝 | 超声波图像取得装置 |
CN101527047A (zh) * | 2008-03-05 | 2009-09-09 | 深圳迈瑞生物医疗电子股份有限公司 | 使用超声图像检测组织边界的方法与装置 |
CN102274046A (zh) * | 2010-06-08 | 2011-12-14 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置以及医用图像诊断装置 |
CN104720851A (zh) * | 2013-12-23 | 2015-06-24 | Ge医疗系统环球技术有限公司 | 计算机断层扫描成像的方法和装置 |
CN104720894A (zh) * | 2015-02-11 | 2015-06-24 | 中山大学附属第一医院 | 一种血管手术方式的合理性分析方法 |
CN105096388A (zh) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | 基于计算流体力学的冠状动脉血流仿真系统和方法 |
CN105190630A (zh) * | 2013-01-15 | 2015-12-23 | 凯瑟沃克斯有限公司 | 计算血流储备分数 |
CN106327487A (zh) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
CN106539622A (zh) * | 2017-01-28 | 2017-03-29 | 北京欣方悦医疗科技有限公司 | 基于血流动力学分析的冠脉虚拟支架植入方法和系统 |
Family Cites Families (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011046511A1 (en) * | 2009-10-13 | 2011-04-21 | Agency For Science, Technology And Research | A method and system for segmenting a liver object in an image |
US8682626B2 (en) | 2010-07-21 | 2014-03-25 | Siemens Aktiengesellschaft | Method and system for comprehensive patient-specific modeling of the heart |
US8315812B2 (en) * | 2010-08-12 | 2012-11-20 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
US9119540B2 (en) | 2010-09-16 | 2015-09-01 | Siemens Aktiengesellschaft | Method and system for non-invasive assessment of coronary artery disease |
US9141763B2 (en) | 2011-02-07 | 2015-09-22 | Siemens Aktiengesellschaft | Method and system for patient-specific computational modeling and simulation for coupled hemodynamic analysis of cerebral vessels |
US9724164B2 (en) | 2011-08-26 | 2017-08-08 | EBM Corporation | Blood-vessel bloodstream simulation system, method therefor, and computer software program |
US10162932B2 (en) | 2011-11-10 | 2018-12-25 | Siemens Healthcare Gmbh | Method and system for multi-scale anatomical and functional modeling of coronary circulation |
US10373700B2 (en) | 2012-03-13 | 2019-08-06 | Siemens Healthcare Gmbh | Non-invasive functional assessment of coronary artery stenosis including simulation of hyperemia by changing resting microvascular resistance |
US9135699B2 (en) * | 2012-03-15 | 2015-09-15 | Siemens Aktiengesellschaft | Method and system for hemodynamic assessment of aortic coarctation from medical image data |
JP2014028124A (ja) | 2012-06-27 | 2014-02-13 | Toshiba Corp | 医用画像処理装置及びx線診断装置 |
US10398386B2 (en) | 2012-09-12 | 2019-09-03 | Heartflow, Inc. | Systems and methods for estimating blood flow characteristics from vessel geometry and physiology |
EP2916735B1 (en) | 2012-11-06 | 2019-09-25 | Koninklijke Philips N.V. | Fractional flow reserve (ffr) index |
JP2014100249A (ja) * | 2012-11-19 | 2014-06-05 | Toshiba Corp | 血管解析装置、医用画像診断装置、血管解析方法、及び血管解析プログラム |
JP6334902B2 (ja) | 2012-11-30 | 2018-05-30 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置 |
JP6091870B2 (ja) | 2012-12-07 | 2017-03-08 | 東芝メディカルシステムズ株式会社 | 血管解析装置、医用画像診断装置、血管解析方法、及び血管解析プログラム |
JP5953438B2 (ja) * | 2012-12-11 | 2016-07-20 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 冠状動脈を通る血流量を決定する方法 |
CN104871161A (zh) | 2012-12-18 | 2015-08-26 | 皇家飞利浦有限公司 | 用于在从估计的心脏射血输出导出的患者特异性边界条件下模拟血流量的方法和装置 |
US9042613B2 (en) | 2013-03-01 | 2015-05-26 | Heartflow, Inc. | Method and system for determining treatments by modifying patient-specific geometrical models |
AU2014296238B2 (en) | 2013-07-30 | 2016-10-06 | Heartflow, Inc. | Method and system for modeling blood flow with boundary conditions for optimized diagnostic performance |
US9700219B2 (en) * | 2013-10-17 | 2017-07-11 | Siemens Healthcare Gmbh | Method and system for machine learning based assessment of fractional flow reserve |
US10595806B2 (en) | 2013-10-22 | 2020-03-24 | Koninklijke Philips N.V. | Fractional flow reserve (FFR) index with adaptive boundary condition parameters |
JP6362853B2 (ja) * | 2013-11-20 | 2018-07-25 | キヤノンメディカルシステムズ株式会社 | 血管解析装置、および血管解析装置の作動方法 |
US9390232B2 (en) * | 2014-03-24 | 2016-07-12 | Heartflow, Inc. | Systems and methods for modeling changes in patient-specific blood vessel geometry and boundary conditions |
EP3140757B1 (en) * | 2014-05-05 | 2020-06-24 | Siemens Healthcare GmbH | Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in models based on diseased patients and hypothetical normal anatomical models |
US10206587B2 (en) * | 2014-05-16 | 2019-02-19 | Toshiba Medical Systems Corporation | Image processing apparatus, image processing method, and storage medium |
US9747525B2 (en) * | 2014-06-16 | 2017-08-29 | Siemens Healthcare Gmbh | Method and system for improved hemodynamic computation in coronary arteries |
WO2016001017A1 (en) | 2014-06-30 | 2016-01-07 | Koninklijke Philips N.V. | Apparatus for determining a fractional flow reserve value |
JP6778174B2 (ja) | 2014-07-18 | 2020-10-28 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 狭窄評価 |
KR20160026808A (ko) | 2014-08-29 | 2016-03-09 | 강원대학교산학협력단 | 환자별 혈관 정보 결정 방법 |
CN105513036B (zh) * | 2014-09-26 | 2019-05-31 | 上海联影医疗科技有限公司 | 三维ct图像的分割方法及装置 |
JP6505444B2 (ja) | 2015-01-16 | 2019-04-24 | キヤノンメディカルシステムズ株式会社 | 観察装置 |
US9928609B2 (en) | 2015-03-10 | 2018-03-27 | Toshiba Medical Systems Corporation | Medical image processing apparatus and medical image processing method |
CN106033603B (zh) * | 2015-03-10 | 2019-06-11 | 东芝医疗系统株式会社 | 医用图像处理装置和医用图像处理方法 |
US10716513B2 (en) | 2015-04-17 | 2020-07-21 | Heartflow, Inc. | Systems and methods for cardiovascular blood flow and musculoskeletal modeling for predicting device failure or clinical events |
DE102015207596A1 (de) * | 2015-04-24 | 2016-10-27 | Siemens Healthcare Gmbh | Verfahren sowie Rechen- und Druckeinheit zum Erstellen einer Gefäßstütze |
CN105654497B (zh) * | 2016-01-20 | 2019-03-01 | 华北电力大学(保定) | 一种血管内光声图像的时间反演重建方法 |
CN106127819B (zh) | 2016-06-30 | 2019-10-08 | 上海联影医疗科技有限公司 | 医学图像中提取血管中心线的方法及其装置 |
-
2017
- 2017-01-23 EP EP17731461.4A patent/EP3375364A4/en active Pending
- 2017-01-23 WO PCT/CN2017/072256 patent/WO2018133118A1/zh active Application Filing
- 2017-04-19 CN CN201710257724.8A patent/CN107123112B/zh active Active
- 2017-04-19 CN CN201710257262.XA patent/CN107115111B/zh active Active
- 2017-04-19 CN CN201710258086.1A patent/CN107123159A/zh active Pending
- 2017-04-26 US US15/498,428 patent/US10360682B2/en active Active
- 2017-08-31 US US15/693,375 patent/US10325369B2/en active Active
- 2017-11-10 US US15/809,975 patent/US9984465B1/en active Active
-
2019
- 2019-07-19 US US16/517,514 patent/US10776922B2/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101347343A (zh) * | 2007-07-17 | 2009-01-21 | 株式会社东芝 | 超声波图像取得装置 |
CN101527047A (zh) * | 2008-03-05 | 2009-09-09 | 深圳迈瑞生物医疗电子股份有限公司 | 使用超声图像检测组织边界的方法与装置 |
CN102274046A (zh) * | 2010-06-08 | 2011-12-14 | 株式会社东芝 | 超声波诊断装置、超声波图像处理装置以及医用图像诊断装置 |
CN105190630A (zh) * | 2013-01-15 | 2015-12-23 | 凯瑟沃克斯有限公司 | 计算血流储备分数 |
CN104720851A (zh) * | 2013-12-23 | 2015-06-24 | Ge医疗系统环球技术有限公司 | 计算机断层扫描成像的方法和装置 |
CN105096388A (zh) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | 基于计算流体力学的冠状动脉血流仿真系统和方法 |
CN104720894A (zh) * | 2015-02-11 | 2015-06-24 | 中山大学附属第一医院 | 一种血管手术方式的合理性分析方法 |
CN106327487A (zh) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
CN106539622A (zh) * | 2017-01-28 | 2017-03-29 | 北京欣方悦医疗科技有限公司 | 基于血流动力学分析的冠脉虚拟支架植入方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107123159A (zh) | 2017-09-01 |
EP3375364A1 (en) | 2018-09-19 |
EP3375364A4 (en) | 2019-01-23 |
WO2018133118A1 (zh) | 2018-07-26 |
US20180211386A1 (en) | 2018-07-26 |
US10776922B2 (en) | 2020-09-15 |
CN107123112A (zh) | 2017-09-01 |
US9984465B1 (en) | 2018-05-29 |
US20180211388A1 (en) | 2018-07-26 |
CN107115111B (zh) | 2020-11-27 |
US10360682B2 (en) | 2019-07-23 |
US20190340764A1 (en) | 2019-11-07 |
CN107115111A (zh) | 2017-09-01 |
US10325369B2 (en) | 2019-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107123112B (zh) | 血流状态分析系统及方法 | |
CN106887000B (zh) | 医学图像的网格化处理方法及其系统 | |
US20210358634A1 (en) | Systems and methods for image processing to determine blood flow | |
US11357469B2 (en) | Systems and methods for assessing the severity of plaque and/or stenotic lesions using contrast distribution predictions and measurements | |
CN102525443B (zh) | 用于基于医学图像数据对心血管循环建模的方法和设备 | |
CN108550388B (zh) | 基于血管造影的计算血管血流参数的装置和系统 | |
WO2018133098A1 (zh) | 血管壁应力应变状态获取方法及系统 | |
US20230310085A1 (en) | Systems and methods for estimation of blood flow using response surface and reduced order modeling | |
Gharleghi et al. | Deep learning for time averaged wall shear stress prediction in left main coronary bifurcations | |
CN114708974A (zh) | 对新冠肺炎患者的住院时长进行预测的方法和相关产品 | |
Roy | Evaluation of accuracy of patient-specific non-linear biomechanical models for predicting intra-operative brain shift |
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 | ||
CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
CP01 | Change in the name or title of a patent holder |