CN108700651A - 成像方法、实施所述方法的装置、计算机程序以及计算机可读存储介质 - Google Patents

成像方法、实施所述方法的装置、计算机程序以及计算机可读存储介质 Download PDF

Info

Publication number
CN108700651A
CN108700651A CN201580085761.4A CN201580085761A CN108700651A CN 108700651 A CN108700651 A CN 108700651A CN 201580085761 A CN201580085761 A CN 201580085761A CN 108700651 A CN108700651 A CN 108700651A
Authority
CN
China
Prior art keywords
beam forming
image
wave
weighing vector
pixel
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
CN201580085761.4A
Other languages
English (en)
Other versions
CN108700651B (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.)
SuperSonic Imagine SA
Original Assignee
SuperSonic Imagine SA
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 SuperSonic Imagine SA filed Critical SuperSonic Imagine SA
Publication of CN108700651A publication Critical patent/CN108700651A/zh
Application granted granted Critical
Publication of CN108700651B publication Critical patent/CN108700651B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/52019—Details of transmitters
    • G01S7/5202—Details of transmitters for pulse systems
    • 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/8959—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
    • 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/8997—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques
    • 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/52049—Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
    • 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/52077—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 with means for elimination of unwanted signals, e.g. noise or interference

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Manufacturing Optical Record Carriers (AREA)
  • Holo Graphy (AREA)
  • Debugging And Monitoring (AREA)
  • Traffic Control Systems (AREA)

Abstract

适用于通过换能器阵列(2)来生成介质内的区域的图像的成像方法，包括在介质内传输多个波的传输步骤(201)、用于采集数据集的接收步骤(202)、根据不同的发送加权向量来提供多个波束成形像素值的波束成形步骤(203)，以及将波束成形像素值组合到图像中每个像素的像素值中的组合步骤(204)。发送加权向量(WTn)互不相同且彼此相互正交。

Description

成像方法、实施所述方法的装置、计算机程序以及计算机可读 存储介质
技术领域
本发明涉及成像方法以及实施所述方法的装置，尤其是用于医学成像的装置。
背景技术
本发明更准确地涉及适用于在介质内生成区域图像的成像方法，更准确地说，超声成像方法。
所述方法通过连接着与所述介质有关的换能器阵列的处理单元来实施。
由响应入射波的反射波所产生的这种类型的图像通常包括斑点噪声，所述斑点噪声与所讨论区域的反射波所产生的随机波动相对应。该斑点噪声会导致图像判读的困难。
斑点噪声可以通过使多个图像平均化来减少，因为每个图像的斑点噪声会不同于所述多个图像之中的另一个图像的斑点噪声。
最后，适用于平均化的图像可以通过从不同角度观察区域(空间复合)或者通过改变信号频率(频率复合)来获得。
最终，为了获得平均化的图像，可以将多个不同的接收加权向量应用于接收波束成形。
然而，这些方法都增加了要进行采集的数量，降低了帧速率，并且有时降低了生成图像的分辨率。
发明内容
本发明的一个目的旨在提供适用于在介质内产生区域图像的另一种成像方法，该方法可减少斑点噪声且没有现有技术方法的缺点。
为此目的，该方法包括以下步骤：
(a)传输步骤，其中，第一个复数波通过换能器在介质中传输，
(b)接收步骤，其中，数据集通过对波响应的所述换能器来采集，
(c)波束成形步骤，其中，数据集通过第二个多个波束成形处理来处理，波束成形处理用于提供至少一部分图像的波束成形像素值，每个波束成形处理或者使用与通过发送加权向量所生成的波相对应的数据集或者使用波束成形像素值运算中的发送加权向量，以及，
(d)组合步骤，其中，组合所述第二个复数波束成形像素值，以便提供图像内像素的像素值，而且其中，发送加权向量是互不相同的且彼此相互正交。
得益于上述这些特征，每个发送加权向量都产生一个不相关的斑点噪声，加权数据的组合使之能够计算斑点噪声减少的区域的图像。
在成像方法的各个不同实施例中，可选择性地包含以下特征中的一项和/或其它项。
根据方法的一方面：
-在传输步骤(a)过程中，每个发送加权向量都用于生成一个波，以及，
-在波束成形步骤(c)过程中，波束成形处理是常规的波束成形，其中，波束成形的像素值可由以下公式计算：
式中：
DS(k，l，m)是数据集的矩阵，
WR是接收加权向量，
τ(x，z，l)是适合于波束成形处理的延迟函数，
l是阵列2中换能器的指数，介于l1至l2之间，构成线波束成形的孔径，以及，
x，z是图像内像素的坐标。
根据方法的一方面，在波束成形步骤(c)中应用发送加权向量，以及，
-在波束成形步骤(c)过程中，波束成形处理是合成波束成形，其中，波束成形的像素值可由以下公式计算：
式中：
DS(k，l，m)是数据集的矩阵，
WR是接收加权向量，
τ(x，z，l)是适合于波束成形处理的延迟函数，
l是阵列2中换能器的指数，介于l1至l2之间，构成线波束成形的孔径，
m是指数，介于1至M之间，M是第一个复数，即介质中传输波的数量，
WTn是发送加权向量，
x，z是图像内像素的坐标。
根据方法的一方面，波束成形处理选自包括合成孔径聚焦技术(SAFT)波束成形处理、虚拟换能器SAFT波束成形处理、空间编码SAFT波束成形处理、圆形波合成波束成形处理、平面波合成波束成形处理的列表中。
根据方法的一方面，图像的像素值在组合步骤(d)过程中可由以下公式计算：
式中：
i是复单位虚数，
n是介于1至N之间的指数，N是第二个复数，即发送加权向量的数量，
|X|是X的模数，
HT{X}是X的希尔伯特变换，
Un(x，z)是所述第二个复数的波束成形像素值。
根据方法的一方面，发送加权向量可由选自包括Riedel-Sidorenko函数、离散椭球函数和哈达玛函数的列表中的正交函数来确定。
根据方法的一方面，换能器是发送或接收超声波的超声换能器，方法生成介质内的区域的超声图像。
根据方法的一方面，进一步包括以下步骤：
-初始成像步骤，其中，区域的第一个图像可由所述处理单元和阵列来确定，
-评估步骤，其中，度量值根据第一个图像的像素来确定，
-成像步骤，其中，如果度量值包含在第一个范围中，则不使用发送加权向量来计算图像中的像素值，如果度量值包含在与所述第一个范围不同的第二个范围中，则使用发送加权向量来计算图像中的像素值，所述发送加权向量不相同并且彼此相互正交。
根据方法的一方面，确定度量值，以便区分出与显著反射信号相对应的第一个图像中的像素和与非显著散斑信号相对应的第一个图像中的像素。
根据方法的一方面，度量值通过自相关函数的运算来确定。
根据方法的一方面，度量值是关于介于10度至30度之间所滞后的自相关函数的平均值。
本发明的另一个目的旨在提供实施所述成像方法的装置。用于生成介质内的区域的图像的所述装置包括连接着与所述介质有关的换能器阵列的处理单元，以及，
其中，阵列和处理单元实施以下步骤：
(a)传输步骤，其中，第一个复数(M)波通过换能器在介质内传输，
(b)接收步骤，其中，数据集通过对波响应的所述换能器来采集，以及，
其中，处理单元实施以下步骤：
(c)波束成形步骤，其中，数据集通过第二个复数个波束成形处理来处理，所述波束成形处理用于提供至少一部分图像的波束成形像素值，每个波束成形处理或者使用与通过发送加权向量所生成的波相对应的数据集或者使用波束成形像素值运算中的发送加权向量，以及，
(d)组合步骤，其中，组合所述第二个复数个波束成形像素值，以便提供图像中每个像素的像素值，以及其中，发送加权向量不相同并且彼此相互正交。
在装置的各个实施例中，可选择性地包含以下特征中的一项和/或其它项。
根据装置的一方面，在波束成形步骤(c)过程中应用发送加权向量，以及，
-在波束成形步骤(c)过程中，波束成形处理是合成波束成形，其中，波束成形像素值可由以下公式计算：
式中：
DS(k，l，m)是数据集的矩阵，
WR是接收加权向量，
τ(x，z，l)是适合于波束成形处理的延迟函数，
l是阵列2中换能器的指数，介于l1至l2之间，构成一条线波束成形的孔径，
m是介于1至M之间的指数，M是第一个复数，即介质内传输波的数量，
WTn是发送加权向量，
x，z是图像内像素的坐标。
根据装置的一方面，合成波束成形处理选自包括合成孔径聚焦技术(SAFT)波束成形处理、虚拟换能器SAFT波束成形处理、空间编码SAFT波束成形处理、圆形波合成波束成形处理、平面波合成波束成形处理的列表中。
根据装置的一方面，进一步实施以下步骤：
-初始成像步骤，其中，区域的第一个图像可由所述处理单元和阵列来确定，以及，
-评估步骤，其中，度量值根据第一个图像的像素来确定，
-成像步骤，其中，如果度量值包含在第一个范围中，则不使用发送加权向量来计算图像中的像素值，如果度量值包含在与所述第一个范围不同的第二个范围中，则使用发送加权向量来计算图像中的像素值，所述发送加权向量不相同并且彼此相互正交。
本发明的另一个目的旨在提供一种计算机程序，在计算机执行所述程序时，所述计算机程序包括实施上述成像方法的步骤的指令。
本发明的另一个目的旨在提供一种计算机可读存储介质，所述计算机程序储存于所述计算机可读存储介质，且在计算机执行所述程序时，所述计算机程序包括执行上述成像方法的步骤的指令。
附图说明
本发明的其它特征和优点将通过以非限制性示例方式所列出的本发明两个实施例的以下详细说明及其参考附图而更为显而易见。在附图中：
-图1是显示了根据本发明一个实施例的超声成像装置的示意图；
-图2是显示了图1所示装置的一部分的方框图；
-图3是显示了根据本发明通过图1所示装置来实施成像方法的第一个示例的图表；
-图4是显示了根据本发明通过图1所示装置来实施成像方法的第二个示例的图表；
-图5a至图5c是彼此相互正交的三个发送加权向量的第一个示例，所述向量与Riedel-Sidorenko函数相对应；
-图6a至图6c是彼此相互正交的三个发送加权向量的第二个示例，所述向量与离散椭球函数相对应；
-图7a至图7c是彼此相互正交的三个发送加权向量的第三个示例，所述向量与哈达玛函数相对应；
-图8是根据本发明的成像方法的第三个示例，所述第三个示例是自适应方法；
-图9显示了针对一个样本计算且与图8所示方法中的度量评价决策相对应的决策图像的示例；
-图10显示了针对与图9中相同的样本计算的图像的示例；以及，
-图11显示了在图8的成像方法的第三个示例中使用的自相关函数的示例。
具体实施方式
在各图中，相同标号指代相同或相似元素。
图1所示的装置适合于区域1的成像，所述区域1例如是活体组织，尤其是患者的人体组织。装置可包括，例如：
-换能器阵列2，例如，通常包括数十个(例如，100至300个)换能器的线性阵列，所述换能器沿着轴X(水平或者阵列方向X)并列，如同使区域1成像的常用探头中已知的那样(因此，阵列2适合执行区域1的二维(2D)成像，但是，阵列2也可以是适合执行区域1的3D成像的二维阵列)；
-电子舱3，用于控制换能器阵列并采集其信号；
-微型计算机4，用于控制电子舱3并查看由电子舱所获得的图像(在一个变体中，单个电子装置可以完成电子舱3和微型计算机4的所有功能)。
图1上的轴Z是垂直于轴X的轴，通常指向由阵列换能器生成的超声束的方向。在本文件中，该方向是指垂直或轴向方向。
换能器阵列2也可以是包括沿着曲线对齐的多个换能器的凸阵。
如图2所示，例如，电子舱3可包括：
-L个模拟/数字转换器5(A/D1-A/DL)，分别连接着换能器阵列2的L个换能器(T1-TL)；
-L个缓冲存储器6(B1-Bn)，分别连接着n个模拟/数字转换器5，
-与缓冲存储器6和微型计算机4相通讯的中央处理单元8(CPU)，
-连接着中央处理单元8的存储器9(MEM)；
-连接着中央处理单元8的数字信号处理器10(DSP)。
在此公开的装置是用于超声成像的装置，换能器是超声换能器，所实施的方法适用于生成区域1的超声图像。
然而，装置可以是利用除超声波以外其它波(波长不同于超声波长的波)的任何成像装置，换能器和电子舱组件可用于所述波。
图3和图4显示了通过图1和图2所示装置实施所述方法的两个示例。方法的步骤主要通过中央处理单元8并最终籍助于数字信号处理器10或者任何其它装置来控制。
方法包括以下主要步骤：
(a)传输步骤(101；201)，其中，通过换能器在介质的区域内传输第一个复数波；
(b)接收步骤(102；202)，其中，通过对波响应的所述换能器来采集数据集；
(c)波束成形步骤(103；203)，其中，通过第二个复数波束成形处理来处理在接收步骤所采集到的数据集，以便提供适用于至少一部分图像的第二个复数波束成形像素值Un(x，z)，以及，
(d)组合步骤(104；204)，其中，组合N个波束成形像素值，以便提供图像中每个像素的像素值。
第一个复数是为了处理图像在介质的区域内所传输的波的数量。即换能器为了生成所述波而连续发射的数量。
第二个复数是方法中所使用的发送加权向量WTn的数量。
根据本发明的方法采用第二个复数(数量N)发送加权向量WTn，所述发送加权向量彼此相互不同。
此外，在方法中所使用的发送加权向量WTn彼此相互正交，换言之：
指数i，j属于1…N，
指数i与指数j不同，则WTi.WTj＝0。
即：
发送加权向量WTn是包括数量L个分量的向量，每个分量都与施加的信号的放大系数相对应，所述信号通常在传输步骤(a)过程中发送到阵列2的换能器。发送加权向量WTn的分量可定义为小于或等于1的值：
WTn(l)<＝1，当l＝1至L。
根据图3所阐释的第一个变体，在传输步骤a)(101)过程中将发送加权向量应用于传输波，从而提供包括与多个发送加权向量相对应的数据的数据集。
然后，在波束成形步骤c)过程中使用该数据集，这个步骤的每个波束成形处理都是用与特定发送加权向量(特定传输波)相对应的数据。
在该第一个变化例中，由于每个发送加权向量的加权效应，所以在介质中传输的振幅波可以较小，因此，所接收到的数据集的信噪比可以较低。此外，该第一个变体需要对波进行MxN次传输和采集，所以，也降低了成像帧频。而且，其生成的图像具有斑点噪声且在整个扫描区域中都不是最优的，因为它只适用于波束的焦点深度。
根据图4所示的第二个变化例中，发送加权向量都应用在波束成形步骤(c)(203)过程中的各个波束成形处理计算中的数据集(未加权的)。
在那种情况下，第一个变体的加权效应可通过在波束成形步骤(c)(203)过程中的适当的计算(求和)来获得，所述计算利用各个波束成形处理的发送-接收波束成形公式。
在该第二个变体中，波通过单位放大系数(与加窗或孔径相关的系数除外)在介质内传输。传输加权系数仅通过计算在接收步骤(b)中所采集到的数据集来应用。
得益于该发送-接收波束成形处理，真正在介质内传输的波的振幅没有减少，而且保存了数字化为数据集的所接收到的信号的信噪比(SNR)。
在方法的所有变体中，将换能器阵列2与待成像的介质(例如，患者的身体)相接触。
传输波的数量可例如介于2至100之间。超声波的频率可例如介于0.5至100MHz之间，例如，介于1至10MHz之间。
发送加权向量WTn的数量N是例如较小的数，例如，介于2至10之间，或，例如，N＝3。可有很多正交函数用于确定该发送加权向量WTn。
图5a至图5c显示了3个发送加权向量WTn的第一个示例。在这些图中，横坐标是换能器指数，纵坐标是各个向量中分量的值。
这些发送加权向量WT1-3是彼此相互正交的Riedel-Sidorenko函数的向量。
图6a至图6c显示了采用离散椭球函数的3个发送加权向量WTn的第二个示例，也是彼此相互正交的。
图7a至图7c显示了采用哈达玛函数的3个发送加权向量WTn的第二个示例，也是彼此相互正交的。
发送加权向量的数量N最好小于或等于传输到介质中的波的数量(第一个复数)。
关于本发明效果的演示
由于采用正交发送加权向量，使图像中的斑点噪声光滑，因为每个传输向量都成产生互不相关的散斑。
可以用数学方法证明这一点。在这部分，我们采用1991年11月《美国声学学会杂志》第90(5)期中由Raoul Mallart和Mathias Fink编著的文件《脉冲回波测量中的vanCittert-Zernike定理》中所用的形式和符号。
在点X1及频率f下的入射压力场可由下式给出：
式中：
O(X)是聚焦孔径的发送孔径函数；以及，
X＝(x，y，0)，其中，x，y是图像中的坐标。
这部分的目的是表明正交传输孔径Ok(X)和Ol(X)产生互不相关的散斑图。由位于再点X1的个别散射体散射的压力场是球面波。
在点X0的接收压力场表达为：
式中：r01＝|X0-X1|。
假设散射介质是不连贯的，换言之，其结构是随机的并且比成像系统所使用的最小波长更细。介质就不能分辨出，其散射函数RM的自相关形式如下：
RM(X1-X2，f)＝χ(X，f)δ(X1-X2)
式中：χ是点坐标(x，z)附近的局部散射系数。
在点X0感测到的响应狄拉克脉冲由整个介质反向散射的压力场可由下式给出：
压力场Pk(X0，f)和P1(X0，f)的互相关Rk1(X0，f)可由下式给出：
Rkl(X0，f)＝E{Pk(X0，f)Pl(X0，f)}
式中：
E{.}表示数学期望值。
可以通过上述方程式把乘积Pk(X0，f).P1(X0，f)表达为：
通过上述方程式，可以推导出：
于是，压力场Pk(X0，f)和P1(X0，f)的互相关可由下式给出：
通过上文引用的Mallart等人的出版物，在点X1的入射压力场大约为：
式中：是相位项。
关于在深度z聚焦的孔径，孔径函数O(x)包含相位项因此是实值。
使则入射场的乘积可表达为：
而且，将最后一个方程注入前一个方程中，可以推导出：
上述方程式表明，正交孔径，即使之从而获得的互不相关的采集数据，即
因此，压力场的互相关为空值，而且斑点噪声是互不相关的。
所以，与现有技术的方法相比，本发明的方法可以减少斑点噪声。实际上，根据所述方法的斑点噪声是最小的。
本发明的方法可应用于任何超声成像方法。以下说明将对其中几项进行解释。
示例1：常规成像
现在对第一个示例进行说明，所述第一个示例与传统成像方法相对应，通常称为常规的聚焦孔径或传统的B型成像。
在传输步骤(a)过程中，通过换能器将数量M个波连续传输到区域中。这些波都是按照大致垂直于阵列方向(X)的垂直方向(Z)在介质内传输的聚焦波束，并且聚焦于在与换能器阵列2相分隔开的焦距(或者焦点深度)处。
M个连续聚焦波束按照侧向彼此横向移动，所述侧向与阵列方向(X)一致，以便扫描区域。
根据本发明，每个聚焦波束都要重复N次，换言之，一个聚焦波束适用于各个指数n的发送加权向量WTn聚焦波束。将发送加权向量WTn直接应用于换能器：将发送加权向量的指数l的各个分量用于放大或减小发送到阵列2的指数l的换能器的信号。
在该第一个示例中，第一个复数(在介质内传输的波的数量)因此等于MxN。
每个波都在区域中传播，在此与扩散粒子相互作用，所述扩散粒子是超声波的反射。于是，波将作为朝向换能器阵列2的反射波(包括回声)的反向散射。
在接收步骤(b)中，如果各个反射波能被换能器所接收且通过模拟数字转换器采集和转换为数据并储存于存储器中。然后，将所有反射波都以数据集或数据组的方式储存于存储器中。
将来自采集波的数据集组成矩阵DS(k，l，m)，式中：
k是超时样本的指数，
l是阵列中换能器的指数，
m是数量MxN个发射或传输波(数量与第一个复数相一致)中的传输波的指数。
在波束成形步骤(c)过程中，通过波束成形处理来计算图像的一行或几行(垂直行或水平行)。这些行都平行于垂直或轴向方向Z(聚焦波束的方向)并且包含在聚焦波束内。
对于N个发送加权向量WTn中的各个而言，行内像素的波束成形像素值Un(x,z)都可由以下与接收波束成形相对应的波束成形公式来计算：
式中：
DS(k,l,m)是储存所有数据集的矩阵，
WR是接收加权向量，
τ(x,z,l)是波束成形处理的延迟函数，即与当前传统的波束成形处理相一致，
l是阵列2中换能器的指数，介于l1至l2之间，构成一行波束成形的孔径，
m是指向指数n的唯一或确定发送加权向量WTn的指数并且与在介质内坐标为(x，z)点附近聚焦的波束相对应，以及，
x，z是图像被像素的坐标。
在组合步骤(d)中，N(第二个复数)个波束成形像素值Un(x,z)可根据各个发送加权向量来计算，并将其组合，以便提供图像内各个像素的像素值I(x,z)。
然后，像素值I(x,z)可由以下求和公式来计算：
式中：
i是复单位虚数。
n是介于1至N之间的指数，N是第二个复数中发送加权向量的数量，
|X|是X的模数，
HT{X}是X的希尔伯特变换，
Un(x,z)是指数n的波束成形像素值。
可惜，上述方法：
-需要MxN次传输和采集波，因此，步骤a)和b)要花费一些时间，
-生成的图像具有斑点噪声且在整个扫面区域中不是最优的，因为它只适用于波束的焦点深度。
示例2合成成像
现在对与合成成像方法相对应的第二个示例进行说明。
作为合成成像方法，应理解为称之为合成方法的任何方法，以及至少以下列表中的所有方法：
1)合成孔径聚焦技术方法(SAFT)，
2)虚拟换能器SAFT方法，
3)空间编码SAFT方法，
4)圆形波合成方法，以及，
5)平面波合成方法。
1)SAFT方法，例如，在2006年的《超声波学》44，e5-e15中J.A.Jensen，S.I.Nikolov，K.L.Gammelmark，M.H.Pedersen的发表文件《合成孔径超声成像》中所详细说明的。
该方法执行：
-传输步骤(a)，其中，至少有一次波的传输是根据阵列2的各个换能器来进行(发射)的：阵列2的各个换能器是一个接一个激发的，以及，
-接收步骤(b)，其中，采集所有换能器的信号，并将其作为数据集记录(储存)在存储器中。
也可以将数据集组合成矩阵DS(k，l，m)，通常称为“完整数据集”，式中：
k是超时样本的指数，
l是阵列中换能器的指数，
m是在发射波数量(与第一个复数相一致的数量)中的传输波的指数。
因此，发射波的数量M通常等于阵列2中换能器的数量L。但是，如果没有使用阵列的所有换能器，则该数量可小于换能器的数量L，或者，如果使其部分作平均化，则可大于换能器的数量L。
于是，波束成形步骤(c)和组合步骤(d)不同于已知的SAFT方法，在下文对其进行说明。
在波束成形步骤(c)过程中，波束成形像素值Un(x,z)可由包括双重求和(一个根据所采用的具有指数l的换能器的数量，另一个根据具有指数m的超声波的数量)的以下波束成形公式来计算：
式中：
DS(k,l,m)是数据集的矩阵，
WR是接收加权向量，
τ(x,z,l)是波束成形处理的延迟函数，即与当前SAFT波束成形处理相对应，
l是阵列2中换能器的指数，介于l1至l2之间，构成线波束成形的孔径，
m是介于1至M之间的指数，M是第一个复数，即介质内传输波的数量，
WTn是指数n的发送加权向量，
x，z是图像内像素的坐标。
在组合步骤(d)过程中，根据数量N(第二个复数)个发送加权向量所计算的波束成形像素值Un(x,z)也进行组合，以便提供图像内每个像素的像素值I(x,z)。
然后，像素值I(x,z)可由上文所述的相同求和公式来计算：
发送加权向量WT通过在波束成形步骤(c)中的计算来应用，该方法无需额外发射(将波发射到介质中)。
通过SAFT方法生成的图像通常在整个扫描区域内具有较高质量：高分辨率及高对比度。信噪比(SNR)不是最佳的，因为每次发射仅使用一个转移的阵列，而且，与其它方法相比，发射能量较低。
得益于正交发送加权向量WTn，上述改进的SAFT方法能够减少图像的斑点噪声。
2)虚拟换能器SAFT方法在信噪比方面改进了SAFT方法。该方法，例如，在2008年IEEE国际超声研讨会会议记录第966-969页，J.Kortbek，J.A.Jensen，K.L.Gammelmark的发表文件《合成孔径的连续波束成形》中公开了。
与以前的SAFT方法相比，发射波并非仅通过阵列的一个换能器来生成，而是由阵列的多个换能器来生成的，使得发射波是聚焦在区域内预定聚焦区的聚焦波束。
于是，在该虚拟换能器SAFT方法中：
-在传输步骤(a)过程中，通过多个换能器来传输分别与区域内聚焦区的聚焦波束相对应的多个发射波，以及，
-在接收步骤(b)过程中，数据集可根据多个聚焦区来采集并且数据集可组成矩阵DS(k，l，m)，通常称之为“完整数据集”，与SAFT方法相等效。
然后，将本发明的波束成形步骤(c)和组合步骤(d)应用于虚拟换能器SAFT方法，并且与上述改进的SAFT方法所应用的方法相同或相似。
得益于该改进的虚拟换能器SAFT方法，利用了第二个复数个发送加权向量WTn，所生成的图像的斑点噪声减少了。
3)空间编码SAFT方法也在信噪比方面改进了SAFT方法。
在该方法中，
-在传输步骤(a)过程中，波是通过将传输矩阵TM应用于SAFT方法的换能器信号来发射：对于M个发射波(第一个复数)中的各个波而言，将换能器的信号乘以传输矩阵TM，所述传输矩阵是可逆的矩阵，以及，
-在接收步骤(b)响应发射波所采集到的数据集储存于接收矩阵RM(k，l，m)中，而且数据集的矩阵DS(k，l，m)可由以下反演公式利用接收矩阵RM和传输矩阵TM来获得：
DS(k,l,m)|k＝constant＝TM-1·RM(k,l,m)|k＝constant，
关于k＝1至K。
数据集的矩阵DS相同于上文在SAFT方法中所述的。
然后，将本发明的波束成形步骤(c)和组合步骤(d)应用于空间编码SAFT方法，并且与上述改进的SAFT方法所应用的方法相似。
得益于该改进的空间编码SAFT方法，利用了发送加权向量WTn的第二个复数，所生成的图像的斑点噪声减少了。
4)圆形波合成方法也在信噪比方面改进了SAFT方法。也与虚拟换能器SAFT方法相似，但区别之处在于聚焦波束聚焦于阵列之后，使得圆形波可在介质内传输。
该方法的详细信息可见于2009年IEEE超声波研讨会会议记录第515-518页，M.Couade，M.Pernot，M.Tanter，E.Messas，A.Bel，M.Ba，A.-A.Hagege，M.Fink的发表文章《利用相控阵列的圆形波合成成像进行心脏超快成像》。
然后，将本发明的波束成形步骤(c)和组合步骤(d)应用于该圆形波合成方法，并且与上述改进的SAFT方法所应用的方法相同或相似。
得益于该改进的圆形波合成方法，利用第二个复数个发送加权向量WTn，所生成的图像的斑点噪声减少了。
5)平面波合成方法也在信噪比方面改进了SAFT方法。
该方法的详细信息可见于已发布专利US6，551，246或已发布专利申请书US2009/0234230。
在该方法中：
-在传输步骤(a)过程中，将第一个复数(M)平面波发射到介质中，以及，
-在接收步骤(b)过程中，采集换能器信号并将其作为数据集的矩阵DS(k，l，m)记录(储存)于存储器，m是传输到介质中的平面波的指数。
因此，该方法与SAFT方法的不同之处在于在介质内传输(射出、发射)数量M个平面波。
然后，将本发明的波束成形步骤(c)和组合步骤(d)应用于该平面波合成方法，而且与上述改进的SAFT方法所应用的方法相同或相似。
在波束成形步骤(c)中，波束成形像素值Un(x,z)可由以下波束成形公式来计算：
按照平面波合成方法，各个平面波可通过不同的发送加权向量WTn来加权。
在以上所有波束成形处理中，延迟函数τ(x,z,l,m)是取决于各种类型的波束成形处理以及各种类型的探针(形状、尺寸)所众所周知的函数。
根据图8阐释的成像方法的第三个变体，成像方法没有把多个发送加权向量用于图像中的所有像素。由于这一特征，所以没有降低这些像素的横向分辨率(没有组合发送加权向量)而且减少了图像中其它像素的斑点噪声(具有发送加权向量组合)。
在该第三个变体中，成像方法进一步包括以下步骤。
-初始成像步骤(301)，其中，区域的第一个图像可由所述处理单元和阵列来确定，
-评估步骤(302)，其中，度量值根据第一个图像内的像素来确定，
-成像步骤(303，304，305)，其中，如果度量值包含在第一个范围内，则图像中的像素值可以直接计算而不使用多个发送加权向量，如果度量值包含在与所述第一个范围不同的第二个范围中，则图像中的像素值要使用多个(第二个复数，N个)发送加权向量来计算，所述发送加权向量互不相同且按照(第二个)复数彼此相互正交。
然后，成像方法通过组合在成像步骤(303，304，305)所计算的像素就可以建立最终图像307。根据上文所述的方法之一，或者在步骤304不使用发送加权向量，或者在步骤305中使用发送加权向量，来计算最终图像得像素值。
有利的是，确定度量值，以便区分出与显著反射信号相对应的第一个图像中的像素和与非显著散斑信号相对应的第一个图像中的像素。来自第一个图像中像素的显著反射信号通常是与具有强反射体要素的介质内的位置相对应的像素。在那种情况下，成像方法不使用降低横向分辨率的发送加权向量。来自第一个图像中像素的非显著散斑信号通常是与不具有强反射体要素的介质内的位置相对应的像素。该位置与图像中的斑点噪声的位置相对应。在那种情况下，成像方法可采用减少所述斑点噪声的发送加权向量，所述这些斑点噪声对用户来说既难以令人满意且又不重要。
图9和图10是通过成像方法的第三个变体所生成的图像的示例。图9是二值图像303，其中，黑色像素是没有传输加权(方法的步骤304)的像素，白色像素是具有传输加权(方法的步骤305)的像素。图10是通过第三个变体的成像方法300所提供的最终图像。
得益于上述混合的自适应方法，该图像307的斑点噪声减少了且具有良好的横向分辨率。
度量值可由自相关函数R(θ，x，z)的运算来确定。
例如，自相关函数可由下式确定：
R(θ,x,z)＝E{Im1(x,z).Im2(x,z)}
式中：E{}是数学期望值，以及，
Im(x，z)＝|Vm(x，z)+iHT{Vm(x，z)}|2
DS(k，l，m)是数据集的矩阵，
k是超时指数，
l是阵列中换能器的指数，
m是关于合成波束成形的传输波(例如，平面波)的指数，
WR是接收加权向量，
τ(x，z，l，m)是波束成形处理的延迟函数，例如，与平面波的波束成形处理相对应。
滞后θ是第一个角度与第二个角度之间相差的角度，所述第一个角度与指数m1的第一次发射相对应，第二个角度与指数m2的第二次发射相对应。第一次发射和第二次发射都是平面波发射，第一个角度和第二个角度是与换能器阵列相关的所述平面波的角度。
图11显示了这种自相关函数的两个示例：在与强反射体相对应的介质内的位置创建的第一个自相关函数的曲线401，以及在与斑点噪声相对应的介质内的位置创建的第二个自相关函数402。
两个自相关函数401，402的区别是10度至30度的滞后范围。该差异可用于区分第一个图像中的两种类型的像素，以及选择(步骤303)是否使用(步骤304，305)发送加权向量来计算成像方法的最终图像的像素。

Claims (17)

1.适用于生成介质内的区域的图像的成像方法，其中，所述方法通过连接着换能器阵列(2)的处理单元(8)来实施，所述换能器阵列与所述介质有关，所述方法包括以下步骤：
(a)传输步骤(101；201)，其中，第一个复数(M)波通过换能器在介质中传输，
(b)接收步骤(102；202)，其中，数据集通过对波响应的所述换能器来采集，
(c)波束成形步骤(103；203)，其中，数据集通过第二个复数波束成形处理来处理，所述波束成形处理用于提供至少一部分图像的波束成形像素值(Un(x，z))，每个波束成形处理或者使用与通过发送加权向量生(WTn)成的波相对应的数据集或者使用波束成形像素值运算中的发送加权向量(WTn)，以及，
(d)组合步骤(104；204)，其中，组合所述第二个复数个波束成形像素值，以便提供图像内像素的像素值(I(x，z))，而且其中，发送加权向量(WTn)互不相同且彼此相互正交。
2.根据权利要求1所述的方法，其特征在于：
-在传输步骤(a)过程中，每个发送加权向量(WTn)都用于生成一个波，以及，
-在波束成形步骤(c)过程中，波束成形处理是常规的波束成形处理，其中，波束成形的像素值(Un(x，z))可由以下公式计算：
式中：
DS(k，l，m)是数据集的矩阵，
WR是接收加权向量，
τ(x，z，l)是适合于波束成形处理的延迟函数，
l是阵列2中换能器的指数，介于l1至l2之间，构成线波束成形的孔径，以及，
x，z是图像内像素的坐标。
3.根据权利要求1所述的方法，其特征在于，所述发送加权向量(WTn)应用于波束成形步骤(c)中，以及，
-在波束成形步骤(c)过程中，波束成形处理是合成波束成形，其中，波束成形的像素值(Un(x，z))可由以下公式计算：
式中：
DS(k，l，m)是数据集的矩阵，
WR是接收加权向量，
τ(x，z，l)是适合于波束成形处理的延迟函数，
l是阵列2中换能器的指数，介于l1至l2之间，构成线波束成形的孔径，
m是指数，介于1至M之间，M是第一个复数，即介质中传输波的数量，
WTn是发送加权向量，
x，z是图像内像素的坐标。
4.根据权利要求3所述的方法，其特征在于，所述合成波束成形处理选自包括合成孔径聚焦技术(SAFT)波束成形处理、虚拟换能器SAFT波束成形处理、空间编码SAFT波束成形处理、圆形波合成波束成形处理、平面波合成波束成形处理的列表中。
5.根据权利要求1至4中任一项所述的方法，在组合步骤(d)过程中，图像的像素值(I(x，z))可由以下公式计算：
式中：
i是复单位虚数。
n是介于1至N之间的指数，N是第二个复数，即发送加权向量的数量，
|X|是X的模数，
HT{X}是X的希尔伯特变换，
Un(x，z)是所述第二个复数的波束成形像素值。
6.根据权利要求1至5中任一项所述的方法，其特征在于，所述发送加权向量(WTn)通过选自包括Riedel-Sidorenko函数、离散椭球函数和哈达玛函数的列表中的正交函数来确定。
7.根据权利要求1至6中任一项所述的方法，其特征在于，所述换能器是适用于发送或接收超声波的超声换能器，以及所述方法生成在介质内的区域的超声图像。
8.根据权利要求3至7中任一项所述的成像方法，进一步包括以下步骤：
-初始成像步骤(301)其中，所述区域的第一个图像通过所述处理单元和阵列来确定，
-评估步骤(302)，其中，度量值根据第一个图像的像素来确定，
-成像步骤(303，304，305)，其中，如果度量值包含在第一个范围中，则不使用发送加权向量来计算图像中的像素值，如果度量值包含在与所述第一个范围不同的第二个范围中，则使用发送加权向量来计算图像中的像素值，所述发送加权向量互不相同且彼此相互正交。
9.根据权利要求8所述的成像方法，其特征在于，确定度量值，以便区分与显著反射信号相对应的第一个图像中的像素和与非显著散斑信号相对应的第一个图像中的像素。
10.根据权利要求8或9所述的成像方法，其特征在于，所述度量值通过自相关函数的运算来确定。
11.根据权利要求10所述的成像方法，其特征在于，所述度量值是关于介于10度至30度之间所滞后的自相关函数的平均值。
12.适用于生成介质内的区域的图像的装置，包括连接着与所述介质有关的换能器阵列(2)的处理单元(8)，以及，
其中，阵列和处理单元实施以下步骤：
(a)传输步骤(101；201)，其中，第一个复数(M)波通过换能器在介质内传输，
(b)接收步骤(102；202)，其中，数据集通过对波响应的所述换能器来采集，以及，
其中，所述处理单元实施以下步骤：
(c)波束成形步骤(103；203)，其中，所述数据集通过第二个复数波束成形处理来处理，所述波束成形处理用于提供至少一部分图像的波束成形像素值(Un(x，z))，每个波束成形处理或者使用与通过发送加权向量(WTn)所生成的波相对应的数据集或者使用波束成形像素值运算中的发送加权向量(WTn)，以及，
(d)组合步骤(104；204)，其中，组合所述第二个复数波束成形像素值，以便提供图像中每个像素的像素值(I(x,z))，以及
其中，所述发送加权向量(WTn)互不相同且彼此相互正交。
13.根据权利要求12所述的装置，其特征在于，所述发送加权向量(WTn)应用于波束成形步骤(c)过程中，以及，
-在波束成形步骤(c)过程中，所述波束成形处理是合成波束成形，其中，波束成形像素值(Un(x，z))可由以下公式计算：
式中：
DS(k，l，m)是数据集的矩阵，
WR是接收加权向量，
τ(x，z，l)是适合于波束成形处理的延迟函数，
l是阵列2中换能器的指数，介于l1至l2之间，构成线波束成形的孔径，m是介于1至M之间的指数，M是第一个复数，即介质内传输波的数量，WTn是发送加权向量，
x，z是图像内像素的坐标。
14.根据权利要求13所述的装置，其特征在于，合成波束成形处理选自包括合成孔径聚焦技术(SAFT)波束成形处理、虚拟换能器SAFT波束成形处理、空间编码SAFT波束成形处理、圆形波合成波束成形处理、平面波合成波束成形处理的列表中。
15.根据权利要求12至14中任一项所述的装置，进一步实施以下步骤：
-初始成像步骤(301)，其中，区域的第一个图像通过所述处理单元和阵列来确定，以及，
-评估步骤(302)，其中，度量值根据第一个图像的像素来确定，
-成像步骤(303，304，305)，其中，如果度量值包含在第一个范围中，则不使用发送加权向量来计算图像中的像素值，如果度量值包含在与所述第一个范围不同的第二个范围中，则使用发送加权向量来计算图像中的像素值，所述发送加权向量互不相同且彼此相互正交。
16.计算机程序，在计算机执行所述程序时，所述计算机程序包括执行根据权利要求1至15中任一项所述方法的步骤的指令。
17.计算机可读存储介质，在所述计算机可读存储介质上储存计算机程序，在计算机执行所述程序时，所述计算机程序包括执行根据权利要求1至15中任一项所述方法的步骤的指令。
CN201580085761.4A 2015-12-01 2015-12-01 成像方法、实施所述方法的装置以及计算机可读存储介质 Active CN108700651B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/IB2015/002652 WO2017093778A1 (en) 2015-12-01 2015-12-01 An imaging method, an apparatus implementing said method, a computer program and a computer-readable storage medium

Publications (2)

Publication Number Publication Date
CN108700651A true CN108700651A (zh) 2018-10-23
CN108700651B CN108700651B (zh) 2022-03-29

Family

ID=56373079

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201580085761.4A Active CN108700651B (zh) 2015-12-01 2015-12-01 成像方法、实施所述方法的装置以及计算机可读存储介质

Country Status (10)

Country Link
US (1) US11487006B2 (zh)
EP (2) EP4235216A3 (zh)
JP (1) JP6734376B2 (zh)
KR (1) KR102452220B1 (zh)
CN (1) CN108700651B (zh)
BR (1) BR112018011282B1 (zh)
CA (1) CA3007284A1 (zh)
ES (1) ES2953404T3 (zh)
IL (1) IL259759B (zh)
WO (1) WO2017093778A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3936891A1 (en) 2020-07-10 2022-01-12 Supersonic Imagine Method and system for estimating an ultrasound attenuation parameter
KR102550262B1 (ko) * 2021-02-22 2023-07-03 광주과학기술원 무작위 간섭을 이용한 초음파 이미징 장치 및 그 방법
GB2624639A (en) * 2022-11-22 2024-05-29 Darkvision Tech Inc Surface texture of manufactured parts
EP4394441A1 (en) 2022-12-29 2024-07-03 Supersonic Imagine A beamforming method and system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6048315A (en) * 1998-09-28 2000-04-11 General Electric Company Method and apparatus for ultrasonic synthetic transmit aperture imaging using orthogonal complementary codes
US6213946B1 (en) * 1998-12-24 2001-04-10 Agilent Technologies, Inc. Methods and apparatus for speckle reduction by orthogonal pulse compounding in medical ultrasound imaging
CN102393520A (zh) * 2011-09-26 2012-03-28 哈尔滨工程大学 基于目标回波多普勒特性的声纳运动目标成像方法
CN102429684A (zh) * 2010-09-28 2012-05-02 深圳迈瑞生物医疗电子股份有限公司 一种多普勒彩色血流成像方法和装置
US20150269417A1 (en) * 2014-03-19 2015-09-24 Samsung Electronics Co., Ltd. Method and apparatus for processing images

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5851187A (en) * 1997-10-01 1998-12-22 General Electric Company Method and apparatus for ultrasonic beamforming with spatially encoded transmits
US6551246B1 (en) 2000-03-06 2003-04-22 Acuson Corporation Method and apparatus for forming medical ultrasound images
JP2009536855A (ja) * 2006-05-12 2009-10-22 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ コヒーレントでない遡及的で動的な送信の焦点合わせ
US7782521B2 (en) * 2007-05-31 2010-08-24 Texas Instruments Incorporated System and method for displaying images
US8641628B2 (en) * 2007-09-26 2014-02-04 Siemens Medical Solutions Usa, Inc. Aperture synthesis using cMUTs
US9117439B2 (en) * 2008-03-13 2015-08-25 Supersonic Imagine Method and apparatus for ultrasound synthetic imagining
US9254116B2 (en) * 2010-04-02 2016-02-09 Duke University Methods, systems and apparatuses for Van-Cittert Zernike imaging
US8743657B1 (en) * 2011-04-22 2014-06-03 The United States Of America As Represented By The Secretary Of The Navy Resolution analysis using vector components of a scattered acoustic intensity field
EP2574956A1 (en) * 2011-09-30 2013-04-03 GE Inspection Technologies Ltd Ultrasound imaging system and method with side lobe suppression via coherency factor weighting
WO2013116807A1 (en) * 2012-02-03 2013-08-08 Los Alamos National Security, Llc Systems and methods for synthetic aperture ultrasound tomography
JP6288998B2 (ja) * 2013-09-17 2018-03-07 キヤノンメディカルシステムズ株式会社 超音波診断装置及び超音波イメージングプログラム
EP3232937B1 (en) * 2014-12-19 2023-09-27 University of Rochester Ultrasound system for high-speed and high resolution imaging applications

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6048315A (en) * 1998-09-28 2000-04-11 General Electric Company Method and apparatus for ultrasonic synthetic transmit aperture imaging using orthogonal complementary codes
US6213946B1 (en) * 1998-12-24 2001-04-10 Agilent Technologies, Inc. Methods and apparatus for speckle reduction by orthogonal pulse compounding in medical ultrasound imaging
CN102429684A (zh) * 2010-09-28 2012-05-02 深圳迈瑞生物医疗电子股份有限公司 一种多普勒彩色血流成像方法和装置
CN102393520A (zh) * 2011-09-26 2012-03-28 哈尔滨工程大学 基于目标回波多普勒特性的声纳运动目标成像方法
US20150269417A1 (en) * 2014-03-19 2015-09-24 Samsung Electronics Co., Ltd. Method and apparatus for processing images

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHANGHAN YOON ET AL.: "Orthogonal Quadratic Chirp Signals for Simultaneous Multi-Zone Focusing in Medical Ultrasound Imaging", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL》 *
邹强鑫: "超声数字波束合成方法的研究与优化设计", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *

Also Published As

Publication number Publication date
US20200256988A1 (en) 2020-08-13
EP4235216A3 (en) 2023-10-25
EP4235216A2 (en) 2023-08-30
ES2953404T3 (es) 2023-11-10
WO2017093778A1 (en) 2017-06-08
JP2018535789A (ja) 2018-12-06
KR102452220B1 (ko) 2022-10-07
BR112018011282B1 (pt) 2024-01-30
BR112018011282A2 (pt) 2019-01-22
US11487006B2 (en) 2022-11-01
CN108700651B (zh) 2022-03-29
EP3384313B1 (en) 2023-07-19
EP3384313A1 (en) 2018-10-10
JP6734376B2 (ja) 2020-08-05
IL259759A (en) 2018-07-31
CA3007284A1 (en) 2017-06-08
IL259759B (en) 2022-08-01
KR20180113501A (ko) 2018-10-16

Similar Documents

Publication Publication Date Title
Kortbek et al. Sequential beamforming for synthetic aperture imaging
US8317712B2 (en) Retrospective dynamic transmit focusing for spatial compounding
JP6408297B2 (ja) ビームフォーミング方法、計測イメージング装置、及び、通信装置
US20130258805A1 (en) Methods and systems for producing compounded ultrasound images
US20100004540A1 (en) Dual path processing for optimal speckle tracking
Tasinkevych et al. Modified synthetic transmit aperture algorithm for ultrasound imaging
CN109363714A (zh) 一种超声成像设备及其超声成像方法
CN108700651A (zh) 成像方法、实施所述方法的装置、计算机程序以及计算机可读存储介质
Zhang et al. Ultrafast ultrasound imaging using combined transmissions with cross-coherence-based reconstruction
KR20220036339A (ko) 매질의 초음파 특성화 방법 및 시스템
Hemmsen et al. Tissue harmonic synthetic aperture ultrasound imaging
KR20220036338A (ko) 매질의 초음파 특성화 방법 및 시스템
WO2021009495A1 (en) Method and apparatus for adaptive beamforming
US8235906B2 (en) System and method for accelerated focused ultrasound imaging
CN108024795A (zh) 频谱解析装置、频谱解析方法以及超声波摄像装置
KR20230088284A (ko) 빔형성된 데이터를 처리하는 방법 및 시스템
CN114027873B (zh) 一种超声成像方法、装置及计算机可读存储介质
CN115227288B (zh) 一种基于滞后的延时乘累加超声虚源波束形成方法
EP4394441A1 (en) A beamforming method and system
CN115607185A (zh) 超声成像方法和超声成像系统
CN115844445A (zh) 用于处理由换能器元件接收的信号组的方法和系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: French Aiang Provence

Applicant after: Sound imaging Co.,Ltd.

Address before: French Aiang Provence

Applicant before: SUPER SONIC IMAGINE

GR01 Patent grant
GR01 Patent grant