专利名称：超声波诊断装置及测量内中膜的厚度的方法
技术领域：
本发明涉及计算例如颈动脉等的血管的内中膜复合体的厚度的超声波诊断装置。
背景技术：
一般利用超声波诊断装置进行动脉硬化的诊断，该超声波诊断装置向被检体收发超声波，并基于来自所述被检体的回波信号获取生物体信息。周知在动脉硬化的诊断中，将颈动脉的内中膜复合体厚作为了解初始的粥样硬化的重要指标。“内中膜复合体”或者“内中膜”是存在于血管内腔与外膜之间的、血管壁的内膜及中膜的总称。并且，内中膜复合体厚或者内中膜厚也被称为“頂T” (Intima-Media Thickness)，在本申请说明书中记为“MT”。图I示出存在于血管内腔I和外膜2之间的内中膜3。在检查中，检测作为该血管内腔I与内膜的边界的内腔边界4、以及作为中膜与外膜2的边界的外膜边界5，并测量其厚度。作为自动测量该頂T的方法，例如已知专利文献I的技术。图2示出在专利文献I的技术中在确定外膜边界及内腔边界时所利用的回波强度分布。在专利文献I中公开了如下的方法，如图2所示那样，在深度方向的回波强度分布中，将强度变化在规定值以上的位置作为边界候选，将较大的一方定为外膜边界，将较小的一方定为内腔边界。现有技术文献专利文献专利文献I JP特开2008-168016号公报
发明内容
发明所要解决的技术问题但是，在专利文献I中，由于根据強度变化大或强度变化小统ー进行处理，因此有时无法测量出正确的MT。例如，图3表示内腔边界的强度变化大于外膜边界的强度变化的例子。是强度变化大的位置未必始终是外膜边界的例子。其结果，在专利文献I所公开的方法中，会将内腔边界作为外膜边界进行检测，故无法测量出正确的MT。本发明是鉴于上述问题而提出的，其目的在于正确地检测内腔边界和外膜边界，并实现正确的MT測量，而不依赖于内腔侧与外膜侧的回波强度变化的大小关系。用于解决技术问题的方法本发明的超声波诊断装置向被检体发送超声波，并基于接收到的来自所述被检体的回波信号获取生物体信息。本发明的某个所述超声波诊断装置具备強度分布生成部，其根据所述回波信号生成与血管壁相关的深度方向的回波強度分布；模板生成部，其基于预 先准备的表示回波強度分布的基准式样的模板，来生成用于边界检测的模板；式样类似度计算部，其是计算所述模板与所述强度分布之间的类似度的式样类似度计算部，即使所述模板在所述深度方向上变化、且使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；和边界决定部，其基于按每个所述变化所计算出的所述类似度，来决定内腔边界及外膜边界。本发明的另ー个超声波诊断装置具备強度分布生成部，其根据所述回波信号生成与血管壁相关的深度方向的回波強度分布；模板生成部，其基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；式样类似度计算部，其是计算所述模板与所述强度分布之间的类似度的式样类似度计算部，即针对在所述深度方向的不同的位置处存在的边界的各候选的组合，使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；边界连续度计算部，其针对关于相邻的多个声音线而确定的所述组合，基于各声音线的回波信号的强度值的差值来计算边界连续度；边界评价部，其针对所述组合，生成综合所述式样类似度和所述边界连续度之后的评价值；和边界决定部，其基于综合之后的所述评价值，来决定内腔边界和外膜边界。也可，所述模板生成部所生成的模板至少包含3个区域，在所述至少3个区域中至少包含I个能够改变所述区域的范围的长度可变区域。也可，所述至少3个区域按照顺序排列，所述长度可变区域是所述至少3个区域之中的中央的区域。也可，所述至少3个区域的一个端部的区域对应于低回波的组织区域，所述至少3个区域的另ー个端部的区域对应于高回波的组织区域。也可，所述低回波的组织区域对应于血流区域，所述高回波的组织区域对应于血管外膜壁及体组织。也可，所述长度可变区域对应于血管的内中膜区域。也可，所述式样类似度计算部计算使所述长度可变区域的长度标准化之后的类似度。也可，所述式样类似度计算部将所述长度可变区域的平均強度用作所述长度可变区域的強度，由此使所述长度可变区域的长度标准化。
也可，所述模板生成部所生成的模板由2个种类的子模板构成，所述2个种类的子模板各自由至少2个区域构成，井能够改变2个种类的模板的距离，所述式样类似度计算部利用所述2个种类的子模板分别计算类似度，并将合起且综合两者之后的值作为类似度。也可，构成所述2个种类的模板的2个区域分别对应于低回波及高回波的组织区域。也可，所述2个种类的模板的低回波及高回波的组织区域分别对应于血流区域及内中膜区域、或者分别对应于内腔边界及规定的组织区域，所述规定的组织区域是外膜边界、动脉外膜壁及体组织。也可，所述模板生成部根据所述強度分布生成部生成的強度分布的对比度，来改变所生成的模板的系数值。也可，所述式样类似度计算部所评价的类似度通过由所述模板生成部生成的模板而指定的边界间的强度差来定义。也可，所述式样类似度计算部所评价的类似度通过所述模板生成部生成的模板与所述强度分布生成部生成的強度分布之间的标准化相关来定义。也可，所述超声波诊断装置还具备后壁检测部，该后壁检测部根据所述強度分布生成部生成的強度分布检测高回波的区域，所述边界决定部缩减为所检测出的后壁范围内来决定边界。也可，所述后壁检测部所检测的高回波区域对应于血管壁。也可，所述边界评价部利用动态规划法计算所述评价值。也可，所述边界评价部利用降低了分辨率的強度分布来计算所述评价值，所述边界决定部基于降低了分辨率的所述强度分布来暂时决定候选，然后在原分辨率的強度分布中从所述临时决定出的候选的周边的区域中决定内腔边界和外膜边界。也可，所述边界评价部利用动态规划法缩减为所述评价值良好的上位候选、或者具有规定的评价值的候选，同时在声音线方向上计算所述评价值。也可，所述边界连续度计算部计算由声音线方向的边界位置的差所定义的边界连续度。也可，所述式样类似度计算部、所述边界连续度计算部、所述边界评价部、及所述边界决定部以内腔边界及外膜边界这两个边界作为对象，来评价所述可取得的边界候选，所述边界决定部同时决定所述内腔边界及外膜边界。也可，所述边界连续度计算部计算边界连续度，该边界连续度由在声音线方向的内腔边界及外膜边界的位置处确定的厚度的差来定义。也可，所述边界决定部在基于确定的顺序来决定内腔边界和外膜边界时，所述式样类似度计算部在另ー个边界位置尚未确定的阶段计算类似度的情况下，使一个边界固定，针对具有与另ー个边界位置对应的多个厚度的模板来计算类似度，将计算出的所述类似度的最大值、或者平均值作为所述ー个边界的类似度。也可，所述边界决定部，在检测出内腔边界之后缩减为所述内腔边界的下方并检测出外膜边界的情况下、和在检测出外膜边界之后缩减为所述外膜边界的上方并检测出内腔边界的情况下，比较在各个检测中用到的评价值，将评价值良好的检测边界作为最終的边界采用。也可，所述超声波诊断装置还具备MT计算部，该MT计算部根据所述内腔边界和外膜边界计算内中膜的厚度。本发明的方法是向被检体发送超声波，并基于接收到的来自所述被检体的回波信号測量血管壁的内中膜的厚度的方法，其中，该方法包括如下步骤根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；基于预先准备的表示回波强度分布的基准式样的模板，来生成用于边界检测的模板；是计算所述模板与所述强度分布之间的类似度的步骤，即使所述模板在所述深度方向上变化、且使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；基于按每个所述变化所计算出的所述类似度，来决定内腔边界及外膜边界；和根据所決定的所述内腔边界及所述外膜边界，计算内中膜的厚度。本发明的方法是向被检体发送超声波，并基于接收到的来自所述被检体的回波信号測量血管壁的内中膜的厚度的方法，其中，该方法包括如下步骤根据所述回波信号生成与血管壁相关的深度方向的回波强度分布；基于预先准备的表示回波强度分布的基准式样、的模板，来生成用于边界检测的模板；是计算所述模板与所述强度分布之间的类似度的步骤，即针对在所述深度方向的不同的位置处存在的边界的各候选的组合，使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；针对关于相邻的多个声音线而确定的所述组合，基于各声音线的回波信号的强度值的差值来计算边界连续度；针对所述组合，生成综合所述式样类似度和所述边界连续度之后的评价值；基于综合之后的所述评价值，来决定内腔边界和外膜边界；和根据所决定出的所述内腔边界及所述外膜边界，计算内中膜的厚度。发明效果根据本发明，并不是独立地检测内腔边界和外膜边界，而是将2个边界作为组来进行检测，即便在内腔边界的强度变化大于外膜边界的强度变化的情况下，也能够正确地測量IMT0
图I是表示存在于血管内腔I与外膜2之间的内中膜3的图。图2是表示在专利文献I的技术中在确定外膜边界及内腔边界时所利用的回波强度分布的图。图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是表示路径成本Cpath最小的内腔边界和外膜边界的对iliM的例子的图。图18是表示针对内腔边界及外膜边界分别进行路径检索的处理的步骤的流程图。
图19是表示具有后壁检测部26的MT测量部16a的结构的框图。
具体实施例方式以下，參照附图，说明本发明的超声波诊断装置的实施方式。(第一实施方式)
图4是本发明的第一实施方式所涉及的超声波诊断装置10的结构框图。
超声波诊断装置10具备探头11、发送部12、接收部13、B图像生成部14、R0I设定部15、IMT測量部16、合成部17、和显示部18。
发送部12生成超声波脉冲，并输出至探头11。探头11将由发送部12输出的超声波脉冲发送至生物体内，并接收从生物体内反射回来的回波信号。探头11将接收到的回波信号输出至接收部13。图5表示用于发送超声波脉冲的、内置于超声波探头11的超声波振动元件组6。超声波振动元件组6 (以下记为“振动器6”。)构成例如各超声波振动元件沿着ー个方向排列的、所谓的ID阵列振动器。振动器6例如由压电体构成，通过驱动压电体来发送超声波，此外，压电体接收超声波并将超声波变换为电信号。振动器6使各超声波振动元件依次驱动来进行超声波的收发，能够扫描规定范围。此外，振动器6使来自多个超声波振动元件的各超声波的相位在规定的位置(焦点位置)处重合，也能够接收在焦点位置处反射的信号。图6表示后者的例子。图6(a)及(b)示意地表示使用沿着X方向排列的多个超声波振动元件形成焦点时的超声波的集束波。超声波的集束波具有图示那样的规定宽度，在z轴方向的规定的深度处具有焦点。其中，该集束波的中心轴如图6(c)及(d)中所示那样也被称为“声音线”。再次參照图4。接收部13使从探头11输出的回波信号成束后，输出至B图像生成部14和MT測量部16。B图像生成部14根据回波信号生成断面图像，将其断面图像的数据输出至合成部17。ROI设定部15从用户接受确定ROI位置的指示，基于该指示在动脉壁的位置处设定关心区域(Region of Interest R0I)。MT测量部16测量MT。更为具体而言，頂T测量部16根据由接收部13输出的回波信号和由ROI设定部15输出的ROI位置的信息，检测内腔边界和外膜边界。进而，IMT測量部16基于检测出的内腔边界及外膜边界的位置来计算MT值。MT測量部16将确定边界位置的信息(边界位置信息)输出至合成部17，将IMT值输出至显示部18。合成部17基于由B图像生成部14输出的B图像的数据和由MT測量部16输出的边界位置信息，合成IMT值和B图像。更为具体而言，合成部17将内腔边界及外膜边界重叠在B图像上，并输出得到的合成图像数据。显示部18将由MT测量部16输出的MT值和由合成部17输出的合成结果显示在监视器上。显示部18将从合成部17输出的合成图像数据显示在例如液晶显示器等的显示设备上。
接下来，详细说明MT測量部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和MT计算部25。IMT计算部25基于所決定的内腔边界及外膜边界来计算MT值，并将计算出的IMT值输出至显示部17。接下来，说明本实施方式的MT计算处理的步骤。图8是表示IMT计算处理的步骤的流程图。首先，在步骤Sll中，強度分布生成部21对所接收的回波信号进行包络线检波、对数压缩，生成回波強度分布。该处理针对ROI内的各点进行。接下来，在步骤S12中，強度分布生成部21为了降低噪声，针对强度分布应用例如3X3的平均值滤波器，进行平滑化。接着，处理转移至式样类似度计算部23。再者，以下的处理按照每个声音线(以线为单位)来进行。首先，在步骤S13中，式样类似度计算部23评价ROI中的强度分布和表示内中膜式样的模板之间的类似度。模板是由模板生成部22生成的。作为模板，使用图9所示的系数式样(内中膜式样)。在图9中，范围I1对应于血管内腔，范围I2和范围I3对应于内中膜，范围I4对应于外膜。图9所示的模板是微分型滤波器，作为基准的系数式样预先在深度方向上交替设置取负值的し、I3的范围、和取正值的12、I4的范围。如图示那样，在模板中至少设置三个区域、即对应于血管内腔的区域(或者式样)、对应于内中膜的区域、及对应于外膜的区域。式样类似度计算部23将该模板作为检测窗，如图10所示那样，ー边在深度方向上移动ー边类似度，搜索合适的深度。此外，式样类似度计算部23还评价相对于内中膜厚的类似度。S卩、如图11所示那样，在模板生成部22中制作将内中膜厚改变成多祥的多个模板，并计算这些的模板与強度分布之间的类似度，捜索合适的厚度。以下，进行更为详细的说明。
图10表示作为对象的声音线(line)的強度分布90、和改变了深度的模板组91。模板组91的各模板是与声音线平行的深度方向的I列式样。在图中示出了深度pl、p2、p3、p4、……、ρ10的模板。再者，图示的各模板仅深度不同，其式样是相同的。式样类似度计算部23首先计算作为对象的声音线(line)的強度分布90、和深度Pl且内中膜厚q时的模板之间的类似度。类似度由模板T与强度分布X之间的内积值给出，在将深度P的强度分布设为Xp，将模板的内中膜厚设为q，将对应于内中膜厚q的模板设为T,，将模板尺寸设为L时，类似度Y通过以下的公式计算。
[公式I]Y{p,q) = YjXpiD ■ Tq{i)
/=1在图10所示的例子中，在最合适的深度p4处，模板中包含的内中膜厚q的式样与对象线的強度分布90中的内中膜厚的式样相一致。但是，这仅是便于理解的一例。在实际中，在最合适的深度处，模板的内中膜厚q的式样与对象线的强度分布90中的内中膜厚的式样未必一致。因此，在本实施方式中，接下来，式样类似度计算部23在深度pi处ー边改变内中膜厚q的式样，一边计算类似度。图11表示作为对象的声音线(line)的強度分布90和改变了厚度的模板组92。式样类似度计算部23 —边改变内中膜厚q，一边求得与各个厚度对应的类似度Y。接着，式样类似度计算部23将深度变更至p2，一边改变内中膜厚q —边求得类似度Y。这种的处理持续至深度plO。在图8的步骤S14中，边界决定部24搜索类似度Y (P、q)最大的深度p_max和内中膜厚q_max，在步骤中S15中，分别决定内腔边界Ii和外膜边界ma。即，若将模板的血管内腔长度设为11，则由以下的公式给出。[公式2]Ii = p_max+llma = li+q_max以上是以线为单位的处理。在上述的处理中，系数值&1’、&2’、&3’、&4’也可以设定为所对应的范围长lpl2、l3、I4的倒数值，以免类似度缩减为特定的内中膜厚。即，在与范围长Ip 12、13、I4对应的各区域中，计算平均強度，成为加权的形式。此时的类似度由以下的公式计算。[公式3]ド(Aが:-f全ん(/)トf.Sw,)-笋1^(ハ + デ·; Σ'(,·)
,1 /=1*2 /W] i l3 /-Z1--Z2114 / /ド/2マ/3]1在此，是权重系数。再有，在改变内中膜厚吋，既可以使I2固定而使I3可变，也可以使12、I3基于规定的比例都可变。
步骤S15是切换至下ー线的切換处理。以下，针对各线进行上述的处理。以线为単位的处理结束之后，在步骤S16中，在线方向上对边界位置进行平滑化。以上是第一实施方式的动作流程。再者，在上述说明中，将类似度设定为強度分布与取得微分型滤波器形式的模板之间的内积值，但也可以是相关值。此时，在模板的系数值中不需要包含负值。此外，如图12所示那样，类似度也可以设定为准备分别与内腔边界和外膜边界匹配的子模板，分别计算类似度，通过加权相加等而组合之后的值。在各子模板中，规定了与边界及除此之外的组织等对应的式样，但这些的式样是可变的。即，与边界对应的式样的位置是可变的。上述的2种模板中的、与边界区域及除此之外的组织等对应的区域，分别对应于各高回波及低回波的组织区域。 此外，模板的系数值也可以基于ROI中的描绘度进行变更，例如，可以计算强度值的方差，在方差大的情况下，增强模板的对比度，相反在方差小的情况下，减弱模板的对比度。最后，在上述说明中，假定直接检测内腔边界和外膜边界，但也可以首先使用微分型滤波器等检测后壁，缩减为后壁内部来检测内腔边界和外膜边界。这在血管中出现多重回波的情况下可有效地发挥作用。以上，在第一实施方式中，确定与对象线的強度分布最类似的模板，以检测内腔边界和外膜边界。模板规定了包含内腔边界及外膜边界的強度分布的式样。这意味着并不是独立地检测内腔边界和外膜边界，而是在考虑内腔边界的强度变化及外膜边界的强度变化的同时来检测内腔边界及外膜边界。由此，在内腔边界的強度大于外膜边界的情况下，也不会发生误检测。(第二实施方式)本实施方式与第一实施方式的不同点在于，本实施方式中的MT测量部的处理。具体而言，本实施方式中的IMT測量部在线方向上进行平滑化之前，除了考虑式样类似度以外，还考虑线方向的连续性来选择边界。如果进行该处理，则特别是即便在对应于内腔边界的强度不明确的情况下，也能够检测稳定的边界。本实施方式的超声波诊断装置的结构及动作除了 MT測量部以外，与图4所示的超声波诊断装置10相同。因此，作为本实施方式的超声波诊断装置的结构援引图4，并且省略IMT测量部以外的说明。图13是表示本实施方式的IMT测量部30的结构的框图。MT測量部30具有強度分布生成部31、模板生成部32、式样类似度计算部33、边界连续度计算部34、边界评价部35、边界决定部36、和MT计算部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基于所決定的内腔边界及外膜边界计算MT值，并将计算出的MT值输出至显示部17。该处理与图7所示的MT计算部25的处理相同。接下来，说明本实施方式的边界检测处理的步骤。下述的处理主要由边界连续度计算部34、边界评价部35及边界决定部36来进行。图14是表示基于动态规划法的边界检测处理的步骤的流程图。在本实施方式中，按照搜索由式样类似性和线方向的边界的连续性所定义的成本最小的线方向的路径(=边界)的形式来实施边界检测处理。在步骤S21中，边界评价部35针对所有的内腔边界和外膜边界的对iliM使路径成本Cpath初始化。初始值由与内中膜式样之间的类似度Cpattwn给出。若将权重设为WpattOTn，则路径成本Cpath由以下的公式给出。
[公式4] Cpath (iline)^pattern Cpattern (土line)在此，类似度Cpatten是与第一实施方式同样的值。但是，主要想关注的是在线内进行标准化之后的值。不过，不需要在所有的线内进行标准化。例如，可以至少仅对与内中膜对应的线内的值进行标准化。所谓“内腔边界和外膜边界的对”相当于内腔边界和外膜边界的组合。若将ROI的采样点设为N点时，则对的数量为N的平方。不过，由于内腔边界处在比外膜边界浅的位置、或者系统规定的最大測量厚等的约束条件，实际计算中使用的对数少于N的平方。接下来，进入至第2线。然后，在步骤S22中，式样类似度计算部33计算与针对前ー线的各内腔边界和外膜边界的对Ut1的内中膜式样之间的类似度。此外，边界连续度计算部34针对前ー线的各内腔边界和外膜边界的对Ut1，分别计算内中膜厚和边界位置的连续度。边界评价部35接受这些计算结果，评价线方向的内中膜厚和边界位置的连续性。与内中膜式样之间的类似度Cpatten,与第一实施方式同样既可以由强度分布和模板的内积值给出，也可以是相关值。
若将对象线和前ー线的内腔边界位置分别设为LIlim、LIlil^1，将外膜边界设为MAline, MAn1，则内中膜厚的连续度Cmdiness由以下的公式给出。[公式5]Cthickness (iine-i) — I (MAline-Llline) - (MAline_1-Llline_1)此外，边界位置的连续度Cdistanra由以下的公式给出。[公式6]Cdis tan ce (iiine-ι) — I Llline_Llline_11 +1MA1 Jne-MA1边界评价部35在线内分别对评价值Cttiekness和Cdistanee标准化之后，如以下的公式那样将它们统一成候选路径的成本Ct()tal。 [公式7] Ctotal 、Iline—Iノ レpath (i
line-1)+'^thickness しthickness (line-1)+Wdis tan ce ^dis tan ce \line-1)在此，权重Wthidiness及Wdistanra是相对于各评价值的权重系数。再者，公式7中包含评价值cThi.ss和Cdista_，但只要包含边界位置的连续度Cdistanre即可，可以没有内中膜厚的连续度
^Thickness°边界评价部35针对预先決定的路径对象的内腔边界和外膜边界的对Ut1实施以上的计算。例如在图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中，搜索成本Cttrtal的最小值和此时的内腔边界和外膜边界的对，作为路径进行选择。此时，记录所选择的路径，并且更新路径成本。更新可按照以下的公式给出。[公式8] ^path ^ line ) < ^pattern ' ^pattern^Hne^tota!^hine~\^针对这些对象线的所有的内腔边界和外膜边界的对iline实施步骤S22和S23的处理。边界评价部35将所得到的结果作为评价值送往边界决定部36。针对所有线完成计算之后，在步骤S26中边界决定部36捜索路径成本Cpath最小的内腔边界和外膜边界的对られ。在此，路径成本Cpath最小的对为ROI的右端线的边界。以后，对于左侧的线的边界，參照预先记录的路径来设定。图17表示路径成本Cpath最小的内腔边界和外膜边界的对iliM的例子。
以上是处理的流程。再者，此后也可以与第一实施方式同样地实施线方向的平滑化。由此，能够得到平滑的边界线。此外，对于基于动态规划法的路径搜索，也可以针对降低了分辨率(特别是深度方向)強度分布来实施以暂时决定候选之后，复原至原始的分辨率，在暂时決定的候选的周边(附近)的区域中进行搜索，以最终决定内腔边界和外膜边界。由此，能够维持检测性能的同时削减处理量。此外，在基于动态规划法进行路径搜索时，也可以缩减为路径成本小的上位对、或者路径成本为规定值以下的对，来进行路径搜索。由 此，能够维持检测性能的同时削减处理量。最后，在上述说明中，以内腔边界和外膜边界为对进行路径搜索，但也可以单独进行路径搜索。例如，图18是表示针对内腔边界及外膜边界分别进行路径检索的处理的步骤的流程图。在该处理中，首先检测内腔边界，然后缩减为比它深的位置处检测外膜边界。接着，按照相反的顺序检测外膜边界和比它浅的位置的内腔边界。然后，选择合计路径成本小的一方。在步骤S31中，边界评价部35基于动态规划法暂时检测内腔边界，并基于检测结果记录路径成本ん。在接下来的步骤S32中，边界评价部35以暂时检测出的内腔边界的下方为对象，基于动态规划法暂时检测外膜边界，并基于检测结果记录路径成本A2。在步骤S33中，边界评价部35计算并记录按照内腔边界及外膜边界的顺序进行暂时检测时的路径成本 A ( = A^A2)。接下来，在步骤S34中，边界评价部35基于动态规划法对外膜边界进行暂时检測，并基于检测结果记录路径成本B”在接下来的步骤S35中，边界评价部35以暂时检测出的外膜边界的上方为对象，基于动态规划法对内腔边界进行暂时检測，并基于检测结果记录路径成本B2。在步骤S33中，边界评价部35计算并记录按照外膜边界及内腔边界的顺序进行暂时检测时的路径成本B ( = B^B2)。然后，边界决定部36利用成本A及成本B之中的、成本较小的一方的暂时检测结果，来决定外膜边界及内腔边界。不过，在检测ー个边界时，另ー个边界尚未确定。因此，也可以如图11所示那样，使用固定一个边界通过改变暂时检测出的另ー个边界以改变内中膜厚的多个模板，来评价合适度，将其评价值的最大值、或者平均值作为类似度。由此，能够维持检测性能的同时削
减处理量。以上，在第二实施方式中，由于考虑了线方向的边界连续性，因此即便在对应于内腔边界的强度较弱的情况下也能够检测稳定的边界。在上述的第一及第ニ实施方式中，也可以进ー步检测后壁，并将其检测结果用于边界的決定。如图I所示，后壁是存在于比外膜边界5深的位置处的区域。图19是表示具有后壁检测部26的MT测量部16a的结构的框图。在此，作为图7中记载的第一实施方式所涉及的IMT測量部16的变形例进行说明。后壁检测部26位于强度分布生成部21的后级，利用由強度分布生成部21输出的強度分布检测后壁，将检测出的后壁范围输出至边界决定部24。更加具体地进行说明。如图I所示，后壁是血管壁，在深度方向上观察，存在于比内腔边界4、内中膜3、外膜边界5更靠里的位置。并且，后壁作为高回波的区域而得到。因此，后壁检测部26根据由强度分布生成部21生成的強度分布来检测高回波的区域。其结果，边界决定部26能够缩减为检测出的后壁范围内，决定外膜边界5。产业上的可利用性本发明所涉及的MT測量装置是精度良好的可自动地測量MT的装置，能够缩短检查时间，例如对于动脉硬化的筛查很有用。符号说明11 探头12发送部13接收部14B图像生成部15R0I 设定部16MT 测量部17合成部18显示部21强度分布生成部22模板生成部23式样类似度计算部24边界决定部25MT 计算部26后壁检测部27对比度计算部28检测窗31强度分布生成部32模板生成部33式样类似度计算部34边界连续度计算部35边界评价部36边界决定部 37MT 计算部
权利要求
1.一种超声波诊断装置，其向被检体发送超声波，并基于接收到的来自所述被检体的回波信号获取生物体信息，其中，所述超声波诊断装置具备 強度分布生成部，其根据所述回波信号生成与血管壁相关的深度方向的回波強度分布； 模板生成部，其基于预先准备的表示回波強度分布的基准式样的模板，来生成用于边界检测的模板； 式样类似度计算部，其是计算所述模板与所述强度分布之间的类似度的式样类似度计算部，即使所述模板在所述深度方向上变化、且使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；和 边界决定部，其基于按每个所述变化所计算出的所述类似度，来决定内腔边界及外膜边界。
2.一种超声波诊断装置，其向被检体发送超声波，并基于接收到的来自所述被检体的回波信号获取生物体信息，其中，所述超声波诊断装置具备 強度分布生成部，其根据所述回波信号生成与血管壁相关的深度方向的回波強度分布； 模板生成部，其基于预先准备的表示回波強度分布的基准式样的模板，来生成用于边界检测的模板； 式样类似度计算部，其是计算所述模板与所述强度分布之间的类似度的式样类似度计算部，即针对在所述深度方向的不同的位置处存在的边界的各候选的组合，使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度； 边界连续度计算部，其针对关于相邻的多个声音线而确定的所述组合，基于各声音线的回波信号的强度值的差值来计算边界连续度； 边界评价部，其针对所述组合，生成综合所述式样类似度和所述边界连续度之后的评价值；和 边界决定部，其基于综合之后的所述评价值，来决定内腔边界和外膜边界。
3.根据权利要求I或2所述的超声波诊断装置，其中， 所述模板生成部所生成的模板至少包含3个区域，在所述至少3个区域中至少包含I个能够改变所述区域的范围的长度可变区域。
4.根据权利要求3所述的超声波诊断装置，其中， 所述至少3个区域按照顺序排列， 所述长度可变区域是所述至少3个区域之中的中央的区域。
5.根据权利要求3所述的超声波诊断装置，其中， 所述至少3个区域的一个端部的区域对应于低回波的组织区域，所述至少3个区域的另ー个端部的区域对应于高回波的组织区域。
6.根据权利要求5所述的超声波诊断装置，其中， 所述低回波的组织区域对应于血流区域，所述高回波的组织区域对应于血管外膜壁及体组织。
7.根据权利要求4所述的超声波诊断装置，其中，所述长度可变区域对应于血管的内中膜区域。
8.根据权利要求3所述的超声波诊断装置，其中， 所述式样类似度计算部计算使所述长度可变区域的长度标准化之后的类似度。
9.根据权利要求8所述的超声波诊断装置，其中， 所述式样类似度计算部将所述长度可变区域的平均強度用作所述长度可变区域的强度，由此使所述长度可变区域的长度标准化。
10.根据权利要求I或2所述的超声波诊断装置，其中， 所述模板生成部所生成的模板由2个种类的子模板构成，所述2个种类的子模板各自由至少2个区域构成，并能够改变2个种类的模板的距离， 所述式样类似度计算部利用所述2个种类的子模板分别计算类似度，并将合起且综合两者之后的值作为类似度。
11.根据权利要求10所述的超声波诊断装置，其中， 构成所述2个种类的模板的2个区域分别对应于低回波及高回波的组织区域。
12.根据权利要求10所述的超声波诊断装置，其中， 所述2个种类的模板的低回波及高回波的组织区域分别对应于血流区域及内中膜区域、或者分别对应于内腔边界及规定的组织区域，所述规定的组织区域是外膜边界、动脉外膜壁及体组织。
13.根据权利要求I或2所述的超声波诊断装置，其中， 所述模板生成部根据所述強度分布生成部生成的強度分布的对比度，来改变所生成的模板的系数值。
14.根据权利要求I或2所述的超声波诊断装置，其中， 所述式样类似度计算部所评价的类似度通过由所述模板生成部生成的模板而指定的边界间的强度差来定义。
15.根据权利要求I或2所述的超声波诊断装置，其中， 所述式样类似度计算部所评价的类似度通过所述模板生成部生成的模板与所述强度分布生成部生成的强度分布之间的标准化相关来定义。
16.根据权利要求I或2所述的超声波诊断装置，其中， 所述超声波诊断装置还具备后壁检测部，该后壁检测部根据所述強度分布生成部生成的强度分布检测高回波的区域， 所述边界决定部缩减为所检测出的后壁范围内来决定边界。
17.根据权利要求16所述的超声波诊断装置，其中， 所述后壁检测部所检测的高回波区域对应于血管壁。
18.根据权利要求2所述的超声波诊断装置，其中， 所述边界评价部利用动态规划法计算所述评价值。
19.根据权利要求18所述的超声波诊断装置，其中， 所述边界评价部利用降低了分辨率的強度分布来计算所述评价值， 所述边界决定部基于降低了分辨率的所述强度分布来暂时决定候选，然后在原分辨率的強度分布中从所述临时决定出的候选的周边的区域中决定内腔边界和外膜边界。
20.根据权利要求18所述的超声波诊断装置，其中，所述边界评价部利用动态规划法缩减为所述评价值良好的上位候选、或者具有规定的评价值的候选，同时在声音线方向上计算所述评价值。
21.根据权利要求18所述的超声波诊断装置，其中， 所述边界连续度计算部计算由声音线方向的边界位置的差所定义的边界连续度。
22.根据权利要求18所述的超声波诊断装置，其中， 所述式样类似度计算部、所述边界连续度计算部、所述边界评价部、及所述边界决定部以内腔边界及外膜边界这两个边界作为对象，来评价所述可取得的边界候选， 所述边界决定部同时决定所述内腔边界及外膜边界。
23.根据权利要求18所述的超声波诊断装置，其中， 所述边界连续度计算部计算边界连续度，该边界连续度由在声音线方向的内腔边界及外膜边界的位置处确定的厚度的差来定义。
24.根据权利要求18所述的超声波诊断装置，其中， 所述边界决定部在基于确定的顺序来决定内腔边界和外膜边界吋， 所述式样类似度计算部在另ー个边界位置尚未确定的阶段计算类似度的情况下，使一个边界固定，针对具有与另ー个边界位置对应的多个厚度的模板来计算类似度，将计算出的所述类似度的最大值、或者平均值作为所述ー个边界的类似度。
25.根据权利要求18所述的超声波诊断装置，其中， 所述边界决定部，在检测出内腔边界之后缩减为所述内腔边界的下方并检测出外膜边界的情况下、和在检测出外膜边界之后缩减为所述外膜边界的上方并检测出内腔边界的情况下，比较在各个检测中用到的评价值，将评价值良好的检测边界作为最終的边界采用。
26.根据权利要求I至25任意一项所述的超声波诊断装置，其中， 所述超声波诊断装置还具备頂T计算部，该MT计算部根据所述内腔边界和外膜边界计算内中膜的厚度。
27.一种测量内中膜的厚度的方法，是向被检体发送超声波，并基于接收到的来自所述被检体的回波信号測量血管壁的内中膜的厚度的方法，其中，该方法包括如下步骤 根据所述回波信号生成与血管壁相关的深度方向的回波強度分布； 基于预先准备的表示回波強度分布的基准式样的模板，来生成用于边界检测的模板；是计算所述模板与所述强度分布之间的类似度的步骤，即使所述模板在所述深度方向上变化、且使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度； 基于按每个所述变化所计算出的所述类似度，来决定内腔边界及外膜边界；和 根据所決定的所述内腔边界及所述外膜边界，计算内中膜的厚度。
28.一种测量内中膜的厚度的方法，是向被检体发送超声波，并基于接收到的来自所述被检体的回波信号測量血管壁的内中膜的厚度的方法，其中，该方法包括如下步骤 根据所述回波信号生成与血管壁相关的深度方向的回波強度分布； 基于预先准备的表示回波強度分布的基准式样的模板，来生成用于边界检测的模板；是计算所述模板与所述强度分布之间的类似度的步骤，即针对在所述深度方向的不同的位置处存在的边界的各候选的组合，使与所述模板的内中膜厚对应的式样变化，从而按每个变化来计算所述模板与所述强度分布之间的类似度；针对关于相邻的多个声音线而确定的所述组合，基于各声音线的回波信号的强度值的差值来计算边界连续度； 针对所述组合，生成综合所述式样类似度和所述边界连续度之后的评价值； 基于综合之后的所述评价值，来决定内腔边界和外膜边界；和 根据所决定出的所述内腔边界及所述外膜边界，计算内中膜的厚度。
全文摘要
本发明提供一种超声波诊断装置及测量内中膜的厚度的方法。超声波诊断装置向被检体发送超声波，并基于接收到的来自被检体的回波信号获取生物体信息。超声波诊断装置具备强度分布生成部，根据回波信号生成与血管壁相关的深度方向的回波强度分布；模板生成部，基于预先准备的表示回波强度分布的基准式样的模板，生成用于边界检测的模板；式样类似度计算部，计算模板与强度分布之间的类似度，即使模板在深度方向上变化、且使与模板的内中膜厚对应的式样变化，从而按每个变化来计算模板与强度分布之间的类似度；和边界决定部，基于按每个变化计算出的类似度来决定内腔边界及外膜边界。在内中膜的厚度测量的内腔边界与外膜边界的检测中，不依赖于内腔侧和外膜侧的回波强度变化的大小关系。
文档编号A61B8/06GK102695457SQ20108006024
公开日2012年9月26日 申请日期2010年11月12日 优先权日2010年2月10日
发明者川端章裕, 福元刚智, 铃木隆夫, 高木一也 申请人:松下电器产业株式会社