CN107582041A - 管状构造解析装置及管状构造解析方法 - Google Patents
管状构造解析装置及管状构造解析方法 Download PDFInfo
- Publication number
- CN107582041A CN107582041A CN201710695287.8A CN201710695287A CN107582041A CN 107582041 A CN107582041 A CN 107582041A CN 201710695287 A CN201710695287 A CN 201710695287A CN 107582041 A CN107582041 A CN 107582041A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- index
- time series
- identification
- 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.)
- Granted
Links
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/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
-
- 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/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/504—Clinical applications involving diagnosis of blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- 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/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
- G06T2207/10136—3D ultrasound image
-
- 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
血管的构造流体解析的精度的提高。图像解析-跟踪处理部(53)将时间序列的医用图像进行图像处理,计算解析对象区域的时间序列的形态指标和时间序列的形状变形指标。力学模型构建部(55)基于时间序列的形态指标、时间序列的形状变形指标和时间序列的医用图像,暂定地构建关于解析对象区域的构造流体解析的力学模型。统计性鉴定部(61)鉴定关于潜在变量鉴定区域的潜在变量,以使基于暂定地构建出的力学模型的血管形态指标的预测值及血液流量指标的预测值的至少一方与血管形态指标的观测值及血液流量指标的观测值的至少一方匹配。计算部(57、59)对已将潜在变量分配给潜在变量鉴定区域的力学模型实施构造解析、流体解析或构造流体复合解析,计算时间序列的力学性指标的预测值和时间序列的血液流量指标的预测值的至少一方。
Description
本申请为2013年11月12日提交的申请号为201380054742.6、发明名称为《血管解析装置、医用图像诊断装置及血管解析方法》的分案申请。
技术领域
本实施方式涉及管状构造解析装置及管状构造解析方法。
背景技术
希望有用来非侵袭或低侵袭地预防、诊断冠状动脉的狭窄、脑动脉瘤或基于作为它们的前兆的颈动脉的斑块的狭窄的技术,所述冠状动脉的狭窄是作为三大疾病之一的心疾病的原因。
冠状动脉的狭窄是导致虚血性心疾病的重大的病变。作为冠状动脉的狭窄诊断,通过导管进行的冠状动脉造影检查(CAG:Coronary Angiography)是主流。作为冠状动脉的器质性病变的诊断指标,有心肌血流储备分数(FFR:Fractional Flow Reserve)。FFR被定义为狭窄存在下的最大冠脉血流相对于狭窄非存在下的最大冠脉血流的比率。FFR与狭窄远位部冠内压相对于狭窄近位部冠内压的比率大致一致。测量设在导管前端的压力传感器。即,FFR的测量需要导管手术。
冠状动脉的狭窄解析如果能够通过心脏CT进行,则与通过导管手术进行的FFR的测量相比,能够低侵袭且减轻患者的负担、以及节约医疗成本。然而,在心脏CT中,仅基于CT图像中包含的斑块区域或狭窄区域的大小的指标能够低侵袭地计测。如果能够基于CT图像通过构造流体解析来计测狭窄前后的压力差等,则能够期待狭窄(或斑块)带来的影响的定量化。
作为冠脉循环的动态评价,在临床上开发导入了包括超高速CT、心血管荧光电影照相、超声波法、SPECT(单光子发射计算机断层成像)及PET(正电子发射断层扫描)的核医学成像、MRI(核磁共振图像法)等,对诊断及治疗法的评价很有帮助。
但是,难以将冠脉微小血管用医用图像诊断装置正确地捕捉。此外,即使血管形状鲜明,在医用图像中包含噪声的情况或在生物体组织的边界的阈值设定中存在含糊性的情况也较多。这样,从医用图像诊断装置得到的血管形状具有不确定性。
在临床应用中利用医用图像诊断装置的情况下,从比冠脉微小血管靠上游的大动脉起始部起、仅将冠状动脉的较粗的区域作为对象进行解析的情况也较多。由于冠状动脉的血流也较大地受冠脉微小血管的紧张性(tonus,强直性)影响,所以适当地设定较粗的区域的冠状动脉的出口处的流量或压力、或者它们的变化率等的流体解析的边界条件成为课题。此外,冠状动脉的血流受到因心脏的搏动带来的机械性因素(因搏动带来的整体性的运动、由局部性的伸缩、扭转、剪切变形带来的强制位移或外力)影响。仅通过流体解析,不能考虑心脏的搏动等的机械性因素的影响,所以不能精确地计测血流的流量分布及内压分布。另一方面,还以在图像中捕捉到的心脏及血管系统为对象,实施考虑到机械性因素的影响的构造-流体复合解析。但是,在进行构造-流体复合解析的情况下,难以正确地设定血液(包括造影剂)的流体解析中的血管的入口及出口的边界条件、血管及斑块的材料模型的情况也较多。此外,在存在图像中未描绘出的微小血管的情况下,还有不能考虑微小血管给血流带来的影响的情况。因此,构造-流体复合解析的解析结果有可能不能将实际的血流及血管变形再现。此外,还有边界条件、负荷条件、材料模型不适当的情况、或在血管伴随着较大的运动的情况下收敛性及解析稳定性上有问题的情况。这样,以往的血管的构造流体解析有需要极大的解析资源和解析时间的情况、或解析结果的误差变大的情况,有在现实在临床的现场中利用上发生问题的情况。
现有技术文献
专利文献
专利文献1:日本特开2008-241432号公报
非专利文献
非专利文献1:门冈等,(ITUジャーナル,心臓シミュレータが拓くテーラーメード医療~世界最先端の心臓シミュレータとその適用例のご紹介~,Vol.41,No.6)
发明内容
有关本实施方式的管状构造解析装置,具备:设定部,对于与对象物的管状构造体有关的时间序列的图像设定解析对象区域,对于上述解析对象区域设定潜在变量鉴定(identification、日语:同定)区域;计算部,计算上述解析对象区域的时间序列的形态指标和时间序列的形状变形指标;构建部,基于上述时间序列的形态指标、上述时间序列的形状变形指标和上述时间序列的图像,构建与上述解析对象区域的构造流体解析有关的力学模型;以及鉴定部,对与上述潜在变量鉴定区域有关的潜在变量进行鉴定,以使基于构建出的上述力学模型的形态指标的预测值及流量指标的预测值中的至少一方与预先计测出的形态指标的观测值及流量指标的观测值中的至少一方匹配。
实施方式的目的是提供一种能够实现与血管(包括血液)的构造流体解析及图像解析-跟踪处理有关的精度的提高的管状构造解析装置及管状构造解析方法。
附图说明
图1是表示有关本实施方式的医用图像诊断装置(X射线计算机断层摄影装置)的概略性的块结构的图。
图2是表示与本实施方式的构造流体解析的对象区域有关的力学模型的一例的图。
图3是表示基于图1的系统控制部的控制而进行的构造流体解析处理的典型的流程的图。
图4是表示图1的图像处理装置的块结构的图。
图5是示意地表示与血管的芯线正交的截面的图。
图6是表示在图4的图像解析-跟踪处理的图像跟踪处理中使用的血管芯线的形态的时间性变化的图。
图7是表示图4的图像解析-跟踪处理的血管芯线的弯曲变形及旋转位移的计算例的图。
图8是用来说明图4的图像解析-跟踪处理的图像跟踪处理的图,是表示时刻t与时刻t+Δt之间的跟踪处理的一例的图。
图9是表示通过图4的力学模型构建部构建的形状模型的与芯线正交的截面的图。
图10是用来说明由图4的力学模型构建部进行的向形状模型的强制位移履历的分配的图。
图11是表示由图4的力学模型构建部进行的向形状模型的强制位移履历的其他分配方法的图,是表示血管形状变形指标为扭转的情况下的分配例的图。
图12是表示由图4的力学模型构建部进行的向形状模型的强制位移履历的其他分配方法的图,是表示血管形状变形指标为弯曲的情况下的分配例的图。
图13是用来说明由图4的统计性鉴定部进行的、关于层次贝叶斯模型及马尔科夫链蒙特卡罗法的负荷条件(血管内的平均压力)的事后分布计算和平均内压的鉴定的图。
图14是用来说明由图4的统计性鉴定部进行的、关于层次贝叶斯模型及马尔科夫链蒙特卡罗法的材料模型参数的事后分布计算和材料模型参数(血管壁的等价弹性率)的鉴定的图。
图15是表示由图1的显示设备进行的、作为力学性指标之一的内压的空间分布的显示例的图。
图16是表示由图1的显示设备进行的、作为血液流量指标之一的流速值的空间分布的显示例的图。
图17是由图1的显示设备显示的与左冠状动脉起始部的血压有关的曲线图。
图18是关于由图1的显示设备显示的LCX与LDA的分叉点附近的血压的曲线图。
图19是由图1的显示设备显示的关于与芯线方向有关的血压变化的曲线图。
图20是表示由图4的力学模型构建部进行的强制位移履历的其他分配例的图。
图21是表示由图4的力学模型构建部进行的强制位移履历的其他分配例的图。
图22是表示由图4的力学模型构建部进行的强制位移履历的其他分配例的图。
图23是表示有关本实施方式的纤维组的变形图。
图24是表示有关本实施方式的壁厚圆筒的力学模型的正交截面的图。
图25是图20的微小扇形要素的放大图。
具体实施方式
以下,参照附图说明有关本实施方式的血管解析装置、医用图像诊断装置及血管解析方法。
有关本实施方式的血管解析装置是用来对在通过医用图像诊断装置产生的医用图像中包含的血管区域进行构造流体解析的计算机装置。有关本实施方式的血管解析装置既可以装入在医用图像诊断装置中,也可以是与医用图像诊断装置分体的工作站等的计算机装置。以下,为了具体地进行说明,假设有关本实施方式的血管解析装置被装入在医用图像诊断装置中。
有关本实施方式的医用图像诊断装置在装备用来将被检体摄像的摄像机构的任何种类的图像诊断装置中都能够应用。作为有关本实施方式的医用图像诊断装置,例如可以适当地使用X射线计算机断层摄影装置(X射线CT装置)、磁共振诊断装置、超声波诊断装置、SPECT装置、PET装置、放射线治疗装置等。以下,为了具体地进行说明,假设有关本实施方式的医用图像诊断装置是X射线计算机断层摄影装置。
图1是有关本实施方式的医用图像诊断装置(X射线计算机断层摄影装置)的概略性的块结构图。如图1所示,X射线计算机断层摄影装置具有CT架台10和控制台20。CT架台10按照来自控制台20的架台控制部23的控制,用X射线将被检体的摄像部位摄像。摄像部位例如是心脏。CT架台10具有X射线管11、X射线检测器13及数据收集装置15。X射线管11和X射线检测器13可绕旋转轴Z旋转地装备在CT架台10上。X射线管11向被注入了造影剂的被检体照射X射线。X射线检测器13检测从X射线管11产生、透过了被检体的X射线,产生与检测到的X射线的强度对应的电信号。数据收集装置15从X射线检测器13读出电信号,变换为数字数据。将每1个视角的数字数据的集合称作原始数据集。与多个扫描时刻有关的时间序列的原始数据集通过非接触数据传送装置(未图示)被向控制台20传送。
控制台20以系统控制部21为中枢,具有架台控制部23、再构成装置25、图像处理装置27、输入设备29、显示设备31及存储装置33。
架台控制部23根据由用户经由输入设备29设定的扫描条件控制控制台20内的各装置。
再构成装置25基于原始数据集而产生关于被检体的CT图像的数据。具体而言,首先,再构成装置25对原始数据集实施前处理,产生投影数据集。作为前处理,包括对数变换及不均匀修正、校准修正等。接着,再构成装置25对投影数据集实施图像再构成处理而产生CT图像的数据。作为图像再构成算法,可以采用FBP(filtered back projection)法等的解析学的图像再构成法、ML-EM(maximum likelihood expectation maximization)法、OS-EM(ordered subset expectation maximization)法等的依次近似图像再构成等的已有的算法。在本实施方式中,再构成装置25基于时间序列的投影数据集而产生时间序列的CT图像的数据。CT图像包含关于通过造影剂造影的血管的像素区域(以下,称作血管区域)。另外,CT图像既可以是表现CT值的2维空间分布的切片数据,也可以是表现CT值的3维空间分布的体数据。以下,假设CT图像是体数据。将时间序列的CT图像的数据向存储装置33存储。
图像处理装置27基于时间序列的CT图像来构建力学模型,执行构造流体解析。关于图像处理装置27的处理的详细情况在后面叙述。
输入设备29受理来自用户的各种指令及信息输入。作为输入设备29,可以使用键盘、鼠标、开关等。
显示设备31将CT图像及构造流体解析结果等的各种信息显示。作为显示设备31,可以适当使用例如CRT显示器、液晶显示器、有机EL显示器、等离子显示器等。
存储装置33由硬盘装置等的各种存储介质构成。存储装置33将时间序列的投影数据及时间序列的CT图像数据等的各种数据存储。例如,存储装置33将时间序列的CT图像数据以依据DICOM(digital imaging and communications in medicine)标准的医用图像文件形式存储。此外,存储装置33也可以将由外部设备收集的医用数据与时间序列的CT图像数据在医用图像文件内建立关联而存储。
系统控制部21具有中央运算处理装置(CPU:central processing unit)及读出专用存储器(ROM:read only memory)、随机访问存储器(RAM:random access memory)。系统控制部21作为X射线计算机断层摄影装置的中枢发挥功能。系统控制部21执行存储在ROM或RAM中的血管解析程序,执行有关本实施方式的血管构造解析处理。
另外,系统控制部21、图像处理装置27、输入设备29、显示设备31及存储装置33构成血管解析装置50。如本实施方式那样,血管解析装置50既可以装入到医用图像诊断装置(X射线计算机断层摄影装置)中,也可以是与医用图像诊断装置分体的计算机装置。在血管解析装置50与医用图像诊断装置分体的情况下,血管解析装置50只要从医用图像诊断装置或PACS(picture archiving and communication systems)经由网络收集时间序列的CT图像等的医用数据就可以。
接着,详细地说明有关本实施方式的动作例。另外,有关本实施方式的血管解析装置、医用图像诊断装置及血管解析方法可以将心脏血管或颈动脉、脑动脉等的人体的所有的部位的血管作为解析对象。但是,以下为了具体地进行说明,假设有关本实施方式的解析对象是心脏的血管。
作为心脏的血管,可以举出例如冠状动脉和大动脉。冠状动脉从大动脉的冠状动脉起始部开始,在心肌表面上行进,从心外膜侧向内膜侧进入。冠状动脉在心肌的内膜中分叉为无数的毛细管。在分叉后,无数的毛细管再次合并而形成大心静脉,连接到冠状窦(日语:冠静脈洞)。冠脉血管系统与其他内脏不同,具有必须在心肌的收缩及松弛的力学性变化中保障灌流的特征。
作为冠脉血流的特征,相比在由心肌收缩带来的机械性血流阻碍作用下冠状动脉起始部的内压变高的收缩期,在左心室舒张期中当灌流压下降时流过较多。因此,正常的冠状动脉血流速波形是收缩期和舒张期的双峰性,舒张期血流是占优的。已知在肥大型心肌症或大动脉瓣狭窄症中,在收缩期中确认逆行性波,在大动脉逆流症中通过收缩期顺行波变大等疾病而呈现特异性的血流波形。此外,舒张期的顺行性波形与左心室舒张功能、特别是左心室松弛有密切的关系。在左心室松弛延迟例中,舒张期波形的尖峰向后偏移,此外有减速脚变平缓的趋向。此外,在这样的症例中,在快速心搏时舒张期的冠脉血流不能充分增大,可以考虑会助长心肌虚血。
在解剖学上,通过在从大动脉起始部开始分叉的左右冠状动脉上作用与大动脉压相等的冠脉灌流压、即通过作用冠状动脉分支的大动脉起始部的压力,产生冠脉血流。对于决定冠脉血流而言,作为大动脉压的驱动压和冠脉血管阻力是重要的。在140~180μm以上的较粗的冠脉血管中存在冠脉血管阻力的20%左右,相对于此,在100~150μm以下的微小血管中可以说较多存在阻力成分的残留。因而,在没有所谓冠脉狭窄等的情况下,阻力值受冠脉微小血管的紧张性(强直性)左右。
作为血管阻力因素,可以举出血管特性、动脉硬化、管狭窄、血液粘性、机械性因素。冠脉微小血管的强直性由血管特性、心肌代谢(心肌氧消耗)、神经体液性因素、机械性因素、作为体液因素的各种血管动作性物质、血液粘性规定,进而,也受到包括心肥大、冠状动脉硬化等在内的各种各样的病变影响,而引起冠脉循环障碍。
冠状动脉血流搏动受到冠状动脉血流的搏动样式、通过心肌收缩进行的心肌内血流的控制、心肌内血管对于机械性刺激的反应影响。作为心肌收缩阻碍血流的机理,可以举出心肌内压的上升、心肌内血管容量的变化、心肌内血管的压迫。在心肌舒张期的血流规定因素中,存在舒张期的冠状动脉压、舒张期的血管外力、心搏数、舒张期在心周期中所占的比例、心肌松弛。
有关本实施方式的血管解析装置50基于时间序列的CT图像来构建力学模型,利用力学模型来执行关于心脏的血管的构造流体解析,高精度地计算血管内的力学性指标及血管流量指标。为了计算精度较高的力学性指标及血管流量指标,需要对力学模型分配精度较高的潜在变量。血管解析装置50在构建力学模型时,对初始的力学模型实施逆解析,将潜在变量统计性地鉴定。由此,血管解析装置50能够决定高精度的潜在变量。力学性指标是指关于血管壁的力学性的指标。作为关于血管壁的力学性指标,例如分类为关于血管壁的位移的指标、关于在血管壁上发生的应力或应变的指标、关于在血管内腔上负荷的内压分布的指标、关于表示血管的硬度等的材料特性的指标等。关于表示血管的硬度等的材料特性的指标可以举出表示血管组织的应力与应变的关系的曲线的平均斜率等。血液流量指标是指与在血管中流动的血液有关的血行动态的指标。作为血管流量指标,例如可以举出血液的流量、血液的流速、血液的粘性等。
潜在变量例如包括血管的材料结构式或血液的材料结构式等的材料模型的参数(例如杨氏模量或泊松比等)、在血管内腔上负荷的内压分布等的负荷条件参数、构造解析或流体解析的边界条件参数、与时间序列的形态指标或形状变形指标的不确定性相关联的偏差分布参数的至少一个。这里,所谓与时间序列的形态指标或形状变形指标的不确定性相关联的偏差分布参数,是考虑在医用图像数据中存在由各CT值的噪声引起的偏差分布、及由生物体组织的边界阈值的含糊性引起的概率分布等、而将各种不确定性表现为概率分布。作为各种不确定性,例如可以举出血管组织及血液的边界坐标及特征点(血管分叉部及造影剂分散配置等)的空间坐标中的不确定性、几何学构造参数(与芯线垂直的截面的内腔半径等)的不确定性、医用图像数据自身(CT值及边界阈值等)的不确定性。
力学模型是用来表现血管及血液的动态的数值模型。力学模型按照构造流体解析的手法而具有不同的类型。例如,将力学模型分类为连续体力学模型和简单力学模型。连续体力学模型例如被用在有限要素法(FEM:finite element method)或边界要素法中。将简单力学模型例如分类为基于材料力学的材料力学模型和基于流动学的流体力学模型。另外,在以下的说明中没有特别言及的情况下,关于力学模型的类型不做特别限定。初始的力学模型是指被分配了与潜在变量的参数有关的采样集合(各参数的组合的集合)的力学模型,所述潜在变量的参数是根据潜在变量的概率分布及变量范围得到的。
图2是表示关于构造流体解析的对象区域(以下称作解析对象区域)的力学模型M1的一例的图。如图2所示,力学模型M1具有大动脉区域R1、右冠状动脉区域R2和左冠状动脉区域R3。血液从大动脉向右冠状动脉或左冠状动脉流动。
如图2所示,力学模型M1的大动脉起始部侧的末端被设定在血流的入口,右冠状动脉区域的末端和左冠状动脉区域的末端被设定在血流的出口。对于入口和出口分别设定边界条件。关于入口的边界条件例如包括入口处的血流的流速或由血流带来的压力或它们的变化率。关于出口的边界条件例如包括出口处的血流的流速或由血流带来的压力或它们的变化率。大动脉、右冠状动脉及左冠状动脉的变形依存于由血流引起的向血管壁的力学作用、由心脏的搏动带来的向血管壁的力学作用(外力)、血管截面边界的负荷条件、血管壁的材料模型、血管的无应力状态及血管壁的几何学形状等的各种各样的因素。这里,由血流引起的向血管壁的力学作用例如包括起因于血流的内压和起因于血流的剪切力。通过起因于血流的内压,在血管半径方向或血管内腔面的垂直方向上变形。由心脏的搏动带来的向血管壁的力学作用和起因于血流的剪切力以与血管芯线方向有关的伸缩及扭转、弯曲等的血管的变形为负荷条件,被分配给大动脉区域R1、右冠状动脉区域R2及左冠状动脉区域R3。具体而言,与血管芯线方向有关的伸缩及扭转、弯曲等的血管的变形通过强制位移(移动矢量或旋转位移)或载荷矢量的时间性变化来表现。此外,基于起因于血流的内压的、与血管半径方向或内腔面的垂直方有关向的变形被作为向血管内腔的压力分布的时间性变化而分配。
对于大动脉区域R1、右冠状动脉区域R2及左冠状动脉区域R3,在构造流体解析中分配由强制位移带来的位移约束条件。由此,能够将构造流体解析中的血管壁的变形自由度缩小,能够实现计算收敛性的稳定化及解析时间的缩短。
此外,例如血管的形状的变形程度依存于血管壁的材料。因此,对大动脉区域R1、右冠状动脉区域R2及左冠状动脉区域R3分别分配材料模型。此外,血管的形状的变形程度也依存于血管的无应力状态。作为负荷条件的初始值,也可以被分配血管的残留应力分布。在对血管解析对象的数值计算用力学模型的空间进行了离散化的节点集合及由节点构成的要素集合中,划分为将材料模型及边界条件、负荷条件等的解析条件进行鉴定的区域和不进行鉴定的区域。对于解析条件的鉴定区域内的节点,分配强制位移履历的位移约束条件,对于材料模型的鉴定区域,仅对血管壁表面(外表面)的节点分配强制位移履历的位移约束条件,对于血管壁内部确保位移自由度而不分配位移约束。由此,能够抑制构造流体解析的变形自由度,能够稳定且有效率地进行解析。但是,也可以在血管壁表面上设置缓冲用的虚拟要素集合,对其表面的节点赋予强制位移。由此,在血管内腔上有由斑块等形成的突起的情况下、或在血管分叉部等由内压带来的负荷给芯线方向的截面外的变形带来影响的情况下,可以不仅参照血管内腔、还参照血管壁的形态指标,来将负荷在血管上的载荷矢量和内压进行分离而鉴定。此外,在生理学上模拟壁表面的脂肪层,另一方面,还有在数值计算上通过对血管壁表面赋予强制位移而避免在血管壁内部局部地发生与实际不同的较高的应力的效果。
与这些材料模型、边界条件及负荷条件等的潜在变量有关的参数通过后述的基于力学模型的逆解析(统计性鉴定处理)来鉴定。将通过逆解析鉴定出的精确的潜在变量向力学模型分配。通过被分配了精确的潜在变量的力学模型,能够执行基于考虑了由解析对象血管区域外的血管或心脏等的外部因素带来的对该解析对象血管区域的影响的构造流体解析或流体解析或构造解析或图像解析的血行动态解析。血管解析装置50关于力学模型的构建,通过由逆解析进行的潜在变量的鉴定,能够解决以往例的以下的4点困难。困难1.冠状动脉的材料模型的鉴定方法。困难2.心脏的形状的变形向冠状动脉的影响的并入。困难3.冠状动脉的边界条件的鉴定方法。困难4.基于具有由医用图像数据的不确定性引起的偏差的血管形状的图像解析或构造流体解析。通过这4点困难的克服,血管解析装置50与不进行基于逆解析的潜在变量的鉴定的以往的血管构造流体解析相比,实现解析精度的提高。
接着,对有关本实施方式的构造流体解析处理的详细情况进行说明。图3是表示基于有关本实施方式的系统控制部21的控制进行的构造流体解析处理的典型的流程的图。图4是表示图像处理装置27的块结构的图。
如图3所示,在构造流体解析处理中,首先,通过系统控制部21从存储装置33将处理对象的医用图像文件读出,向图像处理装置27供给。医用图像文件除了时间序列的CT图像的数据以外,还包括关于该被检体的有关血管内腔的压力值的数据、血液流量指标的观测值的数据及斑块指标。在CT图像以外,也可以是MRI图像或超声波回波图像。时间序列的CT图像的数据是表现时间序列的CT值的3维空间分布的数据。时间序列的CT图像例如在1心搏下包括约20张、即约20心相位的CT图像。
如图3所示,系统控制部21使图像处理装置27进行区域设定处理(步骤S1)。在步骤S1中,图像处理装置27的区域设定部51对在时间序列的CT图像中包含的血管区域设定构造流体解析的解析对象区域。解析对象区域被设定在关于冠状动脉的血管区域的任意的一部分。例如区域设定部51通过用户经由输入设备29的指示或图像处理,对血管区域设定解析对象区域和鉴定对象区域。
这里,参照图5对血管的构造进行说明。图5是示意地表示与血管的芯线正交的截面(以下称作血管截面)的图。如图5所示,血管具有管状的血管壁。将血管壁的中心轴称作芯线。将血管壁的内侧称作内腔。血液流到内腔中。将内腔与血管壁的边界称作血管内壁。心肌等的血管周边组织分布在血管壁的外侧。将血管壁和血管周边组织的边界称作血管外壁。有在血管壁内部发生斑块的情况。如图5所示,将斑块例如分类为石灰化斑块或粥状斑块等。也有将粥状斑块称作不稳定斑块的情况。粥状斑块较柔软,有血管内壁破裂而作为血栓渗出到血管内部的危险性。因而,掌握斑块的性状在临床上是有用的。斑块的性状及存在区域可以通过医用图像文件中包含的斑块指标确定。斑块指标能够通过以例如骨的CT值为基准进行了标准化的CT值的大小等来相对地判别。但是,对血管内部的斑块的变形特性及硬度进行解析并不容易。
如果进行步骤S1,则系统控制部21使图像处理装置27进行图像解析-跟踪处理(步骤S2)。在步骤S2中,图像处理装置27的图像解析-跟踪处理部53对时间序列的CT图像实施图像处理,计算时间序列的血管形态指标和时间序列的血管形状变形指标。具体而言,图像解析-跟踪处理部53通过对时间序列的CT图像实施图像解析处理而计算时间序列的血管形态指标,通过对时间序列的CT图像实施跟踪处理而计算时间序列的血管形状变形指标。
更详细地讲,图像解析-跟踪处理部53在图像解析处理中,从各CT图像中提取血管区域,确定与血管的内腔有关的像素区域(以下称作血管内腔区域)和与血管壁有关的像素区域(以下称作血管壁区域)。图像解析-跟踪处理部53作为血管形态指标而确定垂直于血管的芯线的截面或垂直于血管内腔面的面与血管内腔、血管壁、斑块区域交叉的区域上的多个像素的3维坐标。另外,血管形态指标并不限定于3维坐标,也可以是各种几何学指标。例如,作为有关本实施方式的几何学指标,可以利用与芯线垂直的截面中的每一定角度的血管内腔的半径或直径及0°的方向矢量、或者截面中的与全部角度相对的平均面积或平均半径、或者由与芯线方向垂直的多个截面包围的血管内腔容积、或者由与内腔面垂直的多个截面包围的血管壁容积或斑块容积等。
在跟踪处理中,图像解析-跟踪处理部53通过来自用户的经由输入设备29的指示或图像处理来设定特征点和/或特征形状、代表点、像素等的多个跟踪点。跟踪点具体而言被设定在血管区域和/或血液区域、造影剂区域、质子区域等中。例如,图像解析-跟踪处理部53设定血管分叉部或表面的特征形状等的跟踪点集合。根据通过各时刻(各心相位)的图像解析-跟踪处理部53的跟踪处理而得到的跟踪点集合的位移数据,通过插补处理等来计算力学模型的血管壁表面、血管壁内部、或血管内腔中的节点的位移的时间性变化,将计算出的时间性变化作为强制位移而进行分配。此外,例如图像解析-跟踪处理部53对力学模型定义血管芯线上的节点。图像解析-跟踪处理部53也可以根据力学模型的血管壁表面、血管壁内部或血管内腔中的节点的位移的时间性变化,提取与血管的关于芯线方向的伸缩、扭转、弯曲有关的变形,作为与血管芯线和芯线垂直的截面中的节点的强制位移给予来表现。这样,作为血管形状变形指标,确定力学模型中的各时刻的节点的强制位移数据(强制位移履历)。
以下,参照图6、图7及图8,说明图像解析-跟踪处理。图6是表示血管芯线的形态的时间性变化的图。如图6所示,例如假设时间序列的医用图像每1心搏包含20张CT图像。即,假设从心相位0%到95%每隔5%得到CT图像。由图像解析-跟踪处理部53从各CT.图像中提取血管区域的芯线。如图6所示,芯线的形态随着心相位的经过而变化。
图8是表示时刻t与时刻t+Δt之间的跟踪处理的一例的图。如图8所示,在血管芯线上设定了从P1到P10的力学模型中的节点,与各截面上的血管的力学模型的节点在力学上相连。但是,与血液的力学模型的节点是独立的。基于血管的跟踪点的位移数据,通过插补等的处理而计算血管芯线上的从P1到P20的节点的位移数据,对各节点设定了强制位移。由于对血管形状变形指标和血管形态指标进行说明,所以考虑由节点P13和节点P14规定的局部血管区域RA。假设在时刻t,与芯线方向有关的节点P13与节点P14之间的距离是L,血管区域的半径是r。通过从图像解析-跟踪处理部53提取节点P13和节点P14的血管芯线方向的伸缩、扭转、弯曲等的强制位移,计算节点P13的强制位移(3维空间中的移动位移和芯线方向的旋转位移)和节点P14的强制位移(3维空间中的移动位移和芯线方向的旋转位移)。
图7是表示血管芯线的弯曲变形或旋转位移的计算例的图。如图7所示,例如,扭转角也可以根据由线a及b构成的面的法线方向矢量c的变化来计算。
如图7及图8所示,图像解析-跟踪处理部53基于跟踪点的坐标和移动矢量,计算芯线上的各节点的强制位移(3维空间中的移动位移和芯线方向的旋转位移),计算血管形状变形指标。例如,图像解析-跟踪处理部53计算相邻的两个节点的坐标差的时间变化来作为与芯线方向有关的伸缩距离ΔL。此外,图像解析-跟踪处理部53关于各芯线上的节点,计算该节点与包含该节点在内的血管区域截面上的其他节点(血管内腔或血管壁或斑块区域中的节点)之间的距离的时间变化,作为与半径方向有关的伸缩距离Δr。此外,图像解析-跟踪处理部53关于各跟踪点,基于该跟踪点的附近的多个跟踪点的坐标和移动矢量,计算芯线上的该节点的芯线方向的扭转角度Δθ。
此外,图像解析-跟踪处理部53也可以通过血液区域的造影剂或质子的图像跟踪,作为血液流量指标而计算流速或芯线方向截面的平均流速或平均流量。
将血管形状变形指标作为力学模型中的强制位移使用。以下,将时间序列的血管形态指标称作形状履历,将时间序列的血管形状变形指标称作强制位移履历。
如果进行步骤S2,则系统控制部21使图像处理装置27进行构建处理(步骤S3)。在步骤S3中,图像处理装置27的力学模型构建部55基于形状履历(时间序列的血管形态指标)、强制位移履历(时间序列的血管形状变形指标)和时间序列的医用图像(CT图像、MRI图像、超声波回波图像等的DICOM数据),暂定地构建关于解析对象区域的力学模型。力学模型是与用来进行构造流体解析的解析对象区域有关的数值模型。
以下,对步骤S3详细地说明。力学模型构建部55首先基于医用图像和形状履历,构建用来解明力学模型(数理模型)的形状模型。形状模型是将各时刻的血管区域的几何学构造示意性地表现的模型。形状模型例如被划分为多个离散化区域。将各离散化区域的顶点称作节点。力学模型构建部55既可以基于各时刻的医用图像中包含的血管区域和血管形态指标来构建各时刻的形状模型,也可以基于在特定的时相的医用图像中包含的血管区域和血管形态指标来构建各时刻的形状模型。此外,例如在作为初始的负荷状态而假定在与解析对象区域对应的血管中不存在残留应力的情况下,作为无应力状态的时相而将与解析对象区域对应的血管最收缩的时相假定是无应力状态。
图9是表示与形状模型的芯线正交的截面的图。如图9所示,形状模型从芯线朝向外侧具有血管内腔区域、血管壁区域。在斑块存在的情况下,也可以在血管壁区域中设置斑块区域。此外,在考虑由血管周边组织带来的向血管的影响的情况下,也可以将血管周边组织的虚拟要素设在血管壁区域的外侧。
如果构建形状模型,则力学模型构建部55将根据各潜在变量的概率分布或变量范围得到的关于潜在变量的参数的采样值(例如基于马尔科夫链蒙特卡罗法等的、来自各参数的组合的集合的采样)向力学模型设定。例如,力学模型构建部55如图2所示,对大动脉区域R1的大动脉起始部侧的末端设定与入口有关的边界条件的鉴定对象的区域(以下称作边界条件鉴定区域),对右冠状动脉区域R2的末端和左冠状动脉区域R3的末端设定与出口有关的边界条件鉴定区域。力学模型构建部55向各边界条件鉴定区域分配根据边界条件的概率分布和/或变量范围得到的与边界条件的参数有关的采样值。此外,力学模型构建部55对大动脉区域R1、右冠状动脉区域R2及左冠状动脉区域R3设定材料模型的鉴定对象的区域(以下称作材料模型鉴定区域)及负荷条件的鉴定对象的区域(以下称作负荷条件鉴定区域)。力学模型构建部55向各材料模型鉴定区域分配根据材料模型的概率分布和/或变量范围得到的关于材料模型的参数的采样值,向各负荷条件鉴定区域分配根据负荷条件的概率分布及变量范围得到的与负荷条件的参数有关的采样值。血管据说即使流量是0也存在残留应力。例如,力学模型构建部55也可以将流量为0的情况下的残留应力作为负荷条件的初始值而向解析对象区域分配。此外,力学模型构建部55也可以对在几何学构造上有不确定性部位设定几何学构造的鉴定对象的区域(以下称作几何学构造鉴定区域)。另外,几何学构造的参数是与几何学构造的不确定性相关联的偏差分布参数、或内在于医用图像数据中的偏差分布参数,也可以是各CT值的偏差分布或生物体组织的边界阈值的偏差分布等。详细后述,但力学模型构建部55也可以对斑块区域设定材料模型。关于材料模型的详细情况后述。
如果构建出形状模型,则力学模型构建部55将在步骤S2中计算出的时间序列的血管形状变形指标即强制位移履历向形状模型分配。将被分配了潜在变量及强制位移履历的形状模型称作力学模型。
图10是形状模型M2表示血管及血液的力学模型的一部分、用来说明向力学模型的节点的强制位移履历的分配的图。图10表示形状模型M2的一部分。但是,在图10中表示芯线位于M2内的情况,但也可以是芯线位于M2外的情况。如图10所示,对于形状模型M2设定了多个节点PN(PN1、PN2)。将芯线上的节点称作PN1,将血管及血液的力学模型中的节点称作PN2。形状模型M2被设定在虚拟要素表面、血管外壁、血管内壁、斑块区域表面、斑块区域内部或血液区域。力学模型构建部55按照每一时刻对形状模型M2的各节点PN1分配强制位移、即血管形状变化指标。
具体而言,力学模型构建部55将芯线上的相邻的节点PN1和节点PN1用横梁要素(或刚性要素)EB连结。此外,力学模型构建部55将节点PN1和在穿过该节点PN1的正交截面中包含的其他节点PN2用横梁要素EB连结。力学模型构建部55向节点PN1及横梁要素EB分配与各血管形状变形指标的形状位移方向有关的约束条件。在对材料模型或血管内腔的内压进行鉴定的区域中,作为强制位移,可以举出与芯线方向有关的血管壁(或虚拟要素)表面的伸缩、血管壁(或虚拟要素)表面的扭转、血管壁(或虚拟要素)表面的弯曲变形。例如,在不鉴定材料模型或血管内腔的内压的区域中,不仅是与芯线方向有关的强制位移,将与半径方向有关的血管壁的时间序列的伸缩(位移)也作为强制位移履历分配。此外,在血管内腔有突起的情况或如血管分叉部等那样内压贡献于向芯线方向截面外的变形的情况下,不对该区域分配强制位移履历,而仅对该区域的周边部(例如虚拟要素的表面节点)分配强制位移履历。此外,力学模型构建部55对节点PN1及横梁要素EB分配时间序列的血管形状变形指标作为强制位移履历。由此,表现与血管的整体及局部有关的伸缩变形、扭转变形及弯曲变形。
另外,强制位移履历的分配对象并不限定于芯线上的节点及横梁要素。图11及图12是表示向形状模型的强制位移履历的其他分配方法的图。图11表示血管形状变形指标为扭转的情况下的分配例,图12表示血管形状变形指标为弯曲的情况下的分配例。如图11及图12所示,力学模型构建部55也可以对形状模型的表面或内部的节点直接分配强制位移履历。例如,在步骤S2中,图像解析-跟踪处理部53计算关于特征点的伸缩量、扭转量及弯曲量等的血管形状变形指标。力学模型构建部55通过将计算出的血管形状变形指标向特征点的周围的节点进行插补(内插或外插)而直接分配。
有关本实施方式的图像处理装置27使用在步骤S3中暂定地构建的力学模型实施逆解析,将对力学模型设定的潜在变量统计性地鉴定。统计性鉴定处理在后述的步骤S6中进行。步骤S4及S5为了分别计算在统计性鉴定处理中使用的血管形态指标及血液流量指标而设置。
如果进行步骤3,则系统控制部21使图像处理装置27进行血管应力解析处理(步骤S4)。在步骤S4中,图像处理装置27的血管应力解析部57对现阶段的力学模型实施血管应力解析,计算时间序列的血管形态指标的预测值。血管形态指标是已述的血管形态指标的哪种都可以,例如可以使用与血管芯线方向相关的内腔区域的截面形状指标或血管壁的截面形状指标。具体而言,内腔区域的截面形状指标是内腔区域的关注像素的坐标值、内腔区域的几何学构造参数(内腔区域的半径、内腔区域的直径等)的至少一个。此外,血管壁区域的截面形状指标具体而言是血管壁区域的关注像素的坐标值、血管壁区域的几何学构造参数(血管壁区域的半径、壁区域的直径等)的至少一个。另外,预测值是指对力学模型实施血管应力解析而计算出的血管形态指标的计算值。
此外,如果进行了步骤3,则系统控制部21使图像处理装置27进行血液流体解析处理(步骤S5)。在步骤S5中,图像处理装置27的血液流体解析部59对暂定地构建的力学模型实施血液流体解析而计算时间序列的血液流量指标的预测值。血液流量指标是血流量或流速。或者,血液流量指标也可以是血流量或流速的空间性或时间性的平均值。另外,预测值是指对力学模型实施血液流体解析而计算出的血液流体指标的计算值。
如果进行了步骤S4及S5,则系统控制部21使图像处理装置27进行鉴定处理(步骤S6)。在步骤S6中,图像处理装置27的统计性鉴定部61对力学的模型的潜在变量的参数进行统计性地鉴定,以使在步骤S4中计算出的血管形态指标的预测值和在步骤S5中计算出的血液流量指标的预测值与事前收集到的血管形态指标的观测值和血液流量指标的观测值匹配。
如图4所示,统计性鉴定部61具有第1统计性鉴定部61-1和第2统计性鉴定部61-2。第1统计性鉴定部61-1统计性地鉴定潜在变量的参数,以使血管形态指标的预测值与血管形态指标的观测值匹配。第2统计性鉴定部61-2统计性地鉴定潜在变量的参数,以使血液流量指标的预测值与血液流量指标的观测值匹配。以下,依次说明第1统计性鉴定部61-1和第2统计性鉴定部61-2。
具体而言,在步骤S6中,第1统计性鉴定部61-1设定基于在步骤S4中计算出的血管形态指标的预测值和观测值的数据分布。数据分布例如表示与血管形态指标的预测值与观测值的误差有关的多变量正规分布函数。具体而言,统计性鉴定处理部61-1对于力学模型中的各节点或各要素,计算与预测值与观测值之间的误差有关的正规分布函数值,将各正规分布函数值的积设定为数据分布。数据分布既可以按照每一时刻单独地设定,也可以多个时刻一起设定。接着,第1统计性鉴定部61-1向力学模型的潜在变量分配事前分布(事前概率分布)。具体而言,向与材料模型、边界条件、负荷条件及时间序列的形态指标及形状变形指标的不确定性关联的参数分别分配事前分布。例如,分配与作为负荷条件的参数之一的血管内腔相关的压力值有关的事前分布。压力值的可取的值的范围(设想范围)可以在经验上预先限定。第1统计性鉴定部61-1通过限定于这些设想范围内而执行关于内压值的蒙特卡洛模拟,从而对于各离散化区域计算内压值的概率分布、即事前分布。此外,力学模型构建部55作为事前分布,也可以在观测到芯线方向的压力分布较平滑、此外伴随着时间经过的压力变化较平滑、没有血流的倒流的情况下,将芯线方向的平均的压力变化的斜率为负、例如将通过多变量正规分布函数在数学上表现的概率分布,设定为事前分布。按照设定在这些设想范围内的概率分布,能够执行关于负荷条件的参数的蒙特卡洛模拟,能够得到用来向力学模型设定的负荷条件(潜在变量)的采样值。接着,第1统计性鉴定部61-1关于各潜在变量,通过对事前分布和数据分布实施统计性鉴定处理来计算事后分布(事后概率分布)。统计性鉴定处理例如可以举出层次贝叶斯模型或马尔科夫连锁模型。并且,第1统计性鉴定部61-1关于各潜在变量,根据事后分布的众数或平均值等的统计值来对各潜在变量的参数进行鉴定。例如,在上述例子的情况下,计算关于血管内腔压力值的事后分布,根据该事后分布来计算血管内腔压力值的鉴定值。
图13是用来说明通过层次贝叶斯模型及马尔科夫链蒙特卡罗法进行的关于负荷条件(血管内的平均压力)的事后分布计算和平均内压的鉴定的图。如图13所示,假设对于从血管起始部延伸的血管设定了石灰化斑块区域和粥状斑块区域。石灰化斑块区域被设定为材料模型鉴定区域,粥状斑块区域被设定为材料模型鉴定区域。随着从血管起始部沿着血管芯线方向行进而血管内压下降。沿着血管芯线设定多个节点。在包含各节点的正交截面(节点截面)中计算内腔内压的事后分布,鉴定事后分布的众数。
另外,作为血管形态指标的观测值,例如使用在步骤S2中计算出的血管形态指标。
第2统计性鉴定部61-2的处理只是在数据分布的计算中使用的指标不同,与第1统计性鉴定部61-1的处理是同样的。即,第2统计性鉴定部61-2首先设定基于在步骤S5中计算出的血液流量指标的预测值和观测值的数据分布。接着,第2统计性鉴定部61-2向力学模型的潜在变量分配事前分布。例如,分配与关于血管的材料模型的参数及关于血液的材料模型的参数、关于斑块的材料模型的参数有关的事前分布。作为这些材料模型的参数,例如可以举出弹性率等的材料模型参数、或血液的构成式中的关于粘性的参数。材料模型的参数的设想范围及概率分布可以在经验上预先设定。第2统计性鉴定部61-2关于各离散化区域而设定材料模型的参数的概率分布、即事前分布,能够按照限定于这些设想范围内的概率分布,执行关于材料模型的参数的蒙特卡洛模拟,能够得到用来向力学模型设定的材料模型参数(潜在变量)的采样值。接着,第2统计性鉴定部61-2对于各潜在变量,通过对事前分布和数据分布实施统计性鉴定处理而计算事后分布,根据计算出的事后分布的统计值来鉴定各潜在变量的参数。例如,在上述例子的情况下,计算关于材料模型的参数的事后分布,根据该事后分布来计算材料模型的参数的鉴定值。
图14是用来说明通过层次贝叶斯模型及马尔科夫链蒙特卡罗法进行的关于材料模型参数的事后分布计算和材料模型参数(血管壁的等价弹性率)的鉴定的图。如图14所示,假设血管模型与图13是同样的。限定于材料模型鉴定区域而计算血管壁的材料模型的参数(例如等价弹性率)的事后分布,鉴定事后分布的众数。
另外,血液流量指标的观测值例如假定是向大动脉送出的血流量变化,可以使用血管形态指标的观测值作为从时间序列的CT图像中通过图像处理计测的左心室的容积变化值(CFA)。也可以通过由造影剂的冠状动脉内注入后的造影剂的图像跟踪来计算特征点的移动量的时间性变化,从而计算流速或流量。此外,也可以取得造影剂的血管芯线方向或时间性的特定区域的浓度变化量,根据用各区域的芯线方向距离间隔除以该浓度变化后的值或浓度变化的时间性变化率,来计算流速或流量。在MRI的情况下使用质子的图像跟踪,在超声波回波的情况下通过对比回波图法等来计算流量。
此外,在不以解析对象区域的各像素的坐标值是确定值为前提的情况下,即假定在解析对象区域的几何学构造中有不确定性的情况下,也可以在潜在变量中包含几何学构造。在此情况下,统计性鉴定部61可以将表现出各节点的坐标值的关于芯线方向的规定范围内的变动、或解析对象区域的径的规定范围内的变动的正规分布等的概率分布设定为事前分布。在此情况下,也可以将解析对象区域的形状是平滑的、此外芯线中的节点的顺序是不变的这样的制约设定为事前分布。按照限定在设想范围内的概率分布,能够执行关于几何学构造的参数的蒙特卡洛模拟,能够得到用来向力学模型设定的几何学构造的不确定性参数(潜在变量)的采样值。
另外,在各步骤S6中也可以不进行第1统计性鉴定部61-1的统计性鉴定处理和第2统计性鉴定部61-2的统计性鉴定处理这两者。即,在各步骤S6中,也可以进行第1统计性鉴定部61-1的统计性鉴定处理和第2统计性鉴定部61-2的统计性鉴定处理的某一方。
此外,在上述例子中,第1统计性鉴定部61-1统计性地鉴定潜在变量的参数,以使血管形态指标的预测值与血管形态指标的观测值匹配,第2统计性鉴定部61-2统计性地鉴定潜在变量的参数,以使血液流量指标的预测值与血液流量指标的观测值匹配。但是,统计性鉴定部61也可以基于构造-流体复合解析,统计性地鉴定潜在变量的参数,以使血管形态指标的预测值和血液流量指标的预测值与血管形态指标的观测值和血液流量指标的观测值匹配。关于由统计性鉴定部61进行的统计性鉴定处理的更详细的情况在后面叙述。
如果进行了步骤S6,则图像处理装置27的力学模型构建部55将在步骤S6中计算出的潜在变量的参数向力学模型设定。
并且,系统控制部21判定是否满足鉴定结束条件(步骤S7)。在步骤S7中判定为不满足鉴定结束条件的情况下(步骤S7:NO),系统控制部21反复步骤S4、S5、S6及S7。这里,鉴定结束条件通过用来判定鉴定结束的指标(以下称作鉴定结束指标)是否达到规定值来表现。作为鉴定结束指标,例如可以举出血管形态指标的预测值与观测值的差值。在此情况下,系统控制部21在该差值比既定值大的情况下判定为不满足鉴定结束条件,在差值比既定值小的情况下判定为满足鉴定结束条件。此外,鉴定结束指标例如也可以是蒙特卡罗法的采样点的数量。在此情况下,系统控制部21在该采样点的数量比既定值小的情况下判定为不满足鉴定结束条件,在采样点的数量比既定值大的情况下判定为满足鉴定结束条件。在满足鉴定结束条件的情况下,力学模型构建部55将该时点的最新的力学模型设定为最终的力学模型。
如果构建了最终的力学模型,则力学模型构建部55计算将血管形状变形指标的观测值、最终的力学模型中的负荷条件的参数及材料模型的参数建立了关联的模型(以下称作关联模型)。将关联模型向存储装置33存储。关联模型为了检索的容易性等,可以与患者信息或检查信息等建立关联而存储。另外,血管形态指标或血液流量指标的观测值、最终的力学模型中的负荷条件的参数及材料模型的参数并不一定需要以模型的形态建立关联,例如也可以是表或数据库。
上述的步骤S4、S5、S6及S7既可以通过相同的鉴定法反复,也可以通过不同的鉴定法反复。在通过不同的鉴定法反复的情况下,例如也可以首先利用简单力学模型暂定地鉴定潜在变量,接着利用连续体力学模型正确地鉴定潜在变量。通过这样将统计性鉴定处理用不同的手法分为两阶段进行,能够使潜在变量的参数在短时间收敛。作为利用简单力学模型的方法,对于内压及外压可以举出厚壁圆筒的材料力学式。此外,作为利用简单力学模型的方法,也可以使用哈根-泊肃叶流动及修正伯努利的式子。作为利用连续体力学模型的方法,可以举出FEM构造流体解析。关于利用简单力学模型的鉴定法和利用连续体力学模型的鉴定法的详细情况在后面叙述。
在步骤S7中判定为满足鉴定结束条件的情况下(步骤S7:YES),系统控制部21可以使图像解析-跟踪处理部53进行修正处理(步骤S8)。在步骤S8中,图像解析-跟踪处理部53可以将时间序列的医用图像中包含的血管区域的形状修正,以使基于通过统计性鉴定法的逆解析得到的潜在变量而实施的构造流体解析结果(力学性指标的预测值及血液流体指标的预测值)与观测值(力学性指标的观测值及血液流体指标的观测值)匹配。显示设备31将基于修正后的时间序列的医用图像的诊断结果显示。由此,血管解析装置50能够将考虑到最终的力学模型的诊断结果显示。或者,显示设备31也可以将逆解析的鉴定和通过该构造流体解析与观测结果不匹配的血管部位/区域显示到画面上。例如,血管的动态的运动较快的心相位的图像模糊的情况较多,在通过基于医用图像的图像解析而观测到的血管形状中存在误差较大的部位或区域。血管的动态比较稳定的心相位的医用图像数据噪声较少。可以基于误差分布较小的血管形状数据将误差较大的心相位的血管形状通过使用力学模型来正确地内插。关于误差较大的血管部位或区域,能够将正确地内插的形状与自源数据起的偏差分布一起显示。由此,能够确保血管形状显示的稳定性,并且用户能够认识到形状的不确定性。
如果进行了步骤S8,则系统控制部21使图像处理装置27进行血管应力解析处理(步骤S9)。在步骤S9中,图像处理装置27的血管应力解析部57对最终的力学模型实施血管应力解析,计算时间序列的力学性指标的预测值的空间分布。具体而言,按照每个离散化区域来计算力学性指标的预测值。
此外,如果进行了步骤S8,则系统控制部21使图像处理装置27进行血液流体解析处理(步骤S10)。在步骤S10中,图像处理装置27的血液流体解析部59对暂定地构建的力学模型实施血液流体解析而计算时间序列的血液流量指标的预测值的空间分布。具体而言,按照离散化区域计算血液流量指标的预测值。
另外,作为力学性指标或血液流量指标,也可以计算FFR。
如果进行了步骤S9及S10,则系统控制部21使显示设备31进行显示处理(步骤S11)。在步骤S11中,显示设备31将在步骤S9中计算出的时间序列的力学性指标的预测值和在步骤S10中计算出的时间序列的血液流量指标的预测值显示。例如,显示设备31将时间序列的力学性指标或时间序列的血管流量指标在时间序列的力学模型中以与该预测值对应的颜色动态地显示。因此,显示设备31保持有表示各种预测值与颜色值(例如RGB)的关系的颜色表。显示设备31利用颜色表确定与预测值对应的颜色值,用与所确定的颜色值对应的颜色显示与该预测值对应的离散化区域。
图15是表示作为力学性指标之一的内压的空间分布的显示例的图。如图15所示,显示设备31将构成力学模型的各离散化区域以与关于该离散化区域的内压值相应的颜色,动态地显示。用户通过观察力学模型,能够利用颜色来掌握随着时间经过且空间性变化的力学性指标。
图16是表示作为血液流量指标之一的流速值的空间分布的显示例的图。如图16所示,显示设备31将构成力学模型的各离散化区域以与关于该离散化区域的流速值相应的颜色,动态地显示。用户通过观察力学模型,能够利用颜色来掌握随着时间经过且空间性变化的血液流量指标。
例如,在血管内完全狭窄的情况下,狭窄部位的内压比非狭窄部位的内压小。在作为力学性指标而指定了内压的情况下,用户能够根据力学模型上的局部性的颜色的差异来判断狭窄的有无。此外,狭窄部位的流量比非狭窄部位的流量小。在作为血液流量指标而指定了流量的情况下,用户能够根据力学模型上的局部性的颜色的差异来判断狭窄的有无。
此外,血管应力解析部57也可以基于斑块区域的材料模型参数的鉴定结果,作为力学性指标而计算硬度值的空间分布。在此情况下,显示设备31也可以将关于斑块区域的硬度值的空间分布在力学模型上显示。此外,显示设备31也可以对斑块区域周边的内压分布或应力分布、应变分布进行显示。用户能够将这些显示用于推测斑块的性状和破裂容易度。
力学性指标及血液流量指标的预测值并不仅限定于用力学模型的离散化区域的颜色表现的方法。例如,如图17、图18及图19所示,也可以用曲线图显示。另外,图17是关于左冠状动脉起始部的血压的曲线图。图17的曲线图的纵轴被规定为标准化的血压,横轴被规定为心相位[%]。图18是关于LCX与LDA的分叉点附近的血压的曲线图。图18的曲线图的纵轴被规定为标准化的血压,横轴被规定为心相位[%]。图19是与关于芯线方向的血压变化有关的曲线图。图19的曲线图的纵轴被规定为血压比,横轴被规定为距大动脉的距离[mm]。显示设备31通过将力学性指标及血液流量指标的预测值用曲线图显示,能够使用户简便地掌握这些值。
如果进行了步骤S11,上述构造流体解析处理结束。
另外,在图10中,将强制位移履历对形状模型的芯线部和外壁部设定,但强制位移履历的设定部位并不限定于此。例如,也可以将强制位移履历对芯线部与外壁部之间的血管壁区域设定。
另外,强制位移履历的约束条件的分配对象可以根据是否鉴定边界条件及材料模型来进行划分。图20是表示强制位移履历的其他分配例的图,表示形状模型的截面。例如,如图20(a)所示,在鉴定边界条件及材料模型的情况下,可以仅对形状模型的外壁部OW上的节点PN2分配强制位移履历,不对血管壁区域RV的节点PN3分配强制位移履历。如图20(b)所示,在不鉴定边界条件及材料模型的情况下,可以对形状模型的外壁部的节点PN2和血管壁区域RV的节点PN3这两者分配强制位移履历。在此情况下,对芯线上的节点PN1分配强制位移履历。此外,也可以将外壁部OW的节点PN2和节点PN1用横梁要素EB连结,对横梁要素EB上的节点PN2及PN3也分配强制位移履历。此时,将与半径方向有关的收缩及膨胀用横梁要素EB的伸缩位移表现。另外,也可以不对血管内腔区域RI分配强制位移履历。
图21是表示强制位移履历的分配的其他例的图,表示包含血管周边组织的虚拟要素RD在内的形状模型的截面。如图21所示,将虚拟要素RD设定在血管壁区域RN的外侧。在形状模型包含虚拟要素RD的情况下,除了血管壁区域RN以外,对虚拟要素RD也设定节点PN4。对于节点PN4也分配强制位移履历。力学模型构建部55也可以在鉴定边界条件及材料模型的情况下对血管壁区域RV中包含的节点PN3分配强制位移履历、在不鉴定边界条件及材料模型的情况下不分配强制位移履历。在对节点PN3分配强制位移履历的情况下,除了与内腔区域RI有关的形状指标以外,还参照关于血管壁区域RV的形状指标而进行材料模型的鉴定。
图22是表示强制位移履历的分配的其他例的图,表示包含斑块区域RP在内的形状模型的截面。如图22所示,斑块区域RP包含在血管壁区域RV中。斑块区域RP被设定为材料模型鉴定区域。关于斑块区域RP,考虑内腔形状指标、血管壁形状指标及斑块指标来鉴定材料模型。如已述那样,斑块指标例如是通过超声波诊断装置的组织性状诊断得到的关于斑块的性状的数据。力学模型构建部55根据性状将斑块区域划分为多个部分区域,对多个部分区域单独设定材料模型鉴定区域。可以对各部分区域预先设定与该部分区域的性状相应的参数范围。通过已述的统计性鉴定处理,将关于各部分区域的材料模型参数通过统计性鉴定部61进行鉴定。并且,在步骤S10中通过显示设备31作为力学性指标而显示关于血管的材料特性的指标,用户能够正确且容易地掌握斑块的性状。
接着,对作为潜在变量之一的材料模型详细地说明。作为血管的材料模型,可以使用弹性模型、超弹性模型、各向异性超弹性模型、考虑到粘性特性的超弹性模型等。作为各向异性超弹性模型,例如可以使用由Y.C.Fun提出的数理模型、称作Holzapfel-Gasser结构式的数理模型。每单位参照体积的应变能量用以下的(1)式表示。(1)式的第1项表示关于不包含骨胶原的各向同性基础材料的剪切变形的能量,第2项表示关于不包含骨胶原的各向同性基础材料的体积变形的能量,第3项表示骨胶原纤维各组的贡献(考虑纤维方向的分散)。
[数学式1]
图23是表示纤维组的变形的图。如图23所示,假定圆筒形上的外膜。将由芯线方向z和周向θ规定的面中的平均方向A上的纤维组的变形用下述的(2)式表示。
[数学式2]
作为(1)及(2)式中的关于材料模型的参数,如以下的表1所示,有材料参数和纤维分散参数。作为材料参数,使用C10、D、K1、K2等,作为纤维分散参数,使用Kappa、γ等。各参数的默认值和制约条件是表1所示那样的。
[表1]
血液的材料模型优选的是以下的(3)式那样的Casson结构式和(4)式那样的HB结构式。
[数学式3]
[数学式4]
另外,τ表示剪切应力,τ0表示屈服应力,μ表示粘性系数,表示剪切应变速度。
这些材料模型的参数,在上述步骤S6中由统计性鉴定部61通过使用了血管形态指标及血液流体指标的统计性鉴定处理来鉴定。
接着,说明由统计性鉴定部61进行的统计性鉴定处理的详细情况。
根据时间序列的医用图像计测的血管形态指标及血液流量指标那样的观测变量具有不确定性。统计性鉴定部61作为存在这样的不确定性的状况下的潜在变量的统计性鉴定法,而利用基于层次贝叶斯模型和马尔科夫链蒙特卡罗法的统计性方法。
如已述那样,在步骤S6中,统计性鉴定部61设定基于在步骤S4中计算出的血管形态指标或血液流量指标的预测值和观测值的数据分布。数据分布例如表示与血管形态指标或血液流量指标的预测值和观测值的误差有关的多变量正规分布函数。数据分布既可以按照每个时刻单独地设定,也可以多个时刻一起设定。接着,统计性鉴定部61向形状模型的强制位移和潜在变量分配事前分布。事前分布表示能够取的值的概率分布。接着,统计性鉴定部61执行关于潜在变量的数值模拟的参数勘察,构建表现潜在变量与血管形态指标或血液流量指标的关系的模型。例如,材料模型参数与内压分布参数与血管形态指标或血液流量指标的关系被模型规定。另外,血管形态指标或血液流量指标和潜在变量的关系也可以不是由模型的形态,而是由数据库或表来规定。将这些模型、数据库或表向存储装置33存储。利用该模型、数据库或表,统计性鉴定部61根据事前分布计算血管形态指标或血液流量指标的概率分布。统计性鉴定部61根据通过层次贝叶斯模型和马尔科夫链蒙特卡罗法得到的事后分布,统计性地鉴定潜在变量。
该鉴定问题具体而言,为不满足以下的3个条件的不良设定问题。3个条件是:(1)保证解的存在,(2)解唯一决定,(3)解相对于数据连续地变化,解相对于测量误差稳定。不良设定问题如果在由标准化理论和起扩展构成的框架内掌握,则容易处理。仅通过标准标准化理论解答不良设定问题是不充分的。为了解答不良设定问题,需要检测内部状态的不连续点且使检测出的不连续点对于内部状态的推测起作用的方法论。在这一点上,在有关本实施方式的统计性鉴定处理中,马尔科夫概率场理论是有效的。
在血管形态指标和血液流量指标中存在不确定性的环境下,统计性鉴定部61在适当的制约条件下鉴定潜在变量的概率分布参数。为了决定适当的制约条件,需要事前知道解的性质。统计性鉴定部61基于模拟和观测值产生关于解空间的制约条件的数据库。统计性鉴定部61利用产生的数据库,对超多自由度大规模问题执行基于马尔科夫概率场理论和层次贝叶斯模型的统计性鉴定处理。在作为制约条件的事前分布的设定中,基于许多的数值实验结果,将关于这些因素的参数的概率分布并列地单独构成。统计性鉴定部61通过将多个概率分布进行统一并对数据的缺损进行插补,来鉴定潜在变量的参数。为了该处理,统计性鉴定部61进行通过基于使用马尔科夫概率场理论的模型的层次贝叶斯法的推测。其方法是,根据作为解析对象的构造的变形状态的实测结果,基于鉴定出的中间变量,能够推测某个负荷条件和边界条件下的压力及流量分布。
冠状动脉的构造流体解析中的材料模型、边界条件及负荷条件的鉴定问题被定位为非线形逆解析,不能保证解的唯一性及稳定性的情况较多。生物体组织的材料特性及血压的现实中可取的范围能够作为先验信息设想,所以能够将它们设定为事前分布的概率分布。此外,压力及位移也可以设想在空间时间上是平滑的,所以能够将该信息也作为先验信息设定为事前分布的概率分布。或者,在能够考虑在血液的流动中不发生倒流的事实的情况下,也可以使用血管芯线方向的整体性的压力分布的斜率为负(存在压力下降)来作为制约条件。对于负荷条件(内压分布等)、边界条件及材料模型,能够将基于时间序列的CT图像的血管形状变形指标的观测值与基于力学模型的血管形状变形指标的预测值的方差(日语:2乗誤差)分布设定为数据分布。关于可观测的平均流量的方差分布也可以作为数据分布追加。能够基于这些事前分布和数据分布使用层次贝叶斯模型和蒙特卡罗法来计算事后分布。通过事后分布的发生概率及分散,能够得到潜在变量的参数的鉴定值。发生概率越高、分散越小,可以说是确信程度越高的鉴定值。在事后分布为多峰性分布的情况下,也只要选择多个鉴定值中的分散较小的鉴定值就可以。或者,在可能存在多个鉴定值的情况下,可以以各鉴定条件实施构造流体解析,识别各自的可能性,将鉴定值或解析结果作为诊断或预防的指针信息来利用。由于在时间序列的CT图像中也包含误差,所以在关于力学模型的各节点的血管形态指标中也包含误差。因此,也可以将各血管形态指标例如作为以从时间序列的CT图像计测出的血管形态指标的预测值作为平均值后的正规分布的概率变量来进行处理,在包含保持位置的空间的顺序这样的制约的基础上,设定事前分布。此外,在潜在变量的参数的鉴定中,存在不具有唯一性、可以考虑多个候选的情况。在此情况下,通过检查对于遵循根据时间序列的CT图像计测出的血管形态指标的不确定性的随机数的采样点的、潜在变量的鉴定值的样本集合的变动幅度,判定各鉴定值的候选的健壮性(稳定性)。也可以基于各鉴定值的候选的健壮性来决定最终的鉴定值。
接着,对力学模型的详细情况进行说明。如已述那样,力学模型构建部55可以根据力学模型的种类来构建不同类型的力学模型。在使用基于连续体力学的FEM的情况下,力学模型构建部55构建用于血管壁的应力解析用的形状模型(FEM模型)和用于血液的流体解析用的形状模型(FEM模型)这两者。
在使用基于材料力学的简单的鉴定法的情况下,根据材料力学中的受到内压的厚壁圆筒的式子近似地求出压力、弹性率和位移的关系。在此情况下,关于在芯线方向上排列的多个离散化区域分别使用厚壁圆筒近似来作为形状模型。具体而言,力学模型构建部55确定穿过在芯线上离散地排列的节点的截面上的血管内腔形状、血管壁表面形状和截面中心。接着,力学模型构建部55基于血管内腔形状和血管壁表面形状来计算平均面积、内腔的平均半径及平均壁厚。并且,力学模型构建部55基于平均面积、内腔的平均半径及平均壁厚,对各离散化区域的血管区域实施厚壁圆筒近似而构建形状模型。
在使用基于流动力学的简单的鉴定法的情况下,为了将血流的平均压力和平均流量近似地求出,使用流体力学中的修正伯努利式、或Hagen-Poiseuille Flow(哈根-泊肃叶流动)式。在此情况下,构建用来近似地求出多个离散化区域的各自的血压差与流量的关系的形状模型。这里,血压差是入口的压力与出口的压力的压力差,流量是指每单位时间的入口流量(或流速)和出口流量(或流速)。但是,力学模型构建部55可以基于在步骤S2中计算出的关于芯线方向的伸缩距离和扭转角,使关于各节点的移动矢量和关于与该节点相邻的节点的旋转位移在各时刻处相对应。
接着,对使用连续体力学模型的构造流体解析进行说明。
关于血管及血液(构成血管和血液的物质)的运动学与引起该运动的力是无关的。关于血管及血液的运动学的基本概念是将关于位置、时间、物体、运动及能够变形的物质的集合的直观的概念向数学用语进行了抽象化的概念。对关于血管及血液的变形及运动的局部的解析进行支配的基础的运动学张量是变形梯度张量F和速度梯度张量L。变形梯度张量F规定在运动的血管及血液的物质要素中发生的大小及形状的变化。变形梯度张量F用旋转张量(正规(日语:正格)正交张量)R、拉伸张量(正值对称张量)U、V的积表示。拉伸张量U、V通过首先施以关于由基本形态下的正规正交矢量R决定的方向的拉伸、接着施以由正规正交矢量R提供的刚体旋转而被带来。另外,施以拉伸的顺序和施以刚体旋转的顺序可以相反。速度梯度张量L不依存于基准形态而仅依存于当前形态。速度梯度张量L规定在运动的血管及血液的物质要素中发生的大小及形状的变化发生的速度。速度梯度张量L可以分离为应变速度张量D(对称张量)和自旋张量(反对称张量)。应变速度张量D表示物体刚好经过该当前形态时的拉伸的变化率。自旋张量表示物体刚好经过该当前形态时的旋转的变化率。对与血管的力学模型的外表面节点(或积分点)的一部分有关的变形梯度张量和速度梯度张量分配位移约束条件(包括时间性变化)。鉴定材料模型的参数、负荷条件(力学模型的内腔中的表面力矢量)及边界条件(血管边界中的力矢量),以使与力学模型的内腔节点(或积分点)有关的预测值(变形梯度张量、速度梯度张量或它们的函数值(例如也可以是位移或面积))和观测值(根据观测数据得到的变形梯度张量、速度梯度张量或它们的函数值)匹配。这里,血管内部的应力的初始状态既可以预先假定,也可以鉴定。
基于连续体力学的力学模型以表示运动中的血管及血液的质量、运动学、运动量、角运动量及能量的平衡的方程式为基础。质量、力、热及内部能量的概念是基本。所谓平衡法则,是指运动量的时间全微分等于物体力与接触力之和,此外角运动量的时间全微分等于物体转矩与接触转矩之和、此外运动能量和内部能量的时间性变化等于功率(力学能量)、每单位时间的热供给和热流束之和。根据平衡法则、结构式及跳跃条件,导出变形梯度张量、速度梯度张量及应力张量等场的方程式,能够记述血管及血液的力学模型。此外,血管及血液的应变场满足适合条件。这里,所谓结构式,是提供由密度、内部能量、速度矢量、应力张量、热流束矢量及温度构成的10个标量方程式组的关系的式子。通过场的平衡式中的17个标量场、即密度、内部能量、速度矢量、应力张量及热流束矢量中的场的平衡,提供8个标量关系,提供加上了温度的其余的未知量的关系。但是,物体力b和热源r为已知。提供这些标量场的方程式的参数是材料模型参数。
基于连续体力学的力学模型可以基于通过有限要素法或边界要素法的数值解析法,基于所提供的边界条件、负荷条件及材料模型来计算位移矢量、应力张量、应变张量及速度矢量等的场的近似解。
在构造-流体复合解析中,求解构造及流体的方程式的方法是一体型解法(monolithic method)和分离型解法(partitioned method)的哪种都可以。此外,构造与流体之间的边界面上的复合既可以是弱复合也可以是强复合。此外,在流体解析中,既可以用以ALE法为代表那样的边界面跟踪型的手法处理移动边界,也可以是Immersed Boundarymethod,Immersed Finite Element Method,Fictitious Domain Method等的边界面补充型的手法。
接着,作为简单的力学模型的例子,对受到内压和外压的厚壁圆筒的材料力学的式子和哈根-泊肃叶流动及修正伯努利的式子详细地说明。
首先,参照图24和图25,对厚壁圆筒的材料力学的式子进行说明。图24是表示壁厚圆筒的力学模型的正交截面的图,图25是图24的微小扇形要素的放大图。说明在内半径ra、外半径rb的厚壁圆筒上作用内压pa和外压pb的情况下的应力及应变、位移等的式子。E及ν是材料模型参数。E表示弹性率,ν表示泊松比。在厚壁圆筒中,需要也考虑半径应力σr,也考虑圆周应力σθ的半径方向的分布。以后,轴向的应变εz关于截面的位置及朝向是一样的。此外,对圆筒截面是轴对称的情况进行说明,但也可以是任意形状。在圆筒截面是轴对称的情况下,平衡条件只要仅对任意的截面的半径方向考虑就可以。对于在任意的截面上由半径r、r+dr的同心圆筒和中心角dθ切取的单位厚度1的微小扇形要素,考虑关于半径方向的力的平衡。由于变形也是轴对称的,所以在ab面及bc面上不发生剪切应力,所以仅作用垂直应力。因而,关于半径方向的力的平衡可以如以下的(5)式那样表现。
σr rdθ+2σθr sin(dθ/2)-(σr+(dσr/dr)dr)(r+dr)dθ=0…(5)
这里,由于dr比r小,dσr比σr小,所以如果将(5)式中包含的高次的微小项省略,设为sin(dθ/2)≒dθ/2,则(5)式可以如以下的(6)式那样表现。
rdσr/dr+σr-σθ=0…(6)
如果设半径r的半径方向的位移为u,则由于u+dr中的同方向的位移为u+(du/dr)dr,所以半径r的半径方向的应变εr为εr=du/dr。此外,根据半径方向的位移u,半径r的圆为半径r+u的圆。因而,圆周应变εθ可以如以下的(7)式那样表现。
εθ=(2π(r+u)-2πr)/2πr=u/r…(7)
此外,根据应力与应变的关系式,可得到以下的(8)式或(9)式。
d2u/dr2+(1/r)(du/dr)-u/r2=0…(8)
d2u/dr2+d(u/r)/dr=0…(9)
如果将(8)式或(9)式积分,则能够得到以下的(10)式。
u=c1r+c2/r…(10)
由此,可得到以下的(11)式、(12)式、(13)式。
σr=(E/((1+ν)(1-2ν)))(c1-(1-2ν)(c2/r2)+νεz)…(11)
σθ=(E/((1+ν)(1-2ν)))(c1+(1-2ν)(c2/r2)+νεz)…(12)
σz=(Eν/((1+ν)(1-2ν)))(2c1+((1-ν)/ν)εz)…(13)
(11)式、(12)式、(13)式中的常数c1、c2可以根据周边条件,即在圆筒的内周r=ra根据σr=-pa、在外周r=rb根据σr=-pb决定。通过该周边条件,从(11)式、(12)式及(13)式能够分别得到以下的(14)式、(15)式、(16)式。此外,位移u可以如以下的(17)式那样表现。
σr=(1/(rb 2-ra 2))(ra 2(1-rb 2/r2)pa-rb 2(1-ra 2/r2)pb)…(14)
σθ=(1/(rb 2-ra 2))(ra 2(1+rb 2/r2)pa-rb 2(1+ra 2/r2)pb)…(15)
σz=2ν(ra 2pa-rb 2pb)/(rb 2-ra 2)+Eεz=ν(σr+σθ)+Eεz…(16)
u=((1+ν)(1-2ν)/E)((ra 2pa-rb 2pb)/(rb 2-ra 2))r
+((1+ν)/E)((ra 2rb 2)/((rb 2-ra 2)r))(pa-pb)-νεzr…(17)
上述(16)式及(17)式包含εz的项。因而,(16)式根据作为对象的圆筒的边界处的约束条件而不同。例如在作为对象的圆筒的两端被约束的情况下为εz=0,两端被开放而为σz=0。
rθ面内的剪切应力τr可以如以下的(18)式那样表现。
τr=(1/2)(σr-σθ)=((ra 2rb 2)/((rb 2-ra 2)r2))pb…(18)
θz面内的剪切应力τr’当σz=0时为最大,所以如果σz=0,则可以如以下的(19)式那样表现。
τr’=(1/2)|σθ|=(1/2)((rb 2)/(rb 2-ra 2))(1+ra 2/r2)pb…(19)
接着,对哈根-泊肃叶流动及修正伯努利的式子进行说明。
如果流体向圆管内流入,则随着向下游行进而压力下降,并且流动的速度分布也逐渐变化。如果流动向管内流入,则边界层从管壁发展,随着向下游行进而边界层增加厚度,最终管内的流动被边界层覆盖。因此,速度分布从管入口的大致平坦的分布向下游的抛物线形分布变化,其以后的速度分布不再变化。将该状态称作完全发展的流动,由管摩擦损失带来的压力下降也为一定的比例。将流动从管入口达到发展后的流动的区间称作助跑区间或入口区间,将该区间的长度称作助跑距离或入口长度。由于完全发展后的流动的速度分布不向下游方向变化,所以如以下的(20)式所示,通过管摩擦损失发生的压力损失ΔP的作用力和通过流体的粘性发生的剪切应力τ的摩擦力平衡。在圆管内流动的情况下,当雷诺斯数Re为约2300以下时可以称为层流。
[数学式5]
速度分布u在层流的情况下,如以下的(21)式那样用轴对称的旋转抛物面表示。
[数学式6]
流量Q通过将速度分布u遍及管截面整体进行积分,能够如以下的(22)式那样表现。
[数学式7]
截面平均流速v可以如以下的(23)式那样表现。
[数学式8]
压力梯度向流动方向为一定,压力减小。压力梯度如果使用管长l的压力下降Δp,则可以如以下的(24)式那样表现。如果将(22)式向(24)式代入,则可得到(25)式。(25)式意味着流量Q与压力损失Δp成比例。将满足该关系的流动称作哈根-泊肃叶流动。
[数学式9]
有损失的情况下的修正伯努利的式子可以如以下的(26)式那样表现。
[数学式10]
Δh表示损失压头。损失压头是血管内壁的摩擦损失压头、狭窄部(变宽管或变细管的损失压头)或分叉部、弯曲管部等的损失压头、血管的解析对象区域入口/出口部的入口/出口损失压头、由起因于管内非恒常流动的压力变动的水击带来的损失压头的某种。例如,将管摩擦损失压头用以下的式子表示。
[数学式11]
这里,l是管的长度,d是管内径,v是管内平均速度,λ是管摩擦系数。管摩擦系数λ是在流动为层流的情况下由雷诺斯数Re决定、在紊流的情况下由雷诺斯数Re和表面粗糙度决定的值。根据流体的粘性,必定作用摩擦阻力。该摩擦阻力由于消耗驱动流动的动力或能量,所以成为能量损失。
如果将哈根-泊肃叶流动的式子如达西-威斯巴赫的式子那样变形,则可得到关于压力差和流速及管内径的以下的(28)式。
[数学式12]
以上,结束哈根-泊肃叶流动及修正伯努利的式子的说明。
如上述那样,图像处理装置27可以利用该材料力学的式子和哈根-泊肃叶流动及修正伯努利的式子来鉴定潜在变量。例如,考虑作为力学模型而将血管的变形使用厚壁管的材料力学的式子、将管径变化通过内压变化和弹性率表现的情况。在将无应力状态假定为初始形状(例如血管最收缩的状态)的情况下,如果将血管壁及斑块的弹性率设定为某个值,则能够得到血管内腔的平均半径等的血管形状变形指标的观测值的时间性变化量和内压的变化量的关系式。血管形状变形指标的观测值根据时间序列的CT图像来计测。决定血管的内压分布的时间性变化,以与该血管形状变形指标的观测值的时间性变化量一致。通过在该内压分布下进行血液的流体解析,计测血管流量指标的预测值。在该血管流量指标的预测值与观测值不一致的情况下,图像处理装置27将最初决定的血管壁或斑块的弹性率变更,再进行同样的解析。通过将其反复,图像处理装置27能够决定与血管形状变形指标的观测值和血液流量指标的观测值匹配的血管壁及斑块的弹性率、内压分布、流体解析的压力边界条件等的潜在变量。为了更有效率且稳定地进行该决定方法,也可以使用基于层次贝叶斯模型和马尔科夫链蒙特卡罗法的统计性鉴定方法。
如上述那样,有关本实施方式的血管解析装置50具有存储装置33、区域设定部51、图像解析-跟踪处理部53、力学模型构建部55及统计性鉴定部61。存储装置33将关于被检体的血管的时间序列的医用图像的数据存储。区域设定部51对时间序列的医用图像中包含的血管区域设定解析对象区域。图像解析-跟踪处理部53将时间序列的医用图像进行图像处理,计算解析对象区域的时间序列的形态指标和时间序列的形状变形指标。力学模型构建部55基于时间序列的形态指标、时间序列的形状变形指标和时间序列的医用图像,暂定地构建关于解析对象区域的构造流体解析的力学模型。统计性鉴定部61鉴定关于解析对象区域的力学模型的潜在变量,以使基于暂定地构建出的力学模型的血管形态指标及血液流量指标的预测值与预先计测出的血管形态指标及血液流量指标的观测值一致。
通过上述结构,有关本实施方式的血管解析装置50能够通过使用了血管形状变形指标和血液流量指标的逆解析,来鉴定材料模型、边界条件、负荷条件及几何学构造等的潜在变量。血管解析装置50通过一边将潜在变量变更一边反复进行逆解析,能够鉴定考虑了全部下述4点困难的潜在变量,所述4点困难为:1.冠状动脉的材料模型的鉴定方法、2.心脏的形状的变形对冠状动脉的影响的并入、3.冠状动脉的边界条件的鉴定方法、4.利用了具有不确定性的血管形状的材料模型、负荷条件、边界条件的鉴定方法。因而,血管解析装置50能够执行考虑了由在CT图像中未描绘出的血管或心脏等的外部因素带来的影响的构造流体解析。
这样,根据本实施方式,能够实现血管的构造流体解析的精度的提高。
说明了本发明的一些实施方式,但这些实施方式是作为例子提示的,并不是要限定发明的范围。这些新的实施方式能够以其他各种各样的形态实施,在不脱离发明的主旨的范围内能够进行各种省略、替换、变更。这些实施方式及其变形包含在发明的范围或主旨中,并且包含在权利要求书所记载的发明和其等价的范围中。
标号说明
10 CT架台;11X射线管;13X射线检测器;15数据收集装置;20控制台;21系统控制部;23架台控制部;25再构成装置;27图像处理装置;29输入设备;31显示设备;33存储装置;50血管解析装置;51区域设定部;53图像解析-跟踪处理部;55力学模型构建部;57血管应力解析部;59血液流体解析部;61统计性鉴定部。
Claims (15)
1.一种管状构造解析装置,具备:
设定部,对于与对象物的管状构造体有关的时间序列的图像设定解析对象区域,对于上述解析对象区域设定潜在变量鉴定区域;
计算部,计算上述解析对象区域的时间序列的形态指标和时间序列的形状变形指标;
构建部,基于上述时间序列的形态指标、上述时间序列的形状变形指标和上述时间序列的图像,构建与上述解析对象区域的构造流体解析有关的力学模型;以及
鉴定部,对与上述潜在变量鉴定区域有关的潜在变量进行鉴定,以使基于构建出的上述力学模型的形态指标的预测值及流量指标的预测值中的至少一方与预先计测出的形态指标的观测值及流量指标的观测值中的至少一方匹配。
2.如权利要求1所述的管状构造解析装置,其中,
上述鉴定部作为上述潜在变量而对与上述潜在变量鉴定区域有关的材料模型、与上述潜在变量鉴定区域的流体流入口及流体流出口有关的边界条件、包含上述潜在变量鉴定区域的由血流带来的内压分布在内的负荷条件、以及上述潜在变量鉴定区域的几何学构造的至少其一进行鉴定。
3.如权利要求2所述的管状构造解析装置,其中,
上述鉴定部具有:
应力解析部,对构建出的上述力学模型实施应力解析,推测上述形态指标的预测值;
流体解析部,对构建出的上述力学模型实施流体解析,推测上述流量指标的预测值;
第1统计性鉴定部,统计性地鉴定上述潜在变量,以使上述形态指标的预测值与上述形态指标的观测值匹配;以及
第2统计性鉴定部,统计性地鉴定上述潜在变量,以使上述流量指标的预测值与上述流量指标的观测值匹配。
4.如权利要求1所述的管状构造解析装置,其中,
还具备控制上述构建部和上述鉴定部的控制部;
上述控制部每当由上述鉴定部鉴定潜在变量就判定规定的鉴定结束指标是否达到既定值,在判定为上述鉴定结束指标没有达到上述既定值的情况下,控制上述构建部,利用鉴定出的上述潜在变量再构建力学模型,控制上述鉴定部,鉴定再构建出的上述力学模型的潜在变量,在判定为上述鉴定结束指标达到了上述既定值的情况下,控制上述构建部,将最新的力学模型设定为最终的力学模型。
5.如权利要求4所述的管状构造解析装置,其中,
还具备:
应力解析部,对上述最终的力学模型实施应力解析,计算时间序列的力学性指标;以及
显示部,显示上述时间序列的力学性指标。
6.如权利要求4所述的管状构造解析装置,其中,
还具备:
流体解析部,对上述最终的力学模型实施流体解析,计算时间序列的流体指标;以及
显示部,显示上述时间序列的流体指标。
7.如权利要求4所述的管状构造解析装置,其中,
还具备:
修正部,对上述时间序列的图像进行修正,以使基于上述最终的力学模型的力学性指标和流量指标中的至少一方的预测值与力学性指标和流量指标中的至少一方的观测值匹配。
8.如权利要求2所述的管状构造解析装置,其中,
上述构建部基于上述时间序列的形态指标和上述时间序列的图像来构建上述解析对象区域的形状模型,对上述形状模型分配上述时间序列的形状变形指标作为强制位移,通过分配上述潜在变量来构建上述力学模型。
9.如权利要求8所述的管状构造解析装置,其中,
上述构建部将上述力学模型划分为材料模型鉴定区域和非鉴定区域,对上述材料模型鉴定区域的表面部分配强制位移的约束条件,对上述材料模型鉴定区域的内部不分配强制位移的约束条件,对上述非鉴定区域的表面部及内部分配强制位移的约束条件。
10.如权利要求1所述的管状构造解析装置,其中,
还具备:
关联存储部,将基于上述力学模型的、由液压带来的内压的时间序列的空间分布、材料模型参数及形态指标建立关联而存储。
11.如权利要求1所述的管状构造解析装置,其中,
上述时间序列的形态指标和上述时间序列的形状变形指标以概率分布被提供。
12.如权利要求1所述的管状构造解析装置,其中,
上述管状构造体是被检体的血管,
上述形态指标是对应于上述血管的与血管区域的形态有关的血管形态指标,
上述流量指标是与流过上述血管区域的血液的流量有关的血液流量指标,
上述构建部将构建出的上述力学模型按照来自超声波诊断装置的斑块指标划分为多个部分区域,对上述多个部分区域的每个分配用于材料模型的参数。
13.如权利要求1所述的管状构造解析装置,其中,
上述管状构造体是被检体的血管,
上述形态指标是对应于上述血管的与血管区域的形态有关的血管形态指标,
上述流量指标是与流过上述血管区域的血液的流量有关的血液流量指标,
上述血液流量指标的观测值是基于上述时间序列的图像的血管内腔的容积的时间变化、或基于超声波多普勒计测的血流量的时间变化。
14.如权利要求1~11中任一项所述的管状构造解析装置,其中,
上述管状构造体是被检体的血管,
上述形态指标是对应于上述血管的与血管区域的形态有关的血管形态指标,
上述流量指标是与流过上述血管区域的血液的流量有关的血液流量指标。
15.一种管状构造解析方法,具备如下步骤:
对于与对象物的管状构造体有关的时间序列的图像设定解析对象区域,对于上述解析对象区域设定潜在变量鉴定区域;
计算上述解析对象区域的时间序列的形态指标和时间序列的形状变形指标;
基于上述时间序列的形态指标、上述时间序列的形状变形指标和上述时间序列的图像,构建与上述解析对象区域的构造流体解析有关的力学模型;以及
对与上述潜在变量鉴定区域有关的潜在变量进行鉴定,以使基于构建出的上述力学模型的形态指标的预测值及流量指标的预测值中的至少一方与预先计测出的形态指标的观测值及流量指标的观测值中的至少一方匹配。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012253619A JP2014100249A (ja) | 2012-11-19 | 2012-11-19 | 血管解析装置、医用画像診断装置、血管解析方法、及び血管解析プログラム |
JP2012-253619 | 2012-11-19 | ||
CN201380054742.6A CN104736061A (zh) | 2012-11-19 | 2013-11-12 | 血管解析装置、医用图像诊断装置及血管解析方法 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380054742.6A Division CN104736061A (zh) | 2012-11-19 | 2013-11-12 | 血管解析装置、医用图像诊断装置及血管解析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107582041A true CN107582041A (zh) | 2018-01-16 |
CN107582041B CN107582041B (zh) | 2021-07-06 |
Family
ID=50731144
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380054742.6A Pending CN104736061A (zh) | 2012-11-19 | 2013-11-12 | 血管解析装置、医用图像诊断装置及血管解析方法 |
CN201710695287.8A Active CN107582041B (zh) | 2012-11-19 | 2013-11-12 | 血管解析装置及血管解析方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380054742.6A Pending CN104736061A (zh) | 2012-11-19 | 2013-11-12 | 血管解析装置、医用图像诊断装置及血管解析方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US11357409B2 (zh) |
JP (1) | JP2014100249A (zh) |
CN (2) | CN104736061A (zh) |
WO (1) | WO2014077233A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110559015B (zh) * | 2019-08-26 | 2020-12-22 | 清华大学 | 血管生理参数测量方法、设备、计算机设备和存储介质 |
CN113876304A (zh) * | 2021-09-08 | 2022-01-04 | 深圳市中科微光医疗器械技术有限公司 | 一种基于oct图像和造影图像确定ffr的方法及装置 |
Families Citing this family (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9808222B2 (en) | 2009-10-12 | 2017-11-07 | Acist Medical Systems, Inc. | Intravascular ultrasound system for co-registered imaging |
JP2014100249A (ja) * | 2012-11-19 | 2014-06-05 | Toshiba Corp | 血管解析装置、医用画像診断装置、血管解析方法、及び血管解析プログラム |
JP6091870B2 (ja) | 2012-12-07 | 2017-03-08 | 東芝メディカルシステムズ株式会社 | 血管解析装置、医用画像診断装置、血管解析方法、及び血管解析プログラム |
EP2784748B1 (en) * | 2013-03-28 | 2017-11-01 | Expert Ymaging, SL | A computer implemented method for assessing vascular networks from medical images and uses thereof |
US9986938B2 (en) * | 2014-02-25 | 2018-06-05 | Medis Associated B.V. | Method and device for determining a geometrical parameter of a blood vessel |
NL2012459B1 (en) | 2014-03-18 | 2016-01-08 | Medis Ass B V | Method and device for determining deviation in pressure in a blood vessel. |
US9058692B1 (en) * | 2014-04-16 | 2015-06-16 | Heartflow, Inc. | Systems and methods for image-based object modeling using multiple image acquisitions or reconstructions |
US9754082B2 (en) | 2014-05-30 | 2017-09-05 | Heartflow, Inc. | Systems and methods for reporting blood flow characteristics |
WO2016081517A2 (en) * | 2014-11-17 | 2016-05-26 | Borkholder David A | Pulse wave velocity, arterial compliance, and blood pressure |
WO2016182508A1 (en) * | 2015-05-12 | 2016-11-17 | Singapore Health Services Pte Ltd | Medical image processing methods and systems |
JP2018519047A (ja) | 2015-06-19 | 2018-07-19 | ニューラル アナリティクス、インコーポレイテッド | 頭蓋内ドップラープローブ |
US10653393B2 (en) | 2015-10-08 | 2020-05-19 | Acist Medical Systems, Inc. | Intravascular ultrasound imaging with frequency selective imaging methods and systems |
US10909661B2 (en) | 2015-10-08 | 2021-02-02 | Acist Medical Systems, Inc. | Systems and methods to reduce near-field artifacts |
US11369337B2 (en) * | 2015-12-11 | 2022-06-28 | Acist Medical Systems, Inc. | Detection of disturbed blood flow |
CN108778140A (zh) | 2016-01-05 | 2018-11-09 | 神经系统分析公司 | 用于确定临床指征的系统和方法 |
CN108778141A (zh) | 2016-01-05 | 2018-11-09 | 神经系统分析公司 | 集成化探针结构 |
US11589836B2 (en) | 2016-01-05 | 2023-02-28 | Novasignal Corp. | Systems and methods for detecting neurological conditions |
US10278662B2 (en) | 2016-02-05 | 2019-05-07 | Toshiba Medical Systems Corporation | Image processing apparatus and medical image diagnostic apparatus |
JP2017176202A (ja) * | 2016-03-28 | 2017-10-05 | コニカミノルタ株式会社 | 動態解析システム |
JP6378715B2 (ja) * | 2016-04-21 | 2018-08-22 | ゼネラル・エレクトリック・カンパニイ | 血管検出装置、磁気共鳴イメージング装置、およびプログラム |
JP6929689B2 (ja) * | 2016-04-26 | 2021-09-01 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置及び医用画像診断装置 |
US9947093B2 (en) | 2016-05-03 | 2018-04-17 | Konica Minolta, Inc. | Dynamic analysis apparatus and dynamic analysis system |
US10489919B2 (en) * | 2016-05-16 | 2019-11-26 | Acist Medical Systems, Inc. | Motion-based image segmentation systems and methods |
US10610184B2 (en) | 2016-08-12 | 2020-04-07 | Canon Medical Systems Corporation | Medical-information processing apparatus and X-ray CT apparatus |
DE102017217599A1 (de) * | 2016-10-04 | 2018-04-05 | Toshiba Medical Systems Corporation | Medizinische Informationsverarbeitungsvorrichtung, Röntgen-CT-Vorrichtung und medizinisches Informationsverarbeitungsverfahren |
JP2018068631A (ja) * | 2016-10-28 | 2018-05-10 | キヤノン株式会社 | 放射線撮影システム、放射線表示方法 |
US10163209B2 (en) * | 2016-11-23 | 2018-12-25 | Toshiba Medical Systems Corporation | Medical image processing apparatus, medical image processing method, and X-ray CT apparatus |
DE102017221276A1 (de) * | 2016-11-28 | 2018-05-30 | Toshiba Medical Systems Corporation | Medizinische Bildverarbeitungsvorrichtung, Röntgen-CT Vorrichtung und medizinisches Bildverarbeitungsverfahren |
CN106780527B (zh) * | 2016-11-29 | 2020-09-15 | 上海联影医疗科技有限公司 | 医学图像中血管进出口、边界条件获取方法及处理装置 |
EP3554375A1 (en) * | 2016-12-15 | 2019-10-23 | Koninklijke Philips N.V. | Visualizing collimation errors |
WO2018133118A1 (zh) * | 2017-01-23 | 2018-07-26 | 上海联影医疗科技有限公司 | 血流状态分析系统及方法 |
EP3382583A1 (en) * | 2017-03-29 | 2018-10-03 | Koninklijke Philips N.V. | Hemodynamic simulation of movement inducted vascular deformations |
JP6847453B2 (ja) * | 2017-05-30 | 2021-03-24 | 国立大学法人九州大学 | 地図情報生成方法、作動方法、及びプログラム |
JP7073661B2 (ja) * | 2017-09-27 | 2022-05-24 | コニカミノルタ株式会社 | 動態解析装置及び動態解析システム |
CN111356406A (zh) * | 2017-10-06 | 2020-06-30 | 埃默里大学 | 用于确定一个或多个动脉节段的血液动力学信息的方法和系统 |
JP7170459B2 (ja) * | 2018-08-10 | 2022-11-14 | 株式会社東芝 | 血行検出装置、方法、及びプログラム |
CN112022137B (zh) | 2018-11-30 | 2021-07-13 | 博动医学影像科技(上海)有限公司 | 建立血管截面函数和血管应力的方法及装置 |
EP3685756A1 (en) * | 2019-01-24 | 2020-07-29 | Koninklijke Philips N.V. | Methods and systems for investigating blood vessel characteristics |
JP2019171102A (ja) * | 2019-05-30 | 2019-10-10 | コニカミノルタ株式会社 | 画像処理装置、画像処理方法、およびプログラム |
US11024034B2 (en) | 2019-07-02 | 2021-06-01 | Acist Medical Systems, Inc. | Image segmentation confidence determination |
CN112151180B (zh) | 2019-12-05 | 2024-03-08 | 苏州润迈德医疗科技有限公司 | 具有狭窄病变的血管数学模型的合成方法和装置 |
US20220110530A2 (en) * | 2019-12-09 | 2022-04-14 | Nordsletten David | Method and System for Estimating Pressure Difference in Turbulent Flow |
JP2023061477A (ja) * | 2021-10-20 | 2023-05-02 | 富士通株式会社 | プログラム、データ処理方法及びデータ処理装置 |
CN116705330B (zh) * | 2023-07-31 | 2023-11-10 | 柏意慧心(杭州)网络科技有限公司 | 确定血管壁的弹性特征的方法、计算设备和介质 |
CN117316359B (zh) * | 2023-09-22 | 2024-04-09 | 杭州威灿科技有限公司 | 血液检测过程跟踪方法、装置、设备及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1647067A (zh) * | 2002-04-04 | 2005-07-27 | 石原产业株式会社 | 数据解析装置及方法 |
JP2008514368A (ja) * | 2004-09-30 | 2008-05-08 | ザ リージェント オブ ザ ユニバーシティ オブ カリフォルニア | 特性判定方法、記録装置、評価方法 |
CN101542530A (zh) * | 2006-11-22 | 2009-09-23 | 皇家飞利浦电子股份有限公司 | 基于有限数据集合的图像生成 |
CN101680018A (zh) * | 2007-01-10 | 2010-03-24 | 海莫希尔有限责任公司 | 体外血液动力学的内皮/平滑肌细胞共培养模型在鉴定血管疾病的新型治疗靶标中的应用 |
US20120041739A1 (en) * | 2010-08-12 | 2012-02-16 | Heartflow, Inc. | Method and System for Patient-Specific Modeling of Blood Flow |
WO2012028190A1 (en) * | 2010-09-02 | 2012-03-08 | Pie Medical Imaging Bv | Method and apparatus for quantitative analysis of a tree of recursively splitting tubular organs |
CN104736061A (zh) * | 2012-11-19 | 2015-06-24 | 株式会社东芝 | 血管解析装置、医用图像诊断装置及血管解析方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7117026B2 (en) * | 2002-06-12 | 2006-10-03 | Koninklijke Philips Electronics N.V. | Physiological model based non-rigid image registration |
US20060069318A1 (en) | 2004-09-30 | 2006-03-30 | The Regents Of The University Of California | Method for assessment of the structure-function characteristics of structures in a human or animal body |
EP1960965A2 (en) | 2005-12-09 | 2008-08-27 | Koninklijke Philips Electronics N.V. | Model-based flow analysis and visualization |
JP4439533B2 (ja) | 2007-03-27 | 2010-03-24 | 株式会社東芝 | 負荷算定装置および負荷算定方法 |
US8705819B2 (en) | 2007-06-21 | 2014-04-22 | Koninklijke Philips N.V. | Adjusting acquisition protocols for dynamic medical imaging using dynamic models |
US8675940B2 (en) | 2009-10-27 | 2014-03-18 | Siemens Aktiengesellschaft | Generation of moving vascular models and blood flow analysis from moving vascular models and phase contrast MRI |
US8157742B2 (en) | 2010-08-12 | 2012-04-17 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
JP6091870B2 (ja) * | 2012-12-07 | 2017-03-08 | 東芝メディカルシステムズ株式会社 | 血管解析装置、医用画像診断装置、血管解析方法、及び血管解析プログラム |
-
2012
- 2012-11-19 JP JP2012253619A patent/JP2014100249A/ja not_active Abandoned
-
2013
- 2013-11-12 WO PCT/JP2013/080502 patent/WO2014077233A1/ja active Application Filing
- 2013-11-12 CN CN201380054742.6A patent/CN104736061A/zh active Pending
- 2013-11-12 CN CN201710695287.8A patent/CN107582041B/zh active Active
-
2015
- 2015-05-19 US US14/716,433 patent/US11357409B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1647067A (zh) * | 2002-04-04 | 2005-07-27 | 石原产业株式会社 | 数据解析装置及方法 |
JP2008514368A (ja) * | 2004-09-30 | 2008-05-08 | ザ リージェント オブ ザ ユニバーシティ オブ カリフォルニア | 特性判定方法、記録装置、評価方法 |
CN101542530A (zh) * | 2006-11-22 | 2009-09-23 | 皇家飞利浦电子股份有限公司 | 基于有限数据集合的图像生成 |
CN101680018A (zh) * | 2007-01-10 | 2010-03-24 | 海莫希尔有限责任公司 | 体外血液动力学的内皮/平滑肌细胞共培养模型在鉴定血管疾病的新型治疗靶标中的应用 |
US20120041739A1 (en) * | 2010-08-12 | 2012-02-16 | Heartflow, Inc. | Method and System for Patient-Specific Modeling of Blood Flow |
US20120041318A1 (en) * | 2010-08-12 | 2012-02-16 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
WO2012028190A1 (en) * | 2010-09-02 | 2012-03-08 | Pie Medical Imaging Bv | Method and apparatus for quantitative analysis of a tree of recursively splitting tubular organs |
CN104736061A (zh) * | 2012-11-19 | 2015-06-24 | 株式会社东芝 | 血管解析装置、医用图像诊断装置及血管解析方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110559015B (zh) * | 2019-08-26 | 2020-12-22 | 清华大学 | 血管生理参数测量方法、设备、计算机设备和存储介质 |
CN113876304A (zh) * | 2021-09-08 | 2022-01-04 | 深圳市中科微光医疗器械技术有限公司 | 一种基于oct图像和造影图像确定ffr的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
US11357409B2 (en) | 2022-06-14 |
JP2014100249A (ja) | 2014-06-05 |
CN104736061A (zh) | 2015-06-24 |
WO2014077233A1 (ja) | 2014-05-22 |
CN107582041B (zh) | 2021-07-06 |
US20150245776A1 (en) | 2015-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107582041A (zh) | 管状构造解析装置及管状构造解析方法 | |
US11854704B2 (en) | Systems and methods for anatomical modeling using information obtained from a medical procedure | |
US11887305B1 (en) | Systems and methods for identifying personalized vascular implants from patient-specific anatomic data | |
US10354759B2 (en) | Systems and methods for modeling changes in patient-specific blood vessel geometry and boundary conditions | |
CN104837407B (zh) | 血管解析装置、医用图像诊断装置以及血管解析方法 | |
JP7475376B2 (ja) | Mriのストレイン測定値による心機能の迅速な定量評価 | |
CN104736046B (zh) | 用于数字评估脉管系统的系统和方法 | |
KR101879560B1 (ko) | 환자별 혈류 모델링 방법 및 시스템 | |
CN109036551B (zh) | 一种冠状动脉生理学指标关系建立及应用方法、装置 | |
JP6362853B2 (ja) | 血管解析装置、および血管解析装置の作動方法 | |
CN106659400A (zh) | 用于确定血流储备分数值的装置 | |
WO2013180851A1 (en) | A method for estimating pressure gradients and fractional flow reserve from computed tomography angiography: transluminal attenuation flow encoding | |
JP6362851B2 (ja) | 血管解析装置、血管解析プログラム、及び血管解析装置の作動方法 | |
CN108109698A (zh) | 应用计算流体力学来模拟计算血流储备分数的系统和方法 | |
JP6275797B2 (ja) | 管状構造解析装置、管状構造解析方法及び管状構造解析プログラム |
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 |