超声波诊断装置和超声波诊断方法与流程

本申请是申请号为201480018632.9、申请日为2014年3月27日、发明名称为“超声波诊断装置和超声波诊断方法”的发明专利申请的分案申请。

本发明涉及超声波诊断装置和超声波诊断方法，并且更特别地涉及使用声辐射压力的超声波诊断装置和超声波诊断方法。

背景技术：

在诸如癌症的疾病中，已知疾病的状况与患部的硬度相关。已知活体内的硬度分布在其诊断中是很重要的信息。一种定量获得组织的硬度的方法是根据横波传播的速度来获得弹性分布的方法。在这种方法中，在活体内产生横波并且使用超声波来测量该横波，由此估计传播速度。

最近已经提出使用通过发射超声波来产生的声辐射压力来挤压活体中的组织的方法作为生成横波的方法(例如，参见npl1)。

在npl1中所述的方法中，用于生成横波的大约数百微秒的超声波突发波(推压超声波(pushultrasonicwave))被发射，并且使用声辐射压力来施加压力(在下文中还可以称为振动)。其后，为了观察目标对象内生成的横波的传播，以大约几khz的重复频率来执行超声波图像测量(跟踪超声波(trackultrasonicwave))。使用跟踪超声波的观察一般执行大约几毫秒至数十毫秒的时长。根据获得的图像来测量位移分布并且估计横波传播的速度。

引文列表

非专利文献

npl1：j.bercoff,m.tanter和m.fink，“supersonicshearimaging:anewtechniqueforsofttissueelasticitymapping”，ieeetrans.ultrason.,ferroelectr.,freq.control,pp.396-409,vol.51,no.14,2004/04。

技术实现要素：

技术问题

但是，与一般的诊断成像相比，该方法需要在短时间段内施加很强的能量。当这种超声波被发射时，会出现诸如暴露于超声波的部分的温度升高的缺陷，并因此需要改进。

因此，本发明的目的是提供能够在不导致暴露于超声波的部分的温度升高的情况下使用声辐射压力来进行超声波诊断的超声波诊断装置和超声波诊断方法。

问题的解决方案

为了解决上述问题，根据本发明的超声波诊断装置包含输出编码的加压脉冲信号的第一脉冲输出单元；

输出用于测量的测量脉冲信号的第二脉冲输出单元；

探测器，输出基于加压脉冲的用于在目标对象中产生横波的超声波和基于测量脉冲信号的用于测量的超声波；

接收器，接收用于测量的超声波的回波并输出电信号；以及

弹性模量估计器，对由接收器输出的电信号进行解码并且基于解码的信号来估计目标对象的弹性模量。

由于加压脉冲信号是用于加压的信号，即，由于加压脉冲信号只可以“推压”而无法“拉拽”，因而编码的加压脉冲信号使用单极性码。优选地，编码的加压脉冲信号基于双极性码而生成。优选地，超声波诊断装置还包含位移估计器。该位移估计器通过使用双极性码对由接收器输出的电信号进行解码来估计位移，并且弹性模量估计器基于所估计的位移来估计弹性模量。

优选地，第一脉冲输出单元在编码的加压脉冲信号的位为第一值的情况下输出加压脉冲信号，而在编码的加压脉冲信号的该位为第二值的情况下不输出加压脉冲信号。

编码的加压脉冲信号优选地为单极性码对，并且该单极性码对由双极性码生成。该单极性码对包含第一单极性码和第二单极性码。第一单极性码的位在双极性码的位为“+1”的情况下为“1”，而在双极性码的该位为“-1”的情况下为“0”。第二单极性码的位在双极性码的该位为“+1”的情况下为“0”，而在双极性码的该位为“-1”的情况下为“1”。

优选地，双极性码为巴克码(barkercode)。

优选地，在第一时段中，在双极性码的位数为n并且位移被估计的时间点的数量为m的情况下，第一脉冲输出单元在每个特定周期基于第一单极性码而输出数量为0个以上到n个以下的加压脉冲信号，并且第二脉冲输出单元每个周期在不同于加压脉冲信号的时刻输出测量脉冲信号(n+m)次。在第二时段中，第一脉冲输出单元每个周期基于第二单极性码输出数量为0个以上到n个以下的加压脉冲信号，并且第二脉冲输出单元每个周期在不同于加压脉冲信号的时刻输出测量脉冲信号(n+m)次。

优选地，第一脉冲输出单元在相邻的测量脉冲信号之间输出加压脉冲信号。

优选地，超声波诊断装置包含位移估计器，并且位移估计器包含：第一计算器，该第一计算器基于第一时段中的相邻的时间点的超声波回波来计算表示编码的单极性位移的第一单极性码位移信号，并且基于第二时段中的相邻的时间点的超声波回波来计算表示编码的单极性位移的第二单极性码位移信号；第二计算器，该第二计算器计算表示作为第一单极性码位移信号与第二单极性码位移信号之差的编码的双极性位移的双极性码位移信号；以及第三计算器，该第三计算器通过使用双极性码位移信号与双极性码的相关处理执行解码来计算表示位移的信号。

根据本发明的超声波诊断装置控制方法包括：输出编码的加压脉冲信号的步骤；输出基于加压脉冲信号的用于在目标对象中产生横波的超声波的步骤；输出用于测量的测量脉冲信号的步骤；输出基于该测量脉冲信号的用于测量的超声波的步骤；接收用于测量的超声波的回波信号并输出电信号的步骤；以及对该电信号进行解码并基于解码的电信号来估计弹性模量的步骤。

优选地，超声波诊断装置控制方法还包括：估计目标对象的位移的步骤，其中估计位移的步骤通过对该电信号进行解码来估计位移，以及估计弹性模量的步骤包括基于估计的位移来估计弹性模量。优选地，编码的加压脉冲信号包括以单极性码对而被编码。该单极性码对基于双极性码而生成。估计位移的步骤包括通过使用该双极性码解码基于接收的超声波回波而输出的电信号来估计位移的步骤。

优选地，该单极性码对由双极性码生成。该单极性码对包含第一单极性码和第二单极性码。在双极性码的位数为n并且位移被估计的时间点的数量为m的情况下，输出加压脉冲信号的步骤包括在第一时段中于每个特定周期基于第一单极性码而输出加压脉冲信号的步骤，该加压脉冲信号的数量为0个以上到n个以下，并且包括在第二时段中于每个周期基于第二单极性码而输出加压脉冲信号的步骤，该加压脉冲信号的数量为0个以上到n个以下。输出测量脉冲信号的步骤包括：在第一时段中于每个周期在不同于加压脉冲信号的时刻输出测量脉冲信号(n+m)次的步骤，以及在第二时段中于每个周期在不同于加压脉冲信号的时刻输出测量脉冲信号(n+m)次的步骤。

优选地，输出加压脉冲信号的步骤包括在相邻的测量脉冲信号之间输出加压脉冲信号的步骤。

发明的效果

按照根据本发明的超声波诊断装置和超声波诊断方法，可以在不导致暴露于超声波的部分的温度升高的情况下执行使用声辐射压力进行的超声波诊断。

附图说明

[图1]图1是表示根据本发明的实施例的超声波诊断装置的配置的图。

[图2]图2是表示根据本发明的实施例超声波诊断装置的操作过程的流程图。

[图3]图3是表示图2中的步骤s302中的处理过程的流程图。

[图4]图4是表示图2中的步骤s303中的处理过程的流程图。

[图5]图5是用于描述基于现有技术的推压脉冲和跟踪脉冲的测量的图。

[图6]图6是用于描述基于根据本发明的实施例的推压脉冲和跟踪脉冲的测量的图。

具体实施方式

在下文中，将利用附图来描述本发明的实施例。

图1是表示根据本发明的实施例的超声波诊断装置的配置的示图。

参照图1，超声波诊断装置50包含超声波探测器3和超声波诊断装置主体4。

超声波探测器3包含一维排列的换能器1-1至1-n，以及一维排列的换能器2-1至2-n。

换能器1-1至1-n输出作为用于测量的超声波的超声波平面波，并且接收回波(用于测量的超声波的反射波)。

换能器2-1至2-n输出用于在作为目标对象的诊断目标20中产生横波的超声波，由此向该诊断目标20施加声辐射压力。

超声波诊断装置主体4包含：将作为加压脉冲信号的推压脉冲信号发送到探测器的推压脉冲发送器5、单极性码存储器6、将作为测量脉冲信号的跟踪脉冲信号发送到探测器的跟踪脉冲发送器7、回波接收器(接收器)8、位移估计器90、传播速度估计器12、弹性模量估计器13、显示器14以及双极性码存储器15。位移估计器90包含单极性位移估计器9、双极性位移估计器10和解码器11。

双极性码存储器15存储n位巴克码a[i](i＝0至n-1)。在此，a[i]表示巴克码的第i位。

单极性码存储器6存储表示正部的单极性码ap[i](i＝0至n-1)以及表示负部的an[i](i＝0至n-1)，该单极性码构成了基于巴克码a[i]生成的一个组。使用n位巴克码a[i]，以下列等式(1)和(2)表示ap[n]和an[n]。

当a[i]＝+1时，ap[i]＝1；并且

当a[i]＝-1时，ap[i]＝0…(1)

当a[i]＝+1时，an[i]＝0；并且

当a[i]＝-1时，an[i]＝1…(2)

在n＝5的情况下，{a[i]}＝{+1,+1,+1,-1,+1}，{ap[i]}＝{1,1,1,0,1}，并且{an[i]}＝{0,0,0,1,0}。

推压脉冲发送器5使用单极性码ap和an来生成推压脉冲。因为使用声辐射压力来进行的压力施加只可以“推压”而无法“拉拽”，所以使用单极性码。在第一时段中，推压脉冲发送器5依据ap[0]至ap[n-1]按照特定的间隔向换能器2-i发送第一组编码的推压脉冲，由此驱动换能器2-i。在第二时段中，推压脉冲发送器5依据an[0]至an[n-1]按照特定的间隔向换能器2-i发送第二组编码的推压脉冲，由此驱动换能器2-i。

在第一时段和第二时段中，跟踪脉冲发送器7按照特定的间隔向换能器1-i发送跟踪脉冲，由此驱动换能器1-i。

在第一时段和第二时段中，回波接收器8接收基于从换能器1-i输出的超声波回波(反射波)的信号，执行诸如放大之类的处理，并且输出电信号(在下文中还可以被称为回波电信号)。

位移估计器90包含单极性位移估计器9、双极性位移估计器10和解码器11。位移估计器90通过使用双极性码对电信号进行解码来估计相对一个较早时间点的位移。

根据相邻时间点(帧)中的电信号，单极性位移估计器9通过使用在超声波测量中使用的位移测量方法来计算表示编码的单极性位移的单极性码位移信号。

根据第一时段中的单极性码位移信号以及第二时段中的单极性码位移信号，双极性位移估计器10计算表示编码的双极性位移的双极性码位移信号。

解码器11通过使用巴克码a与双极性码位移信号的相关处理来执行解码(脉冲压缩)从而计算表示位移的双极性位移的信号。在不执行诸如在npl1中描述的编码之类的编码的情况下，可以获得与在大振幅推压超声波于很短的时间段内被施加的情况下获得的信号等同的信号。

基于双极性位移信号，传播速度估计器12估计在目标对象中产生的横波的传播速度。

弹性模量估计器13基于横波传播的速度计算弹性模量(杨氏模量(young'smodulus))。

显示器14显示具有作为像素值的弹性模量的二维图像。

(操作)

图2是表示根据本发明的实施例的超声波诊断装置的操作过程的流程图。

参照图2，在步骤s301中，变量x被设置为“dx/2”。变量x表示诊断目标20的x坐标，如图1所示。

在步骤s302中，推压脉冲发送器5基于单极性码ap[i](i＝0至n-1)输出推压脉冲信号。在接收到推压脉冲信号时，位于线x上的换能器2-i发射用于在线x上产生横波的超声波(推压超声波)。跟踪脉冲发送器7输出(n+m)个跟踪脉冲信号。在接收到跟踪脉冲信号时，换能器1-1至1-n发射用于测量的超声波平面波。在此，n是巴克码a以及单极性码ap和an的位数。m是待显示于显示器14上的图像的帧数，即，估计位移的时间点的数量。

换能器1-1至1-n接收用于测量的超声波的回波(反射波)。回波接收器8基于该超声波回波而生成电信号rp[0,y,xα]至rp[n+m-1,y,xα]。在此，y是基于超声波回波被接收到的时间t以及超声波的速度而计算的诊断目标20的y坐标。此外，y＝0至ys-1。ys是在待显示于显示器14上的二维图像的y方向(垂直方向)上的大小。在此，xα＝x-dx/2至x+dx/2-1。

在步骤s303中，推压脉冲发送器5基于单极性码an[i](i＝0至n-1)而输出推压脉冲信号。在接收到推压脉冲信号时，位于线x上的换能器2-i发射用于在线x上产生横波的超声波。跟踪脉冲发送器7输出(n+m)个跟踪脉冲信号。在接收到跟踪脉冲信号时，换能器1-1至1-n发射用于测量的超声波平面波。

换能器1-1至1-n接收用于测量的超声波的回波。回波接收器8基于超声波回波而生成电信号rn[0,y,xα]至rn[n+m-1,y,xα]。

在步骤s304中，在变量x不等于(xs-dx/2-1)的情况下，处理进行到步骤s305；并且，在变量x变为等于(xs-dx/2-1)的情况下，处理进行到步骤s306。注意，xs是在待显示于显示器14上的二维图像的x方向(水平方向)上的大小。

在步骤s305中，变量x仅增加“dx”，并且处理返回步骤s302。

在步骤s306中，变量k被设置为“0”。

在步骤s307中，单极性位移估计器9根据基于超声波回波的电信号rp[k,y,x]以及基于超声波回波的电信号rp[k+1,y,x]、通过使用超声波测量中使用的位移测量方法而计算单极性码位移信号dp[k,y,x]。

具体地，单极性位移估计器9如等式(1)所指出地在获得位移的点附近定义了尺寸为w的窗口，并且通过执行模板匹配来获得在互相关系数变为最大时的点ν作为位移。在此，△y是表示y方向上的采样间隔的距离。

[公式1]

在步骤s308中，单极性位移估计器9根据基于超声波回波的电信号rn[k,y,x]以及基于超声波回波的电信号rn[k+1,y,x]、通过使用超声波测量中使用的位移测量方法而计算出单极性码位移信号dn[k,y,x]。

具体地，单极性位移估计器9如等式(2)所指出地在获得位移的点附近定义了尺寸为w的窗口，并且通过执行模板匹配来获得在互相关系数变为最大时的点ν作为位移。在此，δy是表示y方向上的采样间隔的距离。

[公式2]

在步骤s309中，双极性位移估计器10通过从单极性码位移信号dp[k,y,x]中减去单极性码位移信号dn[k,y,x]来计算双极性码位移信号db[k,y,x]。

在步骤s310中，在变量k不等于(n+m-2)的情况下，处理进行到步骤s311；并且，在变量k变为等于(n+m-2)的情况下，处理进行到步骤s312。

在步骤s311中，变量k仅增加“1”，并且处理返回步骤s307。

在步骤s312中，变量k被设置为“0”。

在步骤s313中，解码器11针对i＝0至n-1计算巴克码a的第i位a[i]与双极性码位移信号db[k+i,y,x]的乘积，并且计算这些乘积的总和作为双极性位移信号d[k,y,x]。

在步骤s314中，在变量k不等于(m-1)的情况下，处理进行到步骤s315；并且在变量k变为等于(m-1)的情况下，处理进行到步骤s316。

在步骤s315中，变量k仅增加“1”，并且处理返回步骤s313。

在步骤s316中，传播速度估计器12基于双极性位移信号d[k,y,x]通过使用互相关等来估计传播时间，并然后估计横波传播的速度cs[y,x]。注意，y＝0至ys-1，且x＝0至xs-1。在此，ys是在待显示于显示器14上的二维图像的y方向(垂直方向)上的大小。

具体地，传播速度估计器12在获得弹性的点[y,x]附近按照横波传播的方向来定义两个点[y1,x1]和[y2,x2]，并且如下面的等式那样通过使用互相关来估计在这两个点之间的传播时间τ[y,x]。在此，△t是在时间点k与时间点k+1之间的间隔。

[公式3]

传播速度估计器12通过使用下面的等式根据横波传播时间τ[y,x]来获得横波传播速度cs[y,x]。在此，l是在已经获得了传播时间的两个点之间的距离。

[公式4]

在步骤s317中，弹性模量估计器13通过使用下面的等式来获得刚性模量g[y,x]。在此，ρ是密度。由于密度分布一般是未知的，因而使用平均密度。

g[y,x]＝ρ×cs[y,x]²…(4)

弹性模量估计器13按照下面的等式来计算杨氏模量e[y,x]。在此，ν是泊松比，在软生物组织中ν可以被认为基本上为0.5。

e[y,x]＝2×g[y,x]×(1+ν)…(5)

在步骤s318中，显示器14显示其中像素值用作e[y,x]的二维图像。

图3是表示图2中的步骤s302中的处理过程的流程图。

参照图3，在步骤s101中，变量i和j被设置为0。

在步骤s102中，跟踪脉冲发送器7发送跟踪脉冲信号。在接收到跟踪脉冲信号时，换能器1-1至1-n发射用于测量的平面超声波。

在步骤s103中，换能器1-1至1-n接收用于测量的超声波的回波。回波接收器8通过执行波束形成来生成基于超声波回波的电信号rp[i,y,xα]。如上所述，y＝0至ys-1。在此，xα＝x-dx/2至x+dx/2-1。

在步骤s104中，在单极性码ap[i]为“1”的情况下，处理进行到步骤s105；并且，在单极性码ap[i]为“0”的情况下，处理进行到步骤s106。

在步骤s105中，推压脉冲发送器5输出推压脉冲信号。在接收到推压脉冲信号时，位于线x上的换能器2-i发射用于在目标对象的线x上方的部分中产生横波的超声波。

在步骤s106中，推压脉冲发送器5不输出推压脉冲信号，并且位于线x上的换能器2-i不发射用于在目标对象的线x上方的部分中产生横波的超声波。

在步骤s107中，跟踪脉冲发送器7发送跟踪脉冲信号。在接收到跟踪脉冲信号时，换能器1-1至1-n发射用于测量的超声波平面波。

在步骤s108中，换能器1-1至1-n接收用于测量的超声波的回波。回波接收器8通过执行波束形成来生成基于超声波回波的电信号rp[i,y,xα]。如上所述，y＝0至ys-1。在此，xα＝x-dx/2至x+dx/2-1。

在步骤s109中，变量i仅增加“1”。

在步骤s110中，在变量i小于帧数n的情况下，处理返回步骤s104；并且在变量i变为等于帧数n的情况下，处理进行到步骤s111。

在步骤s111中，跟踪脉冲发送器7发送跟踪脉冲信号。在接收到跟踪脉冲信号时，换能器1-1至1-n发射用于测量的平面超声波。

在步骤s112中，换能器1-1至1-n接收用于测量的超声波的回波。回波接收器8通过执行波束形成来生成基于超声波回波的电信号rp[n+1+j,y,xα]。如上所述，y＝0至ys-1。在此，xα＝x-dx/2至x+dx/2-1。

在步骤s113中，变量j仅增加“1”。

在步骤s114中，在变量j小于(m-1)的情况下，处理返回步骤s111；并且在变量j变为等于(m-1)的情况下，处理结束。

作为上述操作的结果，获得基于超声波回波的电信号rp[0,y,xα]至rp[n+m-1,y,xα]。由于y＝0至ys-1且xα＝x-dx/2至x+dx/2-1，因而ys×dx×(n+m)个电信号rp被获得。

图4是表示图2中的步骤s303中的处理过程的流程图。

参照图4，在步骤s201中，变量i和j被设置为0。

在步骤s202中，跟踪脉冲发送器7发送跟踪脉冲信号。在接收到跟踪脉冲信号时，换能器1-1至1-n发射用于测量的超声波平面波。

在步骤s203中，换能器1-1至1-n接收用于测量的超声波的回波。回波接收器8通过执行波束形成来生成基于超声波回波的电信号rn[i,y,xα]。如上所述，y＝0至ys-1。在此，xα＝x-dx/2至x+dx/2-1。

在步骤s204中，在单极性码an[i]为“1”的情况下，处理进行到步骤s205；并且，在单极性码an[i]为“0”的情况下，处理进行到步骤s206。

在步骤s205中，推压脉冲发送器5输出推压脉冲信号。在接收到推压脉冲信号时，位于线x上的换能器2-i发射用于在诊断目标20的线x上方的部分中产生横波的超声波。

在步骤s206中，推压脉冲发送器5不输出推压脉冲信号，并且位于线x上的换能器2-i不发射用于在诊断目标20的线x上方的部分中产生横波的超声波。

在步骤s207中，跟踪脉冲发送器7发送跟踪脉冲信号。在接收到跟踪脉冲信号时，换能器1-1至1-n发射用于测量的超声波平面波。

在步骤s208中，换能器1-1至1-n接收用于测量的超声波的回波。回波接收器8通过执行波束形成来生成基于超声波回波的电信号rn[i,y,xα]。如上所述，y＝0至ys-1。在此，xα＝x-dx/2至x+dx/2-1。

在步骤s209中，变量i仅增加“1”。

在步骤s210中，在变量i小于帧数n的情况下，处理返回步骤s204；并且在变量i变为等于帧数n的情况下，处理进行到步骤s211。

在步骤s211中，跟踪脉冲发送器7发送跟踪脉冲信号。在接收到跟踪脉冲信号时，换能器1-1至1-n发射用于测量的超声波平面波。

在步骤s212中，换能器1-1至1-n接收用于测量的超声波的回波。回波接收器8通过执行波束形成来生成基于超声波回波的电信号rn[n+1+j,y,xα]。如上所述，y＝0至ys-1。在此，xα＝x-dx/2至x+dx/2-1。

在步骤s213中，变量j仅增加“1”。

在步骤s214中，在变量j小于(m-1)的情况下，处理返回到步骤s211；并且，在变量j变为等于(m-1)的情况下，处理结束。

作为上述操作的结果，获得基于用于测量的超声波回波的电信号rn[0,y,xα]至rn[n+m-1,y,xα]。由于y＝0至ys-1且xα＝x-dx/2至x+dx/2-1，因而ys×dx×(n+m)个电信号rn被获得。

(参见现有技术的推压脉冲信号和跟踪脉冲信号(npl1))。

图5是用于描述基于现有技术的推压脉冲信号和跟踪脉冲信号的测量的图。

在现有技术中，测量被多次重复以便获得特定的精度，并且使用作为测量结果的基于测量超声波(跟踪超声波)的回波的电信号的总体平均来执行步骤s316至s318中的处理。

在现有技术中，在第一测量中每周期tc1输出推压脉冲信号p[1,0]、p[1,1]、…、p[1,n-1]，由此发射用于在目标对象中产生横波的超声波。与本实施例中的ap和an不同，推压脉冲信号没有被编码并且每周期tc1都要输出。

在基于推压脉冲信号而发射用于在目标对象中产生横波的超声波之前，跟踪脉冲信号t[1,0]被输出，由此发射超声波平面波。另外，在基于推压脉冲信号而发射用于在目标对象中产生横波的超声波之后，每周期tc2输出跟踪脉冲信号t[1,1]、t[1,2]、…、t[1,m-1]，由此发射用于测量的超声波平面波。因此，基于用于测量的超声波的回波信号的电信号r[1,0]、r[1,1]、…、r[1,m-1]被获得。

在现有技术中，一次测量耗费时间tc1×n+tc2×m。如果完成s次测量来获得总体平均，则在现有技术中需要时间en1＝s×(tc1×n+tc2×m)。

(根据实施例的推压脉冲信号和跟踪脉冲信号)

图6是用于描述基于根据本发明的实施例的推压脉冲信号(加压脉冲信号)和跟踪脉冲信号(测量脉冲信号)的测量的图。

在该实施例中，每周期tc0输出跟踪脉冲信号tp[0]至tn[n+m-1]，由此发射出用于测量的超声波平面波。另外，在相邻的跟踪脉冲信号tp之间的时刻，基于单极性码ap[i](i＝0至n)每周期tc输出或者不输出推压脉冲信号(加压脉冲信号)pp[i]。因此，待输出的推压脉冲信号pp的数量大于或等于0且小于或等于n。在推压脉冲信号pp[i]被输出的情况下，用于在目标对象中产生横波的超声波基于该推压脉冲信号pp[i]被发射。

其后，每周期tc0输出跟踪脉冲信号tn[0]至tn[n+m-1]，由此发射出用于测量的超声波平面波。另外，在相邻的跟踪脉冲信号tn之间的时刻，基于单极性码an[i](i＝0至n)每周期tc输出或者不输出推压脉冲信号(加压脉冲信号)pn[i]。因此，待输出的推压脉冲信号pn的数量大于或等于0且小于或等于n。在推压脉冲信号pn[i]被输出的情况下，用于在目标对象中产生横波的超声波基于该推压脉冲信号pn[i]被发射。

因为推压脉冲信号已经被编码，所以可以使由响应于根据本实施例的推压脉冲信号而发送的超声波(推压超声波)在目标对象中产生的横波的振幅变得显著小于由现有技术的推压超声波产生的横波的振幅。另外，还可以使tc0变得比tc1长。因此，在该实施例中可以使每单位时间施加于诊断目标20的能量变得小于现有技术中每单位时间施加于诊断目标20的能量，由此防止温度升高。

在该实施例中，测量耗费时间en2＝2×tc0×(n+m)。与现有技术的时间en1(＝s×(tc1×n+tc2×m))的这个时间相比，tc0>tc1且tc0>tc2。根据由本发明的申请人所进行的实验，指示需要使现有技术中为获得总体平均而执行的测量的次数s相当大以便获得与本发明的实施例的测量精度等同的测量精度。因此，可以使该实施例中的测量时间en2整体上小于现有技术的测量时间en1。

在不使用编码的加压脉冲信号的情况下仅仅通过减小声辐射压力的振幅或者缩短超声波暴露时间来降低入射能量的情况下，对所产生的横波的观测中的snr(信噪比)降低，由此降低了估计传播波速的精度。相比之下，可以通过多次执行测量来获得总体平均从而提高snr。但是会存在测量时间随为获得均值而执行的测量的次数成比例增加的问题。根据本发明，与这样的技术(在不使用编码的脉冲信号的情况下仅仅通过减小声辐射压力的振幅或者缩短超声波暴露时间并且通过增加测量次数来提高snr的技术)相比，变得可以在不导致暴露于超声波的部分的温度升高的情况下在较短的时间段中使用声辐射压力进行超声波诊断。

尽管在该实施例中巴克码被用作双极性码的实例，但是双极性码并不限于此。还可以使用其他二进制码，诸如格雷码、m系列码(m-seriescode)和勒让德符号(legendresymbol)。

在格雷码的情况下，具体地使用两个格雷码a[i]和b[i](i＝0至n-1)。

超声波诊断装置通过使用与该实施例所描述的方法相同的或类似的方法使用{a[i]}来计算双极性位移信号d1[k,y,x](k＝0至m-1、x＝0至xs-1且y＝0至ys-1)。类似地，超声波诊断装置通过使用与该实施例所描述的方法相同的或类似的方法使用{b[i]}来计算双极性位移信号d2[k,y,x](k＝0至m-1、x＝0至xs-1且y＝0至ys-1)。超声波诊断装置获得极性位移信号d1[k,y,x]与双极性位移信号d2[k,y,x]之和作为目标双极性位移信号d[k,y,x]。

尽管在该实施例中利用使用等式(1)和(2)所示的互相关的模板匹配来计算位移，但是位移的计算并不限于此。例如，可以通过使用根据回波信号的正交检测而获得iq信号，可以根据相邻时间点之间的iq信号来估计rf信号的相位差，并且可以根据该相位差来计算位移。

尽管在该实施例中计算了相对一个较早时间点的位移，但是位移的计算并不限于此。使用参照初始值的位移，可以以类似的算法来执行解码。像相邻的帧之间的位移一样，相对初始值的累加位移通过横波传播。因此，可以使用参照初始值的位移来执行弹性估计。尽管在该实施例中推压脉冲信号的发送周期和跟踪脉冲信号的发送周期被描述为同一周期，但是这些发送周期并不限于此。跟踪脉冲信号的发送周期可以短于推压脉冲信号的发送周期。在这种情况下可以使用更多的跟踪脉冲来执行位移估计，由此进一步提高测量精度。而且，更优选的是，为了以高精度执行解码，跟踪脉冲信号的发送周期为推压脉冲信号的发送周期的整数分之一。

尽管在该实施例中跟踪脉冲信号被描述为以平面波的形式来发送，但是跟踪脉冲信号并不限于平面波。跟踪脉冲信号的发送范围可以被限定于包含用于执行弹性估计的区域的范围之内。例如，可以通过执行所谓的波束形成而形成发送波束来提高反射信号的信噪比，并且通过使用作为其结果而获得的具有高信噪比的信号来提高估计位移的精度。当估计位移的精度变高时，可以提供能够以甚至更高的精度来估计弹性的装置。

尽管在该实施例中通过执行位移估计并进一步计算传播速度来估计弹性模量，但是如果执行数学上等价的处理，则可以在不计算位移或传播速度本身的情况下根据所接收的信号直接估计出弹性模量。

尽管在该实施例中已经描述了对双极性码位移信号执行解码，但是待解码的信号并不限于此。例如，还可以通过对用于测量的超声波的回波信号进行解码或者对通过在不执行解码的情况下估计位移而获得的信号进行解码来实现本发明的有利效果。

本文所公开的实施例在所有方面都是示例性的，并且应当被理解为是非限制性的。本发明的范围并非由上述描述来指示，而是由权利要求书指示，并且意在包括具有与权利要求书的范围等价的含义的且属于权利要求书的范围的所有更改。例如，尽管上述描述讨论了目标对象的位移被估计的情况，但是基于估计的位移来估计横波传播的速度并然后估计弹性模量，其中仅位移估计或传播速度估计被执行且弹性模量估计不被执行的配置同样被包含于本发明的范围中。尽管用于发送推压超声波的元件以及用于发送/接收跟踪超声波的元件被分离地设置于探测器3中，但是可以用一个元件来执行这些操作。在这种情况下，存在可以使探测器变小的优点。另外，用于发送推压超声波的探测器以及用于发送/接收跟踪超声波的探测器可以被分离地设置。在这种情况下，这些探测器可以具有适于推压超声波的配置(形状)以及适于发送/接收跟踪超声波的配置(形状)，从而导致设计优化。

本申请基于在2013年3月28日在日本专利局提交的日本专利申请no.2013-069852并要求其优先权，其全文通过引用并入本文。

附图标记列表

1-1至1-n、2-1至2-n换能器

3超声波探测器

4超声波诊断装置主体

5推压脉冲发射器

7跟踪脉冲发射器

8回波接收器

9单极性位移估计器

10双极性位移估计器

11解码器

12传播速度估计器

13弹性模量估计器

14显示器

50超声波诊断装置

90位移估计器