具体实施方式
下面,根据实施方式来说明本发明。图1中示出本发明一实施方式的超声波诊断装置的结构框图。图2中示出本发明特征部分的超声波诊断装置之处理步骤的一实施方式的流程图。
如图1所示,抵接于被检体1使用的超声波探头2如图3(A)所示,被形成具有排列了与被检体1之间发送及接收超声波的多个振子的超声波发送接收面21。探头2由从发送电路3供给的超声波脉冲驱动,机械式或电子地进行束扫描。发送接收控制电路4控制用于驱动探头2的多个振子的超声波脉冲的发送计时,向设定在被检体1内的焦点形成超声波束。并且,发送接收控制电路4沿探头2的振子的排列方向电子地扫描超声波束。
另一方面,探头2接收从被检体1内产生的反射回声信号,输出至接收电路5。接收电路5根据从发送接收控制电路4输入的计时信号,进行取入反射回声信号后放大等接收处理。由接收电路5接收处理后的反射回声信号控制在相位调整加法电路6中由多个振子接收到的反射回声信号的相位,对一个或多个收敛点形成超声波受波束。将在相位调整加法电路6中相位调整相加后的反射回声信号(下面称为超声波断层数据)输入到信号处理部7,进行增益补正、记录压缩、检波、轮廓强调、滤波处理等信号处理。另外,在相位调整加法电路6中生成的超声波断层数据的高频(RF)信号也可是合成解调后的I、Q信号。这样,由探头2在被检体1体内按一定方向扫描超声波束,得到对应于1个断层像的超声波断层数据。
将由信号处理部7处理后的超声波断层数据导入黑白扫描转换器8,在其中转换成数字信号,同时,转换成对应于超声波束的扫描面的2维断层像数据。即,黑白扫描转换器8是为了以超声波周期取得包含运动组织的被检体1内的RF信号帧数据和将该帧数据转换成图像并显示,而用于视频同步读出的断层扫描部件及进行系统控制用的部件,由如下构成:即,将来自信号处理部7的反射回声信号转换成数字信号的A/D转换器;按时间序列存储由该A/D转换器数字化后的多个断层像数据的帧存储器;和控制这些动作的控制器等。由该信号处理部7和黑白扫描转换器8构成断层像的图像再构成部件。从黑白扫描转换器8输出的断层像数据经切换加法部9供给图像显示器10,显示断层像。
图像显示器10显示由黑白扫描转换器8得到的时间序列的断层像数据,由将经切换加法器9输入的图像数据转换成模拟信号的D/A转换器、和输入来自该D/A转换器的模拟视频信号并显示为图像的彩色视频监视器构成。
另一方面,将从相位调整加法电路6输出的超声波断层数据导入RF信号帧数据选择部11。RF信号帧数据选择部11将超声波束的扫描面(断层面)所对应的RF信号群作为帧数据,选择多个帧大小,存储在存储器等中。位移测量部12依次取入存储在RF信号帧数据取得部11中的、取得时刻不同的多对帧数据,根据取入的一对帧数据,求出断层面中多个测量点的位移矢量,作为位移帧数据,输出至失真/弹性模数运算部13。
失真/弹性模数运算部13根据输入的位移帧数据,求出断层面中多个测量点的失真,作为弹性帧数据,输出至弹性数据处理部14。并且,失真/弹性模数运算部13从压力测量部19取入施加于被检体的压力测量数据,求出被检体各部分的应力分布,根据先求出的失真帧数据和应力分布求出弹性模数,作为弹性帧数据,输出至弹性数据处理部14。压力测量部19如图3(c)所示,取入安装在压迫板22等上的多个压力传感器23的输出信号,测量施加于被检体1体表面的压力,将测量出的压力数据发送至失真/弹性模数运算部13。另外,压力测量如图3(D)所示,可使用安装了在压迫板22的前面安装的参照变形体24的探头2,根据参照变形体24的位移帧数据,测量施加于体表面的压力。该压力测量法在特开2005-13283号公报或特开2005-66041号公报中记述。
弹性数据处理部14对从失真/弹性模数运算部13输入的失真或弹性模数的弹性帧数据实施坐标平面内的平滑处理、对比度最佳化处理或帧间的时间轴方向的平滑处理等各种图像处理,发送至彩色扫描转换器15。
彩色扫描转换器15转换从弹性数据处理部14输出的弹性帧数据,生成彩色的弹性图像,经切换加法部9,在图像显示器10中显示。即,彩色扫描转换器15根据预先设定的弹性(失真或弹性模数)的上限值及下限值的范围,对弹性图像赋予灰度化(例如,256灰度)后的红、绿、蓝等色调代码。例如,弹性帧数据的弹性模数测量得大的硬区域转换成蓝色代码,相反,弹性模数测量得小的柔软区域转换成红色代码。另外,可使用黑白扫描转换器替代彩色扫描转换器15。这时,可通过弹性模数测量得大的硬区域使亮度变亮,相反,弹性模数测量得小的柔软区域使亮度变暗等,来表示弹性模数的分布。
并且,切换加法部9具有:输入从黑白扫描转换器8输出的黑白断层像数据和从彩色扫描转换器15输出的彩色弹性图像数据、切换两个图像和显示任一个的功能;使两图像的一个成为半透明、相加合成后在图像显示器10中重叠显示的功能;和并排显示两个图像的功能。并且,电影存储器部18将从切换加法部9输出的图像数据存储在存储器中,根据来自装置控制接口部17的指令,调出过去的图像数据,在图像显示器10中显示。并且,可将选择出的图像数据传输至MO等记录媒体。
并且,作为本实施方式的特征的直方图评价部20,取入从位移测量部12和失真/弹性模数运算部13输出的位移帧数据或失真帧数据等弹性帧数据,如后所述,在生成直方图的同时,根据直方图来评价统计特征。而且,生成直方图图像数据,经切换加法部9,在图像显示器10中显示直方图图像。
参照图2中示出的流程图来说明这样构成的本实施方式的详细结构及动作。
(步骤S1)
首先,边对被检体1的组织施加压力,边连续地取得超声波断层数据(S1)。这时,使图3(A)示出的探头2的超声波发送接收面21接触被检体1的体表进行超声波发送接收,同时,用手上下移动超声波发送接收面21,压迫被检体1。这时,为了对被检体1的诊断部位有效地供给应力分布,如图3(B)所示,优选地,使面与探头2的超声波发送接收面21吻合而安装压迫板22,用由超声波发送接收面21和压迫板22构成的压迫面压迫被检体1。
(步骤S2)
RF信号帧数据选择部11在帧存储器内依次保存从相位调整加法电路6以超声波诊断装置的帧速率按时间序列输出的超声波断层数据(RF信号帧数据)。这时,保存最新的RF信号帧数据作为N。而且,根据控制命令,从时间上过去的RF信号帧数据N-1、N-2、N-3、…、N-M中选择1个RF信号帧数据X。接着,设最新的RF信号帧数据N和选择出的RF信号帧数据X为1组RF信号帧数据,输出至位移测量部12。
(步骤S3)
位移测量部12对选择出的1组RF信号帧数据N、X进行1维或2维相关处理,测量断层像上各测量点的位移或位移矢量(位移的方向和大小),生成位移帧数据。作为该位移矢量的检测法,例如如特开平5-317313号公报所述,具有块匹配法和梯度(gradient)法。块匹配法是将图像分成例如由n×n像素构成的多个块,从在前帧搜索最近似于当前帧中着眼的块的块,求出位移矢量。
(步骤S4)
失真/弹性模数运算部13根据由位移测量部12及压力测量部19求出的位移帧数据及压力数据,运算断层像上各测量点的失真及弹性模数。然后,生成对应断层像的失真或弹性模数的弹性帧数据,输出至弹性数据处理部14。这里,失真不需压力数据,通过将位移进行空间微分来计算。并且,作为弹性模数的一种,例如杨氏模数Ym的运算,如下式所示,通过根据压力数据求出各测量点的应力(压力),用各测量点的失真量除该应力来求出。这里,i、j的指标表示帧数据的测量点的坐标。
Ymi,j=压力(应力)ij/(失真量i,j)
这里,(i,j=1,2,3,…)
(步骤S5)
弹性数据处理部14对从失真/弹性模数运算部13输入的弹性帧数据实施坐标平面内的平滑处理、对比度最佳化处理、帧间的时间轴方向的平滑处理等图像处理,发送至彩色扫描转换器15。该弹性数据处理部14中的图像处理例如在特开2004-261198号公报中记载了。
彩色扫描转换器15输入从弹性数据处理部14输出的弹性帧数据、和用于指定弹性帧数据的灰度化选择范围的上限值及下限值,生成将弹性帧数据转换成色调信息的弹性图像数据。该上限值及下限值由来自未图示的装置控制部的命令、或来自弹性数据处理部14的弹性帧数据中的灰度化选择范围指定。例如,彩色扫描转换器15对于输入的弹性帧数据,将弹性模数测量得大的区域转换成蓝色代码,相反,将弹性模数测量得小的区域转换成红色代码。
这里,也可使用黑白扫描转换器替代彩色扫描转换器15。这时,失真大的区域可使亮度变亮,相反,失真小的区域使亮度变暗。
这样,将由彩色扫描转换器15生成的弹性图像数据输出至切换加法部9。切换加法部9相加或切换从黑白扫描转换器8输出的黑白的断层像数据、和从彩色扫描转换器15输出的彩色的弹性图像数据,输出至图像显示器10并显示。例如,仅显示黑白的断层像数据、或仅显示彩色的弹性图像数据,或者相加合成两图像数据来显示。并且,例如,如特开2000-60853号公报所述,可同时显示黑白断层像和彩色或黑白弹性图像2个画面。并且,例如,如特开2004-135929号公报所述,可将彩色的弹性图像半透明地重叠在黑白断层像上显示。
并且,在将从切换加法部9输出的显示图像数据发送至图像显示器10的同时,输入至电影存储器部18。电影存储器部18将输入的显示图像数据存储在存储器中。并且,电影存储器部18根据从装置控制接口部17供给的控制信号,调出过去的显示图像数据和在图像显示器10中显示,或将选择出的显示图像数据传输至MO等存储媒体。
(步骤S6、S7)
在步骤S6、S7中,执行本发明特征的直方图评价处理。直方图评价部20使用从失真/弹性模数运算部13输出的弹性帧数据,求出在当前时刻的关心区域内各测量点的失真或弹性模数,从而求出关心区域内的直方图。而且,生成反映了求出的直方图的图像信息,经切换加法部9在图像显示器10中显示,并作为诊断信息提供给检测者。在反映该直方图的图像信息中,可包括包含了着眼部位的关心区域的各部位的弹性的统计特征。
这里,参照图4说明直方图评价部20中的一系列处理的流程。直方图评价部20具备存储部31、直方图运算部32、统计评价部33和图像构筑部34。存储部31暂时存储从失真/弹性模数运算部13输出的弹性帧数据,输出至直方图运算部32。
直方图运算部32输入从存储部31输出的弹性帧数据、和从装置控制接口部17输出的关心区域(ROI)的坐标数据。这里,ROI在步骤S5中在图像显示器10显示出的彩色弹性图像上,从装置控制接口部17输入设定(S6)。
而且,直方图运算部32运算设ROI内分布的弹性帧数据为母集团的直方图,生成由数值数据构成的直方图数据,输出至统计评价部33。即,将ROI中的失真或弹性模数的值分割成多个区间,算出具有对应于各区间的失真或弹性模数的值的测量点个数,例如,将失真或弹性模数的多个区间分配给横轴,将测量点个数分配给纵轴,生成直方图数据。
接着,统计评价部33解析从直方图运算部32输出的直方图数据,抽取直方图数据的统计特征。而且,生成统计特征的统计评价数据,与直方图数据一起输出至图像构筑部34。图像构筑部34构筑反映了从统计评价部33输出的统计评价数据和直方图数据的图像作为直方图图像数据,输出至切换加法部9。
这里,参照图5~图8具体地说明直方图评价部20中的处理。首先,如图5的原理图所示,设在失真/弹性模数运算部13中运算出的时刻t的弹性帧数据之弹性模数的数值数据群为Ei,j(t)。这里,i=1,2,3,…,N的自然数,j=1,2,3,…、M的自然数。即,在图5中,指标i表示弹性图像横轴方向的坐标,指标j表示弹性图像纵轴方向的坐标。
在存储部31中保存弹性帧数据Ei,j(t),输出至直方图运算部32。另外,在显示于图像显示器10中的弹性图像上,检测者通过装置控制接口部17的输入部件(例如跟踪球等设备),例如,如图6所示,设定关心区域ROI。装置控制接口部17在设定的ROI坐标、例如图6所示矩形的ROI的情况下,读取4个角的坐标P1(x1,y1)、P2(x2,y1)、P3(x1,y2)、P4(x2,y2),作为ROI坐标数据而输出至直方图运算部32。
直方图运算部32使用弹性帧数据和ROI坐标数据,运算直方图数据。首先,参照ROI坐标数据,例如,如图7所示,从弹性帧数据中抽取属于ROI内部的测量点之弹性模数的数值数据群。
接着,以属于ROI内部的测量点之弹性模数的数值数据群为母集团,运算直方图数据。运算该直方图数据的处理与通常构筑直方图时的步骤相同。例如,参照图8具体地说明。若直方图的显示范围为弹性模数,设为从0(kPa)至700(kPa),并且设作为直方图的区间的管脚(pin)宽度为25(kPa),则设定28个管脚。即,各管脚的宽度如下。
管脚1:0-24.99…(kPa)
管脚2:25-49.99…(kPa)
管脚3:50-74.99…(kPa)
...
管脚27:650-674.99…(kPa)
管脚28:675-699.99…(kPa)
对应于这样设定的管脚,从母集团中计算具有各管脚范围的弹性模数的测量点个数。即,例如如下计算。
管脚1:17(个)
管脚2:87(个)
管脚3:142(个)
...
管脚27:74(个)
管脚28:63(个)
通过该处理,生成管脚序号(管脚的范围)和测量点个数的关系作为直方图数据。
统计评价部33运算生成的直方图数据的统计特征,形成统计评价数据。作为该统计特征,可从作为母集团的弹性模数数据群的标本的测量点个数、弹性模数数据的平均值、中央值、分散值、标准偏差、峰值的半值宽度、失真度、尖度等中适当选择设定。但是,表示统计特征的当然不限于此。
(步骤8)
图像构筑部34根据由直方图运算部32求出的直方图数据和由统计评价部33求出的统计评价数据,构筑直方图图像数据。构筑的直方图图像数据经切换加法部9在图像显示器10中显示,作为组织鉴别的辅助弹性信息提供给检测者。
另外,在本实施方式中,说明了使探头2接触被检体1的体表面来测量的情况,但本发明不限于此,经直肠探头、经食道探头、术中用探头、血管内探头等任意的超声波探头可同样适用。
下面,说明由直方图评价部20生成的直方图图像数据的具体实施例。另外,在下面的实施例中,以如图6所示设定ROI的情况为例来说明。
实施例1
本实施例1是生成图8所示的直方图图像的具体例。如图6所示,在ROI中包含组织1、组织2、组织3。这些组织1~3的组织固有的硬度的信息与组织1~3的范围显示一起,作为弹性图像在图像显示器10中被实时地显示。由此,可容易地把握
(组织1的硬度)<(组织2的硬度)<(组织3的硬度)的关系。
根据该ROI内的弹性模数数据生成的直方图成为图8所示的直方图图像。即,该图的直方图图像的纵轴与属于各管脚的弹性模数的测量点个数相关。因此,个数越多,管脚的长度越长。由ROI指定的弹性帧数据中的母集团保有组织1~3各自固有的弹性模数数据。因此,形成对应的弹性模数范围的测量点个数变高的直方图。
并且,在统计评价部33中,例如运算属于母集团的测量点个数、母集团的弹性模数的平均值、中央值、分散值,如图8所示,与直方图相匹配,在图像显示器10中实时显示。
并且,切换加法部9例如如图9所示,利用各种显示方式,在图像显示器10中显示图像。图示例组合了从黑白扫描转换器8输出的断层像、从彩色扫描转换器15输出的弹性图像和从直方图评价部20输出的直方图图像,构筑1个显示图像数据,在图像显示器10中实时显示。
由此,根据本实施例,可客观、定量地评价着眼区域与其周围组织以哪种硬度关系结合,可支援确定的组织鉴别。
即,在用于显示亮度(色调)大小或其扩展程度等的弹性图像的情况下,根据依赖于检测者主观的诊断,在检测者间给出不同的诊断结果,但根据本实施例,可提供避免该情况、确立客观普遍的诊断、临床上适用的超声波诊断装置。
尤其是,由于在直方图中表现出组织的集团举动的特征,所以,由此可评价病变部位等着眼部位的组织的硬度均匀性等,可识别作为非常硬的充实性组织的肿瘤、作为柔软的充实性组织的纤维腺肿、或者作为流动的囊胞性组织的囊胞等的组织性状,可提供临床上有用的诊断信息。
实施例2
在实施例1中示出使用弹性模数的实例,但本发明不限于此,作为弹性数据,也可适用失真数据、位移数据、位移矢量。即,可根据在失真/弹性模数运算部13中运算出的失真的弹性帧数据,求出同样的直方图及统计特征数据,并显示图像。并且,可根据在位移测量部12中运算出的位移帧数据,求出同样的直方图及统计特征数据,并显示图像。并且,如果由位移矢量构成位移帧数据,则可根据该矢量数据,求出同样的直方图及统计特征数据,并显示图像。
即,通过根据与组织硬度相关的弹性数据,求出直方图及统计特征数据并显示图像,可取得与实施例1相同的效果。
实施例3
在实施例1中,说明了根据由失真/弹性模数运算部13中运算出的弹性模数的数值构成的弹性帧数据、求出直方图及统计特征数据的实例。但是,本发明不限于此,即便将在弹性数据处理部14中实施了图像处理后的弹性帧数据输入直方图评价部20,求出直方图及统计特征数据,也可取得与实施例1相同的效果。
并且,即便根据在彩色扫描转换器15中转换后的弹性图像数据,求出直方图及统计特征数据,也可取得与实施例1相同的效果。
实施例4
图10示出实施例4的直方图图像的显示例。在实施例1中,说明了根据作为与弹性相关的弹性数据的弹性模数来生成并显示直方图的实例。但是,本发明不限于此,在统计评价部33中,例如可利用双高斯(ダブルガウシアン)或多高斯函数等适当的函数近似(配合(fitting))求出的直方图,从而如图10所示那样进行图像显示。这时,可求出各个高斯函数的平均值、标准偏差,输出其作为统计评价数据。即,如该图所示,可求出配合函数、该函数的峰值位置的弹性模数m1~m3、或者弹性模数的平均值M和标准偏差σ1~σ3等宽度信息,显示作为直方图的特征量数据。
根据本实施例,与实施例1相比,可进一步客观、定量地评价着眼区域与其周围组织以哪种硬度关系结合,可支援确定的组织鉴别。
实施例5
图11中示出实施例5的直方图图像的显示例。本实施例指定实施例1中示出的直方图的任意弹性模数的直方图。显示具体的直方图信息。即,经装置控制接口部17的输入部件,在直方图图像数据上显示多个光标(C1、C2、…),检测者可自由移动光标。而且,检测者指定的光标C1、C2、…位置的直方图的信息、例如个数N1、N2、平均值m1、m2等数值数据与光标的移动同步显示。
根据本实施例,可进一步客观、定量地评价着眼区域与其周围组织以哪种硬度关系结合,可支援确定的组织鉴别。
实施例6
图12示出实施例6的直方图图像的显示例。本实施例是由实施例5示出的光标C1、C2、…指定直方图的范围以及运算该范围的统计特征并图像显示的实例。
即,例如,如图12所示,若在直方图中指定光标C1~C2的范围,则运算对应于该范围的直方图数据的测量点个数N、平均值m、标准偏差σ等统计特征,显示作为数值数据。
并且,可限定于由检测者自由指定的光标C1~C2的范围,实施基于实施例4中说明的近似函数的配合,抽取其统计特征并图像显示。
实施例7
图13中示出实施例7的直方图图像的显示例。本实施例如图13所示,如以网眼示出的区域41那样,显示实施例6中选择出的直方图的光标C1~C2之范围中计算的测量点集团对应于弹性图像的哪个区域。由此,检测者可容易地识别直方图与组织的对应关系。
实施例8
图14中示出实施例8的直方图图像的显示例。本实施例是运算直方图的2个峰值位置C1、C2的弹性模数的比率(m2/m1)并显示的实例。根据本实施例,例如,若设C1是正常组织的区域,则可立即把握C2比正常组织硬几倍。
即,在以前,压迫活体组织,运算内部各区域中产生的失真或弹性模数,构筑并显示灰度化其大小的弹性图像。可是,不能把握着眼区域与其周围组织以哪种硬度关系结合。这点,根据本实施例,可客观、定量地评价着眼区域与其周围组织以哪种硬度关系结合,可进行确定的组织鉴别。
并且,同时,虽然根据弹性图像中的亮度(色调)大小或其扩展程度等依赖于检测者主观的诊断,在检测者间给出不同的诊断结果,但可提供避免该情况、确立客观普遍的诊断、临床上有用的超声波诊断装置。
实施例9
本实施例除实施例8外,还可在直方图评价部20中增加以下功能。
求出由图14的光标C1、C2等指定的峰值部的弹性模数的半值宽度,根据弹性模数的半值宽度的值,评价组织性状,输出显示其结果。即,例如,弹性模数的半值宽度为1~20(kPa)时,评价为作为非常硬的充实性组织的肿瘤。弹性模数的半值宽度为21~70(kPa)时,评价为作为柔软的充实性组织的纤维腺肿。弹性模数的半值宽度为71(kPa)以上时,评价为作为流动的囊胞性组织的囊胞。
并且,求出峰值部的分散值/平均值的比并无维化。而且,例如,在癌症肿瘤的情况下,由于分散小、平均大,所以可通过分散值/平均值的比变小(例如0.1)来评价。纤维腺肿的情况下,由于分散是中间程度、平均也是中间程度,所以通过分散值/平均值之比也变为中间程度(例如0.3)来评价。囊胞情况下,由于分散大、平均小,所以可通过分散值/平均值之比变大(例如0.5)来评价。
这样,根据本实施例,由于可输出显示评价了组织性状的结果,所以可进一步提高组织鉴别的可靠性。
实施例10
图15中示出实施例10的直方图图像的显示例。本实施例的特征在于,追随根据压力而位移的组织的移动,自动变更设定在弹性图像中的ROI,对应显示匹配变更的关心区域中的直方图。由此,由于跟踪同一组织区域来设定ROI,所以可得到适当的直方图信息。
即,通常,在使用了超声波的弹性的图像化中,用手上下移动接触于被检体诊断部位体表的超声波探头、压迫被检体,同时,运算并显示该压迫过程的弹性图像。这时,如图16所示,若固定设定的ROI的坐标数据不变,则由于对应压力而位移的组织远离ROI,或ROI外的组织进入ROI内,ROI内的母集团变化,所以即便原样不变运算ROI内的测量点的直方图,也不能客观地进行对比。因此,即便在时间经过的同时对比压迫力不同的弹性图像,也不能得到准确的鉴别信息。
因此,在本实施例中,响应被检体的压迫动作状况或跳动等活体的运动状况来更新ROI的坐标数据。即,在直方图评价部20中,作为表示由装置控制接口部17设定的ROI区域的代表点,例如,如图6所示,在ROI设定时刻t0,输入P1(t0)(x1,y1)、P2(t0)(x2,y1)、P3(t0)(x1,y2)、P4(t0)(x2,y2)。接着,在时刻t1,从图1的位移测量部12输入时刻t1的位移帧数据Di,j(t1)。对于时刻t1取得的位移帧数据Di,j(t1),可通过参照对应的坐标的位移数据,求出在时刻t0位于ROI坐标P1(t0)(x1,y1)、P2(t0)(x2,y1)、P3(t0)(x1,y2)、P4(t0)(x2,y2)的组织在时刻t1位移至哪个位置。因此,将求出的坐标P1(t1)、P2(t1)、P3(t1)、P4(t1)更新为新的ROI坐标数据。而且,使用分布在更新后的ROI坐标内的母集团,求出直方图和统计特征并显示。
这样,根据本实施例,通过压迫活体,活体内的组织在压迫过程中接近超声波探头,但通过在任意时刻重复本实施例的处理,对应被检体的压迫动作状况或跳动等活体的运动状况,更新ROI坐标数据。由此,例如,如图15所示,ROI的位置和大小因时刻不同而变化,但在任意时刻抽取活体内的同一组织区域。
这样,根据本实施例,由于根据对应于同一组织区域的弹性帧数据构成母集团,运算直方图和统计特征并显示,所以在压迫过程中其他组织不能进入ROI内。结果,虽然母集团的标本数变化,但作为用于区别组织特征的重要信息的平均值或分散值等统计特征在任意时刻都可保持稳定的值,确定地评价组织固有的值。
即,如图16所示,在固定了ROI的坐标时,由于运算出的直方图因时刻不同而不同,并且,表示根据直方图运算的统计特征的指标也因时刻不同而变动,所以妨碍客观、确定的诊断。可是,根据本实施例,由于可对应压迫过程来更新ROI坐标数据,在任意时刻抽取同一组织区域,评价其统计特征,所以可高精度地进行客观、确定的诊断。
实施例11
这里,说明直方图图像的显示例的另一实施例11。本实施例的特征在于,在直方图评价部20中,保持过去求出的直方图数据,构成按时间积累了当前时刻对每个管脚的直方图数据和过去的直方图数据的积累直方图数据。并且,可求出积累直方图数据的统计特征,并图像显示。
并且,也可在直方图评价部20中保持过去求出的直方图数据,构成按时间平滑处理了对每个管脚在当前时刻的直方图数据和过去的直方图数据的平滑直方图数据。这时,也可求出平滑直方图数据的统计特征,并图像显示。
根据本实施例,可稳定地抽取高精度的统计特征量。
实施例12
在上述的各实施例中,以实时运算直方图和统计特征且进行显示的情况为例来进行说明,但本发明不限于此,可使用存储在电影存储器部18中的图像数据来进行图2所需的处理。
上述处理经装置控制接口部17冻结(freeze)装置来进行。电影存储器部18根据装置控制接口部17的控制信号,调出过去的显示图像数据,从任意帧的显示图像数据中指定进行基于直方图的评价的显示图像数据和ROI。而且,将ROI的坐标数据输出至直方图评价部20。并且,将用于指定处理对象的帧的信息输出至失真/弹性模数运算部13,读出失真/弹性模数运算部13保存的过去的弹性帧数据,发送至直方图评价部20,进行直方图和统计特征的运算。由此,可离线实施本发明。
实施例13
上述实施例中说明的直方图的显示与否、或者直方图的显示设定(管脚宽度、纵轴横轴的范围等)可任意地设定、变更。并且,可利用装置控制接口部17自由指定诸如运算并显示统计特征的平均值、标准偏差等统计信息。
实施例14
图17示出本发明使用作为与弹性相关的弹性数据之一的位移矢量来生成直方图的实施例14的显示例。本实施例根据对被检体的组织施加压力测量出的超声波断层数据,求出被检体断层部位的多个测量点中的组织的位移矢量,对应于断层部位的测量点,生成用于表示位移矢量的弹性图像并显示,用直方图显示设定在弹性图像中的关心区域的特定方向即深度方向之位移矢量分量的分布。另外,该图中的直方图的横轴表示位移的深度方向分量。
图17(A)示出设定在ROI内的着眼区域是作为非常硬的充实性组织的癌症肿瘤的情况。这时各测量点的位移矢量如图17(A)的上图用箭头所示那样,面向同一方向。因此,表示ROI内的位移矢量分布的直方图如图17(A)的下图所示,成为集中于狭窄范围的直方图。即,该实例由于组织要素间牢固结合,所以可理解即便受到压迫,整体也不变形地运动。
图17(B)示出设定在ROI内的着眼区域是作为柔软的充实性组织的纤维腺肿的情况。这时各测量点的位移矢量如图17(B)的上图用箭头所示,方向稍微分散。因此,表示ROI内的位移矢量分布的直方图如图17(B)的下图所示,分布的宽度稍微变宽。即,该实例在组织要素间有某种程度结合,可理解在受到压迫时,边相互束缚、边某种程度变形。
图17(C)示出设定在ROI内的着眼区域是作为流动的囊胞性组织的囊胞的情况。这时各测量点的位移矢量如图17(C)的上图用箭头所示,位移的方向分散得大。因此,表示ROI内的位移矢量分布的直方图,如图17(C)的下图所示,分布的宽度扩展到大的范围。即,本实例由于组织要素间未完全结合,所以可理解一旦受到压迫,则各要素群分别自由运动。
由此,如果由位移矢量的深度方向分量或横方向分量构成直方图,则可识别组织性状的差异。尤其是,不仅是硬度的差异,还通过观察分布宽度,可识别组织要素间的结合程度,从而得到组织鉴别有效的辅助弹性信息。
实施例15
图18示出本发明实施例15的弹性信息的显示例。本实施例是通过对应于关心区域ROI内的直方图的弹性数据的宽度来变更弹性图像的色调或亮度的灰度(例如256灰度)分配范围,使关心区域ROI内的弹性差异的对比度显眼的实施例。即,图18(a)显示了视野内的弹性图像和直方图,图18(b)手动或通过交界检测自动在弹性图像中设定关心区域ROI,显示了关心区域ROI内的测量点的直方图。相反,如图18(c)所示,对应于关心区域ROI内的测量点之直方图的弹性数据的宽度,自动变更灰度的分配范围,使对比度显眼,来显示ROI内的弹性差异。
根据本实施例,由于强烈显示ROI内的组织之硬度分布的对比度,所以得到组织鉴别有效的辅助弹性信息。
实施例16
图19中示出本发明实施例16的弹性信息的显示例。在图9示出的实施例1的弹性信息的显示例中,示出在1个画面中并排显示弹性图像和包含直方图的辅助弹性信息的图像(直方图图像)的实例。在本实施例中,如图19所示,通过小窗口显示,在弹性图像的一部分上重叠并显示直方图图像。
即,由于即便直方图图像牺牲分辨率也不会发生看漏,所以不必以与断层像或弹性图像同等的大小对等地显示。因此,如本实施例,可在弹性图像的图像区域内用小窗口显示缩小后的直方图图像。据此,由于可用足够的分辨率观察弹性图像的弹性信息,所以可减少信息的损失。
另外,根据图19的实例,担心因会隐藏重叠显示直方图图像的部分的弹性图像信息、根据情况不同会丢失重要的弹性信息。该情况如图20所示,优选地,由半透明图像构成直方图图像,重叠在弹性图像上进行显示。据此,可避免重叠显示直方图图像导致的弹性信息的看漏。
这里,说明生成彩色半透明图像的、称为阿尔法混合的通常的图像构筑方法。例如,若用三原色(R,G,B)表现图像的彩色信息,则彩色弹性图像的坐标(i,j)的像素的彩色信息、E(i,j)如下式,为3原色的彩色信息Er(i,j)、Eg(i,j)、Eb(i,j)的函数。
E(i,j)=(Er(i,j),Eg(i,j),Eb(i,j))
同样地,彩色直方图图像的坐标(i,j)的像素的彩色信息为H(i,j),为下式的函数。
H(i,j)=(Hr(i,j),Hg(i,j),Hb(i,j))
这里,若设两图像间的重叠比率为α(0≤α≤1),重叠了直方图图像的区域之半透明重叠图像的坐标(i,j)之像素的彩色信息G(i,j)可用下式求出。
C(i,j)=α×E(i,j)+(1-α)×H(i,j)
即,半透明重叠图像的坐标(i,j)之像素的彩色信息按α比率具有弹性图像和直方图图像两者的图像信息,可半透明地同时识别各图像的彩色信息。
这样,根据本实施例,若在限定的监视器显示画面内以同等大小不干扰地配置多个图像,则即便各图像缩小的情况,也可避免降低图像诊断的精度。
实施例17
图21示出本发明实施例17的弹性信息的显示例。图19及图20中,示出了在断层像或弹性图像等诊断图像区域内利用窗口显示缩小后的直方图图像的实例。在本实施例中,如图21所示,根据需要,切换诊断图像和直方图图像的显示区域的大小关系。
例如,可构成为通过操作安装在超声波诊断装置的装置控制接口部17中的规定键或设定在按钮中的图像切换键,来切换诊断图像和直方图图像的显示区域的大小关系。另外,图像切换键不限于装置控制接口部17的键或按钮,可在监视器的图像上显示,并设定成图标。
根据本实施例,在图11、图12的实施例中,可在大的直方图图像上高效地执行直方图图像上的光标移动等操作。
实施例18
图22示出将本发明适用于鉴别甲状腺滤胞癌的诊断时的弹性信息的显示例。在甲状腺的图像诊断中,甲状腺滤胞癌的鉴别非常重要。已知滤胞癌在肿瘤内部的细胞密度有分布,肿瘤边缘部具有高于肿瘤中心部的细胞密度。另外,在利用弹性图像评价滤胞癌时,存在将肿瘤边缘图像显示为比肿瘤中心部硬的区域的倾向。即,认为将细胞密度高的部位作为比细胞密度低的部位硬的部位来测量。
本实施例显示适于评价肿瘤内部的弹性分布特征的情况的弹性信息。例如,即便如甲状腺肿瘤那样、在肿瘤边缘部和中心部之间存在硬度差异的情况,在通常的弹性图像中多数情况不能识别该差异。图22是甲状腺肿瘤的弹性图像和滤胞癌之情况下的直方图图像的显示例。如该图所示,为了客观地鉴别滤胞癌,通过将肿瘤中心部ROI1和环状的肿瘤边缘ROI2作为各自独立的母集团来直方图解析,取得对鉴别有效的信息并显示。
即,在本实施例中,以下面的显示为特征。
(1)在认为滤胞癌的区域的边缘部和中心部各自独立地设定ROI1和ROI2并显示。
(2)对每个ROI分色ROI1和ROI2的直方图,并合成显示。
(3)运算各直方图的统计信息并显示。在图示例中,显示平均、标准偏差、中央值。
考虑各种ROI1、ROI2的设定方法。例如,操作装置控制接口17的鼠标或跟踪球等输入部件,由检测者用线描绘环状的边缘部ROI1的外周和内周并指定。由此,该描绘线间的环状区域可识别为边缘部ROI1,并进行设定。
并且,预先准备椭圆形的多个模板,操作鼠标或跟踪球,指定成使2个椭圆形模板与边缘部ROI1的外周和内周吻合。由此,可设被2个椭圆形模板的外周线夹持的环状区域为边缘部ROI1,将小的椭圆形模板的内侧区域识别为RROI2,并进行设定。
并且,预先准备适于边缘区域指定用的多个圆环(torus)(轮环体)状的ROI,利用鼠标或跟踪球等操作,在图像上配置2个圆环ROI。由此,可设被2个圆环ROI夹持的环状区域为边缘部ROI1,将小的圆环ROI的内侧区域识别为RROI2,并进行设定。
根据本实施例,即便是不能定量地识别肿瘤的边缘部和中心部的硬度差异的弹性图像,也可进行定量客观的滤胞癌诊断。
实施例19
例如,如文献「Ueno E,et al:Dynamic tests in real-time breastechography.Ultrasound Med Biol,14:53-57,1988」所述,在乳腺区域的图像诊断中,已知评价肿瘤的可动性(称为活动性(mobility))的动态测试。所谓动态测试,是评价关心肿瘤相对于其周围组织被限制到何种程度的测试。具体地说,是在使整个乳房大变形时,边实时观察断层像,边评价是表示着眼对象的肿瘤被周围组织的运动拖拉移动,还是未被周围组织的运动拖拉的独立的运动的诊断方法。例如,由于即便乳癌中,浸润性的肿瘤也浸润至周围组织中发育,所以与周围组织的限制度高,可动性低。另外,由于良性的肿瘤不示出浸润性,所以不示出与周围组织结合的举动,认为可动性高。但是,由于该可动性评价依赖于检测者的经验、主观,所以客观性低。
因此,在本实施例中,使用直方图进行可动性评价。图23(a)、(b)中示出本实施例的弹性信息的显示例。图23(a)示出的关心部假定良性肿瘤等中可动性高的肿瘤的情况。这时,由于肿瘤不限制在周围的乳腺部,所以在从外部使整个乳房大变形时,在肿瘤和乳腺的交界部出现滑动(slip)的举动。即,在交界的肿瘤侧的区域和交界的乳腺侧的区域中,因在位移方向和大小上出现差异,看情况沿反方向位移。因此,如图23(a)所示,在包含肿瘤和乳腺的交界的区域设定ROI,在该ROI内设定例如提供0.05mm的深度方向的位移时的压迫操作作为评价的基准。若将这时ROI内的横方向位移大小显示为直方图,则如图所示,在正方向(例如,图像中向右)和反方向(例如,图像中向左)出现2个具有相同程度大小位移的直方图的峰值。该直方图的扩展程度(峰值的半值宽度、标准偏差或分散等)或峰值间的距离(mm)可评价为提供可动性的程度的尺度。
另外,图23(b)中示出的关心部是恶性肿瘤的可动性低的肿瘤的情况。这时,由于肿瘤被强有力地限制在周围的乳腺部,所以在从外部使整个乳房大变形时,肿瘤和乳腺的交界部进行同等的运动。即,在交界的肿瘤侧的区域和交界的乳腺侧的区域中,位移的方向和大小的差异减小,看情况而完全结合和位移。因此,与图23(a)相同,在包含肿瘤和乳腺的交界的区域设定ROI,以该ROI内提供了0.05mm的深度方向的位移时的压迫操作为评价基准,若表示为直方图,则ROI内的横方向的位移的大小变为该图(b)那样。即,可在正方向(例如,图像中向右)或反方向(例如,图像中向左)的任一方出现具有相同程度大小的位移的直方图的峰值,将直方图的扩展程度(峰值的半值宽度、标准偏差或分散等)评价为提供可动性的程度的尺度。
实施例20
在实施例19中,说明了作为好坏性鉴别的1个信息,为了客观地评价与肿瘤对周围的浸润度有关的可动性,提供基于位移信息的直方图作为辅助弹性信息的实例。在本实施例中,作为良恶性鉴别的1个信息,示出作为评价组织间的粘连(syncretio)用的辅助弹性信息的适用直方图的方法。
图24示出在乳腺组织中具有肿瘤时的弹性图像的实例。尤其是,是适于进行摘出在乳腺和大胸筋交界部附近的肿瘤的手术时的辅助弹性信息的实例。在摘出在乳腺和大胸筋交界部附近的肿瘤时,从乳腺与大胸筋之间接近(approach)肿瘤。在进行这种手术时,肿瘤和大胸筋、及乳腺和大胸筋是否不粘连、可剥离非常重要。
因此,在本实施例中,如图24(a)所示,在希望评价乳腺和大胸筋的粘连的区域设定ROI。与图23的情况相同,如果乳腺与大胸筋不粘连,则乳腺和大胸筋产生相对乳房变形各自进行独立的运动、在交界部进行滑动的现象。因此,例如若提供规定大小的深度方向的位移,显示横方向的位移作为直方图图像,则如图24(a)所示,出现2个直方图的峰值。
另外,乳腺和大胸筋强有力地粘连时,如图24(b)所示,由于乳腺和大胸筋结合和运动,所以得到为在一定位移处具有峰值的分布。因此,成为利用ROI中的直方图的扩展程度(峰值的半值宽度、标准偏差或分散等)来评价粘连或剥离的可能性的尺度。
这样,可在是否可剥离的评价中使用基于位移信息的直方图,可利用作为制定肿瘤摘出手术中的手术计划用的重要信息。
另外,粘连的评价不限于摘出手术的情况,在丰胸手术等中在乳腺和大胸筋之间插入袋(生理盐水袋或硅袋等)的情况下,也可有效地用于评价组织间的粘连。在使用内视镜的外科手术中,其效果也相同。
实施例21
参照图25,说明将本发明的、定量的辅助弹性信息适用于高密度焦点式超声波(下面称为HIFU:high-intensity focused ultrasound)的治疗效果判定的本实施例21。
基于HIFU的治疗是从超声波治疗用探针(probe)向被检体内的治疗部位放射超声波,使高密度的超声波聚焦于治疗部位的小区域,将焦点部位加热至例如60~90℃的高温,使病变部(例如癌等的病灶)的组织死亡。而且,使焦点在治疗部位内稍微移动,使病变部的全部组织死亡来治疗。
由于通过该治疗死亡的病变部的组织因热变性而变硬,所以可使其硬度程度和治疗效果相对应。可是,如果在治疗部位设定关心区域ROI、测量ROI内的弹性信息并观察弹性图像,则可判定是否在治疗部位的组织全部区域均匀地实施治疗。可是,根据以对应于硬度的灰度所显示的弹性图像,如背景技术栏中所述,就色灰度近似的区域的微小部位而言,存在不能地明确地识别灰度差、不能明确地区别是否实施治疗的情况。这时,可利用作为本发明的辅助弹性信息的直方图,定量及客观地支援正确的诊断。
图25(a)示出用弹性图像描绘包含治疗前的前列腺癌的前列腺的截面情况。这时,若在前列腺癌的区域中设定ROI,如图所示,可知例如具有以100(kPa)程度为中心的弹性。
下面,在图25(b)中示出对前列腺癌区域实施HIFU治疗的中途过程的情况。充分照射了超声波的前列腺癌区域因热变性而硬化,例如图25(b)所示,利用直方图识别变化成300(kPa)程度的硬度。可是,如从图25(b)的直方图把握的那样,可知具有表示未治疗的前列腺癌的硬度的100(kPa)硬度的组织区域残留在ROI内部,该区域是未充分照射超声波的治疗不完全的区域。
这时,与图13的实施例相同,如果用光标C1、C2选择治疗不完全的区域,执行将该测量点反馈至弹性图像中显示的处理,则可在图像上示出治疗不完全的区域,检测者可识别。
而且,最后,如图25(C)所示,通过对热变性不充分的治疗不完全的区域再次实施HIFU治疗、确认直方图,可确认在ROI内的全部区域中均匀地进行治疗。
另外,作为治疗结束的判断基准,也可设定200(kPa)以上作为阈值。由此,可对治疗效果判定设定共同条件,可实现根据一定治疗基准的可靠治疗。
并且,热烧灼癌组织治疗的方法不限于HIFU,RFA(经皮的放射波烧灼疗法)的情况也可适用本发明。并且,由于即便相反急剧冷却癌组织治疗的低温疗法,也仅对应热变性导致的硬化,所以可同样适用本发明。
并且,还包含与治疗方法无关、例如利用药剂的激素疗法等,在癌组织等的治疗中,对于硬度随着组织变性而变化的治疗法,作为评价是否均匀地进行该治疗处理、或是否均匀地得到治疗效果的方法,可利用本发明。
实施例22
参照图26说明最适于乳癌检诊的本发明的实施例22。在乳癌的检诊中,通常,如图26(a)、(b)所示,将超声波探头触及乳房扫描,实时扫描乳房的断层像,进行乳癌检诊。在这种检诊中,难以对一个被验者分配足够的检查时间,要求高效的诊断。可是,无遗漏地观察被实时扫描的断层像之测量截面内的全部区域是非常需要时间的作业,也大大依赖于经验。即,B模式像等断层像是将活体内的组织结构构筑为图像的信息,将或形成肿瘤像、或后方的回声增强、或边缘部粗杂等组织的形态信息利用于鉴别。即,由于根据非常多的复杂角度检测或鉴别肿瘤,所以为了进行高效的诊断,要求检测者的经验。
可是,由于肿瘤比周围的脂肪组织或乳腺组织硬,所以例如在进行基于弹性模数作为弹性信息的检查时,观察在测量截面内是否存在具有大的弹性模数的区域。并且,由于已知越恶性的组织越硬,所以仅对测量截面中是否包含具有异常硬度的区域集中注意即可,所以即便实时也可高效地识别。
可是,如背景技术栏中所述,在要根据弹性图像判别弹性的不同、鉴别乳癌等肿瘤时,就不同色灰度的交界区域或色灰度近似的区域而言,有时不能明确区别弹性的差异。因此,必需根据弹性图像普遍且更注意观察测量截面的全部区域,实时地扫描测量截面,执行该动作使检测者的负担大。
因此,在本实施例中,通过除弹性图像之外,或者还直方图化弹性信息的分布、显示作为辅助弹性信息来替代弹性图像,实现高效的检诊。即,如图26(c)所示,在截面的扫描中实时测量失真或弹性模数等弹性信息,直方图化该弹性信息的分布后显示。例如,在设怀疑恶性的弹性模数的阈值Eth为200(kPa)时,通过仅观察直方图分布,可立即判断存在于当前弹性图像中的Eth以上弹性模数的测量点分布。即,如果是弹性图像,则可确认至角落,但根据本实施例,没有这种负担。
并且,在对直方图的纵轴(度数)也设定阈值Cth,测量Eth以上的测量点Cth以上时,如果在直方图的显示图像上显示「有弹性异常区域」等警告,则可防止看漏。
并且,直方图的纵轴(度数)标度也可从线性标度变更成记录标度。通过成为记录标度,即便在怀疑恶性的硬的肿瘤非常小、度数少时,也可无遗漏地高灵敏度地检测该肿瘤。
以上,如各实施例中所说明的那样,通过在图像显示器10的显示画面中,除根据位移、失真或弹性模数等的弹性图像之外,或者代替弹性图像,来显示包含弹性分布的直方图的辅助弹性信息,可高效地实现检测者负担少的诊断。