具体实施方式
下面,根据本发明的实施方式进行说明。图1表示本发明的一个实施方式的超声波诊断装置的整体构成框图。如图所示,用来接触被测体1的超声波探头(以下省略为探头)2形成为,包括多个与被测体1之间发送和接收超声波的振子。探头2由发送电路3所提供的超声波脉冲驱动。收发控制电路4,控制驱动探头2的多个振子的超声波脉冲发送定时,对准设于被测体1内的焦点,形成超声波束。此外,收发控制电路4还使超声波束在探头2的振子的排列方向上进行电子式扫描。
另一方面,探头2接收被测体1内产生的反射回声信号,向接收电路5输出,接收电路5按照收发控制电路4输入的定时信号,取得反射回声信号,进行放大等接收处理。接收电路5处理的反射回声信号在调相加法电路6中得到放大,方法是使多个振子接收的反射回声信号相位同步,然后相加。由调相加法电路6进行过调相相加的反射回声信号被输入信号处理部7,进行增益修正、对数压缩、检波、轮廓强调、过滤处理等信号处理。
由信号处理部7处理的反射回声信号,被导入黑白扫描转换器8,转换成与超声波束的扫描面对应的二维断层图像数据(数字数据)。断层图像(B模式图像)的图像重组机构,由这些信号处理部7和黑白扫描转换器8构成。黑白扫描转换器8输出的断层图像数据通过切换加法部9,被输出至图像显示器10,使B模式图像得到显示。
另一方面,由调相加法电路6输出的反射回声信号被导入RF信号帧数据选择部11。RF信号帧数据选择部11,将与超声波束的扫描面对应的反射回声信号群作为帧数据,在取得多个帧份之后保存在存储器等中。位移运算部12,依次获得与RF信号帧数据选择部11所保存的取得时刻不同的多对帧数据,根据获得的一对帧数据,求出断层面上的多个测量点的位移矢量,使其作为位移帧数据输出至形变信息运算部13。
本实施方式的形变信息运算部13构成为:根据位移帧数据,求出各测量点的生物体组织的形变。使得形变信息运算部13所求出的形变分布(帧数据)输出至形变分布修正部14。
形变分布修正部14,通过形变分布修正函数制作部18输出的形变分布修正函数,对从形变信息运算部13输入的形变分布进行修正。然后,对基于修正后的形变分布的形变信息,实施坐标平面内的平滑处理、对比度优化处理、帧与帧之间的时间轴方向的平滑处理等各种图像处理,之后向彩色扫描转换器15送出。
彩色扫描转换器15,获得形变分布修正部14修正后的形变分布,按照设定好的形变的色彩映射(colour map),对每个形变分布的帧数据的像素附加色调代码,生成彩色形变图像。
彩色扫描转换器15所生成的彩色形变图像,通过切换加法部9被显示在图像显示器10上。此外,切换加法部9形成为具有以下功能,即输入黑白扫描转换器8输出的黑白断层图像、彩色扫描转换器15输出的彩色形变图像,切换两种图像来显示任意一方的功能;使两种图像中的一个半透明化,进行相加合成,在图像显示器10上重叠显示的功能;并列显示两种图像的功能。此外,由切换加法部9输出的图像数据,根据装置控制接口部19的控制,会被保存在场景存储器(cine memory)20中。场景存储器20存放的图像数据,根据装置控制接口部19的控制,会被显示在图像显示器10上。
与本实施方式的特征有关的形变分布修正函数制作部18,获得例如加压机构(图1的情况下就是探头2)与被测体1之间的接触面形状、测量对象的视野范围大小(边界条件)、衍射角Φ等条件,它们是从装置控制接口部19输入的加压测量条件。然后,形变分布修正函数制作部18,运算或者选择并设定以下实施例所示的形变分布修正函数。设定的形变分布修正函数被输出至形变分布修正部14。
对如上构成的本实施方式的基本动作进行说明。首先,用探头2压迫被测体1,对被测体1进行超声波束扫描,连续接收来自扫描面的反射回声信号。然后,根据调相加法电路6输出的反射回声信号,由信号处理部7和黑白扫描转换部8对断层图像进行重组,经由切换加法器9,显示在图像显示器10上。
另一方面,RF信号帧数据选择部11取入反射回声信号,与帧率同步,来反复取得帧数据,按照时间顺序保存在内置的帧存储器内。然后,以取得时刻不同的一对帧数据为单位,连续选择多对帧数据,输出至位移运算部12。位移运算部12,对选择的一对帧数据进行一维相关或二维相关的处理,测量扫描面上的多个测量点的位移,生成位移帧数据。作为上述位移矢量的检测方法,已知有例如特开平5-317313号公报等所述的块映射(block mapping)法或渐变法。块映射法是将图像分成块,每个块例如由N×N个像素组成,从前一帧起查找最近似于当前帧中的所着眼的块,由此求出测量点的位移。此外,还可以对一对RF信号帧数据进行在同一区域的自相关计算,算出位移。
形变信息运算部13取入位移帧数据,求出各测量点的形变变化,将形变分布(帧数据)输出至形变分布修正部14。众所周知,形变变化的运算,是对各测量点的位移做空间微分,计算各测量点的形变变化Δε。此外,也可以像非专利文献1等提出的那样,设定视野范围中的感兴趣区ROI和基准区ROI0,求出这些区域中的形变变化Δε、Δε0的平均值,根据其比值(Δε0的平均值/Δε的平均值)的大小,来鉴别组织是良性还是恶性。
形变分布修正部14对输入的形变分布进行平滑处理等处理,然后,使用从形变分布修正函数制作部18输入的形变分布修正函数,对形变分布进行修正,将基于修正后的形变分布的形变信息输出至彩色扫描转换器15。彩色扫描转换器15,根据形变分布生成彩色形变图像。彩色形变图像按照像素单位,使用例如256级的色调灰度,根据帧数据的形变,进行着色。另外,可以使用黑白扫描转换器来取代彩色扫描转换器15。在这种情况下,可以通过对形变大的区域提高亮度,对形变小的区域减小亮度等方法,鉴别是良性还是恶性。
以下,针对使用本实施方式的形变分布的修正,对基于不同加压方法和不同加压测量条件的具体实施例进行说明。另外,各实施例由作为本实施方式特征部分的形变信息运算部13、形变分布修正部14、形变分布修正函数制作部18、装置控制接口部19等实施。
因此,本发明是预先对每个加压测量条件,都进行图4(A)的形变分布的测量,由形变分布修正函数制作部18设定形变分布修正函数,形成接触面上的应力在任意深度上都不衰减的形变分布。然后,形变分布修正部14利用形变分布修正函数,对形变信息运算部13得到的形变分布进行修正,由此获得与距离加压机构的深度和方向无关的准确的形变信息。
实施例1
本实施例1是以图2所示的线形探头21作为加压机构,将该探头21的超声波收发面(接触面)推挤并压迫于被测体1,对该情况下的形变信息进行修正。线形探头21的接触面与被测体1相比,具有足够的硬度,在测量范围的压迫力下,是不会变形的。
此外,如图3所示,设接触面23的x轴方向的长度为2·x0,设y轴方向的长度为2·y0,设接触面23上的应力σ为σ=σ0(z=0)。现在,假定施加在接触面23上的压迫力的弹性波,向压缩方向大范围传播,衍射角为Φ。在弹性波的路径区域内的任意xy平面(z=恒定)上,设定一种模型,其中如下式(1)所示,z方向的应力σ(z)与x、y坐标无关,是恒定值。也就是说,施加在被测体1上的压迫力(外力)=应力×面积是恒定的,与深度无关。在式(1)中,ds是面积微元(微小面积要素)。
∫σ(z)·ds=恒定 (1)
此外,如果设衍射角为Φ,那么深度z下的x方向上的扩展范围Ux(z)由下式(2)表示。
Ux(z)=2(x0+z·tanφ) (2)
同样,如果设衍射角为Φ,那么深度z下的y方向上的扩展范围Uy(z)由下式(3)表示。
Uy(z)=2(y0+z·tanφ) (3)
根据这些公式,下式(4)成立。
σ0·(2x0)·(2y0)=σ(z)·2(x0+z·tanφ)·2(y0+z·tanφ) (4)
根据式(4),z轴的任意位置的应力变为下式(5)。
σ(z)=σ0·(x0·y0)/
{(x0+z·tanφ)·(y0+z·tanφ)} (5)
这里,衍射角Φ依赖于弹性波的频率(反复压迫操作的频率),在例如Φ=π/4的条件下,式(5)就变为式(6)。
σ(z)=σ0·{x0·y0}/{(x0+z)·(y0+z)} (6)
此外,如果是接近接触面的较浅部位,则认为z<<x0、y0,于是σ(z)≒σ0(恒定);如果是深部位,则认为z>>x0、y0,于是σ(z)≒σ0·{x0·y0}/{z·z}。
因此,深部的应力会大大衰减,以1/z2的关系变化。其结果,即便是硬度均匀的生物体组织,应力σ(z)也会随传播衰减,由此,就得到形变值也随之减小的形变分布。
在本实施例1中,形变分布修正部14加入上述的应力衰减,对形变分布进行修正,根据修正后的形变分布(以下为修正形变分布),构筑形变信息。以下详细说明本实施例1的形变分布的具体修正方法。
设:利用图2的探头21,在上述加压测量条件下,进行形变分布测量。此外,假定视野范围24的生物组织的弹性是均匀的。设测量得到的形变分布数据为ε(x,z)。将这时的视野范围置为—x0≦x≦x0,0≦z≦z0。在上述条件下,例如中央线x=0的深度方向上的形变分布ε(0,z),就会变为图4(A)所示的分布,随着应力的衰减,在深度方向上减小。而且,基于该形变分布ε(x,z)的形变信息就会变成图4(B)那样,越往深处的部位形变越小,因此,有可能误认为深部区存在硬的部位。
因此,本实施例加入上式(6),在形变分布修正函数制作部18上用下式(7)定义形变分布修正函数w(z)。形变分布修正函数w(z),是式(6)所示的应力衰减量的倒数。
w(z)={(x0+z)·(y0+z)}/{x0·y0} (7)
进而,形变分布修正部14用下式(8)求出修正形变分布ε’(x,z)。
ε’(x,z)=w(z)·ε(x,z) (8)
形变分布修正部14,将应力的衰减量的倒数与形变分布相乘,从而对形变分布进行修正。也就是说,形变分布修正部14,使用从形变分布修正函数制作部18输出的形变分布修正函数w(z),预测应力的衰减,来对形变分布进行修正。由此,修正过的形变分布就会像图5(A)那样,在深度方向上分布均匀,按照基于修正后的形变分布的形变信息而形成的弹性图像,也会像图5(B)那样,在整个区域没有形变大小的差异,可以防止错误识别。
此外,如果是对具有硬度不相同的区域的生物体组织进行测量,那么使用形变分布修正函数w(z),可以高精度地求出其硬度差异。
另外,虽然是假定衍射角Φ=π/4、用式(6)进行了近似运算,但本发明不限于此,也可以进行可变设定。此外,也可以预先将弹性波衍射角Φ设定为重复压迫操作频率的函数,测量重复压迫操作频率,并使用式(5),推定应力的衰减,来设定形变分布修正函数w(z)。
实施例2
本实施例2是使用图6所示的凸形经直肠探头31来作为加压机构,对安装在经直肠探头31前端的球状气囊33进行膨胀·收缩来压迫被测体,对该情况下的形变信息进行修正。气囊33,被安装成包围凸形的超声波收发面32,是通过连通到气囊33内的流路34,从注射器等注入·排出水来对气囊33进行膨胀·收缩的情况下的一个例子。
如上所述,应力的衰减取决于施加压迫力的接触面的形状,且取决于因弹性波的衍射而带来的大范围传播。也就是说,应力的衰减,在视野范围相对于接触面面积非常广阔的加压测量条件下表现得非常显著。像本实施例2这样的经直肠探头31等体内插入型探头就是如此。作为体内插入型探头,例如还有经阴道探头、经食道探头等。
用图6所示的球状气囊33施加压迫力、测量弹性的方法,在专利文献1中已被提出。本实施例2的情况下,气囊33的膜面与被测体体腔内的表皮接触,水注入气囊33时,膜面膨胀的方向如图7所示,向球面的法线方向压迫生物体组织。对于图7(A)所示的xy平面,按照与实施例1同样的条件,在使用相较于压迫对象而言充分的力、且可以保持球面地实施压迫的前提下,讨论应力的传递。
本实施例的情况也与图2的压迫操作同样,在初始状态下实施压迫,然后反复调整压迫力,压迫生物体组织。现在,设初始状态下的气囊33的曲率半径为r0,设接触面上的应力σ为σ=σ0(r=r0)。然后,如图8所示,xy平面内的测量点坐标用(r,θ)指定。此外,设接触面36上产生的弹性波是球面波,向球面的法线方向传播。这时,与实施例1同样,构筑一种模型,力=应力×面积是恒定的,不依赖于深度。也就是说,根据式(1),在本实施例2的情况下,下式(9)成立。
σ0·4π(r0)2=σ(r)·4π(r)2 (9)
因此,σ(r)就由下式(10)求出。
σ(r)=σ0·(r0/r)2 (10)
由式(10)可知,σ(r)大大衰减,以1/r2的关系变化。因此,在本实施例中,形变分布修正函数制作部18加入式(10),用下式(11)定义形变分布修正函数w(r),形变分布修正函数w(r),是式(10)所示的应力的衰减量的倒数。
w(r)=(r/r0)2 (11)
形变分布修正部14,利用上述的形变分布修正函数w(r),对形变分布进行修正,可以通过下式(12)求出修正形变分布ε’(r,θ)。
ε’(r、θ)=w(r)×ε(r,θ) (12)
根据本实施例,就会与实施例1同样,形变分布修正部14通过将应力的衰减量的倒数与形变分布相乘,从而对形变分布进行修正。可以在修正形变信息的整个区域排除因应力衰减而带来的形变大小的差异,按照基于修正过的形变分布的形变信息来防止误诊。此外,如果是对具有硬度不相同的区域的生物体组织进行测量,那么使用形变分布修正函数w(r),可以高精度地求出其硬度差异。
实施例3
实施例2中,表示了以具有球状膜面的气囊33为加压机构的情况。本实施例3是说明将具有图9(A)、(B)所示的圆筒状膜面的气囊41作为加压机构使用时的形变分布修正函数的例子。就本实施例的气囊41而言,圆筒状膜面与压迫对象接触,在保持圆筒膜面的状态下进行膨胀·收缩,就会在圆筒膜面的法线方向对压迫对象进行压迫。本实施例的情况下气囊41与压迫对象之间的接触面非常大,图9(B)的z轴方向的长度2·z0比视野范围的半径r的尺寸大得多,这时,与与实施例1的浅部的状况同样,应力在yz平面内的传播可以忽略衰减。也就是说,在使用这种足够大的圆筒状的气囊41进行压迫的情况下,只需考虑图8所示的xy平面内的应力衰减即可。根据式(1)的模型条件,设定下式(13)、(14)。
σ0·2πr0=σ(r)·2πr (13)
σ(r)=σ0·(r0/r) (14)
由此可知,本实施例情况下的应力大大衰减,以1/r的关系变化。因此,根据式(14),形变分布修正函数制作部18用下式(15)定义形变分布修正函数w(r)。形变分布修正函数w(r)是式(14)所示的应力的衰减量的倒数。
w(r)=(r/r0) (15)
然后,形变分布修正部14利用形变分布修正函数w(r),对测量的形变分布ε(r,θ)进行修正,由下式(16)求出修正形变分布ε’(r,θ)。
ε’(r,θ)=w(r)×ε(r、θ) (16)
由此,根据本实施例,就会与实施例1、2同样,形变分布修正部14通过将应力的衰减量的倒数与形变分布相乘,对形变分布进行修正。这样,根据基于修正过的形变分布的形变信息,就可以防止误诊。此外,如果是对具有硬度不相同的区域的生物体组织进行测量,那么使用形变分布修正函数w(r),可以高精度地求出其硬度差异。
此外,如果是实施例2的气囊,应力就会按照1/r2的关系衰减,如果是实施例3的气囊,应力就会按照1/r的关系衰减。在加压机构使用气囊的情况下,应力的衰减特性会因膨胀·收缩方式和接触面大小的不同而不同。
此外,本实施例3的形变修正法,也可以应用在将血管管壁因搏动而膨胀·收缩的现象视为加压机构,利用施加在血管壁或周围的生物体组织上的压迫力测量形变的情况。例如,可以应用在利用颈动脉的搏动的甲状腺诊断、利用下肢动脉的脉搏的深部静脉血栓诊断等中。
再有,也可以应用在以下情况,即,取代利用搏动的方法,将球囊扩张导管等、对插入血管内的气囊进行膨胀·收缩的方法作为加压机构,利用施加在血管壁或周围的生物体组织的压迫力测量形变。
实施例4
本实施例4是将凸形探头2本身作为加压机构使用,并对此时的形变分布进行修正的例子。实施例1~3中,加压测量条件都是:在视野范围内,压迫对象在深度方向上的压迫力均匀衰减。本实施例4是以图6所示的经直肠探头31为加压机构,不使用气囊,向压迫对象施加压迫力的情况下的形变修正的例子。
凸形探头2在超声波收发面32的长轴方向上具有曲率,例如像图10所示的那样,将长轴中心在法线方向上移动,来对压迫对象实施压迫。在这种情况下,接触面上的压迫力方向就与扇形的视野范围的深度方向不同,所以,施加在接触面上的y轴方向的压迫力的视野范围的深度方向的分力,就成为往超声波束方向有效的压迫力。因此,根据实施例4的加压测量条件,就会具有如下特性:随着超声波束的方向,接触面上的压迫力变得不均匀。其结果,在视野范围内,应力分布就会不均匀,由此可能带来不均匀的形变分布,导致误诊。
在图10中,若设探头2沿图的y轴方向移动并进行压迫,则被压迫的方向的范围至少为0<θ<π,没被压迫的方向的范围是π≦θ≦2π。在被压迫的方向的范围内,在接触面的任意坐标(r0,θ)上,大小一定的压力σ0都会如图所示地施加在y轴方向。与气囊的情况不同,接触面上的压力的法线方向的成分根据坐标(r0,θ)的不同而不同。因此,法线方向的成分σ0’就会如下式(17)所示,随sinθ的变化而变化。θ是法线与x轴之间的成角。
σ0’(θ)=σ0·sinθ (17)
由此,在对大角度的视野范围进行诊断的情况下,视野范围的中央部(θ=π/2附近)应力大,形变大;视野范围的两侧部(θ=0或π附近)应力变小。当然,视野范围内的形变的测量值也是随应力而变化的值,因此,例如有可能发生误诊,认为两侧部相对于中央部,存在硬组织。所以,在本实施例中,鉴于加压测量条件的不均匀性,形变分布修正函数制作部18按照以下方式设定形变分布修正函数w(θ),形变分布修正部14对形变分布进行修正。
首先,形变分布函数制作部18根据式(17),以下式(18)设定形变分布修正函数w(θ)。形变分布修正函数w(θ)是式(17)表示的应力的衰减量的倒数。
w(θ)=1/(sinθ) (18)
因此,形变分布修正部14用下式(19),对测量的形变分布ε(r,θ)进行修正,求出修正形变分布ε’(r,θ)。
ε’(r,θ)=w(θ)×ε(r、θ) (19)
形变分布修正部14,通过将应力的衰减量的倒数与形变分布相乘,对形变分布进行修正。进而,根据随深度而衰减的应力,使用例如实施例3的形变分布修正函数w(r),可以同样进行修正。也就是说,根据视野范围的大小,特别是根据深度范围的不同,深度方向所对应的应力分布的衰减效果也同时出现。在这种情况下,构筑作为r和θ的函数的形变分布修正函数w(r,θ)。例如,在实施例3的测量条件下,进一步在深度方向上确认式(14)的关系的情况下,作为形变分布修正函数w(r,θ)使用下式(20)。
w(r、θ)=w(r)×w(θ)
=(r/r0)×(1/sinθ) (20)
另外,在本实施例4中,虽然假定施加在两侧部的压力仅为垂直成分,但实际上,探头2的移动会产生被推挤的部位,在该部位会产生横向压迫力成分,所以优选加入这种情况,设定形变分布修正函数w(r,θ)。
实施例5
实施例1~4的形变分布修正函数,由于是仅与视野范围的坐标相应的函数,所以,优选例如按照加压机构的压迫条件,预先对一个帧的每个坐标位置,求出形变分布修正函数的计算值(构成值),并在形变分布修正函数制作部18内的存储器中保存。这样,参照与形变的测量值的坐标相对应的构成值,就可以进行实时修正。
此外,形变分布修正函数不限于实施例1~4所示的(1/r)、(1/r2)的函数,可以使用对数函数或指数函数等,适合对加压测量条件进行解析的函数。
另外,有时会因探头2的形状或压迫对象的形状等加压测量条件的不同,使应力分布传递十分复杂,很难模式化地构筑准确的形变分布修正函数。在这种情况下,例如形变分布修正函数制作部18,可以使用有限要素法等仿真方法,构筑形变分布修正函数。
根据利用了模体(phantom)的实测值,制作形变分布修正函数,将该形变分布修正函数预先保存在形变分布修正函数制作部18内的存储器中,形变分布修正部14就可以按照存储的形变修正函数,修正形变分布。
实施例6
这里,对经由图1的装置控制接口部19来控制形变分布修正函数制作部18的实施例进行说明。
如实施例1~5所示,所使用的形变分布修正函数,需要结合探头2或气囊等加压机构的差异、视野范围的形状、以及衍射角Φ等加压测量条件来进行切换。因此,要设法在装置控制接口部19上设置形变分布修正函数的切换机构以及ON/OFF的切换机构,使检测者能够结合加压测量条件进行切换。具体就是,形变分布修正函数制作部18,根据探头2和加压机构的种类设定形变分布修正函数,并保存在存储器中。
例如在对超声波诊断装置设置线形探头21,并通过装置控制接口部19将测量模式切换为线形扫描的情况下,形变分布修正函数制作部18使用实施例1的形变分布修正函数。此外,在对超声波诊断装置设置凸形经直肠探头31,并通过装置控制接口部19将测量模式切换为凸形扫描的情况下,形变分布修正函数制作部18使用实施例2的形变分布修正函数。这样,对每种探头2的形状都准备一个固定的形变分布修正函数,就可以做到在切换探头2时,自动地随之切换形变分布修正函数。
此外,如果利用了气囊,在B模式图像上,就会观察到没有反射到膜面上的回声信号的区域。通过在B模式图像上检出无回声的层状区域,就可以认为试用了气囊,可以自动切换成气囊用形变分布修正函数。
此外,也可以设置切换机构,对是不是进行形变分布修正处理进行切换。另外,可以设法读出保存在场景存储部中的形变信息,切换各实施例中表示的形变分布修正函数,制作并比较修正形变信息。
此外,通常的超声波诊断装置设有TGC(time·gain·control)和STC(sensitivity·time·control)等功能,它们是根据测量深度调整接收信号感度的功能。通过这些功能,就能借助微调开关,使测量深度的各个位置上的感度得到调整。因此,为了能根据深度来独立调整形变分布修正函数,可以像图11所示的那样,设置形变分布修正函数的微调开关(也包含侧方向(lateral))。而且,例如,在判定为深部的修正强度较弱的情况下,可以调整与该深部对应的开关,加强修正效果。
在这种情况下,形变分布修正函数的微调开关,也可以对当前选择的形变分布修正函数进行微小程度的调整,此外,它还可以限制形变分布修正函数,以免变更巨大。
此外,也可以更换B模式图像用TGC调整开关,替代为形变分布修正函数的调整开关。
另外,不仅可以进行w(r)的深度方向的修正,同样也可以进行w(θ)的角度方向(侧方向)的调整。
此外,虽然上述各实施例示出的是设定形变分布修正函数,来对形变分布进行修正的例子,但本发明不限于此,也可以事先设定位移的修正函数,对用于形变运算的位移分布进行修正,从而求出修正形变分布,这样可以起到与上述各实施例相同的效果。具体而言就是,预先对每个加压测量条件进行位移分布测量,使用位移分布形变函数制作部(未图示)来设定位移分布修正函数,形成使接触面上的应力在任意深度方向上都不衰减的位移分布。然后,位移分布修正部(未图示。),用位移分布修正函数来修正从位移运算部12得到的位移分布,就会得到准确的位移信息,不论距离加压机构的深度和方向如何。然后,形变信息运算部13,根据修正后的位移信息,求出形变分布。