CN103796594A - 超声波成像装置 - Google Patents

超声波成像装置 Download PDF

Info

Publication number
CN103796594A
CN103796594A CN201280042033.1A CN201280042033A CN103796594A CN 103796594 A CN103796594 A CN 103796594A CN 201280042033 A CN201280042033 A CN 201280042033A CN 103796594 A CN103796594 A CN 103796594A
Authority
CN
China
Prior art keywords
similarity
operational part
imaging apparatus
ultrasonic imaging
phasing
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.)
Granted
Application number
CN201280042033.1A
Other languages
English (en)
Other versions
CN103796594B (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.)
Hitachi Ltd
Original Assignee
Hitachi Medical Corp
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 Medical Corp filed Critical Hitachi Medical Corp
Publication of CN103796594A publication Critical patent/CN103796594A/zh
Application granted granted Critical
Publication of CN103796594B publication Critical patent/CN103796594B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046—Techniques for image enhancement involving transmitter or receiver
    • G01S7/52047—Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/18—Methods or devices for transmitting, conducting or directing sound
    • G10K11/26—Sound-focusing or directing, e.g. scanning
    • G10K11/34—Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
    • G10K11/341—Circuits therefor
    • G10K11/346—Circuits therefor using phase variation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/60—Control of cameras or camera modules
    • H04N23/63—Control of cameras or camera modules by using electronic viewfinders
    • H04N23/631—Graphical user interfaces [GUI] specially adapted for controlling image capture or setting capture parameters
    • H04N23/632—Graphical user interfaces [GUI] specially adapted for controlling image capture or setting capture parameters for displaying or modifying preview images prior to image capturing, e.g. variety of image resolutions or capturing parameters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4405—Device being mounted on a trolley
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461—Displaying means of special interest
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G01S15/8927—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array using simultaneously or sequentially two or more subarrays or subapertures

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Health & Medical Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Multimedia (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Signal Processing (AREA)
  • Human Computer Interaction (AREA)
  • Quality & Reliability (AREA)
  • Theoretical Computer Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

用多个元件(401)接收来自被检测体的超声波信号，通过相似度运算部(404)求出接收信号x(n)间的相似度。利用相似度运算部(404)求出的接收信号间的相似度C(n)，通过自适应权重运算部(407)求出相似度所对应的自适应权重w(n)。利用自适应权重w(n)与接收信号x(n)，由定相运算部(408)生成定相输出。图像处理部(108)利用定相输出来生成图像数据。例如，相似度运算部(404)在时间方向上进行相似度的运算。因此，利用自适应信号处理的方法，能够以少的运算量高精度地求出超声波成像装置的接收信号的定相处理所利用的权重值。

Description

超声波成像装置
技术领域
本发明涉及利用超声波对被检测体内的图像进行成像的超声波成像技术。
背景技术
所谓超声波成像技术，是利用超声波(不愿听到的声波、一般而言是20kHz以上的高频声波)以非侵袭的方式对以人体为首的被检测体的内部进行图像化的技术。作为一个例子，对医用超声波成像装置简单地进行说明。超声波探针朝向患者的体内发送超声波并接收从患者体内反射的回波信号。接收信号在超声波探针及超声波成像装置主体的一方或双方中被实施信号处理后交接给图像显示部，从而显示超声波图像。若更详细地说明的话，例如，在超声波成像装置主体的发送波束成形器中，生成发送波束的信号，经过了收发分离电路后被传送至超声波探针。超声波探针发射超声波。超声波探针在接收到来自体内的回波信号后，向成像装置主体传递信号。在成像装置主体中，接收信号经过收发分离电路及接收波束成形器后被传递给图像处理部。在图像处理部中，各种滤波器、扫描变换器等执行各种各样的图像处理。最终在图像显示部上显示超声波图像。
这样，一般的超声波诊断装置通过发送波束成形、接收波束成形及后级的图像处理这3种技术构成。特别是，发送及接收时的波束成形器进行RF(高频)电平下的信号处理，因此波束成形器的算法或安装体系机构决定超声波图像的基本画质。因此，波束成形器是装置的基干部。
接收波束成形器针对构成超声波探针的多个元件的各接收信号(接收数据)，依据焦点位置与元件位置的关系，提供延迟量分布成凹面型的延迟时间，使焦点(focus)对准空间的虚拟的某一点之后，对接收信号数据进行相加。该方法被称为基于延迟相加方式的定相。在该延迟相加方式中，通过将由超声波诊断装置的多个元件接收到的接收数据和诊断装置所蓄积的固定权重向量相乘这样的处理单元，来实现延迟。不仅是在接收波束成形器中这样实现，在发送波束成形器中也是同样的。
另一方面，作为超声波成像装置的基本问题，公知方位方向分辨率的制约。因为超声波的收发是利用具备有限的开口径的阵列来进行的，故发生开口部边缘的衍射的影响。如果准备无限长的阵列，则虽然存在与深度方向同样地可无限提高分辨率的可能性，但在现实中，存在收发阵列的长度等装置设计上的物理性制约，因此会妨碍方位方向分辨率的提高。近年来，以下尝试倍受关注：针对时序的收发数据的每一个，使在波束成形器的延迟相加时为了进行延迟而利用的前述的固定权重向量自适应地变化，从而得到更精细的超声波图像。由此，存在可显著提高波束成形技术中的本质性课题之一的方位方向分辨率的可能性。
特别是，例如，专利文献1公开了将以近年来在移动体通信领域发展的MVDR法(Minimum Variance Distortionless Response；Capon法)为首的自适应信号处理的技术利用于超声波图像处理。通过利用自适应的方法，从而按接收信号在时间方向的每个采样点来求出以往是固定值的权重向量，并将其与接收信号相乘，由此可以进行延迟。
【在先技术文献】
【专利文献】
【专利文献1】美国专利申请公开第10/676777号说明书
发明内容
-发明想要解决的课题-
在超声波诊断装置中，在将以MVDR法为首的自适应信号处理的技术适用于波束成形器的技术(自适应波束成形器)中，若在时间方向上不进行推测，而是逐次(每个样本·每个通道)计算空间方向的相关矩阵，则难以使在时间方向上分散的能量的误差收敛。因此，图像中的点像在时间方向(深度方向)上模糊，产生针对各种噪音处理的不稳定度，成为图像噪声/误差的原因。
另一方面，在通过自适应信号处理进行了时间方向的推测和空间方向的推测这两种推测的情况下，造成庞大的处理负荷，安装成本会增大。这些处理能力和推测处理负荷的折衷对于安装来说成为大的障碍。
本发明的目的在于，利用自适应信号处理的方法且以较少的运算量来高精度地求出超声波成像装置的接收信号的定相处理所利用的权重值。
-解决课题的手段-
为了达成上述目的，根据本发明可提供如下的超声波成像装置。即，一种超声波成像装置，具有：多个元件，接收来自被检测体的超声波信号；相似度运算部，求出多个元件的接收信号间的相似度；自适应权重运算部，利用相似度运算部求出的接收信号间的相似度，求出与该相似度对应的自适应权重；定相运算部，利用自适应权重与接收信号来生成定相输出；以及图像处理部，利用定相输出来生成图像数据。
-发明效果-
在本发明中，预先对接收信号进行相似度运算处理，利用求出的相似度来运算自适应权重，由此可以降低运算量并进行正确的点像的推测。例如，通过进行时间方向的相似度运算处理，从而可以用较少的运算量来补偿时间方向的偏离，能进行更正确的点像的推测。由此，能修正时间方向(深度方向)上的图像模糊，可获得直径小的点像，并且可以稳定地得到伪像或噪音也被降低的超声波图像。
附图说明
图1(a)是表示本实施方式的超声波成像装置的概略构成的立体图，(b)是框图。
图2是表示第1实施方式的接收波束成形器的构成的框图。
图3是表示第1实施方式的接收波束成形器的信号处理的流程图。
图4是表示第1实施方式的接收波束成形器的动作的说明图。
图5(a)是表示通过第1实施方式的延迟电路412的处理使波面一致的接收信号208的说明图、(b)是表示通过相似度运算部404求出的相似度函数、和提取变换部413提取的指标值的例子的说明图、(c)是表示使波面一致的接收信号的采样时间n的输出在通道方向上的分布的图表、(d)是表示相似度函数的采样时间n的指标(φ)在通道方向上的分布的图表。
图6是表示比较例的接收波束成形器的动作的说明图。
图7是表示第2实施方式的接收波束成形器的构成的框图。
图8是表示第2实施方式的接收波束成形器的信号处理的流程图。
图9是表示第3实施方式的接收波束成形器的构成的框图。
图10是表示第4实施方式的接收波束成形器的构成的框图。
图11是表示第5实施方式的接收波束成形器的构成的框图。
图12(a)是表示第6实施方式的接收波束成形器的部分结构的框图、(b)是表示抽取(decimate)运算处理的说明图。
图13是本实施方式的超声波成像装置的中央控制台的立体图。
图14(a)是通过第1实施方式的超声波成像装置得到的图像、(b)是通过现有技术型的延迟相加方式(无相似度运算处理与自适应定相引擎405)而得到的比较例的图像、(c)是仅利用自适应定相引擎在时间方向上进行平均化处理而得到的比较例的图像、(d)是仅利用自适应定相引擎且在时间方向上未进行平均化而得到的比较例的图像。
图15(a)是沿着通道方向描绘图14(a)～(d)的超声波图像的最大亮度点所处的深度的亮度而得到的图表、(b)是沿着深度方向描绘图14(a)～(d)的超声波图像的最大亮度点所处的通道的亮度而得到的图表。
图16是表示第7实施方式的接收波束成形器的构成的框图。
图17是表示第7实施方式的接收波束成形器的信号处理的流程图。
具体实施方式
根据本发明的第1形态，可提供以下的超声波成像装置。即，一种超声波成像装置，具有：多个元件，接收来自被检测体的超声波信号；相似度运算部，求出多个元件的接收信号间的相似度；自适应权重运算部，利用相似度运算部求出的接收信号间的相似度，求出与该相似度对应的自适应权重；定相运算部，利用自适应权重与接收信号来生成定相输出；以及图像处理部，利用定相输出来生成图像数据。通过采取这种构成，从而可以利用自适应信号处理的方法，以少的运算量高精度地求出超声波成像装置的接收信号的定相处理所利用的权重值。
相似度运算部进行相似度运算的方向优选是时间方向。
自适应权重运算部例如构成为：通过利用相似度运算部求出的相似度来进行自适应信号处理，从而求出自适应权重。
优选在多个元件与相似度运算部之间还配置延迟部，该延迟部依据超声波信号的焦点位置，分别使多个元件接收到的接收信号延迟，以使波面一致。由此，相似度运算部可以求出该延迟部延迟后的接收信号的相似度。
再有，例如采取以下构成：在相似度运算部与自适应权重运算部之间配置提取表示相似度的特征的规定的指标值的提取部，自适应权重运算部将提取部提取出的指标值用作相似度。
自适应权重运算部例如构成为根据接收信号间的相似度来生成空间协方差矩阵，进行自适应信号处理，由此求出自适应权重。
可以并排配置接收超声波信号的多个元件。该情况下，相似度运算部采用求出多个元件之中仅相隔了规定数量的2个元件的接收信号间的相似度的构成。
在自适应权重运算部求出的自适应权重的数量比接收信号的数量还少的情况下，定相运算部也可以根据自适应权重的数量进行使多个接收信号衰退的运算，利用衰退后的接收信号与自适应权重来生成定相输出。
根据本发明的第2形态，可提供如下的超声波成像装置。即，一种超声波成像装置，具有：多个元件，接收来自被检测体的超声波信号；相似度运算部，求出多个元件的接收信号间的相似度；自适应权重运算部，利用相似度运算部求出的接收信号间的相似度，求出自适应权重；定相运算部，利用自适应权重与接收信号来生成定相输出；以及图像处理部，利用定相输出来生成图像数据。这种形态的超声波成像装置可以利用自适应信号处理的方法，以少的运算量高精度地求出超声波成像装置的接收信号的定相处理所利用的权重值。
上述第2形态中，可以在相似度运算部与所述自适应权重运算部之间配置：提取部，提取表示相似度的特征的规定的指标值；以及延迟部，基于由提取部提取出的指标值，分别使多个元件的接收信号延迟，并使波面一致。该情况下，自适应权重运算部可以利用延迟部延迟后的接收信号来求出自适应权重。
这样，在本发明中，预先对接收信号进行相似度运算处理，利用求出的相似度、或者根据求出的相似度而进行延迟后的接收信号，运算自适应权重，由此可以降低运算量并进行正确的点像的推测。例如，通过进行时间方向的相似度运算处理，从而能以较少的运算量补偿时间方向的偏离，可以进行更正确的点像的推测。由此，可以修正时间方向(深度方向)的图像模糊，可获得直径小的点像，并且能稳定地得到伪像或噪音也被降低的超声波图像。
对本发明的一实施方式的具体例进行说明。
(第1实施方式)
将上述本发明的第1形态的超声波成像装置作为第1实施方式具体地进行说明。
首先，利用图1(a)、(b)对超声波成像装置的整体构成进行说明。图1(a)是装置的立体图，图1(b)是表示内部的概略构成的框图。
如图1(a)所示，超声波成像装置具备超声波探针101、装置主体102和图像显示部103。如图1(b)所示，在装置主体102内配置着发送波束成形器104、收发分离电路411、接收波束成形器107和图像处理部108。
在发送波束成形器104中生成发送波束用的信号，经过收发分离电路411后被交接给超声波探针101。超声波探针101朝向被检测体100的体内发送超声波，且超声波探针101接收体内反射的回波信号。接收信号经过收发分离电路411后，在接收波束成形器107中被实施定相运算处理等。定相运算处理后的接收信号被交接给图像处理部108，在各种滤波器、扫描变换器等中执行各种各样的图像处理，生成超声波图像。超声波图像被交接给图像显示部103而被显示。
图2是表示接收波束成形器107的构成的框图。该接收波束成形器是将自适应信号处理的技术适用于波束成形器的自适应波束成形器。如图2所示，接收波束成形器107包括延迟电路412、相似度运算部404、采样数量调整部410、提取部变换部413和自适应定相引擎405。这些接收波束成形器107的各部分既能由独立的电路分别构成，也能构成为通过预先保存有程序的存储器和读取该程序并执行的CPU、GPU来实现各部分的动作。
相似度运算部404通过运算来求出构成超声波探针101的多个元件的接收数据(以下也称为接收数据)彼此在时间方向上的相似度，并将基于运算结果而计算出的信息作为自适应定相引擎405的输入。通过预先进行时间方向上的相似度运算处理，从而在自适应定相引擎405中能够以较少的运算量补偿时间方向的偏离，可以进行更正确的点像推测。另外，作为相似度运算部404的前级而配置延迟电路412，依据元件的位置而向构成超声波探针的多个元件的各接收信号提供延迟时间，进行使焦点(焦点)对准空间的虚拟的某一点的处理。
超声波探针101具备排列成阵列状的多个元件(超声波振荡器)400。在本实施方式中采用有效通道(active channel)技术，将接收了1个发送超声波波束所对应的接收回波的超声波探针101之中、元件400中的一部分区域的元件设为有效通道401，利用有效通道401的接收信号，生成超声波传播方向的1个图像数据(1簇)。如图2所示，将元件的位置稍微错开而依次构成多个有效通道402、401、403，针对各个有效通道402、401、403来生成簇，将该结果排列而成的图形就是超声波图像。
在以下的说明中，通过对1次收发所对应的1个有效通道401的各元件的接收数据进行自适应定相处理，从而对生成1簇之际的各部的动作进行说明。图3是表示此时的接收波束成形器的处理的流程。
由有效通道401接收到的多个接收数据通过收发分离电路411后，被输入至接收波束成形器107的延迟电路412。
在步骤31中，如图4所示，延迟电路412依据元件的位置，向构成有效通道401的多个元件400的各接收信号(接收数据)提供以空间内的一点201为中心且延迟量分布成凹面型205的延迟时间，由此假使焦点(focus)对准焦点201(步骤31)。通过将凹面型205的形状例如改变成凹面型206、207，从而可以使焦点分别对准点202、203。由此，根据所期望的焦点，使各元件400的接收信号204延迟，可以得到波面1000一致的接收信号(接收数据)208。
如图5(a)所示，如果构成有效通道401的元件数量(通道数量)为K个，则某一采样时刻n的K个接收数据208可以表示为下述式(1)的向量x(n)。
【数式1】
x(n)=[x1(n)，x2(n)，…，xK(n)]T···(1)
在步骤32中，配置于延迟电路412后级的相似度运算部404从延迟电路412接受由K个接收数据组成的向量x(n)作为输入信号，进行不同通道(元件)的接收信号间的相似度运算。具体地，在相似度运算部404中，计算有效通道(总数K)的通道间的相似度函数，并输出该结果。作为相似度函数，也可以利用马氏距离(Mahalanobis distance)、皮尔逊相似度函数、互相关函数等输出多个信号向量间的相似度的任何函数。
在此，作为相似度运算的一例，对利用了互相关函数的运算进行说明。互相关函数是表示某一信号与其他信号的相似度的方法之一，一般而言可用下式(2)的函数Cp(n)来表示。如式(2)所示，利用某一通道p的接收数据xp(n)、与从通道p相隔了q通道的通道p+q的接收数据xp+q(n)在时间方向上仅浮动τ后取得共轭的信号x* p+q(n+τ)的交迭，表示互相关函数Cp(n)。在此，将通道p与通道p+q设定为相隔几个通道是任意的，只要满足式(3)q就可以是任意值。例如，若q=1，则Cp(n)成为相邻通道间的互相关函数。在总通道数为K时，如式(4)所表示那样，输出K-q个互相关函数。例如，在q=3时取得与相隔了3个通道的通道之间的相关，在(ch.1，ch.4)、(ch2，ch5)…(ch.K-3，ch.K)全部中输出(K-3)组的互相关函数C1(n)～CK-q(n)。再有，在式(2)中，积分区间-r～r表示图5(a)所示的互相关窗1003的区间，针对互相关窗1003内的每个采样时刻的接收数据x(n)进行式(2)的运算。
【数式2】
C p ( n ) = ∫ - r r x p ( n ) x p + q * ( n + τ ) dτ . . . ( 2 )
【数式3】
1≤q≤K-2···(3)
【数式4】
1≤p≤K-q···(4)
图5(b)表示通过式(2)求出的互相关函数C1(n)～CK-q(n)。其中，在图5(b)中，例示了q=1的情况。
另外，互相关窗1003的大小既可以是预先确定的固定值，也能构成为采样点调整部410依据操作者的指示而设定成任意的大小。具体地，在式(2)中，通过使积分区间r的大小变化，从而可以变更互相关窗1003的大小。即，采样点调整部410作为针对接收信号的时间方向设定任意长度的窗的窗长度调整部发挥功能。
向提取变换部413交接通过式(2)运算出的K-q个互相关函数C1(n)～CK-q(n)。在步骤33中，提取变换部413如式(5)及式(6)所示提取表示互相关函数Cp(n)的特征的规定的1个以上的指标(参数)值。按K-q个互相关函数C1(n)～CK-q(n)的每1个提取参数值。作为参数，利用Cp(n)的时间方向的最大值(峰值振幅)ap、取最大值ap的时刻与基准时刻t0的偏离量Δtp(n)、将Δtp(n)变换成相位的φp(n)、式(7)所示的由最大值ap和φp(n)表示的复数(复数数据)ξp(n)、式(8)所示的复数分量Ip与Qp的组、及仅复数分量的实部或仅虚部之中的规定的1个以上的参数。
将提取变换部413提取出的规定的参数值作为自适应定相引擎405的输入。在作为参数而利用了复数数据的情况下，通过利用相位φp(n)与振幅ap双方，从而在自适应定相引擎405中，可以进行高精度的互相关推测。其中，式(7)中的T表示超声波的周期。
【数式5】
a p ( n ) = max [ C p ( n ) ] . . . ( 5 )
【数式6】
Δt p ( n ) = t max [ C p ( n ) ] - t 0 . . . ( 6 )
【数式7】
ξ p ( n ) = a p ( n ) exp ( jφ p ) , φ p = 2 π Δt p T . . . ( 7 )
【数式8】
Ip=αpcosφp(n)，Qp=αpsinφp(n)···(8)
如图2所示，在提取变换部413的后级配置有自适应定相引擎405。自适应定相引擎405是基于多个输入信号来生成定相输出的模块，是自适应波束成形器(接收波束成形器)107的主运算部。如图2所示，自适应定相引擎405具备矩阵运算部406、自适应权重运算部407及定相运算部408。在本实施方式中，作为自适应权重运算部407的算法，以利用MVDR算法的情况为例，对自适应定相引擎405的处理进行说明。
在步骤34中，矩阵运算部406利用式(9)来计算空间协方差矩阵R(n)。利用提取变换部413提取出的ap、Δtp(n)、φp(n)、ξp、及、Ip与Qp的组之中的规定的1个以上的值，运算R(n)。在此，以利用复数数据ξp来求出式(9)的空间协方差矩阵R(n)的情况为例进行说明。如式(9)所示，R(n)是取了用式(10)表现的复数向量ξ(n)和其复数转置向量ξH(n)之积的系综平均而得的值。
【数式9】
Figure BDA0000470637980000102
【数式10】
ξ(n)=[ξ1(n)，ξ2(n)，...，ξK-q(n)]T···(10)
本发明的特征在于：作为式(9)的R(n)的输入(要素)，利用由相似度运算部404运算出的相似度。由于利用相似度，故将接收的有效阵列401的大小设为K，空间协方差矩阵R(n)成为(K-q)×(K-q)的方矩阵。如果是现有技术中的自适应波束成形器，则对空间协方差矩阵的输入使用式(1)的x(n)，因此协方差矩阵成为K×K的方矩阵。
另外，在式(9)中，系综平均数N，作为对象采样点ξ(n)前后各S采样，合计N=2S+1点，如式(9)的最右边所示取一致平均。
在步骤35中，在接受了空间协方差矩阵R(n)的自适应权重运算部407中，利用MVDR法来计算权重向量w(n)。在此，通过式(11)来求出基于MVDR法的权重向量。
【数式11】
w ( n ) = R - 1 ( n ) a H R - 1 ( n ) a . . . ( 11 )
式(11)中，R(n)是由式(9)生成的、时间方向的某一采样点n处的相关矩阵，a是模向量。
例如，在将MVDR法适用于线性扫描的某一时刻的接收信号采样的情况下，将到来方向q=0°作为关心方向而得到的复数权重向量w(n)是使关心方向以外的方向的响应最小化的自适应滤波器，可期待方位方向的分辨率提高在。本实施方式中，如上所述，由于在延迟电路412中实施延迟处理，故输入信号成为使波面1000在q=0°方向一致的数据。因此，式(9)中模向量a可以简单地设为a=[1，1，…，1]T。
在步骤36中，定相运算部408接受由自适应权重运算部407求出的复数权重向量w(n)，如式(12)～式(14)所示运算迂回延迟电路412而接受的接收数据向量x(n)。由此，定相运算部408得到与有效通道401对应的1簇定相输出y(n)。
式(12)、式(13)是在从延迟电路412进行了旁路输入的时间点使作为K个分量的接收数据的向量x(n)向由互相关函数的个数K-q的要素构成的向量z(n)衰退的梯形权重处理。在从K个要素向K-q个要素的变换中，只要可以从K个要素向K-q衰退，任意运算都可以。由此，也能利用与式(12)、式(13)的梯形权重处理不同的线形运算。
【数式12】
x′p(n)=[xp(n)，xp+i(n)，...，xp+K-q-1(n)]T···(12)
【数式13】
z ( n ) = Σ p = 1 q + 1 x p ' ( n ) = [ z 1 ( n ) , z 2 ( n ) , . . . z K - q ( n ) ] T . . . ( 13 )
【数式14】
y(n)=WH(n)z(n)···(14)
伴随于在接收阵列上从有效通道401移位至有效通道402、403，按每1簇将可通过式(14)获得的1簇的定相输出y(n)交接给图像处理部108。图像处理部108中，由扫描变换器排列所有簇，以生成二维图像。再有，进行各种滤波处理等各种各样的后端图像处理。最终，超声波图像被显示于图像显示部103。
这样，在本实施方式中，在自适应方法中，按接收信号x(n)在时间方向上的每个采样点，利用接收信号并通过运算而求出权重向量w(n)，对该w(n)与x(n)进行运算处理，由此可以得到定相输出y(n)。由此，作为比较例，如图6所示，与利用固定值的权重向量w的情况相比较，由于可以使权重向量自适应地变化，故可以得到高精细的超声波图像。
而且，在本实施方式中，由于预先进行时间方向的相似度运算处理，故在自适应定相引擎405中以较少的运算量就能补偿时间方向的偏离，可以进行更正确的点像的推测。由此，可以修正时间方向(深度方向)的图像模糊，可得到更紧凑的点像。还可以稳定地得到伪像或噪音少的超声波图像。
关于进行相似度运算处理的效果，利用图5(c)、(d)进行说明。图5(c)在纵轴上标示了在某一采样时刻n图5(a)的延迟电路412的输出x(n)，在横轴上标示了通道编号。图5(d)在纵轴上标示了在某一采样时刻n提取变换部413提取出的值、即相似度运算部404的运算结果Cp的输出的特征(参数：ap、Δtp(n)、φp(n)、ξp(n)或者、复数分量Ip与Qp组的任一个，在此作为一例标示了φp(n))，在横轴上标示了通道编号。
如图5(c)所示，接收信号在通道方向的描绘直接排列了声波的RF信号(接收信号)，因此有时会噪声多，因微小散射体或噪声的影响在描绘中产生突起1009。通过在时间方向上取平均，从而可削弱突起1009的影响，但这种方法是有界限的。由此，在将图5(a)的RF信号直接取算术平均、或者在时间方向上取算术平均的基础上输入至自适应定相引擎的情况下，得到的权重w(n)受到噪声或突起1009的影响。特别是，在将MVDR法等的自适应波束成形适用于RF信号的情况下，按照使焦点对准数据突起1009的方式生成权重w(n)，在原本物体不存在的情况下，会在突起1009的位置上生成伪像。例如，在不进行相似度运算的情况下(即不进行时间方向推测的情况下)，自适应定相引擎的推测所利用的相位φ不得不利用与图5(d)示出的实际相位1010基本不一致的固定相位1016，因此时间方向的推测精度会下降，最终的超声波图像的画质会下降。
另一方面，在本实施方式中，相似度运算部404进行的相似度运算处理将波视作一个信息包来进行式(2)的匹配滤波处理，如图5(b)所示，可获得向时间方向的平均化作用，可以得到降低了噪声的通道方向分布1010。例如，通过利用根据相似度运算部404的运算结果而由提取变换部413提取出的相位1010，从而存在以下效果：根据时间方向可以进行高精度的推测，将最终的超声波图像的点像设为时间方向(超声波传播方向)的点直径小的紧凑的像。
另外，作为上述矩阵运算部406中的运算算法的其他方法之一，也能进行利用了子阵列矩阵的空间平均运算。利用式(15)、式(16)来表示子阵列矩阵。
【数式15】
R l ~ ( n ) = ξ l ~ ( n ) ξ l ~ H ( n ) . . . ( 15 )
【数式16】
ξ l ~ ( n ) = [ ξ l ( n ) , ξ l + 1 ( n ) , . . . , ξ l + L - 1 ( n ) ] T . . . ( 16 )
若使子阵列矩阵的主对角分量与相关矩阵R(n)的主对角分量一致，1个样本1个样本地错开，则成为(K-q)-L+1个子阵列矩阵的空间平均处理，可获得式(17)的相关矩阵R^(n)。在由自适应权重运算部406运算该相关矩阵R^(n)之际，通过置换成上述式(11)的R(n)并加以利用，从而可以运算权重w^(n)。在定相运算部408中，利用式(18)、式(19)来输出定相输出y(n)。
【数式17】
R ^ ( n ) = 1 N ( K - q - L + 1 ) Σ s = - S s Σ l = 1 K - q - L + 1 R l ~ ( n - s ) . . . ( 17 )
【数式18】
z l ~ ( n ) = [ z l ( n ) , z l + 1 ( n ) , . . . , z l + L - 1 ( n ) ] T . . . ( 18 )
【数式19】
y ( n ) = w ^ H ( n ) Σ l = 1 K - q - L + 1 z l ~ ( n ) . . . ( 19 )
这样，通过在矩阵运算部406中进行子阵列矩阵的空间平均处理，从而可获得抑制超声波接收信号的相关性的噪音的效果。再有，作为空间平均法的一种，例如也能进行非专利文献2所公开的公知的前置/后置空间平均处理。
在上述的本实施方式中，对在相似度运算部404中求出相隔q个的通道的接收信号间的相似度的例子进行了说明。例如，在q=1的相邻通道间求出接收信号的相似度的状况的情况下，相似度运算部404输出K-1个互相关函数，因此在自适应定相引擎405中进行运算的式(10)、式(11)的向量的要素数也为K-1个。因而，在对接收数据向量x(n)和复数权重向量w(n)进行运算处理之前，需要利用式(12)及式(13)，使接收数据的向量x(n)朝由互相关函数的个数K-1的要素构成的向量z(n)衰退。因此，为了避免式(12)、(13)的运算，也可以在相似度运算部404中如式(20)所示求出第K通道的输入数据XK(n)的自相关函数，作为相似度q=1时第K个互相关函数，将其用作CK(n)。
由此，如式(21)、式(22)所示，在保持通道数K不变的情况下能够准备式(10)以后的向量，可以不使旁路输入的接收信号x(n)衰退就加以利用。最终通过式(23)可以得到定相输出。
【数式20】
C K ( n ) = ∫ - r r x K ( n ) x K ( n + τ ) dτ . . . ( 20 )
【数式21】
ξ(n)=[ξ1(n)，ξ2(n)，…，ξK(n)]T···(21)
【数式22】
z(n)=x(n)=[x1(n)，x2(n)，…，xK(n)]T···(22)
【数式23】
y(n)=wH(n)x(n)···(23)
如以上所述，在本实施方式中，其特征在于：作为自适应波束成形的算法，进行来自接收阵列的多个元件的接收信号的相似度运算，利用该运算结果来生成空间协方差矩阵R(n)。由此，在自适应权重运算部407中，进行基于空间协方差矩阵R(n)的波束成形的算法也可以是任何算法。即，除了MVDR法，例如也可以利用MUSIC(Multiple SignalClassification)法、APES(Amplitude and Phase Estimation法、ESPRIT(Estimation of Signal Parameters via Rotational InvarianceTechniques)法、MEM法(Maximum Entropy Method)等。
(第2实施方式)
利用图7，对本发明的第2实施方式的超声波成像装置进行说明。图7是第2实施方式的接收波束成形器107的框图。图8是表示接收波束成形器107的各部分的动作的流程图。
第2实施方式中，自适应权重运算部构成为具有：预先保存有相似度的分布与权重值的多种组合的权重存储部、和权重推测部。权重推测部构成为：基于从相似度运算部接受的多个相似度的分布，选择权重存储部内所保存的相似度的分布与权重值的组合，从而选择与从相似度运算部接受的多个相似度对应的权重值。
具体地，如图7所示，第2实施方式的接收波束成形器107和第1实施方式的图4的接收波束成形器107在自适应定相引擎501的构成上存在不同。其他构成与第1实施方式同样，因此省略说明。
第2实施方式中，如图7所示，自适应定相引擎501构成为具备权重推测部502、权重存储部503和定相运算部504。在第1实施方式的图4的自适应定相引擎405中，按接收数据x(n)的时间方向的每个样本(或者隔几个点间隔剔除)进行相似度运算与自适应权重运算，以进行了逐次自适应定相，但第2实施方式的特征在于，不按每个采样进行权重运算。第2实施方式中，将预先设想的相似度提取值和每个该设想出的相似度提取值的权重向量w(n)建立对应后预先保存在权重存储部503(Look UpTable：LUT)中。权重推测部502根据权重存储部503的权重向量w(n)来推测权重。
具体地，在权重存储部503中通过事先的脱机处理，保存预先设想出的接收信号x(n)的相似度的规定参数(ap、Δtp(n)、φp(n)、ξp(n)、及、复数分量Ip与Qp的组之中规定的1个以上)的值(例如、φp(n))在通道方向的分布(例如n=1～(k-q))。另外，按每个上述参数值的分布，将预先通过运算而求出的权重向量w(n)和参数值的分布建立对应后预先保存于权重存储部503中。该权重向量w(n)的运算是利用第1实施方式的矩阵运算部406及自适应权重运算部407的自适应波束成形器的算法而求出的。
对本实施方式的接收波束成形器107的动作进行说明。图8的步骤31～33中，收发分离电路411、延迟电路412、相似度运算部404及提取变换部413与第1实施方式的图3的步骤31～33同样地，针对利用延迟处理使波面一致了的接收信号x(n)运算出相似度，根据相似度来提取规定的参数(例如φp(n)、n=1～(k-q))。该规定的参数是与权重存储部503保存有通道方向的分布的参数相同的参数。
步骤44中，权重推测部502比较从提取变换部413接受的表示相似度的规定参数在通道方向的分布、和权重存储部503所保存的数据，选择参数在通道方向的分布最接近的保存数据。选择参数在通道方向的分布最接近的保存数据之际，能够利用基于最大概似推测、最小二乘法、多项式近似的递归式拟合算法等已知的曲线拟合算法。权重推测部502将与选择出的参数在通道方向的分布的保存数据对应地保存着的权重向量w(n)交接给定相运算部504，以作为推测权重向量。
步骤45中，定相运算部504利用式(14)、式(19)、式(23)的任一个来运算出推测权重向量w(n)和迂回延迟电路412而接受的接收信号x(n)，得到与有效通道401对应的1簇的定相输出y(n)。
在第2实施方式的构成中，通过利用图7的自适应定相引擎501，从而不计算逐次权重，因此与利用图4的自适应定相引擎405的情况相比可以实现格外的运算成本的降低。
在本实施方式中，是将对接收信号实施了相似度运算的信号输入至自适应定相引擎501的构成。如在第1实施方式中利用图5(d)所说明过的那样，预先实施了相似度运算的信号在通道方向分布上噪音少、且信号的稳定性高。因此，第2实施方式中，即便在作为自适应定相引擎而利用了LUT型的自适应定相引擎501的情况下，也能够进行高精度的权重推测。例如，如图5(b)所示，在利用噪音或突起多的通道分布1006并利用了LUT型的自适应定相引擎501的情况下，虽然权重推测的精度低，但如本实施方式所示，在相似度运算部404中通过利用进行了互相关处理的分布，从而在权重推测部502中能够进行更正确的近似、例如多项式近似，利用了权重存储部503的通道权重推测的精度的得以提高。
(第3实施方式)
第3实施方式中，采用自适应权重运算部具有第1运算部和第2运算部的构成。第1运算部利用相似度运算部求出的相似度进行自适应信号处理来运算出自适应权重。第2运算部具备预先保存有相似度的分布和权重值的多种组合的权重存储部及权重推测部。权重推测部通过选择在权重存储部内所保存的相似度的分布和权重值的组合，从而选择与从相似度运算部接受的多个相似度对应的权重值。再有，自适应权重运算部具备选择性地驱动第1运算部或第2运算部的驱动切换部、和选择性地将第1运算部及第2运算部的输出的一方交接给图像处理部的输出切换部中的至少一方。
以下，利用图9具体地说明本发明的第3实施方式的超声波成像装置。图9是表示第3实施方式的接收波束成形器的一部分的框图。如图9所示，第3实施方式的接收波束成形器是作为自适应定相引擎而具备了第1实施方式的自适应定相引擎405与第2实施方式的定相引擎501的复式定相引擎的接收波束成形器。将第1实施方式的自适应定相引擎405称为逐次型自适应定相引擎，将第2实施方式的自适应定相引擎501称为LUT型自适应定相引擎。
第3实施方式中，具备用于选择性地使2种定相引擎501、405中的任一个工作的切换部603、和控制切换部603的控制部601。切换部603具备将提取变换部413的输出交接给2种定相引擎501、405中的任一个的开关602。
定相引擎501、405各自的结构及动作和第1及第2实施方式中说明过的内容相同，故在此省略说明。
在第3实施方式的构成中，超声波成像装置的做手术的人根据被检测体的超声波图像的对比度比的大小、超声波图像在时间/空间上的变动的大小来判断定相引擎的切换，可以向控制部601指示利用2种定相引擎501、405中的哪一个。由此，能够选择性地利用适于被检测体的状态的定相引擎来提供超声波图像。
(第4实施方式)
第4实施方式构成为：差分部求出第3实施方式中的第1运算部的运算结果与第2运算部的运算结果的差分，由权重变更部依据差分部求出的差分来变更第2运算部的权重存储部内所保存的权重值。
以下，利用图10具体地说明本发明的第4实施方式的超声波成像装置。图10是第4实施方式的接收波束成形器107的框图。
如图10所示，第4实施方式的接收波束成形器107在第3实施方式的图9的复式定相引擎501、405中还附加了反馈环。具体地，具备差分部701，其计算LUT型定相引擎501的定相输出y(n)和逐次型定相引擎405的定相输出y(n)的差分。差分部701的差分输出被输入至权重存储部503。权重变更部703被配置在权重存储部503中，基于来自差分部701的差分输出，对预先保存于权重存储部503中的权重w(n)的值进行变更。具体地，权重变更部703依据预先确定的算法进行运算，由此将差分部701的差分输出作为误差，按照误差变为最小的方式来变更在权重存储部503内预先保存的权重w(n)的值。
权重变更部703的运算算法只要是将误差设为最小的算法即可，可以是任何算法，例如优选利用类似于MMSE(Minimum Mean SquareError：最小误差二乘法)的算法。作为MMSE的算法可以利用基于鞍点法的LMS(Least-Mean Squares)、采样值的直接解法SMI(Sample MatrixInversion)、类似于卡尔曼滤波器的RLS(Recursive Least Square：递归式最小二乘法)等方法中的任一个。
这样，在第4实施方式中，可以进行使通过基于LUT型自适应定相引擎501的权重推测而求出的权重值w(n)靠近通过逐次型自适应定相引擎405运算的权重值w(n)的校准。因而，在拍摄开始时，使LUT型自适应定相引擎501与逐次型自适应定相引擎405都工作，由差分部701进行反馈，以使基于LUT型自适应定相引擎501的权重存储部503内的权重值最佳化。如果进行了最佳化，则使逐次型自适应定相引擎405与差分部701停止，仅进行基于LUT型自适应定相引擎501的权重推测，由此可以利用LUT型自适应定相引擎501来推测与逐次型自适应定相引擎405同样的权重w(n)并进行定相运算。
由此，与第2实施方式相比，虽然计算负荷会有些许上升，但由于可以将LUT型自适应定相引擎501的权重存储部503的权重值变更为最佳的权重值w(n)，故能够提供更适于被检测体的状态的超声波图像。
在本实施方式中，可以将对接收信号实施了相似度运算的噪音少且稳定性高的信号输入至自适应定相引擎405、501。由此，即便如第2实施方式中所描述的那样利用了LUT型的自适应定相引擎501的情况下，通道权重推测的精度也会较高。因此，本实施方式中，在利用了反馈环的复式自适应定相引擎的构成中，用作误差量的差分部的输出的扰动变小，可以稳定地进行反馈处理。进而，在权重变更部703中，在进行变更为通过MMSE等使误差最小的权重值的运算处理之际，可以实现反复计算次数的降低，进一步可降低运算成本。
(第5实施方式)
利用图11对本发明的第5实施方式的超声波成像装置进行说明。图11是表示第5实施方式的接收波束成形器107的一部分的框图。
图11的接收波束成形器107与第4实施方式相同，构成为：在具备基于差分部701的反馈环的复式定相引擎501、405中附加了定相引擎的切换部801、802、803。再有，配置有控制切换部801、802、803的切换动作的切换控制部806和向切换控制部806指示切换定时的定时控制部808。
在图11的构成中，依据来自切换控制部806的控制信号来对切换部801的开关804、805及切换部802、803进行切换，由此可以将逐次型定相引擎405的定相输出、LUT型定相引擎501的定相输出、或经过了基于差分部701的反馈环的来自LUT型定相引擎501的定相输出选择性地传递给图像处理部108。
再有，在图11的构成中，可以由定时控制部808来控制切换控制部806的切换定时。由此，例如在手术中的被检测体(患者)的活动较大的情况下等，使差分部701的反馈环返回，权重存储部503中的权重值的更新次数增多，并追踪运动。另一方面，在被检测体的活动少的情况下，基于差分部701的反馈环不返回，能够灵活地应对降低计算负荷、向图像处理部108中的其他图像处理算法分配计算成本等自适应超声波图像的运算方案或各种各样的应用程序。
(第6实施方式)
利用图12(a)、(b)来说明本发明的第6实施方式。图12(a)是说明本实施方式的相似度运算部404及提取变换部413的框图的图。图12(b)是表示图12(a)的抽取运算部的处理的说明图。
如图12(a)所示，在本实施方式中，在相似度运算部404内具备切换部902。切换部902进行以下切换，即，对从延迟电路403输入至相似度运算部404的信号与第1实施方式等同样地进行相似度运算处理，还是不进行相似度运算而是直接交接给提取变换部413。
再有，提取变换部413具备与第1实施方式等同样地进行规定参数的提取运算的提取运算部905、和抽取运算部904。提取运算部905内配置有切换部(参数(指标)切换部)903，其对从相似度运算部的运算结果Cp(n)中提取的参数种类进行切换。抽取运算部904对提取运算部903的输出进行间隔剔除处理后输出给定相引擎405、501。
再有，切换部902、903连接着对这些部件的动作进行控制的提取参数变更部901。
提取参数变更部901依据操作者(做手术的人)的指示对切换部902进行切换，由此可选择将向相似度运算部404输入的来自延迟电路403的信号直接交接给提取运算部905、交接进行完相似度运算的结果中的个。还有，通过对切换部903进行切换，从而可变更提取运算部905从自相似度运算部404接受的信号中提取哪个参数。即，可以选择并切换设为峰值振幅ap(n)、设为峰值的时刻偏离Δtp(n)、设为相位φ还是设为复数信号(ξp(n)、Ip(n)、Qp(n))。
抽取运算部904在时间方向上对从提取运算部905接受的相似度运算结果或延迟电路403的输出进行间隔剔除后进行输出。具体地，如图12(b)所示，抽取运算部904以规定的间隔将提取运算部905的输出间隔剔除后输出。即，在从采样时刻n到采样时刻n+4为止的期间1014内，继续输出采样时刻n的提取运算部905的输出和采样时刻n+5的输出的平均，在采样时刻n+5～n+9的期间1015内，继续输出采样时刻n+5的输出和采样时刻n+9的输出的平均。这样，在采样时刻的规定期间(例如5个采样)1014、1015内继续输出相同的值，由此可以降低后级的定相引擎405、501的运算处理的负荷。例如，每隔5个采样间隔剔除输出，由此可将定相引擎的处理成本降低至五分之一。再有，作为规定期间1014、1015的输出，利用与下一期间的输出之间的平均，由此可防止间隔剔除采样而导致输出变动较大的状况。
在第6实施方式中，其他构成可以利用第1～第5实施方式的任一个。
(超声波诊断装置的中央控制台)
图13表示本发明的超声波成像装置的中央控制台130。如图13所示，上述实施方式的几个特征可以配置于超声波诊断装置的中央控制台。例如，如图13所示作为中央控制台的开关部1103，可配置用于选择性地使2种定相引擎405、501工作的图9的控制部601、图11的切换控制部806、和图12的提取参数变更部901。由此，做手术的人可以切换定相引擎的种类或提取参数的种类。图11的切换定时控制部808或图2、图7、图10、图12的采样点调整部401作为带有刻度的提钮部1101、1102等而可以配置于中央控制台。由此，例如，在使LUT型定相引擎501运转几次的情况下，做手术的人可以以明示的方式改变使逐次型定相引擎405运转1次的重复定时值或相关采样点数(式(2)中的r)。由此，做手术的人一边观察实际的超声波图像，一边对自适应波束成形器的各种参数进行变更，可以利用患者依存的最佳条件进行诊断。再有，也可构成为将设定完的值显示于图像显示部103的显示区域的一部分1104上。
(上述实施方式的效果的例示)
利用图14(a)～(d)对上述实施方式的效果进行说明。图14(a)～(d)是利用本实施方式及比较例的接收波束成形器等将点散射体的回波信号做成超声波图像的图。在图14(a)～(d)中，对于所有图像，横轴是阵列状元件(超声波振荡器)400的通道方向(口径方向)，纵轴是超声波传播方向，成为朝下远离元件400的方向。
图14(a)的超声波图像1203是利用了相似度运算部404与自适应定相引擎405的第1实施方式而得到的图像，图14(b)的超声波图像1204是作为比较例通过现有技术的延迟相加方式(无相似度运算处理与自适应定相引擎405)而得到的图像。图14(c)的超声波图像1205是作为比较例而从第1实施方式移除了相似度运算部404与提取变换部413，仅利用自适应定相引擎而进行时间方向平均化处理(平均化的点数等同于超声波图像1203时的间隔剔除点数)得到的图像。图14(d)的超声波图像1206是仅利用自适应定相引擎但未进行时间方向平均化而得到的图像。
再有，图15(a)是沿着通道方向描绘了图14(a)～(d)的超声波图像1203～1206的最大亮度点所处的深度的亮度的图，图15(b)是沿着深度方向描绘了最大亮度点所处的通道的亮度的图。图15(a)、(b)中，以分贝显示亮度。
根据图14及图15(a)、(b)可知，第1实施方式得到的超声波图像1203与一点一点计算自适应权重的超声波图像1206相比毫不逊色，是点散射体的像紧凑的图像。再有，与在时间方向上计算了自适应权重的超声波图像1205相比，可知成为了点散射体的像小的点。还有，与现有技术的延迟相加方式的超声波图像1204相比，可知成为了点散射体的像小3倍以上的点，显著地提高了超声波图像的分辨率。
这样，根据图14及图15(a)、(b)可确认：发明与延迟相加方式相比，可显著地提高方位方向分辨率，与其他自适应波束成形方法相比可获得毫不逊色的性能和点直径小的紧凑的点像。
如上所述，在上述实施方式中，通过将相似度运算的结果作为自适应定相引擎的输入，从而可以使时间方向的点像的点直径缩减得较小，还能以低成本稳定地得到伪像或噪音少的超声波图像。由此，可以实现兼具高画质与高稳定性的自适应超声波成像装置。
(第7实施方式)
将上述本发明的第2形态的超声波成像装置作为第7实施方式具体地进行说明。
利用图16对本发明的第7实施方式的超声波成像装置进行说明。图16是第7实施方式的接收波束成形器107的框图。图17是表示接收波束成形器107的各部分的动作的流程图。
如图16所示，第7实施方式的接收波束成形器107和第1实施方式的图2的接收波束成形器107的构成的差异在于：在自适应定相引擎405与提取变换部413之间配置有第2延迟电路1600。再有，接收波束成形器107的动作代替第1实施方式的图3的流程的步骤34而如图17所示那样增加了步骤1700、1701。其他构成及动作和第1实施方式同样，因此省略说明。
第7实施方式中，如图17所示，在提取变换部413中按照根据相似度而计算出的相似度提取值，对接收数据x(n)再次实施延迟，使波面一致(步骤1700)。在第2延迟电路1600中，取代第1实施方式中的、提取出的规定指标(ap、Δtp(n)、φp(n)、ξp(n)、式(8)所示的复数分量Ip与Qp的组等)，而将再次使波形一致的数据x’(n)作为向自适应定相引擎405的输入，由此生成空间协方差矩阵R(n)(步骤1701)。
具体地，如第1实施方式所描述过的，提取变换部413通过互相关处理，可以从式(6)提取与基准时刻的偏离量Δtp(n)。在此，延迟电路1600中，利用该偏离量Δtp，再次对接收信号x(n)实施延迟，使波面一致(步骤1700)。用式(24)来表示通过延迟电路1600使波面一致了的数据x’(n)。在此，式(24)中的Δτp(n)是根据接收数据的采样频率对与基准时刻的偏离量Δtp(n)的刻度进行变换后的延迟采样点数。式(24)表示将第K通道取基准点时的实施例，但成为基准点的通道也可以是1～K的任一通道元件。
【数式24】
x′p(n)=xp(n-Δτp(n))···(24)
1≤p≤K-1
x′K(n)=xK(n)
式(24)表示本实施方式中、特别是作为优选形态而使式(3)、式(4)中的q=1的情况。然而，q不是1时，也可在通道方向对接收数据进行再次采样等在通道间使数据适当地衰退，由此适用第7实施方式。
由第2延迟电路1600使波面一致的数据x’(n)被输入至自适应定相引擎405，在矩阵运算部406中计算空间协方差矩阵(步骤1701)。因而，在第7实施方式中空间协方差矩阵R’(n)成为式(25)，利用该空间协方差矩阵R’(n)进行自适应权重运算、定相运算。以后的构成及信号处理顺序和第1实施方式是同样的。
【数式25】
Figure BDA0000470637980000251
这样，在第7实施方式的构成中，在具有固定参数的延迟电路中可以使并不一致的接收数据的相位一致。如图5(b)所示，偏离量Δtp表示接收元件在通道间的到来时刻的差异，对应于因被检测体内部的波面变形的影响而发生的时刻偏离量。利用再次使波面一致的接收数据x’(n)，由此在自适应定相引擎中的通道权重的推测中，能够利用更正确的值。
由此，根据本实施方式，通过将图16的第2延迟电路1600的输出作为自适应定相引擎的输入来利用，从而提高自适应定相引擎中的通道权重的推测精度。
-符号说明-
100  被检测体
101  超声波探针
102  超声波成像装置
103  图像显示部
104  发送波束成形器
107  接收波束成形器
108  图像处理部
130  中央控制台
201、202、203  空间内的某一点
205、206、207  延迟凹面型
208  使波面一致的接收信号
400  阵列状的元件(超声波振荡器)
401、402、403  有效通道
404  相似度运算部
405  逐次型自适应定相引擎
406  矩阵运算部
407  自适应权重运算部
408  定相运算部
410  采样点调整部
411  收发分离电路
412  延迟电路
413  提取变换部
501  LUT型自适应定相引擎
502  权重推测部
503  权重存储部
504  定相运算部
601  控制部
602  开关
701  差分部
703  权重变更部
801、802、803  切换部
804、805  开关
806  定相引擎切换控制部
808  切换定时控制部
901  提取参数变更部
902、903  切换部
904  抽取运算部
905  提取运算部
1000  波面
1003  互相关窗
1006  通常的接收信号的通道方向分布
1009  描绘中的突起数据
1010  进行完互相关处理的信号的通道方向分布
1011、1015  期间
1016  固定相位
1101、1102  带有刻度的提钮部
1103  开关部
1104  显示区域的一部分
1203  第1实施方式中得到的超声波图像
1204  延迟相加方式的超声波图像
1205  仅利用自适应定相引擎并进行了时间方向平均化时的超声波图像
1206  仅利用自适应定相引擎、但未进行时间方向的平均化时的超声波图像
1600  第2延迟电路

Claims (19)

1.一种超声波成像装置，其特征在于，具有：
多个元件，接收来自被检测体的超声波信号；
相似度运算部，求出所述多个元件的接收信号间的相似度；
自适应权重运算部，利用所述相似度运算部求出的所述接收信号间的相似度，求出与所述相似度对应的自适应权重；
定相运算部，利用所述自适应权重与所述接收信号，生成定相输出；和
图像处理部，利用所述定相输出来生成图像数据。
2.根据权利要求1所述的超声波成像装置，其特征在于，
所述相似度运算部在时间方向上进行所述相似度的运算。
3.根据权利要求1所述的超声波成像装置，其特征在于，
所述自适应权重运算部通过利用所述相似度运算部求出的所述相似度来进行自适应信号处理，从而求出所述自适应权重。
4.根据权利要求1所述的超声波成像装置，其特征在于，
在所述多个元件与所述相似度运算部之间还具有延迟部，该延迟部依据所述超声波信号的焦点位置，分别使所述多个元件接收到的接收信号延迟，从而使波面一致，
所述相似度运算部求出由所述延迟部延迟后的所述接收信号的相似度。
5.根据权利要求1所述的超声波成像装置，其特征在于，
在所述相似度运算部与所述自适应权重运算部之间配置对表示所述相似度的特征的规定的指标值进行提取的提取部，
所述自适应权重运算部将所述提取部提取出的指标值作为所述相似度来加以利用。
6.根据权利要求3所述的超声波成像装置，其特征在于，
所述自适应权重运算部根据所述接收信号间的所述相似度，生成空间协方差矩阵，进行自适应信号处理，由此求出所述自适应权重。
7.根据权利要求1所述的超声波成像装置，其特征在于，
所述多个元件被并排配置，
所述相似度运算部求出所述多个元件之中仅相隔规定数量的2个元件的接收信号间的相似度。
8.根据权利要求1所述的超声波成像装置，其特征在于，
在所述自适应权重运算部求出的所述自适应权重的数量比所述接收信号的数量还少的情况下，所述定相运算部根据所述自适应权重的数量，进行使所述多个接收信号衰退的运算，利用衰退后的接收信号和所述自适应权重来生成所述定相输出。
9.根据权利要求1所述的超声波成像装置，其特征在于，
所述自适应权重运算部具有：
权重存储部，预先保存了相似度的分布和权重值的多种组合；和
权重推测部，基于从所述相似度运算部接受的多个相似度的分布，选择在所述权重存储部内保存的所述相似度的分布与权重值的组合，由此选择与从所述相似度运算部接受的多个相似度相对应的权重值。
10.根据权利要求1所述的超声波成像装置，其特征在于，
所述自适应权重运算部具有第1运算部和第2运算部，
所述第1运算部利用所述相似度运算部求出的所述相似度，进行自适应信号处理，从而运算出所述自适应权重，
所述第2运算部具备：
权重存储部，预先保存了相似度的分布和权重值的多种组合；和
权重推测部，基于从所述相似度运算部接受的多个相似度的分布，选择在所述权重存储部内保存的所述相似度的分布与权重值的组合，由此选择与从所述相似度运算部接受的多个相似度相对应的权重值。
11.根据权利要求10所述的超声波成像装置，其特征在于，
该超声波成像装置还具有：
差分部，求出所述第1运算部的运算结果与所述第2运算部的运算结果的差分；及
权重变更部，其依据所述差分部求出的差分，变更在所述第2运算部的权重存储部内所保存的权重值。
12.根据权利要求10所述的超声波成像装置，其特征在于，
所述自适应权重运算部具有选择性地驱动所述第1运算部或所述第2运算部的驱动切换部、和选择性地将所述第1运算部及所述第2运算部的输出中的一方交接给所述图像处理部的输出切换部的至少一方。
13.根据权利要求12所述的超声波成像装置，其特征在于，
该超声波成像装置还具有：定时控制部，对所述驱动切换部及所述输出切换部的切换定时进行控制。
14.根据权利要求2所述的超声波成像装置，其特征在于，
该超声波成像装置还具有：窗长度调整部，在所述接收信号的时间方向上设定任意长度的窗，
所述相似度运算部针对由所述窗长度调整部设定的所述窗内的接收信号，求出所述接收信号间的相似度。
15.根据权利要求5所述的超声波成像装置，其特征在于，
所述相似度是函数，
表示所述相似度的特征的规定的所述指标值是振幅、相位、利用了所述振幅与所述相位的复数、所述复数的实部、及所述复数的虚部中的1个以上的值。
16.根据权利要求5所述的超声波成像装置，其特征在于，
所述提取部具有：指标切换部，作为表示所述相似度的特征的所述指标值，该指标切换部对振幅、相位、利用了所述振幅与所述相位的复数、所述复数的实部、及所述复数的虚部中的1个以上的值进行任意地切换后进行提取。
17.根据权利要求1所述的超声波成像装置，其特征在于，
所述相似度运算部具有在时间方向上对求出的所述相似度进行间隔剔除后输出的抽取运算部。
18.一种超声波成像装置，其特征在于，具有：
多个元件，接收来自被检测体的超声波信号；
相似度运算部，求出所述多个元件的接收信号间的相似度；
自适应权重运算部，利用所述相似度运算部求出的所述接收信号间的相似度，求出自适应权重；
定相运算部，利用所述自适应权重和所述接收信号，生成定相输出；和
图像处理部，利用所述定相输出来生成图像数据。
19.根据权利要求18所述的超声波成像装置，其特征在于，
在所述相似度运算部与所述自适应权重运算部之间配置：提取部，提取表示所述相似度的特征的规定的指标值；和延迟部，基于所述提取部提取出的指标值，分别使所述多个元件的接收信号延迟，以使波面一致，
所述自适应权重运算部利用由所述延迟部延迟后的所述接收信号，求出所述自适应权重。
CN201280042033.1A 2011-09-15 2012-08-07 超声波成像装置 Expired - Fee Related CN103796594B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2011-202342 2011-09-15
JP2011202342 2011-09-15
PCT/JP2012/070113 WO2013038847A1 (ja) 2011-09-15 2012-08-07 超音波撮像装置

Publications (2)

Publication Number Publication Date
CN103796594A true CN103796594A (zh) 2014-05-14
CN103796594B CN103796594B (zh) 2015-11-25

Family

ID=47883085

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280042033.1A Expired - Fee Related CN103796594B (zh) 2011-09-15 2012-08-07 超声波成像装置

Country Status (5)

Country Link
US (1) US9754185B2 (zh)
EP (1) EP2756805B1 (zh)
JP (1) JP5813776B2 (zh)
CN (1) CN103796594B (zh)
WO (1) WO2013038847A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108778147A (zh) * 2016-03-14 2018-11-09 富士胶片株式会社 超声波诊断装置及超声波诊断装置的控制方法
CN110087301A (zh) * 2018-01-25 2019-08-02 株式会社东芝 通信中继装置、控制方法以及程序

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104114099B (zh) * 2012-02-15 2016-02-17 株式会社日立制作所 超声波摄像装置
JP6143541B2 (ja) * 2013-05-10 2017-06-07 東芝メディカルシステムズ株式会社 超音波診断装置及び制御プログラム
JP5877816B2 (ja) * 2013-05-27 2016-03-08 日立アロカメディカル株式会社 超音波画像撮像装置及び超音波画像撮像方法
JP6197505B2 (ja) * 2013-09-05 2017-09-20 セイコーエプソン株式会社 超音波測定装置、超音波画像装置及び超音波測定方法
JP2015173879A (ja) * 2014-03-17 2015-10-05 古野電気株式会社 音速測定装置、音速測定方法、及び音速測定プログラム
US10175347B2 (en) 2015-12-02 2019-01-08 Butterfly Network, Inc. Ultrasound receiver circuitry and related apparatus and methods
JP6469746B2 (ja) * 2017-03-23 2019-02-13 キヤノンメディカルシステムズ株式会社 超音波診断装置及び制御プログラム
US10945706B2 (en) 2017-05-05 2021-03-16 Biim Ultrasound As Hand held ultrasound probe
EP3622320A1 (en) * 2017-05-11 2020-03-18 Koninklijke Philips N.V. Methods and systems for controlling the generation of a compound ultrasound image
CA3064439A1 (en) 2017-06-19 2018-12-27 Butterfly Network, Inc. Mesh-based digital microbeamforming for ultrasound applications
WO2018236778A1 (en) 2017-06-20 2018-12-27 Butterfly Network, Inc. MULTI-STAGE TRANSMITTED AMPLIFIER (TIA) FOR AN ULTRASONIC DEVICE
WO2018236799A1 (en) 2017-06-20 2018-12-27 Butterfly Network, Inc. ASYMMETRIC TRANS-IMPEDANCE (TIA) AMPLIFIER FOR ULTRASONIC DEVICE
US10857567B2 (en) 2017-06-20 2020-12-08 Butterfly Network, Inc. Analog to digital signal conversion in ultrasound device
EP3641656A4 (en) 2017-06-20 2021-03-17 Butterfly Network, Inc. AMPLIFIER WITH INTEGRATED TIME GAIN COMPENSATION FOR ULTRASOUND APPLICATIONS
EP3424434A1 (en) * 2017-07-07 2019-01-09 Koninklijke Philips N.V. Method and device for processing ultrasound signal data
CN108354626A (zh) * 2018-03-31 2018-08-03 华南理工大学 基于gpu的多种mv高清算法快速医学超声影像系统

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3474278B2 (ja) * 1994-03-16 2003-12-08 フクダ電子株式会社 超音波診断装置
US6193659B1 (en) * 1997-07-15 2001-02-27 Acuson Corporation Medical ultrasonic diagnostic imaging method and apparatus
JP4334032B2 (ja) * 1998-02-23 2009-09-16 株式会社東芝 超音波診断装置
US6135963A (en) * 1998-12-07 2000-10-24 General Electric Company Imaging system with transmit apodization using pulse width variation
JP2000350729A (ja) * 1999-06-10 2000-12-19 Ge Yokogawa Medical Systems Ltd 干渉アーチファクト抑制方法および超音波診断装置
CA2290240C (en) 1999-11-24 2008-03-11 Stergios Stergiopoulos High resolution 3d ultrasound imaging system deploying a multi-dimensional array of sensors and method for multi-dimensional beamforming sensor signals
JP4958348B2 (ja) * 2001-09-06 2012-06-20 株式会社日立メディコ 超音波撮像装置
US7280627B2 (en) * 2002-12-09 2007-10-09 The Johns Hopkins University Constrained data-adaptive signal rejector
EP1491914B1 (en) * 2003-06-25 2006-10-18 Aloka Co. Ltd. Ultrasound diagnosis apparatus comprising a 2D transducer with variable subarray shape pattern
US7207942B2 (en) * 2003-07-25 2007-04-24 Siemens Medical Solutions Usa, Inc. Adaptive grating lobe suppression in ultrasound imaging
US6894642B2 (en) * 2003-10-01 2005-05-17 University Of Florida Research Foundation, Inc. Doubly constrained robust capon beamformer
US7914454B2 (en) * 2004-06-25 2011-03-29 Wilk Ultrasound Of Canada, Inc. Real-time 3D ultrasonic imaging apparatus and method
US7740583B2 (en) * 2004-06-30 2010-06-22 General Electric Company Time delay estimation method and system for use in ultrasound imaging
US7867166B2 (en) * 2004-07-29 2011-01-11 University Of Rochester Statistical estimation of ultrasonic propagation parameters for aberration correction
WO2008051639A2 (en) * 2006-10-25 2008-05-02 Maui Imaging, Inc. Method and apparatus to produce ultrasonic images using multiple apertures
JP4974781B2 (ja) * 2007-06-26 2012-07-11 株式会社日立メディコ 超音波診断装置
US8968205B2 (en) * 2011-02-10 2015-03-03 Siemens Medical Solutions Usa, Inc. Sub-aperture control in high intensity focused ultrasound

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108778147A (zh) * 2016-03-14 2018-11-09 富士胶片株式会社 超声波诊断装置及超声波诊断装置的控制方法
CN108778147B (zh) * 2016-03-14 2021-04-09 富士胶片株式会社 超声波诊断装置及超声波诊断装置的控制方法
CN110087301A (zh) * 2018-01-25 2019-08-02 株式会社东芝 通信中继装置、控制方法以及程序
CN110087301B (zh) * 2018-01-25 2021-05-25 株式会社东芝 通信中继装置、控制方法以及记录介质

Also Published As

Publication number Publication date
WO2013038847A1 (ja) 2013-03-21
JPWO2013038847A1 (ja) 2015-03-26
EP2756805A4 (en) 2015-05-20
EP2756805B1 (en) 2017-11-08
EP2756805A1 (en) 2014-07-23
US9754185B2 (en) 2017-09-05
CN103796594B (zh) 2015-11-25
US20140240482A1 (en) 2014-08-28
JP5813776B2 (ja) 2015-11-17

Similar Documents

Publication Publication Date Title
CN103796594B (zh) 超声波成像装置
EP1620015B1 (en) Efficient ultrasound system for two-dimentional c-scan imaging and related method thereof
US8057392B2 (en) Efficient architecture for 3D and planar ultrasonic imaging—synthetic axial acquisition and method thereof
CN101209211B (zh) 接收孔径可调整的数字化超声波束合成方法及装置
US20090069693A1 (en) Retrospective Dynamic Transmit Focusing for Spatial Compounding
US10989810B2 (en) Systems and methods for beamforming using variable sampling
CN108056791A (zh) 超声波诊断装置的控制装置以及控制方法
CN104422931A (zh) 超声波测量装置、超声波图像装置及超声波测量方法
CN102695456B (zh) 波束形成方法、超声波诊断装置、集成电路
CN106796286A (zh) 超声信号模拟波束形成器/波束形成
CN101305924B (zh) 计算容积超声扫描
Yen Beamforming of sound from two-dimensional arrays using spatial matched filters
Zhou et al. Dynamic-transmit focusing using time-dependent focal zone and center frequency
CN105377142A (zh) 超声波诊断装置
Wygant et al. Beamforming and hardware design for a multichannel front-end integrated circuit for real-time 3D catheter-based ultrasonic imaging
Hechkel et al. Evaluation of transmit/receive apodization scheme for hardware implementation used in B-mode ultrasound imaging
Sreejeesh et al. Beamforming algorithm architectures for medical ultrasound
US20140350406A1 (en) Dynamic Operation for Subarrays in Medical Diagnostic Ultrasound Imaging
Son et al. A Novel Minimum Variance Beamformer and Its Circuit Design for Ultrasound Beamforming
Raghunathan et al. Design of a miniature ultrasound probe for 3D transesophageal echocardiography
JP4130604B2 (ja) 超音波診断装置
Dusa et al. Integrated 16-channel transmit and receive beamforming ASIC for ultrasound imaging
Stuart et al. An architecture and implementation of real-time synthetic aperture compounding with SARUS
JP3886366B2 (ja) 超音波診断装置
Du et al. User parameter free approaches to multistatic adaptive ultrasound imaging

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20170320

Address after: Tokyo, Japan, Japan

Patentee after: Hitachi Ltd.

Address before: Japan Tokyo Chiyoda District Kanda four chome 14 No. 1

Patentee before: Hitachi Medical Corporation

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151125

Termination date: 20200807