CN101732036A - 医用图像解析装置以及图像解析用控制程序 - Google Patents

医用图像解析装置以及图像解析用控制程序 Download PDF

Info

Publication number
CN101732036A
CN101732036A CN200910208856A CN200910208856A CN101732036A CN 101732036 A CN101732036 A CN 101732036A CN 200910208856 A CN200910208856 A CN 200910208856A CN 200910208856 A CN200910208856 A CN 200910208856A CN 101732036 A CN101732036 A CN 101732036A
Authority
CN
China
Prior art keywords
mentioned
view data
parameter
data
image analysis
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
Application number
CN200910208856A
Other languages
English (en)
Other versions
CN101732036B (zh
Inventor
川岸哲也
大内启之
阿部康彦
T·海莱-瓦利
J·A·C·利马
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Johns Hopkins University
Original Assignee
Toshiba Medical Systems Corp
Johns Hopkins University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Medical Systems Corp, Johns Hopkins University filed Critical Toshiba Medical Systems Corp
Publication of CN101732036A publication Critical patent/CN101732036A/zh
Application granted granted Critical
Publication of CN101732036B publication Critical patent/CN101732036B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/246Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • G06T7/248Analysis of motion using feature-based methods, e.g. the tracking of corners or segments involving reference images or patches
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Quality & Reliability (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Multimedia (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明提供一种医用图像解析装置以及图像解析用控制程序,根据从不同的医用图像诊断装置中收集到的图像数据有效评价生物体组织的运动功能。通过对另外设置的医用图像诊断装置供给的被检体的时间序列的图像数据适用共用解析算法,设定部对从上述图像数据中提取的基准图像数据的心肌组织设定多个关心点,跟踪处理部根据通过上述基准图像数据与其后续图像数据之间的跟踪处理取得的上述关心点的心肌组织的位移信息,测量出运动参数。另外,数据生成部根据这些测量结果,作为参数数据生成表示运动参数的2维分布的参数图像数据或表示运动参数的时间变化的参数时间序列数据并显示该参数数据。

Description

医用图像解析装置以及图像解析用控制程序
技术领域
本发明涉及医用图像解析装置以及图像解析用控制程序,特别涉及根据通过解析从被检体中收集到的时间序列的图像数据得到的生物体组织的运动参数进行参数图像数据等的生成的医用图像解析装置以及图像解析用控制程序。
背景技术
医用图像诊断近几年由于随着计算机技术的发展而被实用化的X射线CT装置、MRI装置、超声波诊断装置等的发展而取得了快速的进步,在现在的医疗中已成为必不可少的技术。特别是,在X射线CT装置或MRI装置中,随着生物体信息检测装置或运算处理装置的高速化以及高性能化能够实时显示图像数据。另外,超声波诊断装置只需通过将超声波探头接触体表的简单操作就可以容易地观察实时的2维图像数据。因此,X射线CT装置、MRI装置和超声波诊断装置被广泛应用于生物体内脏器官的形态诊断或功能诊断。
特别是,在近几年的超声波诊断领域中,开发了根据通过解析按照时间序列收集到的B模式图像数据等超声波图像数据得到的心肌组织的位移信息,例如能够进行“应变(strain)”的2维观测的应变成像法。
在进行心脏的功能诊断的应变成像法中,首先,根据通过对被检体的超声波扫描得到的接收信号,按照时间序列收集B模式图像数据,其次,对在时间方向上邻接的超声波图像数据适用基于模式匹配的跟踪处理,测量出心肌组织的各部位的“位移”。然后,通过算出作为每单位长度的位移量被定义的“应变”的2维分布来进行应变图像数据的生成。
另外,也提出了根据通过应用彩色多普勒法2维显示心肌组织的移动速度的TDI(Tissue Doppler Imaging:组织多普勒成像)法取得的上述移动速度的空间倾斜,测量出“应变速度”的2维分布,并通过对该“应变速度”进行时间积分生成应变图像数据的方法(例如,参照专利文献:日本特开2005-130877号公报)。
另一方面,也探讨了使用能够在短时间内收集时间序列的图像数据的X射线CT装置或MRI装置的应变成像法,特别是,因为肋骨或肺等阻碍超声波传播而难以使用超声波诊断装置得到四腔剖面或两腔剖面的应变图像数据,而根据使用X射线CT装置或MRI装置的应变成像法可以容易地生成四腔剖面或两腔剖面的应变图像数据。
即,通过选择与诊断目的或诊断部位对应的应变成像法中使用的医用图像诊断装置,能够生成对诊断有效的应变图像数据,同时通过比较观察根据从不同的多个医用图像诊断装置取得的图像数据生成的各种应变图像数据,可以取得更高诊断精度。
但是,上述那样解析从不同的多个医用图像诊断装置收集到的图像数据来生成应变图像数据时,以往的解析装置或解析算法根据医用图像诊断装置的类型而不同,因此难以有效进行应变图像数据的生成,再者,不能高精度地比较观察各种应变图像数据。
发明内容
本发明鉴于上述问题点而作出,其目的在于:提供一种使用对由不同的医用图像诊断装置的各个生成的生物体组织的时间序列的图像数据适用共用的解析算法而测量出的“应变”等运动参数,能够有效评价生物体组织的运动功能的医用图像解析装置以及图像解析用控制程序。
为了实现上述目的,根据本发明的一个方面,提供一种医用图像解析装置,其特征在于,包括:图像数据收集部,收集由不同的医用图像诊断装置的各个生成的时间序列的图像数据;运动参数测量部,通过对上述时间序列的图像数据进行处理,测量出生物体组织的运动参数;参数数据生成部,根据上述运动参数生成参数数据;和显示部,显示上述参数数据;其中,上述运动参数测量部通过使用共用的解析算法解析从上述医用图像诊断装置的各个收集到的上述图像数据,测量出上述运动参数。
根据本发明的另一个方面,提供一种图像解析用控制程序,其特征在于,使医用图像解析装置执行以下功能:执行收集由不同的医用图像诊断装置的各个生成的时间序列的图像数据的功能;通过对上述时间序列的图像数据进行处理来测量出生物体组织的运动参数的功能;根据上述运动参数生成参数数据的功能;和显示上述参数数据的功能。
根据本发明,通过对由不同的医用图像诊断装置的各个生成的生物体组织的时间序列的图像数据适用预先设定的共用的解析算法,测量出希望的运动参数,可以有效评价生物体组织的运动功能。因此可以提高诊断效率与诊断精度,减轻对操作者的负担。
附图说明
图1为本发明的实施例中的医用图像解析装置的功能框图。
图2为用于说明在该实施例中由医用图像诊断装置收集到的1个心跳周期的时间序列的图像数据的图。
图3A和3B为用于说明该实施例的运动参数测量部进行的跟踪处理的具体例子的图。
图4为表示该实施例的参数数据生成部生成的参数图像数据的具体例子的图。
图5为表示该实施例的参数数据生成部生成的参数时间序列数据的具体例子的图。
图6为表示该实施例的显示部上显示的显示数据的具体例子的图。
图7为该实施例中的参数数据的生成/显示顺序的流程图。
具体实施方式
以下,参照附图说明本发明的实施例。
以下所述的本发明的实施例中的医用图像解析装置,收集另外设置的医用图像诊断装置经由网络等供给的被检体的规定心跳周期的时间序列的图像数据。然后,根据预先设定的共用解析算法(即,对利用不同的多个医用图像诊断装置取得的图像数据都能够适用的解析算法),对从时间序列的多个图像数据中提取的基准图像数据的心肌组织以规定间隔设定多个关心点。然后,根据通过使用上述基准图像数据与该基准图像数据的后续图像数据的跟踪处理取得的上述关心点的心肌组织的位移量或位移方向等,测量出“应变”等运动参数。然后,根据这些测量结果,作为参数数据生成表示运动参数的2维分布的参数图像数据或表示运动参数的时间变化的参数时间序列数据并在显示部上加以显示。
另外,在本实施例中,针对将心肌组织的“应变”作为运动参数进行测量的情况进行了叙述,但是并不只限于此,例如,也可以将心肌组织的“位移”、“旋转”、“扭转(torsion)”、“速度”等作为运动参数进行测量。也可以将表示时间变化的“应变率”、“旋转率”、“扭转率”、“加速度”等作为运动参数,并且,在心肌组织以外的生物体组织中也可以测量上述运动参数。
并且,在本实施例中,针对从医用图像诊断装置预先供给的规定心跳周期的时间序列的图像数据的各个生成参数图像数据,并将取得的规定心跳周期的参数图像数据与上述图像数据重叠并重复显示(循环显示)的情况进行了叙述,但是运动参数的测量中的处理速度比图像数据的收集中的处理速度快时,可以从医用图像诊断装置依次供给的图像数据的各个生成参数图像数据,也可以将取得的参数图像数据与上述图像数据重叠并实时显示。
(装置的构成)
使用图1至图6对本发明的第1实施例中的医用图像解析装置的构成与基本动作进行说明。另外,图1为医用图像解析装置的功能框图。
图1所示的医用图像解析装置100具备:图像数据收集部1,经由网络10收集各种医用图像诊断装置或图像数据服务器等供给的规定心跳周期的心脏的时间序列的图像数据;运动参数测量部2,根据该时间序列的图像数据测量出被检体的心肌组织的运动参数;参数数据生成部3,根据测量出的运动参数生成各种参数数据;显示部4,显示取得的参数数据;输入部5,进行收集图像数据的医用图像诊断装置(医用图像摄影装置:modality)的选择、运动参数的选择、对于后述基准图像数据的心内膜以及心外膜的设定、各种指令信号的输入等;系统控制部6,整体上控制医用图像解析装置100具备的上述各单元。
图像数据收集部1具备未图示的存储电路,另外设置的超声波诊断装置、X射线CT装置、MRI装置等医用图像诊断装置或图像数据服务器经由网络10供给的规定心跳周期的心脏的时间序列的图像数据,与心跳时相信息等附带信息一起保存到上述存储电路中。此时,将在检测出心电波形的R波的时刻生成的图像数据与R波定时(timing)信息一起保存到上述存储电路中。
运动参数测量部2,具备基准图像数据提取部21、关心点设定部22、跟踪处理部23以及跟踪条件保管部24。基准图像数据提取部21读出图像数据收集部1的存储电路中保存的时间序列的图像数据,从这些图像数据中,例如作为基准图像数据提取具有R波定时信息的图像数据并在显示部4上加以显示。
关心点设定部22,对由输入部5对显示部4上显示的上述基准图像数据设定的心内膜以及心外膜包围的心肌组织以规定间隔设定多个关心点。
跟踪处理部23形成将对基准图像数据的心肌组织以规定间隔设定的关心点的各个作为中心的规定尺寸的模板。并且,按照该模板中所包含的多个像素与上述基准图像数据的后续图像数据的像素之间的模式匹配(模板匹配)进行跟踪处理,测量出设定上述关心点的心肌组织的局部位移量以及局部位移方向。
例如,图2所示,跟踪处理部23从图像数据收集部1的存储电路中保存的1个心跳周期T0的时间序列的图像数据A1至AM中提取具有心电波形Ec的R波定时信息的时刻T1的基准图像数据A1,进行对该基准图像数据A1形成的模板与基准图像数据A1的时刻T2(T2=T1+δT)的后续图像数据A2之间的模式匹配,测量出以基准图像数据A1的关心点为中心的心肌组织的图像数据A2的局部位移量以及局部位移方向,进而根据该位移量与位移方向,作为运动参数测量出心肌组织的“应变”。
具体来说,跟踪处理部23对基准图像数据A1形成将对基准图像数据A1的心肌组织以规定间隔设定的关心点的各个作为中心的2维模板,相对于基准图像数据A1的后续图像数据A2使该模板相对移动,从而进行上述模板中包含的多个像素与图像数据A2的像素之间的互相关运算。并且,通过检测出取得最大相关值的模板相对于图像数据A2的相对移动方向以及相对移动距离,测量出以关心点为中心的心肌组织的局部位移量以及局部位移方向,进而作为“应变”测量出每单位长度的位移量。
通过对基准图像数据A1形成的模板与图像数据A2之间的模式匹配,如果测量出设定了上述关心点的心肌组织的局部“应变”,则按照与上述同样的方式,通过使用了图像数据A2与图像数据A3、图像数据A3与图像数据A4.....这样邻接的2个图像数据的模式匹配,依次测量出时刻T3(T3=T1+2δT)、时刻T4(T4=T1+3δT)、.....的各个局部“应变”。
另外,通过先行的图像数据与后续的图像数据(例如,图像数据A2与图像数据A3)之间的模式匹配测量心肌组织的“应变”时,图像数据A2的模板也可以根据关心点设定部22对图像数据A1设定的关心点的位置信息形成,但是也可以根据通过图像数据A1与图像数据A2之间的模式匹配测量出的局部心肌组织的移动信息形成。
其次,针对使用互相关运算的图像数据的跟踪处理,使用图3更详细地进行说明。图3A所示的关心点Cg是对图像数据A1的心肌组织以规定间隔设定的多个关心点之一,如果以该关心点Cg为中心将具有规定尺寸(即,规定像素数N0(N0=Px·Qy))的模板Tg的像素值作为f1(px,qy)、将图像数据A2的像素值作为f2(px,qy),则通过根据下面公式(1)算出互相关系数γ12(k,s),可以测量出设定了关心点Cg的心肌组织的时间δT后的位移量与位移方向。
γ 12 ( k , s ) = 1 No σ 1 σ 2 Σ Px = 1 Px Σ qy = 1 Qy ( f 1 ( px , qy ) - f ‾ 1 ) ( f 2 ( px + k , qy + s ) - f 2 ‾ )
f 1 = ‾ 1 No Σ Px = 1 Px Σ qy = 1 Qy f 1 ( px , qy ) f 2 = ‾ 1 No Σ Px = 1 Px Σ qy = 1 Qy f 2 ( px + k , qy + s )
σ 1 2 = 1 No Σ Px = 1 Px Σ qy = 1 Qy ( f 1 ( px , qy ) - f 1 ‾ ) 2 σ 2 2 = 1 No Σ Px = 1 Px Σ qy = 1 Qy ( f 2 ( px , qy ) - f 2 ‾ ) 2
No=PxQy                                   …(1)
其中,Px以及Qy是模板Tg的px方向以及qy方向上的像素数,通常,对图像数据A1的心肌组织设定的关心点Cg位于模板Tg的大致中央位置。该互相关运算的结果,当k=k1(参照图3B)、s=s1(未图示)、γ12(k,s)具有最大值时,表示图像数据A1的具有关心点Cg的局部心肌组织在图像数据A2中向px方向移位k1个像素,向qy方向移位s1个像素。
对于图像数据A1的心肌组织设定的全部关心点进行这样的跟踪处理,从而测量出设定了这些关心点的图像数据A1的心肌组织的时间δT后的局部位移量与局部位移方向。然后,基于2维测量出的位移量以及位移方向,测量出被定义为每单位长度的位移量的“应变”作为运动参数。
并且,以同样的顺序,对图像数据A2的后续图像数据A3至AM的各个也进行同样的跟踪处理,测量出这些图像数据的心肌组织的局部“应变”。
返回到图1,在运动参数测量部2的跟踪条件保管部24中,针对各医用图像摄影装置保存关心点设定部22对基准图像数据A1的心肌组织设定的关心点之间的间隔或跟踪处理部23在图像数据之间的模式匹配中使用的模板尺寸、观察区域、动态范围、增益、平滑处理等跟踪条件。
其次,参数数据生成部3具备图像数据生成部31与时间序列数据生成部32。
图像数据生成部31通过排列运动参数测量部2的跟踪处理部23从设定了多个关心点的图像数据的心肌组织测量出的运动参数,使其与关心点的位置信息对应,从而生成参数图像数据。并且,根据上述方法,依次生成与图像数据收集部1供给的规定心跳周期的时间序列的图像数据的各个对应的参数图像数据,向显示部4供给作为附带信息附加在上述图像数据中的心跳时相信息。
图4为图像数据生成部31生成的参数图像数据的具体例子,该参数图像数据是在根据医用图像诊断装置生成包含被检体的左心房LA以及左心室LV的剖面的图像数据时取得的参数图像数据。另外,直接使用通过对基准图像数据的关心点的心肌组织进行跟踪处理测量出的图像数据的关心点的运动参数,也可以生成参数图像数据。但是,如图4所示,通过使用在沿着心肌组织预先设定的测量区域R1至R11的各个中所包含的多个关心点上测量出的运动参数的平均值(平均运动参数),可以生成对诊断有效的参数图像数据。
即,在图4中,例如对左心房LA的心肌组织设定测量区域R1至R4、对左心室LV的心肌组织设定测量区域R7至R11,并且,对左心房LA与左心室LV的边界区域设定测量区域R5以及R6,根据在各个测量区域中包含的多个关心点上测量出的运动参数的平均值设定上述测量区域的亮度或色调等。
另外,图1所示的参数数据生成部3的时间序列数据生成部32通过在时间轴方向上排列跟踪处理部23从图像数据收集部1供给的规定心跳周期的时间序列的图像数据的心肌组织测量出的希望部位的局部运动参数或希望测量区域的平均运动参数,生成表示运动参数或平均运动参数的时间变化的参数时间序列数据,向显示部4供给所取得的规定心跳周期的参数时间序列数据。
图5为时间序列数据生成部32生成的参数时间序列数据的具体例子,在与图4所示的参数图像数据一起生成的该参数时间序列数据中,例如针对一个心跳周期描绘出在预先设定的左心房LA的测量区域R1以及R2中测量出的平均运动参数的时间变化TR1以及TR2、以及在预先设定的左心室LV的测量区域R7以及R8中测量出的平均运动参数的时间变化TR7以及TR8。
再次返回到图1,显示部4具备显示数据生成部41与监视器42,显示数据生成部41将图像数据收集部1供给的图像数据与参数数据生成部3供给的参数时间序列数据以及参数图像数据重叠或合成后转换成规定的显示格式,由此生成显示数据并在监视器42上加以显示。此时,参数图像数据通过与上述运动参数或平均运动参数对应的亮度和色调显示出来。
即,显示数据生成部41接收参数数据生成部3的图像数据生成部31供给的时间序列的参数图像数据以及时间序列数据生成部32供给的参数时间序列数据,从图像数据收集部1的存储电路中读出具有与附加在参数图像数据的各个中的心跳时相信息相同的心跳时相信息的图像数据。并且,合成与该图像数据重叠的参数图像数据与上述参数时间序列数据。
图6表示显示部4的监视器42上所显示的显示数据的具体例子。该显示数据,例如具有图像数据显示区域DA1与时间序列数据显示区域DA2。在图像数据显示区域DA1中显示在同一心跳时相与图像数据B1重叠的参数图像数据B2,在时间序列数据显示区域DA2中显示根据规定心跳周期的时间序列的图像数据生成的参数时间序列数据B3和基于附加在上述图像数据中的心跳时相信息的心电波形B4。此时,对参数时间序列数据B3设定表示配置在图像数据显示区域DA1中的图像数据B1以及参数图像数据B2的心跳时相的可动线。
例如,图像数据收集部1重复供给从被检体中收集到的1个心跳周期的图像数据时,在图像数据显示区域DA1中,根据这些图像数据生成的参数时间序列数据B2与上述图像数据一起作为运动画面重复显示(循环显示)出来。此时,表示图像数据显示区域DA1中显示的图像数据B1以及参数图像数据B2的心跳时相的可动线沿着时间序列数据显示区域DA2中的参数时间序列数据B3的时间轴向箭头方向移动。
其次,图1所示的输入部5在操作面板上具备键盘、轨迹球、鼠标、选择按钮、输入按钮等输入装置或显示面板,并且具有选择生成希望的图像数据的医用图像诊断装置(医用图像摄影装置)的医用图像摄影装置选择部51、选择运动参数的运动参数选择部52、对基准图像数据设定心内膜以及心外膜的内外膜设定部53等。另外,使用上述显示面板或输入装置也可以进行被检体信息的输入、图像数据收集条件的设定、参数数据生成条件的设定、参数数据显示条件的设定、各种指令信号的输入等。
另外,作为医用图像摄影装置,包括超声波诊断装置、X射线CT装置、MRI装置、X射线诊断装置、内窥镜装置等,作为运动参数,除了上述“应变”之外还包括“位移”、“旋转”、“扭转”、“速度”、“应变率”、“旋转率”、“扭转率”、“加速度”等。
其次,系统控制部6具备未图示的CPU与存储电路,由输入部5输入/设定/选择的信息被保存到上述存储电路中。并且,上述CPU根据输入部5输入的上述信息或预先保管到自己的存储电路中的信息,整体上控制医用图像解析装置100的各单元,从而进行参数图像数据以及参数时间序列数据的生成与显示。
(参数数据的生成/显示顺序)
其次,针对本实施例中的参数图像数据以及参数时间序列数据的生成/显示的顺序,使用图7的流程图进行说明。
在参数数据生成之前,医用图像解析装置100的操作者在输入部5中输入被检体信息后,作为收集图像数据的医用图像摄影装置,例如选择“X射线CT装置”,作为运动参数选择心肌组织的“应变”。然后,操作者在输入部5中设定图像数据收集条件、参数数据生成条件、参数数据或图像数据的显示条件等。而且,上述初始设定如果结束,则操作者在输入部5中输入图像数据的收集指令(图7的步骤S1)。
从输入部5经由系统控制部6接收上述被检体信息、医用图像摄影装置选择信息、图像数据收集条件信息以及图像数据收集开始指令的图像数据收集部1对经由网络10连接的X射线CT装置发送基于上述信息的指示信号。并且,根据该指示信号从X射线CT装置供给的该被检体的规定心跳周期的时间序列的图像数据(CT图像数据)与作为附带信息的心跳时相信息一起暂时保存到图像数据收集部1的存储电路中(图7的步骤S2)。
图像数据的收集与保存完成后,操作者在输入部5中输入参数数据的生成开始指令,经由系统控制部6接收该指令信号的运动参数测量部2的基准图像数据提取部21,作为基准图像数据从图像数据收集部1的存储电路中保存的时间序列图像数据中提取具有R波定时信息的图像数据,并在显示部4的监视器42上加以显示。
观测显示部4上显示的基准图像数据的操作者使用输入部5的内外膜设定部53对基准图像数据的心肌组织设定心内膜以及心外膜。其次,运动参数测量部2的关心点设定部22读出预先保管在跟踪条件保管部24中的X射线CT装置的跟踪条件,根据该跟踪条件,对由上述心内膜以及心外膜包围的基准图像数据的心肌组织以规定间隔设定多个关心点(图7的步骤S3)。另外,该设定也可以通过手动操作进行。
另外,跟踪处理部23读出图像数据收集部1的存储电路中保存的规定心跳周期的时间序列的图像数据、以及预先保管在跟踪条件保管部24中的X射线CT装置的跟踪条件。并且,根据上述跟踪条件对这些图像数据中包含的基准图像数据设定以上述关心点为中心的规定尺寸的模板。然后,进行该模板中包含的多个像素与上述基准图像数据的后续图像数据的像素之间的模式匹配,测量出心肌组织的局部“位移”,进而作为运动参数测量出被定义为每单位长度的位移量的“应变”(图7的步骤S4)。
其次,参数数据生成部3的图像数据生成部31通过排列运动参数测量部2的跟踪处理部23从时间序列的图像数据的各个测量出的运动参数,使其与关心点的位置信息对应,来生成2维参数图像数据,时间序列数据生成部32通过在时间轴方向依次排列跟踪处理部23从上述图像数据的心肌组织测量出的希望部位的局部运动参数或基于该运动参数的希望区域的平均运动参数,生成表示运动参数或平均运动参数的时间变化的参数时间序列数据(图7的步骤S5)。
并且,向显示部4供给从规定心跳周期的图像数据生成的参数图像数据以及参数时间序列数据。此时,参数图像数据的各个,向显示部4供给作为附加信息附加在与参数图像数据对应的图像数据中的心跳时相信息。
其次,显示部4的显示数据生成部41接收参数数据生成部3的图像数据生成部31供给的时间序列的参数图像数据以及时间序列数据生成部32供给的参数时间序列数据,从图像数据收集部1的存储电路中读出具有与各个参数图像数据的心跳时相信息相同的心跳时相信息的图像数据。并且,对与该图像数据重叠的参数图像数据附加参数时间序列数据,生成显示数据并在监视器42上加以显示(参照图6)(图7的步骤S6)。
按照上述顺序对于图像数据收集部1的存储电路中保存的规定心跳周期的图像数据的参数数据的生成与显示结束时,运动参数测量部2根据系统控制部6供给的控制信号重复读出上述规定心跳周期的时间序列的图像数据,从而测量出运动参数,参数数据生成部3生成基于运动参数测量部2测量出的运动参数的参数图像数据以及参数时间序列数据并在显示部4上加以显示。即,在显示部4上重复显示(循环显示)与规定心跳周期的图像数据对应的参数图像数据以及参数时间序列数据(图7的步骤S4至S6)
根据上述的本发明的实施例,通过对由不同的医用图像诊断装置的各个生成的生物体组织的时间序列的图像数据适用预先设定的共用的解析算法,测量出希望的运动参数,可以有效评价生物体组织的运动功能。因此,可以提高诊断效率与诊断精度,减轻对操作者的负担。
特别是,因为肋骨或肺等阻碍超声波传播而难以从超声波图像数据得到心脏的四腔剖面或两腔剖面的参数数据,而根据CT图像数据生成参数数据时,可以容易地生成心脏的四腔剖面或两腔剖面的参数数据。
另外,根据上述实施例,通过使用在沿着心肌组织设定的测量区域的多个关心点上测量出的多个运动参数的平均值(平均运动参数),可以抑制由噪音等引起的变动,从而可以生成对诊断有效的参数图像数据。并且,通过将参数图像数据与图像数据重叠并加以显示,可以明确相对于心脏的整个区域的参数数据的生成区域。
以上,针对本发明的实施例进行了叙述,但是并不只限于此,能够应变来实施。例如,在上述实施例中,针对通过使用适用了共用的解析算法的共用的医用图像解析装置100解析由不同的医用图像诊断装置生成的图像数据,来生成各种参数数据的情况进行了叙述,但是也可以将共用的解析算法适用于与上述医用图像诊断装置对应的专用的医用图像解析装置来生成参数数据。
另外,在上述实施例中,针对根据从图像数据收集部1的存储电路重复供给的规定心跳周期的图像数据生成的时间序列的参数图像数据在显示部4上循环显示的情况进行了叙述,但是可以设置保存根据上述规定心跳周期的图像数据生成的参数数据的参数数据存储部,可以重复读出该参数数据存储部中保存的参数图像数据并循环显示出来。
并且,运动参数的测量中的处理速度比图像数据的收集中的处理速度快时,可以从医用图像诊断装置依次供给的图像数据生成参数数据,也可以实时显示所取得的参数数据。此时,未必需要图1的图像数据收集部1具备的存储电路或上述参数数据存储部,
另一方面,在上述实施例中,针对根据运动参数测量部2的跟踪条件保管部24中保管的跟踪条件对时间序列的图像数据进行跟踪处理的情况进行了叙述,但是并不只限于此,例如,在输入部5中也可以直接设定适合该图像数据的跟踪处理的跟踪条件。
另外,针对根据医用图像诊断装置供给的2维图像数据生成2维参数图像数据的情况进行了叙述,但是也可以根据3维图像数据生成2维或3维的参数图像数据。另外,针对将附加了R波定时信息的图像数据作为基准图像数据来设定关心点的情况进行了叙述,但是也可以将在任意心跳时相得到的图像数据作为基准图像数据。
并且,在上述实施例中,针对作为运动参数测量出心肌组织的“应变”的情况进行了叙述,但是并不仅限于此,例如,作为运动参数也可以测量出“位移”、“旋转”、“扭转”、“速度”等。作为运动参数也可以测量出表示它们的时间变化的“应变率”、“旋转率”、“扭转率”、“加速度”等。在心肌组织以外的生物体组织中也可以测量上述运动参数。
另外,针对重叠或合成参数图像数据、参数时间序列数据以及图像数据并加以显示的情况进行了叙述,但是也可以单独地显示这些数据。针对经由网络10收集医用图像诊断装置生成的时间序列的图像数据的情况进行了叙述,但是也可以经由存储媒体进行收集,其收集方法没有特别限定。
另外,由图1所示的运动参数测量部2以及参数数据生成部3的各单元进行的跟踪处理或参数数据生成处理也可以通过硬件进行,但是这些处理的全部或其一部分,通常通过基于共用的解析算法的软件进行。

Claims (13)

1.一种医用图像解析装置,其特征在于,包括:
图像数据收集部,收集由不同的医用图像诊断装置的各个生成的时间序列的图像数据;
运动参数测量部,通过对上述时间序列的图像数据进行处理,测量出生物体组织的运动参数;
参数数据生成部,根据上述运动参数生成参数数据;和
显示部,显示上述参数数据;
其中,上述运动参数测量部通过使用共用的解析算法解析从上述医用图像诊断装置的各个收集到的上述图像数据,测量出上述运动参数。
2.根据权利要求1所述的医用图像解析装置,其特征在于:
上述运动参数测量部通过跟踪由上述医用图像诊断装置生成的上述时间序列的图像数据,测量出上述运动参数。
3.根据权利要求2所述的医用图像解析装置,其特征在于:
上述运动参数测量部根据与上述医用图像诊断装置对应的跟踪条件来跟踪上述图像数据。
4.根据权利要求3所述的医用图像解析装置,其特征在于:
上述跟踪条件为模板尺寸、观察区域、动态范围、增益和平滑处理中的至少一个。
5.根据权利要求3所述的医用图像解析装置,其特征在于,还包括:
跟踪条件保管部,预先保管与上述不同的医用图像诊断装置中的各个对应的跟踪条件,
其中,上述运动参数测量部从上述跟踪条件保管部中读出与生成上述时间序列的图像数据的上述医用图像诊断装置对应的上述跟踪条件,根据该跟踪条件跟踪上述图像数据。
6.根据权利要求3所述的医用图像解析装置,其特征在于,还包括:
跟踪条件保管部,预先保管与上述不同的医用图像诊断装置中的各个对应的跟踪条件;和
医用图像摄影装置选择部,选择生成上述时间序列的图像数据的上述医用图像诊断装置,
其中,上述运动参数测量部根据上述医用图像摄影装置选择部供给的医用图像摄影装置选择信息,使用与从上述跟踪条件保管部中读出的上述医用图像诊断装置对应的跟踪条件,跟踪上述图像数据。
7.根据权利要求3所述的医用图像解析装置,其特征在于,还包括:
跟踪条件输入部,输入与上述医用图像诊断装置对应的跟踪条件,
其中,上述运动参数测量部根据由上述跟踪条件输入部输入的上述跟踪条件来跟踪上述图像数据。
8.根据权利要求1所述的医用图像解析装置,其特征在于:
上述参数数据生成部包括根据上述运动参数测量部测量出的上述运动参数生成2维或3维参数图像数据的图像数据生成部。
9.根据权利要求1所述的医用图像解析装置,其特征在于:
上述参数数据生成部包括生成表示上述运动参数测量部测量出的上述运动参数的时间变化的参数时间序列数据的时间序列数据生成部。
10.根据权利要求1所述的医用图像解析装置,其特征在于:
上述显示部实时显示根据上述时间序列的图像数据生成的上述参数图像数据。
11.根据权利要求8所述的医用图像解析装置,其特征在于:
上述显示部循环显示根据上述时间序列的图像数据生成的规定的心跳周期的上述参数图像数据。
12.根据权利要求10或11所述的医用图像解析装置,其特征在于:
上述显示部将上述参数图像数据与上述图像数据重叠或合成,并同步显示。
13.一种图像解析用控制程序,其特征在于,使医用图像解析装置执行以下功能:
收集由不同的医用图像诊断装置的各个生成的时间序列的图像数据的功能;
通过对上述时间序列的图像数据进行处理来测量出生物体组织的运动参数的功能;
根据上述运动参数生成参数数据的功能;和
显示上述参数数据的功能。
CN200910208856.7A 2008-11-05 2009-11-05 医用图像解析装置以及图像解析用控制方法 Active CN101732036B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US12/265,302 US8208703B2 (en) 2008-11-05 2008-11-05 Medical image analysis apparatus and image analysis control program
US12/265,302 2008-11-05

Publications (2)

Publication Number Publication Date
CN101732036A true CN101732036A (zh) 2010-06-16
CN101732036B CN101732036B (zh) 2014-06-25

Family

ID=42131455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910208856.7A Active CN101732036B (zh) 2008-11-05 2009-11-05 医用图像解析装置以及图像解析用控制方法

Country Status (3)

Country Link
US (1) US8208703B2 (zh)
JP (1) JP5491747B2 (zh)
CN (1) CN101732036B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102293669A (zh) * 2010-06-25 2011-12-28 株式会社东芝 超声波诊断装置及超声波图像处理装置
CN103156647A (zh) * 2011-12-08 2013-06-19 株式会社东芝 超声波诊断装置、图像处理装置及图像处理方法
CN104067277A (zh) * 2012-01-27 2014-09-24 皇家飞利浦有限公司 医学选择系统
CN108652670A (zh) * 2018-03-19 2018-10-16 武汉理工大学 左心室功能异常超声成像治疗一体化设备

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4727737B2 (ja) * 2009-02-24 2011-07-20 三菱重工業株式会社 放射線治療装置制御装置および目的部位位置計測方法
DE102009021492B4 (de) * 2009-05-15 2012-01-12 Siemens Aktiengesellschaft Verfahren zur Bestimmung einer MR-Relaxationszeit im Herzmuskel bei einer Magnetresonanzuntersuchung, Magnetresonanzanlage, Computerprogrammprodukt und elektronisch lesbarer Datenträger
US8861830B2 (en) * 2011-11-07 2014-10-14 Paieon Inc. Method and system for detecting and analyzing heart mechanics
US20150042677A1 (en) * 2012-03-23 2015-02-12 Konica Minolta, Inc. Image-generating apparatus
BR112014032101A2 (pt) * 2012-06-27 2017-06-27 Koninklijke Philips Nv computador configurado para estimar um parâmetro de movimentos para o movimento de um objeto móvel, determinador paramétrico de movimento de um objeto, e método automatizado realizado por computador
WO2014054660A1 (ja) * 2012-10-01 2014-04-10 株式会社東芝 画像処理装置及びx線ct装置
KR102094501B1 (ko) * 2012-11-30 2020-03-27 삼성전자주식회사 Mri 다변수 영상의 획득 장치 및 방법
KR20140088434A (ko) 2013-01-02 2014-07-10 삼성전자주식회사 환자별 특성을 기반으로 한 mri 다변수 영상 획득 지원 장치 및 방법
CN106030657B (zh) * 2014-02-19 2019-06-28 皇家飞利浦有限公司 医学4d成像中的运动自适应可视化

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3715580B2 (ja) * 2002-03-19 2005-11-09 株式会社東芝 医用運動解析装置及びその方法
WO2004092766A1 (en) * 2003-04-15 2004-10-28 Koninklijke Philips Electronics N.V. Heart wall strain imaging
JP2005130877A (ja) 2003-10-28 2005-05-26 Aloka Co Ltd 超音波診断装置
US7106189B2 (en) * 2004-04-29 2006-09-12 Tracetech Incorporated Tracking system and methods thereof
JP4733938B2 (ja) * 2004-07-16 2011-07-27 株式会社東芝 超音波診断装置および超音波画像処理装置
JP2006110190A (ja) * 2004-10-15 2006-04-27 Toshiba Corp 医用画像データ解析装置及びその方法
JP4714000B2 (ja) * 2005-10-31 2011-06-29 株式会社東芝 心機能解析装置及びその方法
GB2438461A (en) * 2006-05-23 2007-11-28 Univ Cambridge Tech Weighted phase separation based ultrasonic deformation estimation
JP4745133B2 (ja) * 2006-05-30 2011-08-10 株式会社東芝 超音波診断装置、医用画像処理装置及び医用画像処理プログラム
US8238624B2 (en) * 2007-01-30 2012-08-07 International Business Machines Corporation Hybrid medical image processing
JP5624258B2 (ja) * 2007-04-26 2014-11-12 株式会社東芝 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102293669A (zh) * 2010-06-25 2011-12-28 株式会社东芝 超声波诊断装置及超声波图像处理装置
US9161737B2 (en) 2010-06-25 2015-10-20 Kabushiki Kaisha Toshiba Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus
CN103156647A (zh) * 2011-12-08 2013-06-19 株式会社东芝 超声波诊断装置、图像处理装置及图像处理方法
CN103156647B (zh) * 2011-12-08 2016-03-16 株式会社东芝 超声波诊断装置、图像处理装置及图像处理方法
US9345446B2 (en) 2011-12-08 2016-05-24 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus, image processing apparatus, and image processing method
CN104067277A (zh) * 2012-01-27 2014-09-24 皇家飞利浦有限公司 医学选择系统
CN110085302A (zh) * 2012-01-27 2019-08-02 皇家飞利浦有限公司 医学选择系统
CN108652670A (zh) * 2018-03-19 2018-10-16 武汉理工大学 左心室功能异常超声成像治疗一体化设备

Also Published As

Publication number Publication date
JP2010110612A (ja) 2010-05-20
JP5491747B2 (ja) 2014-05-14
US20100111380A1 (en) 2010-05-06
CN101732036B (zh) 2014-06-25
US8208703B2 (en) 2012-06-26

Similar Documents

Publication Publication Date Title
CN101732036B (zh) 医用图像解析装置以及图像解析用控制方法
KR102269467B1 (ko) 의료 진단 이미징에서의 측정 포인트 결정
CN100393283C (zh) 生物体组织动状态跟踪方法、使用该方法的图像诊断装置
CN100430025C (zh) 超声波诊断装置和图象处理装置
US8172756B2 (en) Ultrasonic imaging system and method
CN101336844B (zh) 医用图像处理装置以及医用图像诊断装置
US6994673B2 (en) Method and apparatus for quantitative myocardial assessment
CN102047140B (zh) 具有引导efov扫描的扩展视野超声成像
CN102348417B (zh) 心脏m模式图片的自动分析
JP5661453B2 (ja) 画像処理装置、超音波診断装置、及び画像処理方法
US6749571B2 (en) Method and apparatus for cardiac elastography
CN101066210B (zh) 超声系统中显示信息的方法
CN107072635B (zh) 用于中间用户反馈的多跳超声心动图采集的质量度量
CN109758178A (zh) 超声成像中的机器辅助工作流
US20100016718A1 (en) Shear Wave Imaging
CN101461723B (zh) 超声波诊断装置、超声波图像解析装置以及超声波图像解析方法
CN101484074A (zh) 超声波诊断装置、医用图像处理装置以及医用图像处理程序
CN102056547A (zh) 医用图像处理装置及医用图像处理方法
CN102939050B (zh) 超声波诊断装置以及超声波收发方法
CN101681516A (zh) 标记超声成像系统的2d显示器上的3d体积图像的系统和方法
DE102017204471A1 (de) Unterstützende Kennzeichnung für Ultraschallbildgebung im Kontrollmodus
US11903760B2 (en) Systems and methods for scan plane prediction in ultrasound images
US20140153358A1 (en) Medical imaging system and method for providing imaging assitance
CN107481259A (zh) 尤其是在超声波空间复合中用于估计图像间运动的方法和系统
CN106030657A (zh) 医学4d成像中的运动自适应可视化

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant