CN102695457A - 超声波诊断装置及测量内中膜的厚度的方法 - Google Patents

Abstract

本发明提供一种超声波诊断装置及测量内中膜的厚度的方法。超声波诊断装置向被检体发送超声波，并基于接收到的来自被检体的回波信号获取生物体信息。超声波诊断装置具备：强度分布生成部，根据回波信号生成与血管壁相关的深度方向的回波强度分布；模板生成部，基于预先准备的表示回波强度分布的基准式样的模板，生成用于边界检测的模板；式样类似度计算部，计算模板与强度分布之间的类似度，即：使模板在深度方向上变化、且使与模板的内中膜厚对应的式样变化，从而按每个变化来计算模板与强度分布之间的类似度；和边界决定部，基于按每个变化计算出的类似度来决定内腔边界及外膜边界。在内中膜的厚度测量的内腔边界与外膜边界的检测中，不依赖于内腔侧和外膜侧的回波强度变化的大小关系。

Description

超声波诊断装置及测量内中膜的厚度的方法

技术领域

本发明涉及计算例如颈动脉等的血管的内中膜复合体的厚度的超声波诊断装置。

背景技术

一般利用超声波诊断装置进行动脉硬化的诊断，该超声波诊断装置向被检体收发超声波，并基于来自所述被检体的回波信号获取生物体信息。周知在动脉硬化的诊断中，将颈动脉的内中膜复合体厚作为了解初始的粥样硬化的重要指标。“内中膜复合体”或者“内中膜”是存在于血管内腔与外膜之间的、血管壁的内膜及中膜的总称。并且，内中膜复合体厚或者内中膜厚也被称为“IMT”(Intima-Media Thickness)，在本申请说明书中记为“IMT”。

图1示出存在于血管内腔1和外膜2之间的内中膜3。在检查中，检测作为该血管内腔1与内膜的边界的内腔边界4、以及作为中膜与外膜2的边界的外膜边界5，并测量其厚度。作为自动测量该IMT的方法，例如已知专利文献1的技术。

图2示出在专利文献1的技术中在确定外膜边界及内腔边界时所利用的回波强度分布。在专利文献1中公开了如下的方法，如图2所示那样，在深度方向的回波强度分布中，将强度变化在规定值以上的位置作为边界候选，将较大的一方定为外膜边界，将较小的一方定为内腔边界。

现有技术文献

专利文献

专利文献1：JP特开2008-168016号公报

发明内容

发明所要解决的技术问题

但是，在专利文献1中，由于根据强度变化大或强度变化小统一进行处理，因此有时无法测量出正确的IMT。

例如，图3表示内腔边界的强度变化大于外膜边界的强度变化的例子。是强度变化大的位置未必始终是外膜边界的例子。

其结果，在专利文献1所公开的方法中，会将内腔边界作为外膜边界进行检测，故无法测量出正确的IMT。

本发明是鉴于上述问题而提出的，其目的在于正确地检测内腔边界和外膜边界，并实现正确的IMT测量，而不依赖于内腔侧与外膜侧的回波强度变化的大小关系。

用于解决技术问题的方法

本发明的超声波诊断装置向被检体发送超声波，并基于接收到的来自所述被检体的回波信号获取生物体信息。本发明的某个所述超声波诊断装置具备：强度分布生成部，其根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；模板生成部，其基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；式样类似度计算部，其是计算所述模板与所述强度分布之间的类似度的式样类似度计算部，即：使所述模板在所述深度方向上变化、且使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；和边界决定部，其基于按每个所述变化所计算出的所述类似度，来决定内腔边界及外膜边界。

本发明的另一个超声波诊断装置具备：强度分布生成部，其根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；模板生成部，其基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；式样类似度计算部，其是计算所述模板与所述强度分布之间的类似度的式样类似度计算部，即：针对在所述深度方向的不同的位置处存在的边界的各候选的组合，使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；边界连续度计算部，其针对关于相邻的多个声音线而确定的所述组合，基于各声音线的回波信号的强度值的差值来计算边界连续度；边界评价部，其针对所述组合，生成综合所述式样类似度和所述边界连续度之后的评价值；和边界决定部，其基于综合之后的所述评价值，来决定内腔边界和外膜边界。

也可，所述模板生成部所生成的模板至少包含3个区域，在所述至少3个区域中至少包含1个能够改变所述区域的范围的长度可变区域。

也可，所述至少3个区域按照顺序排列，所述长度可变区域是所述至少3个区域之中的中央的区域。

也可，所述至少3个区域的一个端部的区域对应于低回波的组织区域，所述至少3个区域的另一个端部的区域对应于高回波的组织区域。

也可，所述低回波的组织区域对应于血流区域，所述高回波的组织区域对应于血管外膜壁及体组织。

也可，所述长度可变区域对应于血管的内中膜区域。

也可，所述式样类似度计算部计算使所述长度可变区域的长度标准化之后的类似度。

也可，所述式样类似度计算部将所述长度可变区域的平均强度用作所述长度可变区域的强度，由此使所述长度可变区域的长度标准化。

也可，所述模板生成部所生成的模板由2个种类的子模板构成，所述2个种类的子模板各自由至少2个区域构成，并能够改变2个种类的模板的距离，所述式样类似度计算部利用所述2个种类的子模板分别计算类似度，并将合起且综合两者之后的值作为类似度。

也可，构成所述2个种类的模板的2个区域分别对应于低回波及高回波的组织区域。

也可，所述2个种类的模板的低回波及高回波的组织区域分别对应于血流区域及内中膜区域、或者分别对应于内腔边界及规定的组织区域，所述规定的组织区域是外膜边界、动脉外膜壁及体组织。

也可，所述模板生成部根据所述强度分布生成部生成的强度分布的对比度，来改变所生成的模板的系数值。

也可，所述式样类似度计算部所评价的类似度通过由所述模板生成部生成的模板而指定的边界间的强度差来定义。

也可，所述式样类似度计算部所评价的类似度通过所述模板生成部生成的模板与所述强度分布生成部生成的强度分布之间的标准化相关来定义。

也可，所述超声波诊断装置还具备后壁检测部，该后壁检测部根据所述强度分布生成部生成的强度分布检测高回波的区域，所述边界决定部缩减为所检测出的后壁范围内来决定边界。

也可，所述后壁检测部所检测的高回波区域对应于血管壁。

也可，所述边界评价部利用动态规划法计算所述评价值。

也可，所述边界评价部利用降低了分辨率的强度分布来计算所述评价值，所述边界决定部基于降低了分辨率的所述强度分布来暂时决定候选，然后在原分辨率的强度分布中从所述临时决定出的候选的周边的区域中决定内腔边界和外膜边界。

也可，所述边界评价部利用动态规划法缩减为所述评价值良好的上位候选、或者具有规定的评价值的候选，同时在声音线方向上计算所述评价值。

也可，所述边界连续度计算部计算由声音线方向的边界位置的差所定义的边界连续度。

也可，所述式样类似度计算部、所述边界连续度计算部、所述边界评价部、及所述边界决定部以内腔边界及外膜边界这两个边界作为对象，来评价所述可取得的边界候选，所述边界决定部同时决定所述内腔边界及外膜边界。

也可，所述边界连续度计算部计算边界连续度，该边界连续度由在声音线方向的内腔边界及外膜边界的位置处确定的厚度的差来定义。

也可，所述边界决定部在基于确定的顺序来决定内腔边界和外膜边界时，所述式样类似度计算部在另一个边界位置尚未确定的阶段计算类似度的情况下，使一个边界固定，针对具有与另一个边界位置对应的多个厚度的模板来计算类似度，将计算出的所述类似度的最大值、或者平均值作为所述一个边界的类似度。

也可，所述边界决定部，在检测出内腔边界之后缩减为所述内腔边界的下方并检测出外膜边界的情况下、和在检测出外膜边界之后缩减为所述外膜边界的上方并检测出内腔边界的情况下，比较在各个检测中用到的评价值，将评价值良好的检测边界作为最终的边界采用。

也可，所述超声波诊断装置还具备IMT计算部，该IMT计算部根据所述内腔边界和外膜边界计算内中膜的厚度。

本发明的方法是向被检体发送超声波，并基于接收到的来自所述被检体的回波信号测量血管壁的内中膜的厚度的方法，其中，该方法包括如下步骤：根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；是计算所述模板与所述强度分布之间的类似度的步骤，即：使所述模板在所述深度方向上变化、且使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；基于按每个所述变化所计算出的所述类似度，来决定内腔边界及外膜边界；和根据所决定的所述内腔边界及所述外膜边界，计算内中膜的厚度。

本发明的方法是向被检体发送超声波，并基于接收到的来自所述被检体的回波信号测量血管壁的内中膜的厚度的方法，其中，该方法包括如下步骤：根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；是计算所述模板与所述强度分布之间的类似度的步骤，即：针对在所述深度方向的不同的位置处存在的边界的各候选的组合，使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；针对关于相邻的多个声音线而确定的所述组合，基于各声音线的回波信号的强度值的差值来计算边界连续度；针对所述组合，生成综合所述式样类似度和所述边界连续度之后的评价值；基于综合之后的所述评价值，来决定内腔边界和外膜边界；和根据所决定出的所述内腔边界及所述外膜边界，计算内中膜的厚度。

发明效果

根据本发明，并不是独立地检测内腔边界和外膜边界，而是将2个边界作为组来进行检测，即便在内腔边界的强度变化大于外膜边界的强度变化的情况下，也能够正确地测量IMT。

附图说明

图1是表示存在于血管内腔1与外膜2之间的内中膜3的图。

图2是表示在专利文献1的技术中在确定外膜边界及内腔边界时所利用的回波强度分布的图。

图3是表示内腔边界的强度变化大于外膜边界的强度变化的例子的图。

图4是本发明的第一实施方式所涉及的超声波诊断装置10的结构框图。

图5是表示用于发送超声波脉冲的、内置于超声波探头11的超声波振动元件组6的图。

图6(a)及(b)是使用沿着x方向排列的多个超声波振动器形成焦点时的超声波集束波的示意图，(c)及(d)是超声波集束波的略图。

图7是表示IMT测量部16的结构的框图。

图8是表示IMT计算处理的步骤的流程图。

图9是表示系数式样(内中膜式样)的图。

图10是表示作为对象的声音线(line)的强度分布90和改变了深度的模板组91的图。

图11是表示作为对象的声音线(line)的强度分布90和改变了厚度的模板组92的图。

图12是表示利用与内腔边界和外膜边界分别匹配的子模板来计算类似度并通过加权相加等来求出类似度的例子的图。

图13是表示本发明的第二实施方式的IMT测量部30的结构的框图。

图14是表示基于动态规划法的边界检测处理的步骤的流程图。

图15是表示关于对象采样130预先规定的路径对象130a～130e的例子的图。

图16(a)～(c)是表示利用了本发明的第二实施方式中的动态规划法的成本计算方法的例子的图。

图17是表示路径成本C_path最小的内腔边界和外膜边界的对i_line的例子的图。

图18是表示针对内腔边界及外膜边界分别进行路径检索的处理的步骤的流程图。

图19是表示具有后壁检测部26的IMT测量部16a的结构的框图。

具体实施方式

以下，参照附图，说明本发明的超声波诊断装置的实施方式。

(第一实施方式)

图4是本发明的第一实施方式所涉及的超声波诊断装置10的结构框图。

超声波诊断装置10具备：探头11、发送部12、接收部13、B图像生成部14、ROI设定部15、IMT测量部16、合成部17、和显示部18。

发送部12生成超声波脉冲，并输出至探头11。

探头11将由发送部12输出的超声波脉冲发送至生物体内，并接收从生物体内反射回来的回波信号。探头11将接收到的回波信号输出至接收部13。

图5表示用于发送超声波脉冲的、内置于超声波探头11的超声波振动元件组6。超声波振动元件组6(以下记为“振动器6”。)构成例如各超声波振动元件沿着一个方向排列的、所谓的1D阵列振动器。

振动器6例如由压电体构成，通过驱动压电体来发送超声波，此外，压电体接收超声波并将超声波变换为电信号。振动器6使各超声波振动元件依次驱动来进行超声波的收发，能够扫描规定范围。此外，振动器6使来自多个超声波振动元件的各超声波的相位在规定的位置(焦点位置)处重合，也能够接收在焦点位置处反射的信号。图6表示后者的例子。

图6(a)及(b)示意地表示使用沿着x方向排列的多个超声波振动元件形成焦点时的超声波的集束波。超声波的集束波具有图示那样的规定宽度，在z轴方向的规定的深度处具有焦点。其中，该集束波的中心轴如图6(c)及(d)中所示那样也被称为“声音线”。

再次参照图4。接收部13使从探头11输出的回波信号成束后，输出至B图像生成部14和IMT测量部16。

B图像生成部14根据回波信号生成断面图像，将其断面图像的数据输出至合成部17。

ROI设定部15从用户接受确定ROI位置的指示，基于该指示在动脉壁的位置处设定关心区域(Region of Interest：ROI)。

IMT测量部16测量IMT。更为具体而言，IMT测量部16根据由接收部13输出的回波信号和由ROI设定部15输出的ROI位置的信息，检测内腔边界和外膜边界。进而，IMT测量部16基于检测出的内腔边界及外膜边界的位置来计算IMT值。IMT测量部16将确定边界位置的信息(边界位置信息)输出至合成部17，将IMT值输出至显示部18。

合成部17基于由B图像生成部14输出的B图像的数据和由IMT测量部16输出的边界位置信息，合成IMT值和B图像。更为具体而言，合成部17将内腔边界及外膜边界重叠在B图像上，并输出得到的合成图像数据。

显示部18将由IMT测量部16输出的IMT值和由合成部17输出的合成结果显示在监视器上。

显示部18将从合成部17输出的合成图像数据显示在例如液晶显示器等的显示设备上。

接下来，详细说明IMT测量部16。

图7是表示IMT测量部16的结构的框图。

IMT测量部16具有：强度分布生成部21、模板生成部22、式样(pattern)类似度计算部23、边界决定部24、IMT计算部25。

强度分布生成部21基于从接收部13输出的回波信号、及从ROI位置设定部15输出的ROI位置的信息，生成回波强度分布。然后，强度分布生成部21将所生成的强度分布输出至式样类似度计算部23。该回波强度分布是与从被检体的表皮向生物体内部的方向(深度方向)相关的回波强度分布。

模板生成部22基于预先准备并设定的、表示回波强度分布的基准式样的模板，生成用于边界检测的模板。模板生成部22将用于边界检测的模板输出至式样类似度计算部23。

式样类似度计算部23计算从强度分布生成部21输出的、接收信号的强度分布和从模板生成部22输出的模板之间的类似度。以下，将该类似度记为“式样类似度”。式样类似度计算部23将计算出的类似度输出至边界决定部24。

边界决定部24基于从式样类似度计算部23输出的式样类似度，决定内腔边界及外膜边界。边界决定部24将所决定的边界位置的信息输出至合成部18和IMT计算部25。

IMT计算部25基于所决定的内腔边界及外膜边界来计算IMT值，并将计算出的IMT值输出至显示部17。

接下来，说明本实施方式的IMT计算处理的步骤。

图8是表示IMT计算处理的步骤的流程图。

首先，在步骤S11中，强度分布生成部21对所接收的回波信号进行包络线检波、对数压缩，生成回波强度分布。该处理针对ROI内的各点进行。

接下来，在步骤S12中，强度分布生成部21为了降低噪声，针对强度分布应用例如3×3的平均值滤波器，进行平滑化。

接着，处理转移至式样类似度计算部23。再者，以下的处理按照每个声音线(以线为单位)来进行。

首先，在步骤S13中，式样类似度计算部23评价ROI中的强度分布和表示内中膜式样的模板之间的类似度。模板是由模板生成部22生成的。

作为模板，使用图9所示的系数式样(内中膜式样)。在图9中，范围l₁对应于血管内腔，范围l₂和范围l₃对应于内中膜，范围l₄对应于外膜。图9所示的模板是微分型滤波器，作为基准的系数式样预先在深度方向上交替设置取负值的l₁、l₃的范围、和取正值的l₂、l₄的范围。如图示那样，在模板中至少设置三个区域、即对应于血管内腔的区域(或者式样)、对应于内中膜的区域、及对应于外膜的区域。

式样类似度计算部23将该模板作为检测窗，如图10所示那样，一边在深度方向上移动一边类似度，搜索合适的深度。

此外，式样类似度计算部23还评价相对于内中膜厚的类似度。即、如图11所示那样，在模板生成部22中制作将内中膜厚改变成多样的多个模板，并计算这些的模板与强度分布之间的类似度，搜索合适的厚度。

以下，进行更为详细的说明。

图10表示作为对象的声音线(line)的强度分布90、和改变了深度的模板组91。模板组91的各模板是与声音线平行的深度方向的1列式样。在图中示出了深度p1、p2、p3、p4、……、p10的模板。再者，图示的各模板仅深度不同，其式样是相同的。

式样类似度计算部23首先计算作为对象的声音线(line)的强度分布90、和深度p1且内中膜厚q时的模板之间的类似度。

类似度由模板T与强度分布X之间的内积值给出，在将深度p的强度分布设为X_p，将模板的内中膜厚设为q，将对应于内中膜厚q的模板设为T_q，将模板尺寸设为L时，类似度Y通过以下的公式计算。

[公式1]

$Y(p,q)=\underset{i=1}{\overset{L}{\mathrm{\Σ}}}{X}_p\left(i\right)\·T_q\left(i\right)$

在图10所示的例子中，在最合适的深度p4处，模板中包含的内中膜厚q的式样与对象线的强度分布90中的内中膜厚的式样相一致。但是，这仅是便于理解的一例。在实际中，在最合适的深度处，模板的内中膜厚q的式样与对象线的强度分布90中的内中膜厚的式样未必一致。

因此，在本实施方式中，接下来，式样类似度计算部23在深度p1处一边改变内中膜厚q的式样，一边计算类似度。图11表示作为对象的声音线(line)的强度分布90和改变了厚度的模板组92。式样类似度计算部23一边改变内中膜厚q，一边求得与各个厚度对应的类似度Y。

接着，式样类似度计算部23将深度变更至p2，一边改变内中膜厚q一边求得类似度Y。

这种的处理持续至深度p10。

在图8的步骤S14中，边界决定部24搜索类似度Y(p、q)最大的深度p_max和内中膜厚q_max，在步骤中S15中，分别决定内腔边界li和外膜边界ma。即，若将模板的血管内腔长度设为l1，则由以下的公式给出。

[公式2]

li＝p_max+l1

ma＝li+q_max

以上是以线为单位的处理。

在上述的处理中，系数值a₁’、a₂’、a₃’、a₄’也可以设定为所对应的范围长l₁、l₂、l₃、l₄的倒数值，以免类似度缩减为特定的内中膜厚。即，在与范围长l₁、l₂、l₃、l₄对应的各区域中，计算平均强度，成为加权的形式。此时的类似度由以下的公式计算。

[公式3]

$Y(p,q)=-\frac{{a_1}^{\′}}{l_1}\·\underset{i=1}{\overset{l_1}{\mathrm{\Σ}}}{X}_p\left(i\right)+\frac{{a_2}^{\′}}{l_2}\·\underset{i=l_1+1}{\overset{l_1+l_2}{\mathrm{\Σ}}}{X}_p\left(i\right)-\frac{{a_3}^{\′}}{l_3}\·\underset{i=l_1+l_2+1}{\overset{l_1+l_2+l_3}{\mathrm{\Σ}}}{X}_p\left(i\right)+\frac{{a_4}^{\′}}{l_4}\·\underset{i=l_1+l_2+l_3+1}{\overset{l_1+l_2+l_3+l_4}{\mathrm{\Σ}}}{X}_p\left(i\right)$

在此，a₁’、a₂’、a₃’、a₄’是权重系数。

再有，在改变内中膜厚时，既可以使l₂固定而使l₃可变，也可以使l₂、l₃基于规定的比例都可变。

步骤S15是切换至下一线的切换处理。以下，针对各线进行上述的处理。

以线为单位的处理结束之后，在步骤S16中，在线方向上对边界位置进行平滑化。

以上是第一实施方式的动作流程。

再者，在上述说明中，将类似度设定为强度分布与取得微分型滤波器形式的模板之间的内积值，但也可以是相关值。此时，在模板的系数值中不需要包含负值。

此外，如图12所示那样，类似度也可以设定为：准备分别与内腔边界和外膜边界匹配的子模板，分别计算类似度，通过加权相加等而组合之后的值。在各子模板中，规定了与边界及除此之外的组织等对应的式样，但这些的式样是可变的。即，与边界对应的式样的位置是可变的。上述的2种模板中的、与边界区域及除此之外的组织等对应的区域，分别对应于各高回波及低回波的组织区域。

此外，模板的系数值也可以基于ROI中的描绘度进行变更，例如，可以计算强度值的方差，在方差大的情况下，增强模板的对比度，相反在方差小的情况下，减弱模板的对比度。

最后，在上述说明中，假定直接检测内腔边界和外膜边界，但也可以首先使用微分型滤波器等检测后壁，缩减为后壁内部来检测内腔边界和外膜边界。这在血管中出现多重回波的情况下可有效地发挥作用。

以上，在第一实施方式中，确定与对象线的强度分布最类似的模板，以检测内腔边界和外膜边界。模板规定了包含内腔边界及外膜边界的强度分布的式样。这意味着并不是独立地检测内腔边界和外膜边界，而是在考虑内腔边界的强度变化及外膜边界的强度变化的同时来检测内腔边界及外膜边界。由此，在内腔边界的强度大于外膜边界的情况下，也不会发生误检测。

(第二实施方式)

本实施方式与第一实施方式的不同点在于，本实施方式中的IMT测量部的处理。具体而言，本实施方式中的IMT测量部在线方向上进行平滑化之前，除了考虑式样类似度以外，还考虑线方向的连续性来选择边界。如果进行该处理，则特别是即便在对应于内腔边界的强度不明确的情况下，也能够检测稳定的边界。

本实施方式的超声波诊断装置的结构及动作除了IMT测量部以外，与图4所示的超声波诊断装置10相同。因此，作为本实施方式的超声波诊断装置的结构援引图4，并且省略IMT测量部以外的说明。

图13是表示本实施方式的IMT测量部30的结构的框图。

IMT测量部30具有：强度分布生成部31、模板生成部32、式样类似度计算部33、边界连续度计算部34、边界评价部35、边界决定部36、和IMT计算部37。

强度分布生成部31基于从接收部13输出的回波信号、及从ROI位置设定部15输出的ROI位置的信息，生成回波强度分布。并且，强度分布生成部31将所生成的强度分布输出至式样类似度计算部23。

模板生成部32基于预先设定的模板(表示基准式样)，生成用于边界检测的模板。模板生成部32将用于边界检测的模板输出至式样类似度计算部33。

式样类似度计算部33接受从强度分布生成部31输出的接收信号的强度分布、和从模板生成部22输出的模板。进而，式样类似度计算部33从边界评价部35接受与内腔边界的位置及外膜边界的位置相关的指示。式样类似度计算部33基于由边界评价部35指定的内腔边界的位置及外膜边界的位置，计算该位置处的接收信号的强度分布、与从模板生成部22输出的模板之间的类似度。式样类似度计算部33将计算出的类似度输出至边界评价部35。

边界连续度计算部34从边界评价部35接受与内腔边界的位置及外膜边界的位置相关的指示。边界连续度计算部34针对所指定的内腔边界的位置及外膜边界的位置，利用对象线和相邻线来评价线方向的边界位置和厚度的连续性。边界连续度计算部34将评价结果作为“连续度”输出至边界评价部35。

边界评价部35基于由式样类似度计算部33输出的式样类似度、和由边界连续度计算部34输出的边界的连续度，针对预先规定的全部边界的组合，评价内腔边界相似度、及外膜边界相似度。评价值的计算基于后述的动态规划法来进行。边界评价部35将计算出的评价值输出至边界决定部36。

边界决定部36基于由边界评价部34输出的评价值，决定内腔边界及外膜边界。边界决定部36将决定出的边界位置输出至合成部18及IMT计算部37。

IMT计算部37基于所决定的内腔边界及外膜边界计算IMT值，并将计算出的IMT值输出至显示部17。该处理与图7所示的IMT计算部25的处理相同。

接下来，说明本实施方式的边界检测处理的步骤。下述的处理主要由边界连续度计算部34、边界评价部35及边界决定部36来进行。

图14是表示基于动态规划法的边界检测处理的步骤的流程图。

在本实施方式中，按照搜索由式样类似性和线方向的边界的连续性所定义的成本最小的线方向的路径(＝边界)的形式来实施边界检测处理。

在步骤S21中，边界评价部35针对所有的内腔边界和外膜边界的对i_line使路径成本C_path初始化。初始值由与内中膜式样之间的类似度C_pattern给出。若将权重设为W_pattern，则路径成本C_path由以下的公式给出。

[公式4]

C_path(i_line)←-W_pattern·C_pattern(i_line)

在此，类似度C_pattern是与第一实施方式同样的值。但是，主要想关注的是在线内进行标准化之后的值。不过，不需要在所有的线内进行标准化。例如，可以至少仅对与内中膜对应的线内的值进行标准化。

所谓“内腔边界和外膜边界的对”相当于内腔边界和外膜边界的组合。若将ROI的采样点设为N点时，则对的数量为N的平方。不过，由于内腔边界处在比外膜边界浅的位置、或者系统规定的最大测量厚等的约束条件，实际计算中使用的对数少于N的平方。

接下来，进入至第2线。

然后，在步骤S22中，式样类似度计算部33计算与针对前一线的各内腔边界和外膜边界的对i_line-1的内中膜式样之间的类似度。此外，边界连续度计算部34针对前一线的各内腔边界和外膜边界的对i_line-1，分别计算内中膜厚和边界位置的连续度。边界评价部35接受这些计算结果，评价线方向的内中膜厚和边界位置的连续性。

与内中膜式样之间的类似度C_pattern，与第一实施方式同样既可以由强度分布和模板的内积值给出，也可以是相关值。

若将对象线和前一线的内腔边界位置分别设为LI_line、LI_line-1，将外膜边界设为MA_line、MA_line-1，则内中膜厚的连续度C_Thickness由以下的公式给出。

[公式5]

C_thickness(i_ine-1)＝|(MA_line-Ll_line)-(MA_line-1-Ll_line-1)|

此外，边界位置的连续度C_distance由以下的公式给出。

[公式6]

C_dis tan ce(i_line-1)＝|Ll_line-Ll_line-1|+|MA_line-MA_line-1|

边界评价部35在线内分别对评价值C_Thickness和C_distance标准化之后，如以下的公式那样将它们统一成候选路径的成本C_total。

[公式7]

C_total(i_line-1)＝C_path(i_line-1)+W_thickness·C_thickness(_line-1)+W_dis tan ce·C_dis tan ce(_line-1)

在此，权重W_thickness及W_distance是相对于各评价值的权重系数。再者，公式7中包含评价值C_Thickness和C_distance，但只要包含边界位置的连续度C_distance即可，可以没有内中膜厚的连续度C_Thickness。

边界评价部35针对预先决定的路径对象的内腔边界和外膜边界的对i_line-1实施以上的计算。

例如在图15中，表示关于对象采样130(例如内腔边界的候选)所预先决定的路径对象130a～130e。边界评价部35针对由对象采样130和各路径对象130a～130e之间的各个组合所确定的路径，计算上述的成本。

图16(a)～(c)表示利用了本实施方式中的动态规划法的成本计算方法的例子。在图16(a)中，在路径对象130a～130e的内部记载的数值是路径成本的累积值。在将对象采样130和各路径对象130a～130e连接起来的线段的附近记载的数值是距离的成本。并且，在对象采样130的右侧沿着深度方向记载的数值是式样的成本。

如图16(b)所示，边界评价部35在各路径对象130a～130e的路径成本的累积值上加上距离的成本之后的结果中、采用提供最小值(图16(b)中为“66”)的路径。边界评价部35预先记录成本最小的路径。

进而，边界评价部35在上述的最小值上加上式样的成本，决定为对象采样130的值。

再次参照图14。在步骤S23中，搜索成本C_total的最小值和此时的内腔边界和外膜边界的对，作为路径进行选择。此时，记录所选择的路径，并且更新路径成本。更新可按照以下的公式给出。

[公式8]

$C_{\mathrm{path}}\left(i_{\mathrm{line}}\right)\←-W_{\mathrm{pattern}}\·C_{\mathrm{pattern}}\left(i_{\mathrm{line}}\right)+\underset{i_{\mathrm{line}-1}}{\mathrm{MIN}}\left[C_{\mathrm{total}}\left(i_{\mathrm{line}-1}\right)\right]$

针对这些对象线的所有的内腔边界和外膜边界的对i_line实施步骤S22和S23的处理。边界评价部35将所得到的结果作为评价值送往边界决定部36。

针对所有线完成计算之后，在步骤S26中边界决定部36搜索路径成本C_path最小的内腔边界和外膜边界的对i_line。在此，路径成本C_path最小的对为ROI的右端线的边界。以后，对于左侧的线的边界，参照预先记录的路径来设定。

图17表示路径成本C_path最小的内腔边界和外膜边界的对i_line的例子。

以上是处理的流程。

再者，此后也可以与第一实施方式同样地实施线方向的平滑化。由此，能够得到平滑的边界线。

此外，对于基于动态规划法的路径搜索，也可以针对降低了分辨率(特别是深度方向)强度分布来实施以暂时决定候选之后，复原至原始的分辨率，在暂时决定的候选的周边(附近)的区域中进行搜索，以最终决定内腔边界和外膜边界。由此，能够维持检测性能的同时削减处理量。

此外，在基于动态规划法进行路径搜索时，也可以缩减为路径成本小的上位对、或者路径成本为规定值以下的对，来进行路径搜索。由此，能够维持检测性能的同时削减处理量。

最后，在上述说明中，以内腔边界和外膜边界为对进行路径搜索，但也可以单独进行路径搜索。

例如，图18是表示针对内腔边界及外膜边界分别进行路径检索的处理的步骤的流程图。在该处理中，首先检测内腔边界，然后缩减为比它深的位置处检测外膜边界。接着，按照相反的顺序检测外膜边界和比它浅的位置的内腔边界。然后，选择合计路径成本小的一方。

在步骤S31中，边界评价部35基于动态规划法暂时检测内腔边界，并基于检测结果记录路径成本A₁。在接下来的步骤S32中，边界评价部35以暂时检测出的内腔边界的下方为对象，基于动态规划法暂时检测外膜边界，并基于检测结果记录路径成本A₂。在步骤S33中，边界评价部35计算并记录按照内腔边界及外膜边界的顺序进行暂时检测时的路径成本A(＝A₁+A₂)。

接下来，在步骤S34中，边界评价部35基于动态规划法对外膜边界进行暂时检测，并基于检测结果记录路径成本B₁。在接下来的步骤S35中，边界评价部35以暂时检测出的外膜边界的上方为对象，基于动态规划法对内腔边界进行暂时检测，并基于检测结果记录路径成本B₂。在步骤S33中，边界评价部35计算并记录按照外膜边界及内腔边界的顺序进行暂时检测时的路径成本B(＝B₁+B₂)。

然后，边界决定部36利用成本A及成本B之中的、成本较小的一方的暂时检测结果，来决定外膜边界及内腔边界。

不过，在检测一个边界时，另一个边界尚未确定。因此，也可以如图11所示那样，使用固定一个边界通过改变暂时检测出的另一个边界以改变内中膜厚的多个模板，来评价合适度，将其评价值的最大值、或者平均值作为类似度。由此，能够维持检测性能的同时削减处理量。

以上，在第二实施方式中，由于考虑了线方向的边界连续性，因此即便在对应于内腔边界的强度较弱的情况下也能够检测稳定的边界。

在上述的第一及第二实施方式中，也可以进一步检测后壁，并将其检测结果用于边界的决定。如图1所示，后壁是存在于比外膜边界5深的位置处的区域。

图19是表示具有后壁检测部26的IMT测量部16a的结构的框图。在此，作为图7中记载的第一实施方式所涉及的IMT测量部16的变形例进行说明。

后壁检测部26位于强度分布生成部21的后级，利用由强度分布生成部21输出的强度分布检测后壁，将检测出的后壁范围输出至边界决定部24。

更加具体地进行说明。如图1所示，后壁是血管壁，在深度方向上观察，存在于比内腔边界4、内中膜3、外膜边界5更靠里的位置。并且，后壁作为高回波的区域而得到。因此，后壁检测部26根据由强度分布生成部21生成的强度分布来检测高回波的区域。其结果，边界决定部26能够缩减为检测出的后壁范围内，决定外膜边界5。

产业上的可利用性

本发明所涉及的IMT测量装置是精度良好的可自动地测量IMT的装置，能够缩短检查时间，例如对于动脉硬化的筛查很有用。

符号说明

11探头

12发送部

13接收部

14B图像生成部

15ROI设定部

16IMT测量部

17合成部

18显示部

21强度分布生成部

22模板生成部

23式样类似度计算部

24边界决定部

25IMT计算部

26后壁检测部

27对比度计算部

28检测窗

31强度分布生成部

32模板生成部

33式样类似度计算部

34边界连续度计算部

35边界评价部

36边界决定部

37IMT计算部

Claims (28)

1.一种超声波诊断装置，其向被检体发送超声波，并基于接收到的来自所述被检体的回波信号获取生物体信息，其中，所述超声波诊断装置具备：

强度分布生成部，其根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；

模板生成部，其基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；

式样类似度计算部，其是计算所述模板与所述强度分布之间的类似度的式样类似度计算部，即：使所述模板在所述深度方向上变化、且使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；和

边界决定部，其基于按每个所述变化所计算出的所述类似度，来决定内腔边界及外膜边界。

2.一种超声波诊断装置，其向被检体发送超声波，并基于接收到的来自所述被检体的回波信号获取生物体信息，其中，所述超声波诊断装置具备：

强度分布生成部，其根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；

模板生成部，其基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；

式样类似度计算部，其是计算所述模板与所述强度分布之间的类似度的式样类似度计算部，即：针对在所述深度方向的不同的位置处存在的边界的各候选的组合，使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；

边界连续度计算部，其针对关于相邻的多个声音线而确定的所述组合，基于各声音线的回波信号的强度值的差值来计算边界连续度；

边界评价部，其针对所述组合，生成综合所述式样类似度和所述边界连续度之后的评价值；和

边界决定部，其基于综合之后的所述评价值，来决定内腔边界和外膜边界。

3.根据权利要求1或2所述的超声波诊断装置，其中，

所述模板生成部所生成的模板至少包含3个区域，在所述至少3个区域中至少包含1个能够改变所述区域的范围的长度可变区域。

4.根据权利要求3所述的超声波诊断装置，其中，

所述至少3个区域按照顺序排列，

所述长度可变区域是所述至少3个区域之中的中央的区域。

5.根据权利要求3所述的超声波诊断装置，其中，

所述至少3个区域的一个端部的区域对应于低回波的组织区域，所述至少3个区域的另一个端部的区域对应于高回波的组织区域。

6.根据权利要求5所述的超声波诊断装置，其中，

所述低回波的组织区域对应于血流区域，所述高回波的组织区域对应于血管外膜壁及体组织。

7.根据权利要求4所述的超声波诊断装置，其中，

所述长度可变区域对应于血管的内中膜区域。

8.根据权利要求3所述的超声波诊断装置，其中，

所述式样类似度计算部计算使所述长度可变区域的长度标准化之后的类似度。

9.根据权利要求8所述的超声波诊断装置，其中，

所述式样类似度计算部将所述长度可变区域的平均强度用作所述长度可变区域的强度，由此使所述长度可变区域的长度标准化。

10.根据权利要求1或2所述的超声波诊断装置，其中，

所述模板生成部所生成的模板由2个种类的子模板构成，所述2个种类的子模板各自由至少2个区域构成，并能够改变2个种类的模板的距离，

所述式样类似度计算部利用所述2个种类的子模板分别计算类似度，并将合起且综合两者之后的值作为类似度。

11.根据权利要求10所述的超声波诊断装置，其中，

构成所述2个种类的模板的2个区域分别对应于低回波及高回波的组织区域。

12.根据权利要求10所述的超声波诊断装置，其中，

所述2个种类的模板的低回波及高回波的组织区域分别对应于血流区域及内中膜区域、或者分别对应于内腔边界及规定的组织区域，所述规定的组织区域是外膜边界、动脉外膜壁及体组织。

13.根据权利要求1或2所述的超声波诊断装置，其中，

所述模板生成部根据所述强度分布生成部生成的强度分布的对比度，来改变所生成的模板的系数值。

14.根据权利要求1或2所述的超声波诊断装置，其中，

所述式样类似度计算部所评价的类似度通过由所述模板生成部生成的模板而指定的边界间的强度差来定义。

15.根据权利要求1或2所述的超声波诊断装置，其中，

所述式样类似度计算部所评价的类似度通过所述模板生成部生成的模板与所述强度分布生成部生成的强度分布之间的标准化相关来定义。

16.根据权利要求1或2所述的超声波诊断装置，其中，

所述超声波诊断装置还具备后壁检测部，该后壁检测部根据所述强度分布生成部生成的强度分布检测高回波的区域，

所述边界决定部缩减为所检测出的后壁范围内来决定边界。

17.根据权利要求16所述的超声波诊断装置，其中，

所述后壁检测部所检测的高回波区域对应于血管壁。

18.根据权利要求2所述的超声波诊断装置，其中，

所述边界评价部利用动态规划法计算所述评价值。

19.根据权利要求18所述的超声波诊断装置，其中，

所述边界评价部利用降低了分辨率的强度分布来计算所述评价值，

所述边界决定部基于降低了分辨率的所述强度分布来暂时决定候选，然后在原分辨率的强度分布中从所述临时决定出的候选的周边的区域中决定内腔边界和外膜边界。

20.根据权利要求18所述的超声波诊断装置，其中，

所述边界评价部利用动态规划法缩减为所述评价值良好的上位候选、或者具有规定的评价值的候选，同时在声音线方向上计算所述评价值。

21.根据权利要求18所述的超声波诊断装置，其中，

所述边界连续度计算部计算由声音线方向的边界位置的差所定义的边界连续度。

22.根据权利要求18所述的超声波诊断装置，其中，

所述式样类似度计算部、所述边界连续度计算部、所述边界评价部、及所述边界决定部以内腔边界及外膜边界这两个边界作为对象，来评价所述可取得的边界候选，

所述边界决定部同时决定所述内腔边界及外膜边界。

23.根据权利要求18所述的超声波诊断装置，其中，

所述边界连续度计算部计算边界连续度，该边界连续度由在声音线方向的内腔边界及外膜边界的位置处确定的厚度的差来定义。

24.根据权利要求18所述的超声波诊断装置，其中，

所述边界决定部在基于确定的顺序来决定内腔边界和外膜边界时，

所述式样类似度计算部在另一个边界位置尚未确定的阶段计算类似度的情况下，使一个边界固定，针对具有与另一个边界位置对应的多个厚度的模板来计算类似度，将计算出的所述类似度的最大值、或者平均值作为所述一个边界的类似度。

25.根据权利要求18所述的超声波诊断装置，其中，

所述边界决定部，在检测出内腔边界之后缩减为所述内腔边界的下方并检测出外膜边界的情况下、和在检测出外膜边界之后缩减为所述外膜边界的上方并检测出内腔边界的情况下，比较在各个检测中用到的评价值，将评价值良好的检测边界作为最终的边界采用。

26.根据权利要求1至25任意一项所述的超声波诊断装置，其中，

所述超声波诊断装置还具备IMT计算部，该IMT计算部根据所述内腔边界和外膜边界计算内中膜的厚度。

27.一种测量内中膜的厚度的方法，是向被检体发送超声波，并基于接收到的来自所述被检体的回波信号测量血管壁的内中膜的厚度的方法，其中，该方法包括如下步骤：

根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；

基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；

是计算所述模板与所述强度分布之间的类似度的步骤，即：使所述模板在所述深度方向上变化、且使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；

基于按每个所述变化所计算出的所述类似度，来决定内腔边界及外膜边界；和

根据所决定的所述内腔边界及所述外膜边界，计算内中膜的厚度。

28.一种测量内中膜的厚度的方法，是向被检体发送超声波，并基于接收到的来自所述被检体的回波信号测量血管壁的内中膜的厚度的方法，其中，该方法包括如下步骤：

根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；

基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；

是计算所述模板与所述强度分布之间的类似度的步骤，即：针对在所述深度方向的不同的位置处存在的边界的各候选的组合，使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；

针对关于相邻的多个声音线而确定的所述组合，基于各声音线的回波信号的强度值的差值来计算边界连续度；

针对所述组合，生成综合所述式样类似度和所述边界连续度之后的评价值；

基于综合之后的所述评价值，来决定内腔边界和外膜边界；和

根据所决定出的所述内腔边界及所述外膜边界，计算内中膜的厚度。

Images