CN102293669B - 超声波诊断装置及超声波图像处理装置 - Google Patents

超声波诊断装置及超声波图像处理装置 Download PDF

Info

Publication number
CN102293669B
CN102293669B CN2011101726375A CN201110172637A CN102293669B CN 102293669 B CN102293669 B CN 102293669B CN 2011101726375 A CN2011101726375 A CN 2011101726375A CN 201110172637 A CN201110172637 A CN 201110172637A CN 102293669 B CN102293669 B CN 102293669B
Authority
CN
China
Prior art keywords
mentioned
value
unit
time changing
relevant
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN2011101726375A
Other languages
English (en)
Other versions
CN102293669A (zh
Inventor
阿部康彦
川岸哲也
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
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
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 Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of CN102293669A publication Critical patent/CN102293669A/zh
Application granted granted Critical
Publication of CN102293669B publication Critical patent/CN102293669B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0883Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Cardiology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明提供一种超声波诊断装置,给出对于CRT应用判断有用的与非同步性的程度相关的指标及非同步性的评价法。该超声波诊断装置具备:数据取得单元,将包括被检体的心脏的至少一部分的二维区域或三维区域作为被扫描区域进行超声波扫描,取得与该被扫描区域相关的超声波数据;时间变化曲线取得单元,使用与上述被扫描区域相关的超声波数据,取得与心脏的多个部位相关的壁运动参数的时间变化曲线;指标值取得单元,使用与上述心脏的多个部位中的任意两个部位相关的时间变化曲线,取得用来评价心脏的壁运动的能力的指标值;以及显示单元,将上述指标值以规定的形态显示。

Description

超声波诊断装置及超声波图像处理装置
本申请基于2010年6月25日提出的日本专利申请第2010-145519号并主张其优先权,这里引用其全部内容。
技术领域
本发明涉及超声波诊断装置及超声波图像处理装置。
背景技术
关于心肌等的生物体组织,客观且定量地评价其功能对于该组织的诊断是非常重要的。例如,近年来,对于重度心脏功能不全患者的心脏再同步疗法(CRT:Cardiac Resynchronization Therapy)受到关注,为了向CRT的应用的事前判断及治疗的效果判断,通过心回波进行的定量评价变得重要。CRT是能够改善归并到重度心脏功能不全患者中的情况较多的心壁运动的收缩不同步(dyssynchorny)的治疗法。在该治疗法有效的患者(Responder)中,可以确认戏剧性的症状的改善。另一方面,虽然是心脏功能不全症,但CRT无效的患者(Non Responder)在2005年时存在约3成而成为问题。
CRT无效的患者是心脏功能不全中的不是收缩不同步的患者。到目前为止,用心电图的QRS宽度>130msec、射血分数(Ejection Fraction、即EF)≤35%判断CRT的应用。但是,在该基准中,也包括虽然心脏功能不全但不是收缩不同步的人,产生了CRT无效的患者。
所以,在世界中开始想要通过使用了心回波的定量评价法仅提取收缩不同步的尝试,提出了各种使用了速度的到达时间成像、位移或应变的到达时间(峰值时间或重心时间)成像等的方法。这些方法以能够将局部的心肌的收缩定时的差异作为彩色图像输出、容易进行收缩不同步的有无及收缩定时异常的部位的掌握为目标。这里,作为收缩定时的差异,最一般地采用所谓速度、位移及应变(strain)的“关于壁运动参数的时间变化的峰值(peak)时间差”。
但是,PROSPECT study(Circulation.2008;117:2608-16)的结果报告出,通过速度及位移的心回波多普勒法的指标即峰值时间差,不能显著地识别出CRT有效的患者和CRT无效的患者。
另一方面,不仅是峰值时间,还提出了通过重心时间的非同步性的评价。即,鉴于如果是峰值则在速度及位移的峰值时相附近混入一些噪声的情况、或存在多个峰值的情况下的非同步性的评价不稳定而提出的是重心时间。但是,即使是重心时间,结果与上述峰值时间差同样,在检测对于CRT应用有效的非同步性(异常)的方面存在制约。此外,一般而言,重心时间与峰值时间相比正常与异常的时间差较小,所以正常与异常的分离较困难。
因而,通过以往的关于非同步性的指标,不能显著地识别CRT的有效的患者和无效的患者。
专利文献1:日本特许第3187008号公报
专利文献2:日本特愿2004-173589号
专利文献3:日本特愿2005-073788号
发明内容
本发明是鉴于上述情况而做出的,目的是提供一种给出对于CRT应用判断有用的关于非同步性的程度的指标及非同步性的评价法的超声波诊断装置及超声波图像处理装置。
有关技术方案的超声波诊断装置具备:数据取得单元,将包括被检体的心脏的至少一部分的二维区域或三维区域作为被扫描区域进行超声波扫描,取得与该被扫描区域相关的超声波数据;时间变化曲线取得单元,使用与上述被扫描区域相关的超声波数据,取得与心脏的多个部位相关的壁运动参数的时间变化曲线;指标值取得单元,使用与上述心脏的多个部位中的任意两个部位相关的时间变化曲线,取得用来评价心脏的壁运动的能力的指标值;以及显示单元,将上述指标值以规定的形态显示。
发明效果
以上,根据本发明,能够实现给出对于CRT应用判断有用的关于非同步性的程度的指标及非同步性的评价法的超声波诊断装置及超声波图像处理装置。
附图说明
图1是有关实施方式的超声波诊断装置1的结构图。
图2A~图2D是表示不同的两个部位的壁运动参数的时间变化曲线f(t)、g(t)的4种情况的图。
图3是表示遵循有关第1实施方式的非同步性评价支援功能的处理(非同步性评价支援处理)的流程的流程图。
图4是表示对于CRT有效的左束支传导阻滞病例的本非同步性评价支援功能的应用例的图。
图5是表示对健康例的本非同步性评价支援功能的应用例的图。
图6是表示对健康例的本非同步性评价支援功能的应用例的图。
图7是表示非同步性评价支援信息的显示形态的另一例的图。
图8是表示有关第2实施方式的非同步性评价支援处理的流程的流程图。
标记说明
1超声波诊断装置,11超声波探头,13发送单元,15接收单元,17B模式处理单元,19移动矢量处理单元,21图像生成单元,23显示单元,31控制单元(CPU),37壁运动参数运算单元,38DI值计算单元,39存储单元,41操作单元,43收发单元
具体实施方式
有关实施方式的超声波诊断装置具备:数据取得单元,将包括被检体的心脏的至少一部分的二维区域或三维区域作为被扫描区域进行超声波扫描,取得与该被扫描区域相关的超声波数据;时间变化曲线取得单元,使用与上述被扫描区域相关的超声波数据,取得与心脏的多个部位相关的壁运动参数的时间变化曲线;指标值取得单元,使用与上述心脏的多个部位中的任意两个部位相关的时间变化曲线,取得用来评价心脏的壁运动的能力的指标值;以及显示单元,将上述指标值以规定的形态显示。
以下,按照附图说明本发明的第1实施方式及第2实施方式。另外,在以下的说明中,对于具有大致相同的功能及结构的构成单元赋予相同的标记,只在需要的情况下进行重复的说明。
另外,在以下的实施方式中,在超声波诊断装置中针对应用例进行说明。但是,并不局限于此,也可以对工作站、个人计算机等的超声波处理装置应用。
(第1实施方式)
图1是有关第1实施方式的超声波诊断装置1的结构图。本超声波诊断装置1具备超声波探头11、发送单元13、接收单元15、B模式处理单元17、移动矢量处理单元19、图像生成单元21、显示单元23、控制单元(CPU)31、壁运动参数运算单元37、DI值计算单元38、存储单元39、操作单元41、收发单元43。另外,在应用在超声波图像处理装置中的情况下,图1的虚线内为其构成单元。
超声波探头11具有基于来自发送单元13的驱动信号产生超声波,将来自被检体的反射波变换为电信号的多个压电振子、设在该压电振子上的调整层、以及防止从该压电振子向后方的超声波的传播的垫材(バッキング材)等。如果从该超声波探头11向被检体发送超声波,则通过生物体组织的非线性等,随着超声波的传播而发生各种谐波成分。构成发送超声波的基波和谐波成分因体内组织的声阻抗的边界、微小散射等而被向后方散射,作为反射波(回波)被超声波探头11接收。
发送单元13具有未图示的延迟电路及脉冲发生器电路等。在脉冲发生器电路中,以规定的频率f Hz(周期为1/f秒)反复产生用来形成发送超声波的速率脉冲。此外,在延迟电路中,按照每个通道将超声波集束为束状,并且对各速率脉冲施加决定发送指向性所需要的延迟时间。发送单元13以基于该速率脉冲的定时以朝向规定的扫描线形成超声波束的方式对每个振子施加驱动脉冲。
接收单元15具有未图示的放大电路、A/D转换器、加法器等。在放大电路中,将经由探头11取入的回波信号按照每个通道放大。在A/D转换器中,对放大后的回波信号赋予决定接收指向性所需要的延迟时间,然后在加法器中进行相加处理。通过该相加,生成对应于规定的扫描线的超声波回波信号。
B模式处理单元17通过对从接收单元15接受到的超声波回波信号实施包络线检波处理,生成对应于超声波回波的振幅强度的B模式信号。
移动矢量处理单元19在时相不同的两个超声波数据间使用图案匹配处理检测组织的移动位置,基于该移动位置求出各组织的移动量(或速度)。具体而言,对于一个超声波数据内的关注区域,求出类似性最高的另一个超声波数据内的对应区域。通过求出该关注区域与对应区域之间的距离,能够求出组织的移动量。此外,通过将该移动量用超声波数据间的时间差(帧速率或体速率(ボリュームレート))除,能够求出组织的移动速度。通过将该处理在超声波数据上的各位置进行,能够取得各组织的位移(移动矢量)或关于组织的位移的时空间分布数据。另外,这里将超声波数据定义为具有二维或三维的位置信息的接收信号的集合(即具有空间信息的接收信号的集合)。
图像生成单元21生成表示有关B模式信号的规定断层的二维分布的B模式超声波像。此外,图像生成单元21基于运算出的壁运动参数,使用表面绘制、极映射(Polar-Mapping)等的方法生成映射了该壁运动参数的二维图像或三维图像。
显示部23基于来自图像生成单元21的视频信号,将超声波图像、在对应的各位置上映射了壁运动参数的壁运动参数图像、各部位的壁运动参数的时间变化曲线等用规定的形态显示。此外,显示部23将通过分段处理分割后的各部位按照后述的非同步评价支援功能标记,以规定的形态显示。
控制单元(CPU)31具有作为信息处理装置(计算机)的功能,静态或动态控制本超声波诊断装置主体的动作。特别是,控制单元31通过将存储在存储单元39中的专用程序展开到未图示的存储器中,实现后述的非同步评价支援功能。
壁运动参数运算单元37基于移动矢量处理单元19输出的时空间分布数据,按照时相生成壁运动参数。这里,所谓壁运动参数,例如是心壁等的规定组织的与规定方向相关的位移、应变、应变率、速度、扭转、扭转率、及关于其他组织运动能够取得的物理信息。
DI值计算单元38按照后述的非同步评价支援功能,计算用来在非同步性的评价中使用的DI值。
存储单元39是磁盘(软盘(注册商标)、硬盘等)、光盘(CD-ROM、DVD等)、半导体存储器等的记录介质、以及将记录在这些介质中的信息读出的装置。在该存储单元37中,存储收发条件、规定的扫描顺序、对应于各时相的原始数据及超声波图像数据(例如通过组织多普勒模式、B模式等摄影的组织图像数据)、预先生成的各时相的超声波数据、关于移动矢量的时空间分布数据、用来实现后述的非同步评价支援功能的程序、诊断信息(患者ID、医生的所见等)、诊断草案、身体标记生成程序等。
操作单元41连接在装置主体上,具有用于进行来自操作者的各种指示、关注区域(ROI)的设定指示、各种画质条件设定指示、后述的非同步评价支援功能中的任意的壁运动参数及任意的解析期间、心时相的选择等的鼠标或跟踪球、模式切换开关、键盘等。
收发单元43是经由网络与其他装置进行信息的收发的装置。能够将在本超声波诊断装置1中得到的超声波图像等的数据及解析结果等通过收发单元43经由网络向其他装置转送。
(理论背景)
CRT为有效的类型(Responder)的心壁运动的非同步被指出部位间的峰值时间差相同时,部位间的壁运动参数值的差异更大。例如,如果存在因心肌梗塞造成的局部性的壁运动下降(akinesis)或心室壁的纤维化,则在该部位和其他健康的部位中容易发生壁运动参数的峰值时间差,但在这样的病变部位中,即使通过CRT起搏,也因为没有心肌的余力而发生心室收缩不全,结果成为CRT无效的患者。在该例中,仅通过峰值时间差并不清楚部位间的壁运动参数值的差异的程度是本质性的问题。
图2A~图2D是表示不同的两个部位的壁运动参数的时间变化曲线(例如应变波形)f(t)、g(t)的4种情况的图(纵轴表示应变,横轴表示时间)。在该图中,如果将例如图2A的情况与图2B的情况比较,则2部位间的峰值时间差相同,但在图2A的情况比图2B的情况CRT更有效的观点中,非同步性更高。因而,仅通过峰值时间差并不能区别CRT相对有效的图2A的情况和图2B的情况。同样,如果将图2C的情况与图2D的情况比较,则在图2C的情况比图2D的情况CRT更有效的观点中,非同步性较高。这里,在图2C的情况中,两部位间的重心时间差相等,预先附加了通过重心时间不能检测到本例的非同步性的情况。
将以往的以峰值时间差及重心时间差为代表的部位间的非同步性称作局部心肌的收缩不同步(dyssynchrony)。另一方面,将考虑到这样的壁运动参数值的部位间差异的非同步性的概念称作“左心室的失调(discoordination)”。并且,认为左心室的失调的程度越大则CRT越有效。并且,部位间的壁运动参数值的差异如果用时间变化曲线观察,则相当于由对应于各部位的两个时间变化曲线包围的区域的面积值(面积差分值),可以认为反映左心室的失调。
(非同步性评价支援功能)
接着,对本超声波诊断装置1具有的非同步性评价支援功能进行说明。本功能作为定量地反映左心室的失调的指标而导入用来评价心脏的壁运动的能力的指标值。作为该指标值的代表性的例子,计算将不同的两个部位(例如通过分段将心壁分割的情况下的不同的两个段)间的面积差分值用规定的值标准化而得到的DI值(Discoordinatin Index)。通过将该DI值作为非同步性评价支援信息以规定的形态显示,能够支援心壁运动的非同步性的评价。
另外,在以下的说明中,所谓段,与局部部位是同义的。此外,在以下的说明中,以将包括左心室的二维区域作为被扫描区域而应用该非同步性评价支援功能的情况为例。但是,并不局限于该例,也可以将例如包括左心室的三维区域作为被扫描区域而应用该非同步性评价支援功能。
图3是表示遵循有关本实施方式的非同步性评价支援功能的处理(非同步性评价支援处理)的流程的流程图。以下,按照该图,对非同步性评价支援处理的各步骤的处理的内容进行说明。
[持续规定期间的时间序列体数据组的取得:步骤S1]
首先,将至少包括被检体的左心室的二维区域作为被扫描区域,取得持续至少一个心搏量以上的规定期间的时间序列的超声波数据(以下称作“时间序列超声波数据组”)(步骤S1)。
[接受壁运动参数的选择:步骤S2]
接着,控制单元31从操作者经由操作单元41受理壁运动参数的选择(步骤S2)。另外,在本实施方式中,为了具体地进行说明,假设在该步骤S2中选择了应变(strain)作为壁运动参数。但是,并不局限于该例,例如可以采用应变率(strain rate)、位移、速度等作为壁运动参数。
[壁运动参数组的生成:步骤S2]
接着,移动矢量处理单元19在构成所取得的持续规定期间的时间序列超声波数据组的1个心搏以上的对应于各时相的超声波数据中的、规定的时相的超声波数据中,基于来自用户的指示等提取心肌部位,将所提取的局部的心肌部位通过二维的图案匹配处理执行在时间上追踪的波谱跟踪处理,由此,运算时空间的移动矢量信息。此外,壁运动参数运算单元37使用运算出的时空间的移动矢量信息,运算由操作者经由操作单元41指定的壁运动参数,生成由1个心搏以上的壁运动参数构成的二维的壁运动参数组(步骤S3)。
[分段处理/时间变化曲线的计算:步骤S4]
接着,DI值计算单元38将在各时相的超声波数据中提取的心肌部位分割为N个局部部位(段)(分段处理),按照各局部部位计算壁运动参数的时间变化曲线(步骤S4)。另外,接着作为局部部位的数量N,例如在作为分割对象的数据是二维数据的情况下6左右是适当的。此外,在作为分割对象的数据是体数据的情况下,遵循ASE的16或遵循AHA的17是适当的。在本实施方式中,为了具体地进行说明,以N=6的情况为例。
[关于全部组合的DI值计算:步骤S5]
接着,DI值计算单元38关于不同的两个局部部位的全部组合,计算DI值(步骤S5)。即,DI值计算单元38在将任意的两个局部部位Si、Sj(其中,i≠j,i及j是满足1≤i、j≤6(分段数)的整数)的时间变化曲线分别设为f(t)、g(t)的情况下,按照以下的式(1)计算关于壁运动参数的DI值。
DI值=(Σ|f(t)-g(t)|/norm)·dt  (1)
另外,Σ意味着解析期间(例如1个心周期)中的时间和的运算。此外,按照该式(1)关于不同的两个局部部位的全部组合计算DI值的结果是,得到的DI值为
Figure GDA00003531890600081
个(在此情况下是
Figure GDA00003531890600082
个)。
在上述式(1)中,Σ|f(t)-g(t)|只能是不同的两个局部部位间的壁运动参数f(t)、g(t)的曲线的面积差分值。此外,norm是为了吸收局部部位间的时间变化曲线的值的波动(例如如果是健康例则strain等的壁运动参数值在段间具有显著误差)及因病态带来的曲线绝对值的波动(例如在健康例中壁运动参数的绝对值相对较大,但在心脏功能不全等的病例中绝对值变小)等的影响而需要的标准化值。该标准化值优选的是使用N段整体的平均曲线(全局(global))的峰值作为具有与分子相同单位的量。如果使用该值,则在前面的标准化的效果以外,还具有如下效果:对于相对于整体具有相对较小的参数值的段的DI值的加权相对减小,不易反映难以对CRT反应的壁运动降低部位的影响。作为其他标准化值,可以举出f(t)或g(t)的峰值。它们都具有用成对的一个参数值将另一个标准化的意义。进而,上述式(1)的dt意味着规定的时间间距。乘以时间间距dt意味着将曲线的面积差分的大小换算(时间积分)为时间值。本项并不是必须的结构。但是,如果将DI值的单位换算为时间,则具有能够与以往的峰值时间等比较的优点。
另外,标准化的方法并不局限于上述例子,例如如果使用能够按照每个壁运动参数设定的规定的值也能够执行。
此外,作为将上述式(1)的面积差分积分的时间区间(解析期间),在壁运动参数如位移或应变那样在健康例中为单峰性的物理量中,如上述那样使用1个心周期是优选的。另一方面,在壁运动参数如速度或应变速率那样在健康例中在收缩期和扩张期中具有极性不同的峰值的物理量中,优选的是将收缩期和扩张期分开设为解析期间。
[最大DI值及对应于DI值的部位的标识显示:步骤S6]
接着,DI值计算单元38从关于计算出的全部组合的DI值之中决定最大DI值(代表值)。显示单元23将与所决定的最大DI值对应的两个局部部位中的至少一个、优选的是将标识化了两者的壁运动参数图像、时间变化曲线与最大DI值一起作为非同步性评价支援而用规定的形态显示(步骤S6)。另外,代表值并不局限于最大DI值,也可以是其他规定的统计值。
以下,表示本非同步性评价支援功能的几个应用例。哪个例子都是根据二维的短轴图像用二维的波谱跟踪的方法求出的作为壁运动参数的应变计算DI值的。此外,作为计算DI值时的norm,使用N段整体的平均曲线的峰值。在这些各应用例中,对于段1~6,1:inf、2:pst、3:lat、4:ant、5:sp、6:asp的心脏的短轴像的解剖学的部位(通过分段的局部部位)分别对应。因而,DI值只分为不同的两个局部部位的组合1-2、1-3、1-4、……、5-6的共计15对计算。在本实施方式的DI值的定义中,由于非同步性的程度越大、DI值也越大,所以寻找DI值为最大的对。在各应用例中,作为优选的一个显示形态,在得到的最大DI值上同时标记显示以往的峰值时间差。
(应用例1)
图4是表示对于CRT有效的左束支传导阻滞病例的本非同步性评价支援功能的应用例的图。壁运动参数使用作为关于向壁厚方向的变化的壁运动参数的径向应变(Radial-Strain)。由该图可知,相对于最大峰值时间差是asp-inf部位间的414[msec],表述为最大面积差分时间的DI值在asp-lat部位间为997[msec]。此外,将asp和lat的部位在图像上用彩色的椭圆显示标识。
通过观察这样提供的非同步性评价支援信息,操作者能够根据曲线图的结果容易且迅速地掌握asp部位提前收缩、lat部位最具有非同步性的情况,可知其非同步性的程度比以往的峰值时间差更大地表现。此外,能够将这些各部位在超声波图像上容易地辨识。在CRT中,将左心室侧的起搏导线植入到非同步性最大的部位中可以说是有效的。因而,能够与通过DI值的非同步性的程度同时也同时掌握应植入起搏导线的位置的信息。
(应用例2)
图5是表示对健康例的本非同步性评价支援功能的应用例的图。壁运动参数使用作为关于向壁厚方向的变化的壁运动参数的径向应变。最大峰值时间差是1-5部位间的70[msec],最大重心时间差是1-4部位间的86[msec],相对于此,表述为最大面积差分时间的DI值在2-5部位间为72[msec]。
此外,在本应用例中,将局部部位2与局部部位5的凡例在曲线图上用彩色的椭圆显示,并且将对应的波形通过强调显示而标识。DI值具有以往的指标同等的程度,并且根据对于左束支传导阻滞病例的上述结果例,启发为是比以往的指标高灵敏度的非同步性的指标。
(应用例3)
图6是表示对健康例的本非同步性评价支援功能的其他的应用例的图。壁运动参数使用关于向圆周方向的变化的壁运动参数圆周应变(Circumferential-Strain)。最大峰值时间差是3-5部位间的94[msec]、最大重心时间差是3-4部位间的87[msec],相对于此,表述为最大面积差分时间的DI值在3-6部位间为75[msec]。
此外,在本应用例中,将局部部位3与局部部位6的凡例在曲线图上用彩色的椭圆显示,并且将对应的波形通过强调显示而标识。DI值具有与以往的指标同样的程度,并且是最小的。在本例的情况下,启发为DI值是比以往的指标更高特异度的非同步性的指标。
(变形例)
图7是表示例如对于图4所示的非同步性评价支援信息、改变了得到最大DI值的部位的显示方式的例子的图。在按照该例显示的非同步性评价支援信息中,在asp-lat部位间作为最大DI值而得到997[ms],但将asp和lat部位的位置分别在图像上用强调的框线标识显示。
根据本超声波诊断装置,作为反映左心室的失调的指标,对心脏的不同的两个部位(局部部位)的全部组合,计算将心脏的不同的两个部位间的面积差分值用规定的值标准化而得到的DI值,从其中判断成为最大的DI值。由于认为局部部位间的壁运动参数的面积差分值反映左心室的失调,所以得到的该最大DI值可以成为反映左心室的非同步性的指标。因而,如果使用该DI指标,则能够包罗对于CRT应用判断有用的非同步性的类型,与通过以往的指标进行的非同步性的评价相比,能够实现更高精度的CRT应用判断。
此外,根据本超声波诊断装置,与得到的最大DI值一起,将标识了与该对大DI值对应的两个局部部位的至少一个的壁运动参数图像、时间变化曲线与最大DI值一起作为非同步性评价支援而以规定的形态显示。因而,操作者通过观察显示的非同步性评价支援信息,能够容易且迅速地辨识出心脏的哪两个局部部位最具有非同步性等。此外,通过根据需要而将最大峰值时间差等一起显示,能够一次辨识许多诊断信息,能够实现自由度更高的超声波图像诊断。
(第2实施方式)
接着,对有关第2实施方式的超声波诊断装置进行说明。
在第1实施方式中,对于不同的两个局部部位的全部组合计算DI值,从其中使用最大的DI值生成非同步性评价支援信息并显示。相对于此,有关本实施方式的超声波诊断装置是操作者选择任意的不同的两个局部部位,计算关于该选择的两个局部部位的DI值,使用它生成非同步性评价支援信息并显示的。
图8是表示有关第2实施方式的非同步性评价支援处理的流程的流程图。在该图中,步骤S1'~步骤S4'的各处理分别与图3所示的步骤S1~步骤S4的各处理实质上是相同的。以下,对步骤S5'、步骤S6'的各处理的内容进行说明。
[关于任意的两个部位的DI值计算:步骤S5']
接着,控制单元31从操作者经由操作单元41接受任意的两个部位的指定。对于指定方法并没有特别限定。例如,可以举出对分割为各部位的壁运动参数图像或超声波图像指定希望的两个部位,或者指定对应于各部位的时间变化曲线中的、对应于希望的两个部位的两个时间变化曲线等的方法。DI值计算单元38计算关于所指定的两个部位的DI值(步骤S5')。计算方法是在第1实施方式中叙述那样的。
[DI值及指定的两个部位的标识显示:步骤S6']
显示单元23将标识了所指定的两个局部部位的至少一个的壁运动参数图像、时间变化曲线与最大DI值一起作为非同步性评价支援而以规定的形态显示(步骤S6')。
根据以上所述的超声波诊断装置,作为反映左心室的失调的指标,对希望的两个部位(局部部位),计算将心脏的不同的两个部位间的面积差分值用规定的值标准化而得到的DI值。此外,可以与关于得到的任意的两个部位的DI值一起,将标识了该两个局部部位的至少一个的壁运动参数图像、时间变化曲线与最大DI值一起作为非同步性评价支援而用规定的形态显示。
另外,本发明并不原样限定于上述实施方式,在实施阶段中在不脱离其主旨的范围内能够将构成要素变形而具体化。作为具体的变形例,例如有如下的形态。
(1)有关本实施方式的各功能也可以通过将执行该处理的程序安装到工作站等的计算机中、将它们在存储器上展开来实现。此时,能够使计算机执行该方法的程序也可以保存在磁盘(软盘(注册商标)、硬盘等)、光盘(CD-ROM、DVD等)、半导体存储器等的记录介质中而发布。
(2)在上述各实施方式中,表示了通过B模式摄像取得关于心脏组织的超声波数据的例子。但是,并不局限于该例,也可以通过组织多普勒法取得关于心脏组织的超声波数据。
此外,通过在上述实施方式中公开的多个构成要素的适当的组合,能够形成各种发明。例如,也可以从实施方式所示的全部构成要素中删除某几个构成要素。进而,也可以将跨越不同的实施方式的构成要素适当组合。
已经说明了一些实施方式,这些实施方式只是例示而并不限定本发明的范围。事实上,这里描述的新的方法和系统可以通过各种其他形式实施,进而,在不脱离本发明的主旨的范围内能够对这里描述的方法和系统进行各种省略、替代及变更。权利要求书和其等价物涵盖这些在本发明的主旨范围内的形式及变更。

Claims (20)

1.一种超声波诊断装置,具备:数据取得单元,将包括被检体的心脏的至少一部分的二维区域或三维区域作为被扫描区域进行超声波扫描,取得与该被扫描区域相关的超声波数据;以及时间变化曲线取得单元,使用与上述被扫描区域相关的超声波数据,取得与心脏的多个部位相关的壁运动参数的时间变化曲线,该超声波诊断装置的特征在于,还具备:
指标值取得单元,使用与上述心脏的多个部位中的任意两个部位相关的时间变化曲线,计算差分值,并且计算上述差分值在解析期间中的总和,从而取得用来评价心脏的壁运动的能力的指标值;以及
显示单元,将上述指标值以规定的形态显示。
2.如权利要求1所述的超声波诊断装置,其特征在于,
上述指标值取得单元通过计算解析期间的各时相中的上述任意两个部位的上述壁运动参数的差分值、并且计算上述差分值在上述解析期间中的总和并用规定的值标准化,由此取得上述指标值。
3.如权利要求2所述的超声波诊断装置,其特征在于,
上述指标值取得单元计算与能够从上述多个部位选择的两个组合的全部相关的上述指标值,并且从与上述全部的组合相关的上述多个指标值中决定代表值;
上述显示单元将上述代表值及得到该代表值的两个部位中的至少一个进行显示。
4.如权利要求3所述的超声波诊断装置,其特征在于,
上述代表值是上述指标值的最大值。
5.如权利要求2所述的超声波诊断装置,其特征在于,
上述时间变化曲线取得单元使用与上述多个部位相关的各时间变化曲线取得平均的时间变化曲线;
上述指标值取得单元将上述平均的时间变化曲线在上述解析期间内的峰值作为上述规定的值而执行上述标准化。
6.如权利要求2所述的超声波诊断装置,其特征在于,
上述指标值取得单元将与上述两个部位的一个部位相关的上述时间变化曲线在上述解析期间内的峰值作为上述规定的值而执行上述标准化。
7.如权利要求2所述的超声波诊断装置,其特征在于,
上述解析期间是一个心搏周期、收缩期、扩张期中的某个。
8.如权利要求2所述的超声波诊断装置,其特征在于,
还具备选择单元,该选择单元用于选择应变、应变率、位移、速度、扭转、扭转率、局部面积变化、局部面积变化率中的某个作为上述壁运动参数;
上述时间变化曲线取得单元取得所选择的某个上述壁运动参数的时间变化曲线。
9.如权利要求2所述的超声波诊断装置,其特征在于,
上述指标值取得单元使用能够对每个上述壁运动参数设定的规定的值,执行上述标准化。
10.如权利要求2所述的超声波诊断装置,其特征在于,
上述指标值取得单元取得具有时间的维度的单位的上述指标值而作为上述指标值。
11.一种超声波图像处理装置,具备时间变化曲线取得单元,该时间变化曲线取得单元使用将包括被检体的心脏的至少一部分的二维区域或三维区域作为被扫描区域进行超声波扫描而取得的、与上述被扫描区域相关的超声波数据,取得与心脏的多个部位相关的壁运动参数的时间变化曲线,该超声波图像处理装置的特征在于,还具备:
指标值取得单元,使用与上述心脏的多个部位中的任意两个部位相关的时间变化曲线,计算差分值,并且计算上述差分值在解析期间中的总和,从而取得用来评价心脏的壁运动的能力的指标值;以及
显示单元,将上述指标值以规定的形态显示。
12.如权利要求11所述的超声波图像处理装置,其特征在于,
上述指标值取得单元通过计算解析期间的各时相中的上述任意两个部位的上述壁运动参数的差分值、并且计算上述差分值在上述解析期间中的总和并用规定的值标准化,由此取得上述指标值。
13.如权利要求12所述的超声波图像处理装置,其特征在于,
上述指标值取得单元计算与能够从上述多个部位选择的两个组合的全部相关的上述指标值,并且从与上述全部的组合相关的上述多个指标值中决定代表值;
上述显示单元将上述代表值及得到该代表值的两个部位中的至少一个进行显示。
14.如权利要求13所述的超声波图像处理装置,其特征在于,
上述代表值是上述指标值的最大值。
15.如权利要求12所述的超声波图像处理装置,其特征在于,
上述时间变化曲线取得单元使用与上述多个部位相关的各时间变化曲线取得平均的时间变化曲线;
上述指标值取得单元将上述平均的时间变化曲线在上述解析期间内的峰值作为上述规定的值而执行上述标准化。
16.如权利要求12所述的超声波图像处理装置,其特征在于,
上述指标值取得单元将与上述两个部位的一个部位相关的上述时间变化曲线在上述解析期间内的峰值作为上述规定的值而执行上述标准化。
17.如权利要求12所述的超声波图像处理装置,其特征在于,
上述解析期间是一个心搏周期、收缩期、扩张期中的某个。
18.如权利要求12所述的超声波图像处理装置,其特征在于,
还具备选择单元,该选择单元用于选择应变、应变率、位移、速度、扭转、扭转率、局部面积变化、局部面积变化率中的某个作为上述壁运动参数;
上述时间变化曲线取得单元取得所选择的某个上述壁运动参数的时间变化曲线。
19.如权利要求12所述的超声波图像处理装置,其特征在于,
上述指标值取得单元使用能够对每个上述壁运动参数设定的规定的值,执行上述标准化。
20.如权利要求12所述的超声波图像处理装置,其特征在于,
上述指标值取得单元取得具有时间的维度的单位的上述指标值而作为上述指标值。
CN2011101726375A 2010-06-25 2011-06-24 超声波诊断装置及超声波图像处理装置 Active CN102293669B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2010145519A JP5597455B2 (ja) 2010-06-25 2010-06-25 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP145519/2010 2010-06-25

Publications (2)

Publication Number Publication Date
CN102293669A CN102293669A (zh) 2011-12-28
CN102293669B true CN102293669B (zh) 2013-12-18

Family

ID=45353189

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101726375A Active CN102293669B (zh) 2010-06-25 2011-06-24 超声波诊断装置及超声波图像处理装置

Country Status (3)

Country Link
US (1) US9161737B2 (zh)
JP (1) JP5597455B2 (zh)
CN (1) CN102293669B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130211256A1 (en) * 2010-10-26 2013-08-15 Oslo Universitetssykehus Hf Method for myocardial segment work analysis
CN104640506B (zh) * 2013-09-18 2017-06-30 东芝医疗系统株式会社 超声波诊断装置、医用图像处理装置和医用图像处理方法
KR20150086717A (ko) * 2014-01-20 2015-07-29 삼성메디슨 주식회사 의료 영상 표시 방법 및 장치
RU2017105599A (ru) * 2014-08-05 2018-09-06 Иновакор Аб Система для определения состояния сердца
JP7019287B2 (ja) * 2016-03-03 2022-02-15 キヤノンメディカルシステムズ株式会社 超音波診断装置及び画像処理プログラム
US11304681B2 (en) 2016-03-03 2022-04-19 Canon Medical Systems Corporation Ultrasonic diagnostic apparatus and image processing method
JP6748866B2 (ja) * 2016-07-26 2020-09-02 国立大学法人九州大学 画像処理のためのコンピュータプログラム、画像処理装置及び方法
JP7271126B2 (ja) * 2018-10-15 2023-05-11 キヤノンメディカルシステムズ株式会社 超音波診断装置及び医用画像処理装置
JP7236312B2 (ja) * 2019-04-04 2023-03-09 富士フイルムヘルスケア株式会社 超音波診断装置、信号処理装置、および、プログラム
CN110720948B (zh) * 2019-11-12 2021-02-02 无锡海斯凯尔医学技术有限公司 基于超声检测系统的生物体征检测方法
US12076182B1 (en) * 2019-11-25 2024-09-03 Analog Devices, Inc. Chest patch for measuring displacement of ventricular walls of subject's heart

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2952913B2 (ja) 1989-12-15 1999-09-27 ソニー株式会社 記録再生装置
JPH09313486A (ja) * 1996-06-02 1997-12-09 Noritoshi Nakabachi 超音波生体計測装置
DE69732511T2 (de) 1996-10-29 2006-01-12 Koninklijke Philips Electronics N.V. Verarbeitungsverfahren für Signale von Objekten mit sich bewegenden Teilen und Echographie-Vorrichtung dafür
EP1421905B1 (en) * 2001-08-20 2011-06-08 Japan Science and Technology Agency Tissue identifying method in ultrasonography and ultrasonograph
JP3761510B2 (ja) 2002-11-27 2006-03-29 日興製薬株式会社 花粉症を改善または予防する食品
US20040116810A1 (en) 2002-12-17 2004-06-17 Bjorn Olstad Ultrasound location of anatomical landmarks
CN100383554C (zh) 2003-04-15 2008-04-23 皇家飞利浦电子股份有限公司 心壁应变成像
JP4256234B2 (ja) 2003-08-28 2009-04-22 株式会社岡村製作所 積み重ね椅子における肘掛け装置
JP4805669B2 (ja) 2005-12-27 2011-11-02 株式会社東芝 超音波画像処理装置及び超音波画像処理装置の制御プログラム
JP4966108B2 (ja) * 2007-06-25 2012-07-04 株式会社東芝 超音波診断装置、超音波画像処理装置、及び超音波画像処理プログラム
JP5238201B2 (ja) 2007-08-10 2013-07-17 株式会社東芝 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP5454844B2 (ja) 2008-08-13 2014-03-26 株式会社東芝 超音波診断装置、超音波画像表示装置及び超音波画像表示プログラム
JP5438936B2 (ja) * 2008-08-29 2014-03-12 株式会社東芝 超音波診断装置、画像処理装置及び画像処理プログラム
US8208703B2 (en) 2008-11-05 2012-06-26 Toshiba Medical Systems Corporation Medical image analysis apparatus and image analysis control program

Also Published As

Publication number Publication date
US9161737B2 (en) 2015-10-20
JP5597455B2 (ja) 2014-10-01
US20110319761A1 (en) 2011-12-29
JP2012005708A (ja) 2012-01-12
CN102293669A (zh) 2011-12-28

Similar Documents

Publication Publication Date Title
CN102293669B (zh) 超声波诊断装置及超声波图像处理装置
CN101361665B (zh) 超声波诊断装置、超声波图像处理装置以及方法
Hashemi et al. Global time-delay estimation in ultrasound elastography
Pernot et al. Real-time assessment of myocardial contractility using shear wave imaging
JP6640922B2 (ja) 超音波診断装置及び画像処理装置
WO2014080833A1 (ja) 超音波診断装置、画像処理装置及び画像処理方法
CN101292882B (zh) 超声波图像处理装置以及超声波图像处理方法
Pollick et al. Echocardiographic and cardiac Doppler assessment of mice
US7374538B2 (en) Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing
CN101791230B (zh) 超声波诊断装置、超声波图像处理装置、医用图像诊断装置
CN104688266A (zh) 三维弹性超声成像中的运动修正
CN101647717A (zh) 超声波诊断装置、超声波图像显示装置及医用图像诊断装置
CN101563626A (zh) 用于最佳散斑跟踪的双路径处理
WO2003077765A1 (fr) Systeme d'echographie
CN107427278A (zh) 通过超声量化材料的弹性的方法
JP5944633B2 (ja) 超音波診断装置、画像処理装置及びプログラム
Salles et al. 3D myocardial mechanical wave measurements: toward in vivo 3D myocardial elasticity mapping
Ma et al. Comparison of cardiac displacement and strain imaging using ultrasound radiofrequency and envelope signals
US20070161895A1 (en) Ultrasound system and method of displaying ultrasound image
JP7142115B2 (ja) 超音波診断装置および超音波診断装置の制御方法
CN103153197A (zh) 超声波诊断装置、图像处理装置以及图像处理方法
CN102958452A (zh) 超声波诊断装置、医用图像处理装置、医用图像处理方法
Lokesh et al. Diverging beam with synthetic aperture technique for rotation elastography: preliminary experimental results
RU2356494C1 (ru) Способ диагностики множественного поражения коронарных артерий у больных ибс
CN114503214A (zh) 用于对心脏系统进行建模的方法和系统

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
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160802

Address after: Japan Tochigi

Patentee after: Toshiba Medical System Co., Ltd.

Address before: Tokyo, Japan, Japan

Patentee before: Toshiba Corp

Patentee before: Toshiba Medical System Co., Ltd.