具体实施方式
以下,说明应用本发明而成的超声波诊断装置、超声波图像处理装置、超声波图像处理程序、以及超声波图像生成方法的实施方式。
(第一实施方式)
参照附图来说明应用本发明而成的超声波诊断装置的第1实施方式。图1是表示第1实施方式的超声波诊断装置的整体构成的框图。
如图1所示,本实施方式的超声波诊断装置100具备如下部件而构成:与被检体1之间收发超声波的超声波探头2;对超声波探头2提供驱动信号的发送部4;接收由超声波探头2接收到的反射回声信号的接收部6;进行发送部4以及接收部6的收发控制的超声波收发控制部8;对由接收部6接收到的反射回声信号进行相位累加来生成RF信号的相位累加部10。
此外,还具备:对从相位累加部10输出的RF信号(超声波信号),进行对数压缩、滤波处理、图像处理等各种处理来生成断层像数据的断层图像处理部12;对被检体1的多个断层面的断层像数据进行坐标变换,作成各轴正交的三维断层像体数据的断层像体数据生成部14;进行三维断层像体数据的各体素的亮度值的偏置运算的作为本实施方式的特征部分的偏置运算部16;和针对进行了偏置运算的三维断层像体数据利用体绘制技术来生成三维断层图像的断层像体绘制部18。
另一方面,还具备:基于来自断层图像处理部12的输出数据来生成三维断层像体数据的至少1个断层面上的二维断层图像的断层像二维坐标变换部20;对由断层像体绘制部18生成的三维断层图像和由断层像二维坐标变换部20生成的二维断层图像进行合成的图像合成部22;和作为显示从图像合成部22输出的图像的显示部的显示部24。
此外,还具备:作为接受来自检查者的指令等输入的输入接口的控制面板26;和基于来自控制面板26的输入指令来控制构成超声波诊断装置100的上述各构成要素的控制部28。例如,检查者通过控制面板26指定了三维断层像体数据的任意的剖面后,所指定的剖面位置信息被送往断层像二维坐标变换部20,生成其剖面位置上的二维断层图像。
在此,具体说明上述各构成要素。超声波探头2在超声波探头的长轴方向上排列了1~m信道(channel)的量的振子元件。在此,在短轴方向上也被切断为k个从而排列了1~k信道的量的振子元件的情况下,通过改变提供给短轴方向的各振子元件(1~k信道)的延迟时间,在短轴方向上也进行发送波或接收波的聚焦(focus)。此外,通过改变提供给短轴方向的各振子元件的超声波发送信号的振幅,来进行发送波加权,通过改变来自短轴方向的各振子元件的超声波接收信号的放大度或衰减度,来进行接收波加权。并且,通过开启、关闭短轴方向的各个振子元件,能够进行口径控制。
该超声波探头2一边根据来自控制部28的控制信号以电动机驱动或者手动方式在短轴方向上移动,一边进行扫描,能够进行被检体1的多个断层面的三维数据收集。此外,在短轴方向上也被切断为k个从而排列了1~k信道的量的振子元件的情况下,通过沿着探头头部的曲率的短轴方向、或通过电子聚焦而生成的短轴方向的超声波束能够收集三维超声波数据。
另外,该超声波探头2,根据重叠在从发送部4提供的驱动信号上而施加的偏置电压的大小,超声波收发灵敏度即机电耦合系数发生变化,例如可以应用cMUT(Capacitive Micromachined Ultrasonic Transducer:IEEETrans.Ultrason.Ferroelect.Freq.Contr.Vol45 pp.678-690May 1998等)。cMUT是通过半导体微加工工艺(例如,LPCVD:Low Pressure ChemicalVapor Deposition)而制造的超微电容式超声波振子。
发送部4以及接收部6对超声波探头2提供发送信号并处理接收到的反射回声信号,在其内部,具有对超声波探头2进行控制使其发出超声波束的发送波电路、和接收该所发出的超声波束的来自被检体1内的反射回声信号并收集生物体信息的接收波电路,通过超声波收发控制部8来控制。
相位累加部10对从接收部6输出的反射回声信号的相位进行控制,针对一个或多个收敛点形成超声波接收波束。另外,在相位累加部10中生成的RF信号也可以是复合解调后的I、Q信号。
断层图像处理部12对在相位累加部10中进行了相位累加的反射回声信号进行处理,具有如下部分而构成:基于依次输入的反射回声信号来进行对数压缩、滤波处理、图像处理的信号处理电路;和由存储超声波图像的磁盘装置以及RAM构成的存储装置。
断层像体数据生成部14基于由断层图像处理部12处理后的多个断层面的断层像数据来生成三维断层像体数据。偏置运算部16对三维断层像体数据的各体素进行亮度值的偏置运算。其详细内容见后述。
断层像体绘制部18对进行了偏置运算的三维断层像体数据进行体绘制等二维投影处理,生成三维断层图像并传输到图像合成部22。
断层像二维坐标变换部20基于来自断层图像处理部12的输出数据,通过重新取样以及插值处理来生成三维断层像体数据的至少1个断层面上的二维断层图像并传送到图像合成部22。图像合成部22对由断层像体绘制部18生成的三维断层图像和由断层像二维坐标变换部20生成的二维断层图像进行合成并传送到显示部24。显示部24输入由图像合成部22生成的图像,并显示为超声波图像,例如由CRT监视器、液晶监视器形成。
以下,更具体地说明本实施方式的超声波诊断装置100的断层像体绘制部18。断层像体绘制部18基于三维断层像体数据,在二维投影面上生成至少从一个视线方向来看的投影图像(三维断层图像),通过对三维断层像体数据中的视线方向的亮度值,乘以由控制部28传送来的每个亮度的透明度值并相加,来构成三维断层图像。在此,在本实施方式中使用的公知的体绘制方式的式子定义如下。
(式1)Cout=Cout-1+(1-Aout-1)·Ai·Ci
(式2)Aout=Aout-1+(1-Aout-1)·Ai
式1中,Ci是从所作成的二维投影面上的某点来观察三维断层像体数据的情况下的、在视线上的第i个体素亮度值,在视线上排列了N个体素的数据时,累计到i=0~N-1为止的值Cout成为最终输出像素值。Cout-1表示第i-1个为止的累计值。
此外,式1、2中的Ai是视线上的第i个体素值的不透明度,取0.0~1.0的值。Cout、Aout都将0作为初始值,如式2所示那样,Aout每当通过体素时被累计(累积相加)并收敛于1.0。因此,如式1所示,在第i-1个为止的体素的不透明度的累计值Aout-1成为约1.0的情况下,第i个体素值Ci不被反映于输出图像中。另外,假设不透明度为Ai时的透明度用1-Ai来表示,不透明度和透明度是互补关系。因此,在本说明书中,将透明度和不透明度的概念适当记载为透明度/不透明度。此外,例如在列举不透明度来进行了说明的情况下,同时也包含透明度的概念进行了说明。
图2是表示体素值和不透明度的关系的图。如图2所示,体素值和不透明度的关系一般表示为在横轴取亮度、在纵轴取不透明度的不透明度表(不透明度(opacity)表),根据体素的亮度值来参照不透明度。
如上所述,在本实施方式的体绘制处理中,能够捕捉不透明度较高的体素作为表面,从而立体地显示三维断层像数据。此外,作为将非表面的内部构造物透过性地可视化的绘制手法,一般使用只显示关心区域中的高亮度构造物的最大值投影法(Maximum intensity projection)、只描绘低亮度构造物的最小值投影法(Minimum intensity projection)、显示视线方向的体素值的累计图像的方法(Ray summation)等。
接下来,对本实施方式的超声波诊断装置的特征部分进行说明。本实施方式的超声波诊断装置具备如下部件而构成:与被检体1之间收发超声波的超声波探头2;基于由超声波探头2测量到的被检体1的多个断层面的反射回声信号,来生成三维断层像体数据的断层像体数据生成部14;基于三维断层像体数据,在二维投影面上生成至少从一个视线方向来看的三维断层图像的断层像体绘制部18;和显示三维断层图像的显示部24。具备根据三维断层像体数据的各体素的亮度值来对该各体素的亮度值进行增减的偏置运算部16,偏置运算部16的各体素的亮度值的增减量能够通过输入接口(控制面板)26来调整,断层像体绘制部18基于通过偏置运算部16对亮度值进行了偏置后的三维断层像体数据来生成三维断层图像。图3是表示在偏置运算部16中使用的偏置表的设定画面的一例的图。图3(a)(b)分别是对偏置表进行设定的一例。
如图3所示,偏置表设定有将三维断层像体数据的各体素可能取得的亮度值的范围(0~255)划分为多个区域(N:在本实施方式中为5)的边界值(P0~P5)、和所划分出的多个区域的每一个的亮度值的增减量(O1~O5)。
例如,利用图3(a)所示的图表31,通过控制面板26将P0~P5所示的边界值的设定指针向左右移动,由此能够设定区域的边界。同样,01~05是表示各区域的亮度值的增减量(偏置值)的指针,通过向上下移动能够设定各区域的偏置值。图表31的上限、下限值分别设为100,-100,但这是示出一例,只要在输入以及输出数据的可能取得的范围内则为多少都可以。另外,移动P0或P5,对于成为偏置的对象外的输入亮度值不进行偏置处理。此外,P0以及P5也可以固定于各体素可能取得的亮度值的范围的两端。此外,边界值(P0~P5)、亮度值的增减量(O1~O5)的设定可以利用切换器(toggle)、编码器或液晶面板上的调整按钮等来进行。
此外,如图3(b)所示,边界值以及亮度值的增减量(偏置值)也可以经由控制面板26,通过下拉菜单或直接数值输入来分别设定显示窗口32内的01~05、P0~P5。此外,如图3(b)所示,也可以将通过边界值(P0~P5)划分出的有效的区域数作为亮度调整区域数33(N=5)来进行显示。
图4是在偏置运算部16中进行的亮度值的变换处理的数据流程图。如图4所示,在偏置运算部16中开始偏置调整处理后,计数器i被初始化为0(步骤41)。接下来,读入三维断层像体数据的各体素的亮度值m(i)(步骤42)。判定i是否大于全部数据数(体素数)(步骤43),在大于的情况(在步骤43中为“是”的情况)下结束偏置调整处理。另一方面,在不大于、即没有对三维断层像体数据的所有的体素进行完偏置处理的情况下,进入到步骤44,区域选择用计数器n被初始化为1(步骤44)。
判断区域选择用计数器n是否小于亮度调整区域数N(步骤45),在小于的情况下(在步骤45中为“是”),判定被读入的三维断层像体数据的第i个亮度值m(i)是否处于P(n-1)和P(n)的范围内(步骤46)。在步骤46中为“是”的情况下,在步骤48中加上亮度值的增减量(偏置值)O(n)(步骤48)。另一方面,在步骤46中为“否”的情况下,对区域选择用计数器n进行更新后返回到步骤45(步骤47)。即、通过步骤45~步骤47的循环,被读入的三维断层像体数据的第i个亮度值m(i)被分类于由边界值P(0)~P(N)划分的任意一个区域,在步骤48中,加上与所划分的区域相应的亮度值的增减量(偏置值)O(n)。在步骤48结束后,或者在步骤45中为“否”的情况下,对计数器i进行更新后返回到步骤42(步骤49)。
通过以上的处理,在三维断层像体数据的所有体素的亮度值上,加上由边界值P(0)~P(N)分割出的每个区域的亮度值的增减量(偏置值)O(n)。另外,在本实施方式中,示出了利用偏置表来进行偏置处理的情况,但不限于此。例如,可以在偏置运算部中设定根据三维断层像体数据的各体素的亮度值的输入来输出亮度值的增减量的函数,基于该函数来对各体素的亮度值进行增减。
接着,对本实施方式的偏置运算部16向被检体1的应用例和效果进行说明。图5是用256个灰度的亮度值表示的胎儿的头部的断层图像的示意图。在断层图像51中,胎儿表皮52的亮度值为128。此外,胎儿头盖53的亮度值为200,头盖内腔区域54的亮度值为32。并且,脑实质55的亮度值为80,脑室56的亮度值为16。另外关于胎儿头部的周边的背景,为了简便地说明效果,而不予考虑。
在图6中示出在针对由这种断层图像的集合构成的三维断层像体数据,通过体绘制来构筑三维断层图像的情况下,将亮度较高的部位设定为不透明,将亮度较低的部位设定为透明这种一般的不透明度设定时的三维断层图像。如图6所示,关于三维断层图像61,构筑了由高亮度的胎儿表皮52、胎儿头盖53决定了大部分的输出亮度,绘制了胎儿头部的表面的三维表面图像。
对此,利用本实施方式的偏置运算部16,在想要提高作为检查者所希望的特定组织的脑室56的可视性的情况下,如图7所示那样来设定偏置表。即、偏置表成为图7(a)所示的图表71以及图7(b)所示的显示窗口72那样。此时,调整区域数N为2,P0、P2为输入亮度的最小、最大值且为固定。此外,为了只显示希望的特定目标部位,使得加在目标部位上的亮度值的增减量(偏置值)O1=255,加在目标部位以外的部位上的亮度值的增减量(偏置值)O2=-255,是最极端地将目标部位三维化的方法。边界值P1一边确认图像一边进行设定,因此能够由检查者通过控制面板26来设定。
图8是根据基于图7所示的偏置表通过偏置运算部16进行了偏置运算后的三维断层像体数据而生成的三维断层像的示意图。脑室56的亮度值为16,处于P0(0)~P1(17)的范围内,因此进行+255的偏置,其他组织处于P1(17)~P2(255)的范围内,因此进行-255的偏置。因此,如图8所示,三维断层图像81被生成并显示为提高了脑室56的可视性的三维断层图像。另外,在如本例这样,想要对从最小的亮度值到任意的亮度值进行强调显示的情况下,通过利用切换器、编码器或液晶面板上的调整按钮等来任意地设定参数“显示阈值”、以及亮度值的增减量(偏置值),能够以简便的操作生成图8所示的三维断层图像。
另一方面,在想要提高作为检查者所希望的特定组织的脑实质55的可视性的情况下,如图9所示那样来设定偏置表。即,如图5所示,脑实质55在内部存在低亮度的脑室56,在周围同样存在低亮度的头盖内腔区域54,并且在周围存在胎儿头盖53,因此想要提高可视性的亮度范围是处于中间亮度的一个区域。因此,偏置表如图9(a)所示的图表91以及图9(b)所示的显示窗口92那样,将调整区域数N设为3,P0、P3固定为输入亮度的最小、最大值,亮度值的增减量(偏置值)01=-255,02=255,O3=-255。边界值P1、P2可以一边确认图像一边通过控制面板26任意地设定。
图10是根据基于图9所示的偏置表通过偏置运算部16进行了偏置运算后的三维断层像体数据而生成的三维断层像的示意图。脑实质55的亮度值为80,处于P1(75)~P2(85)的范围内,因此进行+255的偏置,其他组织处于P0(0)~P1(75)或P2(85)~P3(255)的范围内,因此进行-255的偏置。因此,如图10所示,三维断层图像101被生成并显示为提高了脑实质55的可视性的三维断层图像。另外,在如本例这样,想要对处于中间亮度的任意的亮度范围进行强调显示的情况下,通过利用切换器、编码器或液晶面板上的调整按钮等任意地设定参数“显示开始阈值”、“显示结束阈值”以及亮度值的增减量(偏置值),能够以简便的操作生成图10所示的三维断层图像。
另外,在存在2个以上的希望的特定组织(目标部位)的情况下,用切换器或编码器很难设定,因此优选直接编辑显示画面上的设定画面。在控制面板26到显示器上的液晶面板上,通过控制面板26上的轨迹球(trackball)操作、或显示部24上的触摸面板,能够设定亮度调整区域数N的增减、和边界值P0-Pn、偏置00-0n。此外,通过使用显示窗口32内的下拉菜单(上下滑动条(bar)),也能够通过控制面板26来进行设定。
此外,在本实施方式中示出了检查者通过控制面板26来调整偏置表的例子,但不限于此。例如,也可以默认地准备多个偏置表如提高亮度值较低的组织的可视性的偏置表、提高亮度值处于中间的组织的可视性的偏置表、提高亮度值较高的组织的可视性的偏置表等,让检查者选择使用哪个偏置表。由此,检查者能够以简便的操作,得到提高了所希望的特定组织的可视性的三维断层图像。
接下来,利用图11~图13,对作为设定偏置表时的辅助画面的二维断层图像进行说明。图11是表示二维断层图像以及亮度-不透明度映射图的显示例的图。图11(a)示出二维断层图像,图11(b)示出亮度-不透明度映射图。如图11(b)所示,亮度-不透明度映射图111设定有使用了在纵轴、横轴分别取亮度值、不透明度的二维表的颜色代码。即、根据亮度值以及透明度/不透明度而设定了颜色代码。此外,如图11(a)所示,二维断层图像112根据亮度-不透明度映射图111的颜色代码而被着色。亮度-不透明度映射图111中的虚线相当于示出了图2所示的体素值和不透明度的关系的线,输入亮度值对应的不透明度是唯一确定的,在二维断层图像112上显示的像素的颜色确定为处于亮度-不透明度映射图111中的虚线上的1点。
作为二维断层图像112的基础的断层像数据通过断层像二维坐标变换部20而生成。断层像数据是检查者通过控制面板26指定的三维断层像体数据的任意的剖面的断层像数据,没有通过偏置运算部16进行偏置处理。该断层像数据被输入到图像合成部22。图像合成部22通过控制部28取得偏置表的边界值、亮度的增减量(偏置值)、以及亮度-不透明度映射图111。
图像合成部22基于偏置表,对从断层像二维坐标变换部20输入的断层像数据的亮度值进行偏置,针对偏置后的断层像数据的各体素,基于亮度-不透明度映射图111来变换颜色,生成二维断层图像。例如,参照亮度值被进行偏置后的断层像数据的各体素的亮度值、不透明度来变换为具有RGB或YUV等颜色信息的数值。在这种颜色代码变换处理中,通过像亮度-不透明度映射图111那样在纵轴设定亮度值、在横轴设定不透明度越低则与输入数据无关的彩度越高的颜色数据,能够分级显示在三维断层图像中透过的区域、不透过的区域、其中间的区域。
例如,在图11中,在二维断层图像112中,用与亮度不同的色调来对如下区域进行编码并显示:对三维断层图像起支配作用的不透明区域113;对三维断层图像产生很大影响的透明度较低的中间区域114;对三维断层图像不太产生影响的透明度较高的中间区域115;对三维断层图像不产生影响的透明区域116。其结果,检查者能够直观地掌握边界值以及亮度值的增减量(偏置值)对三维断层图像产生的影响,能够提高检查效率。即,检查者能够掌握二维断层图像112上的哪个区域被反映到三维断层图像中,哪个部位没有被反映。例如,因为透明区域116的区域在当前的偏置表的设定中没有被反映到三维断层图像中,因此可知,若想要使其反映到三维断层图像中,则必须调整偏置表的设定。
图12是表示三维断层图像、二维断层图像、和亮度-不透明度映射图的显示的一例的图。如图12所示,三维断层图像121、二维断层图像122、亮度-不透明度映射图123被并列显示。检查者通过参照图12,可以了解不透明区域113被描绘为不透明,透明区域116完全没有被描绘,中间区域114、115被描绘为半透明。此外,通过强调不透明区域113,能够清楚地将目标部位可视化。
图13是表示三维断层图像、二维断层图像、和亮度-不透明度映射图的显示的一例的图。如图13所示,也可以将三维断层图像131、X-Y断层面上的二维断层图像132、X-Z断层面上的二维断层图像133、Y-Z断层面上的二维断层图像134、和亮度-不透明度映射图135并列显示。通过像这样同时显示来自3个方向的剖面,能够更明显地使检查者了解可视化区域。
接下来,利用图14~图16,说明在上述颜色代码变换处理中使用的颜色变换表(亮度-不透明度映射图)的作成方法的例子。图14是表示颜色变换表的作成流程的图。在图14中记载的记号P、O是指上述的边界值P(n)、亮度值的增减量(偏置值)O(n)。
此外,在图15、16中,在本实施方式中使用的颜色变换表OUTMAP[i]表示合成后颜色映射表,BWMAP[i]表示断层像输出用RGB变换表,COLMAP[i]表示透明度设定用RGB颜色表。它们是针对输入i,通过还原作为三原色R、G、B的三要素来进行向自然色的变换的表,分别是能够用*.R,*.G,*.B来访问RGB3要素的颜色映射表。i表示基于输入亮度的灰度,多为0-255的256个灰度,但可以使用任意的数。
透明设定用RGB表COLMAP例如可以对所有的灰度要素i设定*.R=0,*.G=255,*.B=0,且具有随着变为透明而变为灰色的设定,也可以通过对灰度要素i设定不同的色调、彩度、明度来提供任意的视觉效果。
如图14所示,开始颜色映射变换处理后,由计数器i来表现的输入灰度值(表示输入亮度值、弹性值、或其他输入)被初始化为0(步骤141)。接下来,进行计数器i和Mmap(映射图所具有的灰度数)的大小比较(步骤142),若i小于Mmap(在步骤142中为“否”),则进入到步骤143。在步骤143中,区域选择用计数器n被初始化为1,所选择的输入灰度值i通过步骤144以及步骤145,被分类于由边界值P(0)~P(N)分割出的任意一个区域,在步骤146中利用将所选择的亮度值的增减量(偏置值)O(n)加i后所得的值来参照不透明度表OPQ[n],求出不透明度mm(i),利用该不透明度mm(i)在步骤148中作成合成后颜色映射表OUTMAP[i]。
在步骤146中所使用的常数v是提供不透明度颜色的影响度的参数,通常为1.0。步骤143、步骤147是区域选择用计数器n的初始化处理以及更新处理,步骤149是计数器i的更新处理。
图15是表示在步骤148中作成合成后颜色映射表OUTMAP[i]的处理的一例的图。如图15所示,步骤151计算出对BWMAP[i]、OUTMAP[i]分别乘以系数mm(i)、(1.0-mm(i))并相加后得到的值作为OUTMAP[i]的各要素。
这表示作成了如下的颜色映射表OUTMAP[i]:例如在透明设定用RGB表COLMAP针对所有的灰度要素i,设定了*.R=0,*.G=255,*.B=0的情况下,在具有不透明的输入灰度值i的情况下,直接输出断层像亮度BWMAP[i],在具有透明的输入灰度值i的情况下,随着不透明度的减少,逐渐输出绿色的断层像图像。
即,作为不透明的三维体素而被图像化的断层像的区域直接显示黑白、或者所设定的色调的断层像,作为透明的三维体素而被图像化的断层像的区域用所设定的色调以外的透明度设定用颜色来表示显示,随着从不透明变为透明,透明度设定用颜色的比例较强地显示在断层像上。
此时,在完全不透明的输入灰度值i的情况下,完全发出绿色,但通过将在步骤146中使用的常数v设为1.0以下(例如0.9或0.8),能够调整半透明等发色的状态。由此,完全不透明的输入灰度值i的区域并非成为完全的绿色而是成为半透明的绿色,因此检查者能够看到该区域的形态等。
图16是表示在步骤148中作成合成后颜色映射表OUTMAP[i]的处理的一例的图。在步骤161中,不准备透明度设定用RGB颜色表,而是利用不透明度系数mm(i)来计算与固定的透明度设定用颜色值的合成色,并简易地作成合成后颜色映射表OUTMAP。
上述的两例是利用不透明度系数mm(i)来合成RGB颜色映射图的各要素的方法的一例,但也可以采用保持彩度、明度、色调作为参数的HSV形式来保持颜色信息,并利用通过mm(i)而使这些参数变化后的信息来进行合成,输出形式不限定于RGB形式,也可以为YUV方式。
(第二实施方式)
参照附图来说明应用本发明而成的超声波诊断装置的第2实施方式。图17是表示第2实施方式的超声波诊断装置的构成的框图。对于与第1实施方式同样的构成省略说明。
如图17所示,在超声波诊断装置100中具备:存储从相位累加部10输出的RF信号帧数据,并至少选择2枚帧数据的RF信号帧数据选择部171;测量被检体的生物体组织的位移的位移测量部172;根据由位移测量部172测量出的位移信息来求出变形或弹性模量的弹性信息运算部173;根据由弹性信息运算部173运算出的变形或弹性模量来构成彩色弹性图像的弹性图像处理部174;和将来自弹性图像处理部174的输出信号按照适合显示部24的显示的方式进行变换的弹性像二维坐标变换部177。
此外,还具备:基于来自弹性图像处理部174的输出信号来进行向三维弹性数据的坐标变换,并作成弹性像体数据的弹性像体数据生成部175;和基于来自弹性像体数据生成部175的输出信号,对存在于二维投影面上的各像素的视线方向上的弹性像体数据,进行体绘制、最大值、最小值投影、或者平均处理的弹性像体绘制部176。
RF信号帧数据选择部171保存来自相位累加部10的多个RF信号帧数据,并从所保存的RF信号帧数据群中选择1组即2个RF信号帧数据。例如,从相位累加部10将基于时间序列、即图像的帧速率而生成的RF信号帧数据依次存储到RF信号帧数据选择部171中,并选择所存储的RF信号帧数据(N)作为第1数据,同时,从在时间上过去存储的RF信号帧数据群(N-1,N-2,N-3...N-M)中选择1个RF信号帧数据(X)。另外,在此N、M、X是附加给RF信号帧数据的索引编号,是自然数。
然后,位移测量部172根据所选择的1组数据、即RF信号帧数据(N)以及RF信号帧数据(X)来进行一维或二维相关处理,并求出与断层图像的各点对应的生物体组织上的位移或移动向量即位移的方向和大小所相关的一维或二维位移分布。在此,在移动向量的检测中,使用块匹配(block matching)法。块匹配法进行如下处理:将图像分为例如由N×N个像素构成的模块,着眼于关心区域内的模块,从前面的帧中搜索与着眼的模块最近似的模块,并参照此模块来进行预测编码即根据差分决定采样值。
弹性信息运算部173根据从位移测量部172输出的测量值例如移动向量、和从压力测量部178输出的压力值来计算与断层图像上的各点对应的生物体组织的变形或弹性模量,并基于该变形或弹性模量来生成弹性图像信号即弹性帧数据。
此时,变形的数据通过对生物体组织的移动量例如位移进行空间微分来计算。此外,在具备压力测量部178所示的压力测量功能的构成中,能够推算弹性模量,也可以使用弹性模量作为弹性数据。弹性模量的数据是通过将压力的变化除以变形的变化来计算的。例如,若假设由位移测量部172测量出的位移为L(X),由压力测量部178测量出的压力为P(X),则因为变形ΔS(X)能够通过对L(X)进行空间微分来算出,所以利用ΔS(X)=ΔL(X)/ΔX这个式子能够求出。此外,弹性模量数据的杨氏模量Ym(X)通过Ym=(ΔP(X))/ΔS(X)这个式子来计算。根据该杨氏模量Ym能够求出与断层图像的各点相当的生物体组织的弹性模量,因此能够连续地得到二维的弹性图像数据。另外,杨氏模量是指,施加于物体上的单纯拉伸应力相对于与拉伸平行地产生的变形之比。
弹性图像处理部174包括帧存储器和图像处理部而构成,将从弹性信息运算部173按照时间序列输出的弹性帧数据确保在帧存储器中,并对所确保的帧数据进行图像处理。此外,在弹性图像处理部174中根据RF信号帧数据选择部171或位移测量部172或弹性信息运算部173的输出信息来评价弹性图像的错误,并进行输出图像的遮盖(masking)。
弹性像二维坐标变换部177根据监视器来对来自弹性图像处理部174的弹性帧数据进行坐标变换。另外,超声波探头2能够根据来自控制部28的控制信号以电动机驱动、或手动方式在短轴方向上扫描。此外,即使为手动方式,在具备磁传感器179的构成的情况下,通过使用来自磁传感器179的位置信息,也能够检测压迫量和短轴位置。
弹性像体数据生成部175通过弹性图像处理部174来进行向三维弹性数据的坐标变换,弹性像体绘制部176对存在于所输出的二维投影面上的各像素的视线方向上的体数据,进行体绘制、最大值、最小值投影、或者平均处理。
图像合成部22对由断层像以及弹性像体绘制部作成的断层数据和变形/弹性数据进行合成。合成图像的各像素的亮度信息以及色调信息,将黑白断层图像和彩色弹性图像的各信息按照合成比例进行相加,并进行RGB变换来作成显示于显示部24的图像。
在此,偏置运算部16将由断层像体数据生成部14输出的三维断层像体数据利用三维弹性像数据进行亮度变换,并传送到断层像体绘制部18。换言之,本实施方式的偏置运算部16,根据从弹性像体数据生成部175输出的三维弹性像体数据的各体素的弹性值,来对从断层像体数据生成部14输出的三维断层像体数据的对应的各体素的亮度值进行增减。
另外,弹性像数据是通过上述方法计算出的变形值、杨氏模量、纵弹性模量、横弹性模量等表示硬度的弹性参数的总称,不是指特定的1个参数。
接下来,对本实施方式的超声波诊断装置的特征部分进行说明。图18是表示在偏置运算部16中使用的偏置表的设定画面的一例的图。图18(a)(b)分别是对偏置表进行设定的一例。
如图18所示,偏置表设置有将三维弹性像体数据的各体素可能取得的弹性值的范围(0~255)划分为多个区域(N:在本实施方式中为5)的边界值(E0~E5)、和所划分出的多个区域中的每一个的亮度值的增减量(O1~O5)。
例如,利用图18(a)所示的图表181,通过控制面板26将E0~E5所表示的边界值的设定指针向左右移动,由此能够设定区域的边界。同样,01~05是表示各区域的亮度值的增减量(偏置值)的指针,通过向上下移动能够设定各区域的偏置值。图表181的上限、下限值分别为100,-100,但是这是示出的一例,只要在输入以及输出数据的可能取得的范围内则为多少都可以。另外,移动E0或E5,对于成为偏置的对象外的输入弹性值不进行偏置处理。此外,E0以及E5也可以固定于各体素可能取得的弹性值的范围的两端。此外,边界值(E0~E5)、亮度值的增减量(O1~O5)的设定,可以利用切换器、编码器或液晶面板上的调整按钮等来进行。
此外,如图18(b)所示,边界值以及亮度值的增减量(偏置值),也可以经由控制面板26,通过下拉菜单或直接数值输入来分别设定显示窗口182内的01~05、E0~E5。此外,如图18(b)所示,也可以将通过边界值(E0~E5)划分出的有效的区域数作为亮度调整区域数183(N=5)来进行显示。
图19是在偏置运算部16中进行的亮度值的变换处理的数据流程图。如图19所示,在偏置运算部16中开始偏置调整处理后,计数器i被初始化为0(步骤191)。接下来,读入三维断层像体数据的各体素的亮度值m(i)以及三维弹性像体数据的各体素的弹性值e(i)(步骤192)。判定i是否大于全部数据数(体素数)(步骤193),在大于的情况(在步骤193中为“是”的情况)下结束偏置调整处理。另一方面,在不大于、即没有对三维断层像体数据的所有的体素进行完偏置处理的情况下,进入到步骤194,区域选择用计数器n被初始化为1(步骤194)。
判断区域选择用计数器n是否小于亮度调整区域数N(步骤195),在小于的情况下(在步骤195中为“是”),判定所读入的三维弹性像体数据的第i个弹性值e(i)是否处于E(n-1)和E(n)的范围内(步骤196)。在步骤196中为“是”的情况下,在步骤198中亮度值的增减量(偏置值)O(n)被加在三维断层像体数据的对应的体素的亮度值上(步骤198)。另一方面,在步骤196中为“否”的情况下,更新区域选择用计数器n后返回到步骤195(步骤197)。
即、通过步骤195~步骤197的循环,三维弹性像体数据的被读入的第i个弹性值e(i)被分类于由边界值E(0)~E(N)划分出的任意一个区域,在步骤198中,与划分出的区域相应的亮度值的增减量(偏置值)O(n),被加在三维断层像体数据的对应的体素的亮度值上。在步骤198结束后、或者在步骤195中为“否”的情况下,更新计数器i后返回到步骤192(步骤199)。
通过以上的处理,在三维断层像体数据的所有体素的亮度值上,加上了由边界值E(0)~E(N)分割出的每个区域的亮度值的增减量(偏置值)O(n)。另外,在本实施方式中示出了利用偏置表来进行偏置处理的情况,但不限于此。例如,可以在偏置运算部中,设定根据三维弹性像体数据的各体素的弹性值的输入,来输出三维断层像体数据的对应的体素的亮度值的增减量的函数,并基于该函数来对三维断层像体数据的对应的体素的亮度值进行增减。
图像合成部22与第1实施方式相同,基于来自断层像二维坐标变换部20的输出和来自断层像体绘制部18的输出,来显示二维断层图像和三维断层图像。在本实施方式中,除此之外,还可以显示在图像合成部22中基于来自弹性像体绘制部176的输出,将三维弹性像体数据投影于二维投影面后得到的三维弹性图像。
根据本实施方式,检查者能够着眼于所希望的特定组织的弹性值(组织的硬度或柔软度),来生成提高了具有特定弹性值的组织的可视性的三维断层图像。例如,只要掌握了所希望的特定组织为怎样程度的弹性值,则能够预先增大与弹性值附近对应的亮度值的增量(偏置量)。另一方面,在想要观察硬(或软)组织的情况下,只要增大与高(或低)弹性值对应的亮度值的增量(偏置量),就会显示硬(或软)组织的可视性较高的三维断层图像,因此硬(或软)组织的观察变得容易。
(第三实施方式)
参照附图来说明应用本发明而成的超声波诊断装置的第3实施方式。图20是表示第3实施方式的超声波诊断装置的构成的框图。如图20所示,本实施方式相对于第2实施方式的超声波诊断装置,在如下这一点上不同:将弹性像体数据生成部175的输出输入到断层像体绘制部18,使得断层像体绘制部18能够参照三维弹性像体数据。省略其他的与第2实施方式同样的部分的说明。本实施方式具有利用三维断层像体数据和三维弹性像体数据这两者来进行检查者所需要的三维断层图像的生成的功能。
在本实施方式中,偏置运算部16与第1实施方式相同地按照根据三维断层像体数据的各体素的亮度值来对该各体素的亮度值进行增减的方式构成。
此外,控制面板26具有设定弹性值的阈值的功能,通过控制面板26设定的弹性值阈值信息经由控制部28被设定于断层像体绘制部18。本实施方式中的断层像体绘制部18,在利用前述的式(1)、(2)来进行断层像的体绘制时,在各体素的弹性值超过弹性值阈值(或者低于)的情况下,针对该体素将式(1)中的不透明度Ai设为零。
由此,即使输入亮度值通过偏置运算处理进行了亮度变换,并对照了不透明度表,结果成为了不透明、即应被显示的亮度值,但在弹性值不处于目标区域中的情况下,也将透明化。其结果,能够仅对处于目标弹性值、和目标亮度值的范围内的体素数据(亮度数据)进行可视化。
图21是通过本实施方式的处理,仅将处于目标弹性值、和目标亮度值的范围内的体素数据(亮度数据)可视化的情况下的显示例。如图21所示,能够并列显示如下图像:仅将处于目标弹性值、和目标亮度值的范围内的体素数据(亮度数据)可视化后得到的三维断层图像211、和重叠了三维断层图像211的任意剖面上的断层图像和弹性图像的图像212~213。
这些图像212~214由断层像二维坐标变换部20、弹性像二维坐标变换部177生成,在图像合成部22中进行颜色代码化处理后,利用α混合等手法来重叠。另外,这些图像212~213不仅能够如图示那样只同时显示3枚,而且还能够同时显示任意数量的切片图像。
弹性像体绘制部176能够利用通过控制面板26设定的弹性值阈值,将想要显示的范围的弹性值设定为不透明度1.0,将不想显示的范围的不透明度设定为0.0,从而仅将有效显示范围的图像的表面可视化。此外,能够参照断层像体数据生成部14的输入亮度,将输入亮度的不透明度不为零的区域的弹性值设定为不透明度1.0,将输入亮度的不透明度为零的范围的不透明度设定为0.0,从而仅将有效显示范围的图像的表面可视化。
图22是通过本实施方式的处理,仅将处于目标弹性值、和目标亮度值的范围内的体素数据(亮度数据)可视化的情况下的显示例。如图22所示,能够并列显示三维断层弹性图像221、和重叠了三维断层弹性图像221的任意剖面上的断层图像和弹性图像的图像222~224。三维断层弹性图像221是在图像合成部22中利用α混合等手法将利用上述方法将图像表面的弹性值三维化后得到的图像、和将处于目标弹性值和目标亮度值的范围内的体素数据(亮度数据)三维化后得到的图像重叠后的图像。与三维断层图像211相比,进行了弹性图像的着色,因此能够直观地判别与硬度相关的信息。
图23是表示利用了在上述倒数第5段中说明了的弹性值的阈值的三维断层图像的显示例的图。此时,在实施例中示出了虽具备偏置运算处理,但没有进行基于亮度的偏置调整的设定的情况。图23(a)、(b)、(c)分别表示:无阈值(全部显示)的情况下的三维断层图像231、将阈值设为中间值(使一部分软组织不显示)的情况下的三维断层图像232、将阈值设为高值(只显示硬部组织)的情况下的三维断层图像233。
例如,在三维断层图像231中,形成于影像(phantom)内部的内部呈低亮度的圆筒状的硬部组织的进深方向的形状难以看清。对此,通过像三维断层图像232、233那样变更弹性值的阈值的设定,能够逐渐使得形成于影像内部的圆筒状的低亮度区域容易看清。
图24是表示基于弹性值的阈值来改变想要显示的范围的弹性值的情况下的三维弹性图像的显示例的图。图24(a)、(b)、(c)分别表示:无阈值(全部显示)的情况下的三维弹性图像241、将阈值设为中间值(使一部分软组织不显示)的情况下的三维弹性图像242、将阈值设为高值(只显示硬部组织)的情况下的三维弹性图像243。
本例的情况也是同样,在三维弹性图像241中,形成于影像内部的圆筒状的低亮度区域难以看清,但通过像三维弹性图像242、243那样变更弹性值的阈值的设定,能够逐渐使得形成于影像内部的圆筒状的低亮度区域容易看清。
图25是将图23、图24所示的三维断层图像和三维弹性图像的RGB变换后的RGB各要素利用加权相加的α混合而重叠后的图像例。即,图25(a)、(b)、(c)分别表示:无阈值的(全部显示)的情况下的三维断层弹性图像251;将阈值设为中间值(使一部分软组织不显示)的情况下的三维断层弹性图像252;将阈值设为高值(只显示硬部组织)的情况下的三维断层弹性图像253。
本例的情况也是同样,在三维断层弹性图像251中,形成于影像内部的圆筒状的低亮度区域难以看清,但通过像三维断层弹性图像252、253那样来变更弹性值的阈值的设定,能够逐渐使得形成于影像内部的圆筒状的低亮度区域容易看清。此外,通过重叠显示三维断层图像和三维弹性图像,能够对检查者提供更加适合诊断的图像。
接下来,在图26、图27中示出利用在上述段落中说明了的弹性值的阈值处理,并且进行了基于亮度的偏置调整的设定的情况下的图像显示例。对于偏置,按照强调低亮度的目标观察部位,抑制除此之外的高亮度部位的方式进行了设定。
图26是表示利用在上述段落中说明了的弹性值的阈值,并且进行了基于亮度的偏置调整的设定的情况下的三维断层图像的显示例的图。图26(a)、(b)、(c)分别表示:无阈值(全部显示)的情况下的三维断层图像261;将阈值设为中间值(使一部分软组织不显示)的情况下的三维断层图像262;将阈值设为高值(只显示硬部组织)的情况下的三维断层图像263。
例如,在三维断层图像261中,形成于影像内部的内部呈低亮度的圆筒状的硬部组织的进深方向的形状难以看清。对此,通过像三维断层图像262、263那样变更弹性值的阈值的设定,能够逐渐使得形成于影像内部的圆筒状的低亮度区域容易看清。
图27是将图26所示的三维断层图像、和基于弹性值的阈值并且进行了亮度的偏置调整的设定的情况下的三维弹性图像的RGB变换后的RGB各要素,利用加权相加的α混合而重叠后的图像例。即,图27(a)、(b)、(c)分别表示:无阈值(全部显示)的情况下的三维断层弹性图像271;将阈值设为中间值(使一部分软组织不显示)的情况下的三维断层弹性图像272;将阈值设为高值(只显示硬部组织)的情况下的三维断层弹性图像273。
本例的情况也是同样,在三维断层弹性图像271中,形成于影像内部的圆筒状的低亮度区域不易看清,但通过像三维断层弹性图像272、273那样变更弹性值的阈值的设定,能够逐渐使得形成于影像内部的圆筒状的低亮度区域容易看清。此外,通过将三维断层图像和三维弹性图像重叠显示,能够对检查者提供更加适合诊断的图像。
此外,如图26、图27所示,在基于弹性值的阈值处理并且进行了基于亮度的偏置调整的设定的情况下,通过偏置处理,低亮度的目标观察部位被强调并高亮度地显示,抑制了高亮度的周边部位的一部分。另一方面,阈值处理按照显示硬部组织、除去软组织的方式进行了设定,硬的目标观察部位被显示,柔软的周边部位的一部分被除去。结果,用偏置处理以及阈值处理的任意一方无法除去的无用的部位从三维图像中被除去,与图23、图25相比能够更加强调目标观察部位地进行显示。
另外,上述各实施方式,主要对超声波诊断装置、以及超声波图像生成方法进行了说明,但本发明不限定于此。例如,可以应用于针对预先在超声波诊断装置中生成的被检体的三维断层像体数据脱机(off line)地生成三维断层图像的PC等超声波图像处理装置。此外,也可以将本发明作为嵌入到PC等医用图像处理装置中的超声波图像处理程序。
即,应用本发明而成的PC等超声波图像处理装置,可以具备如下部件而构成:保存基于由超声波探头测量到的被检体的多个断层面的反射回声信号而生成的三维断层像体数据的存储器;基于在该存储器中保存的三维断层像体数据,在二维投影面上生成至少从一个视线方向来看的三维断层图像的断层像体绘制部;和显示三维断层图像的显示部。
并且,可以构成为,具备根据三维断层像体数据的各体素的亮度值来对该各体素的亮度值进行增减的偏置运算部,偏置运算部的各体素的亮度值的增减量能够通过输入接口来调整,断层像体绘制部基于通过偏置运算部对亮度值进行了偏置后的三维断层像体数据来生成三维断层图像。
此外,在该超声波图像处理装置中,可以构成为,在存储器中保存了基于被检体的多个断层面的反射回声信号而生成的三维弹性像体数据的情况下,偏置运算部根据三维弹性像体数据的各体素的弹性值,来对三维断层像体数据的对应的各体素的亮度值进行增减。
检查者例如将由超声波诊断装置等生成的被检体的三维断层像体数据保存在USB、CD-ROM等信息记录介质中,并经由图像输入部输入到PC等超声波图像处理装置。或者也可以不利用信息记录介质而经由网络来输入被检体的三维断层像体数据。
于是,超声波图像处理装置,作为超声波图像处理程序,执行如下步骤:根据三维断层像体数据的各体素的亮度值来对该各体素的亮度值进行增减的步骤;基于亮度值被进行偏置后的三维断层像体数据,在二维投影面上生成至少从一个视线方向来看的三维断层图像的步骤;和显示生成的三维断层图像的步骤。
此外,在超声波图像处理程序中嵌入了基于被检体的多个断层面的反射回声信号来生成三维弹性像体数据的步骤的情况下,对各体素的亮度值进行增减的步骤,可以按照根据三维弹性像体数据的各体素的弹性值,来对三维断层像体数据的对应的各体素的亮度值进行增减的方式构成。
由此,检查者能够脱机地通过超声波处理装置的输入接口来对偏置运算部的各体素的亮度值的增减量进行调整,因此能够生成只强调了具有特定亮度值的组织的三维断层图像。因此,通过对所希望的特定组织的亮度值进行增量调整,或者对所希望的特定组织以外的组织的亮度值进行减量调整,或者进行这两者的调整,能够生成提高了所希望的特定组织的可视性的三维断层图像。例如,只要掌握了所希望的特定组织为怎样程度的亮度值,便能够预先增大与该亮度值附近对应的亮度值的增量(偏置量)。此外,例如在观察所生成的三维断层图像时,所希望的特定组织隐藏在高亮度组织的背面从而难以看清的情况下,只要增大与高亮度值附近对应的亮度值的减量(偏置量),则成为障碍的高亮度组织不易反映到三维断层图像中,因此能够提高所希望的特定组织的可视性。
此外,检查者能够着眼于所希望的特定组织的弹性值(组织的硬度或柔软度),脱机地通过超声波处理装置的输入接口来对偏置运算部的各体素的亮度值的增减量进行调整,因此能够生成提高了具有特定弹性值的组织的可视性的三维断层图像。例如,只要掌握了所希望的特定组织为怎样程度的弹性值,则能够预先增大与该弹性值附近对应的亮度值的增量(偏置量)。另一方面,在想要观察硬(或软)组织的情况下,只要增大与高(或低)弹性值对应的亮度值的增量(偏置量),则能够搜索并观察硬(或软)组织。
(符号说明)
1被检体;2超声波探头;4发送部;6接收部;8超声波收发控制部;10相位累加部;12断层图像处理部;14断层像体数据生成部;16偏置运算部;18断层像体绘制部;20断层像二维坐标变换部;22图像合成部;24显示部;26控制面板;28控制部;111、123、135亮度-不透明度映射图;171RF信号帧数据选择部;172位移测量部;173弹性信息运算部;174弹性图像处理部;175弹性像体数据生成部;176弹性像体绘制部;177弹性像二维坐标变换部