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
Energy of the second spectrum
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<f
L(k)<f
H(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, PK1 is the demodulation filter bandwidth range [ f
L(k),f
H(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
Energy of the second spectrum
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<f
L(k)<f
H(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, PK1 is the demodulation filter bandwidth range [ f
L(k),f
H(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.
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
Energy of the second spectrum
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<f
L(k)<f
H(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, PK1 is the demodulation filter bandwidth range [ f
L(k),f
H(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
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
Energy of the second spectrum
WLUT () is a window function, C is a constant, fs is a sampling rate of an input signal of the demodulation filter, 0<f
L(k)<f
H(k)<fs/2, Pk (f) is the frequency response of the demodulation filter at depth k, PK1 is the demodulation filter bandwidth range [ f
L(k),f
H(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
Energy of the second spectrum
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.