阅读说明：本技术 超声信号的处理方法和装置 (Ultrasonic signal processing method and device ) 是由 马克涛 王桂成 于琦 于 2018-06-19 设计创作，主要内容包括：本发明实施例提供一种超声信号的处理方法和装置。本发明的超声信号的处理方法,包括：根据所述初始化的第一迭代变量A<Sub>m</Sub>(0)和所述初始化的第一中间变量B<Sub>m</Sub>(0)之间的大小关系,对第一迭代变量A<Sub>m</Sub>(n)和第二迭代变量a<Sub>m</Sub>(n)进行迭代,确定阵元m发射的超声信号的回波信号对不同的探测距离的延时k<Sub>m</Sub>(n),对第一迭代变量A<Sub>m</Sub>(n)进行插值,分别根据插值变量<Image he="122" wi="276" file="DDA0001700354610000011.GIF" imgContent="drawing" imgFormat="GIF" orientation="portrait" inline="no"></Image>和<Image he="108" wi="299" file="DDA0001700354610000012.GIF" imgContent="drawing" imgFormat="GIF" orientation="portrait" inline="no"></Image>与第一中间变量B<Sub>m</Sub>(n)之间的大小关系,对延时k<Sub>m</Sub>(n+1)进行修正。本发明实施例可以实现降低波束合成过程中的延时计算所需的运算量,并且可以提升延时计算的精度值,提升波束合成质量。(The embodiment of the invention provides a method and a device for processing ultrasonic signals. The ultrasonic signal processing method of the invention comprises the following steps: a first iteration variable A according to said initialization m (0) And said initialized first intermediate variable B m (0) The magnitude relation between them, for the first iteration variable A m (n) and a second iteration variable a m (n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances m (n) for the first iteration variable A m (n) performing interpolation according to the interpolation variables And and a first intermediate variable B m (n) magnitude relation between, to delay k m (n +1) correction is performed. The embodiment of the invention can realize the reduction of the operation required by the delay calculation in the beam forming processAnd the accuracy value of delay calculation can be improved, and the beam synthesis quality is improved.)

1. A method of processing an ultrasound signal, comprising:

time delay variable k of initialization array element m_m(n) obtaining k of array element m_m(0)；

According to the k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0)；

A first iteration variable A according to said initialization_m(0) And said initialized first intermediate variable B_m(0) The magnitude relation between them, for the first iteration variable A_m(n) and a second iteration variable a_m(n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances_m(n)；

For the first iteration variable A_m(n) performing interpolation according to the interpolation variablesAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_m(n +1) correcting;

wherein, different values of m respectively correspond to different array elements, and n respectively takes 0 to X.

2. According to claim1, the method according to which said first iteration variable a is initialized_m(0) And said initialized first intermediate variable B_m(0) The magnitude relation between them, for the first iteration variable A_m(n) and a second iteration variable a_m(n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances_m(n) comprising:

according to the formula delta_m(n)＝B_m(n)-A_m(n) calculating delta_m(n)；

If delta_m(n) < 0, then a_m(n+1)＝a_m(n)，A_m(n+1)＝A_m(n)；

If delta_m(n) is greater than or equal to 0, then a_m(n+1)＝a_m(n)+2，A_m(n+1)＝A_m(n)+a_m(n)；

k_m(n)＝(a_m(n)-1-2n)/2。

3. The method of claim 2, wherein the pair of first iteration variables A_m(n) performing interpolation, comprising:

according to the formulaComputing

According to the formulaComputing

4. Method according to claim 3, characterized in that said respective function is based on an interpolated variableAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_m(n +1) performing a correction including:

when delta_m(n) < 0, andthen, k (n +1) + corr is set to k according to the formula k' (n +1)_m(n +1) correcting;

when delta_m(n) < 0, andthen, k is added to k according to the formula k' (n +1) ═ k (n +1)_m(n +1) correcting;

when delta_m(n) is not less than 0, andthen, k (n +1) + corr is set to k according to the formula k' (n +1)_m(n +1) correcting;

when delta_m(n) is not less than 0, andthen, k is added to k according to the formula k' (n +1) ═ k (n +1)_m(n +1) correcting;

wherein corr is 0.5.

5. Method according to any of claims 1 to 4, wherein the delay variable k of the initialization array element m_m(n) obtaining k of array element m_m(0) The method comprises the following steps:

acquiring type information of a probe;

when the type information of the probe is a linear type, the method is based on a formulaDelay variable k to array element m_m(n) initializing;

when the type information of the probe is a convex array type, the type information is determined according to a formulaDelay variable k to array element m_m(n) initializing;

wherein the content of the first and second substances,indicating a rounding down.

6. The method according to claim 5, wherein when the type information of the probe is a linear type, the k is based on the k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0) The method comprises the following steps:

according to formula A_m(0)＝k_m ²(0) Delay variable A for array element m_m(n) initializing;

according to the formula a_m(0)＝2k_m(0) +1 delay variable a for array element m_m(n) initializing;

according to formula B_m(0)＝x_m ²/d²Delay variable B for array element m_m(n) initializing;

according to formula b_m(0)＝2x_msin theta/d +1 time delay variable b of array element m_m(n) initializing.

7. The method according to claim 5, wherein when the type information of the probe is a convex array type, the k is based on the k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And a second of the initializationIntermediate variable b_m(0) The method comprises the following steps:

according to the formulaDelay variable A for array element m_m(n) initializing;

according to the formula a_m(0)＝2k_m(0) +1 delay variable a for array element m_m(n) initializing;

according to the formulaDelay variable B for array element m_m(n) initializing;

according to formula b_m(0)＝2a_v(m) sin (theta + gamma (m))/d +1 delay variable b of array element m_m(n) initializing.

8. Method according to any of claims 1 to 4, wherein the delay variable k of the initialization array element m_m(n) obtaining k of array element m_m(0) The method comprises the following steps:

acquiring k of the array element m corresponding to the deflection angle according to a preset data table and the deflection angle of the focus point_m(0)；

And the preset data table comprises initialized values corresponding to different array element positions and deflection angles.

9. The method according to any one of claims 1 to 4, further comprising:

determining the range of the effective array elements according to the range of an included angle between the central line and the normal;

and determining the value range of m according to the range of the effective array elements.

10. The method of claim 9, wherein determining the range of valid array elements from the range of angles between the centerline and the normal comprises:

when the type of probe isWhen the information is of linear type, according to the formulaComputingAccording to the formulaComputingThe range of the effective array elements is

When the type information of the probe is a convex array type, according to a formulaComputingAccording to the formulaComputingThe range of the effective array elements is

And theta is an included angle between the central line and the normal line.

11. An apparatus for processing an ultrasound signal, comprising:

a memory for storing a computer program;

a processor for executing the computer program to implement the method of any one of claims 1 to 10.

12. A computer storage medium, comprising: the computer storage medium is for storing a computer program which when executed is for implementing the method of any one of claims 1 to 10.

Technical Field

The embodiment of the invention relates to a signal processing technology, in particular to a method and a device for processing ultrasonic signals.

Background

With the continuous development of scientific technology, ultrasound systems are used more and more widely in medicine. The most important part is beam forming, and the key of beam forming is the calculation of delay on different detection distances of different array elements.

A plurality of array elements of the ultrasonic probe array emit detection sound waves, because the positions of different array elements relative to the same tissue point are different, the length of a path for the infrasonic waves to reach tissues is different, the reflection of different tissues to the sound waves is also different, the sound waves are finally reflected back, the time lengths for reaching the position of a focus point are different, different sound wave beams are superposed to form a focus point, and a focus image is formed by the focus points. Since the sound waves for forming one focusing point do not return at the same time, the respective delays need to be calculated, and the sound waves belonging to the feedback signal of the current sound wave are found to form the focusing point. The beam synthesis is to delay and apodize the sound wave to form a focused image.

For the beam forming, focusing is performed for each detection distance (beam forming), and calculating the delay value results in a large amount of calculation and low calculation accuracy.

Disclosure of Invention

The embodiment of the invention provides a method and a device for processing ultrasonic signals, which are used for reducing the computation amount required by delay calculation in a beam forming process and improving the accuracy value of the delay calculation.

In a first aspect, an embodiment of the present invention provides a method for processing an ultrasound signal, including:

time delay variable k of initialization array element m_m(n) obtaining k of array element m_m(0)；

According to the k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0)；

A first iteration variable A according to said initialization_m(0) And said initialized first intermediate variable B_m(0) The magnitude relation between them, for the first iteration variable A_m(n) and a second iteration variable a_m(n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances_m(n)；

For the first iteration variable A_m(n) performing interpolation according to the interpolation variablesAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_m(n +1) correcting;

wherein, different values of m respectively correspond to different array elements, and n respectively takes 0 to X.

In a second aspect, an embodiment of the present invention provides an apparatus for processing an ultrasound signal, including:

a memory for storing a computer program;

a processor for executing the computer program to implement the method according to the first aspect.

In a third aspect, an embodiment of the present invention provides a computer storage medium, including: the computer storage medium is for storing a computer program which, when executed, is for implementing the method as described in the first aspect.

The ultrasonic signal processing method and device provided by the embodiment of the invention initialize the delay variable k of the array element m_m(n) obtaining k of array element m_m(0) According to said k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0) According to the initialized first iteration variable A_m(0) And said initialized first intermediate variable B_m(0) The magnitude relation between them, for the first iteration variable A_m(n) and a second iteration variable a_m(n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances_m(n) for the first iteration variable A_m(n) performing interpolation according to the interpolation variablesAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_mAnd (n +1) correcting, delaying the echo signal of the ultrasonic signal by the corrected delay value so as to perform space composite imaging subsequently, thereby reducing the computation required by delay calculation in the beam synthesis process, improving the precision value of delay calculation and improving the beam synthesis quality.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.

Fig. 1 is a schematic view of an application scenario of the processing method of an ultrasonic signal according to the present invention;

FIG. 2 is a flow chart of a first embodiment of a method for processing an ultrasound signal according to the present invention;

FIG. 3 is a flowchart of a second embodiment of the method for processing ultrasonic signals according to the present invention;

FIG. 4 is a schematic diagram of the delay calculation when the type information of the probe is a linear type;

FIG. 5 is a schematic diagram of the delay calculation when the type information of the probe is a convex array type;

FIG. 6 is a flow chart of a third embodiment of the method for processing ultrasonic signals according to the present invention;

FIG. 7 is a schematic structural diagram of a first embodiment of an apparatus for processing ultrasonic signals according to the present invention;

FIG. 8 is a schematic structural diagram of a second embodiment of an apparatus for processing ultrasonic signals according to the present invention;

fig. 9 is a schematic structural diagram of a third embodiment of an apparatus for processing an ultrasonic signal according to the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

Fig. 1 is a schematic view of an application scenario of the method for processing an ultrasound signal according to the present invention, as shown in fig. 1, the application scenario includes a probe and a processing apparatus for an ultrasound signal. The ultrasonic signal processing device is connected with the probe.

The ultrasonic signal processing device can execute the ultrasonic signal processing method of the invention to reduce the computation amount required by delay calculation in the beam forming process, and can improve the precision value of the delay calculation and improve the beam forming quality.

The probe is used for transmitting and receiving ultrasonic signals, namely echo signals of the ultrasonic signals can be received. The ultrasonic signal processing device can perform time delay calculation on the echo signal of the ultrasonic signal so as to perform space composite imaging subsequently and output an image.

The processing device of the ultrasonic signal may be a chip, or may be a functional module in the chip, and the processing device of the ultrasonic signal may be disposed separately from the display, or may be disposed together with the display.

The method for processing ultrasonic signals according to the present invention will be specifically explained below using several specific examples.

Fig. 2 is a flowchart of a first embodiment of a method for processing an ultrasound signal according to the present invention, and as shown in fig. 2, the method of this embodiment may include:

step 101, initializing a delay variable k of an array element m_m(n) obtaining k of array element m_m(0)。

Wherein, different values of m respectively correspond to different array elements, n respectively takes 0 to X, and n can indicate different detection distances. Taking the array element m of a linear type probe as an example, thex_mIs the coordinate of the array element m. Taking array element m of a convex array type probe as an example, the methoda_vAnd (m) is the curvature radius of the array element m. Wherein the content of the first and second substances,indicating a rounding down. X is any integer.

Step 102, according to the k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0)。

Wherein A is_m(0) For the initialized first iteration variable, A_m(n) is a first iteration variable, A_m(0) Is A_m(n) a value where n is 0_mThe specific expression of (n) is the following formula (1). a is_m(0) For initialising the second iteration variable, a_m(n) is a second iteration variable, a_m(0) Is a_m(n) a value where n is 0_mThe specific expression of (n) is the following formula (2). B is_m(0) For the initialized first intermediate variable, B_m(n) is a first intermediate variable, B_m(0) Is B_m(n) a value where n is 0, B_mThe specific expression of (n) is the following formula (3). b_m(0) For initialising the second intermediate variable, b_m(n) is a second intermediate variable，b_m(0) Is b_m(n) a value where n is 0, b_mThe specific expression of (n) is the following formula (4).

A_m(n)＝n²+2k(n)n+k²(n) (1)

B_m(n)＝n²+α_mn+β_m (2)

a_m(n)＝2n+2k(n)+1 (3)

b_m(n)＝2n+1+α_m (4)

Wherein, taking the array element m of the linear probe as an example,taking the array element m of the convex array type probe as an example,

can be based on the above equations (1) and k_m(0) Determining an initialized first iteration variable A_m(0) According to the formulae (3) and k_m(0) Determining an initialized second iteration variable a_m(0) According to the formulae (2) and k_m(0) Determining an initialized first intermediate variable B_m(0) According to the formulae (4) and k_m(0) Determining an initialized second intermediate variable b_m(0)。

103, according to the initialized first iteration variable A_m(0) And said initialized first intermediate variable B_m(0) The magnitude relation between them, for the first iteration variable A_m(n) and a second iteration variable a_m(n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances_m(n)。

Specifically, in the embodiment of the present invention, the first iteration variable a is processed in an iteration manner_m(n) and a second iteration variable a_m(n) iterating, i.e. using A_m(0) And a_m(0) Determination of A_m(1) And a_m(1) Using A_m(1) And a_m(1) Determination of A_m(2) And a_m(2) And so on until determining A_m(X) and a_m(X). And k of n at different values can be determined according to the formula (3) in the iterative process_m(n)。

104, for the first iteration variable A_m(n) performing interpolation according to the interpolation variablesAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_m(n +1) correction is performed.

Specifically, this embodiment may also apply to the first iteration variable a_m(n) performing interpolation to determine an interpolation variableAndby comparing A_m(n) and B_m(n)、And B_m(n)、And B_m(n) magnitude relation between, to delay k_mAnd (n +1) correcting, and delaying the echo signal of the ultrasonic signal by the corrected delay value so as to perform space composite imaging subsequently.

In this embodiment, the delay variable k of the array element m is initialized_m(n) obtaining k of array element m_m(0) According to said k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0) According to the initialized first iteration variable A_m(0) And said initialized first intermediate variable B_m(0) Magnitude relationship between, for the firstIteration variable A_m(n) and a second iteration variable a_m(n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances_m(n) for the first iteration variable A_m(n) performing interpolation according to the interpolation variablesAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_mAnd (n +1) correcting, delaying the echo signal of the ultrasonic signal by the corrected delay value so as to perform space composite imaging subsequently, thereby reducing the computation required by delay calculation in the beam synthesis process, improving the precision value of delay calculation and improving the beam synthesis quality.

The following describes in detail the technical solution of the embodiment of the method shown in fig. 2, using several specific embodiments.

Fig. 3 is a flowchart of a second embodiment of the method for processing an ultrasound signal according to the present invention, and as shown in fig. 3, the method of this embodiment may include:

step 201, initializing delay variable k of array element m_m(n) obtaining k of array element m_m(0)。

Step 202, according to the k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0)。

Step 2031, according to the formula delta_m(n)＝B_m(n)-A_m(n) calculating delta_m(n)。

Step 2032, if δ_m(n)<0, then a is determined using equation (5)_m(n +1), determining A according to formula (6)_m(n+1)。

a_m(n+1)＝a_m(n) (5)

A_m(n+1)＝A_m(n) (6)

Step 2033, if δ_m(n) ≧ 0, a is determined using equation (7)_m(n +1), determining A according to formula (8)_m(n+1)。

a_m(n+1)＝a_m(n)+2 (7)

A_m(n+1)＝A_m(n)+a_m(n) (8)

Step 2034, calculating k according to equation (9)_m(n)。

k_m(n)＝(a_m(n)-1-2n)/2 (9)

Step 2041, calculate according to formula (10)Calculation according to equation (11)

Step 2042, when delta_m(n)<0, andthen, k is calculated according to the formula (12)_m(n +1) correction is performed.

k′(n+1)＝k(n+1)+corr (12)

Step 2043, when delta_m(n)<0, andthen, k is calculated according to the formula (13)_m(n +1) correction is performed.

k′(n+1)＝k(n+1) (13)

Step 2044, when delta_m(n) is not less than 0, andthen, k is calculated according to the formula (12)_m(n +1) correction is performed.

Step 2045, when delta_m(n) is not less than 0, andthen, k is calculated according to the formula (13)_m(n +1) correction is performed.

Wherein corr is 0.5. k' (n +1) is a corrected value of k (n + 1).

It should be noted that, by performing step 2041 to step 2045 once, the delay calculation accuracy can be improved by 2 times, and specifically, the delay calculation accuracy before correction isT_uIs the sampling period of the ultrasonic signal. By performing step 2041 to step 2045 once, the delay calculation accuracy can be improved

Optionally, after performing step 2041 to step 2045 once, the formula may be further used(when is (n)<0 time) or the formula (when δ (n) ≧ 0), for δ_m(n) updating. According to the formulaTo a_m(n) updating. According to the formulaThe corr is updated. The updated δ may be used_m(n)、a_m(n) is selected fromAnd corr, by the way as above step 2041 to step 2044 for k_m(n +1) performing a re-correction. The delay calculation precision is improved toN is the number of times step 2041 is performed.

In this embodiment, the delay variable k of the array element m is initialized_m(n) obtaining k of array element m_m(0) According to said k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0) According to the initialized first iteration variable A_m(0) And said initialized first intermediate variable B_m(0) The magnitude relation between them, for the first iteration variable A_m(n) and a second iteration variable a_m(n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances_m(n) for the first iteration variable A_m(n) performing interpolation according to the interpolation variablesAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_mAnd (n +1) correcting, delaying the echo signal of the ultrasonic signal by the corrected delay value so as to perform space composite imaging subsequently, thereby reducing the computation required by delay calculation in the beam synthesis process, improving the precision value of delay calculation and improving the beam synthesis quality.

Fig. 4 is a schematic diagram of delay calculation when the type information of the probe is a linear type, and as shown in fig. 4, the probe includes a plurality of array elements, each array element is represented by a small square, and the array elements are arranged along a straight line. Taking 0 point as the position of a receiving central line and theta as the deflection angle of the receiving central line, array elements on two sides simultaneously receive echo signals of ultrasonic signals,for x⁺For the vibration element(s), the deflection angle (i.e. the included angle between the scanning beam and the vibration element) is larger than 90 degrees; for x^-The deflection angle is less than 90 deg.

T_uIs the sampling period of the ultrasonic signal, c₀Is the propagation velocity of the ultrasonic signal, theta is the receive centerline declination (as shown).

Parameters such as r, l, etc. referred to in fig. 4 are all quantified using d (d represents the distance length of half of the round-trip acoustic wave in one sampling period, being the distance length of the round-trip or the return only), where n-r (nT) is r (nT)_u) D is an integer; k (n) ═ l (nT)_u) And/d, quantified using d.

The above embodiment is made when the type information of the probe is a linear type as shown in fig. 4

And according to formula A_m(0)＝k_m ²(0) Delay variable A for array element m_m(n) initializing according to the formula a_m(0)＝2k_m(0) +1 delay variable a for array element m_m(n) initializing according to formula B_m(0)＝x_m ²/d²Delay variable B for array element m_m(n) initializing according to formula b_m(0)＝2x_msin theta/d +1 time delay variable b of array element m_m(n) initializing.

It should be noted that, in another implementation manner, according to a preset data table and a deflection angle of a focus point, k of an array element m corresponding to the deflection angle is obtained_m(0) (ii) a And the preset data table comprises initialized values corresponding to different array element positions and deflection angles. I.e. k can be determined by means of a table look-up_m(0)。

FIG. 5 is a schematic diagram of delay calculation when the type information of the probe is a convex array type, as shown in FIG. 5, a plurality of array elements of the probe are not in the same arrayOn the horizontal line, there is an angle difference γ, so θ in the above linear type becomes θ + γ, and x becomes a_vWherein a is_vExpressed as values of the properties of the probe. For different array elements m (m E [ -K, K)]) There are different θ + γ (m) and a_v(m)。m>0 corresponds to x in the linear array⁺Part of (a), m<0 corresponds to x in the linear array^-Part (c) of (a).

When the type information of the probe is a convex array type as shown in fig. 5, the above embodiment

According to the formulaDelay variable A for array element m_m(n) initializing according to the formula a_m(0)＝2k_m(0) +1 delay variable a for array element m_m(n) initializing according to the formulaDelay variable B for array element m_m(n) initializing according to formula b_m(0)＝2a_v(m) sin (theta + gamma (m))/d +1 delay variable b of array element m_m(n) initializing.

It should be noted that, in another implementation manner, according to a preset data table and a deflection angle of a focus point, k of an array element m corresponding to the deflection angle is obtained_m(0) (ii) a And the preset data table comprises initialized values corresponding to different array element positions and deflection angles. I.e. k can be determined by means of a table look-up_m(0)。

Fig. 6 is a flowchart of a third embodiment of a processing method of an ultrasound signal according to the present invention, as shown in fig. 6, on the basis of any one of the above embodiments, the method of the present embodiment may further include:

and 301, determining the range of the effective array element according to the range of the included angle between the central line and the normal line.

In particular, when the type of probe is informationIn the linear type, according to the formulaComputingAccording to the formulaComputingThe range of the effective array elements is

When the type information of the probe is a convex array type, according to a formulaComputingAccording to the formulaComputingThe range of the effective array elements is

Where θ is an angle between the center line and the normal, which may also be referred to as a receiving center line declination.

And 302, determining the value range of m according to the range of the effective array elements.

Specifically, x in the above embodiment_mIs required to be atWithin the range. A in the above-mentioned embodiment_v(m) is required to be inWithin the range.

In this embodiment, the range of the effective array element is determined according to the included angle range between the central line and the normal line, and the value range of m is determined according to the range of the effective array element, so that the calculation of the invalid array element can be removed on the basis of not influencing the subsequent imaging, and the operation efficiency can be further improved.

Fig. 7 is a schematic structural diagram of a first embodiment of an apparatus for processing an ultrasonic signal according to the present invention, and as shown in fig. 7, the apparatus of the present embodiment may include: an initialization module 11, an iterative computation module 12 and a modification module 13, wherein the initialization module 11 is configured to initialize a delay variable k of an array element m_m(n) obtaining k of array element m_m(0) (ii) a The initialization module 11 is further configured to initialize the k_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0) (ii) a The iterative computation module 12 is configured to compute the first iteration variable a according to the initialization_m(0) And said initialized first intermediate variable B_m(0) The magnitude relation between them, for the first iteration variable A_m(n) and a second iteration variable a_m(n) iteration is carried out to determine the time delay k of the echo signal of the ultrasonic signal emitted by the array element m to different detection distances_m(n); the modification module 13 is used for modifying the first iteration variable A_m(n) performing interpolation according to the interpolation variablesAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_m(n +1) correcting; wherein, different values of m respectively correspond to different array elements, and n respectively takes 0 to X.

The iterative calculation module 12 is used for calculating the formulaδ_m(n)＝B_m(n)-A_m(n) calculating delta_m(n)；

If delta_m(n)<0, then a_m(n+1)＝a_m(n),A_m(n+1)＝A_m(n)；

If delta_m(n) is greater than or equal to 0, then a_m(n+1)＝a_m(n)+2,A_m(n+1)＝A_m(n)+a_m(n)；

k_m(n)＝(a_m(n)-1-2n)/2。

The modification module 13 is used for modifying the first iteration variable A_m(n) performing interpolation, comprising: according to the formulaComputingAccording to the formula Computing

The correction module 13 is used for respectively according to the interpolation variablesAndand a first intermediate variable B_m(n) magnitude relation between, to delay k_m(n +1) performing a correction including:

when delta_m(n)<0, andthen, k (n +1) + corr is set to k according to the formula k' (n +1)_m(n +1) correcting;

when delta_m(n)<0, andthen, k is added to k according to the formula k' (n +1) ═ k (n +1)_m(n +1) correcting;

when delta_m(n) is not less than 0, andthen, k (n +1) + corr is set to k according to the formula k' (n +1)_m(n +1) correcting;

when delta_m(n) is not less than 0, andthen, k is added to k according to the formula k' (n +1) ═ k (n +1)_m(n +1) correcting;

wherein corr is 0.5.

The initialization module 11 is used for initializing a delay variable k of an array element m_m(n) obtaining k of array element m_m(0) The method comprises the following steps: acquiring type information of a probe; when the type information of the probe is a linear type, the method is based on a formulaDelay variable k to array element m_m(n) initializing; when the type information of the probe is a convex array type, the type information is determined according to a formulaDelay variable k to array element m_m(n) initializing; wherein the content of the first and second substances,indicating a rounding down.

The initialization module 11 is configured to initialize the probe according to the k when the type information of the probe is a linear type_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0) The method comprises the following steps:

according to formula A_m(0)＝k_m ²(0) Delay variable A for array element m_m(n) initializing;

according to the formula a_m(0)＝2k_m(0) +1 delay variable a for array element m_m(n) initializing;

according to formula B_m(0)＝x_m ²/d²Delay variable B for array element m_m(n) initializing;

according to formula b_m(0)＝2x_msin theta/d +1 time delay variable b of array element m_m(n) initializing.

When the type information of the probe is a convex array type, the initialization module 11 is used for determining the k according to the_m(0) Determining an initialized first iteration variable A_m(0) Initialized second iteration variable a_m(0) Initialized first intermediate variable B_m(0) And an initialized second intermediate variable b_m(0) The method comprises the following steps:

according to the formulaDelay variable A for array element m_m(n) initializing;

according to the formula a_m(0)＝2k_m(0) +1 delay variable a for array element m_m(n) initializing;

according to the formulaDelay variable B for array element m_m(n) initializing;

according to formula b_m(0)＝2a_v(m) sin (theta + gamma (m))/d +1 delay variable b of array element m_m(n) initializing.

The delay variable k of the initialization array element m_m(n) obtaining k of array element m_m(0) The method comprises the following steps:

acquiring k of the array element m corresponding to the deflection angle according to a preset data table and the deflection angle of the focus point_m(0)；

And the preset data table comprises initialized values corresponding to different array element positions and deflection angles.

The apparatus of this embodiment may be configured to implement the technical solutions of the above method embodiments, and the implementation principles and technical effects are similar, which are not described herein again.

Fig. 8 is a schematic structural diagram of a second embodiment of an apparatus for processing an ultrasonic signal according to the present invention, and as shown in fig. 8, the apparatus of the present embodiment may further include, on the basis of the apparatus structure shown in fig. 7: an effective array element determining module 14, wherein the effective array element determining module 14 is configured to determine a range of an effective array element according to an included angle range between the central line and the normal; and determining the value range of m according to the range of the effective array elements.

Wherein, the range of the effective array element is determined according to the range of the included angle between the central line and the normal line, and the method comprises the following steps:

when the type information of the probe is a linear type, according to a formulaComputingAccording to the formulaComputingThe range of the effective array elements is

When the type information of the probe is a convex array type, according to a formulaComputingAccording to the formulaComputingThe range of the effective array elements is

And theta is an included angle between the central line and the normal line.

The apparatus of this embodiment may be configured to implement the technical solutions of the above method embodiments, and the implementation principles and technical effects are similar, which are not described herein again.

Fig. 9 is a schematic structural diagram of a third embodiment of an apparatus for processing an ultrasound signal according to the present invention, and as shown in fig. 9, the apparatus of the present embodiment may include: a memory 21 and a processor 22, wherein the memory 21 is used for storing a computer program, and the processor 22 is used for executing the computer program to implement the method described in the above embodiments.

The apparatus of this embodiment may be configured to implement the technical solutions of the above method embodiments, and the implementation principles and technical effects are similar, which are not described herein again.

An embodiment of the present invention further provides a computer storage medium, where the computer storage medium is used to store a computer program, and the computer program is used to implement the method according to any of the above embodiments when executed. The implementation principle and the technical effect are similar, and the detailed description is omitted here.

Those of ordinary skill in the art will understand that: all or a portion of the steps of implementing the above-described method embodiments may be performed by hardware associated with program instructions. The program may be stored in a computer-readable storage medium. When executed, the program performs steps comprising the method embodiments described above; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.

Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.