CN104039235B - 超声波诊断装置 - Google Patents

超声波诊断装置 Download PDF

Info

Publication number
CN104039235B
CN104039235B CN201280066935.9A CN201280066935A CN104039235B CN 104039235 B CN104039235 B CN 104039235B CN 201280066935 A CN201280066935 A CN 201280066935A CN 104039235 B CN104039235 B CN 104039235B
Authority
CN
China
Prior art keywords
frame
subregion
synthesis
frame data
sub area
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201280066935.9A
Other languages
English (en)
Other versions
CN104039235A (zh
Inventor
植田充
宇野隆也
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Healthcare Corp
Original Assignee
Hitachi Aloka Medical Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Aloka Medical Ltd filed Critical Hitachi Aloka Medical Ltd
Publication of CN104039235A publication Critical patent/CN104039235A/zh
Application granted granted Critical
Publication of CN104039235B publication Critical patent/CN104039235B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5246Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
    • A61B8/5253Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode combining overlapping images, e.g. spatial compounding
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8995Combining images from different aspect angles, e.g. spatial compounding
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54Control of the diagnostic device

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

在超声波诊断装置中,在具有由三种类型的帧A、B和C所获得的帧组作为单元的帧序列上连续地执行移动平均处理,从而连续地产生合成的帧。每个合成的帧具有中心子区域和多个外围子区域。为了在产生各合成的帧时匹配多个子区域之间的增益,在个别外围子区域应用特殊的合成条件。具体地,执行将相同数据值冗余地代换合成计算方程中的多个元素的处理。使用该空间合成,能够消除或减小在多个子区域之间的增益差异。

Description

超声波诊断装置
技术领域
本发明涉及一种超声波诊断装置,并且尤其涉及多个帧数据组的合成处理;也就是空间合成处理。
背景技术
超声波诊断装置正被用于医学领域。超声波诊断装置发送超声波到活体并且接收来自活体的超声波,并且基于由超声波的发送和接收所获得的信息来形成超声波图像。例如,当使用超声波诊断装置用B模式图像(二维黑白断层图像)观测颈动脉、乳腺等时,使用电动线性扫描探针。在电动线性扫描中,用于形成波束的发送和接收开口被设定在探针中的阵列换能器上,并且扫描发送和接收开口。通过这种处理,超声波波束被线性地扫描,并且形成波束扫描平面。通常,超声波波束形成在垂直于阵列换能器的方向上;也就是,0度的波束偏转角的方向上。当超声波波束的线性扫描在这种条件下执行时,将形成长方形波束扫描平面。当利用凸形探针时,在0度的波束偏转角的条件下,将会形成扇形的波束扫描平面。
与相关技术的超声波诊断相比,已知一种空间合成的方法。在空间合成中,多个波束扫描平面被重复地形成,同时以循环的方式切换波束偏转条件,并且因此获得以时序次序布置的帧数据阵列。在帧数据阵列中,随后用预定数量的帧作为单元来顺序地执行合成处理,以顺序地生成合成帧数据。根据空间合成,能够基于通过在关于某组织的多个方向上执行超声波的发送和接收而获得的数据来形成图像,产生组织的优等图像。另外,这种方法能够减少诸如在空间中呈现的伪影的与信号相关的噪声。
在上述空间合成中,在二维合成空间中,产生了多个子区域;也就是,基于多个帧数据组的重叠方式的不同而限定了多个子区域;也就是,重叠结构的不同。例如,当在0度的波束偏转角的条件下获得的长方形形状的基本帧数据组、在+15度的波束偏转角的条件下获得的平行四边形形状的正侧偏转帧数据组、以及在-15度的波束偏转角的条件下获得的平行四边形形状的负侧偏转帧数据组被合成时,产生三个帧数据组重叠的中心子区域(I)、两个帧数据组重叠的两外围子区域(II)和(III)、以及仅基本帧数据组存在的深入的子区域(IV)。深入的子区域(IV)是通常不被显示的区域。
基本上,空间合成为在多个帧之间的各区域中的坐标执行加性平均处理。在上述实例结构中,在中心子区域(I)中,执行加性平均计算以将由三个帧数据组获得的三个回波值之和除以3,并且在两外围子区域(II)和(III)中,执行加性平均计算以将由两个帧数据组获得的两个回波值之和除以2。
然而,当比较超声波波束没有被偏转的情况(0度的波束偏转角的情况)和超声波波束被偏转的情况时,通常,后者的增益低于前者的增益。这是因为,即使在阵列换能器上的物理接收开口的尺寸不变,接收开口的有效尺寸随波束偏转角的增加而减小,并且,相似地,由于波束偏转角增加,使个别换能器元件的有效面积减小。从阵列换能器到在某一深度处的观察点的超声波传播距离还取决于波束偏转角而改变,并且因此,除非发送条件和接收条件根据波束偏转角精确地改变,否则很难匹配在具有不同偏转条件的多组数据之中的增益。
在上文所述的环境下,趋于在合成处理之后的合成帧数据中的子区域之间产生增益差异,并且因而,具有在子区域之间的边界突出的问题。专利文献1公开了一种应用加权加性处理的技术,其中权值在子区域之间的边界周围平滑地改变。虽然该处理本身是有效的,但是没有根本地消除在子区域之间的增益差异。专利文献2还公开了一种加权加性处理。专利文献2还公开了校正信号振幅的变化。然而,为了根据波束偏转角精确地调整增益,将会需要非常复杂的控制。
相关技术参考
专利文献
[专利文献1]JP2003-61955A
[专利文献2]JP2002-526229A
发明内容
技术问题
如所述,在空间合成处理之后,趋于在子区域之间产生增益差异。由于这种增益差异,具有在图像中作为障碍物的边界线趋于在超声波图像上的子区域之间出现的问题。特别地,该问题在空间合成处理被应用于谐波分量时变得特别显著。
本发明的优点在于提高通过空间合成而产生的超声波图像的图像质量。本发明的另一个优点在于不产生在合成的帧数据中的子区域之间的增益差异,或尽管产生此增益差异,也减小增益差异。本发明的又一个优点还在于通过简单的结构来消除或减小增益差异。
问题解决方案
根据本发明的一个方案,提供一种超声波诊断装置,包括:发送和接收单元,其在切换波束偏转条件的同时顺序地形成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是显示根据本发明的优选的实施例的超声波诊断装置的框图。
图2是显示合成处理的第一实例的概念图。
图3是显示第一比较例的图,并且其显示各区域的合成处理的内容。
图4是显示第一实例的图,并且其显示各区域的合成处理的内容。
图5是显示图1中所示的合成单元的实例结构的图。
图6是显示合成处理的第二实例的概念图。
图7是显示第二比较例的图,并且其显示各区域的合成处理的内容。
图8是显示第二实例的图,并且其显示各区域的合成处理的内容。
具体实施方式
将会参照附图详细描述本发明的优选的实施例。
图1显示根据本发明的优选的实施例的超声波诊断装置。图1是示意地显示超声波诊断装置的整体结构的框图。该超声波诊断装置用于医学领域,并且其是通过向活体发送超声波以及从活体接收超声波而形成超声波图像的装置。
在图1中,阵列换能器10被设置于探针(未显示)中,并且通过阵列换能器10来发送和接收超声波。在图1的实例结构中,阵列换能器10包括以线性形状布置的多个换能器元件10a。也就是,在本实施例中,利用了线性探针。可替换地,可以利用凸起型探针等。
开口12被设定在阵列换能器10上。例如,开口12是发送和接收开口。通过使用开口12发送超声波,形成发送波束,并且使用开口12来接收反射波。基于通过接收产生的多个接收信号,通过将在下文描述的整相求和处理(phase-alignmentandsummingprocess)来形成接收波束。附图标记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为用于提取在活体中产生的谐波分量的模块。例如,可以使用诸如带通滤波、脉冲反相(pulseinversion)、脉冲调制等各种方法来提取包括在接收信号中的谐波分量。例如,当应用脉冲反相时,执行两次发送同时对于各波束地址反转相位。
合成单元34为用于执行空间合成处理的模块。具体地,在本实施例中,在以时序次序输入的帧阵列中,执行用3帧作为单元的合成帧的处理,以产生合成帧。合成单元34输出合成的帧阵列(合成的帧数据阵列)。稍后将会参照图2等详细描述合成单元34的具体功能。
数字扫描转换器(DSC)36为基于输入帧形成B模式图像(二维黑白断层图像)的模块。DSC36具有插值函数、坐标变换函数等。产生的图像数据通过显示处理器(未显示)被发送到显示单元38。在显示单元38中,B模式图像作为视频图像或作为静止图像被显示。在本实施例中,能够显示出应用了空间合成的B模式图像;特别是,在帧内的增益差异被补偿的B模式图像。由于此,和相关技术相比,B模式图像的图像质量能够被提升。当执行如上文所述的谐波分量的提取时,谐波图像作为B模式图像被显示。在此情况中,在活体中的组织中产生的谐波或在活体中的造影剂中产生的谐波分量被显示在屏幕上。
控制器40执行如图1所示的结构的操作控制。输入单元42被连接到控制器40。输入单元42由操作面板形成,并且该操作面板具有诸如键盘和轨迹球的各种输入装置。例如,在图1的实例结构中,电影存储器(cinememory)被设置在信号处理器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个帧组作为一单元实施合成处理,但是可替换地,可以用偶数数量的帧作为单元实施合成处理。另外,虽然可期望的是扫描平面的类型的数量和用于移动平均处理的帧的数量一致,但是这些数量可以各自不同。

Claims (6)

1.一种超声波诊断装置,包括:
发送和接收单元,其在切换波束偏转条件的同时顺序地形成n个波束扫描平面,其中n是大于或等于2的整数,以产生n个帧数据组;以及
合成单元,其通过所述n个帧数据组的合成处理来产生合成的帧数据,其中
当所述n个帧数据组空间上重叠时,基于重叠结构的差异来定义多个子区域,
所述多个子区域包括中心子区域和多个外围子区域,在所述中心子区域中合成全部的所述n个帧数据组,在所述外围子区域中合成所述n个帧数据组的部分,
在各所述外围子区域的合成处理中,应用对在所述外围子区域中的重叠结构相对于在所述中心子区域中的重叠结构的不足进行补偿的条件,使得在所述中心子区域中的合成处理之后的增益和在各所述外围子区域中的合成处理之后的增益相匹配,并且
在各所述外围子区域的合成处理中,通过覆盖所述外围子区域的相同帧数据的复制使用来对在外围子区域中的重叠结构相对于在所述中心子区域中的重叠结构的不足进行补偿,其中,所述相同帧数据具有关于正负号对称的偏转角或具有与关于正负号对称的偏转角相近的偏转角。
2.根据权利要求1所述的超声波诊断装置,其中
在各所述外围子区域的合成处理中,覆盖所述外围子区域的所述相同帧数据被参照两次。
3.根据权利要求1所述的超声波诊断装置,其中
数字n是大于或等于3的奇数,
所述n个波束扫描平面包括由0度的波束偏转角形成的波束扫描平面,以及由关于正负号对称的波束偏转角形成的至少一个波束扫描平面对,并且
相应于所述波束扫描平面对的帧数据对具有相同的增益关系。
4.根据权利要求1所述的超声波诊断装置,其中
所述合成单元包括:
合成函数,其为所述多个子区域之间的共有函数,并且其在分母中具有值n并且在分子中具有n个加性元素;以及
对应表,其为对于各子区域记录n个帧数据标识符的表,所述n个帧数据标识符用于赋给待代换所述n个加性元素的n个回波值,并且其中以复制的方式记录相同帧数据识别符,从而对各外围子区域的重叠结构的不足进行补偿,并且
所述合成单元在参照各子区域的所述对应表的同时将n个回波值代入合成函数,并且执行所述合成函数。
5.根据权利要求4所述的超声波诊断装置,其中
所述合成单元进一步包括内插单元,其通过基于所述n个帧数据组的全部或部分的内插处理来产生待代换所述n个加性元素的所述n个回波值的全部或部分。
6.根据权利要求1所述的超声波诊断装置,其中
所述n个帧数据组由在活体中产生的谐波回波分量形成。
CN201280066935.9A 2012-01-11 2012-11-15 超声波诊断装置 Active CN104039235B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2012-002796 2012-01-11
JP2012002796A JP5272084B2 (ja) 2012-01-11 2012-01-11 超音波診断装置
PCT/JP2012/079637 WO2013105339A1 (ja) 2012-01-11 2012-11-15 超音波診断装置

Publications (2)

Publication Number Publication Date
CN104039235A CN104039235A (zh) 2014-09-10
CN104039235B true CN104039235B (zh) 2015-11-25

Family

ID=48781303

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280066935.9A Active CN104039235B (zh) 2012-01-11 2012-11-15 超声波诊断装置

Country Status (4)

Country Link
US (1) US20140358003A1 (zh)
JP (1) JP5272084B2 (zh)
CN (1) CN104039235B (zh)
WO (1) WO2013105339A1 (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101555259B1 (ko) 2013-11-29 2015-09-24 알피니언메디칼시스템 주식회사 초음파 영상 합성 방법 및 장치
JP5919311B2 (ja) 2014-01-16 2016-05-18 富士フイルム株式会社 超音波診断装置および超音波画像生成方法
CN104287777B (zh) * 2014-10-17 2017-09-26 苏州佳世达电通有限公司 超音波扫描方法及超音波扫描系统
CN104546008B (zh) * 2015-02-02 2016-09-21 声泰特(成都)科技有限公司 一种基波/谐波融合与空间复合相结合的成像方法
EP3547923B1 (en) 2016-11-29 2021-01-06 Koninklijke Philips N.V. Ultrasound imaging system and method
JP6729798B2 (ja) * 2017-04-13 2020-07-22 株式会社島津製作所 X線撮影装置
EP3513735B1 (en) * 2017-04-25 2022-05-25 Sogang University Research Foundation Device and method for generating ultrasound vector doppler image using plane wave synthesis
KR102655278B1 (ko) * 2018-07-24 2024-04-08 삼성메디슨 주식회사 초음파 영상 장치 및 그 표시 방법
CN109350115A (zh) * 2018-12-04 2019-02-19 珠海医凯电子科技有限公司 交错扫描双波束方法
JP2022037259A (ja) * 2018-12-28 2022-03-09 株式会社Lily MedTech 超音波撮像システム及び超音波撮像方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6224552B1 (en) * 1998-10-01 2001-05-01 Atl Ultrasound Ultrasonic diagnostic imaging system with reduced spatial compounding seam artifacts
CN101199430A (zh) * 2006-12-15 2008-06-18 深圳迈瑞生物医疗电子股份有限公司 空间复合成像方法、设备及其超声成像系统
CN101866480A (zh) * 2009-04-15 2010-10-20 深圳迈瑞生物医疗电子股份有限公司 一种空间复合成像中的帧相关处理方法及系统

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5908389A (en) * 1996-09-27 1999-06-01 Atl Ultrasound, Inc. Ultrasonic diagnostic imaging of harmonic frequencies with speckle reduction processing
US6547732B2 (en) * 1998-10-01 2003-04-15 Koninklijke Philips Electronics N.V. Adaptive image processing for spatial compounding
JP2004290393A (ja) * 2003-03-26 2004-10-21 Matsushita Electric Ind Co Ltd 超音波診断装置
US20050124886A1 (en) * 2003-11-21 2005-06-09 Koninklijke Philips Electronics N.V. System and method for generating ultrasound images having variable spatial compounding
JP2007513727A (ja) * 2003-12-16 2007-05-31 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 空間合成による超音波診断造影画像
US20070049825A1 (en) * 2005-08-22 2007-03-01 Medison Co., Ltd. System and method of forming an ultrasound spatial compound image
KR101120820B1 (ko) * 2009-11-19 2012-03-22 삼성메디슨 주식회사 초음파 공간 합성 영상을 제공하는 초음파 시스템 및 방법
KR101183003B1 (ko) * 2010-01-25 2012-09-18 삼성메디슨 주식회사 마스크에 기초하여 초음파 공간 합성 영상을 제공하는 초음파 시스템 및 방법
KR101386098B1 (ko) * 2010-11-10 2014-04-29 삼성메디슨 주식회사 빔 프로파일에 기초하여 초음파 공간 합성 영상의 화질을 개선시키는 초음파 시스템 및 방법

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6224552B1 (en) * 1998-10-01 2001-05-01 Atl Ultrasound Ultrasonic diagnostic imaging system with reduced spatial compounding seam artifacts
CN101199430A (zh) * 2006-12-15 2008-06-18 深圳迈瑞生物医疗电子股份有限公司 空间复合成像方法、设备及其超声成像系统
CN101866480A (zh) * 2009-04-15 2010-10-20 深圳迈瑞生物医疗电子股份有限公司 一种空间复合成像中的帧相关处理方法及系统

Also Published As

Publication number Publication date
US20140358003A1 (en) 2014-12-04
WO2013105339A1 (ja) 2013-07-18
JP5272084B2 (ja) 2013-08-28
CN104039235A (zh) 2014-09-10
JP2013141519A (ja) 2013-07-22

Similar Documents

Publication Publication Date Title
CN104039235B (zh) 超声波诊断装置
CN100518656C (zh) 超声波摄像装置及其摄像数据处理方法以及处理程序
JP6393703B2 (ja) サブボリュームの連続適応強化超音波イメージング
US6988991B2 (en) Three-dimensional ultrasound imaging method and apparatus using lateral distance correlation function
CN106469461B (zh) 视图方向自适应体积超声成像
US20200341142A1 (en) Motion-compensated processing
KR102014504B1 (ko) 초음파 이미징에서의 섀도우 억제
KR101175398B1 (ko) 합성영상을 제공하는 초음파 시스템 및 방법
KR100954988B1 (ko) 초음파 영상을 형성하는 초음파 시스템 및 방법
US9402600B2 (en) 3-dimensional elastic image generation method and ultrasonic diagnostic apparatus
US20200375574A1 (en) Ultrasound imaging spatial compounding method and system
US20180325492A1 (en) Color Doppler Imaging with Line Artifact Reduction
JP5525930B2 (ja) 三次元超音波画像を生成表示する超音波診断装置
JP2016514530A (ja) 超音波診断イメージングシステム。
CN102429686A (zh) 超声诊断仪的空间复合成像方法
US10561402B2 (en) Motion imaging with multiple parallel receive beams
KR101121267B1 (ko) 스티어링 각도를 고려하여 합성영상을 제공하는 초음파 시스템 및 방법
CN109982646B (zh) 超声波诊断装置和图像处理方法
CN102423265B (zh) 超声诊断仪的复合成像方法
JP5950291B1 (ja) 超音波診断装置及びプログラム
JP2017195913A (ja) 超音波診断装置
JP2000271118A (ja) 超音波診断装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20170216

Address after: Tokyo, Japan, Japan

Patentee after: Hitachi Ltd.

Address before: Tokyo, Japan

Patentee before: Hitachi Aloka Medical Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211109

Address after: Chiba County, Japan

Patentee after: Fujifilm medical health Co.,Ltd.

Address before: Tokyo, Japan

Patentee before: Hitachi, Ltd.