WO2024027529A1 - 一种抑制杂波的超声波束合成方法、系统及存储介质 - Google Patents

一种抑制杂波的超声波束合成方法、系统及存储介质 Download PDF

Info

Publication number
WO2024027529A1
WO2024027529A1 PCT/CN2023/109079 CN2023109079W WO2024027529A1 WO 2024027529 A1 WO2024027529 A1 WO 2024027529A1 CN 2023109079 W CN2023109079 W CN 2023109079W WO 2024027529 A1 WO2024027529 A1 WO 2024027529A1
Authority
WO
WIPO (PCT)
Prior art keywords
window function
random
data
signal
ultrasonic
Prior art date
Application number
PCT/CN2023/109079
Other languages
English (en)
French (fr)
Inventor
王丽
郭建军
尹大军
Original Assignee
飞依诺科技股份有限公司
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 飞依诺科技股份有限公司 filed Critical 飞依诺科技股份有限公司
Publication of WO2024027529A1 publication Critical patent/WO2024027529A1/zh

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/56—Details of data transmission or power supply

Definitions

  • the present invention relates to the field of medical ultrasound systems, and in particular, to an ultrasonic beam synthesis method, system and storage medium for suppressing clutter.
  • DAS Delay and Sum
  • the object of the present invention is to provide an ultrasonic beam synthesis method, system and storage medium that suppress clutter.
  • the present invention provides an ultrasonic beam synthesis method that suppresses clutter.
  • the method includes the steps of: obtaining channel data of the ultrasonic echo signal; applying a traditional window function to the channel data of the ultrasonic echo signal.
  • the ultrasonic beam is synthesized according to the processed coordinate information to form image data.
  • the method further includes: mapping the random window function to the traditional window function, obtaining the projection of the random window function on the traditional window function, and subtracting the projected random value from the random window function.
  • the window function obtains the orthogonalized random window function.
  • the method further includes: the traditional window function is at least one of a Hamming window function, a Hanning window function or a Gausswin window function; the random The window function is a uniformly distributed and linearly independent random matrix generated by the linear congruence method LCG or RAND method.
  • the method also includes: the mean value of the orthogonalized random window function is 0; the orthogonalized random window function used at each point in the ultrasonic echo signal channel data The random window functions are all random, but the random window function used between each frame of image data is consistent.
  • the method further includes: the orthogonalization method adopts Schmidt orthogonal method.
  • the "calculating the processed coordinate information of each point in the ultrasonic echo signal channel data based on the correlation coefficient and the signal data set" specifically includes: obtaining the signal in the coordinate system The signal data component of the data in the vertical axis direction; the signal data component is multiplied by the correlation coefficient in this direction to obtain The processed coordinate information of each point in the longitudinal axis direction of the ultrasonic echo signal channel data.
  • Ultrasound beam to form image data specifically includes: summarizing the processed coordinate information of each point in the vertical axis direction according to the horizontal axis direction, synthesizing the ultrasonic beam to form image data.
  • envelope detection and log compression are performed on the image data.
  • the invention also provides an ultrasonic imaging system that suppresses clutter.
  • the system includes: a probe array element module, an ultrasonic echo receiving module, a channel data processing module, an image data conversion module, and an image data processing module.
  • the system is used to execute Any of the above ultrasonic beam synthesis methods to suppress clutter.
  • the present invention also provides a storage medium that stores a computer program.
  • the computer program is executed by a processor, the steps in any one of the above ultrasonic beam synthesis methods for suppressing clutter are implemented.
  • the ultrasonic beam synthesis method for suppressing clutter provided by the present invention not only receives the traditional apodization window function, but also introduces a random window as the receiving window function.
  • the random window has wide-band characteristics and is different from the traditional apodization window function.
  • the window function can effectively suppress side lobes and improve the contrast and spatial resolution of the image.
  • due to the random phase caused by the random function when the traditional window function and the random window act on the signal data and noise data generated by the channel data respectively, after calculating the correlation coefficient between the two data, they will act on the signal data This affects the phase distribution of the signal data, weakens the impact of speckle noise on the image, and improves the image quality.
  • Figure 1 is a flow chart of an ultrasonic beam synthesis method for clutter suppression in an embodiment of the present invention.
  • Figure 2 is a schematic diagram of an ultrasonic imaging system that suppresses clutter in an embodiment of the present invention.
  • Figure 3 is a schematic diagram of a channel data processing module in an embodiment of the present invention.
  • This application discloses an embodiment of an ultrasonic beam synthesis method for clutter suppression.
  • this application provides the method operation steps as described in the following embodiment or flow chart 1, based on routine or without creative labor, the above Among the steps of the method that do not logically have a necessary causal relationship, the execution order of these steps is not limited to the execution order provided in the embodiments of this application.
  • an embodiment of the present invention provides an ultrasonic beam synthesis method for clutter suppression. The method includes the following steps. The method and each step are described below:
  • Step S1 obtain the channel data of the ultrasonic echo signal
  • the probe array element can be configured according to different applications to emit ultrasonic signals, and the receiving end can obtain channel data based on the corresponding echo signals.
  • the ultrasonic imaging system shown in Figure 2 a certain number of probe array element modules 100 are provided to transmit ultrasonic waves at different spatial positions and angles.
  • each channel has M sampling points, and one echo signal can obtain M*N channel sampling data.
  • the maximum N does not exceed the number of channels of the probe array element.
  • a value between 1 and the number of channels of the probe can be selected according to the application. It is preferred to use all channels.
  • Step S2 the traditional window function is applied to the channel data to obtain a signal data set
  • the ultrasonic echo receiving module 200 may use multi-channel signals reflected from the detection object to form a receiving beam.
  • the multi-channel signal reflected from the detection object may be a radio frequency (RF) signal.
  • the radio frequency (RF) signal is reflected by the detection object, and then received by one or more transducers (not shown).
  • multiple channel sampling data of the echo signal are obtained. Different methods are applied to the obtained multiple channel sampling data.
  • the traditional window function generates multiple signal data respectively, and summarizes them to obtain a signal data set.
  • the traditional window function may include a Hamming window function, a Hanning window function or a Gausswin window function.
  • Step S3 the orthogonalized random window function is applied to the channel data to obtain a noise data set
  • the pure noise vector can be obtained by subtracting the projection vector from the noise vector, that is, Orthogonalized random window function, for example, the random window function is a uniformly distributed and linearly independent random matrix that can be generated by the linear congruence method LCG or RAND method.
  • the traditional window function uses the Hamming window function
  • the random window function uses the rand() function in Matlab.
  • the orthogonalized random window function can be calculated using formula (1).
  • w1 is the traditional window function
  • w2 is the orthogonalized random window function
  • N is the number of channels.
  • a signal data generator (not shown) in the channel data processing module 300 that can generate two data sets.
  • the traditional window function w1 is applied to N channel data respectively to obtain N signals.
  • the N signal data are aggregated to obtain the signal data set RF1;
  • the orthogonalized random window function w2 is applied to the N channel data to obtain N noise data, and the N noise data are aggregated to obtain the noise data set.
  • RF2, RF1 and RF2 can be calculated using formulas (2) and (3).
  • x(n) represents the nth channel data
  • w1(n) represents the Hamming function applied to the nth channel data
  • w2(n) represents the orthogonalized random function applied to the nth channel data. window function.
  • the sampling points of the orthogonalized random window applied to each channel data are random, but in each frame of image
  • the random windows used between them are consistent.
  • the sampling points of channel data are P1 and P2.
  • the orthogonalized random windows applied to P1 and P2 are random, assuming they are f1 and f2 respectively, but when generating other frame images , the random windows used by the sampling points P1 and P2 corresponding to other frames are consistent with those of the first frame, that is, f1 and f2 are also used. That is to say, when the first frame of the image is calculated through each sampling point, the orthogonalized random windows used are random, but for the sampling points corresponding to other frames, the orthogonalized random windows used The window is consistent with the first frame.
  • Step S4 Calculate the correlation coefficients of the signal data set and the noise data set in two directions in the coordinate system respectively;
  • the correlation coefficient is at It will be high in the signal area, but low in areas with side lobes or other clutter.
  • the correlation coefficient can be calculated using formula (4).
  • n the vertical axis direction
  • line represents the horizontal axis direction
  • Step S5 Calculate the processed coordinate information of each point in the channel data according to the correlation coefficient and the signal data set;
  • the signal data component of the signal data set in the vertical axis direction is obtained, and multiplied by the correlation coefficient in the vertical axis direction to obtain coordinate information after channel sampling data processing.
  • the processed coordinate information can be calculated using formula (5).
  • RF(n,line) ⁇ (n,line)RF1(n,line).............(5)
  • the correlation coefficient ranges from -1 to 1. If the correlation coefficient is 1, the two signals are identical; if the correlation coefficient is close to or below 0, the two signals are considered uncorrelated; if the correlation coefficient is between 0 and 1, the two signals are considered There is a certain relationship. This coefficient acts on the signal data set, affecting the phase distribution of the RF data, effectively suppressing clutter in the vertical axis direction while retaining the main lobe.
  • Step S6 synthesize the ultrasonic beam according to the coordinate information to form image data
  • step S5 the processed coordinate information of each channel sampling data is obtained, and then summarized according to the horizontal axis direction, and the ultrasonic beam is synthesized to form image data.
  • An embodiment of the present invention provides an ultrasonic imaging system that suppresses clutter, as shown in Figure 2, including: an ultrasonic array element module 100, which is used to send and receive ultrasonic signals; and an ultrasonic echo receiving module 200, which is used to receive and transmit The echo signal acquired by the ultrasonic signal; the channel data processing module 300 is used to process multiple channel data of the ultrasonic echo signal. During the processing process, it is necessary to use traditional window functions and random window functions to calculate the beam synthesis signal data to form an image. data; the image data conversion module 400, which performs envelope detection and Log compression on the image data to obtain processed image data; the image processing module 500, which is used for post-processing the image data and generating images.
  • Figure 3 shows an example of a channel data processing module 300.
  • the channel data processing module 300 uses a traditional window function to apply to the multiple channel data respectively to obtain multiple signal data, and summarizes them to obtain the data signal set RF1;
  • the orthogonalized random window function to apply to multiple channel data to obtain multiple noise data, and summarize them to obtain the noise data set RF2; in the coordinate system, calculate the noise data of RF1 and RF2 in two directions respectively.
  • Correlation coefficient the correlation coefficient is multiplied by the signal data set to obtain the processed coordinate information of each point in the ultrasonic echo signal channel data; according to the processed coordinate information, the ultrasonic beam is synthesized according to the horizontal axis direction to form image data.
  • Embodiments of the present invention also provide a storage medium that stores a computer program.
  • the computer program is executed by a processor, the steps in the ultrasonic beam synthesis method for suppressing clutter as described in any one of the above are implemented. That is to say, the steps in any one of the technical solutions in the ultrasonic beam synthesis method for suppressing clutter are implemented.
  • the ultrasonic beam synthesis method, system and storage medium for clutter suppression provided by the present invention first obtain the ultrasonic echo Channel data of the signal; select two window functions, a traditional window function and a random window function, use Schmidt orthogonalization to calculate the orthogonal basis functions of the two window functions, and obtain the orthogonalized random window function;
  • the traditional window function and the orthogonalized random window function are applied to the channel data respectively to obtain the signal data set RF1 and the noise data set RF2, and calculate the correlation coefficient between the two data sets; according to the correlation coefficient and the signal
  • the data set is synthesized by ultrasound beams to form image data.
  • the present invention utilizes the random window to have wide-band characteristics.
  • the correlation coefficient When used in conjunction with the traditional window function, it can retain the main lobe, better suppress the side lobes, and improve the contrast and spatial resolution of the image. At the same time, due to the random phase brought by the random window function The reason is that after the correlation coefficient is calculated, it will be applied to the signal data to obtain the beam synthesis data, which affects the phase distribution of the beam synthesis data, weakens the impact of speckle noise on the image, and improves the image quality.

Landscapes

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

Abstract

一种抑制杂波的超声波束合成方法、系统及存储介质，该方法包括通过将传统窗函数和随机窗函数分别作用到通道数据产生的信号数据和噪声数据时，计算这两种数据之间的相关性系数之后并作用到信号数据上，以改变信号数据的相位分布，减弱了散斑噪声对图像的影响，提高了图像质量。

Description

一种抑制杂波的超声波束合成方法、系统及存储介质 技术领域
本发明涉及医学超声系统领域，尤其涉及一种抑制杂波的超声波束合成方法、系统及存储介质。
背景技术
传统的波束合成是将接收到的通道数据，根据发射焦点和接收焦点与超声探头阵元之间的几何关系计算出声波到达的时间关系，将通道数据在时间上对齐，同时使用一个接收变迹窗函数，再叠加起来，即Delay and Sum(DAS)。
在实际的超声图像扫查中，人体组织由于运动干扰和复杂的声阻抗特性，加上声波的传播特性复杂，会产生散射、多重反射，会出现各种干扰噪声，使用传统的DAS方案不能去除各种杂波信号，使图像不能达到较好的对比度和空间分辨率。
发明内容
本发明的目的在于提供一种抑制杂波的超声波束合成方法、系统及存储介质。
为了实现上述发明目的，本发明提供一种抑制杂波的超声波束合成方法，所述方法包括步骤：获取超声回波信号的通道数据；将传统窗函数作用到所述超声回波信号的通道数据中，获得信号数据集，以及将随机窗函数作用到所述超声回波信号的通道数据中，获得噪声数据集，其中，所述随机窗函数为所述传统窗函数正交化后得到的；在坐标系中，分别计算所述信号数据集和所述噪声数据集在两个方向的相关性系数；根据所述相关性系数以及信号数据集，计算超声回波信号通道数据中每个点处理后的坐标信息；根据所述处理后的坐标信息合成所述超声波束，形成图像数据。
作为本发明的进一步改进，所述方法还包括：将所述随机窗函数映射到所述传统窗函数上，获得随机窗函数在传统窗函数上的投影，用随机窗函数减去投影后的随机窗函数得到正交化后的随机窗函数。
作为本发明的进一步改进，所述方法还包括：所述传统窗函数为汉明(Hamming)窗函数、汉宁(Hanning)窗函数或者高斯(Gausswin)窗函数中的至少一种；所述随机窗函数是通过线性同余法LCG或RAND方法产生的均匀分布且线性无关的随机矩阵。作为本发明的进一步改进，所述方法还包括：所述正交化后的随机窗函数的均值为0；在所述超声回波信号通道数据中每个点使用的所述正交化后的随机窗函数均是随机的，但每帧图像数据之间使用的该随机窗函数是一致的。
作为本发明的进一步改进，所述方法还包括：所述正交化方法采用的是施密特正交法。
作为本发明的进一步改进，所述“根据所述相关性系数以及信号数据集，计算超声回波信号通道数据中每个点处理后的坐标信息”具体包括：在坐标系中，获取所述信号数据在纵轴方向的信号数据分量；所述信号数据分量乘以该方向的相关性系数得 到所述超声回波信号通道数据在纵轴方向的每个点处理后的坐标信息。
作为本发明的进一步改进，所述“根据所述处理后的坐标信息合成所述
超声波束，形成图像数据”具体包括：将所述纵轴方向的每个点处理后的坐标信息按照横轴方向进行汇总，合成超声波束，形成图像数据。
作为本发明的进一步改进，对所述图像数据执行包络检测和日志压缩。
本发明还提供一种抑制杂波的超声成像系统，该系统包括：探头阵元模块、超声回波接收模块、通道数据处理模块、图像数据转换模块、图像数据处理模块，所述系统用于执行上述任意一项抑制杂波的超声波束合成方法。
本发明还提供一种存储介质，所述存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上任意一项抑制杂波的超声波束合成方法中的步骤。
与现有技术相比，本发明提供的抑制杂波的超声波束合成方法，在接收传统的变迹窗函数的同时，引入了随机窗作为接收窗函数，随机窗具有宽频带特性，与传统的窗函数作用可以很好的抑制旁瓣，提高图像的对比度和空间分辨率。同时由于随机函数带来的随机相位原因，在传统窗函数和随机窗分别作用到通道数据产生的信号数据和噪声数据时，计算这两种数据之间的相关性系数之后，会作用到信号数据上，影响了信号数据的相位分布，减弱了散斑噪声对图像的影响，提高了图像质量。
附图说明
图1是本发明实施例中抑制杂波的超声波束合成方法的流程图。
图2是本发明实施例中抑制杂波的超声成像系统的示意图。
图3是本发明实施例中通道数据处理模块的示意图。
具体实施方式
以下将结合附图所示的具体实施方式对本发明进行详细描述。但这些实施方式并不限制本发明，本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本发明的保护范围内。
需要说明的是，术语“包括”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。此外，术语“第一”、“第二”等仅用于描述目的，而不能理解为指示或暗示相对重要性。
本申请公开了一个实施方式的一种抑制杂波的超声波束合成方法，虽然本申请提供了如下述实施方式或流程图1所述的方法操作步骤，但是基于常规或者无需创造性的劳动，所述方法在逻辑性上不存在必要因果关系的步骤中，这些步骤的执行顺序不限于本申请实施方式中所提供的执行顺序。如图1所示，本发明实施例提供了一种抑制杂波的超声波束合成方法，所述方法包括以下若干步骤，下面对所述方法及每个步骤分别进行说明：
步骤S1，获取超声回波信号的通道数据；
例如，可以根据不同应用配置，探头阵元发射超声波信号，接收端根据对应的回波信号获取通道数据。在图2所示的超声成像系统中设置一定数目探头阵元模块100进行不同的空间位置和角度的超声波发射。对于N个通道，每个通道M个采样点，一个回波信号可获得M*N个通道采样数据。其中N最高不超过探头阵元的通道个数，可根据应用选择1到探头的通道个数之间的值，优选是全通道使用。
步骤S2，传统窗函数作用到所述通道数据上获得信号数据集；
具体地，超声回波接收模块200可使用从探测对象反射的多通道信号来形成接收波束。例如，从探测对象反射的多通道信号可以是射频(RF)信号。射频(RF)信号被探测对象反射，再被一个或多个换能器(未标出)接收，最后得到回波信号的多个通道采样数据，对所得到的多个通道采样数据应用了不同的传统窗函数，分别产生多个信号数据，对其进行汇总得到信号数据集。例如，传统窗函数可包括汉明(Hamming)窗函数、汉宁(Hanning)窗函数或者高斯(Gausswin)窗函数。
步骤S3，正交化后的随机窗函数作用到所述通道数据获得噪声数据集；
具体地，在获取回波信号的多个通道采样数据后，选用传统窗函数和随机窗函数作用在通道数据上，可以获得两路信号空间向量，一路称之为信号空间向量，另一路称之为噪声空间向量。这两个窗函数并不一定是正交的，为了得到更加纯净的噪声空间，可将两个窗函数进一步优化，即正交化。可以采用施密特正交化的方法，将噪声空间向量映射到信号空间向量中，获得噪声空间向量在信号空间向量上的投影，用噪声向量减去投影向量即可得到纯净的噪声向量，即正交化后的随机窗函数，例如，随机窗函数是可以通过线性同余法LCG或RAND方法产生的均匀分布且线性无关的随机矩阵。
于一具体实施例中，所述传统窗函数选用汉明窗函数，随机窗函数选用Matlab中的rand()函数，例如，汉明窗函数和随机窗函数定义分别如下：w1＝hamming(N),w2＝rand(N)，采用施密特正交化方法，得到正交化后的随机窗函数，该正交化后的随机窗函数可使用公式(1)来计算。
其中w1为传统窗函数，w2为正交化后的随机窗函数，N为通道数。
于一具体实施例中，通道数据处理模块300中有个信号数据产生器(未标出)可产生两个数据集，例如，传统窗函数w1分别作用到N个通道数据上，得到N个信号数据，汇总这N个信号数据得到信号数据集RF1；同样的，正交化后的随机窗函数w2分别作用到N个通道数据上得到N个噪声数据，汇总这N个噪声数据得到噪声数据集RF2，RF1和RF2可使用公式(2)和(3)来计算。

其中x(n)表示第n个通道数据，w1(n)表示作用到第n个通道数据上的汉明函数，w2(n)表示作用到第n个通道数据上的正交化后的随机窗函数。
在坐标系中，为了确保正交化后的随机窗在横向和纵向都是线性无关的，正交化后的随机窗作用到各个通道数据的采样点上都是随机的，但是在每帧图像之间采用的随机窗又是保持一致的。例如，通道数据的采样点有P1和P2，在产生第一帧图像时，作用到P1和P2上的正交化随机窗是随机的，假设分别是f1和f2，但是在产生其它帧图像时，其它帧对应的采样点P1和P2使用的随机窗和第一帧保持一致，即也是使用f1和f2。也就是说，通过各个采样点计算出的第一帧图像时，采用的正交化后的随机窗都是随机的，但是对其它帧对应的采样点而言，使用的正交化后的随机窗和第一帧是保持一致的。
步骤S4，分别计算信号数据集和噪声数据集在坐标系中两个方向的相关性系数；
具体地，在坐标系中，根据获取到的信号数据集RF1和噪声数据集RF2，计算RF1和RF2在每个横轴点处的纵轴方向各个数据点的相关性系数，该相关性系数在信号区域内会高，但在有旁瓣或者其他杂波的区域会低。
作为示例，所述相关性系数可使用公式(4)来计算。
其中，在坐标系中，n表示纵轴方向，line表示横轴方向，K表示纵轴方向上参与相关性系数计算的点的个数。例如，当K＝1时，表示计算纵轴方向的n-1、n和n+1三个点的相关性系数。
步骤S5，根据所述相关性系数和信号数据集，计算通道数据中每个点处理后的坐标信息；
具体地，在坐标系中，获取所述信号数据集在纵轴方向的信号数据分量，乘以纵轴方向的相关性系数，得到通道采样数据处理后的坐标信息。
作为示例，所述处理后的坐标信息可使用公式(5)来计算。
RF(n,line)＝ρ(n,line)RF1(n,line)………..(5)
该相关性系数取值范围在-1到1之间。如果相关性系数为1，则两个信号是相同的；如果相关性系数接近或低于0，则认为这两个信号不相关；如果相关性系数在0和1之间则认为这两个信号有一定的关系。该系数作用到信号数据集上，影响了RF数据的相位分布，在保留主瓣的同时有效的抑制纵轴方向的杂波。
步骤S6，根据坐标信息合成超声波束，形成图像数据；
具体地，基于步骤S5获得各个通道采样数据处理后的坐标信息，再按照横轴方向进行汇总，合成超声波束，形成图像数据。
本发明一实施例提供了一种抑制杂波的超声成像系统，如图2所示，包括：超声阵元模块100，用户发送并接收超声信号；超声回波接收模块200，用于接收通过发送超声信号所获取的回波信号；通道数据处理模块300，用于处理超声回波信号的多个通道数据，在处理过程中需要使用传统窗函数和随机窗函数计算出波束合成信号数据，形成图像数据；图像数据转换模块400，对所述图像数据进行包络检测和Log压缩获得处理后的图像数据；图像处理模块500，用于图像数据后处理以及产生图像。
图3示出通道数据处理模块300的示例。参照图3，通道数据处理模块300在获取多个超声回波信号的通道数据后，使用传统的窗函数分别作用到多个通道数据上得到多个信号数据，对其汇总得到数据信号集RF1；同样的，使用正交化后的随机窗函数分别作用到多个通道数据上得到多个噪声数据，对其汇总得到噪声数据集RF2；在坐标系中，分别计算RF1和RF2在两个方向的相关性系数；所述相关性系数乘以信号数据集得到超声回波信号通道数据中每个点处理后的坐标信息；根据所述处理后的坐标信息按照横轴方向进行汇总合成超声波束，形成图像数据。
本发明实施例还提供了一种存储介质，所述存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上任意一项所述的抑制杂波的超声波束合成方法中的步骤，也就是说，实现所述抑制杂波的超声波束合成方法中任意一个技术方案中的步骤。
综上，本发明提供的杂波抑制的超声波束合成方法、系统及存储介质，先获取超声回波 信号的通道数据；选取两个窗函数，传统窗函数和随机窗函数，采用施密特正交化计算出所述两个窗函数的正交基函数，得到正交化后的随机窗函数；传统窗函数和正交化后的随机窗函数分别作到通道数据上，获得信号数据集RF1和噪声数据集RF2，计算所述两个数据集之间的相关性系数；根据相关性系数和信号数据集得到超声波束合成，形成图像数据。本发明利用随机窗具有宽频带特性，与传统的窗函数配合使用可以保留主瓣的同时，更好的抑制旁瓣，提高图像的对比度和空间分辨率，同时由于随机窗函数带来的随机相位的原因，在计算出相关性系数后会作用到信号数据上得到波束合成数据，影响了波束合成数据的相位分布，减弱了散斑噪声对图像的影响，提高了图像质量。
应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施方式中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。

Claims (10)

  1. 一种抑制杂波的超声波束合成方法，其特征在于，所述方法包括：
    获取超声回波信号的通道数据；
    将传统窗函数作用到所述超声回波信号的通道数据中，获得信号数据集，以及将随机窗函数作用到所述超声回波信号的通道数据中，获得噪声数据集；其中，所述随机窗函数为所述传统窗函数正交化后得到的；
    在坐标系中，分别计算所述信号数据集和所述噪声数据集在两个方向的相关性系数；
    根据所述相关性系数以及信号数据集，计算超声回波信号通道数据中每个点处理后的坐标信息；
    根据所述处理后的坐标信息合成所述超声波束，形成图像数据。
  2. 根据权利要求1所述的方法，其特征在于，所述方法还包括：
    将所述随机窗函数映射到所述传统窗函数上，获得随机窗函数在传统窗函数上的投影，用随机窗函数减去投影后的随机窗函数得到正交化后的随机窗函数。
  3. 根据权利要求1所述的方法，其特征在于，所述方法还包括：
    所述传统窗函数为汉明窗函数、汉宁窗函数或者高斯窗函数中的至少一种；
    所述随机窗函数是通过线性同余法LCG或RAND方法产生的均匀分布且线性无关的随机矩阵。
  4. 根据权利要求1所述的方法，其特征在于，所述方法还包括：
    所述正交化后的随机窗函数的均值为0；
    在所述超声回波信号通道数据中每个点使用的所述正交化后的随机窗函数均是随机的，但每帧图像数据之间使用的该随机窗函数是一致的。
  5. 根据权利要求1所述的方法，其特征在于，所述方法还包括：
    所述正交化方法采用的是施密特正交法。
  6. 根据权利要求1所述的方法，其特征在于，所述“根据所述相关性系数以及信号数据集，计算超声回波信号通道数据中每个点处理后的坐标信息”具体包括：
    在坐标系中，获取所述信号数据在纵轴方向的信号数据分量；
    所述信号数据分量乘以该方向的相关性系数得到所述超声回波信号通道数据在纵轴方向的每个点处理后的坐标信息。
  7. 根据权利要求1所述的方法，其特征在于，所述“根据所述处理后的坐标信息合成所述超声波束，形成图像数据”具体包括：
    将所述纵轴方向的每个点处理后的坐标信息按照横轴方向进行汇总，合成超声波束，形成图像数据。
  8. 根据权利要求7所述的方法，其特征在于，所述方法还包括：对所述图像数据 执行包络检测和日志压缩。
  9. 一种抑制杂波的超声成像系统，其特征在于，该系统包括：探头阵元模块、超声回波接收模块、通道数据处理模块、图像数据转换模块、图像数据处理模块，所述系统用于执行如权利要求1-8任一项所述的方法。
  10. 一种存储介质，所述存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时候实现如权利要求1-8任意一项所述的抑制杂波的超声波束合成方法中的步骤。
PCT/CN2023/109079 2022-08-05 2023-07-25 一种抑制杂波的超声波束合成方法、系统及存储介质 WO2024027529A1 (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202210936540.5A CN115137392A (zh) 2022-08-05 2022-08-05 一种抑制杂波的超声波束合成方法、系统及存储介质
CN202210936540.5 2022-08-05

Publications (1)

Publication Number Publication Date
WO2024027529A1 true WO2024027529A1 (zh) 2024-02-08

Family

ID=83413846

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2023/109079 WO2024027529A1 (zh) 2022-08-05 2023-07-25 一种抑制杂波的超声波束合成方法、系统及存储介质

Country Status (2)

Country Link
CN (1) CN115137392A (zh)
WO (1) WO2024027529A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115137392A (zh) * 2022-08-05 2022-10-04 飞依诺科技股份有限公司 一种抑制杂波的超声波束合成方法、系统及存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63154165A (ja) * 1986-12-19 1988-06-27 富士通株式会社 超音波ドプラ血流計
JPH06319735A (ja) * 1993-05-17 1994-11-22 Yokogawa Medical Syst Ltd 超音波診断装置
JP2003079624A (ja) * 2001-09-14 2003-03-18 Ge Medical Systems Global Technology Co Llc ドプラ音発生方法および超音波診断装置
JP2007298503A (ja) * 2006-04-03 2007-11-15 Mitsubishi Electric Corp 伝搬遅延時間測定装置およびレーダ装置
CN106236147A (zh) * 2016-10-18 2016-12-21 湖南中聚光电科技发展有限公司 一种超声回波数字解调方法及装置
CN113687356A (zh) * 2021-09-16 2021-11-23 中国科学院空天信息创新研究院 一种机载多通道圆迹sar运动目标检测与估计方法
KR20220037835A (ko) * 2020-09-18 2022-03-25 엘아이지넥스원 주식회사 표적 식별 향상을 위한 부엽 감소 방법 및 장치
CN115137392A (zh) * 2022-08-05 2022-10-04 飞依诺科技股份有限公司 一种抑制杂波的超声波束合成方法、系统及存储介质

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63154165A (ja) * 1986-12-19 1988-06-27 富士通株式会社 超音波ドプラ血流計
JPH06319735A (ja) * 1993-05-17 1994-11-22 Yokogawa Medical Syst Ltd 超音波診断装置
JP2003079624A (ja) * 2001-09-14 2003-03-18 Ge Medical Systems Global Technology Co Llc ドプラ音発生方法および超音波診断装置
JP2007298503A (ja) * 2006-04-03 2007-11-15 Mitsubishi Electric Corp 伝搬遅延時間測定装置およびレーダ装置
CN106236147A (zh) * 2016-10-18 2016-12-21 湖南中聚光电科技发展有限公司 一种超声回波数字解调方法及装置
KR20220037835A (ko) * 2020-09-18 2022-03-25 엘아이지넥스원 주식회사 표적 식별 향상을 위한 부엽 감소 방법 및 장치
CN113687356A (zh) * 2021-09-16 2021-11-23 中国科学院空天信息创新研究院 一种机载多通道圆迹sar运动目标检测与估计方法
CN115137392A (zh) * 2022-08-05 2022-10-04 飞依诺科技股份有限公司 一种抑制杂波的超声波束合成方法、系统及存储介质

Also Published As

Publication number Publication date
CN115137392A (zh) 2022-10-04

Similar Documents

Publication Publication Date Title
US8506487B2 (en) Ultrasound imaging device
US20130258805A1 (en) Methods and systems for producing compounded ultrasound images
JP5483905B2 (ja) 超音波装置
US11372094B2 (en) Reverberation artifact cancellation in ultrasonic diagnostic images
US20060241429A1 (en) Aberration correction with broad transmit beams in medical ultrasound
JP2006346459A (ja) 超音波空間複合イメージングのリアルタイム運動補正のための方法及び装置
CN102727255B (zh) 一种超声图像空间复合成像的方法及装置
WO2024027529A1 (zh) 一种抑制杂波的超声波束合成方法、系统及存储介质
CN110101411B (zh) 超声成像空间复合方法及系统
US11650300B2 (en) Ultrasound system and method for suppressing noise using per-channel weighting
KR101792589B1 (ko) 빔포머, 진단시스템, 의료영상시스템 및 진단영상을 표시하는 방법
KR101552427B1 (ko) 초음파 영상의 스페클 저감장치
US20070083109A1 (en) Adaptive line synthesis for ultrasound
US20150327840A1 (en) Ultrasonic diagnostic device and correction method
WO2015166869A1 (ja) 超音波撮像装置
CN110536647B (zh) 用于在存在运动的情况下处理和复合超声图像的系统、方法和计算机可读介质
US11484292B2 (en) Ultrasound signal processing device that uses synthetic aperture method and delay and sum method
US9304191B2 (en) Subject information obtaining apparatus, subject information obtaining method, and program
CN111278363B (zh) 超声成像设备、系统及其超声造影成像的图像增强方法
CN112890855B (zh) 多波束p次根压缩相干滤波波束合成方法及装置
US9289190B2 (en) Ultrasound strain imaging via pixel frame and window correlation
JP6863817B2 (ja) 超音波撮像装置
KR102115446B1 (ko) 의료용 초음파 영상에서 질량중심 위치의 가중치 연산을 이용하여 부엽을 억제하는 방법
JP2824673B2 (ja) 超音波診断装置
US20240142620A1 (en) Ultrasound imaging method and ultrasound imaging system

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23849237

Country of ref document: EP

Kind code of ref document: A1