具体实施方式
使用图1来说明应用本发明的超声波诊断装置100。
如图1所示,在超声波诊断装置100中具备:与被检体101抵接来使用的超声波探头102;经由超声波探头102每隔一定的时间间隔就对被检体101反复发送超声波的发送部105;接收从被检体101反射的回波信号的接收部106;对发送部105和接收部106进行控制的收发控制部107;和对由接收部106接收到的反射回波进行整相相加的整相加法运算部108。
超声波探头102配设多个振子而形成,具有经由振子101来对被检体101进行超声波收发的功能。超声波探头102由形成为矩形或扇形的多个振子构成,在与多个振子的排列方向正交的方向上使振子机械地振动,并能三维地进行超声波的收发。另外,超声波探头102也可以二维排列多个振子,电子地控制超声波的收发。
收发部105驱动超声波探头102的振子来生成用于产生超声波的波发送脉冲。发送部105具有将所发送的超声波的收敛点设定于某一深度的功能。另外,接收部106对由超声波探头102接收到的反射回波信号以规定的增益进行放大,并生成RF信号即接收信号。超声波收发控制部107用于控制发送部105、接收部106。
整相加法运算部108对由接收部106放大后的RF信号的相位进行控制,对于1个或多个收敛点,形成超声波束从而生成RF信号帧数据(相当于RAW数据)。
进而,在超声波诊断装置100中,具备:存储由整相加法运算部108生成的RF信号帧数据的数据存储部109;基于由数据存储部109存储的RF信号帧数据来构成二维断层图像的二维断层图像构成部113;针对由二维断层图像构成部113构成的二维断层图像,基于二维断层图像的取得位置来进行三维坐标变换,并生成断层容积数据的断层容积数据生成部115;从断层容积数据来构成三维断层图像的三维断层图像构成部117;基于由数据存储部109存储的多个RF信号帧数据来构成二维弹性图像的二维弹性图像构成部114;针对由二维弹性图像构成部114构成的二维弹性图像,基于二维弹性图像的取得位置来进行三维坐标变换,并生成弹性容积数据的弹性容积数据生成部116;从坐标变换后的弹性容积数据来构成三维弹性图像的三维弹性图像构成部118;合成二维断层图像与二维弹性图像、或合成三维断层图像与三维弹性图像的合成处理部110;和对由合成处理部110合成后的合成图像、二维断层图像或二维断层图像等进行显示的显示部111。
另外,在超声波诊断装置100中具备:控制上述各构成要素的控制部103、和对控制部103进行各种输入的操作部104。操作部104具备键盘、轨迹球等。
二维断层图像构成部113基于控制部103中的设定条件,输入从数据存储部109输出的RF信号帧数据,并进行增益校正、日志压缩、检波、轮廓增强、滤波处理等信号处理,来构成二维断层图像。
超声波探头2能与超声波的收发同时地计测收发方向(θ、),断层容积数据生成部115基于相当于二维断层图像的取得位置的收发方向(θ、),对多个二维断层图像进行三维变换,来生成断层容积数据。
三维断层图像构成部117基于断层容积数据的亮度和不透明度来进行容积再现。具体地,三维断层图像构成部117如图2(a)所示那样,由断层不透明度表创建部200和断层再现运算部202构成。
断层不透明度表创建部200创建取体素值(voxel value)即断层图像的亮度为横轴、取不透明度为纵轴的断层不透明度表。断层再现运算部202基于所创建的断层不透明度表,使用根据断层容积数据来构成三维断层图像的下式来进行容积再现。
[数1]
Cout(i)=Cout(i-1)+(a-Aout(i-1))·A(i)·C(i)·S(i) (1)
Aout(i)=Aout(i-1)+(1-Aout(i-1))·A(i) (2)
A(i)=Opacity[C(i)] (3)
C(i)是在从所创建的二维投影面上的某点来观察三维断层图像的情况下存在于视线上第i个的体素的亮度值。Cout(i)是输出的像素值。例如,在视线上N个体素的亮度值排列时,对i=0~N-1进行了累积的亮度值Cout(N-1)成为最终输出的像素值。Cout(i)表示直到第(i-1)个为止的累积值。
另外,A(i)是存在于视线上第i个的亮度值的不透明度,如上述(3)所示那样,是取0~1.0的断层不透明度表(断层opacity table)。断层不透明度表由断层不透明度表创建部200创建,通过根据亮度值来参照不透明度,来决定对输出的二维投影面(三维断层图像)上的贡献率。
S(i)是用于通过根据亮度C(i)及其周边的像素值求得的斜率而推算的阴影添加的权重分量,例如,指示如下等的增强效果:在光源与以体素i为中心的面的法线一致的情况下,反射最强,因此赋予1.0,在光源与法线正交的情况下,赋予0.0。
Cout(i)、Aout(i)均以0为初始值。如上述(2)所示,Aout(i)每当通过体素时就进行累积,并收敛于1.0。因而,如上述(1)所示那样,在直到第(i-1)个为止的不透明度的累积值Aout(i-1)成为约1.0的情况下,第i个以后的亮度值C(i)将不会反映在输出图像上。
断层不透明度表创建部200例如通过将断层不透明度表中的低亮度的不透明度设为零,从而能不对亮度低的断层容积数据进行图像化,而是对亮度高的断层容积数据进行图像化。断层再现运算部202基于所创建的断层不透明度表来进行容积再现,并将所构成的三维断层图像输出到合成处理部110。
在本实施例中,能变更断层不透明度表。通过基于操作部104的操作的控制部103的控制,使断层不透明度表创建部200变更断层不透明度表,能与变更后的断层不透明度表相应地来对断层容积数据进行图像化。
二维弹性图像构成部114根据存储于数据存储部109中的多个RF信号帧数据来对位移进行计测。然后,二维弹性图像构成部114基于计测出的位移来运算弹性值,从而构成二维弹性图像。在弹性值中含有失真、弹性率、位移、粘性、失真比等弹性信息。
弹性容积数据生成部116基于相当于二维弹性图像的取得位置的收发方向(θ、)来对多个二维弹性图像进行三维变换,生成弹性容积数据。按照弹性值来对弹性容积数据赋予色彩值(蓝色、淡蓝色、绿色、黄色、红色等)。例如,对硬的部位赋予蓝色,对软的部位赋予红色,从而进行颜色区分。
三维弹性图像构成部118基于弹性值来设定弹性容积数据的容积再现中的不透明度,从而构成三维弹性图像。通过设定了不透明度与弹性值的关系的弹性不透明度表来设定不透明度。具体地,三维弹性图像构成部118如图2(b)所示那样,由弹性不透明度表创建部204和弹性再现运算部206构成。
弹性不透明度表创建部204创建设定了不透明度与弹性值的关系的弹性不透明度表。创建根据弹性值来设定容积再现中的不透明度的弹性不透明度表。在图3、图4中示出由弹性不透明度表创建部204创建的弹性不透明度表。通过弹性不透明度表来决定弹性值与不透明度的关系,设定显示部111中的对二维投影面上的弹性容积数据的贡献率。
弹性不透明度表由用于根据弹性值来设定不透明度的的不透明度线300、和进行用于设定不透明度线300的引导的引导线302构成。另外,在图3、图4中也示出了与不透明度表的弹性值对应的色带304。
如图3所示,弹性不透明度表根据不透明度线300来设定不透明度,以使根据弹性值而被赋予了蓝色的关注部位(硬的部位)得以增强来进行显示。基于不透明度线300来设定不透明度,以使得根据弹性值而决定的关注部位(硬的部位)的不透明度比被赋予了绿色或红色的关注部位以外的部位(标准部位、软的部位)的不透明度高。成为关注部位以外的部位的标准部位是指在弹性容积数据中具有平均的硬度的部位,是主要被赋予绿色的部位。另外,关注部位以外的部位是主要包含标准部位的部位。
在弹性不透明度表中,若设赋予了蓝色的关注部位(硬的部位)的在不透明度线300中的不透明度为T1,则T1被设定为接近1.0的值(例如0.7≤T1≤1.0)。另外,若设赋予了绿色或红色的关注部位以外的部位(标准部位、软的部位)的在不透明度线300中的不透明度为T2,则将T2设定为接近于0的值(例如,0<T2≤0.5)。另外,对不透明度线300的不透明度进行设定,以使得关注部位以外的部位的不透明度T2不成为零。这是因为要与赋予了蓝色的关注部位(硬的部位)进行对比显示。另外,若将赋予了淡蓝色的介于关注部位和关注部位以外的部位之间的部位的在不透明度线300上的不透明度设为T3,则将T3设定为T2和T1之间的值(例如,T2≤T3≤T1)。
在本实施例中,能进一步变更弹性不透明度表的不透明度。通过基于操作部104的操作的控制部103的控制,使弹性不透明度表创建部204变更弹性不透明度表的不透明度。
弹性不透明度表创建部204变更与赋予了绿色等的关注部位以外的部位(包含标准部位)相应的不透明度线300的不透明度。
具体地,若操作者将图3的左图所示的设定为关注部位以外的部位的不透明度线300朝向上方向移动,则弹性不透明度表创建部204设定弹性不透明度表的不透明度,从而在如图3右图所示,关注部位以外的部位的不透明度提高。反之,若操作者将设定为关注部位以外的部位的不透明度线300朝向下方向移动,则弹性不透明度表创建部204设定弹性不透明度表的不透明度,从而关注部位以外的部位的不透明度降低。进行设定不透明度线300的引导的引导线302表示不透明度线300的可移动范围(例如0<T2≤0.5)。如此,操作者能设定被设定为关注部位以外的部位(包含标准部位的部位)的在不透明度线300上的不透明度。
另外,在本实施例中,以直线示出了不透明度线300,但也可以是曲线。在不透明度线300是曲线的情况下,弹性不透明度表创建部204设定关注部位以外的部位(包含标准部位的部位)中的曲线的不透明度。
弹性再现运算部206使用基于由弹性不透明度表创建部204创建的不透明度线300的弹性不透明度表以及下述视线上来进行容积再现。
[数2]
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)=Opacity[E(i)] (6)
E(i)是在从所创建的的二维投影面上的某点观察三维弹性图像的情况下的存在于视线上的第i个的弹性值。Eout(i)是所输出的像素值。例如,在视线上排列N个体素的弹性值时,i=0~N-1为止对弹性值进行累积而得到的累积值Eout(N-1)成为最终输出的像素值。Eout(i-1)表示到第i-1为止的累积值。另外,A(i)是存在于视线上的第i个的弹性值的不透明度,如式(6)所示,是由弹性不透明度表创建部204创建的弹性不透明度表。
S(i)是用于通过根据由弹性值E(i)及其周边的弹性值求得的斜率而推算的阴影添加的权重分量,例如,指示如下等的增强效果:在光源与以体素i为中心的面的法线一致的情况下,反射最强,因此赋予1.0,在光源与法线正交的情况下,赋予0.0。
Cout(i)、Aout(i)均以0为初始值。如上述(5)所示,Aout(i)每当通过体素时就进行累积,并收敛于1.0。因而,如上述(4)所示那样,在直到第(i-1)个为止的体素的不透明度的累积值Aout(i-1)成为约1.0的情况下,第i个以后的体素值E(i)将不会反映在输出图像上。
合成处理部110分别读取如上所述而构成的三维断层图像、三维弹性图像,并将它们与在二维投影面上由操作部104设定的透过率相乘,然后累加,来创建三维合成图像。显示部111对由合成处理部110创建的三维合成图像进行显示。
另外,弹性不透明度表创建部204如图4(a)所示,能够根据弹性不透明度表的不透明度线300来设定不透明度,以使根据弹性值而赋予了红色的关注部位(软的部位)得以增强来进行显示。弹性不透明度表按照根据弹性值而决定的关注部位的不透明度比关注部位以外的部位的不透明度高的方式来进行设定。具体地,根据不透明度线300来设定透明度,以使得根据弹性值而决定的关注部位(软的部位)的不透明度比赋予了绿色或蓝色的关注部位以外的部位(标准部位、硬的部位)的不透明度高。
弹性不透明度表创建部204能设定不透明度线400的不透明度,不透明度线400与通过基于操作部104的操作的控制部103的控制而赋予了绿色或红色等的关注部位以外的部位(标准部位、硬的部位)相应。若将如图4(a)的左图所示的设定为关注部位以外的部位的不透明度线400朝向上方向移动,则如图4(a)的右图所示那样,关注部位以外的部位的不透明度被较高地设定。反之,若将设定为关注部位以外的部位的不透明度线400朝向下方向移动,则关注部位以外的部位的不透明度被较低地设定。进行用于设定不透明度线400的引导的引导线402表示不透明度线400的可移动范围(例如0<T2≤0.5)。
另外,弹性不透明度表创建部204如图4(b)所示,根据弹性不透明度表的不透明度线404来设定不透明度,以使得根据弹性值而赋予了蓝色、红色的关注部位(硬的部位、软的部位)得以增强来进行显示。根据不透明度线404来设定不透明度,以使得根据弹性值而决定的关注部位(硬的部位、软的部位)的不透明度比赋予了绿色等的关注部位以外的部位(标准部位)的不透明度高。
弹性不透明度表创建部204能设定不透明度线404的不透明度,不透明度线404与通过基于操作部104的操作的控制部103的控制而赋予了绿色等的关注部位以外的部位(标准部位)相应。若将如图4(b)的左图所示的设定为关注部位以外的部位的不透明度线404朝向上方向移动,则如图4(b)的右图所示那样,关注部位以外的部位的不透明度被较高地设定。反之,若将设定为关注部位以外的部位的不透明度线404朝向下方向移动,则关注部位以外的部位的不透明度被较低地设定。进行用于设定不透明度线404的引导的引导线406表示不透明度线404的可移动范围(例如0<T2≤0.5)。
图5表示显示部111的一个显示形态。在同一画面上对图像500、不透明度表选择部502、不透明度设定部504、以及弹性不透明度表506进行显示。
图像500例如是合成了三维断层图像和三维弹性图像而得到的三维合成图像或三维弹性图像。
不透明度表选择部502从预先存储的多个弹性不透明度表中选择弹性不透明度表。具备存储多个上面说明的图3、图4(a)、(b)的弹性不透明度表的存储部(未图示)。在显示部111中显示多个弹性不透明度表,能通过基于操作部104的操作的控制部103的控制来选择弹性不透明度表。将由操作部104选择的弹性不透明度表506予以显示。在此,选择图4(b)的凹型的弹性不透明度表,如506那样显示。
不透明度设定部504能通过相当于操作部104的数值设定部或滑动条来设定与关注部位以外的部位(标准部位)相应的不透明度线404的不透明度。在此,将不透明度线404的不透明度设定为“0.5”。另外,数值或滑动条也能表示不透明度线400的不透明度的可设定范围。
由于在同一个画面上显示了画面500、不透明度表选择部502、不透明度设定部504、以及不透明度表,因此,操作者能一边确认图像500,一边选择不透明度表的图案,或者设定弹性不透明度表的不透明度,以使得成为最佳的图像。
以上,根据实施例1,在视线方向上进行容积再现从而构成三维弹性图像时,即使在三维弹性图像的视线上,在比想要观察的关注部位更靠近前存在具有厚度的关注部位以外的部位(包含标准部位的部位),也能通过按照弹性值来设定弹性不透明度表的不透明度,从而不使关注部位被关注部位以外的部位(包含标准部位的部位)隐藏地来显示关注部位。
例如,在关注部位位于视线方向的深处的情况下,通过按照使关注部位以外的部位(包含标准部位的部位)的不透明度变得更低的方式来进行变更,能透过近前的关注部位以外的部位(包含标准部位的部位)来观察关注部位。
即,能适当地显示三维弹性图像中的希望观察的关注部位。
在此,使用图1~图6来说明实施例2。与实施例1不同的点在于由弹性不透明度表创建部204来设定不包含标准部位的部位的不透明度。关于弹性不透明度表创建部204以外的构成,由于与实施例1相同,因此省略弹性不透明度表创建部204以外的构成的说明。
弹性不透明度表创建部204如图6的左图所示,根据弹性不透明度表的不透明度线600来设定不透明度,以使得根据弹性值而赋予了蓝色、红色的关注部位(硬的部位、软的部位)得以增强来进行显示。
由此,例如,在进一步增强赋予了蓝色的关注部位(硬的部位)的情况下,弹性不透明度表创建部204设定与不包含标准部位的赋予了红色的部位相应的不透明度线600的不透明度。具体地,若将设定为图6的左图所示的不包含标准部位的部位的不透明度的不透明度线600朝向下方向移动,则如图6的右图所示那样,不包含标准部位的部位的不透明度被较低地设定。
另外,在本实施例中,在进一步增强赋予了红色的关注部位(软的部位)的情况下,弹性不透明度表创建部204也能设定与不包含标准部位的赋予了蓝色的部位相应的不透明度线600的不透明度。
以上,根据本实施例,在视线方向上进行容积再现来构成三维弹性图像时,能进一步增强特别希望着眼的关注部位来进行显示。
使用图1~图8来说明实施例3。与实施例1、2不同的点在于,基于进行容积再现的弹性容积数据的体素数来按每个视线(方向)而设定不透明度表的不透明度。即,弹性不透明度表基于进行容积再现的弹性容积数据的体素数来设定不透明度。
如图7所示,三维弹性图像构成部118由以下构成:对弹性容积数据设定视线(方向)的视线设定部700;对在弹性容积数据中设定的视线(方向)上的体素数进行计数的体素计数部702;基于计数出的体素数来创建弹性不透明度表的弹性不透明度表创建部704;和使用所创建的弹性不透明度表来进行弹性容积数据的容积再现的弹性再现运算部706。
使用图8来具体说明本实施例中的三维弹性图像构成部118。在图8(a)的形态中,视线设定部700通过基于操作部104的操作的控制部103的控制,来对弹性容积数据800在对角线上设定有视线802。显示部111能对弹性容积数据800和视线802进行显示。体素计数部702对设定在弹性容积数据800中的视线802上的体素数进行计数。计数出的体素数用线806的长度来表征。另外,体素计数部702也可以进行用于搜索位于视线的最浅处或者最深处的位置的构造物的处理,来推算直到作为目的地的构造物为止的体素数。
弹性不透明度表创建部704首先如实施例1所示那样,根据不透明度线808来设定弹性不透明度表的不透明度,以使得根据弹性值而决定的赋予了蓝色、红色(硬的部位、软的部位)的不透明度比赋予了绿色等的关注部位以外的部位(标准部位)的不透明度高。
在图8(a)的形态中,由于对弹性容积数据800在对角线上设定了视线802,因此进行容积再现的体素数最多。因此,弹性不透明度创建部704根据进行容积再现的体素数多的这一情况,按照体素数来设定弹性不透明度表的不透明度,以使得关注部位(硬的部位、软的部位)不会被关注部位以外的部位(标准部位)隐藏。
具体地,若将赋予了绿色的关注部位以外的部位(标准部位)的不透明度线808的不透明度设为T2,则弹性不透明度创建部704按照关注部位以外的部位的不透明度T2降低的方式(例如0<T2≤0.1)来设定不透明度线808。另外,引导线810表示不透明度线808的可移动范围(例如0<T2≤0.5)。
弹性再现运算部706使用基于由弹性不透明度表创建部704创建的不透明度线808的弹性不透明度表以及上述数学式(数2),来进行容积再现。
另一方面,在图8(b)的形态中,视线设定部700通过基于操作部104的操作的控制部103的控制,按照不与弹性容积数据800的一边正交的方式来设定视线802。体素计数部702对设定于弹性容积数据800中的视线802上的体素数进行计数。计数出的体素数用线814的长度来表征。
弹性不透明度表创建部704首先与实施例1所述那样,根据不透明度线808来设定弹性不透明度表的不透明度,以使得根据弹性值而决定的赋予了蓝色、红色的关注部位(硬的部位、软的部位)的不透明度比赋予了绿色的关注部位以外的部位(标准部位)的不透明度高。
在图8(b)的形态中,按照不与弹性容积数据800的一边正交的方式来设定视线802,进行容积再现的体素数比图8(a)的形态要少。弹性不透明度创建部704按照体素数来设定弹性不透明度表的不透明度,以使得关注部位(硬的部位、软的部位)不会被关注部位以外的部位(标准部位)所隐藏。
具体地,若将赋予了绿色的关注部位以外的部位(标准部位)的在不太明显808上的不透明度设为T2,则弹性不透明度创建部704按照关注部位以外的部位的不透明度T2变高的方式(例如0.1<T2≤0.3)的方式来设定不透明度线808。
弹性再现运算部706使用基于由弹性不透明度表创建部704创建的不透明度线808的弹性不透明度表以及上述的数学式(数2)来进行容积再现。
在实施例3中,在进行容积再现的体素数多的情况下,弹性不透明度表创建部704按照使关注部位以外的部位的不透明度降低的方式来设定不透明度线808。另一方面,在进行容积再现的体素数少的情况下,弹性不透明度表创建部704按照使关注部位以外的部位的不透明度提高的方式来设定不透明度线808。因而,即使改变视线(方向)来构成三维弹性图像,也由于弹性不透明度表创建部704按每个视线(方向)对体素数进行计数,按照体素数来按每个视线(方向)适当地设定不透明度,因此,能显示不管从哪个视线(方向)都均匀的三维弹性图像。
使用图1~图10来说明实施例4。与实施例1~3不同点在于,分析视线(方向)上的弹性容积数据的累积不透明度,基于分析出的累积不透明度来设定弹性不透明度表的不透明度。
如图9所示,三维弹性图像构成部118由以下构成:对弹性容积数据设定视线(方向)的视线设定部900;对相对于所设定的视线(方向)上的视线方向距离(体素数)的累积不透明度进行分析的累计不透明度分析部902;基于分析出的累积不透明度来创建弹性不透明度表的弹性不透明度表创建部904;使用所创建的弹性不透明度表来进行弹性容积数据的容积再现的弹性再现运算部906。
使用图10来具体说明本实施例中的三维弹性图像构成部118。图10是表示相对于由视线设定部900设定的视线(方向)上的视线距离(体素数)的累积不透明度的图。在本实施例中,由于累积不透明度曲线1004整体上平衡良好地设定不透明度,因此设为最佳累积不透明度曲线。
累积不透明度分析部902对相对于所设定的视线(方向)上的视线距离(体素数)的累积不透明度进行分析。
在此,设由累积不透明度分析部902对视线(方向)上的累积不透明度进行分析的结果是得到了累积不透明度曲线1002。累积不透明度曲线1002相对于视线距离(体素数)的累积不透明度高,在视线距离(体素数)达到P1的时间点,累积不透明度达到1.0。
于是,弹性不透明度表创建部904按照视线(方向)上的累积不透明度曲线1002成为最佳累积不透明度曲线1004的方式来设定弹性不透明度表中的不透明度。具体地,弹性不透明度表创建部904例如如图4(b)所示那样,将设定为赋予了绿色等的关注部位以外的部位的不透明度线404朝向下方向移动,将关注部位以外的部位的不透明度设定得更低。
另外,设由累积不透明度分析部902对视线(方向)上的累积不透明度进行分析的结果是得到了累积不透明度曲线1006。累积不透明度曲线1006的相对于视线距离(体素数)的累积不透明度低,在视线距离(体素数)达到P2的时间点,累积不透明度达到1.0。另外,累积不透明度曲线1006在视线距离(体素数)P1的位置上,累积不透明度成为0.5。这种情况下,所输出的弹性值成为弹性容积数据的弹性值的二分之一,存在诊断为比实际更硬的可能性。
于是,弹性不透明度表创建部904按照视线(方向)上的累积不透明度曲线1006成为最佳累积不透明度曲线1004的方式来设定弹性不透明度表中的不透明度。具体地,弹性不透明度表创建部904例如如图4(b)所示那样,将设定为赋予了绿色等的关注部位以外的部位的不透明度线404朝向上方向移动,将关注部位以外的部位的不透明度设定得更高。
弹性再现运算部906使用由弹性不透明度表创建部904创建的弹性不透明度表以及上述数学式(数2)来进行容积再现。
在实施例4中,对视线(方向)上的累积不透明度进行分析,并基于分析出的累积不透明度,来设定弹性不透明度表的不透明度,由此,能显示不管从哪个视线(方向)都均匀的三维弹性图像。
符号的说明
100 超声波诊断装置
102 超声波探头
130 控制部
104 操作部
105 发送部
106 接收部
107 收发控制部
108 整相加法运算部
109 数据存储部
110 合成处理部
111 显示部
113 二维断层图像构成部
115 断层图像三维坐标变换部
117 三维断层图像构成部
114 二维弹性图像构成部
116 弹性图像三维坐标变换部
118 三维弹性图像构成部