超声波诊断装置以及方法.pdf

上传人:v**** 文档编号:5234512 上传时间:2018-12-28 格式:PDF 页数:30 大小:1.64MB
返回 下载 相关 举报
摘要
申请专利号:

CN201410637106.2

申请日:

2014.11.06

公开号:

CN104622502A

公开日:

2015.05.20

当前法律状态:

授权

有效性:

有权

法律详情:

授权|||专利申请权的转移IPC(主分类):A61B 8/00登记生效日:20160722变更事项:申请人变更前权利人:日立阿洛卡医疗株式会社变更后权利人:株式会社日立制作所变更事项:地址变更前权利人:日本国东京都变更后权利人:日本东京都|||实质审查的生效IPC(主分类):A61B 8/00申请日:20141106|||公开

IPC分类号:

A61B8/00

主分类号:

A61B8/00

申请人:

日立阿洛卡医疗株式会社

发明人:

田原麻梨江; 吉川秀树; 田中宏树

地址:

日本国东京都

优先权:

2013-232539 2013.11.08 JP

专利代理机构:

中科专利商标代理有限责任公司11021

代理人:

李国华

PDF下载: PDF下载
内容摘要

本发明涉及超声波诊断装置以及方法。在生物体等对象组织为粘弹性体的情况下,剪切波的传播速度随着剪切波的频率而发生变化。超声波诊断装置通过来自位移生成部(10)的信号从超声波探头(1)照射位移生成用送波波束而产生剪切波之后,通过位移检测部(30)对生物体组织照射脉冲波,检测剪切波的粒子速度,根据由粘弹性解析部(34)检测得到的剪切波的粒子速度的波形的时间广度来估计粘性参数,并在显示部(7)中显示估计出的粘弹性。

权利要求书

权利要求书
1.  一种超声波诊断装置,其特征在于,具有:
位移生成部,其对被检测体内施加声音放射压来使被检测体内产生位 移;
位移检测部,其通过在上述被检测体的多个位置上发送接收超声波来 检测在上述被检测体内产生的剪切波的粒子速度;和
粘弹性解析部,其根据上述粒子速度来估计粘性参数,
上述粘弹性解析部根据上述粒子速度波形的时间广度来估计上述粘 性参数。

2.  根据权利要求1所述的超声波诊断装置,其特征在于,
上述粘弹性解析部根据上述粒子速度成为正和负的两个时刻来算出 上述粒子速度波形的时间广度。

3.  根据权利要求1所述的超声波诊断装置,其特征在于,
上述粘弹性解析部根据上述粒子速度成为过零值的时刻、或者成为上 述过零值的时刻前后的上述粒子速度的积分值,来算出上述粒子速度波形 的时间广度。

4.  根据权利要求1所述的超声波诊断装置,其特征在于,
上述粘弹性解析部基于在上述多个位置上算出的上述粒子速度波形 的时间广度,通过拟合来估计上述粘性参数。

5.  根据权利要求4所述的超声波诊断装置,其特征在于,
上述粘弹性解析部采用指数曲线来进行上述拟合。

6.  根据权利要求1所述的超声波诊断装置,其特征在于,
上述粘弹性解析部算出上述粒子速度的中心频率。

7.  根据权利要求1所述的超声波诊断装置,其特征在于,
上述粒子速度通过一个参考信号与由上述位移检测部在不同的时刻 所接收到的多个信号的相关运算而被计算出。

8.  根据权利要求1所述的超声波诊断装置,其特征在于,
上述粒子速度通过由上述位移检测部在不同的时刻所接收到的多个 信号彼此的相关运算而被计算出。

9.  根据权利要求1所述的超声波诊断装置,其特征在于,
还具有反馈参数决定部,
上述反馈参数决定部基于通过第一测量得到的、根据上述粒子速度波 形的时间广度所估计的上述粘性参数,对接下来的测量即第二测量中的发 送条件、以及测量区域当中的至少一者进行调整。

10.  根据权利要求9所述的超声波诊断装置,其特征在于,
上述反馈参数决定部使进行调整的上述测量区域的大小对应于进行 测量的位置而发生变化。

11.  根据权利要求10所述的超声波诊断装置,其特征在于,
上述测量区域的形状为梯形。

12.  根据权利要求9所述的超声波诊断装置,其特征在于,
具有上述反馈控制的接通/断开的切换开关。

13.  一种超声波诊断装置中的粘性参数估计方法,其特征在于,
对被检测体内施加声音放射压来使被检测体内产生位移,
通过在上述被检测体的多个位置上发送接收超声波来检测在上述被 检测体内产生的剪切波的粒子速度,
根据检测出的上述粒子速度波形的时间广度来估计粘性参数。

14.  根据权利要求13所述的粘性参数估计方法,其特征在于,
根据上述粒子速度成为正和负的两个时刻、或者上述粒子速度成为过 零值的时刻、或者成为上述过零值的时刻前后的上述粒子速度的积分值, 来算出上述粒子速度波形的时间广度。

15.  根据权利要求13所述的粘性参数估计方法,其特征在于,
基于通过第一测量得到的、根据上述粒子速度波形的时间广度所估计 的上述粘性参数,对接下来的测量即第二测量中的发送条件、以及测量区 域当中的至少一者进行调整。

说明书

说明书超声波诊断装置以及方法
技术领域
本发明涉及通过超声波收发来进行与被检测体的粘弹性相关的测量 的技术。
背景技术
作为乳腺癌、肝硬化、血管病等的诊断方法,存在根据超声波回波信 号诊断被检测体内部的硬度的方法(弹性成像技术),来代替医生的触诊。 在基于弹性成像技术的硬度的诊断中,操作者将超声波探针压入被检测体 即生物体表面并进行加压,使被检测体内部的组织产生位移。根据加压所 引起的生物体组织的压缩前后的回波信号来估计压缩方向的位移,求得位 移的空间微分量即变形。进而,根据变形和应力算出与硬度相关的值,例 如杨氏模量。在该方法中,具有如下课题:作为拍摄对象,局限于容易从 体表加压之处所存在的内脏器官。例如,作为在体表与肝脏之间的间介层, 存在滑动面,因此使足够的位移产生那样的加压是困难的。
因而,存在如下技术:采用位移生成用超声波会聚波束、即,推力波 束(push beam),或者推力脉冲来对被检测体内部施加声音放射压,抑制 间介层的影响并使对象组织位移,由此进行硬度的诊断。在这些技术中, 对在会聚波束的前进方向上产生的组织的位移量进行图像化,或根据伴随 着焦点处的组织位移而在与会聚波束的前进方向垂直的方向上产生的剪 切波的传播速度的估计来算出刚性率或杨氏模量之类的弹性率。如果采用 该技术,则期待减小上述的滑动面等的间介层的影响的效果等。
进而,近年来,作为慢性肝疾患或肿瘤的诊断信息在临床现场利用弹 性成像技术不断扩大中,想要实现包括粘性的粘弹性测量的技术正在不断 被开发。本来、由于体组织为具有粘性和弹性这两种特性的粘弹性体,因 此更加期待基于实际情况的高精度的粘弹性评估技术。
在粘弹性的评估中,除了根据剪切波的群速度算出刚性率和杨氏模量 之外,还需要粘性系数等粘性参数的估计。在粘弹性的评估中,如美国公 开US2007/0038095号公报和美国公开US2011/0063950号公报中所公开的 那样,具有采用粒子速度(Particle Velocity)波形的手法。
发明内容
在粘弹性的评估中,在对象组织为粘弹性体的情况下,剪切波的传播 速度随着剪切波的频率而产生变化,但采用上述粒子速度波形的粘弹性评 估具有难以受到体动等的剪切波以外的频率分量的影响之类的优点。
另一方面,在美国公开US2007/0038095号公报中,将经AM调制的 推力波束照射到对象组织。由于使具有100~1kHz程度的频率的剪切波产 生,因此在最大15ms期间照射推力波束。此外,在美国公开 US2011/0063950号公报中,每一次的推力波束的照射时间为1ms,但需要 每隔10ms照射多次推力波束(tone burst)。如上那样,在现有的方法中, 由于总计的推力波束的照射时间长,因此存在在生物体内部的局部的温度 上升变大之虞。为了抑制这种局部的温度上升,需要例如采用一次1ms以 内的更少的照射次数/照射时间下的推力波束的照射来评估粘性。
粘性对于剪切波的影响在于包括频率依赖性在内的相位速度或衰减 率的不同。由于该影响伴随着剪切波的传播而放大,因此基本上传播距离 越长越适于粘弹性测量。即、粘弹性测量的优选条件之一为使高振幅的剪 切波产生,这意味着与上述的温度上升所产生的限制条件同时成立变得困 难。
本发明的目的在于,提供一种解决上述的课题,将生物体组织的温度 上升抑制到最小限度,能高精度地解析粘性参数的超声波诊断装置以及粘 性参数估计方法。
为了实现上述的目的,在本发明中提供以下的超声波诊断装置,即具 有:位移生成部,其对被检测体内施加声音放射压来使被检测体内产生位 移;位移检测部,其通过在被检测体的多个位置上发送接收超声波来检测 在被检测体内产生的剪切波的粒子速度;以及粘弹性解析部,其通过粒子 速度来估计粘性参数,粘弹性解析部根据粒子速度波形的时间广度来估计 粘性参数。
此外,为了实现上述的目的,在本发明中,提供以下的超声波诊断装 置中的粘性参数估计方法,即对被检测体内施加声音放射压来使被检测体 内产生位移,通过在被检测体的多个位置上发送接收超声波来检测在被检 测体内产生的剪切波的粒子速度,并根据检测出的粒子速度波形的时间广 度来估计粘性参数。
发明效果
根据本发明,能将生物体组织的温度上升抑制到最小限度,并且能高 精度地评估粘性参数。
附图说明
图1为表示实施例1相关的超声波诊断装置的一系统结构的图。
图2为说明实施例1相关的超声波探头所产生的位移生成的图。
图3为说明实施例1相关的超声波探头所产生的超声波的波束形成的 图。
图4为说明实施例1相关的粘性参数的计算时序的图。
图5A为表示实施例1相关的剪切波的粒子速度检测处理的一例的图。
图5B为表示实施例1相关的粒子速度检测处理的另一例的图。
图6A为表示实施例1相关的剪切波的位移的时间波形的图。
图6B为表示实施例1相关的剪切波的粒子速度的时间波形的图。
图7A为表示实施例1相关的粒子速度的时间波形的一例的图。
图7B为表示实施例1相关的粒子速度的时间波形的另一例的图。
图8为说明粘弹性模型的图。
图9A为表示剪切波的传播速度的频率依赖性的图。
图9B为表示没有粘性的情况下的剪切波的粒子速度波形的一例的 图。
图9C为表示有粘性的情况下的剪切波的粒子速度波形的一例的图。
图10A为表示实施例1相关的具有粘性的情况下的粒子速度的时间波 形的图。
图10B为说明实施例1相关的具有粘性的情况下的粒子速度的距离衰 减特性的图。
图11A为说明实施例1相关的根据粒子速度的时间波形的正的峰值和 负的峰值来估计粘性参数的处理的图。
图11B为表示实施例1相关的用于根据粒子速度的时间波形来估计粘 性参数的处理的模块图。
图12为说明实施例1相关的粒子速度波形和距离衰减特性的图。
图13A为说明实施例1相关的剪切波的距离衰减特性的图。
图13B为说明实施例1相关的用于根据剪切波的距离衰减特性来估计 粘性参数的处理的图。
图14A为表示实施例1相关的粘性参数的表示例的图。
图14B为表示实施例1相关的粘性参数的表示例的图。
图15为表示实施例2相关的超声波诊断装置的一系统结构的图。
图16为表示实施例2中的粘性参数的计算时序的图。
图17A为说明实施例2中的ROI重新设定处理的图。
图17B为表示实施例2中的ROI重新设定处理的标示例的图。
图18为表示实施例2中的测量区域的决定过程的图。
标号说明
1   超声波探头
2   收发切换开关
3   中央控制部
4   彩色DSC
5   黑白DSC
6   合成部
7   显示部
10  位移生成部
12  送波波束条件设定部
13  位移生成用送波波束生成部
20  第一超声波收发部
30  位移检测部
31  第二超声波收发部
32  位移运算部
33  测定范围设定部
34  粘弹性解析部
38  反馈参数决定部
50  彩色标度设定部
51、52  相关运算
80  Voigt模型
100  超声波探头1的各元件
101  正的峰值
102  负的峰值
110  t(xn)
111、131  FFT
112、132  粘性参数α的评估
113、133  拟合
114、134  η输出
130  p(xn)或者p(xn)’
141  标绘显示
142  2D彩色地图
143  彩色标度
171  测量区域ROI的重新设定
172  测量区域
181  粘性参数
182  ROI控制
具体实施方式
以下,参照附图对本发明的实施方式进行说明。在各附图中,相同的 标号表示同一构成要素。
实施例1
第1实施例为以下的超声波诊断装置以及粘性参数估计方法的实施 例,即,该超声波诊断装置具有对被检测体内施加声音放射压来在被检测 体内位移的位移生成部10;通过对被检测体的至少两个以上的多个位置发 送接收超声波来检测在被检测体内生成的剪切波的粒子速度的位移检测 部30;以及根据粒子速度来估计粘性参数的粘弹性解析部34,粘弹性解 析部根据粒子速度波形的时间广度来估计粘性参数。
图1中表示实施例1的装置整体结构的一例。用于与未图示的被检测 体的表皮接触的超声波探头1形成为具有对在与被检测体之间发送以及接 收超声波的多个振荡器进行排列的超声波收发面。作为装置的构成要素的 各模块、即包括用于在被检测体内产生位移的位移生成部10、收发切换开 关2、第一超声波收发部20、第二超声波收发部31、检测位移的位移检测 部30、粘弹性解析部34、彩色标度设定部50被中央控制部3控制。
超声波探头1经由收发切换开关2与位移生成用送波波束生成部13、 第一超声波收发部20、以及第二超声波收发部31连接。发送接收切换开 关2被中央控制部3控制为切断或者放开超声波探头1与位移生成用送波 波束生成部13、第一超声波收发部20、以及第二超声波收发部31的连接。
第一超声波收发部20被中央控制部3控制,以使采用由省略了图示 的位移检测用送波波形生成部所作成的波形,对超声波探头1的元件的送 波信号给予延迟时间和权重,在被检测体所期望的位置即焦点集中位移检 测用超声波波束、即,推力波束、或者推力脉冲。在被检测体内反射并返 回到探头的回波信号在超声波探头1中被变换为电信号,并发送到第一超 声波收发部20。第一超声波收发部20包括将回波信号整相相加并进行包 线检波和log压缩和带通滤波、增益控制等的信号处理电路。从第一超声 波收发部20发出的输出信号被输入到黑白数字扫描变换器(Digital Scan  Converter:DSC)5。在黑白DSC5中能够形成由黑白构成的表示亮度的断 层像、即B模式像的信息。
对位移生成部10进行说明。在位移生成部10的送波波束条件设定部 12中,设定会聚位置、振幅值、用于送波的超声波探头1的各元件的元件 数、载波的频率等。位移生成用送波波形生成部11中,采用由送波波束 条件设定部12所设定的载波的频率来生成波形。位移生成用送波波束生 成部13被中央控制部3控制,以使采用由位移生成用送波波形生成部11 所作成的波形以及由送波波束条件设定部12所设定的设定条件,对超声 波探头1的各元件的送波信号给予延迟时间和权重,在被检测体内部会聚 超声波波束。来自位移生成用送波波束生成部13的电信号在超声波探头1 中被变换为超声波信号,朝向未图示的被检测体,照射位移生成用超声波 波束。
对位移检测部30进行说明。位移检测部30内的测定范围设定部33 中,设定对位移进行测定的范围。第二超声波收发部31被中央控制部3 控制,以使采用由省略了图示的位移检测用送波波形生成部作成的波形、 以及由测定范围设定部33所设定的测定范围,对超声波探头1的各元件 的送波信号给予延迟时间和权重,在未图示的被检测体的期望的位置、焦 点会聚位移检测用超声波波束。在被检测体内反射并返回到超声波探头的 回波信号在超声波探头1中被变换为电信号,并发送到第二超声波收发部 31。
位移检测部30的第二超声波收发部30包括对回波信号进行整相相 加、进行包线检波和log压缩和带通滤波、增益控制等的信号处理电路。 来自第二超声波收发部31的输出信号被输入到位移运算部32。此外,位 移运算部32中,通过后面详细叙述的相关运算来运算各部位的位移、或 者粒子速度。此外,从位移运算部32输出的位移或者粒子速度的信息被 输入到粘弹性解析部34。
在粘弹性解析部34中,基于所输入的信息,对剪切波的传播速度、 进而粘性参数进行解析、估计。由粘弹性解析部34解析/估计的粘性参数 的值被发送到彩色DSC4。在彩色DSC4中采用该粘性参数的值,形成表 示粘弹性的断层像的信息作为彩色信息。来自彩色DSC4的输出在合成部 6中被合成,并被显示于显示部7中与显示B模式像相同的画面上。
进而,采用粘弹性解析部34所计算的粘性参数等的粘弹性信息经由 中央控制部3被发送到彩色标度设定部50中。在彩色标度设定部50中, 基于粘性参数等粘弹性的值,作成与粘弹性彩色图像相对应的彩色标度。 如后面所说明的那样,该彩色标度优选相邻于与B模式像重叠等显示的表 示粘弹性的大小的彩色图像并显示于显示部7中。
另外,在该图中由模块表示的构成要素的一部分即中央控制部3、粘 弹性解析部34等能够通过中央处理部(Central Processing Unit:CPU)所进 行的程序执行来实现。作为该CPU,能够采用具备通常的CPU和存储部 (存储器)和输入输出部的计算机。这种情况下,显示部7能利用计算机 的输入输出部的显示器等。
在本实施例中,对使图2中所示的由多个元件100构成的线性阵列型 的超声波探头1与被检测体的体表面接触,在体内的目的断层面会聚位移 生成用超声波波束的情况进行说明。在此,对在期望的断层面内,位移生 成用超声波波束的传播方向处于相对体表垂直的方向的情况进行说明。超 声波的波束形成如图3所示那样,通过以下来实现,即求得各焦点与超声 波探头1的各元件100的位置之间的距离,对每个元件给予将通过元件间 的距离差除以对象物的音速所算出的延迟时间并进行送波。该图的下层示 意性地表示与上层相比作为焦点的集点在上侧被扫描的情况。
如图2所示那样,如果在焦点(图2中的集点F)上照射会聚波束, 则按照伴随着传播的超声波的吸收或散射而产生辐射压。通常情况下焦点 处辐射压成为最大,在焦点区域的生物体组织中产生位移。此外,如果会 聚波束的照射被停止,则位移量被缓和。通过该辐射压的生成,如图2所 示那样以集点F为起点在与被检测体表面相平行的方向上产生剪切波。位 移的方向为相对体表垂直的方向。
接下来,采用图4对实施例1中的粘性参数等的粘弹性测量的处理流 程进行说明。首先,在步骤S00中开始粘弹性测量。开始的信号使用省略 了图示的输入设备、例如计算机的键盘和鼠标被输入。粘弹性测量开始前, 通常的B模式像被显示于显示部7中。步骤10中,输入粘弹性测量中的 测定信息。测定信息为位移产生用的超声波波束距体表面的深度(即,集 点F),F值(=焦点距离/开口直径),载波的频率等的值。位移生成用超 声波波束的焦点的位置、F值、载波的频率等被输入到中央控制部3,进 而,焦点的位置、F值、载波的频率经由中央控制部3被输入到送波波束 条件设定部12。测定信息也可由未图示的输入设备输入,或者,在中央控 制部3中,从省略了图示的存储部的存储介质读出适于乳腺或前列腺、血 管、肝脏等的部位的值。
在步骤20中,决定粘弹性的测量区域(ROI:Region ofInterest)。ROI 采用输入设备被输入,经由中央控制部3被输入到测定范围设定部33中。 或者,适于乳腺和前列腺、血管、肝脏等部位的值在中央控制部3中,从 存储介质读出并被输入到测定范围设定部33中。测量区域ROI的大小采 用例如深度方向的宽度为10mm程度,与体表水平的方向的宽度为5~ 20mm程度的长方形。如果设定焦点的位置和ROI,则决定在剪切波传播 的位移(通常、数μm~数十μm)的检波中采用的光栅和光栅上的采样点。 在各光栅中,位移检测用波束的接收的PRF(Pulse Repetition Frequency: 反复发送的脉冲的频率)设定为相对于假设的剪切波的频率满足奈奎斯特 定理。例如,在光栅与剪切波的位移的方向相同的情况下,设PRF为剪切 波的频率的2倍以上。
在步骤30中,采用在步骤10中所决定的测定条件并从超声波探头1 照射位移生成用超声波波束,在由步骤20所决定的测量区域内检测剪切 波。而且,剪切波的位移、或者粒子速度由位移运算部32计算。该位移 运算部32作为位移检测部30的一部分进行了图示,但也能通过上述的 CPU中的程序处理来实现。步骤40中,采用进行了检测的剪切波的位移 或者粒子速度,由粘弹性解析部34进行粘性参数的解析、估计。在步骤 50中,粘性参数的解析、估计结果被显示于显示部7中,本实施例的结构 中的粘弹性测量结束。
在此,采用图5A、图5B对图4中的步骤30的剪切波的检测的处理 流程进行详细说明。在此,对使线性阵列型的超声波探头1与被检测体的 体表面接触,在体内的目的断层面会聚位移生成用超声波波束、即,推力 波束,或者,推力脉冲的情况进行说明。此外,对在期望的断层面内,位 移生成用超声波波束的传播方向处于相对体表垂直的方向的情况进行说 明。在本实施例的结构中,检测对剪切波的位移进行时间微分的剪切波的 粒子速度。
在图5A的情况下,首先,从第二超声波收发部31照射位移检测用超 声波波束,取得参考信号。接下来,将位移生成用超声波波束即推力波束、 推力脉冲照射到焦点。接下来,由于检测通过推力波束、推力脉冲的照射 所生成的剪切波,因此位移检测用超声波波束被作为跟踪脉冲照射。在此, 各光栅中对位移检测用超声波波束进行照射的重复频率为PRF,跟踪脉冲 的检测个数的总计为N个(N为任意的整数)。接着,进行参考信号与已 进行跟踪的RF数据的相关运算51,算出剪切波的位移。即、图5A中, 粒子速度通过由位移检测部30等在不同的时刻所接收的多个信号之间的 相关运算而被计算出。相对于采用相关运算51所算出的位移,作为边缘 强调滤波器在时间方向上实施HPF(High Pass Filter)处理,算出粒子速 度。另外,也可代替HPF,采用用于进行时间微分的边缘强调滤波器来得 到粒子速度。
图5B中表示其他的粒子速度的算出方法。图5B中,在取得作为位移 检测用超声波波束的跟踪脉冲1,2,…,N为止总计N个之后,对第i 个和第i+1个(i为1~N-1为止的整数)的高频(RF)数据进行相关运算 52。即、图5B中,粒子速度通过由位移检测部30等在不同的时刻所接收 的多个信号之间的相关运算而被计算。在图5B的方法中,由于能不采用 HPF等的滤波器而直接算出粒子速度,因此能够缩短计算时间。
在此,对检测剪切波的粒子速度的优点进行说明。图6A、图6B表示 在剪切波的位移或者粒子速度的时间波形中叠加低频噪声、例如1Hz的正 弦波的情况下的时间波形。所谓低频噪声为因手术执行者的手抖动或体动 而产生的噪声。将不同的多个测量位置、即,光栅中的位移的时间波形的 峰值表示为图6A的P1、P2、P3。可知为了粘弹性测量,需要对峰值高精 度地进行检测,但随着与P1、P2、P3成为峰值的时刻延迟,难以识别位 移的峰值位置。
另一方面,在图6B中粒子速度的时间波形的峰值表示为P22、P23、 P21’、P22’、P23’。粒子速度的时间波形针对各光栅,检测一个或者两个 峰值。可以明确图6A中所示的位移的时间波形的峰值不同,在图6B中 所示的粒子速度的时间波形中,即使在低频噪声叠加的情况下也能稳定地 检测峰值,这是对粒子速度进行检测的优点。
接下来,对步骤40中的粘性参数的解析估计方法的详细情况进行说 明。图7A、图7B中表示与图6B相对应的粒子速度的时间波形。针对不 同的测量位置上的粒子速度的波形W21、W22、W23,存在成为1个或者 2个峰值的时刻。图7A的T22、T23是成为W22与W23的负峰值的时刻。 此外,图7A的T21’、T22’、T23’是成为W21、W22、W23的正峰值的时 刻。另一方面,图7B的T21”、T22”、T23”是W21、W22、W23与零相 交叉的、即,过零的时刻。
接下来,采用以上所说明的粒子速度的时间波形的峰值,或者,成为 峰值的时刻来对粘性参数进行估计。首先,对粘性与峰值、以及成为峰值 的时刻之间的关系进行说明。图8表示作为示出生物体的粘弹性的模型之 一的Voigt模型(model)80。Voigt模型80由剪切弹性系数μ和剪切粘 性系数η这两个要素构成。除了Voigt模型80之外,也可考虑Zener model 等这3个要素的粘弹性模型。
如图9A所示那样,在生物体由图8所示的Voigt模型表示时,剪切 波的传播速度具有频率分散性。图9B和图9C,表示在没有剪切粘性系数 η的情况下,以及粘性系数η为1.5Pa·s的情况下的剪切波的粒子速度波 形的一例。时间波形为从推动位置、即,从焦点起沿着剪切波的传播方向 的距离x=4、7、10mm的地点处的波形。如图9B所示那样,在没有粘性 的情况下,即使测量位置距推动位置远离4、7、10mm,粒子速度在时间 方向的波形也不会变化。另一方面,如图9C所示那样,在有粘性的情况 下,随着测量位置远离推动位置,成为正的峰值和负的峰值的时刻扩大。 此外,可知如果具有粘性则伴随着传播的峰值的减小程度变大。
即、如图10A示意地所示那样,粘性越大,成为正的峰值101与负的 峰值102的时刻之差越大。此外,如图10B所示那样,可知粘性越大,表 示峰值的振幅与距离的距离衰减特性的倾斜度越大。
接下来,采用图11A、图11B,对本实施例中的粘性参数的估计方法 进行说明。图11A表示测量位置xn的粒子速度的时间波形wn、以及粒子 速度波形的时间广度。图11A为示意性表示n=2时的例子、即在两个位置 检测出的两个粒子速度波形的时间广度的图。在该图中,将各时间波形 w1、w2下的成为正的峰值与负的峰值的时间之差设为t(x1)和t(x2)。 根据该时间的差t(x1)和t(x2),能够算出粒子速度波形的时间广度。
图11B中表示作为剪切波的粒子速度波形的时间广度的一例的、采用 正的峰值和负的峰值的时间差t(xn)、用于粘性参数的估计解析的处理流 程。首先,作为粒子速度波形的时间广度,输入至少两个以上的多个t(xn) (n为整数)110。接下来,通过FFT111所执行的运算,算出粒子速度波 形的中心频率ω。即、在粘弹性解析部34中,算出粒子速度的中心频率。 位置信息xn、成为正的峰值和负的峰值的时间的差t(xn)相对于中心频 率ω1,ω2,ω3,...,在粘弹性模型中以特有的函数形式提供。例如,函 数以表示指数曲线的、A×exp(α·xn)的形式提供。因而,基于所检测的 中心频率ω、xn、t(xn)的信息,进行与函数的拟合(fitting),执行粘性 参数α的评估112。通过粘性参数α的评估112下的拟合,估计参数A和 粘性参数α。
即、在粘弹性解析部34中,基于在多个位置上算出的粒子速度波形 的时间广度,采用指数曲线等并通过拟合来估计粘性参数。
而且,输出例如通过与中心频率ω拟合所求得的粘性参数α。在此所 利用的函数从省略了图示的记录介质经由中央控制部3,被读入到粘弹性 解析部34。该函数除了exp之外,还可采用多项式函数等公知的函数。进 而,基于作为粘性参数α与剪切波生成的超声波波束的条件的、F值、焦 点距离、载波的频率等的参数,与粘弹性模型拟合113,估计voigt model 的粘性系数η。在此,粘性系数η与通过拟合所求得的粘性参数α的关系 也可作为查询表(Look Up Table:LUT)或者函数,经由中央控制部3从省 略了图示的存储部的存储介质读出。而且,剪切粘性系数η被输出114。
接下来,对本实施例中的其他粘性参数的决定方法进行说明。图12 表示测量位置x1与x2处的正的峰值p(x1)、p(x2)以及负的峰值p(x1)’、 p(x2)’。在该决定方法中,代替图11B的t(xn),采用粒子速度的p(xn) 或者p(xn)’。而且,代替粘性参数α和参数A,通过与表示指数曲线的、 B×exp(β·xn)的函数相拟合,从而估计参数B和粘性参数β。
图13A、图13B表示以p(xn)为输入130进行拟合时的例子。通过 FFT131中的中心频率ω的检测、以及粘性参数β的评估132,来输出粘 性参数β和中心频率ω作为粘性参数的输出值。进而,基于粘性参数β和 作为剪切波生成的超声波波束的条件的F值、焦点距离、载波的频率等的 参数,通过拟合133,估计图8所示的Voigt Model的粘性系数η。在此, 粘性系数η与通过拟合所求得的粘性参数β的关系也可作为LUT或者函 数,经由中央控制部3从省略了图示的存储部即存储介质读出。而且,剪 切粘性系数η被输出134。
另外,作为以上说明的用于算出粒子速度波形的时间广度的t(xn)、 p(xn)、p(xn)’以外的输入参数,也可采用粒子速度波形的积分值,成 为过零值的时刻前后的粒子速度的积分值、比成为零的时刻晚的时刻与比 成为零的时刻早的时刻的粒子速度波形的积分值之差、比任意的阈值大或 者小的粒子速度的积分量、成为正的峰值的时刻与成为零的时刻之差,成 为负的峰值的时刻与成为零的时刻之差等。即、也可代替根据粒子速度为 正和负的两个时刻来算出剪切波的粒子速度波形的时间广度,而根据成为 粒子速度的过零值的时刻、或者成为过零值的时刻前后的粒子速度的积分 值等算出剪切波的粒子速度波形的时间广度。此外,在图9A~图12中, 负的峰值在比正的峰值早的时刻被检测出,但在相反的情况下也同样能决 定粘性参数。
图14A、图14B中表示在本实施例的超声波诊断装置的显示部7中显 示的粘性参数α的例子。例如,如图14A所示那样能够对粘性参数α与中 心频率ω处的剪切波的传播速度Vs(ω)进行标绘显示141。传播速度 Vs(ω),在粘弹性解析部34中,采用公知的方法被算出。除了粘性参数 α之外,还能显示粘性参数β和剪切粘性系数η。或者,也可如图14B那 样与在显示部7中显示的B像叠加,进行2D彩色地图142显示。此时, 将粘性参数α,β和剪切粘性系数η的彩色标度143同时显示于显示部7 是有效的这一点如先前所说明的那样。
实施例2
第2实施例为涉及以下的超声波诊断装置以及方法的实施例,即超声 波诊断装置具有:对被检测体内施加声音放射压来使被检测体内产生位移 的位移生成部10;通过对被检测体的至少两个以上的多个位置发送接收超 声波从而检测在被检测体内产生的剪切波的粒子速度的位移检测部30;根 据粒子速度估计粘性参数的粘弹性解析部34;以及反馈参数决定部38, 粘弹性解析部34根据粒子速度波形的时间广度来估计粘性参数,反馈参 数决定部38基于通过第一测量(第一次测量)得到的、根据粒子速度波 形的时间广度所估计的粘性参数,对接下来的测量即第二测量中的发送条 件、或者测量区域中的至少一个进行调整。
以下,对与使用反馈参数决定部38进行反馈控制的超声波诊断装置 相关的实施例2进行说明。图15中表示实施例2的超声波诊断装置的整 体结构例,但在图1所示的结构的超声波诊断装置中新加入反馈参数决定 部38。该反馈参数决定部38与粘弹性解析部34同样,能够通过CPU中 的程序执行来实现。由于其他构成要素的模块与图1所示的实施例1的超 声波诊断装置的结构相同,因此省略重复的说明。
图16中表示实施例2的超声波诊断装置中的粘性参数的计算时序。 在实施例2中,将步骤S40中解析估计出的粘性参数的值反馈到步骤S20 中的剪切波的测量区域决定。进行反馈的值由反馈参数决定部38算出。 所算出的值经由中央控制部3而被输入至位移检测部30的测定范围决定 部33。
图17A中表示包括本实施例的反馈控制的、粘弹性测量的处理流程的 概要。例如如图5B中所说明的那样照射推力脉冲1和跟踪脉冲1,在剪 切波粒子速度的检测、粘性参数的估计之后,基于所估计的粘性参数的值, 进行测量区域ROI的重新设定171。由于粘性越大,伴随着传播的剪切波 的振幅值也越小,因此能测量的区域也越小。因而,进行重新设定以使粘 性的值越大,测量区域ROI的大小越小。接下来,照射推力脉冲2和跟踪 脉冲2,在进行了重新设定的测量区域ROI内再次进行粘弹性测量。
如图17B所示那样,叠加于B像的测量区域172的大小的设定例如采 用梯形的角度Θ来进行设定。即、在本实施例中,测量区域ROI的形状 的最佳的例子为梯形,粘性越大,越增大角度Θ的值。除了梯形的角度Θ 之外,也可设为测量区域ROI的横向的宽度、深度方向的宽度。这些值针 对焦点F而被设定。即、反馈参数决定部38使进行调整的测量区域ROI 的大小对应于进行测量的位置而产生变化。另外,该测量区域ROI的重新 设定值能基于所估计的粘性参数,从省略了图示的存储部经由中央控制部 3读出。
图18中示意地表示本实施例的反馈参数决定部38中的测量区域ROI 的决定过程。采用Vs值、中心频率、粘性参数α、β、粘性系数η等参数 181来进行ROI控制182,由此决定角度Θ。在决定角度Θ时,也可还采 用F值、载波频率、焦点距离等的位移生成用超声波波束的照射条件、即、 由虚线所示的推力脉冲送信条件、推力波束的照射条件,来决定角度Θ。 此时,如图17B所示那样,随着B像的测量深度,测量区域ROI的形状 发生变化。因此,通过将测量区域ROI的形状显示于显示部7,从而手术 执行者能够在视觉上确认所决定的测量区域ROI。
根据本实施例,通过进行反馈控制能实现粘性参数估计、测量的高精 度化。上述的测量区域ROI的形状除了为梯形之外,也可为任意的几何形 状。此外,进行反馈的值除了测量区域ROI之外,也可设为推力脉冲的照 射条件即振幅值、F值、载波的频率等。通过这些值的反馈,能够选择减 小了振幅值和测量区域ROI后的重视安全性的粘弹性测量、或增大了振幅 值和测量区域ROI后的重视测量区域的粘弹性测量等。
进而,实施例2相关的反馈控制的接通/断开也可经由省略了图示的输 入设备,来切换手术执行者。此外,能够设为具有进行反馈参数决定部38 所执行的反馈控制的接通/断开的切换的开关的构成。进而,手术执行者也 可采用省略了图示的输入设备来进行反馈控制的参数的选择。或者,也可 采用测定对象、测量时刻的信息,由中央控制部3自动地进行设定。
尽管在上述所有的实施例中,采用了超声波会聚波束生成辐射压并产 生剪切波,但除了辐射压之外,也可采用由DC电动机、振动泵等机械驱 动、手动加压、电脉冲所产生的加压等公知的方法。
此外,也可代替线性阵列型的探头,而采用凸状类型或扇形类型的探 头,也可采用2维探头。线性阵列型的探头优选用于乳房或腱/肌肉等中。 此外,凸状类型优选用于肝脏,扇形类型优选用于心脏。此外,超声波探 头1的各元件采用例如以陶瓷、高分子作为材料的压电元件或以硅作为材 料的静电电容型的振荡器,但各元件并不限于此。如果采用2维探头,则 能三维地检测剪切波,能更明确地显示剪切波的传播。由此,期待评估性 能提高或者评估时间缩短的效果。
另外,在对粘性进行测量的区域内存在多个空间性的弹性的情况下, 所测量的粘性参数或粘性系数也成为被平均化的值。而且,在进行粘性检 测的区域,如本申请发明者等在先前申请的国际公开WO2012/105152号 公报中所记载的那样,也可选定变形信息均匀的区域。在变形信息均匀的 区域中,由于弹性率均匀,因此能高精度地估计粘性参数和粘性系数。测 定对象例如可举出肝脏、胸部、血管、前列腺、韧带、腱、肌肉等。
根据以上详述的本发明,由于能评估粘性参数,因此能客观地评估肿 瘤等组织性状或者判定治疗效果。此外,由于生成采用了脉冲波的声音放 射压,因此在每一次测量时所需的声音放射压的施加时间为0.5ms~1ms, 能将生物体组织的温度上升抑制到最小限度。此外,通过采用粒子速度波 形,从而能在体动等中进行鲁棒性测量。进而,基于粘性参数,重新设定 对剪切波进行测量的区域,能够再次评估粘性参数,能进行高精度的测量。
另外,本发明并不限于上述的实施例,包括各种变形例。例如,为了 更好地理解本发明而详细地说明了上述的实施例,但不必限定于具备说明 的所有的结构的内容。此外,能将某实施例的结构的一部分置换为其他实 施例的结构,此外,能在某实施例的结构中追加其他实施例的结构。此外, 关于各实施例的一部分结构,能进行其他结构的追加/删除/置换。
进而,上述的各结构、功能、处理部等采用作成实现其中的一部分或 者全部的程序的例子进行了说明,但当然也可通过采用例如集成电路进行 设计等以硬件来实现其中的一部分或者全部。

超声波诊断装置以及方法.pdf_第1页
第1页 / 共30页
超声波诊断装置以及方法.pdf_第2页
第2页 / 共30页
超声波诊断装置以及方法.pdf_第3页
第3页 / 共30页
点击查看更多>>
资源描述

《超声波诊断装置以及方法.pdf》由会员分享,可在线阅读,更多相关《超声波诊断装置以及方法.pdf(30页珍藏版)》请在专利查询网上搜索。

本发明涉及超声波诊断装置以及方法。在生物体等对象组织为粘弹性体的情况下,剪切波的传播速度随着剪切波的频率而发生变化。超声波诊断装置通过来自位移生成部(10)的信号从超声波探头(1)照射位移生成用送波波束而产生剪切波之后,通过位移检测部(30)对生物体组织照射脉冲波,检测剪切波的粒子速度,根据由粘弹性解析部(34)检测得到的剪切波的粒子速度的波形的时间广度来估计粘性参数,并在显示部(7)中显示估计出。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 人类生活必需 > 医学或兽医学;卫生学


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1