CN103118600A - 超声波诊断装置以及超声波图像的显示方法 - Google Patents

Abstract

本发明提供一种超声波诊断装置及超声波图像的显示方法，为了清晰地描绘出应关注的硬度部位的三维弹性像，针对构成体数据的弹性值数据，按弹性值的每个大小对频度进行计数，并进行使基于频度确定的规定范围的弹性值的值变大、使规定范围外的弹性值的值变小的变换。通过对变换后的体数据进行绘制，从而应关注的硬度部位的大的值被积算，生成三维弹性像，从而能够清晰地描绘出应关注的硬度部位。

Description

超声波诊断装置以及超声波图像的显示方法

技术领域

本发明涉及利用超声波对被检体内的诊断部位显示超声波图像的超声波诊断装置，尤其涉及能够将表示被检体的生物体组织的硬度的弹性图像作为三维图像来显示的超声波诊断装置。

背景技术

以往，如专利文献1和2所示，超声波诊断装置通过将超声波传送到被检体内，从其反射波接收生物体组织的超声波反射回波信号，并实施信号处理，从而生成并显示将超声波反射率作为亮度的诊断部位的黑白断层像(B模式像)。

此外，在专利文献2中，公开了如下方法：根据黑白断层像和其取得位置来进行三维坐标变换，从而获得将多个断层像数据进行了三维配置的三维体数据(volume data)，并对其进行体绘制(rendering)，从而获得从任意的视线方向看的诊断部位的三维浓淡像。作为绘制的方法，公开了如下方法：根据构成三维体数据的体素(voxel)的亮度的值来对每个体素赋予不透明度，且直到在视线上排列的体素的不透明度的累积值成为1为止，将其视线上的体素的亮度值依次累积的值作为二维投影面上的像素(pixel)值。

此外，在专利文献1和2中，还公开了求出诊断部位的弹性像的方法。首先，选择2个帧的黑白断层像，通过模块匹配(bloc matching)等求出图像上的各点在2个帧之间位移了多少，并对求出的位移实施已知的运算，从而求出关于图像上的各个点的表示硬度的弹性值(变形或弹性率等)。将求出的弹性值的大小根据色变换表来变换为色调信息，从而获得用色调表示弹性值的二维弹性像。

此外，在专利文献2中，公开了如下方法：为了在三维浓淡像中，解决被外侧的不透明度高的体素所隐藏而得不到内侧的病变部的三维像的问题，对从黑白断层像获得的三维体数据的各个体素，根据弹性值的大小来赋予不透明度，进行体绘制。由此，获得将硬的组织表示成更加不透明、将软的组织表示成更加透明的三维浓淡像。

此外，在专利文献2中，对多个二维弹性像进行坐标变换而获得三维体数据，并对其进行体绘制，从而生成三维弹性像。此时，由于在将体素值简单累积的通常的绘制方法中，产生弹性值被累积、生物体组织的弹性特性被歪曲的三维弹性像的问题，所以使用黑白断层像的亮度数据来求出在视线上贡献率最高的体素，将其体素的弹性值作为二维投影面的弹性值来生成三维弹性像。

作为显示方法，在专利文献1中，重叠显示黑白断层像(二维)和二维弹性像。在专利文献2中，示出了生成将黑白三维像和三维弹性像重叠的图像，并将其与三个方向的黑白断层图像在一个画面上排列后进行显示。

此外，在专利文献3中，公开了在一个画面上排列显示黑白断层像和二维弹性像或重叠显示黑白断层像和二维弹性像的方法。

现有技术文献

专利文献

专利文献1：国际公开WO2005/048847号公报

专利文献2：JP特开2008-259605号公报

专利文献3：JP特开2000-60853号公报

发明内容

(发明要解决的课题)

在以往的求出被检体的弹性像的方法中，如专利文献1和2所记载的那样，选择2个帧的黑白断层像，且对被检体施加压力，从而求出图像上的各个点在2个帧之间位移了多少，并基于求出的位移来求出弹性值(变形或弹性率等)。因此，越是容易产生较大位移的软的组织，越会成为弹性值越大的值，越是难以产生位移的硬的组织，越会成为弹性值越小的值。

操作者想要更加清晰地看到的组织大多数为肿瘤等硬的组织，但硬的组织的弹性值小，包围硬的组织的软的组织的弹性值大。在以往的生成三维浓淡像的体绘制方法中，根据构成三维体数据的体素的亮度的值来对每个体素赋予不透明度，且直到在视线上排列的体素的不透明度的累积值成为1为止，将依次累积了其视线上的体素的亮度值的值作为二维投影面上的像素(pixel)值。因此，若将其用于弹性值的体数据的绘制，则弹性值的数值大的软的体素对二维投影面的贡献率变高，难以描绘出位于软的组织内侧的应关注的硬的组织。

假设将弹性值的数值大的软的组织的体素的数据通过附加过滤器等而除去了的情况下，能够获得硬的组织的体素的三维像，但由于硬的组织的体素的弹性值的数值小，所以通过绘制而进行了积算之后的数值也小，难以获得对比度高的清晰的三维像。

在专利文献2中，为了获得内侧的病变部的三维像，对从黑白断层像获得的三维体数据的各个体素，根据弹性值的大小来赋予不透明度后进行绘制，但由于在该方法中进行积算的是黑白断层像的亮度数据(超声波反射率)，所以得不到三维弹性像。

本发明的目的在于，提供一种能够清晰地描绘出应关注的硬度的部位的三维弹性像的超声波诊断装置。

(用于解决课题的手段)

为了达到上述目的，根据本发明，提供如下的超声波诊断装置。即，该超声波诊断装置包括：弹性像体数据生成部，其使用向被检体内发送超声波而接收到的信号，生成由表示弹性的弹性值数据构成的体数据；弹性值变换部，其针对构成体数据的弹性值数据，按弹性值的每个大小对频度进行计数，并进行使基于频度确定的规定范围的弹性值的值变大、使规定范围外的弹性值的值变小的变换；以及体绘制部，其通过对弹性值变换部进行变换后的体数据进行绘制，从而生成三维弹性像。

此外，根据本发明的其他方式，提供超声波图像的显示方法。即，该超声波图像的显示方法使用向被检体内发送超声波而接收到的信号，生成由表示弹性的弹性值数据构成的体数据，针对构成体数据的弹性值数据，按弹性值的每个大小对频度进行计数，并进行使基于频度确定的规定范围的弹性值的值变大、使规定范围外的弹性值的值变小的变换，通过对变换后的体数据进行绘制，从而生成三维弹性像来进行显示。

(发明效果)

根据本发明，基于频度确定的应关注的硬度部位的弹性值的数值增大，从而在通过体绘制进行了积算的情况下，应关注的硬度部位的值增大，能够生成应关注的硬度部位的对比度大的三维弹性像。

附图说明

图1是表示实施方式的超声波诊断装置的结构的框图。

图2是表示实施方式1的弹性值变换部的动作的流程图。

图3(a)是表示实施方式1的弹性值变换部生成的弹性值与频度的关系的曲线图，(b)是表示在实施方式1中弹性值变换部进行变换后的弹性值(相对值)与变换前的弹性值(相对值)的关系的曲线图。

图4(a)是表示在实施方式1中获得的三维弹性像、以及显示了二维弹性像与黑白断层像的合成图像的画面例的说明图，(b)是图(a)的色彩图15的放大图，(c)是表示比较例的三维弹性像、以及显示了二维弹性像与黑白断层像的合成图像的画面例的说明图。

图5(a)是表示实施方式2的弹性值变换部生成的弹性值与频度的关系的曲线图，(b)是表示在实施方式2中弹性值变换部进行变换后的弹性值(相对值)与变换前的弹性值(相对值)的关系的曲线图。

图6是表示在实施方式2中获得的三维弹性像的说明图。

图7(a)是表示实施方式3的弹性值变换部生成的弹性值与频度的关系的曲线图，(b)是表示在实施方式3中弹性值变换部进行变换后的弹性值(相对值)与变换前的弹性值(相对值)的关系的曲线图。

图8(a)是表示实施方式4的弹性值变换部生成的弹性值与频度的关系的曲线图，(b)是表示在实施方式4中弹性值变换部生成的权重w和变换前的弹性值(相对值)的关系的曲线图。

具体实施方式

以下，基于附图说明本发明的一实施方式的超声波诊断装置。在本实施方式中，将表示被检体的规定截面中的组织的超声波反射率的分布的B模式像称为黑白断层像，将对由黑白断层像的数据构成的体数据进行了绘制的二维投影像称为三维像，将表示对被检体的组织的弹性进行表示的弹性值的规定截面中的二维分布的像称为二维弹性像，将对由二维弹性像的数据构成的体数据进行了绘制的二维投影像称为三维弹性像。

如图1所示，在超声波诊断装置中，包括：探头2，使其与被检体1抵接后使用；发送部3，经由探头2向被检体1内的诊断部位以规定的时间间隔重复发送超声波；接收部4，以时间序列接收由被检体1反射的反射回波信号；超声波接收控制部5；以及调相加法部6，对接收到的反射回波进行调相加法运算。

探头2包括以列状或者扇形排列的多个振子，利用振子对被检体1进行超声波的收发。探头2具有在与排列方向正交的方向(短轴方向)上对多个阵子进行机械扫描的功能，能够以三维方式收发超声波。另外，作为探头2，也可以通过使用将多个振子二维排列的探头，从而不需要使振子机械振动，就能够以三维方式收发超声波。

发送部3驱动探头2来生成用于产生超声波的发送波脉冲。此时，控制传递给探头2的各个振子的发送波信号的相位，将发送的超声波的收敛点设定为某一深度。此外，接收部4以规定的增益对探头2的各个振子接收到的反射回波信号进行放大，从而生成RF信号、即接收波信号。超声波收发控制部5控制发送部3和接收部4。调相加法部6通过对RF信号的相位进行调整之后进行加法运算，从而形成相对于一个点或者多个收敛点收敛的超声波束，由此生成RF信号帧数据(相当于RAW数据)。

此外，在超声波诊断装置中，包括：界面部43，接受来自操作者的设定；图像显示部13；切换加法部12，切换在图像显示部13中显示的图像的种类；图像系统控制部44；以及短轴扫描位置控制部46。短轴扫描位置控制部46通过控制探头2在与排列方向正交的方向上对多个振子进行机械扫描的动作，从而进行规定范围的三维收发。

此外，作为用于从RF帧数据生成以超声波反射率作为亮度的诊断部位的黑白断层像以及三维像的结构，超声波诊断装置包括断层像构成部7、断层像存储部35、断层像的体数据生成部36、体绘制部38以及多帧构成部(断层像)46。除此之外，作为用于从RF信号帧数据生成二维弹性像以及三维弹性像的结构，包括RF信号帧数据存储部27、RF信号帧数据选择部28、位移测量部30、弹性信息运算部32、弹性图像构成部34、二维弹性像存储部39、弹性像的体数据生成部40、体绘制部42、多帧构成部(弹性像)48以及弹性值变换部51。

断层像构成部7对由调相加法部6生成的RF信号帧数据实施增益修正、日志(log)压缩、检波、轮廓强调以及过滤处理等，并生成以反射率作为亮度(浓淡)的诊断部位的黑白断层像(B模式像)。此时，图像系统控制部44经由界面部43从操作者接受黑白断层像的生成条件，控制断层像构成部7。

断层像存储部35使在断层像构成部7中构成的黑白断层像和其取得位置建立对应之后进行存储。这里所指的取得位置是基于短轴扫描位置控制部46的控制的短轴方向的移动量。体数据生成部36对在断层像存储部35中存储的多个黑白断层像(与一个立体相应的量)进行根据短轴方向移动量来再次配置的坐标变换，生成三维体数据。

体绘制部38使用下式(1)～(3)，基于亮度和不透明度对体数据生成部36生成的三维体数据进行绘制，从而构成被检体的诊断部位的三维像(三维体数据的二维投影像)。投影方向(视线方向)是图像系统控制部44经由界面43从操作者接受的。

Cout(i)＝Cout(i-1)+(1-Aout(i-1))·A(i)·C(i)·S(i)···式(1)

Aout(i)＝Aout(i-1)+(1-Aout(i-1))·A(i)···式(2)

A(i)＝Bopacity[C(i)]···式(3)

在式(1)中，Cout(i)是作为二维投影面的像素值来输出的值。C(i)是从二维投影面上的某一点看三维像时在视线上第i个(其中，i＝0～N-1)存在的体素的亮度值。另外，这里所称的体素是指构成三维体数据的各个亮度数据的位置。在视线上排列了N个体素时，将i＝0～N-1为止的体素的亮度值根据式(1)进行积算而得到的亮度值Cout(N-1)成为最终输出的像素值。Cout(i-1)表示体素第i-1个为止的积算值。

式(1)的A(i)为在视线上第i个存在的体素的不透明度，为0～1.0的值。如式(3)所示，通过参照确定了预先决定的亮度值C(i)与不透明度(Opacity)的关系的表(Bopacity[C(i)])，或者向确定了预先决定的亮度值C(i)与不透明度的关系的函数(Bopacity[C(i)])代入亮度值C(i)，从而根据体素的亮度值的大小而确定不透明度A(i)。例如，对亮度值大的体素赋予大的不透明度。由此，通过根据亮度值来赋予不透明度，从而决定对要输出的二维投影面的亮度值Cout(N-1)的、其体素的亮度值C(i)的贡献率。

式(2)的Aout(i)为，将通过式(3)赋予的不透明度A(i)直到第i个体素为止根据式(2)的右边进行了积算的值。在式(1)中，使用如式(2)所示那样计算出的第i-1个体素为止的不透明度的积算值Aout(i-1)。根据式(2)可知，Aout(i)在每次通过体素时被积算且收敛为1.0。因此，如上述(1)所示那样在第i-1个为止的不透明度的积算值Aout(i-1)约为1.0的情况下，式(1)的右边第2项成为0，第i个以后的亮度值C(i)不会被反映到所输出的二维投影像(三维)中。另外，Cout(i)、Aout(i)都将0作为初始值。

式(1)的S(i)为用于赋予阴影的加权成分，且利用根据亮度值C(i)和其周边的亮度值求出的、亮度值的倾斜度来计算出。例如，在以第i个体素为中心的面(亮度值的倾斜度)的法线与预先确定的光源的光轴一致的情况下，为了进一步强烈反射光，基于预先确定的表或函数，对体素i作为S(i)而赋予1.0，在光源和法线正交的情况下，作为S(i)而赋予0.0。由此，对获得的二维投影像施加阴影，获得增强效果。

多帧构成部(断层像)46根据在体数据生成部36中生成的三维体数据，生成任意截面的黑白断层像。任意截面的位置是界面部43从操作者接受并经由图像系统控制部44设定在多帧构成部(断层像)46中的。另外，能够设定多个任意截面位置，多帧构成部46对多个任意截面位置分别生成黑白断层像。

通过这样的结构，形成被检体1的诊断部位的三维像和任意截面的黑白断层像。

另一方面，在RF信号帧数据存储部27中依次存储调相加法部6生成的RF信号帧数据。例如，RF信号帧数据存储部27在帧存储器内依次存储由调相加法部6基于时间序列即图像的帧率而生成的RF信号数据。RF信号帧数据选择部28从存储在RF信号帧数据存储部27中的多个RF信号帧数据中选出一组、即两个RF信号帧数据。例如，RF信号帧数据选择部28根据来自图像系统控制部44的指令，将最新存储的RF信号帧数据(N)选作第一数据，并且从在时间上过去存储的RF信号帧数据组(N-1、N-2、N-3、…、N-M)中选择一个RF信号帧数据(X)。另外，这里，N、M、X为对RF信号帧数据赋予的索引号，且是自然数。

位移测量部30从一组RF信号帧数据中求出生物体组织的位移等。

例如，位移测量部30对由RF信号帧数据选择部28选择出的一组数据、即RF信号帧数据(N)以及RF信号帧数据(X)进行一维或者二维相关处理，求出与断层图像(二维反射率像)的各个点对应的生物体组织中的位移或与移动向量(位移的方向和大小)有关的一维或者二维位移分布。这里，在移动向量的检测中使用模块匹配法。模块匹配法是如下方法：将断层图像分割为例如由N×N像素构成的模块，关注关心区域内的模块，并从之前的帧开始寻找与关注模块的反射率分布最近似的模块，参照这个且通过预测编码即差分来决定标本值。

弹性信息运算部32基于位移测量部30求出的位移或移动向量，进行规定的运算，从而运算弹性值，并作为时间序列的弹性帧数据来输出。这里所指的弹性值只要是表示被检体1的组织的弹性的值即可，作为一例，可列举变形、弹性率、位移、粘性、变形比等。在作为弹性值而使用变形的情况下，能够通过对生物体组织的移动量例如位移进行空间微分来计算出。

二维弹性像构成部34包括帧存储器和图像处理部，通过将从弹性信息运算部32以时间序列输出的弹性帧数据存储在帧存储器中，并在图像处理部中处理所存储的帧数据，从而生成表示被检体的诊断部位中的弹性值的二维分布的二维弹性像。二维弹性像为基于预先确定的颜色变换表将弹性值的值变换为色调信息的彩色像。例如，与弹性值从规定的小值变化为大值相对应地，赋予从蓝(B)到绿(G)、红(R)以255灰度(1～255)依次变化的色调。最硬的部位的弹性值为1，最软的部位的弹性值为255。

二维弹性像存储部39使由弹性像构成部34生成的二维弹性像及其取得位置建立对应后进行存储。体数据生成部40基于存储在二维弹性像存储部39中的二维弹性像及其取得位置来进行三维坐标变换，从而生成对在空间上连续的多个二维弹性像进行三维配置的三维体数据。

弹性值变换部51为了获得能够明确掌握操作者关注的硬度的组织、例如弹性值小的硬的范围的组织的形状的三维像，变换构成三维体数据的弹性值。由此，生成操作者关注的硬度组织的弹性值大、软的组织的弹性值小的三维体数据。关于变换方法，将在后述的实施方式1～4中进行具体说明。

体绘制部42使用下式(4)～(6)，对通过弹性值变换部51变换的弹性体数据进行绘制。投影方向(视线方向)是图像系统控制部44经由界面43从操作者接受的。由此，由于能够获得操作者关注的硬度的组织(变换后的弹性值大的范围)的体素的三维弹性像，所以能够掌握位于软的组织内侧的硬的组织的三维形状。

Eout(i)＝Eout(i-1)+(1-Aout(i-1))·A(i)·E(i)·S(i)···式(4)

Aout(i)＝Aout(i-1)+(1-Aout(i-1))·A(i)···式(5)

A(i)＝Eopacity[E(i)]···式(6)

在式(4)中，Eout(i)为作为投影面的像素值而输出的值。E(i)为从二维投影面上的某一点看三维弹性图像时的在视线上第i个(其中，i＝0～N-1)存在的体素的变换后的弹性值(通过弹性值变换部51变换后的弹性值)。在视线上排列了N个体素的弹性值时，将i＝0～N-1为止的弹性值根据式(4)进行了积算的积算值Eout(N-1)成为最终输出的像素值。Eout(i-1)表示第i-1个为止的积算值。

式(4)的A(i)为在视线上第i个存在的体素的不透明度，为0～1.0的值。如式(6)所示，通过参照确定了预先决定的变换后的弹性值E(i)和不透明度(Opacity)的关系的表(Eopacity[E(i)])，或者通过向确定了预先决定的弹性值E(i)和不透明度的关系的函数(Eopacity[E(i)])代入弹性值E(i)，从而根据体素的弹性值的大小来确定不透明度A(i)。

式(5)的Aout(i)是，将通过式(6)赋予的不透明度A(i)直到第i个体素为止根据式(5)的右边进行了积算的值。在式(4)中，使用如式(5)所示那样计算出的第i-1个体素为止的不透明度的积算值Aout(i-1)。从式(5)可知，Aout(i)在每次通过体素时被积算且收敛为1.0。因此，如上述(4)所示那样第i-1个为止的不透明度的积算值Aout(i-1)约为1.0的情况下，式(4)的右边第2项成为0，第i个以后的弹性值E(i)不会被反映到所输出的二维投影像(三维)中。另外，Eout(i)、Aout(i)都将0作为初始值。

式(4)的S(i)为用于赋予阴影的加权成分，且利用根据变换后的弹性值(i)及其周边的变换后的弹性值求出的、弹性值的倾斜度来计算出。例如，在以第i个体素为中心的面(弹性值的倾斜度)的法线与预先确定的光源的光轴一致的情况下，为了进一步强烈反射光，基于预先确定的表或函数，对体素i作为S(i)而赋予1.0，在光源和法线正交的情况下，作为S(i)而赋予0.0。由此，对获得的二维投影像施加阴影，给出增强效果。另外，以第i个体素为中心的面的倾斜度(法线方向)是根据第i个体素及其周边的体素的弹性值或位置求出的。

将在后面叙述对于三维弹性像的色调的赋予方法。

多帧构成部(弹性像)48根据弹性值变换部51进行变换之前的弹性体数据，生成任意截面中的二维弹性像。对于任意截面的指定而言，界面部43从操作者接受，并经由图像系统控制部44传递给多帧构成部48。另外，能够设定多个任意截面位置，多帧构成部48对多个任意截面位置分别生成二维弹性像。

切换加法部12是具备帧存储器、图像处理部、图像选择部的结构。帧存储器分别保存弹性像构成部7及多帧构成部46生成的黑白断层像、二维弹性像构成部34及多帧构成部48生成的彩色二维弹性像、体绘制部38生成的三维像、以及体绘制部48生成的三维弹性像。

切换加法部12根据操作者的指示，能够生成在黑白断层像上以规定的比例相加彩色二维弹性像而得到的合成图像。作为黑白断层像，使用断层像构成部7生成的黑白断层像或者多帧构成部46生成的任意截面的黑白断层像。作为彩色二维弹性像，使用二维弹性像构成部34生成的二维弹性像、或者多帧构成部48生成的任意截面的二维弹性像。关于合成方法，使用专利文献1等记载的公知的合成图像的生成方法，能够使用在黑白断层像的各个像素的输出值上相加对与彩色二维弹性值对应的各个像素的输出值乘以规定的比例的值的方法。生成的合成图像保存在帧存储器中。

切换加法部12根据经由界面部43接受到的来自操作者的指示，在帧存储器中保存的黑白断层像、彩色二维弹性像、三维像、三维弹性像以及合成图像中，选择显示在图像显示部13中的图像并进行传递。图像显示部13以规定的排列在画面上显示传递到的一个以上的图像。

图像系统控制部44控制与图像生成相关的各个部分。界面部43还从操作者接受弹性图像的色调(色彩图的色调)、ROI(关心区域)的设定、帧率等的设定。图像系统控制部44在图像显示部13的画面上的规定位置上还显示为了生成三维弹性像而设定的弹性值的范围、关于弹性图像的色彩图、所设定的各种参数的值等。

<实施方式1>

使用图2～图4说明实施方式1的弹性值变换部51的动作。

弹性值变换部51通过读入并执行内置的存储器所存储的程序，从而如图2所示的流程图那样进行弹性值变换处理。

首先，弹性值变换部51在图2的步骤61中，获取体数据生成部36所生成的弹性值的体数据。该变换前的体数据成为如下结构：对弹性值最大的软的体素分配相对值255、对弹性值小的硬的体素分配相对值1，从而相对值与弹性值成比例且逐渐变化。即，弹性值的相对值整体上为255个灰度。另外，在因为噪声或错误而没有弹性信息的情况下，对弹性数据分配相对值0，表示没有弹性信息。

在步骤62中，弹性值变换部51对构成所获取的体数据的各个体素的弹性值(变换前的相对值)及其个数进行计数，并生成如图3(a)所示那样表示弹性值(变换前的相对值)和频度的关系的曲线图100。在图3(a)中，曲线图100的线101表示在体数据中包含的最硬的弹性值(变换前的相对值)。线102表示体数据所包含的最软的弹性值(变换前的相对值)。线103表示体数据所包含的弹性值(变换前的相对值)中频度最高的值(众数P)。

在步骤63中，如图3(b)所示，弹性值变换部51进行如下变换处理：对以众数P为中心的规定范围的弹性值(变换前的相对值)，作为变换后的相对值而分配最大值255，对其他范围的弹性值(变换前的相对值)，作为变换后的相对值而分配最小值。另外，作为最小值，既可以分配作为弹性值的最小值的1，也可以分配表示没有弹性信息的0，但在此，以下说明作为最小值而分配0的例子。

具体地说，如图3(a)所示，弹性值变换部51计算出比众数P小预先确定的值k的值P-k的弹性值(变换前相对值)104、和比众数P大k’的值P+k’的弹性值(变换前相对值)105。如图3(b)所示，弹性值变换部51进行如下变换：对P-k的弹性值(变换前相对值)104至P+k’的弹性值(变换前相对值)105的范围，作为变换后的相对值而分配最大值255，对这以外的范围，作为变换后的相对值而分配最小值0。例如，设为k＝10、k’＝10。

由此，生成向与出现频度高的规定硬度的弹性值对应的体素赋予了大的值(最大值255)、向这以外的频度低的弹性值的体素赋予了小的值(最小值0)的体数据。在步骤64中，弹性值变换部51将变换后的体数据传递给体绘制部42。

体绘制部42使用上述的式(4)～(6)对变换后的体数据进行绘制，通过求出Eout(i)而获得三维弹性像。获得的三维弹性像参照图4(b)的三维弹性像用色彩图15而进行着色。

三维弹性像用色彩图15是，向弹性值变换部51赋予了最大值255的弹性值范围(以众数P为中心的规定范围的弹性值)的变换前的弹性值范围赋予对应的色调(一色)，且根据Eout(i)的值的大小来改变亮度(明度)的图。通过式(4)设定为：由于用于附加阴影的加权成分S(i)的贡献越大(S(i)小)则Eout(i)越小，所以接近黑色(低亮度)。通过设置成式(4)的S(i)的贡献越大、即Eout(i)的值越小则越接近黑色(低亮度)，从而获得阴影的效果，成为立体的三维弹性像。

对于与变换前的弹性值范围对应的色调(一色)而言，例如，若以众数P为中心的规定范围的弹性值的、变换前的弹性值的平均值、最大值或者最小值中预先确定的值小于弹性值整体的平均值127，则设为蓝色，若大于127，则设为红色。或者，也可以是将以众数P为中心的规定范围的弹性值的、变换前的弹性值的平均值、最大值或者最小值中预先确定的值作为在二维弹性像的色彩图14中对应的色调。

在本实施方式1中，由于软的组织的相对值为0，所以在三维弹性像中未描绘弹性值变换部51变换后的体数据，能够获得操作者期望的频度高的硬的组织(变换后的弹性值大的范围)的三维弹性像。并且，由于频度高的硬的组织的体素的变换后的相对值大，所以积算后的数值大，能够获得对比度高的清晰的三维像。

经由切换加法部12，在图像显示部13中显示获得的三维像。图4(a)为在图像显示部13中显示的图像的一例。在画面10的右侧区域，显示有三维弹性像9。三维弹性像9是变换为大的相对值的频度高的硬的组织的像。在三维弹性像9的右侧，显示有表示三维弹性像9的颜色和弹性值的关系的三维弹性像用色彩图15。

在画面10的左侧区域，在此显示有黑白断层像11和二维弹性像8的合成图像。由于仅对操作者设定的ROI生成二维弹性像8，所以仅在黑白断层像11的显示区域的中央部的区域中进行合成。在二维弹性像8的左侧，显示有表示二维弹性像8的色调和弹性值的关系的色彩图14。由于根据弹性值变换部51进行变换之前的数据来生成二维弹性像8，所以色彩图14分别被分配将弹性值大的软的区域设为红色调、将弹性值小的硬的区域设为蓝色调、将中间的区域设为绿色调，且各个颜色随着弹性值的值而阶段性地产生变化，整体为256个灰度。另外，在画面10的下部的区域16中显示有从操作者接受到的各种参数的值。

在本实施方式1中，包括：弹性像体数据生成部40，使用向被检体1内发送超声波而接收到的信号，生成由表示弹性的弹性值数据构成的体数据；弹性值变换部51，对于构成体数据的弹性值数据，按弹性值的每个大小对频度进行计数，并进行使基于频度确定的规定范围的弹性值的值变大、使规定范围外的弹性值的值变小的变换；以及体绘制部42，通过对弹性值变换部51进行变换后的体数据进行绘制，从而生成三维弹性像。

弹性值变换部51通过将基于频度确定的规定范围的弹性值变换为第一常数、将规定范围外的弹性值变换为比第一常数小的第二常数，从而进行二值化。可将上述的基于频度确定的规定范围的弹性值设为包含频度最高的弹性值的弹性值范围。

根据图4(a)可知，能够清晰地描绘出被软的组织包围的频度高的硬度组织的三维弹性像9。因此，由于操作者通过确认黑白断层像11和二维弹性像8的合成图像，从而能够掌握在软的组织中存在硬的组织的情况，能够通过同时观察三维弹性像9来掌握频度高的硬度组织的三维形状，所以能够视觉辨认操作者所期望的肿瘤等硬的组织的三维弹性像。

作为比较例，图4(c)表示将弹性值变换部51进行变换之前的弹性值体数据通过式(4)～(6)进行绘制而得到的三维弹性像19。由于在体数据中，硬的组织的弹性值(相对值)小，包围硬的组织的软的组织的弹性值(相对值)大，所以若通过绘制来对视线上的体素的弹性值依次进行积算，则弹性值大的软的体素对二维投影面的贡献率变高，如图4(c)所示，获得描绘出了软的组织的外观的三维像19。因此，在图4(c)中，不能视觉辨认存在于软的组织内侧的应关注的硬的组织的像。

另外，在本实施方式1中，以众数P为中心，对P-k到P+k’的范围赋予了最大值，但赋予最大值的范围也可以是操作者经由界面部而设定范围的结构。此外，也可以使用如方差那样的统计值。此外，也可以以众数P以外的、平均值、最大值、最小值等统计值为中心设定范围。也可以构成为：基于频度确定的规定范围的弹性值包括使用根据弹性值的每个大小的频度求出的、构成体数据的弹性值的众数、平均值、最大值以及最小值中的其中一个来确定的值。例如，基于频度确定的规定范围的弹性值是根据弹性值的每个大小的频度而求出的、构成体数据的弹性值的最小值，弹性值变换部51向该弹性值的最小值分配规定的大的值，且随着构成体数据的弹性值增大，逐渐分配比规定的大的值小的值。

此外，在本实施方式中，对规定范围赋予最大值255，对除此之外的范围赋予最小值0，从而进行了二值化，但所赋予的两个值也可以不是0、255，只要是充分存在差异的两个值即可。

此外，也能够将体数据的弹性值预先分离为多个区域，按分离之后的每个区域，通过本实施方式1进行弹性值的变换来进行绘制。

<实施方式2>

接着，使用图5和图6说明实施方式2的弹性值变换部51的动作。在实施方式2中，弹性值变换部51进行如下变换：向弹性值中的众数P分配最大的相对值，向以众数P为中心的其以上的弹性值范围以及其以下的弹性值范围，以相对值线性减小的方式分配相对值。即，弹性值变换部51将基于频度确定的规定范围的弹性值变换为第一常数，将规定范围外的弹性值变换为随着弹性值远离规定范围的弹性值而逐渐成为比第一常数小的值。

首先，弹性值变换部51进行实施方式1的图2的步骤61、62，针对从弹性像体数据生成部40输出的体数据，生成如图5(a)所示那样表示弹性值(变换前的相对值)和频度的关系的曲线图200。曲线图200所示的线201表示包含在体数据中的最硬的弹性值(变换前的相对值)。线202表示包含在体数据中的最软的弹性值(变换前的相对值)。线203表示包含在体数据中的弹性值(变换前的相对值)中频度最高的值(众数P)。

在实施方式2中，在步骤63中进行如图5(b)所示那样变换相对值的处理。即，弹性值变换部51向线203所表示的众数P，作为变换后的相对值而分配255，对线201所表示的最硬的弹性值(变换前的相对值)分配0。向线203至线201之间的弹性值，分配以随着变换前的弹性值(相对值)减小而变换后的相对值逐渐减小的方式变换后的相对值。同样地，对线203至线202之间的弹性值，分配以随着变换前的弹性值(相对值)增大而变换后的相对值逐渐减小的方式变换后的相对值。

由此，向与体数据的出现频度最高的硬的弹性值对应的体素赋予最大的相对值255，向以众数为中心的除此之外的弹性值的体素赋予随着远离众数的弹性值而减小的相对值，弹性值的最小值和最大值成为最小的相对值0。

在步骤64中，弹性值变换部51将变换后的体数据传递给体绘制部42。体绘制部42使用上述的式(4)～(6)对变换后的体数据进行绘制，获得三维弹性像。

图6表示获得的三维弹性像的一例。在本实施方式2中，由于向弹性值的频度高的组织分配大的相对值，以众数为中心而相对值逐渐变化的方式进行变换，所以能够以大的对比度清晰地描绘出操作者所期望描绘出的频度高的硬的组织(变换后的弹性值大的范围)的三维弹性像400。并且，由于频度低的软的组织的变换后的相对值也不是0，所以通过绘制，如图6所示那样在硬的组织的三维弹性像400的周围描绘出频度低的软的组织的三维弹性像402。由于频度低的软的组织的相对值小，所以基于绘制的积算后的数值也小，频度低的软的组织的三维弹性像402成为对比度低的像。因此，如图6所示那样，频度低的软的组织的三维弹性像402不会覆盖隐藏频度高的硬的组织的三维弹性像400，能够确认硬的组织的三维形状，且还能够掌握软的组织的三维形状。

由此，根据实施方式2，能够掌握操作者所期望的肿瘤等硬的组织的三维弹性像，还能够确认软的组织的三维弹性像。

另外，在实施方式2中，按照以众数为中心而变换后的相对值逐渐减小的方式进行了弹性值的变换，但中心值并不限于众数，还可以是以平均值、最大值、最小值等统计值为基准来计算出中心值的结构。

在实施方式2中，将变换后的相对值的最大值设为255，但也可以是比255小的规定值。

此外，如图5所示，在实施方式2中，变换后的相对值被变换成以中心值为中心而线性(直线)减小，但不需要一定是线性(直线)变化，也可以是曲线或阶段状等变换后的相对值逐渐减小的变化。

此外，也能够将体数据的弹性值预先分离为多个区域，对分离之后的每个区域，通过本实施方式2进行弹性值的变换来进行绘制。

<实施方式3>

接着，使用图7说明实施方式3的弹性值变换部51的动作。在实施方式3中，弹性值变换部51进行如下的变换：向构成体数据的弹性值(变换前相对值)中的最小的弹性值(最硬的组织)分配最大的相对值255，向最大的弹性值(最软的组织)分配最小的相对值0，向其之间的弹性值范围以相对值根据弹性值而线性减小的方式分配相对值。

弹性值变换部51进行实施方式1的图2的步骤61、62，对变换前的体数据，如图7(a)所示那样生成表示弹性值(变换前的相对值)和频度的关系的曲线图300。线301表示包含在弹性体数据中的最小(硬)的弹性值。线302表示包含在弹性体数据中的最大(软)的弹性值。

在步骤63中，弹性值变换部51提取线301、302的弹性值(变换前的相对值)，且如图7(b)所示，向最小(硬)的弹性值(线301)，作为变换后的相对值而赋予最大值255。向最大(软)的弹性值(线302)，作为变换后的相对值而赋予最小值0。线301和线302之间的弹性值线性变化。其他结构与实施方式1相同。

由此，由于能够将硬的组织变换为大的弹性值(变换后的相对值)、将软的组织变换为小的弹性值(变换后的相对值)，且能够向最硬的组织赋予最大值255、向最软的组织赋予最小值0，所以能够最大限度地发挥相对值的灰度。

在本实施方式3中获得的三维弹性像中，由于对操作者期望描绘出的肿瘤等硬的组织赋予了大的弹性值(变换后)，所以能够以大的对比度清晰地描绘出硬的组织的三维弹性像。并且，由于软的组织的变换后的相对值也不是0，所以虽然贡献度低但能够描绘出三维弹性像。因此，能够确认硬的组织的三维形状，并且还能够掌握软的组织的三维形状。

另外，也能够将体数据的弹性值预先分离为多个区域，对分离之后的每个区域，通过本实施方式3进行弹性值的变换来进行绘制。

<实施方式4>

接着，使用图8说明实施方式4的弹性值变换部51的动作。在实施方式4中，通过对在构成体数据的弹性值(变换前相对值)中规定的小的弹性值(硬的组织)乘以权重w(w＞1)，从而进行增大弹性值的变换。

通过对大的弹性值(软的组织)乘以权重w(w＜1)，从而进行减小弹性值的变换。

弹性值变换部51进行实施方式1的图2的步骤61、62，对变换前的体数据，如图8(a)所示那样生成表示弹性值(变换前的相对值)和频度的关系的曲线图300。线401表示包含在弹性体数据中的最小(硬)的弹性值。线402表示包含在弹性体数据中的最大(软)的弹性值。

在步骤63中，弹性值变换部51提取线401、402的弹性值(变换前的相对值)，并且求出构成体数据的全部弹性值的平均值。图8(b)的线403表示求出的平均值。

如图8(b)所示，弹性值变换部51生成弹性值的函数、即权重w。如下生成w：对于弹性值的平均值而言成为1，在从最小(硬)的弹性值(线401)到平均值(线403)的范围内成为w＞1，在从最大(软)的弹性值(线402)到平均(线403)的范围内成为w＜1。即，弹性值变换部51求出根据弹性值的每个大小的频度而求出的、构成体数据的弹性值的平均值、最大值以及最小值，且通过对平均值的弹性值乘以权重1、对从最小值到平均值为止的弹性值乘以大于1的权重、对从平均值到最大值为止的弹性值乘以小于1的权重，从而改变弹性值。

例如，能够将w确定为在对最小值(最硬)的弹性值(线401)乘以了权重w的情况下成为255。此外，能够将最大值(最软)的弹性值(线402)的权重w设为0。这样，能够将连接平均值(线403)的权重w＝1的直线或者曲线设为预先确定的形状。

弹性值变换部51通过将生成的权重w的函数乘以构成体数据的弹性值(变换前)，从而生成变换后的体数据。由此，硬的组织的弹性值变大、软的组织的弹性值变小。此外，由于能够对最硬的组织赋予最大值255、对最软的组织赋予最小值0，所以能够最大限度地发挥相对值的灰度。

在本实施方式4中获得的三维弹性像中，由于对操作者期望描绘出的频度高的硬的组织赋予大的弹性值(变换后)，所以能够以大的对比度清晰地描绘出频度高的硬的组织的三维弹性像。并且，由于频度低的软的组织的变换后的相对值也不是0，所以虽然贡献度低但能够描绘出三维弹性像。因此，能够确认硬的组织的三维形状，并且还能够掌握软的组织的三维形状。

另外，也能够将体数据的弹性值预先分离为多个区域，对分离之后的每个区域，体绘制部42通过本实施方式3进行弹性值的变换而进行绘制。

另外，在上述的实施方式1～4中，说明了变换弹性体数据的整体弹性值的情况，但也可以仅对在弹性体数据中的一部分区域进行弹性值变换。

此外，在本实施方式1～4中，说明了在生成体数据之后改变弹性值的例子，但也可以对生成体数据之前的二维弹性像构成部34的输出实施本实施方式的弹性值的变换，之后生成弹性体数据。

符号说明：

1被检体、2探头、3发送部、4接收部、5超声波收发控制部、6调相加法部、7断层图像构成部、8彩色二维弹性像、9三维弹性像、10画面、11黑白断层像、12切换加法部、13图像显示部、14色彩图(二维弹性像)、15色彩图(三维弹性像)、16参数显示区域、27RF信号帧数据存储部、28RF信号帧数据选择部、30位移测量部、32弹性信息运算部、34二维弹性像构成部、36断层像的体数据生成部、38体绘制部、39二维弹性像存储部、40弹性像的体数据生成部、42体绘制部、43界面部、44图像系统控制部、46短轴扫描位置控制部、46多帧构成部(断层像)、48多帧构成部(弹性像)、51弹性值变换部

Claims (13)

1.一种超声波诊断装置，其特征在于，包括： 

弹性像体数据生成部，其使用向被检体内发送超声波而接收到的信号，生成由表示弹性的弹性值数据构成的体数据； 

弹性值变换部，其针对构成所述体数据的弹性值数据，按弹性值的每个大小对频度进行计数，并进行使基于所述频度确定的规定范围的弹性值的值变大、使所述规定范围外的弹性值的值变小的变换；以及 

体绘制部，其通过对所述弹性值变换部进行变换后的体数据进行绘制，从而生成三维弹性像。 

2.根据权利要求1所述的超声波诊断装置，其特征在于， 

所述弹性值变换部将基于所述频度确定的规定范围的弹性值变换为第一常数，将所述规定范围外的弹性值变换为比所述第一常数小的第二常数，从而进行二值化。 

3.根据权利要求1所述的超声波诊断装置，其特征在于， 

所述弹性值变换部将基于所述频度确定的规定范围的弹性值变换为第一常数，将所述规定范围外的弹性值变换为随着弹性值远离所述规定范围的弹性值而逐渐成为比第一常数小的值。 

4.根据权利要求1所述的超声波诊断装置，其特征在于， 

基于所述频度确定的规定范围的弹性值是包括频度最高的弹性值的弹性值范围。 

5.根据权利要求1所述的超声波诊断装置，其特征在于， 

基于所述频度确定的规定范围的弹性值包括使用根据所述弹性值的每个大小的频度而求出的、构成所述体数据的弹性值的众数、平均值、最大值以及最小值中的其中一个而确定的值。 

6.根据权利要求1所述的超声波诊断装置，其特征在于， 

基于所述频度确定的规定范围的弹性值是根据所述弹性值的每个大小的频度而求出的、构成所述体数据的弹性值的最小值，所述弹性值变换部向该弹性值的最小值分配规定的大的值，且随着构成所述体数据的弹性值增大而分配比所述规定的大的值逐渐变小的值。 

7.根据权利要求1所述的超声波诊断装置，其特征在于， 

所述弹性值变换部求出根据所述弹性值的每个大小的频度而求出的、构成所述体数据的弹性值的平均值、最大值以及最小值，且通过对所述平均值的弹性值乘以权重1、对从所述最小值到平均值为止的弹性值乘以大于1的权重、对从所述平均值到最大值为止的弹性值乘以小于1的权重，从而改变弹性值。 

8.根据权利要求1所述的超声波诊断装置，其特征在于， 

所述弹性值变换部进行变换后的体数据的最大值以及最小值分别为体绘制部能够绘制的最大值以及最小值。 

9.根据权利要求1所述的超声波诊断装置，其特征在于， 

将所述体数据的弹性值预先分离为多个区域，所述体绘制部对分离后的每个区域进行所述弹性值的变换，从而进行绘制。 

10.根据权利要求1所述的超声波诊断装置，其特征在于， 

所述弹性值变换部通过对构成所述体数据的弹性值中的规定的小的弹性值乘以权重，从而进行使弹性值增大的变换。 

11.根据权利要求1所述的超声波诊断装置，其特征在于， 

所述弹性值变换部进行对所述弹性值中的众数分配最大的相对值、对众数以上的弹性值范围以及众数以下的弹性值范围以相对值线性减小的方式分配相对值的变换。 

12.根据权利要求1所述的超声波诊断装置，其特征在于， 

所述弹性值变换部进行对构成所述体数据的弹性值中的最小的弹性值分配最大的相对值、对最大的弹性值分配最小的相对值、对其间的弹性值范围以相对值根据弹性值而线性减小的方式分配相对值的变换。 

13.一种超声波图像的显示方法，其特征在于， 

使用向被检体内发送超声波而接收到的信号，生成由表示弹性的弹性值数据构成的体数据， 

对于构成所述体数据的弹性值数据，按弹性值的每个大小对频度进行计数，并进行使基于所述频度确定的规定范围的弹性值的值变大、使所述规定范围外的弹性值的值变小的变换， 

通过对所述变换后的体数据进行绘制，从而生成三维弹性像来进行显示。 

Images