发明内容
技术问题
如所述,在空间合成处理之后,趋于在子区域之间产生增益差异。由于这种增益差异,具有在图像中作为障碍物的边界线趋于在超声波图像上的子区域之间出现的问题。特别地,该问题在空间合成处理被应用于谐波分量时变得特别显著。
本发明的优点在于提高通过空间合成而产生的超声波图像的图像质量。本发明的另一个优点在于不产生在合成的帧数据中的子区域之间的增益差异,或尽管产生此增益差异,也减小增益差异。本发明的又一个优点还在于通过简单的结构来消除或减小增益差异。
问题解决方案
根据本发明的一个方案,提供一种超声波诊断装置,包括:发送和接收单元,其在切换波束偏转条件的同时顺序地形成n个波束扫描平面(其中n是大于或等于2的整数),以产生n个帧数据组;以及合成单元,其通过所述n个帧数据组的合成处理来产生合成的帧数据,其中当所述n个帧数据组空间上重叠时,基于重叠结构的差异来定义多个子区域,所述多个子区域包括中心子区域和多个外围子区域,在所述中心子区域中合成全部的所述n个帧数据组,在所述外围子区域中合成所述n个帧数据组的部分,并且,在各所述外围子区域的合成处理中,应用对在外围子区域中的重叠结构相对于在中心子区域中的重叠结构的不足进行补偿的条件,使得在所述中心子区域中的合成处理之后的增益和在各所述外围子区域中的合成处理之后的增益相匹配。
根据上文所述的结构,由所述发送和接收单元来形成n个波束扫描平面。所述n个波束扫描平面在各自不同的波束偏转条件下形成。所述合成单元使用相应于所述n个波束扫描平面的所述n个帧数据组来产生合成的帧数据。在二维合成空间中,基于所述n个帧数据组的重叠方式的不同定义了多个子区域;也就是,重叠结构。在此情况下,理论上,将产生具有数量为n的重叠的子区域到具有数量为1的重叠的子区域。然而,在显示区域外的所述子区域或部分可能从处理目标中忽略。由于实际重叠结构在所述多个子区域中不是均匀的,当在各子区域中执行简单的加性平均处理时,将导致增益差异。因此,在上文所述的结构中,各子区域呈现出具有相同的重叠结构。期望地,重叠结构的不足被在合成处理中进行补偿,使得在各外围子区域(部分帧数据组重叠的子区域)的合成处理中产生与所述中心子区域(全部帧数据组重叠的子区域)的重叠结构相似的重叠结构。在此情况中,例如,可以重复地参照覆盖所述外围子区域的帧数据,以增加在合成处理中待使用的数据的呈现。根据此方法,可以在所述多个子区域之中利用用于合成处理的共有计算方程。可替换地,代替重复地参照相同帧数据,数据可以乘以相应于参照物的数量的系数以获得相似的结果。
如所述的,上文所述的方法集中于子区域之中的实际重叠结构的差异,并且应用必要的补偿(在计算上操作),使得重叠结构呈现为在用于合成处理的子区域之中相同。换句话说,在所述多个子区域之间匹配合成处理的条件(期望地,加性平均条件)。通过该结构,消除或减小了在所述子区域之间的增益差异。此结构使得呈现在超声波图像上的并且构成障碍的边界线消失或减少,并且实现了整个超声波图像的均匀的感觉或亮度。
对于不足的补偿,理想地,使用具有与丢失帧数据相同的增益关系的帧数据(具有关于正负号对称的偏转角)。可替换地,可以使用具有相近增益关系(具有与关于正负号对称的偏转角相近的偏转角)的帧数据。期望地,参照所述中心子区域的合成处理条件来匹配所述外围子区域的合成处理条件。可替换地,可以将特定的外围子区域的合成处理条件设定为基准,并且可以将包括所述中心子区域的其它子区域的合成处理条件与该基准相匹配。
可替换地,除如上文所述的增益差异补偿之外,可以同时应用如专利文献1所述的加权处理。根据该组合,当增益差异即使在如上文所述的增益差异补偿之后仍然保持时,能够设定边界线不突出。
根据本发明的另一个方案,优选地,对所述外围子区域中的重叠结构的不足进行补偿的条件是以复制的方式来参照覆盖所述外围子区域并且具有与不覆盖所述外围子区域的帧数据相同的增益关系的帧数据的条件。通常,在关于正负号对称的波束偏转角下获得的帧数据对之间满足相同的增益关系。期望地,判定多个波束偏转条件以使得满足一个或多个这样的关系。
根据本发明的另一个方案,优选地,数字n是大于或等于3的奇数,所述n个波束扫描平面包括由0度的波束偏转角形成的波束扫描平面,以及由关于正负号对称的波束偏转角形成的至少一个波束扫描平面对,并且相应于所述一对波束扫描平面的帧数据对具有相同的增益关系。
根据本发明的另一个方案,优选地,所述合成单元包括:合成函数,其为所述多个子区域之间的共有函数,并且其在分母中具有值n并且在分子中具有n个加性元素;以及对应表,其为其中为各子区域记录n个帧数据标识符的表,所述n个帧数据标识符用于赋给待代入所述n个加性元素的n个回波值,并且其中以复制的方式记录相同帧数据标识符,从而对各外围子区域的重叠结构的不足进行补偿,并且所述合成单元在参照各子区域的所述对应表的同时将n个回波值代入所述合成函数,并且执行所述合成函数。
根据上文所述的结构,利用多个子区域之间的共有合成函数。也就是说,不再需要为各子区域准备专用合成函数。所述合成函数是期望的应用加性平均处理的函数,并且分母为与加项的数量相应的值n。分子由n个加性元素形成。将n个加性元素相加,且实际值被代入加性元素。对应表指示该代换。在表中,以复制的方式记录用于某一外围子区域的相同帧数据标识符,使得能够应用上文所述的补偿。即使在波束偏转角(绝对值)改变时,只要保持正负对称,表的内容也不需要改变。当如专利文件1所述的加权处理将被结合地应用时,可以将加权系数应用到分子的个别加性元素上,或可以在回波值在合成函数中被代换之前对所述回波值应用加权处理。
根据本发明的另一个实施例,优选地,所述合成单元进一步包括内插单元,其基于n个帧数据的全部或部分通过内插处理来产生待被代入n个加性元素的n个回波值的全部或部分。当在所述帧数据之间的采样布置(坐标系)不同时,通过基于外围数据的内插处理来产生所需的数据。
根据本发明的另一个方案,优选地,用在活体中产生的谐波回波分量来形成所述n个帧数据组。由于谐波分量与基波分量相比是显著衰减的,趋于产生由于传播距离等的差异导致的增益差异,并且因此,上述方法的应用是可期望的。
可替换地,所述多个帧数据组的合成处理能够通过诸如计算机的信息处理装置执行。在此情况中,由超声波系统提供的帧数据阵列在线或离线地被输入到信息处理装置。
具体实施方式
将会参照附图详细描述本发明的优选的实施例。
图1显示根据本发明的优选的实施例的超声波诊断装置。图1是示意地显示超声波诊断装置的整体结构的框图。该超声波诊断装置用于医学领域,并且其是通过向活体发送超声波以及从活体接收超声波而形成超声波图像的装置。
在图1中,阵列换能器10被设置于探针(未显示)中,并且通过阵列换能器10来发送和接收超声波。在图1的实例结构中,阵列换能器10包括以线性形状布置的多个换能器元件10a。也就是,在本实施例中,利用了线性探针。可替换地,可以利用凸起型探针等。
开口12被设定在阵列换能器10上。例如,开口12是发送和接收开口。通过使用开口12发送超声波,形成发送波束,并且使用开口12来接收反射波。基于通过接收产生的多个接收信号,通过将在下文描述的整相求和处理(phase-alignment and summing process)来形成接收波束。附图标记14表示实例的超声波波束,并且该超声波波束14与发送/接收总波束(发送/接收联合波束)相关。超声波波束14形成在与阵列换能器10的布置方向垂直的方向上。也就是说,超声波波束14的波束偏转角是0度。通过沿布置方向扫描开口12;也就是,通过执行电动线性扫描,超声波波束14被线性地扫描,并且形成扫描平面16。扫描平面16与二维数据捕捉区域相对应。扫描平面16具有长方形形状。
在本实施例中,为了执行空间合成方法,扫描平面20和24在扫描平面16之后顺序地形成。也就是说,三个扫描平面16,20和24以循环的方式产生。更具体地,使用开口12,在正侧的+α度的波速偏转角的方向上形成超声波波束18,并且超声波波束18被电扫描的同时保持波束偏转角,从而形成扫描平面20。扫描平面20可以被认为是正侧偏转扫描平面。相似地,使用开口20设定在负侧的-α度的波速偏转角,并且在该方向上形成超声波波束22。超声波波束22被电扫描的同时保持波束偏转角,从而形成扫描平面24。扫描平面24可以被认为是负侧偏转扫描平面。扫描平面20和24具有相互成反向定向的平行四边形形状。
在本实施例中,被设定为用于形成扫描平面20的+α度的波束偏转角以及被设定为用于形成扫描平面24的-α度的波束偏转角是相同的角度(按照绝对值)。也就是,扫描平面20和扫描平面24形成一对。通过在正侧和负侧均设定相同的波束偏转角,能够匹配增益。相似地,在使用凸起型探针的情况中,可以应用与上文所述相似的发送和接收序列。
发送单元26为发送波束成形器,并且在发送期间提供平行于多个换能器元件的多个发送信号。借助该结构,超声波被发送到活体的内部。在活体中的各反射点处被反射的回波(反射波)被阵列换能器10接收,并且多个接收信号从形成接收开口的多个换能器元件平行地输出,并且被输入到接收单元28。接收单元28为接收波束成形器,并且对多个接收的信号执行整相求和处理,以在整相求和之后输出作为接收信号的波束数据。
为一个扫描平面输出一个帧数据组(此后简单地提及为“帧”)。一帧由布置在扫描方向上的多个波束数据组形成,并且每个波束数据组由布置在深度方向上的多个回波值(回波数据)形成。实际上,按时序次序发送回波数据阵列。
信号处理器30为对波束数据执行各种信号处理的模块。可替换地,信号处理器30可以基本上通过软件功能来实现。信号处理器30一般包括检波器、增益控制器、对数压缩器、帧相关装置等。在本实施例中,信号处理器30进一步包括谐波分量提取器32和合成单元34。谐波分量提取器32为用于提取在活体中产生的谐波分量的模块。例如,可以使用诸如带通滤波、脉冲反相(pulse inversion)、脉冲调制等各种方法来提取包括在接收信号中的谐波分量。例如,当应用脉冲反相时,执行两次发送同时对于各波束地址反转相位。
合成单元34为用于执行空间合成处理的模块。具体地,在本实施例中,在以时序次序输入的帧阵列中,执行用3帧作为单元的合成帧的处理,以产生合成帧。合成单元34输出合成的帧阵列(合成的帧数据阵列)。稍后将会参照图2等详细描述合成单元34的具体功能。
数字扫描转换器(DSC)36为基于输入帧形成B模式图像(二维黑白断层图像)的模块。DSC36具有插值函数、坐标变换函数等。产生的图像数据通过显示处理器(未显示)被发送到显示单元38。在显示单元38中,B模式图像作为视频图像或作为静止图像被显示。在本实施例中,能够显示出应用了空间合成的B模式图像;特别是,在帧内的增益差异被补偿的B模式图像。由于此,和相关技术相比,B模式图像的图像质量能够被提升。当执行如上文所述的谐波分量的提取时,谐波图像作为B模式图像被显示。在此情况中,在活体中的组织中产生的谐波或在活体中的造影剂中产生的谐波分量被显示在屏幕上。
控制器40执行如图1所示的结构的操作控制。输入单元42被连接到控制器40。输入单元42由操作面板形成,并且该操作面板具有诸如键盘和轨迹球的各种输入装置。例如,在图1的实例结构中,电影存储器(cine memory)被设置在信号处理器30和DSC36之间。该电影处理器具有存储帧阵列的环形缓冲区的结构。可选地,这种电影处理器可以设置在DSC36的下游。此外,合成单元34可以设置在DSC36的下游。在此情况下,对显示帧阵列应用空间合成处理。在这种情况下,同样优选地,应用将在下文描述的增益补偿条件。
图2显示通过图1所示的合成单元执行的合成处理的第一配置。附图标记44显示输入到合成单元的帧阵列(帧数据阵列)。帧阵列44由布置在时间轴方向上的多个帧构成,并且更具体地,一般由3种类型的帧A、B和C构成。帧A相应于图1所示的扫描平面16;也就是,帧A为在0度的波束偏转角的条件下获得的帧。帧B是相应于图1所示的扫描平面20的帧;也就是,帧B为在正侧的+α度的波束偏转角的条件下获得的帧。帧C是相应于图1所示的扫描平面24的帧;也就是,帧C为在负侧的-α度的波束偏转角的条件下获得的帧。以循环的方式获得这3种类型的帧。
在图2中,附图标记46显示移动平均处理。在移动平均处理46中,在作为单元的帧阵列44的三个帧上重复地执行加性平均处理。在图2中,用括号显示加性平均处理的各单元。以此方式,空间合成处理;也就是,联合处理,在3帧的单元上顺序地执行。
由于该处理,产生合成的帧阵列(合成的帧数据阵列)48。合成的帧阵列48由以时序次序布置的多个合成的帧50构成。通过使帧A、B和C在二维合成空间上重叠而产生各合成的帧50中的每一个,并且在右侧和左侧伸出的部分被切除。也就是说,这些部分不是合成处理的目标。由此,在合成帧50中的区域(整个区域)基于三个帧的重叠的方式的差异,也就是,重叠结构的差异,而被分为4个区域(子区域)。更具体地,如图2所示,定义了区域I、区域II、区域III和区域IV。
在此,区域I为A、B和C三帧重叠的区域,区域II为A、B两帧重叠的区域,区域III为A、C两帧重叠的区域,并且区域IV为仅被帧A覆盖的区域。也就是说,在区域IV中重叠的数量是1。以此方式,基于3个帧的形式的差异,重叠的数量以及重叠结构在合成帧50的局部区域中不同,并且由此,定义了4种区域。实际上,仅显示了探针附近的区域。该区域在图2中显示为显示区域52。实际上,原则上不显示区域IV。然而,显示区域52在深度方向上的尺寸;也就是,测量范围,能够变化并且能够通过用户设定。
如所述的,由于重叠结构或重叠的数量取决于局部区域而不同,如果简单地执行加性平均处理而不考虑该差异,将会导致在区域之间的增益差异的上述问题。将会参照图3和图4描述该问题和解决该问题的方法。
图3显示第一比较例。图3是显示各区域的合成处理内容的图。之后将描述的图4、图7和图8也是显示相似内容的图。
在图3中,附图标记54显示一区域。附图标记56显示应用于各区域的计算方程。附图标记58显示将被赋给包括在计算方程中的加性元素X、Y和Z的输入数据。更具体地,显示了待赋给X、Y和Z的回波值所属的帧类型的标识符A、B和C。例如,在区域I中,将由A、B和C标识的3个回波值代入X、Y和Z,这些值的和除以相加元素的数量(也就是3),并且合成的回波数据值作为加性平均值而被产生。合成的回波数据值形成合成帧中的一个数据值,并且相似地,在区域II中,计算两个回波值的平均值,并且在区域III中,计算两个回波值的平均值。在此,应该注意的是,在这些计算中使用的计算方程的分母是2。也就是说,对于形成在与阵列换能器正交的方向上的波束,在自此波束倾斜的波束中,增益相对减小,并且增益在帧A、B和C之间不匹配。更具体地,相比于帧A,对于帧B和C而言增益较低。在该环境下,如果如图3所示的用于各区域的计算方程被应用,则相比于区域I,对于区域II和区域III而言,在计算结果中增益变得较高。也就是说,在区域II和区域III中,由于具有更低增益的回波值的比率的分子增大,引起了上文所述的增益差异的问题。附带地,区域IV具有最大的增益。
图4显示消除如上文所述的增益差异的第一配置。在该第一结构中,用非常简单的结构消除了增益差异。更具体地,如下文中将描述的,使用区域共有的计算方程。现在将详细描述该配置。
与上文相似,附图标记54显示一区域,并且附图标记56显示计算方程。在第一配置中,如上文所述,共有计算方程在多个区域中被使用。也就是说,仅使用单个计算方程。该计算方程是D=(X+Y+Z)/3。在此,D是合成回波值。
由于在分子中的3个加性元素的代换中应用了对在重叠结构中的差异进行补偿的特殊操作,因此可以使用该单个计算方程。具体地,参照显示为附图标记58的输入数据,3种类型的帧A、B和C与区域I中的加性元素X、Y和Z相关,在图3中显示的比较例也与此相同。然而,在第一配置中,在区域II和III中,相同帧标识符被赋给加性元素Y和Z;也就是参照相同的帧两次。在计算方程中,对于分子,由相同的帧获得的相同的回波值被代换两次。
在更进一步的具体说明中,与区域I相比,在区域II中,缺失来自帧C的贡献,并且帧B被以复制的方式被参照从而补偿该缺失的部分。相似地,在区域III中,缺失来自帧B的贡献,并且帧C被参照两次从而补偿不足。由此,在全部区域中,参照具有基本增益的一帧以及具有稍小于基本帧的增益的两帧,并且能够消除区域之间的增益差异。在本实施例中,在帧B和帧C中,应用了在0度的波束偏转角周围的关于正负号对称的波束偏转角,并且由此,帧B和帧C具有对称关系或成对关系,并且因此,能够利用上文所述的复制的参照。
可替换地,可以操纵计算方程以将用于调整的系数赋给用于区域II或区域III的分子,并且以消除增益差异。可替换地,区域II和区域III的增益可以被设定为标准增益,并且区域I的增益可以被减小;更具体地,Y/2可以被代换为Y,并且Z/2可以被代换为Z,从而消除增益差异。
根据上文所述的第一配置,单个计算方程能够不需进一步操纵而被使用,通过由附图标记60显示的表来操纵待赋给计算方程的回波值,使得能够获得的优点在于,能够用简单的结构来实现增益差异补偿。换句话说,由于增益差异能够通过仅记录表60的内容而被补偿,因此结构非常简单。特别地,当图1中显示的合成单元以软件功能实现时,优选地使用单一计算方程和表的联合。
在第一配置中,帧B和帧C满足相同偏转角的条件,但是并不需要是这种情况。更具体地,当帧对不具有完全相同的偏转角、但是满足接近角条件时,能够通过应用上文所述的方法来减小增益差异。
图5显示为图1中显示的合成单元的实例结构的框图。合成单元34实际上被以软件功能实现,并且图5以框图显示该功能以易于理解。
在图5中,显示了在实例结构中内插处理器66和区域判定器68被赋予地址数据。地址数据显示当前输入的回波值的属性;也就是,地址数据具有诸如回波值所属的帧的类型、帧中的坐标等的信息。另一方面,输入回波值被输入到内插处理器66以及帧存储器62。也就是说,帧阵列被输入到内插处理器66以及帧存储器62。帧存储器62为存储先于当前帧一帧的帧的存储器,并且从帧存储器62输出的帧被存储在帧存储器64中。帧存储器64为存储先于当前帧两帧的帧的存储器。当前帧被直接地输入到内插处理器66。
区域判定器68参照当前输入的回波值(所关注的回波值)的地址数据,以判定回波值所属的区域。也就是,区域判定器68判定所关注的回波值属于区域I、区域II、区域III和区域IV中的哪一个。该判定结果被输出到表70。该表70为具有通过图4中显示的附图标记60标识的内容的表。更具体地,表70具有用于标识3个帧的信息,该3个帧用于产生将被赋给计算方程的分子的3个加性元素的3个回波值。从表70中输出用于执行相应于该所关注的回波值的函数计算所需的信息。该信息标识了将被代入加性元素的回波值所属的帧。在实施例中,内插处理器66为参照帧A的地址空间来执行内插处理以将来自帧B和帧C的回波值赋给帧A的坐标的模块。也就是说,在本实施例中,在帧B和帧C上应用内插处理,以提取回波值作为插值。这是因为,在3个类型的帧之间坐标系相互不同。可替换地,可以对帧A应用内插处理。以任一方法,当由附图标记74显示的计算方程在函数计算器72处被执行时,所需的回波值从表70的输出被赋给函数计算器72,并且所需的回波值从内插处理器66被赋给函数计算器72。例如,对于区域I,从帧A、B和C产生的3个回波值被代入计算方程的分子。在此情况下,从帧B和帧C中产生的回波值为插值。在区域II和区域III中,从帧A中提取的回波值,以及作为从帧B或帧C中产生的插值的回波值被代入计算方程的分子。在此情况下,如图4所示,根据加性元素Y和Z的需要,以复制的方式赋予相同的回波值。通过该处理,使用共有计算方程产生合成的回波值,同时补偿回不足的波值。
下面,将会参照图6至图8描述合成处理的第二配置。
图6是第二配置的概念图。附图标记76显示被输入到合成单元的帧阵列,并且在该实例结构中,以循环的方式输入5种类型的帧。在图6中以A,B,C,D和E标识5种类型的帧。应用到帧的波束偏转条件彼此不同;也就是,用5个波束偏转角产生该5种类型的帧。帧A是基本帧;也就是,在0度的波束偏转角条件下产生的帧。帧B和帧E是以相同偏转角但是不同符号而产生的帧,并且形成帧对。相似地,帧C和帧D是在正负号对称的条件下产生的帧对。
附图标记78显示移动平均处理的内容。在本实施例中,用排列在时间轴上作为一单元的5个帧顺序地执行加性平均处理。附图标记80显示合成帧阵列,其由多个合成的帧82构成。通过重叠5种类型的帧,也就是通过合成处理,产生各合成的帧82。基于重叠状态的差异,也就是重叠结构的差异,定义了在合成的帧82上的9种类型的区域。因此,呈现了从区域I到区域IX的9个区域。在各区域中,重叠结构不同于其它区域,并且同时,重叠数量不同于其它区域。
图7显示了第二比较例。在图6中显示的5种类型的帧的合成被执行的情况中,如果在图7显示的条件下执行合成处理,如上文所述,将在区域之间引起增益差异。也就是说,各区域之间的计算方程不同;特别地,各区域之间的分母不同,并且必须为区域准备各自不同的计算方程。另外,由于在各区域中执行简单加性平均处理,其问题在于根据帧类型的增益差异也将出现在合成结果中。
图8显示解决上述问题的第二配置。如图8所示,为9个区域设定共有计算方程,并且相应于区域I(也就是基本区域或中心区域)的重叠数量5显示于分母中。计算方程的分子由5个加性元素构成;也就是,X、Y、Z、V和W。对于各加性元素,输入数据相关,如附图标记58所示。在附图中,A-E显示帧类型,也就是帧标识符。从标识的帧参照回波值,并且将回波值代入相应的加性元素。
例如,在区域II中,帧C被参照两次,并且通过复制的记录来实现缺失的帧的补偿。这类似地应用于区域III,其中帧D被以复制的方式记录。在区域IV中,参照区域I的重叠结构,帧D和帧E被各自参照两次从而匹配区域I的重叠结构。而且,在其它区域中,执行用于匹配区域I的重叠结构的复制的记录,并且通过此处理在区域间匹配增益。
附带地,在区域IX中,帧A与5个加性元素相关,并且增益将变高。然而,区域IX是最深的区域,并且是实际不显示的部分。因此,该配置不会引起问题。
如所述的,还在第二配置中,管理赋给加性元素的回波值的属性使得产生与基本区域相似的重叠结构。也就是说,应用需要的复制记录。借助此处理,能够消除或减小区域之间的增益差异,并且因此,超声波图像的质量能够被显著提高。
在上文所述的实施例中,当需要更加模糊区域之间的边界时,优选地用在边界周围的平滑加权变化来实施加权相加。在此情况下,在图4和图8中显示的计算方程中,各加性元素能够被乘以加权系数,并且加权系数可以变化,使得加权处理与如上文所述的增益补偿同时实现。
在上述实施例中,用包括参照帧的3个帧组或5个帧组作为一单元实施合成处理,但是可替换地,可以用偶数数量的帧作为单元实施合成处理。另外,虽然可期望的是扫描平面的类型的数量和用于移动平均处理的帧的数量一致,但是这些数量可以各自不同。