专利名称：超声波诊断装置、医用图像处理装置以及医用图像处理程序的制作方法
技术领域：
本发明涉及一种通过超声波探头向被检查体的内部发送超声波 从而根据反射波得到被检查体的医用图像的超声波诊断装置、用于处 理该超声波诊断装置得到的医用图像的医用图像处理装置和医用图 像处理程序。特别是，本发明涉及用于评价生物组织的运动功能的超 声波诊断装置、医用图像处理装置和医用图像处理程序。
背景技术：
超声波诊断装置具有通过仅将超声波探头与体表接触的简单操 作就可以当场观察图像的优点，因此超声波诊断装置被广泛用于检查 生物组织的功能或形状。近年来，心脏壁的运动等生物组织的运动功 能的评价、特别是三维运动功能的评价已经引起人们的注意。
专利文献1公开了一种用于评价生物组织的运动功能的超声波 诊断装置。该超声波诊断装置从生物组织的二维图像中提取特征点， 并根据特征点来指定测量点。然后，设定具有以每个指定的测量点为 中心包含预定数量以上的特征点的大小的模板，使用该模板来进行两 个图像之间的互相关处理，以测量组织的位移或速度等物理参数，从 而提高测量精度。
专利文献1:日本未审查专利申请公开No. 2004 - 31329
发明内容
生物组织通常三维地移动。例如，心脏进行壁厚方向上的运动(壁 厚的变化、壁厚的变形等)和与其垂直的方向上的变形或伸缩等运动。 由于通过测量生物组织的二维速度或位移来进行通常的运动功能评价，因此在评价三维运动方面有限制。
同时，当使用生物组织的三维图像来评价三维运动时，即，当指 定三维分布的多个测量点并分析每个测量点的三维位移以评价生物 组织的三维运动时，要处理的数据量增加，从而产生处理时间变长的 问题。
另外，当使用三维图像来进行运动功能评价时，出现难以在图像 上指定运动测量对象区域的问题。例如，当需要测量心脏壁厚的变化 时，必须从心脏的三维图像中找出对应于内膜和外膜的区域并且指定 这些区域。但是，为了准确地执行这样的操作，要求熟练的手部技巧 和相当多的操作时间。
本发明为了解决上述问题而设计，其目的在于提供一种可以在短 时间内测量生物组织的三维运动的超声波诊断装置、医用图像处理装 置和图像处理程序。
另外，本发明的另 一个目的是提供一种可以容易地指定测量生物 组织的运动的区域的超声波诊断装置、医用图像处理装置和医用图像 处理程序。
为了实现上述目的，根据本发明的第一方面，提供一种超声波诊
断装置，包括超声波探头；收发器，向超声波探头发送和从超声波 探头接收超声波；图像生成部分，根据作为发送和接收超声波的结果 而得到的接收信号，针对生物组织的两个以上的断面位置中的每一个 断面位置生成按时间序列的一系列断层图像的图像数据；显示部分， 针对所述两个以上的断面位置中的每一个断面位置，根据所生成的图 像数据，显示所述一系列断层图像中的一个断层图像；指定部分，针 对所述两个以上的断面位置中的每一个断面位置，在所显示的所述一 个断层图像上指定测量图像区域；和计算器，针对所述两个以上的断 面位置中的每 一 个断面位置，根据所述 一 系列断层图像的图像数据， 计算表示所指定的测量图像区域中的生物组织的运动状态的局部运 动信息，并且根据针对所述两个以上的断面位置中的每一个断面位置 计算的局部运动信息，计算表示生物组织的运动状态的运动信息，其
ii中，显示部分显示由计算器计算的运动信息。
另外，根据本发明的第二方面，提供一种超声波诊断装置，包括 超声波探头；收发器，向超声波探头发送和从超声波探头接收超声波; 图像生成部分，根据作为发送和接收超声波的结果而得到的接收信 号，针对生物组织的一个以上的断面位置中的每一个断面位置生成按 时间序列的一系列断层图像的图像数据；显示部分，针对所述一个以 上的断面位置中的每一个断面位置，根据所生成的图像数据，显示所 述一系列断层图像中的一个断层图像；指定部分，针对所述一个以上 的断面位置中的每一个断面位置，在所显示的所述一个断层图像上指 定测量图像区域；和计算器，针对所述一个以上的断面位置中的每一 个断面位置，根据所述一系列断层图像的图像数据，计算表示所指定 的测量图像区域中的生物组织的运动状态的局部运动信息，其中，显 示部分显示由计算器计算的局部运动信息。
另外，根据本发明的第三方面，提供一种医用图像处理装置，处 理由超声波诊断装置得到的生物组织的医用图像的图像数据，所述医 用图像处理装置包括存储器，存储生物组织的两个以上的断面位置 处的、按时间序列的一系列断层图像的图像数据；显示部分，针对所 述两个以上的断面位置中的每一个断面位置，根据所存储的图像数 据，显示所述一系列断层图像中的一个断层图像；指定部分，针对所 述两个以上的断面位置中的每一个断面位置，在所显示的所述一个断 层图像上指定测量图像区域；和计算器，针对所述两个以上的断面位 置中的每一个断面位置，根据所述一系列断层图像的图像数据，计算 表示所指定的测量图像区域中的生物组织的运动状态的局部运动信 息，并且根据针对所述两个以上的断面位置中的每一个断面位置计算 的局部运动信息，计算表示生物组织的运动状态的运动信息，其中， 显示部分显示由计算器计算的运动信息。
另外，根据本发明的第四方面，提供一种医用图像处理装置，处 理由超声波诊断装置得到的生物组织的医用图像的图像数据，所述医 用图像处理装置包括存储器，存储生物组织的一个以上的断面位置中的每一个断面位置处的、按时间序列的一系列断层图像的图像数
据；显示部分，针对所述一个以上的断面位置中的每一个断面位置， 根据所存储的图像数据，显示所述一系列断层图像中的一个断层图 像；指定部分，针对所述一个以上的断面位置中的每一个断面位置， 在所显示的所述一个断层图像上指定测量图像区域；和计算器，针对 所述一个以上的断面位置中的每一个断面位置，根据所述一系列断层 图像的图像数据，计算表示所指定的测量图像区域中的生物组织的运 动状态的局部运动信息，其中，显示部分显示由计算器计算的局部运 动信息。
另外，根据本发明的第五方面，提供一种医用图像处理程序，使 具有存储器和显示部分的计算机执行以下功能，其中该存储器存储生 物组织的两个以上的断面位置中的每一个断面位置处的、按时间序列 的一系列断层图像的图像数据针对所述两个以上的断面位置中的每 一个断面位置，根据所存储的图像数据，在所述显示部分上显示所述 一系列断层图像中的一个断层图像；针对所述两个以上的断面位置中 的每一个断面位置，响应于在所显示的所述一个断层图像上指定了测 量区域，根据所述一系列断层图像的图像数据，计算表示所指定的测 量图像区域中的生物组织的运动状态的局部运动信息；根据针对所述 两个以上的断面位置中的每一个断面位置计算的局部运动信息，计算 表示生物组织的运动状态的运动信息；在显示部分上显示所计算的运 动信息。
另外，根据本发明的第六方面，提供一种医用图像处理程序，使 具有存储器和显示部分的计算机执行以下功能，其中该存储器存储生 物组织的一个以上的断面位置中的每一个断面位置处的、按时间序列 的一系列断层图像的图像数据针对所述一个以上的断面位置中的每 一个断面位置，根据所存储的图像数据，在所述显示部分上显示所述 一系列断层图像中的一个断层图像；针对所述一个以上的断面位置中 的每一个断面位置，响应于在所显示的所述一个断层图像上指定了测 量区域，根据所述一系列断层图像的图像数据，计算表示所指定的测
13量图像区域中的生物组织的运动状态的局部运动信息；在显示部分上 显示所计算的局部运动信息。
在第一、第三和第五方面的发明中，针对生物组织的两个以上的 断面位置中的每一个断面位置显示一个断层图像，响应于在所显示的 断层图像上指定了测量图像区域，计算表示测量图像区域中的生物组 织的运动状态的局部运动信息，根据所计算的两个以上的断面位置中 的每一个断面位置处的局部运动信息，计算生物组织的运动信息，并 且显示运动信息。
根据本发明，可以通过根据生物组织的两个以上的断面位置中的 每一个断面位置处的局部运动信息来得到运动信息，从而测量生物组 织的三维运动。另外，由于仅考虑到两个以上的断面位置处的局部运 动信息来得到运动信息，因此不必计算断面位置之间的部分处的位 移，从而可以在短时间内得到运动信息。
而且，由于在所显示的断层图像上指定测量图像区域，因此可以 容易地指定用于得到运动信息的测量图像区域。
在第二、第四和第六方面的发明中，针对生物组织的一个以上的
断面位置中的每一个断面位置显示一个断层图像，响应于在所显示的 断层图像上指定了测量图像区域，计算表示测量图像区域中的生物组 织的运动状态的局部运动信息，并且显示局部运动信息。
根据本发明，由于在所显示的断层图像上指定测量图像区域，因 此可以容易地指定用于得到局部运动信息的测量图像区域。
图1是示出本发明的优选实施方式的超声波诊断装置的整体结 构的例子的示意性框图。
图2是示出本发明的优选实施方式的超声波诊断装置的超声波 扫描形态的例子的示意性说明图。
图3是示出本发明的优选实施方式的超声波诊断装置的超声波 扫描形态的例子的示意性说明图。图4是示出本发明的优选实施方式的超声波诊断装置的操作状 况的例子的流程图。
图5是示出本发明的优选实施方式的超声波诊断装置的显示画 面的显示状况的例子的示意图。
图6是示出本发明的优选实施方式的超声波诊断装置的显示画 面的显示状况的例子的示意图。
图7是示出本发明的优选实施方式的超声波诊断装置的处理的 例子的流程图。
图8是示出本发明的优选实施方式的超声波诊断装置的处理的 示意性说明图。
图9是示出本发明的优选实施方式的超声波诊断装置的处理的 例子的流程图。
图IO是示出本发明的优选实施方式的超声波诊断装置的处理的 示意性说明图。
图11是示出本发明的优选实施方式的超声波诊断装置的处理的 例子的流程图。
图12是示出本发明的优选实施方式的超声波诊断装置的处理的 示意性说明图。
图13是示出本发明的优选实施方式的超声波诊断装置的处理的 例子的流程图。
图14是示出本发明的优选实施方式的超声波诊断装置的处理的 示意性说明图。
图15是示出本发明的优选实施方式的超声波诊断装置的变形的 整体结构的例子的示意性框图。
图16是示出本发明的优选实施方式的超声波诊断装置的变形的 处理的示意性说明图。
图17是示出本发明的优选实施方式的超声波诊断装置的整体结 构的例子的示意性框图。
附图标记说明1:超声波诊断装置
2: 二维超声波探头
3:收发器
4:信号处理器
41: B模式处理部分
5:图像处理器
51:体数据生成部分
52: MPR处理部分
6:计算处理器
61:位移计算部分
62:运动信息计算部分
7: 存储器
71、 72:医用图像处理程序
8:用户接口
81:显示部分
82:操作部分
9:控制器
P: 二维扫描面
R:三维扫描区域
X:主扫描方向
Y:副扫描方向
ml~m6:内膜位置图像
M1~M6:外膜位置图像
1000:医用图像处理装置
具体实施例方式
以下参照附图详细说明本发明的优选实施方式的超声波诊断装 置、医用图像处理装置和医用图像处理程序。 <第一实施方式>说明本发明实施方式的超声波诊断装置。图l示出本发明的超声 波诊断装置的整体结构的例子。图l所示的超声波诊断装置l是用于 取得例如示出心脏等生物组织的形状的图像或示出血流状态的图像
的装置，包括二维超声波探头2、收发器3、信号处理器4、图像处理 器5、计算处理器6、存储器7、用户接口8和控制器9。以下说明构 成超声波诊断装置1的各个部分的详细例子。 {存储器、用户接口和控制器}
首先说明存储器7、用户接口8和控制器9。存储器7例如由硬 盘驱动器等存储设备构成。在存储器7中，预先存储用于使超声波诊 断装置l执行本发明的特征性操作的医用图像处理程序71。另外，在 存储器7中，存储超声波图像的图像数据和该图像数据附带的附带信 息(DICOM ( Digital Imaging and Communications in Medicine, 医 用数字成像和通信))等各种数据。
用户接口 8上设置有显示部分81和操作部分82。显示部分81 对应于本发明的"显示器"的例子，由液晶显示器或CRT(阴极射线管) 显示器等任意的显示设备构成。显示部分81显示由超声波诊断装置1 取得的超声波图像等图像和DICOM附带信息等信息等。
操作部分82由鼠标、轨迹球、操纵杆、控制面板和键盘等任意 的操作设备或输入设备构成。
特别是，操作部分82作为用于在显示部分81上显示的超声波图 像(断层图像)上指定测量图像区域的"指定部分"。测量图像区域是 用作测量生物组织的运动状态的基准的、断层图像上的区域(实际上 是对应于该区域的图像数据)。超声波诊断装置l进行操作，以便通 过利用超声波三维地扫描生物组织并重复三维扫描操作来生成按时 间序列的一 系列体数据，并且根据该一 系列体数据中的一个来生成生 物组织的断层图像的图像数据，详细内容将后述。操作部分82用于 在断层图像上指定测量图像区域。超声波诊断装置l进行操作，以便 通过分析所指定的测量图像区域如何按时间序列变化来测量生物组 织的运动状态。
17控制器9包括CPU等微处理器，并且根据医用图像处理程序71 控制超声波诊断装置1的各个部分。特别是，控制器9执行用于在显 示部分81上显示图像和画面的处理。另外，控制器9执行用于使超 声波诊断装置1响应于从操作部分82输出的操作信号来执行操作的 处理。
{二维超声波探头}
以往，二维超声波探头2 (可以简称为超声波探头2)具有多个 二维(例如矩阵状(格子状))设置的超声波换能器(未示出)。这 多个超声波换能器由后述的收发器3个别地驱动。
图2和图3示出由二维超声波探头2进行的超声波扫描形态。如 图2A所示，超声波探头2通过在主扫描方向X上利用从超声波换能 器的排列面输出的超声波(束)进行扫描，形成放射状(扇形形状) 的二维扫描面P。另外，如图2B所示，超声波探头2通过在垂直于 主扫描方向X的副扫描方向Y上利用超声波进行扫描，顺序形成在 副扫描方向Y上排列的扇形形状的二维扫描面Pl、 P2、 ...Pn。从而， 如图3所示，形成三维扫描区域R。
副扫描方向Y对应于本发明的"预定方向"，主扫描方向X对应 于"与预定方向垂直的方向"。二维扫描面PI ~ Pn形成在"沿预定方向 的多个(n个)位置"上。
{收发器}
收发器3包括通过向超声波探头2提供电信号来产生超声波的发 送器和用于接收从接收到超声波的反射波的超声波探头2输出的回声 信号(接收信号)的接收器(未图示)。
收发器3的发送器包括未图示的时钟产生电路、发送延迟电路和 脉冲发生器电路。时钟产生电路是用于产生决定超声波的发送定时和 发送频率的时钟信号的电路。发送延迟电路是通过在超声波发送时施 加延迟来执行发送聚焦的电路。脉冲发生器电路将对应数量的脉冲发 生器内置到各个超声波换能器的各个路径(信道)上并且进行操作， 以便在延迟后的发送定时产生驱动脉冲并提供给超声波探头2的各个
18超声波换能器。
另外，收发器3的接收器包括未示出的前置放大器电路、A/D转 换电路和接收延迟/加法电路。前置放大器电路针对每个接收信道放大 从超声波探头2的每个超声波换能器输出的回声信号。A/D转换电路 执行放大后的回声信号的A (模拟)/D (数字)转换。接收延迟/加法 电路向经过A/D转换的回声信号提供用于决定接收指向性的延迟时 间，并使其相加。通过加法处理，来自与接收指向性一致的方向的反 射分量被强调。可将进行了加法处理的信号称为"RF数据(或原始数 据)"等。收发器3将所取得的RF信号输入信号处理器4。
{信号处理器}
信号处理器4根据从收发器3输入的RF数据执行用于使回声信 号的振幅信号可视化的信号处理。由信号处理器4生成的数据被送往 控制器9，以在用户接口 8的显示部分81上显示，或者被输入图像处 理器5。信号处理器4主要包括B模式处理部分41、多普勒处理部分 42和CMF处理部分43。 (B模式处理部分)
B (亮度)模式处理部分41根据RF数据生成B模式超声波光 栅数据。更具体地说，B模式处理部分41对RF数据执行带通滤波处 理并检测输出信号的包络，对所检测的数据执行利用对数变换的压缩 处理。从而，针对二维扫描面Pl Pn中的每一个生成将信号强度表 现为辉度的亮度的断层图像的图像数据。B模式处理部分41对应于 本发明的"第一断层图像生成部分"的例子。 (多普勒处理部分)
多普勒处理部分42例如通过脉冲多普勒方法(PW多普勒方法) 或连续波多普勒方法(CW多普勒方法)来生成生物组织中的血流信 每
在脉冲多普勒方法中，可以使用脉冲波来检测由特定深度(离超
声波探头2的距离)处的血流的多普勒效应引起的超声波频率位移(多 普勒位移频率分量)。这样，由于脉冲多普勒方法具有好的距离分辨率，因此脉沖多普勒方法适合用于血流或特定部位的组织的深度测
量。在实施脉冲多普勒方法时，多普勒处理部分42通过对具有指定 大小的血流观察区域中的信号进行相位检测，从自收发器3输入的 RF数据中提取多普勒位移频率分量，并且执行FFT (快速傅立叶变 换)处理，来生成示出多普勒频率分布的数据，该多普勒频率分布表 示血流观察区域中的血流速度。
另外，与脉冲多普勒方法不同，连续波多普勒方法使用连续波来 得到叠加了超声波发送/接收方向(具有图2A所示的扇形形状的二维 扫描面P上的直径方向)上的所有部位的多普勒位移频率分量的信号， 即，示出超声波路径上的所有血流状态的信号。连续波多普勒方法具 有测量速度优良的优点。在实施连续波多普勒方法时，多普勒处理部 分42通过对在血流观察的采样线上接收的信号进行相位检测，从自 收发器3输入的RF数据中提取多普勒位移频率分量，并且执行FFT 处理，从而生成示出多普勒频率分布的数据，该多普勒频率分布表示 釆样线上的血流速度。 (CFM处理器)
CFM ( Color Flow Mapping,彩色血流成像)处理器43在执行 彩色血流成像方法时进行操作，彩色血流成像方法用彩色将生物组织 的血流信息叠加到单色B模式图像上，并实时地进行显示。所显示的 血流信息的例子包括血流的速度、分散、能力等。血流信息作为二进 制信息而得到。更具体地说，CFM处理器43包括相位检测电路、 MTI ( Moving Target Indication,运动目标指示)滤波器、自相关器、 流速/分散运算器等。CFM处理器43执行高通滤波处理(MTI滤波 处理)以得到表示生物组织的形状的形态信号和表示血流的血流信 号，从而通过自相关处理在多个位置处得到血流信息，例如血流速度、 血流分散和血流能力。另外，为了减少形态信号，可以执行非线性处 理等。
{图像处理器}
图像处理器5根据信号处理器4生成的数据来执行各种图像处理。例如，图像处理器5包括DSC (Digital Scan Converter,数字扫 描转换器)并执行将与超声波扫描同步的、由信号处理器4生成的数 据转换成用于显示的数据(电视扫描方式的数据)的处理，即扫描转 换处理。
另外，图像处理器5包括后述的体数据生成部分51和MPR处 理部分52。
(体数据生成部分)
体数据生成部分51对信号处理器4的B模式处理部分41生成 的二维扫描面Pl~Pn中的每一个的图像数据执行插值处理，生成体 数据(体素数据)。体数据生成部分51对应于"体数据生成部分"的 例子，例如包括DSC或微处理器。
当显示了基于体数据的伪三维图像时，图像处理器5对体数据执 行体呈现(rendering )处理、MIP ( Maximum Intensity Projection, 最大强度投影)处理等。 (MPR处理部分)
MPR ( Multi Plannar Reconstruction,多平面重建)处理器52 根据体数据生成部分51生成的体数据执行断面转换处理，以生成任 意断面中的断层图像的图像数据。MPR处理部分52对应于"第二断 层图像生成部分"的例子，例如包括DSC、微处理器等。另外，MPR 处理部分52、体数据生成部分51和信号处理器4的B模式处理部分 41作为本发明的"图像生成部分"的例子。
{计算处理器}
计算处理器6被用于根据图像处理器5的MPR处理部分52生 成的断层图像的图像数据来计算表示生物组织的局部运动状态的局
部运动信息或表示更大范围的运动状态的运动信息，并且作为本发明 的"计算器"的例子。
例如，当生物组织是心脏时，由计算处理器6得到的局部运动信 息的例子包括心脏壁厚的变化、变化的速度、心脏壁的运动变形、变 形率、心脏壁的内膜或外膜的旋转角度、旋转角度的速度(旋转速度)、
21内膜和外膜的相对旋转角度等(详情后述)。
另外，运动信息的例子例如包括心脏壁的扭转运动、该扭转运动 的速度、伸/缩(缩短)、伸/缩的速度、心脏壁运动的变形、变形率、 相对旋转梯度等(详情后述)。
计算处理器6包括CPU等微处理器。计算处理器6设置有位移 计算部分61和运动信息计算部分62。 (位移计算部分)
位移计算部分61按照时间序列跟踪由操作部分82在断层图像上 指定的测量图像区域，以便计算该测量图像区域(中的生物组织)的 按时间序列的位移。测量图像区域的位移对应于本发明的"局部运动 信息"的例子。
可以通过用图像之间的时间间隔(帧间隔)来除按时间序列的二 维或三维图像之间的位移，得到位移速度。相反，可以通过用图像之 间的时间间隔来乘测量图像区域的位移速度，得到图像之间的位移。 即，当给出图像之间的时间间隔时，可以认为位移与速度彼此同义。 在该意义上，可以认为在本发明中位移和速度彼此同一。
以下更详细地说明位移计算部分61的操作。如上所述，超声波 诊断装置1生成生物组织的按时间序列的一系列体数据，并且根据一 个体数据(特定时间(特定相位)的体数据)生成断层图像的图像数 据。然后，在断层图像上指定测量图像区域。位移计算部分61按时 间序列追踪测量图像区域在不同时间(时间相位)的体数据中位移多 少。这样，可以将按时间序列追踪体数据中的测量图像区域的位移的 处理称为"三维跟踪"。
另外，位移计算部分61能够在与断层图像相同的断面位置处的 不同时间(时间相位)的断层图像中按时间序列跟踪在断层图像上指 定的测量图像区域位移多少。可以将这样的追踪处理称为"二维跟 踪"。二维跟踪以下述方式执行根据按时间序列的体数据生成断面 位置处的按时间序列的断层图像的图像数据，并在按时间序列的断层 图像的图像数据中跟踪位移。另外，可以通过跟踪按时间序列的体数
22据的断面位置处的位移来进行二维跟踪。
这样的跟踪处理可以利用已知方法来进行。例如，与上述专利文
献l所公开的方法相同，可以如下进行二维跟踪从在按时间序列的 多个断层图像上指定的各个测量图像区域提取特征点，并且根据特征 点来指定测量点。然后，根据测量点设定具有包括预定数量以上的特 征点的大小的模板，使用该模板进行两个断层图像(测量图像区域) 的相关处理(图案匹配处理)，从而进行每个测量点的位移的计算。
另外，在三维跟踪中，同样地根据体数据来设定三维模板，使用 该三维模板来进行对不同时间(时间相位)的两个体数据的图案匹配 处理，以计算每个测量点的位移。
另外，可以仅在测量图像区域(例如后述的图5所示的内膜位置 图像ml)中指定特征点和测量点，或者可以在测量图像区域以外的 区域中指定，例如测量图像区域的边界之间的区域(例如与由内膜位 置图像ml和外膜位置图像Ml包围的心脏壁的断面对应的图像区 域)、测量图像区域附近的区域等。在任何情况下，都可以应用任意 跟踪方法，从而可以计算所指定的测量图像区域的按时间序列的位 移。
通过二维跟踪或三维跟踪得到的测量点的位移可以直接作为测 量图像区域的位移。另外，可以根据测量点的位移来计算测量图像区 域的边界(例如后述的图6所示的内膜位置图像ml ~ m3或外膜位置 图像Ml ~ M3 )的位移，并且将边界的位移用作测量图像区域的位移。
这样，位移计算部分61在一个断层图像中所指定的测量图像区 域中指定多个测量点，针对每帧的断层图像分别得到测量点的位置。 然后，位移计算部分根据每帧的断层图像的测量点的位置，计算每个 测量点的按时间序列的位移，并且根据每个测量点的位移计算所指定 的测量图像区域的位移。
(运动信息计算部分)
运动信息计算部分6 2执行根据位移计算部分61计算的测量图像 区域的位移来计算表示生物组织的(更大范围的)运动状态的运动信
23息的处理。运动信息计算部分62计算运动信息的处理的详细例子将 后述。
参照图4~14说明本实施方式的具有上述结构的超声波诊断装 置l的操作状况的例子。这里说明评价心脏运动状态的情况。以下参 照图4~6说明超声波图像取得操作和测量图像区域指定操作，然后 参照图7 ~ 14说明测量图像区域跟踪处理和运动信息计算处理。
{超声波图像的取得和测量图像区域的指定}
首先取得作为评价运动状态的对象的心脏的超声波图像。为了进 行该取得，当在超声波探头被放置在被检查人的心脏附近(一般在心 脏的心尖部的附近)的体表上的状态下执行预定操作时，收发器3在 控制器9的控制下控制超声波探头2,以便执行心脏的三维超声波扫 描(图2和3所示的超声波扫描)(S01)。
超声波诊断装置l重复地执行三维超声波扫描。此时，最好将三 维超声波扫描重复执行不少于一个心跳周期(心跳的一个周期，例如 在心电图中从当前的R波到下一个R波)的时间。
接着，收发器3将从超声波探头2依次输出的回声信号转换成 RF数据，并依次输入信号处理器4。信号处理器4的B模式处理部 分41根据RF数据依次生成每个二维扫描面PI ~ Pn中的断层图像的 图像数据(S02)。将生成的图像数据输入图像处理器5。当事先得到 了体数据时，可以通过对体数据执行MPR处理等来生成断层图像的
接着，图像处理器5的体数据生成部分51依次对二维扫描面 Pl~Pn中的断层图像的图像数据执行插值处理，依次生成对应于一 个三维超声波扫描的三维扫描区域R中的体数据(S03)。这样得到 按时间序列的多个体数据。
然后，MPR处理部分52根据所生成的多个体数据中的一个生成 心脏的预定断面位置处的断层图像的图像数据。在本实施方式中，生 成心尖四腔像(可以称为心尖四腔断层)的图像数据和心尖二腔像(可以称为心尖二腔断层)的图像数据(S04)。这里，心尖四腔像和心 尖二腔像都对应于沿心脏长轴方向的断面位置处的断面图像，并且它 们的断面位置彼此正交。
控制器9控制显示部分81以显示基于在步骤S04中生成的图像 数据的断层图像(心尖四腔像和心尖二腔像)(S05)。图5示出此 时的显示状况的例子。在图5中，显示部分81 (其显示画面)设置有 基于在步骤S04中生成的图像数据的断面位置指定图像显示部分81A 和81B。在图5中，在断面位置指定图像显示部分81A上显示心尖四 腔像，在断面位置指定图像显示部分81B上显示心尖二腔像。
另外，在利用心电同步的超声波诊断的情况下，在显示部分81 的心电图显示部分81F上显示心电图。在心电图显示部分81F上显示 表示取得断面位置指定图像显示部分81A和81B上显示的断层图像的 时间(时间相位和时刻)的时间光标T。在图5中，时间光标T位于 心电图的R波的时间相位。这里，可以构成为时间光标T可以沿心电 图的时间方向(水平方向)移动(例如拖放)，从体数据生成时间光 标T的目的地时间(时间相位)的断层图〗象，并显示在断面位置指定 图^象显示部分81A和81B上。
断面位置指定光标C1、 C2和C3设置在断面位置指定图像显示 部分81B的水平位置上。例如，当观察心尖四腔像或心尖二腔像时， 用户操作操作部分82的鼠标，沿垂直方向(心脏的长轴方向)拖放 断面位置指定光标C1 C3，从而指定断面位置(S06)。
控制器9例如根据心尖四腔像(和/或心尖二腔像)的显示画面 中的坐标和断面位置指定光标C1 C3的显示画面中的坐标，决定心 尖四腔像的图像数据上的、指定的断面位置的坐标，并将其发送给图 像处理器5。
MPR处理部分52根据从控制器9发送的坐标信息和在步骤S03 中生成的多个体数据，在每个指定的断面位置上生成在心脏的短轴方 向上具有断面的断层图像的图像数据(S07)。此时，MPR处理部分 52针对多个体数据中的每一个，生成每个指定的断面位置上的断层图像的图像数据。
控制器9控制显示部分81的断层图像显示部分81C、81D和81E， 以显示基于从与在步骤S04中生成的体数据相同的体数据生成的图像 数据的断层图像(S08)。在步骤S07中，首先生成在步骤S08的显 示处理中涉及的图像数据，其它图像数据的生成处理可以在步骤S08 的处理之后在后台进行。
在图5所示的例子中，将断面位置指定光标C1 C3分别放置到 心脏的心尖部水平、乳头肌水平和心基部水平。在这种情况下，利用 断面位置指定光标Cl指定的心尖部水平的断层图像(心尖部短轴图 像Gl)显示在断层图像显示部分81C上。利用断面位置指定光标C2 指定的乳头肌水平的断层图像(乳头肌短轴图像G2)显示在断层图 像显示部分81D上。利用断面位置指定光标C3指定的心基部水平的 断层图像(心基部短轴图像G3)显示在断层图像显示部分81E上。
用户操作操作部分82,在断层图像显示部分81C、 81D和81E 中的每一个上显示的断层图像上指定测量图像区域(S09)。该操作 例如可以通过拖动鼠标向断层图像输入表示测量图像区域的边界来 进行。
图6示出在图5所示的心尖部短轴图像G1、乳头肌短轴图像G2 和心基部短轴图像G3上指定测量图像区域的状况的例子。在各短轴 图像G1 G3的内侧示出的虚线表示内膜，在外侧示出的虚线表示外 膜。在观察所显示的短轴图像G1~G3时，用户操作操作部分82以 描绘内膜和外膜，从而输入表示测量图像区域的线。
这样，如图6所示，在心尖部短轴图像G1中，输入通过描绘心 脏壁的内膜而得到的内膜位置图像ml和通过描绘外膜而得到的外膜 位置图像M1。在乳头肌短轴图像G2中，输入通过描绘心脏壁的内 膜而得到的内膜位置图像m2和通过描绘外膜而得到的外膜位置图像 M2。在心基部短轴图像G3中，输入通过描绘心脏壁的内膜而得到的 内膜位置图像m3和通过描绘外膜而得到的外膜位置图像M3。
从而结束超声波图像取得操作和测量图像区域指定操作。
26{测量图像区域跟踪处理和运动信息计算处理}
接着，参照图7 ~ 14说明测量图像区域跟踪处理(位移计算处理) 和运动信息计算处理。以下针对要取得的每个运动信息来进行说明。 以下说明的每个处理在图4所示的流程图的步骤S09之后执行。以下 说明的多个处理可以依次执行，也可以并行执行。而且，不必个别地 执行在不同处理中执行的相同过程。 (运动信息扭转运动)
首先参照图7和图8说明取得心脏壁的扭转运动状态作为运动信 息的处理。当在断层图像显示部分81C、 81D和81E中的每一个上显 示的断层图像上指定了测量图像区域时(步骤S09)，位移计算部分 61根据针对按时间序列的多个体数据中的每一个生成的断层图像的 图像数据(S07)，执行内膜位置图像ml的二维跟踪，从而计算内膜 位置图像ml绕着与心尖部短轴图像Gl的断面正交的方向(心脏的 长轴方向)上的轴的旋转角度(局部运动信息)。同样，在乳头肌短 轴图像G2和心基部短轴图像G3中，位移计算部分计算内膜位置图 像m2和m3绕着心脏的长轴方向上的轴的旋转角度(局部运动信息) (S11)。代替内膜位置图像ml、 m2和m3,也可以计算外膜位置图 ^f象Ml、 M2、 M3的S走转角度。
此时，位移计算部分61例如针对每个时间相位计算内膜位置图 像ml、 m2、 m3的旋转角度，作为在步骤S09中输入内膜位置图像 ml等的时间相位(基准时间相位)的旋转角度。位移计算部分可以 依次计算按时间序列相邻的帧(即连续帧)中的内膜位置图像ml、 m2和m3的旋转角度。
运动信息计算部分62计算内膜位置图像ml的旋转角度与内膜 位置图像m2的旋转角度之差(相对旋转角度)(S12)。同样，运动 信息计算部分计算内膜位置图像m2的旋转角度与内膜位置图像m3 的旋转角度之差(相对旋转角度)(S13)。相对旋转角度对应于本 发明的"差信息"的例子。
参照图8详细说明步骤S12和S13中的处理。在断层图像显示部分81C、 81D和81E中，例如将逆时4十方向定义为正S走转方向(+
e方向)。另外，内膜位置图像mi的旋转角度用ei指示，内膜位置
图像m2的旋转角度用62指示，内膜位置图像m3的旋转角度用03 指示。
此时，利用A612 = 91-02 (或02-ei)得到在步骤S12中计算的
— L 丄丄 rfr * 八j 一.， rTV ia一為 八一八> / ^u、 八a a 、 乂曰f，， ，t_ tTfftf
7日^]"》疋将用厌AW1Z。
甲J乂干J = O乂-"J 、-人0J-WZ J S丄^
中计算的相对旋转角度A023。
在步骤S12中得到的相对旋转角度A012是表示心尖部短轴图像 Gl与乳头肌短轴图像G2的断面位置之间的心脏壁的扭转运动的状
态(大小)的信息。即，当相对旋转角度Aei2 = o (ei-e2)时，可 以认为，在断面位置之间的任意位置上，心脏壁在相同方向上旋转相 同角度，并且在旋转方向上没有扭转。
而在IAei2i邦的情况下，在断面位置之间旋转角度存在差异，并 且心脏壁在旋转角度方向上扭转。随着相对旋转角度Aei2的绝对值 变大，心脏壁的扭转变大。例如，当ei和e2的符号彼此不同时，即，
当内膜位置图像ml的旋转方向与内膜位置图像m2的旋转方向彼此 相反时，相对旋转角度A612的绝对值相对较大。
同样，在步骤S13中得到的相对旋转角度A623是表示乳头肌短 轴图像G2和心基部短轴图像G3的断面位置之间的心脏壁的扭转运 动的大小的信息。
控制器9控制显示部分81，以显示在步骤S12和S13中计算的 相对i走转角度和A023，作为表示心脏壁的扭转运动的大小的运 动信息(S14)。用户可以通过参照所显示的相对旋转角度A912和 A023，来把握心脏壁的扭转运动的大小。这里，还可以分别计算心脏 壁的内膜和外膜的相对旋转角度，并且根据两个相对旋转角度评价扭 转运动的大小(例如通过得到两个相对旋转角度的平均)。
通过利用时间对相对旋转角度A612进行微分，可以得到内膜位 置图像ml和m2之间的心脏壁的扭转运动的速度。同样，通过利用 时间对相对旋转角度A023进行微分，可以得到内膜位置图像m2和
28m3之间的心脏壁的扭转运动的速度。然后，可以构成为在显示部分 81上显示速度。这里，"微分"包括用得到相对旋转角度的帧之间的时 间间隔来除相对旋转角度的处理以及通常的微分计算。 (运动信息相对旋转梯度)
参照图9和图IO说明取得心脏壁的相对旋转梯度作为运动信息
,.t lA im丄rt -丄iA丄丄l乂 r& 曰 士一— 、、 riA Pi J"A 、二 -丄-ter rit ,— i /一 A
wv^c王。々曰"力疋ww厌疋尽w"胆^:wv件^"w^渭揮厌wv'i各,&。
首先，与图7所示的步骤S11同样，位移计算部分61计算心尖 部短轴图像G1的内膜位置图像ml的旋转角度ei、乳头肌短轴图像 G2的内膜位置图像m2的旋转角度e2、心基部短轴图像G3的内膜位 置闺像m3的旋转角度03 (S21)。
接着，与步骤S12和S13同样，运动信息计算部分62计算内膜 位置图像ml的旋转角度ei与内膜位置图像m2的旋转角度02的相 对旋转角度A012 (S22)，并且计算内膜位置图像m2的旋转角度e2 与内膜位置图像m3的旋转角度G3的相对旋转角度A023 ( S23 )。
接着，运动信息计算部分62计算心尖部短轴图像Gl和乳头肌 短轴图像G2之间的距离d12 (S24)，并且计算乳头肌短轴图像G2 和心基部短轴图像G3之间的距离d23(S25)。距离d12和d23例如 可以根据在步骤S06之后由控制器9得到的心尖部短轴图像Gl、乳 头肌短轴图像G2和心基部短轴图像G3的断面位置的坐标来计算。
另夕卜，如图IO所示，运动信息计算部分62利用在步骤S24中得 到的距离d12来除在步骤S22中得到的相对旋转角度Aei2，以计算 内膜位置图像ml和内膜位置图像m2之间的相对旋转梯度(8612-A012/dl2) (S26)。同样，运动信息计算部分62利用在步骤S25中 得到的距离d23来除在步骤S23中得到的相对旋转角度A023,以计 算内膜位置图像m2和内膜位置图像m3之间的相对旋转梯度( =A023/d23) (S27)。
控制器9控制显示部分81，以显示在步骤S26和S27中计算的 相对旋转梯度8012和8623，作为表示心脏壁的扭转运动的程度的运 动信息(S28)。相对旋转梯度3912表示心尖部水平和乳头肌水平的内膜之间的 每单位距离的扭转的大小。相对旋转梯度3623表示乳头肌水平和心 基部水平的内膜之间的每单位距离的扭转的大小。即，相对旋转梯度
和洲23是表示心脏壁(内膜)的扭转程度的运动信息。用户可 以通过参照所显示的相对旋转梯度8912和8023来掌握心脏壁的扭转 运动的程度。还可以针对心脏壁的内膜和外膜来计算相对旋转梯度， 并根据两个相对旋转梯度来评价扭转运动的程度(例如从两个相对旋 转梯度得到平均)。
(运动信息长轴方向的伸缩)
参照图11和图12说明取得心脏壁的长轴方向的伸缩(缩短)作 为运动信息的处理。
首先，位移计算部分61根据按时间序列的多个体数据(S03)分 别对心尖部短轴图像G1的内膜位置图像ml、乳头肌短轴图像G2的 内膜位置图像m2和心基部短轴图像G3的内膜位置图像m3执行三 维跟踪，从而分别针对内膜位置图像ml、内膜位置图像m2、内膜位 置图像m3计算指定了测量图像区域的断层图像的三维位移(Axl、 Ayl、 Azl) 、 (Ax2、 Ay2、 Az2 ) 、 ( Ax3、 Ay3、 Az3 ) ( S31)。位 移对应于"局部运动信息"的例子。位移计算部分可以计算外膜位置图 像M1、 M2和M3的三维位移，来代替内膜位置图像ml、 m2和m3。
此时，例如针对每个时间相位，位移计算部分61计算内膜位置 图像ml、 m2和m3的三维位移，作为相对于在步骤S09中输入内膜 位置图像ml等的基准时间相位的三维位移。位移计算部分可以依次 计算连续帧中的内膜位置图像ml、 m2和m3的三维位移。
另外，Ax和Ay分别指示图2和图3所示的X方向(方向的一 侧为+X方向)和Y方向上的位移。包括X方向和Y方向的平面平 行于心尖部短轴图像Gl、乳头肌短轴图像G2和心基部短轴图像G3 的断面。另外，Az指示与X方向和Y方向正交的Z方向(例如当从 乳头肌水平看时，心尖部方向是-Z方向，心基部方向是+Z方向) 上的位移。Z方向平行于心脏的长轴方向。另外，位移计算部分61从三维位移(Axl、 Ayl、 Azl) 、 (Ax2、Ay2、 Az2) 、 (Ax3、 Ay3、 Az3)中提取Z方向(长轴方向)上的位移Azl、 Az2和Az3 (S32)。
这里，假定心尖部短轴图像G1等的断面平行于XY平面。但是，即使当心尖部短轴图像Gl等的断面不平行于XY平面时，也可以通过在Z方向上投影三维位移(向量)而容易地计算Z方向上的位移。
接着，如图12所示，运动信息计算部分62计算内膜位置图像ml在Z方向上的位移Azl与内膜位置图像m2在Z方向上的位移Az2之差，以得到心尖部水平与乳头肌水平之间的心脏壁的伸缩Azl2-Azl - Az2 (或Az2 - Azl) ( S33 )。同样，运动信息计算部分62计算内膜位置图像m2在Z方向上的位移Az2与内膜位置图像m3在Z方向上的位移Az3之差，以得到乳头肌水平与心基部水平之间的心脏壁的伸缩Az23 = Az2 - Az3 (或Az3 - Az2 ) ( S34 )。伸缩Azl2和Az23对应于"差信息"的例子。
控制器9控制显示部分81，以显示在步骤S33和S34中计算的心脏壁的伸缩Azl2和Az23，作为表示心脏壁的伸缩的大小的运动信息(S35)。用户可以通过参照所显示的心脏壁的伸缩Azl2和Az23，来把握心脏壁的伸缩的大小。
还可以针对内膜和外膜中的每一个分别计算心脏壁的伸缩，并根据两个伸缩的值来评价伸缩的大小(例如从两个伸缩值得到平均)。
另外，通过利用时间对伸缩Azl2进行微分，可以得到内膜位置图像ml和m2之间的心脏壁的伸缩运动的速度。同样，通过利用时间对相对旋转角度Az23进行微分，可以得到内膜位置图像m2和m3之间的心脏壁的伸缩运动的速度。然后，可以构成为在显示部分81上显示速度。这里，"微分"具有与上述相同的含义。(运动信息长轴方向的变形)
参照图13和图14说明取得心脏壁的长轴方向的变形作为运动信息的处理。变形是表示心脏壁的变形大小的程度的信息，表示心脏壁的变形状态。
31首先，与取得长轴方向的收缩的情况同样，位移计算部分61分别针对内膜位置图像ml、内膜位置图像m2、内膜位置图像m3计算在其上指定了测量图像区域的断层图像的三维位移(Axl、 Ay 1 、 Azl)、(Ax2、 Ay2、 Az2) 、 (Ax3、 Ay3、 Az3) ( S41)，并且从三维位移中提取Z方向(长轴方向)上的位移Azl、 Az2和Az3 ( S42 )。
接着，以与上述同样的方式，运动信息计算部分62计算心尖部水平与乳头肌水平之间的心脏壁的伸缩Azl2 = Azl-Az2 (S43)，并且计算乳头肌水平与心基部水平之间的心脏壁的伸缩Az23 = Az2 -Az3 ( S44 )。
与图9所示的步骤S24和S25同样，针对指定了测量图像区域的心尖部短轴图像Gl、乳头肌短轴图像G2和心基部短轴图像G3，运动信息计算部分62计算心尖部短轴图像Gl和乳头肌短轴图像G2之间的距离d12 (S45),并且计算乳头肌短轴图像G2和心基部短轴图像G3之间的距离d23 (S46)。
而且，运动信息计算部分62利用在步骤S45中计算的距离d12来除在步骤S43中计算的伸缩Azl2，以计算心尖部水平和乳头肌水平之间的长轴方向的变形(3zl2 = Azl2/dl2 ) (S47)。另外，运动信息计算部分62利用在步骤S46中计算的距离d23来除在步骤S44中计算的伸缩Az23，以计算乳头肌水平和心基部水平之间的长轴方向的变形(Sz23 = Az23/d23 ) ( S48 )。
控制器9控制显示部分81，以显示在步骤S47和S48中计算的心脏壁的变形Szl2和6z23，作为表示心脏壁的变形的大小的运动信息(S49)。用户可以通过参照所显示的心脏壁的变形Szl2和Sz23，来把握心脏壁的变形的大小。
还可以分别计算心脏壁的内膜和心脏壁的外膜的变形，并根据两个变形值来评价变形的大小(例如从两个变形值得到平均)。(运动信息长轴方向的变形率)
说明取得心脏壁的长轴方向的变形率作为运动信息的处理。变形率是表示心脏壁的变形的时间变化率的信息，表示心脏壁的变形状
32态。
在取得变形率的情况下，通过执行与图13所示的流程图的步骤S41 ~ S48同样的处理，来计算心尖部水平与乳头肌水平之间的长轴方向的变形8zl2和乳头肌水平与心基部水平之间的长轴方向的变形Sz23。
这里，针对两个时间相位tl和t2 (tl弁2)的心尖部短轴图4象Gl、乳头肌短轴图像G2和心基部短轴图像G3来计算变形8zl2和Sz23。运动信息计算部分62利用时间间隔At = I tl - t2l来除变形8zl2，
以计算心尖部水平和乳头肌水平之间的长轴方向的变形率。另外，运动信息计算部分利用时间间隔At来除变形Sz23,以计算乳头肌水平与心基部水平之间的长轴方向的变形率。运动信息计算部分可以执行通常的微分运算来从变形计算变形率。
控制器9控制显示部分81,以显示所计算的心脏壁的变形率8zl2/At和5z23/At，作为表示心脏壁的变形的时间变化率的运动信息。用户可以通过参照所显示的心脏壁的变形率，来把握心脏壁的变形的时间变化率。
还可以分别计算心脏壁的内膜和心脏壁的外膜的变形率，并根据两个变形率的值来评价变形的时间变化率(例如得到两个变形率的值的平均)。
上述实施方式的超声波诊断装置1构成为由用户操作操作部分 82以指定测量图像区域。本变形自动地进行测量图像区域指定操作。
图15示出本变形的超声波诊断装置的例子。图15的超声波诊断 装置100具有与上述实施方式的超声波诊断装置l基本相同的结构， 不同点在于，图像处理器5设置有测量图像区域指定部分53。另外， 为了执行本变形的特征性处理，设置有不同于上述实施方式的医用图 像处理程序72。
测量图像区域指定部分53分析由MPR处理部分52产生的断层 图像(例如心尖部、乳头肌、心基部的短轴图像)的图像数据，并且 在断层图像上指定测量图像区域。
更具体地，测量图像区域指定部分53例如实施通常的边界提取 处理，以分析图像数据的每个像素的像素值，从而提取生物组织的边 界部分(例如心脏壁的内膜和外膜)。然后，测量图像区域指定部分 53指定所提取的边界部分作为断层图像的测量图像区域。测量图像区 域指定部分53对应于本发明的"指定部分"的例子，并且包括例如根 据医用图像处理程序72进行操作的微处理器。
由测量图像区域指定部分53执行的处理对应于图4的流程图的步骤S09的过程。当测量图像区域指定部分53进行测量图像区域指 定处理时，可以在显示部分81上显示断层图像，也可以不显示。
另外，为了使用户能够检查自动指定的测量图像区域，最好在显 示部分81上显示输入了测量图像区域的断层图像。此时，最好构成 为可以使用操作部分82适当改变测量图像区域。
由于在本变形的结构中自动指定测量图像区域，因此可以容易地 指定测量图像区域(实际上，当原样使用自动指定的测量图像区域时， 不必执行测量图像区域指定操作)。
在上述实施方式的超声波诊断装置1中，使用将超声波换能器二 维排列的二维超声波换能器来电子地执行三维超声波扫描操作。本变 形涉及超声波换能器装备有在一维方向上排列的一维超声波探头的 超声波诊断装置。
当使用一维超声波探头时，仅在一维方向(图2和图3所示的主 扫描方向X)上电子地执行超声波扫描操作。因此，在三维超声波扫 描操作中，可以在副扫描方向Y上手动或机械地执行扫描操作。
即使使用这样的一维超声波探头，也可以构成为生成基于三维超 声波扫描操作的体数据，根据体数据生成断层图像的图像数据，显示 该断层图像并且指定测量图像区域。
[第四变形
上述实施方式的超声波诊断装置1计算在心尖部水平和乳头肌 水平等两个以上的断面位置处的断层图像上指定的测量图像区域的
按时间序列的位移，并且根据两个以上断面位置处的测量图像区域的 位移来计算生物组织的运动信息。即，在上述实施方式中，仅根据所 指定的测量图像区域的位移来计算运动信息。
本变形自动地另外指定所指定的测量图像区域以外的测量图像 区域，并且通过包括自动指定的测量图像区域的位移在内来计算移动 信息。以下说明本变形的超声波诊断装置。与图15所示的超声波诊 断装置100同样，本变形的超声波诊断装置包括测量图像区域指定部分53。
说明本变形的超声波诊断装置的操作。与上述实施方式同样，作 为在显示部分81上显示的心尖部短轴图像Gl、乳头肌短轴图像G2 和心基部短轴图像G3的测量图像区域，用户例如指定内膜位置图像 ml、 m2、 m3和外膜位置图像Ml、 M2、 M3 (见图5和图6)。可 以通过第二变形的自动指定处理来指定测量图像区域。
测量图像区域指定部分53指定所指定的测量图像区域ml ~ m3 和M1~M3以外的测量图像区域。例如，如图16所示，测量图像区 域指定部分53在心尖部短轴图像G1外侧的位置处指定内膜位置图像 m4和外膜位置图像M4，在心尖部短轴图像Gl和乳头肌短轴图像 G2之间的位置处指定内膜位置图像m5和外膜位置图像M5,在乳头 肌短轴图像G2和心基部短轴图像G3之间的位置处指定内膜位置图 像m6和外膜位置图像M6。
具体说明测量图像区域指定处理。测量图像区域指定部分53首 先根据由用户等指定了测量图像区域的心尖部短轴图像Gl、乳头肌 短轴图像G2和心基部短轴图像G3的断面位置(坐标)，确定用于 指定新的测量图像区域的断面位置(坐标)。该处理例如如下执行 得到在与乳头肌短轴图像G2相反的方向上从心尖部短轴图像Gl离 开预定距离的位置的坐标(第一Z坐标)，根据心尖部短轴图像G1 的Z坐标和乳头肌短轴图像G2的Z坐标，得到断面中心的坐标(第 二 Z坐标)，然后得到乳头肌短轴图像G2的Z坐标和心基部短轴图 像G3的Z坐标的中心的坐标(第三Z坐标)。
接着，MPR处理部分52根据体数据在第一 第三Z坐标生成 平行于心尖部短轴图像G1等的断层图像的图像数据。
测量图像区域指定部分53分析第一 Z坐标处的断层图像的图像 数据的像素值，以提取心脏壁的边界部分，从而将内侧边界部分指定 为内膜位置图像m4 ，并且将外侧边界部分指定为外膜位置图像M4 。 同样，测量图像区域指定部分分析第二Z坐标处的断层图像的图像数 据的像素值，以提取心脏壁的边界部分，从而将内侧边界部分指定为内膜位置图像m5，并且将外侧边界部分指定为外膜位置图像M5。测 量图像区域指定部分分析第三Z坐标处的断层图像的图像数据的像 素值，以提取心脏壁的边界部分，从而将内侧边界部分指定为内膜位 置图像m6，并且将外侧边界部分指定为外膜位置图像M6。然后，将 内膜位置图像m4~m6和外膜位置图像M4 ~ M6分别指定为新的测 量图像区域。
本变形的位移计算部分61对用户等指定的测量图像区域ml~ m3和M1~M3、由测量图像区域指定部分53指定的测量图像区域 m4 ~ m6和M4 ~ M6分别计算按时间序列的位移。
运动信息计算部分62根据所有测量图像区域ml ~ m6和Ml ~ M6的位移来计算运动信息。例如，当评价心脏壁的扭转运动时，运 动信息计算部分计算内膜位置图像m4和内膜位置图像ml的相对旋 转角度、内膜位置图像ml和内膜位置图像m5的相对旋转角度、内 膜位置图像m5和内膜位置图像m2的相对旋转角度、内膜位置图像 m2和内膜位置图像m6的相对旋转角度、内膜位置图像m6和内膜位 置图像m3的相对旋转角度。同样，运动信息计算部分能够计算外膜 位置图像M1 M6的相对旋转角度。
控制器9控制显示部分81，以显示基于所计算的相对旋转角度 的运动信息。这样，根据本变形，可以得到比上述实施方式更精确的 运动信息。
[第五变形
在上述实施方式中，由用户指定用于指定测量图像区域的断面位 置(参见图5的断面位置指定光标C1~C3)。但是，可以构成为自 动指定断面位置。
例如，象术前术后观察或经过观察那样，当存在过去的同一生物 组织的超声波图像的图像数据时，存储所指定的断面位置并且读出过 去指定的断面位置的信息，从而可以自动指定当前的断面位置。
另外，可以预先设定心尖部水平、乳头肌水平和心基部水平等典 型的断面位置，并且例如根据B模式图像的体数据或图像数据来确定
40典型的断面位置。
另外，可以构成为根据B模式图像的图像数据等分析生物组织 的大小(例如心脏的长轴方向的长度)，并且将把大小分成多个部分 的断面位置指定为用于指定测量图像区域的断面位置。
[其它
在上述实施方式中，得到生物组织的位移并且通过对位移进行微 分(用时间间隔来除)而得到速度。但是，可以构成为首先得到速度， 通过对速度进行积分来得到位移。
另外，在上述实施方式中，得到并显示表示生物组织的(更大范
围的)运动状态的运动信息，还可以构成为显示每个测量图像区域的
局部运动信息。在这种情况下，当不计算或显示运动信息时，可以应 用具有计算和显示局部运动信息所需的结构的超声波诊断装置(对于
后述的医用图像处理装置和医用图像处理程序也同样)。 <第二实施方式>
说明本发明的医用图像处理装置。医用图像处理装置例如包括与 超声波诊断装置连接的计算机、以及与用于存储超声波图像的图像数 据的PACS ( Picture Archiving and Communication System, 医学图 像归档和通讯系统)等数据库连接的计算机。图17示出本发明的医 学图像处理装置的例子。在图17中，对与第一实施方式相同的部件 附加相同的附图标记。
图17所示的医用图像处理装置1000包括与图l所示的超声波诊 断装置1同样的图像处理器5、计算处理器6、存储器7、用户接口8 和控制器9。存储器7作为本发明的"存储器"的例子。医用图像处理 装置1000经由LAN ( Local Area Network,局域网)等网络N与超 声波诊断装置2000和医用图像数据库3000连接。本实施方式的控制 器9包括用于经由网络N执行数据通信的网络适配器。
说明医用图像处理装置1000的操作。从超声波诊断装置2000 或医用图像数据库3000向医用图像处理装置1000输入超声波图像的当所输入的图像数据是B模式图像的图像数据时，体数据生成 部分51生成基于图像数据的体数据。与第一实施方式同样，MPR处 理部分52根据体数据生成按时间序列的一 系列断层图像的图像数据。 一系列断层图像的图像数据被存储在存储器7中。
另外，当输入的图像数据是体数据时，MPR处理部分52根据体 数据生成按时间序列的一系列断层图像的图像数据。 一系列断层图像 的图像数据被存储在存储器7中。
医用图像处理装置1000根据存储在存储器7中的一系列断层图 像的图像数据(以及体数据)执行与第一实施方式同样的处理(参见 图4~图14)。从而，可以在短时间内测量生物组织的三维运动。另 外，可以容易地指定测量生物组织的运动的区域。
作为本实施方式的变形，可以应用设置有测量图像区域指定部分 53的图像处理器5(参见图15，基于医用图像处理程序72的处理)。 从而可以执行与第一实施方式的第二变形和第四变形同样的处理。
<第三实施方式>
说明本发明的医用图像处理程序。在第一和第二实施方式中说明 的医用图像处理程序71和72是本发明的医用图像处理程序的例子。 医用图像处理程序用于使计算机执行在第 一 实施方式及其变形中说 明的处理。医用图像处理程序可以预先存储在嵌入计算机中的硬盘驱 动器等存储设备中，或者可以预先存储在LAN等网络上的服务器中， 从而计算机读出并执行该程序。
可以将医用图像处理程序存储在任意的存储介质中，以便可由计 算机读出。存储介质的例子例如包括光盘、光磁盘(CD-ROM、 DVD-RAM、 DVD-ROM、 MO等)、磁存储介质(硬盘、软盘、ZIP 等)以及半导体存储器。
4权利要求
1. 一种超声波诊断装置，包括超声波探头；收发器，向超声波探头发送和从超声波探头接收超声波；图像生成部分，根据作为发送和接收超声波的结果而得到的接收信号，针对生物组织的两个以上的断面位置中的每一个断面位置生成按时间序列的一系列断层图像的图像数据；显示部分，针对所述两个以上的断面位置中的每一个断面位置，根据所生成的图像数据，显示所述一系列断层图像中的一个断层图像；指定部分，针对所述两个以上的断面位置中的每一个断面位置，在所显示的所述一个断层图像上指定测量图像区域；和计算器，针对所述两个以上的断面位置中的每一个断面位置，根据所述一系列断层图像的图像数据，计算表示所指定的测量图像区域中的生物组织的运动状态的局部运动信息，并且根据针对所述两个以上的断面位置中的每一个断面位置计算的局部运动信息，计算表示生物组织的运动状态的运动信息，其中，显示部分显示由计算器计算的运动信息。
2. 如权利要求l所述的超声波诊断装置，其中，超声波探头利用收发器，在沿着与预定方向正交的方向的超声波 接收和发送方向上对沿着该预定方向的多个位置中的每一个位置进 行扫描，并且重复沿着所述预定方向和所述正交的方向的超声波接收 和发送操作，图像生成部分包括第一断层图像生成部分，根据在所述正交的方向的所述超声波接 收和发送方向上进行扫描时依次得到的接收信号，生成生物组织的包 括所述发送和接收方向以及所述正交的方向的断面中的第一断层图 像的图像数据；体数据生成部分，根据针对所述多个位置中的每一个位置生成的第一断层图像的图像数据，生成生物组织的体数据；和第二断层图像生成部分，根据所生成的体数据，针对所述两个以上的断面位置中的每一个断面位置生成第二断层图像的图像数据， 其中，通过针对所述两个以上的断面位置中的每一个断面位置重复进行超声波接收和发送操作而得到第二断层图像的多个图像数据被用作所述按时间序列的 一 系列断层图像的图像数据。
3. 如权利要求l所述的超声波诊断装置，其中， 计算器通过计算表示所述两个以上的断面位置中的两个断面位置处的局部运动信息之差的差信息，来得到运动信息。
4. 如权利要求3所述的超声波诊断装置，其中， 所述两个断面位置中的一个断面位置处的所述一系列断层图像的断面与另一个断面位置处的所述一系列断层图像的断面平行，计算器针对所述两个断面位置中的每一个断面位置计算测量图像区域中的生物组织绕着与所述一系列断层图像的断面正交的方向 的旋转角度，来作为局部运动信息，并且计算针对所述两个断面位置 计算的旋转角度之差，来作为差信息，从而使用所计算的旋转角度之 差作为表示生物组织的扭转运动的运动信息。
5. 如权利要求3所述的超声波诊断装置，其中， 所述两个断面位置中的一个断面位置处的所述一系列断层图像的断面与另一个断面位置处的所述一系列断层图像的断面平行，计算器针对所述两个断面位置中的每一个断面位置计算测量图像区域中的生物组织绕着与所述一系列断层图像的断面正交的方向 的旋转角度，来作为局部运动信息，并且计算针对所述两个断面位置 计算的旋转角度之差，来作为差信息，利用时间对所计算的旋转角度 之差进行微分，从而计算表示生物组织的扭转运动的速度的运动信 每
6. 如权利要求3所述的超声波诊断装置，其中， 所述两个断面位置中的一个断面位置处的所述一系列断层图像的断面与另一个断面位置处的所述一系列断层图像的断面平行，计算器针对所述两个断面位置中的每一个断面位置计算测量图像区域中的生物组织绕着与所述一系列断层图像的断面正交的方向 的旋转角度，来作为局部运动信息，并且计算针对所述两个断面位置 计算的旋转角度之差，来作为差信息，利用所述两个断面位置之间的 距离来除所计算的旋转角度之差，从而得到表示生物组织的相对旋转 梯度的运动信息。
7. 如权利要求2所述的超声波诊断装置，其中， 计算器针对所述两个以上的断面位置中的每一个断面位置，根据体数据计算表示由指定部分指定的测量图像区域中的生物组织的三维运动状态的局部运动倌息，并且根据针对所述两个以上的断面位置 中的每一个断面位置计算的表示三维运动状态的局部运动信息，来计 算运动信息。
8. 如权利要求7所述的超声波诊断装置，其中， 计算器通过计算表示差的差信息来得到运动信息，其中该差是表示所述两个以上的断面位置中的两个断面位置处的三维运动状态的 局部运动信息之差。
9. 如权利要求8所述的超声波诊断装置，其中， 所述两个断面位置中的一个断面位置处的所述一系列断层图像的断面与另一个断面位置处的所述一系列断层图像的断面平行，计算器针对所述两个断面位置中的每一个断面位置，根据所计算的表示三维运动状态的局部运动信息，计算测量图像区域中的生物组 织在与所述 一 系列断层图像的断面正交的方向上的位移，并且计算针 对所述两个断面位置中的每一个断面位置计算的所述正交的方向上 的位移之差，来作为差信息，从而使用所计算的位移之差作为表示生 物组织在所述正交的方向上的伸缩运动的运动信息。
10. 如权利要求8所述的超声波诊断装置，其中， 所述两个断面位置中的一个断面位置处的所述一系列断层图像的断面与另一个断面位置处的所述一系列断层图像的断面平行，计算器针对所述两个断面位置中的每一个断面位置，根据所计算 的表示三维运动状态的局部运动信息，计算测量图像区域中的生物组 织在与所述一 系列断层图像的断面正交的方向上的位移，并且计算针 对所述两个断面位置中的每一个断面位置计算的所述正交的方向上 的位移之差，来作为差信息，利用时间对所计算的位移之差进行微分， 从而计算表示生物组织在所述正交的方向上的伸缩运动的速度的运 动信息。
11. 如权利要求8所述的超声波诊断装置，其中， 所述两个断面位置中的一个断面位置处的所述一系列断层图像的断面与另一个断面位置处的所述一系列断层图像的断面平行，计算器针对所述两个断面位置中的每一个断面位置，根据所计算的表示三维运动状态的局部运动信息，计算测量图像区域中的生物组 织在与所述一系列断层图像的断面正交的方向上的位移，并且计算针 对所述两个断面位置中的每一个断面位置计算的所述正交的方向上 的位移之差，来作为差信息，利用所述两个断面位置之间的原始距离 来除所计算的位移之差，从而使用其商来作为表示生物组织在所述正 交的方向上的变形的运动信息。
12. 如权利要求8所述的超声波诊断装置，其中， 所述两个断面位置中的一个断面位置处的所述一系列断层图像的断面与另一个断面位置处的所述一系列断层图像的断面平行，计算器针对所述两个断面位置中的每一个断面位置，根据所计算 的表示三维运动状态的局部运动信息，计算测量图像区域中的生物组 织在与所述一系列断层图像的断面正交的方向上的位移，并且计算针 对所述两个断面位置中的每一个断面位置计算的所述正交的方向上 的位移之差，来作为差信息，利用所述两个断面位置之间的原始距离 来除所计算的位移之差，并且利用时间对其商进行微分，从而计算表 示生物组织在所述正交的方向上的变形率的运动信息。
13. 如权利要求1所述的超声波诊断装置，其中， 指定部分包括自动指定部分，该自动指定部分根据在显示部分显示的所述一个断层图像上指定的测量图像区域，在与所述两个以上的断面位置不同的断面位置处指定新的测量图像区域，计算器计算表示新指定的测量图像区域中的生物组织的运动状态的局部运动信息，并且根据针对所述两个以上的断面位置和所述不 同的断面位置计算的局部运动信息，来计算运动信息。
14. 如权利要求1 13所述的超声波诊断装置，其中， 所述生物组织是心脏，由指定部分指定的测量图像区域是在所述显示部分显示的所述 一个断层图像中与心脏的心脏壁对应的图像区域。
15. —种超声波诊断装置，包括 超声波探头；收发器，向超声波探头发送和从超声波探头接收超声波； 图像生成部分，根据作为发送和接收超声波的结果而得到的接收信号，针对生物组织的一个以上的断面位置中的每一个断面位置生成按时间序列的一系列断层图像的图像数据；显示部分，针对所述一个以上的断面位置中的每一个断面位置，根据所生成的图像数据，显示所述一系列断层图像中的一个断层图像；指定部分，针对所述一个以上的断面位置中的每一个断面位置， 在所显示的所述一个断层图像上指定测量图像区域；和计算器，针对所述一个以上的断面位置中的每一个断面位置，根 据所述一系列断层图像的图像数据，计算表示所指定的测量图像区域 中的生物组织的运动状态的局部运动信息，其中，显示部分显示由计算器计算的局部运动信息。
16. 如权利要求15所述的超声波诊断装置，其中， 所述生物组织是心脏，所指定的测量图像区域是所显示的所述一个断层图像中的心脏壁，计算器计算以下信息中的至少一个作为局部运动信息心脏壁厚的变化、心脏壁厚的变化的速度、心脏壁在厚度方向上的变形、心脏 壁在厚度方向上的变形率、心脏壁的内膜或外膜在绕着与所述一系列 断层图像正交的方向的旋转方向上的旋转角度、心脏壁在所述旋转方 向上的变形、心脏壁在所述旋转方向上的变形率、心脏壁的内膜或外 膜在所述旋转方向上的相对旋转角度。
17. 如权利要求1或15所述的超声波诊断装置，其中， 计算器在所显示的所述一个断层图像上指定的测量图像区域中指定多个测量点，得到所述多个测量点在所述 一 系列断层图像中的每 一个断层图像中的位置，根据所得到的位置计算所述多个测量点中的 每一个测量点的按时间序列的位移或位移速度，并且根据针对所述多个测量点中的每一个测量点计算的位移或位移速度来计算测量图像 区域的位移。
18. 如权利要求1或15所述的超声波诊断装置，其中， 超声波探头包括二维排列的多个超声波换能器， 收发器个别地或者以两个以上的组来控制所述多个超声波换能器，以便在预定方向和与所述预定方向正交的方向上利用超声波进行 扫描。
19. 一种医用图像处理装置，处理由超声波诊断装置得到的生物 组织的医用图像的图像数据，所述医用图像处理装置包括存储器，存储生物组织的两个以上的断面位置处的、按时间序列 的一系列断层图像的图像数据；显示部分，针对所述两个以上的断面位置中的每一个断面位置， 根据所存储的图像数据，显示所述一系列断层图像中的一个断层图 像；指定部分，针对所述两个以上的断面位置中的每一个断面位置， 在所显示的所述一个断层图像上指定测量图像区域；和计算器，针对所述两个以上的断面位置中的每一个断面位置，根 据所述一系列断层图像的图像数据，计算表示所指定的测量图像区域 中的生物组织的运动状态的局部运动信息，并且根据针对所述两个以上的断面位置中的每一个断面位置计算的局部运动信息，计算表示生 物组织的运动状态的运动信息，其中，显示部分显示由计算器计算的运动信息。
20. —种医用图像处理装置，处理由超声波诊断装置得到的生物 组织的医用图像的图像数据，所述医用图像处理装置包括存储器，存储生物組织的一个以上的断面位置中的每一个断面位 置处的、按时间序列的一系列断层图像的图像数据；显示部分，针对所述一个以上的断面位置中的每一个断面位置， 根据所存储的图像数据，显示所述一系列断层图像中的一个断层图像；指定部分，针对所迷一个以上的断面位置中的每一个断面位置， 在所显示的所述一个断层图像上指定测量图像区域；和计算器，针对所述一个以上的断面位置中的每一个断面位置，根 据所述一系列断层图像的图像数据，计算表示所指定的测量图像区域 中的生物组织的运动状态的局部运动信息，其中，显示部分显示由计算器计算的局部运动信息。
21. —种医用图像处理程序，使具有存储器和显示部分的计算机 执行以下功能，其中该存储器存储生物组织的两个以上的断面位置中 的每一个断面位置处的、按时间序列的一系列断层图像的图像数据针对所述两个以上的断面位置中的每一个断面位置，根据所存储 的图像数据，在所述显示部分上显示所述一系列断层图像中的一个断 层图像；针对所述两个以上的断面位置中的每一个断面位置，响应于在所 显示的所述一个断层图像上指定了测量区域，根据所述一 系列断层图 像的图像数据，计算表示所指定的测量图像区域中的生物组织的运动 状态的局部运动信息；根据针对所述两个以上的断面位置中的每一个断面位置计算的 局部运动信息，计算表示生物组织的运动状态的运动信息；在显示部分上显示所计算的运动信息。
22. —种医用图像处理程序，使具有存储器和显示部分的计算机 执行以下功能，其中该存储器存储生物组织的一个以上的断面位置中 的每一个断面位置处的、按时间序列的一系列断层图像的图像数据针对所述一个以上的断面位置中的每一个断面位置，根据所存储 的图像数据，在所述显示部分上显示所述一系列断层图像中的一个断 层图像；针对所述一个以上的断面位置中的每一个断面位置，响应于在所 显示的所述一个断层图像上指定了测量区域，根据所述一 系列断层图 像的图像数据，计算表示所指定的测量图像区域中的生物组织的运动 状态的局部运动信息；在显示部分上显示所计算的局部运动信息。
全文摘要
本发明提供一种可以在短时间内测量生物组织的三维运动的超声波诊断装置。图像处理器(5)根据生物组织的B模式图像的图像数据生成体数据，并且根据体数据并针对两个以上的断面位置中的每一个断面位置分别生成按时间序列的一系列断层图像的图像数据；控制器(9)在显示部分(81)上针对每个断面位置显示一个断层图像；用户操作操作部分(82)，在所显示的断层图像上指定测量图像区域；位移计算部分(61)针对每个断面位置计算所指定的测量图像区域的按时间序列的位移；运动信息计算部分(62)根据针对每个断面位置计算的测量图像区域的位移，计算生物组织的运动信息。
文档编号A61B8/08GK101484074SQ200780024890
公开日2009年7月15日 申请日期2007年5月30日 优先权日2006年5月30日
发明者川岸哲也, 阿部康彦 申请人:株式会社东芝;东芝医疗系统株式会社