CN107123159A - 血流状态分析***及方法 - Google Patents
血流状态分析***及方法 Download PDFInfo
- Publication number
- CN107123159A CN107123159A CN201710258086.1A CN201710258086A CN107123159A CN 107123159 A CN107123159 A CN 107123159A CN 201710258086 A CN201710258086 A CN 201710258086A CN 107123159 A CN107123159 A CN 107123159A
- Authority
- CN
- China
- Prior art keywords
- blood
- vessel
- blood vessel
- state
- phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000017531 blood circulation Effects 0.000 title claims abstract description 83
- 238000004458 analytical method Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims description 106
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 185
- 230000002792 vascular Effects 0.000 claims abstract description 74
- 210000004351 coronary vessel Anatomy 0.000 claims description 28
- 210000001367 artery Anatomy 0.000 claims description 25
- 230000036772 blood pressure Effects 0.000 claims description 21
- 210000003141 lower extremity Anatomy 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 8
- 230000003187 abdominal effect Effects 0.000 claims description 7
- 208000000058 Anaplasia Diseases 0.000 claims description 4
- 238000012545 processing Methods 0.000 description 83
- 210000004369 blood Anatomy 0.000 description 35
- 239000008280 blood Substances 0.000 description 35
- 230000002452 interceptive effect Effects 0.000 description 30
- 230000000004 hemodynamic effect Effects 0.000 description 29
- 239000012530 fluid Substances 0.000 description 28
- 230000002159 abnormal effect Effects 0.000 description 21
- 238000004364 calculation method Methods 0.000 description 21
- 210000001519 tissue Anatomy 0.000 description 20
- 238000010586 diagram Methods 0.000 description 17
- 210000000056 organ Anatomy 0.000 description 16
- 230000006870 function Effects 0.000 description 12
- 210000002216 heart Anatomy 0.000 description 12
- 230000008569 process Effects 0.000 description 11
- 230000000747 cardiac effect Effects 0.000 description 8
- 238000004891 communication Methods 0.000 description 8
- 230000002107 myocardial effect Effects 0.000 description 8
- 230000009467 reduction Effects 0.000 description 7
- 230000011218 segmentation 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
- 238000005192 partition Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000003709 image segmentation Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000007935 neutral effect Effects 0.000 description 4
- 230000033764 rhythmic process Effects 0.000 description 4
- 101000746134 Homo sapiens DNA endonuclease RBBP8 Proteins 0.000 description 3
- 101000969031 Homo sapiens Nuclear protein 1 Proteins 0.000 description 3
- 102100021133 Nuclear protein 1 Human genes 0.000 description 3
- 210000001015 abdomen Anatomy 0.000 description 3
- 230000000712 assembly Effects 0.000 description 3
- 238000000429 assembly Methods 0.000 description 3
- 238000013481 data capture Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 3
- 210000003462 vein Anatomy 0.000 description 3
- 206010002329 Aneurysm Diseases 0.000 description 2
- 208000007536 Thrombosis Diseases 0.000 description 2
- 230000001594 aberrant effect Effects 0.000 description 2
- 230000033115 angiogenesis Effects 0.000 description 2
- 230000037396 body weight Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000002526 effect on cardiovascular system Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 201000011066 hemangioma Diseases 0.000 description 2
- 238000003706 image smoothing Methods 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 230000009545 invasion Effects 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 208000002251 Dissecting Aneurysm Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000000746 body region Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 210000005242 cardiac chamber Anatomy 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 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
- 238000002845 discoloration Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 230000001435 haemodynamic effect Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 208000030613 peripheral artery disease Diseases 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 208000037804 stenosis Diseases 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
- 210000002385 vertebral artery Anatomy 0.000 description 1
Classifications
-
- 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
-
- 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
- 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
-
- 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
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Cardiology (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Hematology (AREA)
- Data Mining & Analysis (AREA)
- Vascular Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Databases & Information Systems (AREA)
- Mathematical Physics (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Computer Graphics (AREA)
Abstract
本发明提供一种血流状态分析方法和***。在确定血流状态的过程中,获取多个时相的图像,建立与多个时相分别对应的多个血管模型。关联所述多个血管模型,根据关联结果分别设定所述多个血管模型的边界条件,确定所述多个血管模型的血管状态,根据所述第一血管模型中所述血管在所述第一时相的状态和所述第二血管模型中所述血管在所述第二时相的状态,生成所述血管的状态与时间的关系。
Description
【技术领域】
本发明涉及一种血流状态分析***和方法,尤其是一种利用计算流体力 学获取多时相血流动力学参数的***及方法。
【背景技术】
CTA、MRA成像技术广泛应用在外周血管疾病诊断中,尤其在血管狭窄(椎 动脉狭窄)、动脉瘤、夹层动脉瘤、肿瘤供血动脉等血管疾病等方面。血管 分析软件应用提供了血管分析的工具,用于细小四级血管的精确提取、完整 去骨、快速自动测量等。医学图像分析软件***在针对血管的分析中一般采 用图像分割技术,以及图像显示技术对相应对象的血管进行形态学上的三维 模拟重建。医生通过针对血管形态学上的指标(血管狭窄程度、血管瘤膨胀 程度等)对病变程度进行分析处理,可利用信息单一。
计算流体力学(Computational Fluid Dynamics/CFD)是一个20世纪50 年代以来,随着计算机的发展而产生的一个介于数学、流体力学和计算机之 间的交叉学科。其主要研究内容是通过计算机和数值方法来求解流体力学的 控制方程,对流体力学问题进行模拟和分析。血管模型或血流模型是一个利 用计算流体力学的新兴应用。利用计算流体力学对单时相数据进行分析难以 全面的体现分析区域的实际情况和变化规律,也会由于选择时相不准确而使 结果出现偏差。
【发明内容】
本发明解决的问题是提供一种血管流体分析方法,血管流体分析方法,执 行于至少一个包含处理器和存储器的设备,所述方法包括:获取一个第一时 相的图像和一个第二时相的图像;在所述第一时相的图像中选择一个第一血 管区域,其中,所述第一血管区域包含一个血管;在所述第二时相的图像中 选择一个第二血管区域,其中,血流状态的分析方法,其特征在于:所述方 法包括:获取一个第一时相的图像和一个第二时相的图像;在所述第一时相 的图像中选择一个第一血管区域,其中,所述第一血管区域包含一个血管; 在所述第二时相的图像中选择一个第二血管区域,其中,所述第二血管区域 包含至少一部分所述血管;建立一个第一血管模型,其中,所述第一血管模 型与所述第一血管区域相对应;建立一个第二血管模型,其中,所述第二血 管模型与所述第二血管区域相对应;设定所述第一血管模型的边界条件和所 述第二血管模型的边界条件;根据所述第一血管模型的边界条件,确定所述 第一血管模型中所述血管在所述第一时相的状态;基于所述血管在所述第一 时相的状态,关联所述第一血管模型和所述第二血管模型;根据所述关联结 果以及所述第二血管模型的边界条件,确定所述第二血管模型中所述血管在 所述第二时相的状态;以及根据所述第一血管模型中所述血管在所述第一时 相的状态和所述第二血管模型中所述血管在所述第二时相的状态,生成所述 血管的状态与时间的关系。
可选的,所述分析方法进一步包括根据所述血管的状态与时间的所述关 系,生成所述血管在第三时相的状态。
可选的,一种分析方法,被至少一个包含处理器和存储器的设备执行,所 述分析方法包括:获取多个时相的血管图像,包括一个第一时相的血管图像 和一个第二时相的血管图像,所述多个时相的血管图像分别对应于同一个血 管或其一部分;建立多个血管模型,所述多个血管模型分别对应于所述多个 时相的血管图像;根据所述多个血管模型,获得所述血管或其一部分的多个 状态,包括一个第一血管状态和一个第二血管状态,其中,所述第一血管状 态对应于所述第一时相的血管图像,所述第二血管状态对应于所述第二时相 的血管图像;根据所述血管或其一部分的所述多个状态,获得所述血管或其 一部分的状态随时间变化的关系;以及根据所述关系,获得所述血管或其一 部分的第三血管状态。
可选的,所述血管包括冠状动脉、腹部动脉、大脑动脉或下肢动脉中至少 一种。
可选的,所述血管或其一部分的状态包括血液流速、血液压力、血管壁应 力、血管壁切应力、血流储备系数,以及冠状动脉血流储备中至少一种的数 值。
可选的,所述第三血管状态包括一个平均血流储备系数的值。
可选的,一个血流状态的分析***,包括至少一个处理器和存储设备,其 特征在于:所述***包括:一个接收模块,被配置为:获取多个时相的血管 图像,包括一个第一时相的血管图像和一个第二时相的血管图像,所述多个 时相的血管图像分别对应于同一个血管或其一部分;以及一个多时相特征生 成模块,被配置为:建立多个血管模型,所述多个血管模型分别对应于所述 多个时相的血管图像;根据所述多个血管模型,获得所述血管或其一部分的 多个状态,包括一个第一血管状态和一个第二血管状态,其中,所述第一血 管状态对应于所述第一时相的血管图像,所述第二血管状态对应于所述第二 时相的血管图像;根据所述血管或其一部分的所述多个状态,获得所述血管 或其一部分的状态随时间变化的关系;以及根据所述关系,获得所述血管或 其一部分的第三血管状态。
可选的,一种血流状态的分析***,包括:存储器,用来存储指令,所述 指令被所述至少一个处理器执行时,导致所述***实现的操作包括:获取多 个时相的血管图像,包括一个第一时相的血管图像和一个第二时相的血管图 像,所述多个时相的血管图像分别对应于同一个血管或其一部分;建立多个 血管模型,所述多个血管模型分别对应于所述多个时相的血管图像;根据所 述多个血管模型,获得所述血管或其一部分的多个状态,包括一个第一血管 状态和一个第二血管状态,其中,第一血管状态对应于所述第一时相的血管图像,第二血管状态对应于所述第二时相的血管图像;根据所述血管或其一 部分的所述多个状态,获得所述血管或其一部分的状态随时间变化的关系; 以及根据所述关系,获得所述血管或其一部分的第三血管状态。
【附图说明】
图1A和图1B是根据本申请的一些实施例所示的包括血管状态分析***;
图2是根据本申请的一些实施例所示的一个计算设备的结构,该计算设 备可以实施本申请中披露的特定***;
图3是根据本申请的一些实施例所示的一个移动设备的结构示意图,该 移动设备可以实施本申请中披露的特定***;
图4A是根据本申请的一些实施例所示的处理设备的示例性模块示意图;
图4B是根据本申请的一些实施例所示的多时相特征处理的示例性流程 图;
图5是根据本申请的一些实施例所示的多时相特征生成模块的示例性模 块示意图;
图6是根据本申请的一些实施例所示的获取多时相特征的示例性流程图;
图7是根据本申请的一些实施例所示的获取多时相特征的示意图;
图8是根据本申请的一些实施例所示的设定边界条件的示例性流程图;
图9是根据本申请的一些实施例所示的血流模型的示意图;
图10是根据本申请的一些实施例所示的获得任意位点参数的示例性流程 图;以及
图11是根据本申请的一些实施例所示的获得位点血流动力学参数的示意 图。
【具体实施方式】
在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发 明能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不 违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施 的限制。
本申请涉及一种获取血流状态的方法和***。在确定血流状态的过程中, 获取多个时相的图像,建立与多个时相分别对应的多个血管模型。关联所述 多个血管模型,根据关联结果分别设定所述多个血管模型的边界条件,确定 所述多个血管模型的血管状态。所述关联所述多个血管模型包括配准多个血 管模型。
图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及其描述。在一些实施例中,所述分析和计算的结果可以 包括血管、组织、或器官的物理状态和相关系数或参数。例如,对冠状动脉 模型进行分析和计算的结果可以包括冠状动脉的血流动力学参数,如血流速 率、血液压力、血管壁应力、血管壁切应力、血流储备系数或血流储备分数 (Fractional Flow Reserve,FFR)、冠状动脉血流储备(CoronaryFlow Reserve,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-accessmemory,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可以进一步生成任意时相的分析结果的估计值。具体的后处理的流 程和结果可以参考图10及其描述。在一些实施例中,所述多时相特征处理模 块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可 以在模型的边界区域(入口、出口等)生成二维的网格,而在模型的其他区 域生成三维的网格。所述三维网格可以是基于所述二维网格而建立的。
匹配单元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进行网格化处理后的结果。
在610中,可以设定多时相的参数和边界条件。在一些实施例中,所述 参数的设置可以包括选定的模型中参数的设置,如流速u、密度ρ、血压P、 截面积S等。所述边界条件的设置可以包括对边界区域进行初始条件或限制 条件的设置。所述边界区域指的是感兴趣区域的边缘区域。例如,边界区域 可以是血管的出口、入口、血管壁等。设置的边界条件可以包括边界区域的 血压、血流量、血流速率、血管流阻、压强、应力等。在一些实施例中,可以根据感兴趣的区域将一个低阶耦合模型作为边界条件。所述低阶耦合模型 可以利用与感兴趣区域相耦合的区域或组织的经验性模型作为边界条件。所 述低阶耦合模型可以是二阶模型,一阶模型,零阶模型(集中参数模型)或 这些低阶模型的组合形式。。如图7所示,712可以是模型的选择和参数的设 置的一个实施例。714可以是边界条件的设定的一个实施例。
在612中,可以选择一个当前时相(在第一次选择的时候也可称作初始 时相)。在一些实施例中,初始时相可以依据一些特定的规则进行选择出来。 例如,对于冠状动脉模型,可以选择一个模型变化较慢或较平缓的时相作为 初始时相,例如,离心脏收缩前期或者舒张末期最近的时相。初始时相可以 由机器(例如,多时相特征生成模块430)自行判断选择,或由用户选择。如 果机器和用户均没有或无法选择一个初始时相,可以任意选择一个时相或者 选择多时相特征生成模块430第一个接收到的时相作为初始时相。
在614中,可以对当前时相(在第一次执行时为初始时相)进行分析, 例如,计算流体力学(CFD)分析。根据设定的模型、边界条件以及参数可以 对三维血管模型的血流动力学进行求解。所述求解的控制方程可以基于Euler 方程组、纳维尔-斯托克斯方程组(Navier-Stokes Equations)或格子玻尔 兹曼方程(Lattice Boltzmann Method)。求解计算的离散化方法可以包括有 限差分法、有限体积法、有限元法、边界元方法、谱方法、格子玻尔兹曼法、 无网格法等一种或多种组合。求解计算流场的流体可以是无粘的或有粘性的, 可以是可压缩流体也可以是不可压缩流体,可以是层流也可以是湍流,可以 是定常流动也可以是非定常流动。可以根据所模拟流体的物理特性,相应地 选择相应的控制方程以及模拟方法。例如,针对无粘流体的流场计算可以选 用Euler方程组或者格子玻尔兹曼方程,针对有粘性流体的流场计算可以选 用N-S方程组或者玻尔兹曼方程。例如,对于冠状动脉进行的计算流体力学 (CFD)计算可以使用纳维尔-斯托克斯方程(Navier-StokesEquations):
和
其中,ρ代表血液密度,u代表血流的速度,t代表时间,以及σ代表血流应力 项(由血压p和血液粘性力项决定)。在一些实施例中,在初始时相的计算中, 可以将所述模型中的流速的初始值设置为0。在后续时相的计算中,所述模型 中流速的初始值不再设置为0,而是匹配前后时相的网格,并根据上一时相的 计算结果对当前时相的对应网格进行赋值,作为当前时相流速的初始值。
分析结果可以包括当前时相该模型任意一个区域或者位点的物理状态和 相关系数或参数。例如,对冠状动脉模型进行分析的结果可以包括冠状动脉 任意一个区域或者位点的血流动力学参数,如血流速率、血液压力、血管壁 应力、血管壁切应力、血流储备系数(FFR)、冠状动脉血流储备(CFR)等中 的一种或多种的组合。如图7所示,718显示一种血流动力学的分析计算。
在616中,可以判断是否已遍历所有时相。在判断出已遍历所有时相后, 可以执行618。在判断出未遍历所有时相后,可以执行620。
在618中,可以将分析结果进行输出。例如,将分析结果发送至血流状态 分析***100其他模块或设备中。在一些实施例中,可以对分析结果进行后 处理。所述后处理可以包括生成分析结果与时相的关系曲线或对照表。所述 后处理也可以包括根据所述关系曲线或对照表进一步输出任意时相的分析结 果的预估值。如图7所示,716显示一种后处理的结果。在一些实施例中,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是描述获得一个位点对应的血流动力学 参数的示例性流程图。在一些实施例中,获得一个位点对应的血流动力学参 数的流程可以通过多时相特征处理模块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值)。
以上的描述仅仅是本发明的具体实施例,不应被视为是唯一的实施例。在 一些实施例中,在进行拟合位点的血流动力学参数曲线之前,多时相特征处 理模块440可以获得为点额外的血流动力学参数。所述额外的血流动力学参 数可以是通过插值方法处理后获得,也可以是用户通过,交互设备140,输入 后获得。
根据本申请的一些实施例,图11是描述获得位点对应的血流动力学参数 的示意图。图片1402显示包含心脏区域和腹部区域的多时相图像数据。图片 1406显示一个时相的图像中对应的冠脉以及建立的冠脉模型。图像1408显示 所述冠脉在不同时相的FFR分布。所述不同时相的血液流动状态可以根据所 述获取多时相特征的示例性流程获得。图片1404显示对象的特异性临床数据, 包括对象的主动脉压随时间的变化曲线,相位冠脉血流量随时间的变化曲线。 图片1410显示对象的FFR随时间变化的曲线。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任 何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的 方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱 离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何 简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
Claims (10)
1.一种血流状态分析方法,其特征在于:所述方法包括:
获取一个第一时相的图像和一个第二时相的图像;
在所述第一时相的图像中选择一个第一血管区域,其中,所述第一血管区域包含一个血管;
在所述第二时相的图像中选择一个第二血管区域,其中,所述第二血管区域包含至少一部分所述血管;
建立一个第一血管模型,其中,所述第一血管模型与所述第一血管区域相对应;
建立一个第二血管模型,其中,所述第二血管模型与所述第二血管区域相对应;
设定所述第一血管模型的边界条件和所述第二血管模型的边界条件;
根据所述第一血管模型的边界条件,确定所述第一血管模型中所述血管在所述第一时相的状态;
基于所述血管在所述第一时相的状态,关联所述第一血管模型和所述第二血管模型;
根据所述关联结果以及所述第二血管模型的边界条件,确定所述第二血管模型中所述血管在所述第二时相的状态;以及
根据所述第一血管模型中所述血管在所述第一时相的状态和所述第二血管模型中所述血管在所述第二时相的状态,生成所述血管的状态与时间的关系。
2.如权利要求1所述的方法,其特征在于:所述方法进一步包括根据所述血管的状态与时间的所述关系,生成所述血管在第三时相的状态。
3.一种方法,被至少一个包含处理器和存储器的设备执行,所述方法包括:
获取多个时相的血管图像,包括一个第一时相的血管图像和一个第二时相的血管图像,所述多个时相的血管图像分别对应于同一个血管或其一部分;
建立多个血管模型,所述多个血管模型分别对应于所述多个时相的血管图像;
根据所述多个血管模型,获得所述血管或其一部分的多个状态,包括一个第一血管状态和一个第二血管状态,其中,所述第一血管状态对应于所述第一时相的血管图像,所述第二血管状态对应于所述第二时相的血管图像;
根据所述血管或其一部分的所述多个状态,获得所述血管或其一部分的状态随时间变化的关系;以及
根据所述关系,获得所述血管或其一部分的第三血管状态。
4.如权利要求3所述的方法,其特征在于:所述血管包括冠状动脉、腹部动脉、大脑动脉或下肢动脉中至少一种。
5.如权利要求3所述的方法,其特征在于:所述血管或其一部分的状态包括血液流速、血液压力、血管壁应力、血管壁切应力、血流储备系数,以及冠状动脉血流储备中至少一种的数值。
6.如权利要求3所述的方法,其特征在于:所述第三血管状态包括一个平均血流储备系数的值。
7.一个血流状态分析***,包括至少一个处理器和存储设备,其特征在于:所述***包括:
一个接收模块,被配置为:
获取多个时相的血管图像,包括一个第一时相的血管图像和一个第二时相的血管图像,所述多个时相的血管图像分别对应于同一个血管或其一部分;以及
一个多时相特征生成模块,被配置为:
建立多个血管模型,所述多个血管模型分别对应于所述多个时相的血管图像;
根据所述多个血管模型,获得所述血管或其一部分的多个状态,包括一个第一血管状态和一个第二血管状态,其中,所述第一血管状态对应于所述第一时相的血管图像,所述第二血管状态对应于所述第二时相的血管图像;
根据所述血管或其一部分的所述多个状态,获得所述血管或其一部分的状态随时间变化的关系;以及
根据所述关系,获得所述血管或其一部分的第三血管状态。
8.如权利要求7所述的***,其特征在于:所述血管或其一部分的状态包括血液流速、血液压力、血管壁应力、血管壁切应力、血流储备系数,以及冠状动脉血流储备中至少一种的数值。
9.如权利要求7所述的***,其特征在于:所述第三血管状态包括一个平均血流储备系数的值。
10.一种血流状态分析***,包括:
存储器,用来存储指令,所述指令被所述至少一个处理器执行时,导致所述***实现的操作包括:
获取多个时相的血管图像,包括一个第一时相的血管图像和一个第二时相的血管图像,所述多个时相的血管图像分别对应于同一个血管或其一部分;
建立多个血管模型,所述多个血管模型分别对应于所述多个时相的血管图像;
根据所述多个血管模型,获得所述血管或其一部分的多个状态,包括一个第一血管状态和一个第二血管状态,其中,第一血管状态对应于所述第一时相的血管图像,第二血管状态对应于所述第二时相的血管图像;
根据所述血管或其一部分的所述多个状态,获得所述血管或其一部分的状态随时间变化的关系;以及
根据所述关系,获得所述血管或其一部分的第三血管状态。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2017/072256 WO2018133118A1 (zh) | 2017-01-23 | 2017-01-23 | 血流状态分析***及方法 |
CNPCT/CN2017/072256 | 2017-01-23 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107123159A true CN107123159A (zh) | 2017-09-01 |
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 Before (2)
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 | 血流状态分析***及方法 |
Country Status (4)
Country | Link |
---|---|
US (4) | US10360682B2 (zh) |
EP (2) | EP4393386A2 (zh) |
CN (3) | CN107123112B (zh) |
WO (1) | WO2018133118A1 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109620199A (zh) * | 2018-11-30 | 2019-04-16 | 博动医学影像科技(上海)有限公司 | 建立血管截面函数、血管压力差和血管应力的方法及装置 |
WO2020093971A1 (zh) * | 2018-11-06 | 2020-05-14 | 北京三普威盛科技有限公司 | 用于冠脉狭窄评估的方法,装置,存储介质及电子设备 |
CN113925598A (zh) * | 2020-07-13 | 2022-01-14 | 通用电气精准医疗有限责任公司 | 利用超声成像***对组织进行热监测的方法和*** |
CN113963001A (zh) * | 2021-10-25 | 2022-01-21 | 中国人民解放军总医院第一医学中心 | 一种基于多时相cta的侧支循环血管定量方法 |
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 |
US11087459B2 (en) | 2015-08-14 | 2021-08-10 | Elucid Bioimaging Inc. | Quantitative imaging for fractional flow reserve (FFR) |
US11113812B2 (en) | 2015-08-14 | 2021-09-07 | Elucid Bioimaging Inc. | Quantitative imaging for detecting vulnerable plaque |
US11071501B2 (en) | 2015-08-14 | 2021-07-27 | Elucid Bioiwaging Inc. | Quantitative imaging for determining time to adverse event (TTE) |
US12026868B2 (en) | 2015-08-14 | 2024-07-02 | Elucid Bioimaging Inc. | Quantitative imaging for detecting histopathologically defined plaque erosion non-invasively |
US11094058B2 (en) | 2015-08-14 | 2021-08-17 | Elucid Bioimaging Inc. | Systems and method for computer-aided phenotyping (CAP) using radiologic images |
US12008751B2 (en) | 2015-08-14 | 2024-06-11 | Elucid Bioimaging Inc. | Quantitative imaging for detecting histopathologically defined plaque fissure non-invasively |
US10176408B2 (en) | 2015-08-14 | 2019-01-08 | Elucid Bioimaging Inc. | Systems and methods for analyzing pathologies utilizing quantitative imaging |
US11676359B2 (en) | 2015-08-14 | 2023-06-13 | Elucid Bioimaging Inc. | Non-invasive quantitative imaging biomarkers of atherosclerotic plaque biology |
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 | 深圳科亚医疗科技有限公司 | 预测血管树的血流储备分数的设备、***以及非暂时性可读存储介质 |
WO2020044523A1 (ja) * | 2018-08-30 | 2020-03-05 | オリンパス株式会社 | 記録装置、画像観察装置、観察システム、観察システムの制御方法、及び観察システムの作動プログラム |
CN110428420B (zh) * | 2018-09-05 | 2022-05-17 | 深圳科亚医疗科技有限公司 | 基于患者的冠状动脉ct血管造影图像来确定冠状动脉的流动信息的方法、装置和介质 |
CN110367965B (zh) * | 2018-09-19 | 2022-03-08 | 苏州润迈德医疗科技有限公司 | 便捷测量冠状动脉血管评定参数的方法、装置及*** |
CN112384138B (zh) * | 2018-10-25 | 2023-07-25 | 苏州润迈德医疗科技有限公司 | 获取心表大动脉的血流量的方法、装置、***及存储介质 |
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 | 合肥工业大学 | 用于研究颅内动脉瘤生长过程以及探究生长机制的方法 |
WO2021026125A1 (en) | 2019-08-05 | 2021-02-11 | Elucid Bioimaging Inc. | Combined assessment of morphological and perivascular disease markers |
US12039685B2 (en) | 2019-09-23 | 2024-07-16 | Cathworks Ltd. | Methods, apparatus, and system for synchronization between a three-dimensional vascular model and an imaging device |
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 | 杭州脉流科技有限公司 | 基于冠脉造影图像获取血流储备分数的方法、装置、设备、***及可读存储介质 |
CN111369519B (zh) * | 2020-03-02 | 2021-09-17 | 博动医学影像科技(上海)有限公司 | 冠状动脉的血流速度的计算方法、装置及电子设备 |
CN111815598B (zh) * | 2020-06-30 | 2024-04-26 | 上海联影医疗科技股份有限公司 | 一种血管参数的计算方法、装置、设备及存储介质 |
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 | 上海联影医疗科技股份有限公司 | 一种血管参数的计算方法、装置、设备及存储介质 |
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 | 北京圣点云信息技术有限公司 | 一种基于汇流累计量信息的体内静脉识别方法及装置 |
IL311466A (en) * | 2021-09-14 | 2024-05-01 | United Therapeutics Corp | Tissue evaluation systems and methods |
CN113706539A (zh) * | 2021-10-29 | 2021-11-26 | 南京裕隆生物医学发展有限公司 | 一种用于识别肿瘤的人工智能辅助*** |
US20230252632A1 (en) * | 2022-02-10 | 2023-08-10 | 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 (12)
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 |
EP2845537A2 (en) * | 2010-08-12 | 2015-03-11 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
JP2015097759A (ja) * | 2013-11-20 | 2015-05-28 | 株式会社東芝 | 血管解析装置、および血管解析方法 |
CN104736061A (zh) * | 2012-11-19 | 2015-06-24 | 株式会社东芝 | 血管解析装置、医用图像诊断装置及血管解析方法 |
CN104871211A (zh) * | 2012-12-11 | 2015-08-26 | 皇家飞利浦有限公司 | 确定通过冠状动脉的血液流量的方法 |
CN105096388A (zh) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | 基于计算流体力学的冠状动脉血流仿真***和方法 |
US20150359601A1 (en) * | 2014-06-16 | 2015-12-17 | Siemens Aktiengesellschaft | Method and System for Improved Hemodynamic Computation in Coronary Arteries |
CN105513036A (zh) * | 2014-09-26 | 2016-04-20 | 上海联影医疗科技有限公司 | 三维ct图像的分割方法及装置 |
CN105654497A (zh) * | 2016-01-20 | 2016-06-08 | 华北电力大学(保定) | 一种血管内光声图像的时间反演重建方法 |
CN106033603A (zh) * | 2015-03-10 | 2016-10-19 | 东芝医疗***株式会社 | 医用图像处理装置和医用图像处理方法 |
CN106063734A (zh) * | 2015-04-24 | 2016-11-02 | 西门子股份公司 | 用于创建支架移植物的方法以及计算和打印单元 |
CN106327487A (zh) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
Family Cites Families (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4960162B2 (ja) | 2007-07-17 | 2012-06-27 | 株式会社東芝 | 超音波診断装置 |
CN101527047B (zh) | 2008-03-05 | 2013-02-13 | 深圳迈瑞生物医疗电子股份有限公司 | 使用超声图像检测组织边界的方法与装置 |
JP5632203B2 (ja) | 2010-06-08 | 2014-11-26 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
US8682626B2 (en) | 2010-07-21 | 2014-03-25 | Siemens Aktiengesellschaft | Method and system for comprehensive patient-specific modeling of the heart |
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 |
CA2850476A1 (en) | 2011-08-26 | 2013-03-07 | EBM Corporation | Bloodstream simulation system for simulating blood vessel treatment effect, 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 |
CN103635138B (zh) * | 2012-06-27 | 2016-01-20 | 株式会社东芝 | X射线诊断装置 |
US10398386B2 (en) | 2012-09-12 | 2019-09-03 | Heartflow, Inc. | Systems and methods for estimating blood flow characteristics from vessel geometry and physiology |
US9858387B2 (en) | 2013-01-15 | 2018-01-02 | CathWorks, LTD. | Vascular flow assessment |
US9757073B2 (en) | 2012-11-06 | 2017-09-12 | Koninklijke Philips N.V. | Fractional flow reserve (FFR) index |
WO2014084382A1 (ja) | 2012-11-30 | 2014-06-05 | 株式会社 東芝 | 医用画像処理装置 |
JP6091870B2 (ja) | 2012-12-07 | 2017-03-08 | 東芝メディカルシステムズ株式会社 | 血管解析装置、医用画像診断装置、血管解析方法、及び血管解析プログラム |
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 |
KR102380008B1 (ko) | 2013-07-30 | 2022-03-29 | 하트플로우, 인크. | 최적화된 진단적 성과를 위한 경계 조건에서 혈류를 모형화하기 위한 방법과 시스템 |
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 |
CN104720851B (zh) * | 2013-12-23 | 2019-07-16 | Ge医疗***环球技术有限公司 | 计算机断层扫描成像的方法和装置 |
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 |
US10803995B2 (en) * | 2014-05-05 | 2020-10-13 | Siemens Healthcare Gmbh | Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in diseased 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 |
EP3160335B1 (en) | 2014-06-30 | 2023-09-20 | Koninklijke Philips N.V. | Apparatus for determining a fractional flow reserve value |
CN106572824B (zh) | 2014-07-18 | 2021-09-07 | 皇家飞利浦有限公司 | 狭窄评估 |
JP6377856B2 (ja) | 2014-08-29 | 2018-08-22 | ケーエヌユー−インダストリー コーポレーション ファウンデーション | 患者別の心血管情報を決定する方法 |
JP6505444B2 (ja) | 2015-01-16 | 2019-04-24 | キヤノンメディカルシステムズ株式会社 | 観察装置 |
CN104720894A (zh) * | 2015-02-11 | 2015-06-24 | 中山大学附属第一医院 | 一种血管手术方式的合理性分析方法 |
US9928609B2 (en) | 2015-03-10 | 2018-03-27 | Toshiba Medical Systems Corporation | Medical image processing apparatus and medical image processing method |
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 |
CN106127819B (zh) | 2016-06-30 | 2019-10-08 | 上海联影医疗科技有限公司 | 医学图像中提取血管中心线的方法及其装置 |
CN106539622B (zh) * | 2017-01-28 | 2019-04-05 | 北京欣方悦医疗科技有限公司 | 基于血流动力学分析的冠脉虚拟支架植入*** |
-
2017
- 2017-01-23 WO PCT/CN2017/072256 patent/WO2018133118A1/zh active Application Filing
- 2017-01-23 EP EP24176082.6A patent/EP4393386A2/en active Pending
- 2017-01-23 EP EP17731461.4A patent/EP3375364A4/en active Pending
- 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 (13)
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 |
CN102573638A (zh) * | 2009-10-13 | 2012-07-11 | 新加坡科技研究局 | 一种用于分割图像中的肝脏对象的方法和*** |
EP2845537A2 (en) * | 2010-08-12 | 2015-03-11 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
CN104736061A (zh) * | 2012-11-19 | 2015-06-24 | 株式会社东芝 | 血管解析装置、医用图像诊断装置及血管解析方法 |
CN104871211A (zh) * | 2012-12-11 | 2015-08-26 | 皇家飞利浦有限公司 | 确定通过冠状动脉的血液流量的方法 |
JP2015097759A (ja) * | 2013-11-20 | 2015-05-28 | 株式会社東芝 | 血管解析装置、および血管解析方法 |
CN105096388A (zh) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | 基于计算流体力学的冠状动脉血流仿真***和方法 |
US20150359601A1 (en) * | 2014-06-16 | 2015-12-17 | Siemens Aktiengesellschaft | Method and System for Improved Hemodynamic Computation in Coronary Arteries |
CN105513036A (zh) * | 2014-09-26 | 2016-04-20 | 上海联影医疗科技有限公司 | 三维ct图像的分割方法及装置 |
CN106033603A (zh) * | 2015-03-10 | 2016-10-19 | 东芝医疗***株式会社 | 医用图像处理装置和医用图像处理方法 |
CN106063734A (zh) * | 2015-04-24 | 2016-11-02 | 西门子股份公司 | 用于创建支架移植物的方法以及计算和打印单元 |
CN105654497A (zh) * | 2016-01-20 | 2016-06-08 | 华北电力大学(保定) | 一种血管内光声图像的时间反演重建方法 |
CN106327487A (zh) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020093971A1 (zh) * | 2018-11-06 | 2020-05-14 | 北京三普威盛科技有限公司 | 用于冠脉狭窄评估的方法,装置,存储介质及电子设备 |
CN109620199A (zh) * | 2018-11-30 | 2019-04-16 | 博动医学影像科技(上海)有限公司 | 建立血管截面函数、血管压力差和血管应力的方法及装置 |
US11445923B2 (en) | 2018-11-30 | 2022-09-20 | Shanghai Pulse Medical Technology, Inc. | Method and device for establishing blood vessel cross-section function, blood stress vessel pressure difference and blood vessel stress |
CN113925598A (zh) * | 2020-07-13 | 2022-01-14 | 通用电气精准医疗有限责任公司 | 利用超声成像***对组织进行热监测的方法和*** |
CN113925598B (zh) * | 2020-07-13 | 2024-06-11 | 通用电气精准医疗有限责任公司 | 利用超声成像***对组织进行热监测的方法和*** |
CN113963001A (zh) * | 2021-10-25 | 2022-01-21 | 中国人民解放军总医院第一医学中心 | 一种基于多时相cta的侧支循环血管定量方法 |
Also Published As
Publication number | Publication date |
---|---|
US20190340764A1 (en) | 2019-11-07 |
WO2018133118A1 (zh) | 2018-07-26 |
CN107115111B (zh) | 2020-11-27 |
US10325369B2 (en) | 2019-06-18 |
CN107123112A (zh) | 2017-09-01 |
EP4393386A2 (en) | 2024-07-03 |
EP3375364A4 (en) | 2019-01-23 |
US20180211386A1 (en) | 2018-07-26 |
US9984465B1 (en) | 2018-05-29 |
US10776922B2 (en) | 2020-09-15 |
EP3375364A1 (en) | 2018-09-19 |
CN107123112B (zh) | 2020-08-04 |
US20180211388A1 (en) | 2018-07-26 |
CN107115111A (zh) | 2017-09-01 |
US10360682B2 (en) | 2019-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107123159A (zh) | 血流状态分析***及方法 | |
AU2020244586B2 (en) | Method and system for sensitivity analysis in modeling blood flow characteristics | |
US12004841B2 (en) | Systems and methods for estimation of blood flow characteristics using reduced order model and/or machine learning | |
CN102525443B (zh) | 用于基于医学图像数据对心血管循环建模的方法和设备 | |
CN110853029B (zh) | 用于基于医学图像自动预测血流特征的方法、***和介质 | |
CN105249954B (zh) | 预测支架术后血液动力学指标的方法和*** | |
CN104854592B (zh) | 用于根据血管几何形状和生理学估计血流特性的***和方法 | |
CN108109698B (zh) | 计算血流储备分数的***和设置边界条件的方法 | |
CN105380598B (zh) | 用于针对动脉狭窄的自动治疗规划的方法和*** | |
US10622092B2 (en) | Method and system for facilitating physiological computations | |
CN106537392A (zh) | 用于冠状动脉中的血液动力学计算的方法和*** | |
CN107427268A (zh) | 用于基于纯几何机器学习的血流储备分数的方法和*** | |
CN108122616B (zh) | 个体特异性的心血管模型的生成方法及其应用 | |
CN108665449B (zh) | 预测血流矢量路径上的血流特征的深度学习模型及装置 | |
JP2019512131A (ja) | 画像ベースの患者固有の血流力学的モデルにおいて解像されていない血管を特定してモデル化するためのシステム及び方法 | |
CN106659399A (zh) | 使用患病和假想正常解剖学模型中的流计算的冠状动脉狭窄的非侵入功能评价的方法和*** | |
CN107665737A (zh) | 血管壁应力应变状态获取方法、计算机可读介质及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
CB02 | Change of applicant information | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170901 |
|
RJ01 | Rejection of invention patent application after publication |