CN108113703B - Method and device for generating apodization value for ultrasonic scanning - Google Patents

Method and device for generating apodization value for ultrasonic scanning Download PDF

Info

Publication number
CN108113703B
CN108113703B CN201711397019.4A CN201711397019A CN108113703B CN 108113703 B CN108113703 B CN 108113703B CN 201711397019 A CN201711397019 A CN 201711397019A CN 108113703 B CN108113703 B CN 108113703B
Authority
CN
China
Prior art keywords
array element
apodization
value
depth
section
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711397019.4A
Other languages
Chinese (zh)
Other versions
CN108113703A (en
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.)
Feiyinuo Technology Co ltd
Original Assignee
Vinno Technology Suzhou Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Vinno Technology Suzhou Co Ltd filed Critical Vinno Technology Suzhou Co Ltd
Priority to CN201711397019.4A priority Critical patent/CN108113703B/en
Publication of CN108113703A publication Critical patent/CN108113703A/en
Application granted granted Critical
Publication of CN108113703B publication Critical patent/CN108113703B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4411Device being modular

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)
  • Veterinary Medicine (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)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

The invention provides a method for generating an apodization value for ultrasonic scanning, wherein the depth of the ultrasonic scanning is divided into N adjacent segments in sequence, and N is an integer greater than 1; the method comprises the following steps: acquiring a first apodization value of each array element in the kth section of depth value, wherein k is an integer and is more than or equal to 1 and less than k < N; carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value; and obtaining a second apodization value of each array element in the k-th section of the depth value based on the first synthesized beam. Thereby enabling the generation of a suitable apodization value.

Description

Method and device for generating apodization value for ultrasonic scanning
Technical Field
The present invention relates to the field of ultrasound scanning technologies, and in particular, to a method and an apparatus for generating an apodization value for ultrasound scanning.
Background
The beam forming technology has been a hot technique studied in medical ultrasonic imaging systems for a long time, and the quality of the beam forming technology can directly affect the quality of final imaging. The beam synthesis technology comprises a focusing technology, an apodization technology and a dynamic aperture technology. Apodization is one of the more critical techniques.
In an ultrasound scanning apparatus, an ultrasound probe is generally provided with a plurality of array elements, and the ultrasound scanning apparatus obtains a composite signal according to ultrasound echoes received from the plurality of array elements, and it is understood that each array element contributes differently to the composite signal. It can be understood that if the signals of the respective array elements are simply added in equal amplitude, the side lobe of the ultrasonic echo is large, which results in increased artifact and reduced image contrast in the final ultrasonic imaging; therefore, in order to suppress the side lobe and improve the contrast of the image, the echo signals of different array elements are multiplied by different coefficients and then added, and the process is called apodization. It is generally considered that the array elements in the center of the aperture contribute greatly to the beamformed signals, and the array elements on both sides of the aperture contribute little to the beamformed signals. The echo signal of the array element with large contribution is multiplied by a larger coefficient, and the echo signal of the array element with small contribution is multiplied by a smaller coefficient. Apodization can reduce the sidelobe level of the beamformed beam and improve the contrast of the image. However, the use of apodization increases the main lobe width and reduces the lateral resolution of the image.
Therefore, how to obtain a suitable apodization value becomes a problem to be solved urgently.
Disclosure of Invention
The invention aims to provide a method and a device for generating apodization values for ultrasonic scanning.
In order to achieve one of the above objects, an embodiment of the present invention provides a method for generating an apodization value for an ultrasound scan, a depth of the ultrasound scan being sequentially divided into N adjacent segments, N being an integer greater than 1; the method comprises the following steps:
acquiring a first apodization value of each array element in the kth section of depth value, wherein k is an integer and is more than or equal to 1 and less than N;
carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value;
and obtaining a second apodization value of each array element in the k +1 th section of depth value based on the first synthesized beam.
As a further improvement of an embodiment of the present invention, the obtaining, based on the first synthesized beam, a second apodization value of each array element at the k +1 th section of the depth value includes:
obtaining a demodulation filter type F (k), a demodulation bandwidth [ F [ ]L(k),FH(k)],FL(k) High-pass cut-off frequency, F, for the k-th segment depth valueH(k) A low-pass cut-off frequency which is the k-th section depth value;
obtaining a first spectral energy
Figure BDA0001518728460000021
Energy of the second spectrum
Figure BDA0001518728460000022
Figure BDA0001518728460000023
Figure BDA0001518728460000024
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<fL(k)<fH(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, PK1 is the demodulation filter bandwidth range [ fL(k),fH(k)]PK2 is a demodulation filterIn the bandwidth range 0, fs/2]The spectral energy of (a).
As a further improvement of an embodiment of the present invention, the method further comprises the following steps: and demodulating and filtering the first synthesized beam to obtain a baseband signal.
As a further improvement of an embodiment of the present invention, the delaying, superimposing, and apodizing each array element includes: carrying out time delay processing on each array element; for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and superposing the first echo signals multiplied by all the array elements.
As a further improvement of an embodiment of the present invention, the delaying, superimposing, and apodizing each array element includes: for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and performing delay superposition on the first echo signals multiplied by all the array elements.
An embodiment of the present invention provides an apparatus for generating an apodization value for ultrasonic scanning, wherein a depth of the ultrasonic scanning is sequentially divided into N adjacent segments, N being an integer greater than 1; the system comprises the following modules: the first initialization module is used for acquiring a first apodization value of each array element in the kth section of the depth value, wherein k is an integer and is more than or equal to 1 and less than N;
the first synthesis beam generation module is used for multiplying the first echo signal of the array element by the first dynamic apodization value of the array element for each array element; carrying out time delay superposition on the first echo signals multiplied by all the array elements to obtain a first synthesized beam of a k-th section of depth value;
and the first processing module is used for obtaining a second apodization value of each array element at the k +1 th section of depth value based on the first synthesized beam.
One embodiment of the present invention provides a method for generating an apodization value for ultrasound scanning, wherein the depth of the ultrasound scanning is sequentially divided into N adjacent segments, N being an integer greater than 1; the method comprises the following steps:
acquiring a first apodization value of each array element in the kth section of depth value, wherein k is an integer and is more than or equal to 1 and less than N;
carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value;
and obtaining a second apodization value of each array element in the k-th section of the depth value based on the first synthesized beam.
As a further improvement of an embodiment of the present invention, the obtaining, based on the first synthesized beam, a second apodization value of each array element at the k-th section depth value includes:
obtaining a demodulation filter type F (k), a demodulation bandwidth [ F [ ]L(k),FH(k)],FL(k) High-pass cut-off frequency, F, for the k-th segment depth valueH(k) A low-pass cut-off frequency which is the k-th section depth value;
obtaining a first spectral energy
Figure BDA0001518728460000031
Energy of the second spectrum
Figure BDA0001518728460000032
Figure BDA0001518728460000033
Figure BDA0001518728460000034
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<fL(k)<fH(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, PK1 is the demodulation filter bandwidth range [ fL(k),fH(k)]PK2 is a demodulation filter with a bandwidth range of [0, fs/2 ]]The spectral energy of (a).
As a further improvement of an embodiment of the present invention, the method further comprises the following steps: and demodulating and filtering the first synthesized beam to obtain a baseband signal.
As a further improvement of an embodiment of the present invention, the delaying, superimposing, and apodizing each array element includes: carrying out time delay processing on each array element; for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and superposing the first echo signals multiplied by all the array elements.
As a further improvement of an embodiment of the present invention, the delaying, superimposing, and apodizing each array element includes: for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and performing delay superposition on the first echo signals multiplied by all the array elements.
An embodiment of the present invention provides an apparatus for generating an apodization value for ultrasonic scanning, wherein a depth of the ultrasonic scanning is sequentially divided into N adjacent segments, N being an integer greater than 1; the system comprises the following modules:
the second initialization module is used for acquiring a first apodization value of each array element in the kth section of the depth value, wherein k is an integer and is more than or equal to 1 and less than N;
the second synthetic beam generating module is used for carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value;
and the second processing module is used for obtaining a second apodization value of each array element in the k-th section of the depth value based on the first synthesized beam.
Compared with the prior art, the invention has the technical effects that: the invention provides a method for generating an apodization value for ultrasonic scanning, wherein the depth of the ultrasonic scanning is divided into N adjacent segments in sequence, and N is an integer greater than 1; the method comprises the following steps: acquiring a first apodization value of each array element in the kth section of depth value, wherein k is an integer and is more than or equal to 1 and less than N; carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value; and obtaining a second apodization value of each array element in the k +1 th section of depth value based on the first synthesized beam. Thereby enabling the generation of a suitable apodization value.
Drawings
FIG. 1 is a flow chart illustrating a method for generating apodization values according to a first embodiment of the invention;
FIG. 2 is a signal control flow diagram according to a first embodiment of the present invention;
fig. 3a is a first demodulation bandwidth setting in the first embodiment of the present invention;
FIG. 3b is a second demodulation bandwidth setting according to the first embodiment of the present invention;
fig. 4 is a flowchart illustrating a method for generating an apodization value according to a second embodiment of the present invention.
Detailed Description
The present invention will be described in detail below with reference to embodiments shown in the drawings. These embodiments are not intended to limit the present invention, and structural, methodological, or functional changes made by those skilled in the art according to these embodiments are included in the scope of the present invention.
The embodiment one of the present invention provides a method for generating an apodization value for ultrasonic scanning, where a depth of ultrasonic scanning is sequentially divided into N adjacent segments, where N is an integer greater than 1, and optionally, the depth of the scanning may be sequentially divided into N equal segments; as shown in fig. 1, the method comprises the following steps:
step 101: acquiring a first apodization value of each array element in the kth section of depth value, wherein k is an integer and is more than or equal to 1 and less than N;
step 102: carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value;
step 103: and obtaining a second apodization value of each array element in the k +1 th section of depth value based on the first synthesized beam.
Preferably, the obtaining a second apodization value of each array element at the k +1 th section of the depth value based on the first synthesized beam includes:
obtaining a demodulation filter type F (k), a demodulation bandwidth [ F [ ]L(k),FH(k)],FL(k) High-pass cut-off frequency, F, for the k-th segment depth valueH(k) A low-pass cut-off frequency which is the k-th section depth value;
obtaining a first spectral energy
Figure BDA0001518728460000051
Energy of the second spectrum
Figure BDA0001518728460000052
Figure BDA0001518728460000053
Figure BDA0001518728460000054
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<fL(k)<fH(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, PK1 is the demodulation filter bandwidth range [ fL(k),fH(k)]PK2 is a demodulation filter with a bandwidth range of [0, fs/2 ]]The spectral energy of (a).
Preferably, the method further comprises the following steps: and demodulating and filtering the first synthesized beam to obtain a baseband signal.
Preferably, each array element is subjected to delay, superposition and apodization processing, including: carrying out time delay processing on each array element; for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and superposing the first echo signals multiplied by all the array elements.
Preferably, each array element is subjected to delay, superposition and apodization processing, including: for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and performing delay superposition on the first echo signals multiplied by all the array elements.
Here, the method for generating the apodization value may be performed using a control system in the ultrasound scanning apparatus, and therefore, a software module may be provided at the control system for performing the method for generating the apodization value, and fig. 2 shows a signal control flow diagram between the software modules, which has a workflow: (1) receiving echo signals of the kth section of depth value of each array element acquired by a probe, entering a beam synthesis module, calculating the delay of each array element by a focusing delay calculation module, and performing delay processing on the signals of each array element by a signal delay module according to the calculated delay; (2) apodization coefficient calculation modeThe block generates an apodization weighting coefficient according to the aperture calculated by the aperture calculating module and a preset apodization window, the coefficient multiplying module multiplies and adds signals of each array element after time delay according to the calculated apodization coefficient, and outputs the signals to the demodulating and filtering module to generate baseband signals; the baseband signal can be processed differently according to different imaging modes and finally converted into an image; (3) demodulation bandwidth of demodulation filter fL(k),fH(k)]For demodulation bandwidth f dynamically varying with depth according to different applicationsL(k),fH(k)]The choice of (d) also depends on the frequency band of interest in the echo signal for different applications, k characterizing the depth range, fL(k) High-pass cut-off frequency, f, of demodulation filter for k-th depth valueH(k) The low-pass cut-off frequency of the demodulation filter corresponding to the k-th section depth value. Different demodulation filters and different bandwidth settings will result in completely different image qualities. In the prior art, the apodization coefficient calculating module of the beam synthesis is before the demodulation filter, the apodization coefficient is generated according to the full bandwidth signal without considering the perceptive demodulation bandwidth set in the demodulation filter, so that the apodization coefficient cannot be completely matched with the perceptive demodulation bandwidth.
Fig. 3 shows a frequency response diagram of a typical wall filter. Fig. 3(a) and 3(b) show different settings of the demodulation bandwidth for two different filter types: 1. the filter arrangement of fig. 3(a) is of interest for the high frequency part of the signal and the filter has a steeper channel, indicating that the signal-to-noise ratio of the signal is not too strong and a filter with a steep passband is required to obtain contrast; 2. the filter arrangement of fig. 3(b) is of interest for low frequency parts in the signal and the filter passband is flat, indicating that the signal to noise ratio is strong and good resolution is obtained with a flat filter. In addition, when the bandwidth of the demodulation filter is set narrow, it is for obtaining a good contrast. When the bandwidth of the demodulation filter is set wide, a good resolution is obtained. In the apodization control of the medical ultrasonic imaging system in the prior art, the setting of the demodulation filter is not considered, so that the apodization coefficient generated is not necessarily optimally matched with the perceptual bandwidth of the signal, and the optimal image quality cannot be achieved.
The demodulation filter parameters are used for controlling the generation of the apodization coefficients of the apodization coefficient calculation module in the beam forming so as to realize the optimal matching of apodization and demodulation filtering and improve the image quality. Filter type f (k) and demodulation bandwidth f of demodulation filterL(k),fH(k)]The selection of the perceived bandwidth of the signal at different depths is determined. f. ofL(k),fH(k) The range of depth k is consistent with the range of depths for which apodization coefficients are generated, as opposed to different application settings for different probes.
Where k represents depth and N is the total number of segments of depth. The apodization coefficients are dynamically generated, and the apodization window function can be a single window function or can be different from the apodization window function selected at different depths. The window function at depth k is WLUT () and, considering the demodulation filter settings, the n-based apodization coefficient is
Figure BDA0001518728460000071
Figure BDA0001518728460000072
n is 0, 1, 2, …, M-1, n is the index of the open array element, M is the number of open array elements, C is a constant, C is generally in the range of 0-2, but is not limited to this range and can be determined according to the actual image optimization, fs is the sampling rate of the demodulation filter input signal, 0<fL(k)<fH(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, Pk1 is the demodulation filter bandwidth range [ fL(k),fH(k)]Pk2 is a demodulation filter with a bandwidth range [0, fs/2 ]]The spectral energy of (a). From the above equation, it can be seen that when the frequency response of the demodulation filter is at fL(k),fH(k)]The flatter the more apodization coefficients are generated for better resolution. Conversely, when the frequency response of the demodulation filter is at [ f ]L(k),fH(k)]The steeper the apodization coefficient is produced, the better the contrast is obtained. When [ fL(k),fH(k)]The larger the range of (a), the larger the apodization coefficient generated to obtain better resolution. On the contrary, when [ fL(k),fH(k)]The smaller the range of (a) to obtain better contrast. When [ fL(k),fH(k)]The higher the frequency component is, the smaller the apodization coefficient is generated to obtain better contrast. On the contrary, when [ fL(k),fH(k)]The more the low frequency components are biased, the larger the apodization coefficients are generated to obtain better resolution. The apodization coefficient generated by the method can be dynamically adjusted according to the parameter of the demodulation filter, and the optimal balance is made in the resolution and the contrast of the image. The parameters of the demodulation filter are optimized in advance and stored in the internal engineering file, so that the parameters of Pk (f), fL(k)、fH(k) May be known in advance or calculated. So the rest of the equation can be calculated in advance except that n is dynamically determined by the aperture calculation module. Therefore, the real-time calculation amount is small, and the engineering implementation is easy. Especially in a plane wave imaging system, because emission is unfocused, energy is weak, a full array element works, and beam forming computation is large, optimal balance between image contrast and resolution needs to be better made. Therefore, the apodization coefficient in the beam synthesis is dynamically generated according to the setting of the demodulation filter, and better image quality can be obtained under the condition of not increasing real-time computing resources.
The embodiment of the invention provides a device for generating an apodization value for ultrasonic scanning, wherein the depth of the ultrasonic scanning is sequentially divided into N adjacent sections, and N is an integer greater than 1; the system comprises the following modules:
the first initialization module is used for acquiring a first apodization value of each array element in the kth section of the depth value, wherein k is an integer and is more than or equal to 1 and less than N;
the first synthesis beam generation module is used for multiplying the first echo signal of the array element by the first dynamic apodization value of the array element for each array element; carrying out time delay superposition on the first echo signals multiplied by all the array elements to obtain a first synthesized beam of a k-th section of depth value;
and the first processing module is used for obtaining a second apodization value of each array element at the k +1 th section of depth value based on the first synthesized beam.
Preferably, the processing module is specifically configured to: obtaining a demodulation filter type F (k), a demodulation bandwidth [ F [ ]L(k),FH(k)],FL(k) High-pass cut-off frequency, F, for the k-th segment depth valueH(k) A low-pass cut-off frequency which is the k-th section depth value;
obtaining a first spectral energy
Figure BDA0001518728460000081
Energy of the second spectrum
Figure BDA0001518728460000082
Figure BDA0001518728460000091
Figure BDA0001518728460000092
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<fL(k)<fH(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, PK1 is the demodulation filter bandwidth range [ fL(k),fH(k)]PK2 is a demodulation filter with a bandwidth range of [0, fs/2 ]]The spectral energy of (a).
Preferably, the following modules are also included: and the baseband signal acquisition module is used for demodulating and filtering the first synthesized beam to obtain a baseband signal.
Preferably, the synthetic beam generating module is specifically configured to: carrying out time delay processing on each array element; for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and superposing the first echo signals multiplied by all the array elements.
Preferably, the synthetic beam generating module is specifically configured to: for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and performing delay superposition on the first echo signals multiplied by all the array elements.
The second embodiment of the invention provides a method for generating an apodization value for ultrasonic scanning, wherein the depth of the ultrasonic scanning is sequentially divided into N adjacent segments, and N is an integer greater than 1; as shown in fig. 4, the method comprises the following steps:
step 401: acquiring a first apodization value of each array element in the kth section of depth value, wherein k is an integer and is more than or equal to 1 and less than N;
step 402: carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value;
step 403: and obtaining a second apodization value of each array element in the k-th section of the depth value based on the first synthesized beam.
Preferably, the obtaining a second apodization value of each array element at the k-th section of the depth value based on the first synthesized beam includes: obtaining a demodulation filter type F (k), a demodulation bandwidth [ F [ ]L(k),FE(k)],FL(k) High-pass cut-off frequency, F, for the k-th segment depth valueE(k) A low-pass cut-off frequency which is the k-th section depth value;
obtaining a first spectral energy
Figure BDA0001518728460000093
Energy of the second spectrum
Figure BDA0001518728460000094
Figure BDA0001518728460000101
Figure BDA0001518728460000102
WLUT () is a window function.
Preferably, the method further comprises the following steps: and demodulating and filtering the first synthesized beam to obtain a baseband signal.
Preferably, each array element is subjected to delay, superposition and apodization processing, including: carrying out time delay processing on each array element; for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and superposing the first echo signals multiplied by all the array elements.
Preferably, each array element is subjected to delay, superposition and apodization processing, including: for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and performing delay superposition on the first echo signals multiplied by all the array elements.
The embodiment of the invention also provides a device for generating the apodization value for ultrasonic scanning, wherein the depth of the ultrasonic scanning is sequentially divided into N adjacent sections, and N is an integer greater than 1; the system comprises the following modules:
the second initialization module is used for acquiring a first apodization value of each array element in the kth section of the depth value, wherein k is an integer and is more than or equal to 1 and less than N;
the second synthetic beam generating module is used for carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value;
and the second processing module is used for obtaining a second apodization value of each array element in the k-th section of the depth value based on the first synthesized beam.
It should be understood that although the present description refers to embodiments, not every embodiment contains only a single technical solution, and such description is for clarity only, and those skilled in the art should make the description as a whole, and the technical solutions in the embodiments can also be combined appropriately to form other embodiments understood by those skilled in the art.
The above-listed detailed description is only a specific description of a possible embodiment of the present invention, and they are not intended to limit the scope of the present invention, and equivalent embodiments or modifications made without departing from the technical spirit of the present invention should be included in the scope of the present invention.

Claims (4)

1. A method of generating an apodization value for an ultrasound scan, the depth of the ultrasound scan being divided sequentially into N adjacent segments, N being an integer greater than 1; the method is characterized by comprising the following steps:
acquiring a first apodization value of each array element in the kth section of depth value, wherein k is an integer and is more than or equal to 1 and less than k < N;
carrying out time delay, superposition and apodization processing on each array element so as to obtain a first synthetic beam of a k-th section of depth value;
obtaining a second apodization value of each array element at the k-th section of depth value based on the first synthesized beam;
the obtaining a second apodization value of each array element at the k-th section of the depth value based on the first synthesized beam includes:
obtaining a demodulation filter type F (k), a demodulation bandwidth [ f [ ]L(k),fH(k)],fL(k) High-pass cut-off frequency, f, for the k-th segment depth valueH(k) A low-pass cut-off frequency which is the k-th section depth value;
obtaining a first spectral energy
Figure FDA0002633192800000011
Energy of the second spectrum
Figure FDA0002633192800000012
Figure FDA0002633192800000015
Figure FDA0002633192800000014
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<fL(k)<fH(k)<fs/2,Pk(f) For demodulating the frequency response, P, of the filter over depth kk1For demodulation filters in the bandwidth range fL(k),fH(k)]Spectral energy of (P)k2For demodulation filters in the bandwidth range 0, fs/2]The spectral energy of (a).
2. The method of generating an apodization value according to claim 1, further comprising the steps of:
and demodulating and filtering the first synthesized beam to obtain a baseband signal.
3. The method of generating apodization values according to claim 1, wherein delaying, superimposing, and apodizing each array element comprises:
carrying out time delay processing on each array element; for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and superposing the first echo signals multiplied by all the array elements.
4. The method of generating apodization values according to claim 1, wherein delaying, superimposing, and apodizing each array element comprises:
for each array element, multiplying a first echo signal of the array element by a first dynamic apodization value of the array element; and performing delay superposition on the first echo signals multiplied by all the array elements.
CN201711397019.4A 2017-12-21 2017-12-21 Method and device for generating apodization value for ultrasonic scanning Active CN108113703B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711397019.4A CN108113703B (en) 2017-12-21 2017-12-21 Method and device for generating apodization value for ultrasonic scanning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711397019.4A CN108113703B (en) 2017-12-21 2017-12-21 Method and device for generating apodization value for ultrasonic scanning

Publications (2)

Publication Number Publication Date
CN108113703A CN108113703A (en) 2018-06-05
CN108113703B true CN108113703B (en) 2021-02-09

Family

ID=62231172

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711397019.4A Active CN108113703B (en) 2017-12-21 2017-12-21 Method and device for generating apodization value for ultrasonic scanning

Country Status (1)

Country Link
CN (1) CN108113703B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111329517B (en) * 2020-03-13 2023-05-12 京东方科技集团股份有限公司 Ultrasonic imaging method and device and storage medium
CN111882515B (en) * 2020-09-28 2020-12-29 深圳华声医疗技术股份有限公司 Ultrasonic signal processing method, ultrasonic signal processing apparatus, and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5068833A (en) * 1990-04-11 1991-11-26 Hewlett-Packard Company Dynamic control circuit for multichannel system
US6138513A (en) * 1999-01-09 2000-10-31 Barabash; Leonid S. Method and apparatus for fast acquisition of ultrasound images
CN101116622A (en) * 2006-08-02 2008-02-06 深圳迈瑞生物医疗电子股份有限公司 Beam unitized real-time computing technique for receiving apodized parameters and the device thereof

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100406099B1 (en) * 2001-09-05 2003-11-14 주식회사 메디슨 Ultrasound image forming apparatus and method using multi-level pulse compressor
CN1190753C (en) * 2003-02-28 2005-02-23 清华大学 Ultrasonic dynamic receiving apodization method based on delta-sigma transformation
US8744155B2 (en) * 2008-02-16 2014-06-03 University Of Virginia Patent Foundation Imaging or communications system utilizing multisample apodization and method
CN101940479B (en) * 2010-07-15 2012-10-03 重庆大学 Method and device for synthesizing ultrasonic imaging wave beam based on sectional dynamic apodization coefficient

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5068833A (en) * 1990-04-11 1991-11-26 Hewlett-Packard Company Dynamic control circuit for multichannel system
US6138513A (en) * 1999-01-09 2000-10-31 Barabash; Leonid S. Method and apparatus for fast acquisition of ultrasound images
CN101116622A (en) * 2006-08-02 2008-02-06 深圳迈瑞生物医疗电子股份有限公司 Beam unitized real-time computing technique for receiving apodized parameters and the device thereof

Also Published As

Publication number Publication date
CN108113703A (en) 2018-06-05

Similar Documents

Publication Publication Date Title
DE10244411B4 (en) Medical ultrasound imaging method and medical ultrasound imaging device
Ranganathan et al. Direct sampled I/Q beamforming for compact and very low-cost ultrasound imaging
JP2017534358A (en) Ultrasonic signal processing circuit and related apparatus and method
CN108324319B (en) System and method for distortion-free multibeam ultrasound receive beamforming
CN107180414B (en) System and method for reducing ultrasound speckle using harmonic compounding
US9107602B2 (en) Adaptive clutter filtering in an ultrasound system
CN108113703B (en) Method and device for generating apodization value for ultrasonic scanning
JP2003190157A (en) Ultrasonic diagnostic system
US10993701B2 (en) Ultrasonic imaging device
US6704437B1 (en) Noise estimation method and apparatus for noise adaptive ultrasonic image processing
DE10246353A1 (en) Pulse sequence responsive information generation method for ultrasound systems, involves filtering two beamformed signals with two tap filters and combining them in response to two sequential pulses
Ponnle et al. Suppression of grating lobe artifacts in ultrasound images formed from diverging transmitting beams by modulation of receiving beams
JP2008054800A (en) Ultrasonic diagnostic apparatus
Polichetti et al. A computationally efficient nonlinear beamformer based on p-th root signal compression for enhanced ultrasound b-mode imaging
CN107797111A (en) Robust multi-channel SAR signal reconstruction method under non-uniform scattering coefficient scene
Varray et al. A multi-frequency approach to increase the native resolution of ultrasound images
EP2641544A2 (en) Color flow imaging method, and ultrasonic device therefor
KR101646627B1 (en) Phase rotation beamformer and beamforming method for envelope error compensating
CN111012379B (en) Method and system for performing ultrasound imaging
JP7130008B2 (en) Method and apparatus for processing ultrasound signal data
CN111766306A (en) Full-focus imaging method and device
Wang et al. Generalized sidelobe canceller beamforming method for ultrasound imaging
Izadi et al. Weighted Capon beamformer combined with coded excitation in ultrasound imaging
Shen et al. A Flexible Speckle Reduction Strategy using Thomson’s Multitaper in High-order DMAS Beamforming
JP2020534903A5 (en)

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
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee after: Feiyinuo Technology Co.,Ltd.

Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee before: Feiyinuo Technology (Suzhou) Co.,Ltd.

Address after: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee after: Feiyinuo Technology (Suzhou) Co.,Ltd.

Address before: 215123 5th floor, building a, 4th floor, building C, No. 27, Xinfa Road, Suzhou Industrial Park, Jiangsu Province

Patentee before: VINNO TECHNOLOGY (SUZHOU) Co.,Ltd.