Phase correction table optimization method, device, ultrasonic equipment and storage medium

Abstract

The present invention discloses a phase correction table optimization method, device, ultrasonic equipment and storage medium. The method includes: determining the effective focal area corresponding to this transmission based on a preset focus; optimizing the effective focal area corresponding to this transmission using the coherent area function corresponding to the last beam synthesis, and determining the target focal area corresponding to this transmission; using any focal point in the target focal area as a virtual array element, and determining the phase correction table corresponding to this transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point. This scheme optimizes the effective focal area corresponding to this transmission using the coherent area function corresponding to the last beam synthesis, and can modulate a better focal area distribution result, and optimize the process of updating the phase correction table, without the need for advance simulation, and does not occupy storage resources, which can solve the shortcomings of ideal assumptions in the human body environment, and achieve seamless connection and image quality improvement in conventional ultrasound processes.

Description

Technical Field

[0001] The present invention relates to the field of ultrasonic imaging technology, and in particular to a phase correction table optimization method, device, ultrasonic equipment and storage medium. Background Art

[0002] The principle of medical ultrasound imaging is very mature. The multi-element probe is used to receive the echo of the transmitted ultrasound after it is reflected by the tissue in the human body. The flight time of the ultrasound propagation is calculated based on the geometric relationship to obtain the echo intensity at the corresponding position for imaging, thereby obtaining the characterization of the tissue structure. However, in this process, there are many specific technical details that will have a significant impact on the imaging quality, among which the accuracy of the flight time is the most challenging. Because in general ultrasound equipment, for the sake of system simplicity, the flight time is estimated based on the geometric distance between the probe position and the position to be imaged. However, the wave nature of ultrasound determines that its propagation time cannot strictly follow the geometric relationship, and the non-uniform characteristics of normal human tissue will bring many disturbances to the propagation of ultrasound, thereby further reducing the accuracy of the flight time calculated by the geometric relationship.

[0003] On top of this, modern ultrasonic equipment widely uses multi-beam transmission technology, which means that each time an ultrasonic wave is transmitted, multiple echo lines can be obtained through beam synthesis, and the multiple echo lines are often arranged symmetrically with the transmission center line of the ultrasonic wave transmission position. After one transmission, the transmission position of the next ultrasonic wave will be moved, and then multiple new echo lines will be obtained. Since the offset of each transmission position is very small, some echo line positions in the multiple groups of echo lines obtained by multiple transmissions often overlap. In beam synthesis, the echo lines at the same position obtained by these multiple transmissions can be coherently superimposed to improve the signal-to-noise ratio of the echo lines. However, there is a very important influencing factor in this superposition process, that is, whether the phase of the same echo line obtained by multiple transmissions of ultrasonic waves is consistent? Due to the coherence of ultrasonic imaging, if the signal phases of different transmissions at the same position are not aligned, the result of simple superposition may be signal cancellation, thereby failing to achieve the signal-to-noise ratio improvement effect of the application. Therefore, it is urgent to provide a phase correction table for phase correction so that the signal phases of different transmissions at the same position are aligned.

[0004] The existing phase correction table either uses ultrasonic simulation software (such as Field II) to simulate the sound field under multiple transmission positions, and directly obtains the phase correction value under an ideal uniform medium through sound field analysis to determine the phase correction table. This method is too ideal and the effect is not robust enough; or it estimates a phase correction value by performing coherence analysis on the echo signal in beamforming, so that the phase correction table can be directly calculated and determined in beamforming. This method is too cumbersome and occupies a large amount of system resources. Summary of the invention

[0005] The embodiments of the present invention provide a phase correction table optimization method, device, ultrasonic equipment and storage medium to solve the problems that the existing phase correction table determination method is too ideal, the effect is not robust enough, or the process is too complicated and occupies a large amount of system resources.

[0006] A phase correction table optimization method, comprising:

[0007] Based on the preset focus, determine the effective focus area corresponding to this launch;

[0008] The coherent area function corresponding to the last beam synthesis is used to optimize the effective focal area corresponding to this transmission, and determine the target focal area corresponding to this transmission;

[0009] Any focal point in the target focal area is used as a virtual array element, and a phase correction table corresponding to this transmission is determined based on the geometric relationship between the virtual array element, the imaging point and the correction point.

[0010] Preferably, the determining of the effective focus area corresponding to the current emission based on the preset focus includes:

[0011] Based on the focal area extension range corresponding to the previous beamforming and the theoretical focal area extension range, the focal area extension range corresponding to the current launch is determined;

[0012] Based on the preset focus and the focal zone extension range corresponding to this emission, the effective focal area corresponding to this emission is determined.

[0013] Preferably, determining the phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point includes:

[0014] Determining a first distance between the virtual array element and an imaging point;

[0015] Determining a second distance between the virtual array element and the calibration point;

[0016] Determine an original sound path difference corresponding to the virtual array element based on a first distance and a second distance corresponding to the virtual array element;

[0017] In the target focus area, integrating the original sound path difference corresponding to the virtual array element to determine the target sound path difference corresponding to the virtual array element;

[0018] Based on the target sound path differences corresponding to all the virtual array elements, a phase correction table corresponding to this transmission is determined.

[0019] Preferably, after determining the phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point, the phase correction table optimization method further includes:

[0020] Determine a neighboring area corresponding to a preset focus;

[0021] Performing time-of-flight sampling based on the original echo signal corresponding to the adjacent area to determine the full-channel echo signal corresponding to the adjacent area;

[0022] Performing coherence analysis on the full-channel echo signal corresponding to the adjacent area to determine a target coherence factor function corresponding to the adjacent area;

[0023] Based on the target coherence factor function corresponding to the adjacent area, the focal area extension range and the coherence area function corresponding to the current beamforming are determined.

[0024] Preferably, performing coherence analysis on the full-channel echo signal corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area includes:

[0025] Calculating the coherence factor of the full-channel echo signal corresponding to the adjacent area to determine the original coherence factor function corresponding to the adjacent area;

[0026] The original coherence factor function corresponding to the adjacent area is smoothed and filtered to determine the target coherence factor function corresponding to the adjacent area.

[0027] Preferably, determining the focal area extension range and the coherent area function corresponding to the current beamforming based on the target coherence factor function corresponding to the adjacent area includes:

[0028] Determining a neighborhood focal area corresponding to the current beamforming based on a target coherence factor function corresponding to the neighborhood area and a preset proportional coefficient;

[0029] Determine the focal zone extension range corresponding to this beamforming based on the real coordinates and theoretical center coordinates corresponding to the neighborhood focal area;

[0030] Function fitting is performed on the target coherence factor function corresponding to the neighborhood focal area to determine the coherence area function of this beamforming.

[0031] Preferably, after determining the phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point, the phase correction table optimization method further includes:

[0032] Obtaining the original echo signal;

[0033] Using the phase correction table corresponding to this transmission, the original echo signal is corrected for the flight time, and the target echo signal is determined based on the corrected flight time;

[0034] Performing channel compounding and multi-beam compounding on the target echo signal to determine a target composite signal;

[0035] A target ultrasound image is determined based on the target composite signal.

[0036] A phase correction table optimization device, comprising:

[0037] An effective focus area determination module, used to determine the effective focus area corresponding to this emission based on a preset focus;

[0038] A target focal area determination module is used to optimize the effective focal area corresponding to this transmission by using the coherent area function corresponding to the previous beamforming, and determine the target focal area corresponding to this transmission;

[0039] The phase correction table determination module is used to take any focal point in the target focal area as a virtual array element, and determine the phase correction table corresponding to this transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point.

[0040] Preferably, the effective focus area determination module includes:

[0041] A first zone extension range determination unit, configured to determine a focal zone extension range corresponding to this transmission based on a focal zone extension range corresponding to the previous beamforming and a theoretical focal zone extension range;

[0042] The effective focus area determination unit is used to determine the effective focus area corresponding to this emission based on the preset focus and the focal area extension range corresponding to this emission.

[0043] Preferably, the phase correction table determination module comprises:

[0044] A first distance determining unit, used to determine a first distance between the virtual array element and the imaging point;

[0045] A second distance determining unit, used to determine a second distance between the virtual array element and the calibration point;

[0046] an original sound path difference determining unit, configured to determine the original sound path difference corresponding to the virtual array element based on the first distance and the second distance corresponding to the virtual array element;

[0047] A target sound path difference determination unit, configured to perform integration processing on the original sound path difference corresponding to the virtual array element in the target focus area to determine the target sound path difference corresponding to the virtual array element;

[0048] The phase correction table determining unit is used to determine the phase correction table corresponding to the current transmission based on the target sound path differences corresponding to all the virtual array elements.

[0049] Preferably, the phase correction table optimization device further comprises:

[0050] A neighboring area determination module, used to determine a neighboring area corresponding to a preset focus;

[0051] A full-channel echo signal determination module, configured to perform time-of-flight sampling based on the original echo signal corresponding to the adjacent area to determine the full-channel echo signal corresponding to the adjacent area;

[0052] A target coherence factor function determination module is used to perform coherence analysis on the full-channel echo signal corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area;

[0053] The current beamforming information determination module is used to determine the focal area extension range and the coherent area function corresponding to the current beamforming based on the target coherence factor function corresponding to the adjacent area.

[0054] Preferably, the target coherence factor function determination module comprises:

[0055] An original coherence factor function determining unit, configured to calculate the coherence factor of the full-channel echo signal corresponding to the adjacent area, and determine the original coherence factor function corresponding to the adjacent area;

[0056] The target coherence factor function determination unit is used to perform smoothing filtering on the original coherence factor function corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area.

[0057] Preferably, determining the focal area extension range and the coherent area function corresponding to the current beamforming based on the target coherence factor function corresponding to the adjacent area includes:

[0058] A neighborhood focal area determination unit, configured to determine a neighborhood focal area corresponding to the current beamforming based on a target coherence factor function corresponding to the neighborhood area and a preset proportional coefficient;

[0059] A second focal zone extension range determination unit is used to determine the focal zone extension range corresponding to the current beamforming based on the real coordinates and theoretical center coordinates corresponding to the neighborhood focal area;

[0060] The coherent region function determination unit is used to perform function fitting on the target coherence factor function corresponding to the neighborhood focal region to determine the coherent region function of this beamforming.

[0061] Preferably, the phase correction table optimization device further comprises:

[0062] An original echo signal module, used for obtaining an original echo signal;

[0063] A target echo signal determination module, configured to perform flight time correction on the original echo signal using a phase correction table corresponding to the current transmission, and determine the target echo signal based on the corrected flight time;

[0064] A target composite signal determination module, used for performing channel composite and multi-beam composite on the target echo signal to determine the target composite signal;

[0065] The target ultrasonic image determination module is used to determine the target ultrasonic image based on the target composite signal.

[0066] An ultrasonic device comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned phase correction table optimization method when executing the computer program.

[0067] A computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the phase correction table optimization method is implemented.

[0068] The above-mentioned phase correction table optimization method, device, ultrasonic equipment and storage medium determine the effective focal area corresponding to this transmission based on the preset focus, and then use the coherent area function corresponding to the previous beam synthesis to optimize the effective focal area corresponding to this transmission, so as to combine the actual coherence characteristics with the effective focal area, thereby modulating a better focal area distribution result. Any focal point in the target focal area is used as a virtual array element, and the phase correction table corresponding to this transmission is determined based on the geometric relationship between the virtual array element and the imaging point and the correction point, so as to realize theoretical calculation based on the virtual array element, optimize and update the phase correction table, without the need for advance simulation and without occupying storage resources; and this real-time optimization method can solve the shortcomings of ideal assumptions in a real human body environment, and can achieve seamless connection and image quality improvement in conventional ultrasound processes. BRIEF DESCRIPTION OF THE DRAWINGS

[0069] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments of the present invention will be briefly introduced below. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative labor.

[0070] Figure 1is a first flow chart of a method for optimizing a phase correction table in one embodiment of the present invention;

[0071] Figure 2 yes Figure 1 A flow chart of step S101;

[0072] Figure 3 yes Figure 1 A flow chart of step S103;

[0073] Figure 4 is a schematic diagram of the geometric relationship between virtual array elements, imaging points and correction points in one embodiment of the present invention;

[0074] Figure 5 is a second flow chart of a method for optimizing a phase correction table in one embodiment of the present invention;

[0075] Figure 6 yes Figure 5 A flow chart of step S503;

[0076] Figure 7 yes Figure 5 A flow chart of step S504;

[0077] Figure 8 is a third flow chart of a phase correction table optimization method in one embodiment of the present invention;

[0078] Fig. 9 is a schematic diagram of a phase correction table optimization device in one embodiment of the present invention. DETAILED DESCRIPTION

[0079] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0080] The phase correction table optimization method provided in an embodiment of the present invention can be applied in ultrasonic equipment to iteratively optimize its phase correction table before ultrasonic wave transmission, so that after ultrasonic wave transmission, the optimized phase correction table can be used to perform phase correction on its echo signal.

[0081] The working process of ultrasonic equipment generally includes the following steps:

[0082] (1) The emission condition determination step is a step that all conventional ultrasound equipment will have. When the equipment is turned on, there will be a set of default ultrasound emission imaging conditions, including emission frequency, acoustic power, line density, imaging depth, scanning width, etc. These conditions are user-adjustable. These conditions must be determined before any imaging begins.

[0083] (2) Parameter calculation step: After the launch conditions are determined, the system will perform complex internal parameter calculations based on the default conditions or the launch conditions specified by the user to control the ultrasonic emission and reception of the entire system. Generally speaking, this step is relatively routine. After calculation once, these parameters will not be recalculated as long as the launch conditions are not changed.

[0084] (3) The ultrasonic emission step, i.e., controlling the ultrasonic probe to emit ultrasonic waves according to the emission conditions determined in (1), is no different from a general ultrasonic system and will not be described in detail here.

[0085] (4) The echo receiving step is to configure the echo receiving according to the transmitting conditions in (1) and perform sampling. This is no different from the general ultrasonic system and will not be described in detail here.

[0086] (5) Signal processing steps, including analog-to-digital conversion, low-noise amplification and other common signal processing steps. This step is no different from the general ultrasound system and will not be elaborated here.

[0087] (6) A beamforming step is used to receive the echo signal sent by the ultrasound probe, perform beamforming on the echo signals of one or more channels, and obtain one or more echo synthesis signals.

[0088] (7) Image processing steps, including but not limited to modulus, logarithmic compression, spatial compounding, gain, dynamic transformation, image enhancement, scanning digital transformation, etc. These steps are essentially the same as those of general ultrasound systems and will not be elaborated here.

[0089] (8) The display step is to display the processed image on a monitor. This solution is no different from a general ultrasound system and will not be described in detail here.

[0090] In the above parameter calculation step, except for the need to calculate common control parameters such as time delay, echo line position of each emission and other conventional parameters, there is no difference from the general ultrasound system, which will not be described here. The main improvement is that in the parameter calculation step, it is also necessary to iteratively optimize and update its phase correction table according to actual conditions.

[0091] In one embodiment, if Figure 1 As shown, a phase correction table optimization method is provided, which is described by taking the application of ultrasonic equipment as an example, and includes the following steps:

[0092] S101: Determine the effective focus area corresponding to this emission based on the preset focus;

[0093] S102: using the coherent region function corresponding to the last beamforming, optimizing the effective focal region corresponding to the current transmission, and determining the target focal region corresponding to the current transmission;

[0094] S103: Taking any focal point in the target focal area as a virtual array element, and determining a phase correction table corresponding to this transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point.

[0095] The preset focus is a preset focus, which can be represented by (fx, fz), that is, the horizontal coordinate of the preset focus is fx, and the vertical coordinate is fz. The effective focus area refers to the area including the preset focus.

[0096] Generally speaking, when the ultrasound device is turned on, the default or user-specified ultrasound emission imaging conditions are determined, including but not limited to emission frequency, acoustic power, line density, imaging depth, scanning width, etc., which are all user-adjustable. After the emission conditions are determined, the ultrasound device will perform complex internal parameter calculations based on the emission conditions to control the ultrasound emission and reception of the entire ultrasound device. The internal parameters determined include but are not limited to the preset focus.

[0097] As an example, in step S101, after determining the preset focus, the ultrasound device may, according to actual conditions, use different focal zone extension ranges to expand the horizontal coordinate fx and the vertical coordinate fz of the preset focus to determine the effective focal area corresponding to this emission.

[0098] Among them, the coherent area function is a function related to the coherent area, and the coherent area function is characterized by a value greater than 0 and less than 1, representing the coherence characteristics of the focal point in the area; the coherent area here is the area with the strongest coherence.

[0099] As an example, in step 102, the ultrasound device uses the coherent area function corresponding to the last beamforming in the current parameter calculation step to optimize the effective focal area corresponding to the current transmission and determine the target focal area corresponding to the current transmission. The specific representation is as follows:

[0100]

[0101] is the coherent region function, The coherence region is a function related to the lateral distance x and the longitudinal distance z. Its value is a number greater than 0 and less than 1, which represents the coherence feature. When this transmission is the first transmission, the ultrasound device cannot obtain the coherence region function of the last beam synthesis. In this case, the coherence feature defaults to all 1s, that is, .

[0102] In this example, fx and fz are the horizontal and vertical coordinates of the preset focus, respectively. and They are the horizontal expansion area and the vertical expansion area of ​​the focal area corresponding to this launch. is the focal area function, is the effective focal area corresponding to this emission. Generally speaking, the coherence is strongest at the focus. Therefore, the calculation of the coherent area can be used to determine the diffusion area range of the real focus during actual imaging (because the real focus cannot be an ideal point due to the non-ideality of the medium and the non-ideality of the imaging system, but will have an extended range). Because the most essential feature of the focus is complete coherence, the coherence distribution near the preset focus actually reflects the regional distribution of the real focus. Weighting the coherent area function with the effective focal area corresponding to this emission can be understood as combining the actual coherence characteristics with the expanded effective focal area, thereby modulating a better focal area distribution result.

[0103] The focal point is any point in the target focal area.

[0104] As an example, in step S103, during conventional phase correction calculation, due to the coherent emission of multiple array elements, the calculation before the actual ultrasonic wave emission will be very complicated and not easy to analyze. According to the time inversion theory in ultrasound, although the focus is the convergence point of the ultrasonic waves emitted by each array element in the probe, it can also be considered that the sound field at all positions is caused by the spherical wave emitted from the focus, but the positive and negative signs are different before and after the focus. Therefore, when we calculate the phase correction, we can only rely on the geometric relationship between the focus, the imaging point of the correction value to be calculated, and the correction point to calculate the sound path difference, so as to obtain the phase correction value at the imaging point. Under this consideration, any focal point in the target focal area is used as a virtual array element to transmit a single array element, and the transmission center line passing through the focal point is determined. The imaging point of the correction value to be calculated for an echo line to be imaged in this transmission is perpendicular to the transmission center line and is determined as its correction point. Finally, the geometric relationship between the focal point (i.e., the virtual array element), the imaging point, and the correction point is determined. Finally, based on this several-bit relationship, the difference between the phase of the imaging point and the phase of the correction point is dynamically determined, and it is determined as the phase difference that needs to be corrected. Based on all phase differences, the phase correction table is dynamically updated.

[0105] In this embodiment, based on the preset focus, the effective focal area corresponding to this transmission is determined, and then the coherent area function corresponding to the previous beam synthesis is used to optimize the effective focal area corresponding to this transmission, so as to combine the actual coherence characteristics with the effective focal area, thereby modulating a better focal area distribution result. Any focal point in the target focal area is used as a virtual array element, and based on the geometric relationship between the virtual array element and the imaging point and the correction point, the phase correction table corresponding to this transmission is determined to achieve theoretical calculation based on the virtual array element, optimize and update the phase correction table, without the need for advance simulation and without occupying storage resources; and this real-time optimization method can solve the shortcomings of ideal assumptions in a real human environment, and can achieve seamless connection and image quality improvement in conventional ultrasound processes.

[0106] In one embodiment, if Figure 2 As shown, step S101, i.e., determining the effective focus area corresponding to this emission based on the preset focus, includes:

[0107] S201: determining a focal area extension range corresponding to this transmission based on a focal area extension range corresponding to the last beamforming and a theoretical focal area extension range;

[0108] S202: Determine the effective focal area corresponding to this emission based on the preset focus and the focal zone extension range corresponding to this emission.

[0109] Among them, the focal area expansion range corresponding to the last beam synthesis includes the focal area lateral expansion range corresponding to the last beam synthesis The longitudinal extension range of the focal area corresponding to the last beam synthesis .

[0110] The theoretical focal zone expansion range refers to the focal zone expansion range determined by theoretical calculation, including the theoretical focal zone lateral expansion range. And the theoretical focal zone longitudinal expansion range As an example, the theoretical focal zone extension range may adopt but is not limited to the focal zone extension range determined by Bohr theory. Since the expression is derived based on idealized assumptions and the form is based on a simplified calculation method, the expression method is not unique, and therefore the theoretical focal zone extension range is not precisely described.

[0111] As an example, in step S201, the ultrasound device may expand the theoretical focal zone range (i.e. and ) is determined as the focal area extension range corresponding to this launch, that is , Accordingly, before the first transmission, the ultrasound device can determine the focal zone extension range corresponding to the last beamforming (i.e. and ), the focal area corresponding to the last beam synthesis can be used to expand the range (i.e. and ), the theoretical focal area expansion range (i.e. and ) to update, the update method is as follows:

[0112]

[0113] in, is the horizontal expansion area of ​​the focal zone corresponding to this launch, is the longitudinal expansion area of ​​the focal zone corresponding to this launch. The above formula can be used to indicate that the focal zone expansion range actually used during this launch is a smooth expansion range.

[0114] As an example, in step S202, the ultrasound device may use a focal area function to process the preset focus and the focal area extension range corresponding to the current transmission. Specifically, the focal area function is used to process the horizontal coordinate fx of the preset focus, the vertical coordinate fz of the preset focus, and the horizontal extension range of the focal area. And focal area vertical expansion range Processing is performed to determine the effective focal area corresponding to this emission, that is, its effective focal area is .

[0115] In this example, the focal area extension range corresponding to this transmission can be dynamically updated based on the focal area extension range corresponding to the previous beam synthesis and the theoretical focal area extension range, so that the preset focus can be expanded using the focal area extension range corresponding to this transmission, and the effective focal area corresponding to this transmission can be determined to ensure that the effective focal area is more accurate and effective.

[0116] In one embodiment, if Figure 3 As shown, step S103, i.e., determining a phase correction table corresponding to this transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point, includes:

[0117] S301: Determine a first distance between a virtual array element and an imaging point;

[0118] S302: Determine a second distance between the virtual array element and the calibration point;

[0119] S303: determining an original sound path difference corresponding to the virtual array element based on the first distance and the second distance corresponding to the virtual array element;

[0120] S304: Integrate the original sound path difference corresponding to the virtual array element in the target focus area to determine the target sound path difference corresponding to the virtual array element;

[0121] S305: Determine a phase correction table corresponding to this transmission based on the target sound path differences corresponding to all virtual array elements.

[0122] The first distance is the geometric distance between any focal point (virtual array element) and the imaging point, which can be expressed as The second distance is the distance between any focal point (virtual array element) and the calibration point, which can be expressed as express.

[0123] As an example, in step S301, after determining the target focal area corresponding to the current transmission, the ultrasound device may determine any focal point of the target focal area as a virtual array element, and determine the geometric distance between the virtual array element and any imaging point of the imaging system (i.e., a sampling point on an echo line to be imaged in the current transmission) as the first distance .

[0124] As an example, in step S302, after determining the target focal area corresponding to the current transmission, the ultrasound device may determine any focal point of the target focal area as a virtual array element, transmit with the virtual array element, determine its transmission center line, determine any imaging point of the imaging system perpendicular to its transmission center line as a calibration point, and determine the geometric distance between any virtual array element and the calibration point as the second distance .

[0125] As an example, in step S303, the ultrasound device determines the first distance corresponding to the virtual array element. and the second distance After that, the first distance and the second distance The difference is determined as the original sound path difference ,Right now .

[0126] As an example, in step S304, the ultrasonic device determines the original sound path difference After that, the original sound path difference d can be processed by two-dimensional integration in the target focal area to determine the target sound path difference , that is, the target sound path difference is the result of integrating the original sound path difference d, which is characterized as a function related to the position of the imaging point.

[0127] As an example, in step S305, the ultrasonic device determines the target sound path difference Then, based on the target sound path difference , determine the phase correction table corresponding to this transmission. It is characterized as a function related to the position of the imaging point, and there is a corresponding relationship between the sound path difference and the phase difference , is the target sound path difference, is the wavelength, is the phase difference. Based on the above correspondence, the phase difference that needs to be corrected for each imaging point can be determined. Based on the phase differences that need to be corrected for all imaging points, the phase correction table corresponding to this transmission is determined.

[0128] like Figure 4 The figure shows the relationship between any focal point A (i.e., virtual array element), imaging point B, and correction point C in the target focal area during conventional non-deflected transmission under linear array. The figure assumes the conventional non-deflected transmission of the linear array probe, any virtual array element A is assumed to be the focus of the transmission, CA is the transmission center line of the transmission, B is a sampling point on an echo line to be imaged for the transmission, i.e., the imaging point, and BC is perpendicular to AC. Point C is the correction point, which means that the phase of point B is corrected to the phase of point C, and the specific correction needs to be based on the calculated phase difference.

[0129] In this embodiment, any focal point in the target focal area is determined as a virtual array element, and its original acoustic path difference is determined based on the geometric relationship between the virtual array element and the imaging point and the correction point; the original acoustic path difference is then integrated to determine the target acoustic path difference; according to the corresponding relationship between the acoustic path difference and the phase difference, the corresponding phase correction table is dynamically determined. This process performs theoretical calculations based on the virtual array element to optimize and update the phase correction table without the need for advance simulation and without occupying storage resources; and this real-time optimization method can solve the shortcomings of ideal assumptions in a real human body environment, and can achieve seamless connection and image quality improvement in conventional ultrasound processes.

[0130] In one embodiment, if Figure 5 As shown, after step S103, that is, after determining the phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point, the phase correction table optimization method further includes:

[0131] S501: Determine a neighboring area corresponding to a preset focus;

[0132] S502: performing flight time sampling based on the original echo signal corresponding to the adjacent area to determine the full-channel echo signal corresponding to the adjacent area;

[0133] S503: performing coherence analysis on the full-channel echo signals corresponding to the adjacent area to determine a target coherence factor function corresponding to the adjacent area;

[0134] S504: Determine a focal zone extension range and a coherent region function corresponding to the current beamforming based on a target coherence factor function corresponding to the adjacent region.

[0135] Among them, the adjacent area is based on the preset focus The original echo signal refers to the unprocessed echo signal.

[0136] As an example, in step S501, the ultrasound device uses the preset focus of this transmission. A neighboring area (ROI area) is determined as the center. The neighboring area can be a matrix or a circle. The area size can be within the range of 10 wavelengths. The specific value depends on the experiment.

[0137] As an example, in step S502, after determining the neighboring area corresponding to the preset focus, the ultrasonic device can obtain the original echo signal corresponding to each neighborhood point in the neighboring area, and perform flight time (i.e., time of flight, hereinafter referred to as TOF) on all original echo signals, specifically including the time taken for the ultrasonic wave to be transmitted to the imaging position plus the time taken for the echo reflected from the imaging position to each array element of the probe, i.e., transmission time + reception time. After determining the flight time corresponding to each neighborhood point, the ultrasonic device can interpolate the original echo signal corresponding to it according to the flight time to determine the full-channel echo signal corresponding to the neighboring area. The neighborhood point here is any point in the neighboring area. The full-channel echo signal here can be understood as the echo signal formed after the original echo signal is interpolated.

[0138] In this example, after determining the adjacent area corresponding to the preset focus, the ultrasonic device can calculate the TOF corresponding to the adjacent area. The information of each point of the TOF is not for imaging, but simply for focal area analysis. It can be specifically understood as performing TOF sampling on specific sampling points near the adjacent area of ​​the preset focus to determine the TOF of each specific sampling point. The specific sampling points here can be understood as sampling points for which TOF needs to be calculated to determine the focal area. The sampling method can be determined based on system resources. For example, one-quarter of the wavelength can be used as a sampling distance for equal-interval analysis. If the system resources are abundant, it can be more intensive. In theory, the denser the better.

[0139] As an example, in step S503, the ultrasonic device performs a coherence analysis on the full-channel echo signal corresponding to the adjacent area. Specifically, coherence factors such as but not limited to CF coherence factor, SLSC coherence factor or DORT coherence factor may be used to evaluate each neighborhood point in the adjacent area and determine the coherence factor corresponding to each neighborhood point. The coherence factor here can be any coherence factor as long as it can evaluate the relative ratio of the coherence degree of the full-channel echo signal of the neighborhood point. Then, based on the coherence factors corresponding to all neighborhood points in the adjacent area, the target coherence factor function corresponding to the adjacent area is determined. The target coherence factor function here can be understood as the coherence factor function finally determined by performing a coherence analysis on the full-channel echo signal.

[0140] As an example, in step S504, the ultrasound device may radially intercept the target coherence factor function corresponding to each neighborhood point in the neighborhood area, and specifically may use the coherence factor corresponding to the ideal focal position as an evaluation criterion, and determine the area that matches the coherence factor corresponding to the ideal focal position as the target focal area corresponding to this emission. After determining the target focal area corresponding to this emission, the ultrasound device may determine the focal area expansion range corresponding to this beamforming based on the coordinate information of the target focal area, that is, determine the focal area lateral expansion range corresponding to this beamforming. The longitudinal extension range of the focal area corresponding to this beamforming , so that the subsequent iterations can update the horizontal expansion area of ​​the focal zone corresponding to the next launch Vertical extension area of ​​the focal zone In addition, the ultrasound device may also perform function fitting on the target coherence factor function corresponding to the target focal area to determine the coherence area function corresponding to this emission so as to iteratively update the target focal area corresponding to the next emission.

[0141] In this example, before each ultrasonic transmission, the ultrasonic device needs to iteratively update the phase correction table corresponding to the current transmission based on the focal zone extension range and coherent area function corresponding to the previous beamforming; after the phase correction table is updated, in the beamforming step, the phase correction table is used to perform TOF correction on the original echo signal, and it is necessary to further obtain the adjacent area corresponding to the preset focus, perform flight time sampling based on the original echo signal corresponding to the adjacent area to determine the full-channel echo signal, perform coherence analysis on the full-channel echo signal, determine the target coherence factor function corresponding to the adjacent area, and then determine the focal zone extension range and coherent area function corresponding to this beamforming based on the target coherence factor function, so as to iteratively update the corresponding phase correction table before the next ultrasonic transmission.

[0142] In one embodiment, if Figure 6 As shown, step S503, i.e., performing coherence analysis on the full-channel echo signal corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area, includes:

[0143] S601: Calculate the coherence factor of the full-channel echo signal corresponding to the adjacent area to determine the original coherence factor function corresponding to the adjacent area;

[0144] S602: Perform smoothing filtering on the original coherence factor function corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area.

[0145] As an example, in step S601, the ultrasound device may directly calculate the corresponding CF coherence factor or SLSC coherence factor according to the full channel echo signal corresponding to each neighborhood point in the neighborhood area. For example, , where is the echo signal of the ith channel, that is, the value obtained by summing and then finding the modulus divided by the value obtained by summing and finding the modulus is its corresponding CF coherence factor. Alternatively, the ultrasonic device may perform a short-time Fourier transform on the full-channel echo signal corresponding to any neighborhood point in the adjacent area and then calculate its coherence to determine its DORT coherence factor. Alternatively, the ultrasonic device may perform coherence calculation on the full-channel echo signal corresponding to any neighborhood point in the adjacent area with the full-channel echo signals corresponding to other neighborhood points to determine the corresponding local coherence factor. In this example, after calculating the coherence factors such as the CF coherence factor, SLSC coherence factor, DORT coherence factor or local coherence factor corresponding to each neighborhood point in the adjacent area, the ultrasonic device may construct the original coherence factor function corresponding to the adjacent area based on the coherence factors corresponding to all neighborhood points.

[0146] As an example, in step S602, after obtaining the original coherence factor function corresponding to the adjacent area, the ultrasonic device may perform two-dimensional smoothing filtering on the original coherence factor function to determine the target coherence factor function, so as to reduce the interference of noise or clutter. The filter function used to implement the two-dimensional smoothing filter is not limited here, and conventional low-pass filtering can be implemented. In this solution, Gaussian filtering is generally used. In addition, mean filtering is also a commonly used optional solution.

[0147] In one embodiment, if Figure 7 As shown, step S504, i.e., determining the focal area extension range and the coherent area function corresponding to the current beamforming based on the target coherence factor function corresponding to the adjacent area, includes:

[0148] S701: Determine a neighborhood focal area corresponding to the current beamforming based on a target coherence factor function corresponding to the neighboring area and a preset proportional coefficient;

[0149] S702: Determine the focal area extension range corresponding to this beamforming based on the real coordinates and theoretical center coordinates corresponding to the neighborhood focal area;

[0150] S703: Perform function fitting on the target coherence factor function corresponding to the neighborhood focal area to determine the coherence area function for this beamforming.

[0151] The preset proportionality coefficient is a preset proportionality coefficient that is set in advance, and can be specifically understood as a coherence factor corresponding to an ideal focal position.

[0152] As an example, in step S701, the ultrasonic device may radially intercept the target coherence factor function corresponding to the neighboring area. Specifically, the coherence factor corresponding to each neighborhood point in the neighboring area may be compared with a preset proportional coefficient, and the area formed by all neighborhood points less than the preset proportional coefficient σ is identified as a non-focal area; conversely, the area formed by all neighborhood points not less than the preset proportional coefficient σ is identified as a neighborhood focal area.

[0153] As an example, in step S702, after determining the neighborhood focal area corresponding to the current beamforming, the ultrasound device may determine the focal zone extension range corresponding to the current beamforming based on the real coordinates and theoretical center coordinates of the neighborhood focal area, that is, determine the focal zone lateral extension range corresponding to the current beamforming. The longitudinal extension range of the focal area corresponding to this beamforming , so that the subsequent iterations can update the horizontal expansion area of ​​the focal zone corresponding to the next launch Vertical extension area of ​​the focal zone .

[0154] In this example, the ultrasound device can determine the minimum value of the horizontal coordinate, the maximum value of the horizontal coordinate, the minimum value of the vertical coordinate, and the maximum value of the vertical coordinate corresponding to the neighborhood focal area according to the real coordinates of the neighborhood focal area; the horizontal expansion range of the focal area corresponding to this beam synthesis can be determined by subtracting the horizontal coordinate of the theoretical center from the minimum value of the horizontal coordinate and the maximum value of the horizontal coordinate. ; Subtract the ordinate of the theoretical center from the minimum and maximum ordinate values ​​to determine the longitudinal extension range of the focal area corresponding to this beamforming. , so that the subsequent iterations can update the horizontal expansion area of ​​the focal zone corresponding to the next launch Vertical extension area of ​​the focal zone .

[0155] As an example, in step S703, the ultrasound device can perform function fitting on the target coherence factor function corresponding to the neighborhood focal area based on the neighborhood focal area to determine the coherence area function corresponding to this beamforming. In this example, the ultrasound device can use but is not limited to the Gaussian function to perform function fitting on the target coherence factor function corresponding to the focal area to obtain the Gaussian function parameters, thereby obtaining a coherence area function required to be transferred in the previous parameter calculation. , since it is a Gaussian function, the two-dimensional integral can be calculated more conveniently in the focal zone calculation step.

[0156] In one embodiment, if Figure 8 As shown, after step S103, that is, after determining the phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point, the phase correction table optimization method further includes:

[0157] S801: Acquire original echo signal;

[0158] S802: using the phase correction table corresponding to this transmission, performing flight time correction on the original echo signal, and determining the target echo signal based on the corrected flight time;

[0159] S803: performing channel compounding and multi-beam compounding on the target echo signal to determine a target composite signal;

[0160] S804: Determine a target ultrasound image based on the target composite signal.

[0161] As an example, in step S801, the ultrasound device may acquire an original echo signal, where the original echo signal refers to an unprocessed echo signal.

[0162] As an example, in step S802, the ultrasound device may use the phase correction table corresponding to this transmission to perform flight time correction (i.e., TOF correction) on the original echo signal, and determine the target echo signal based on the corrected flight time. In this example, the ultrasound device needs to perform flight time correction for each imaging sampling point corresponding to the original echo signal, where the imaging sampling point is the sampling point that needs to be imaged and displayed. In this example, since the essence of phase correction is the correction of the transmission time, it will not affect the receiving time, and the flight time = transmission time + receiving time, therefore, the target sound path difference can be added to the conventionally calculated transmission distance. , in order to correct its launch time, which can be determined by the following formula:

[0163]

[0164] is the launch distance in conventional calculations, is the receiving distance in conventional ultrasonic calculation, c is the speed of sound, and the distance divided by the speed is the transmission time and the receiving time respectively. The sum of the two is the total flight time. When we have the phase correction table, we need to make corrections on the transmission distance, that is, use the transmission distance in conventional calculation Minus the correction distance , and then divided by the speed of sound c to get the corrected launch time, and then determine the corrected flight time Then, based on the corrected flight time, the original echo signal is interpolated to determine the target echo signal.

[0165] As an example, in step S803, after determining the target echo signal, the ultrasound device may perform channel compounding on the target echo signal, specifically weighted compounding of the probe multi-element echo signals belonging to the same sampling point reflection signal. This step is also a standard operation in ultrasound beam synthesis and will not be described in detail. Then, the echo signal after channel compounding is multi-beam compounded, specifically weighted compounding of the echo line signals at the same position obtained during multiple transmissions to determine the target composite signal. This step is also a standard operation in ultrasound imaging and will not be described in detail.

[0166] As an example, in step S804, after determining the target composite signal, the ultrasound device may perform image processing on the target composite signal, including but not limited to modulus finding, logarithmic compression, spatial compounding, gain, dynamic transformation, image enhancement, scanning digital conversion, etc. These steps are not the focus of this solution and are not essentially different from general ultrasound systems, so they will not be elaborated here.

[0167] It should be understood that the order of execution of the steps in the above embodiment does not necessarily mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiment of the present invention.

[0168] In one embodiment, a phase correction table optimization device is provided, and the phase correction table optimization device corresponds one-to-one to the phase correction table optimization method in the above embodiment. Fig. 9 As shown, the phase correction table optimization device includes an effective focus area determination module 901, a target focus area determination module 902 and a phase correction table determination module 903. The detailed description of each functional module is as follows:

[0169] An effective focus area determination module 901 is used to determine the effective focus area corresponding to this emission based on a preset focus;

[0170] The target focal region determination module 902 is used to optimize the effective focal region corresponding to the current transmission by using the coherent region function corresponding to the previous beamforming, and determine the target focal region corresponding to the current transmission;

[0171] The phase correction table determination module 903 is used to take any focal point in the target focal area as a virtual array element, and determine the phase correction table corresponding to this transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point.

[0172] Preferably, the effective focus area determination module 901 includes:

[0173] A first zone extension range determination unit, configured to determine a focal zone extension range corresponding to this transmission based on a focal zone extension range corresponding to the previous beamforming and a theoretical focal zone extension range;

[0174] The effective focus area determination unit is used to determine the effective focus area corresponding to this emission based on the preset focus and the focal area extension range corresponding to this emission.

[0175] Preferably, the phase correction table determination module 903 includes:

[0176] A first distance determining unit, used to determine a first distance between a virtual array element and an imaging point;

[0177] A second distance determining unit, used to determine a second distance between the virtual array element and the calibration point;

[0178] an original sound path difference determining unit, configured to determine an original sound path difference corresponding to the virtual array element based on a first distance and a second distance corresponding to the virtual array element;

[0179] A target sound path difference determination unit is used to perform integration processing on the original sound path difference corresponding to the virtual array element in the target focus area to determine the target sound path difference corresponding to the virtual array element;

[0180] The phase correction table determination unit is used to determine the phase correction table corresponding to the current transmission based on the target sound path differences corresponding to all virtual array elements.

[0181] Preferably, the phase correction table optimization device further comprises:

[0182] A neighboring area determination module, used to determine a neighboring area corresponding to a preset focus;

[0183] A full-channel echo signal determination module is used to perform flight time sampling based on the original echo signal corresponding to the adjacent area to determine the full-channel echo signal corresponding to the adjacent area;

[0184] A target coherence factor function determination module is used to perform coherence analysis on the full-channel echo signal corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area;

[0185] The current beamforming information determination module is used to determine the focal area extension range and coherent area function corresponding to the current beamforming based on the target coherence factor function corresponding to the adjacent area.

[0186] Preferably, the target coherence factor function determination module comprises:

[0187] The original coherence factor function determination unit is used to calculate the coherence factor of the full-channel echo signal corresponding to the adjacent area and determine the original coherence factor function corresponding to the adjacent area;

[0188] The target coherence factor function determination unit is used to perform smoothing filtering on the original coherence factor function corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area.

[0189] Preferably, based on the target coherence factor function corresponding to the adjacent area, determining the focal area extension range and the coherence area function corresponding to the current beamforming includes:

[0190] A neighborhood focal region determining unit, configured to determine a neighborhood focal region corresponding to the current beamforming based on a target coherence factor function corresponding to the neighborhood region and a preset proportional coefficient;

[0191] A second focal zone extension range determination unit is used to determine the focal zone extension range corresponding to this beamforming based on the real coordinates and theoretical center coordinates corresponding to the neighborhood focal area;

[0192] The coherent region function determination unit is used to perform function fitting on the target coherence factor function corresponding to the neighborhood focal region to determine the coherent region function of this beamforming.

[0193] Preferably, the phase correction table optimization device further comprises:

[0194] An original echo signal module, used for obtaining an original echo signal;

[0195] A target echo signal determination module is used to perform flight time correction on the original echo signal using the phase correction table corresponding to this transmission, and determine the target echo signal based on the corrected flight time;

[0196] A target composite signal determination module is used to perform channel composite and multi-beam composite on the target echo signal to determine the target composite signal;

[0197] The target ultrasonic image determination module is used to determine the target ultrasonic image based on the target composite signal.

[0198] For the specific definition of the phase correction table optimization device, please refer to the definition of the phase correction table optimization method above, which will not be repeated here. Each module in the above-mentioned phase correction table optimization device can be implemented in whole or in part by software, hardware and a combination thereof. The above-mentioned modules can be embedded in or independent of the processor in the ultrasonic device in the form of hardware, or can be stored in the memory of the ultrasonic device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

[0199] In one embodiment, an ultrasound device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method for optimizing the phase correction table in the above embodiment is implemented, for example Figure 1 S101-S103 as shown, or Figure 2 , Figure 3 , Figures 5 to 8Alternatively, when the processor executes the computer program, the functions of each module / unit in the embodiment of the phase correction table optimization device are realized, for example Fig. 9 The functions of the phase correction table determination module 903 are not described here in detail to avoid repetition.

[0200] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the method for optimizing the phase correction table in the above embodiment is implemented. For example, Figure 1 S101-S103 as shown, or Figure 2 , Figure 3 , Figures 5 to 8 Alternatively, when the computer program is executed by a processor, the functions of each module / unit in the above-mentioned phase correction table optimization device are realized, for example Fig. 9 The functions of the phase correction table determination module 903 are not described here in detail to avoid repetition.

[0201] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0202] Those skilled in the art can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example. In actual applications, the above-mentioned functions can be distributed and completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above.

[0203] The embodiments described above are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features may be replaced by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included in the protection scope of the present invention.

Claims

1. A phase correction table optimization method, characterized in that: include: Based on the preset focus, determine the effective focus area corresponding to this launch; The coherent area function corresponding to the last beam synthesis is used to optimize the effective focal area corresponding to this transmission, and determine the target focal area corresponding to this transmission. The coherent area function is a function related to the coherent area. The coherent area function is characterized by a value greater than 0 and less than 1, representing the coherence characteristics of the focal point in the area. The coherent area is the area with the strongest coherence. Taking any focal point in the target focal area as a virtual array element, and determining a phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point; The optimization process of the effective focal area corresponding to the current transmission using the coherent area function corresponding to the previous beamforming is to weight the coherent area function and the effective focal area corresponding to the current transmission; The determining of the phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point includes: Determining a first distance between the virtual array element and an imaging point; Determining a second distance between the virtual array element and the calibration point; Determine an original sound path difference corresponding to the virtual array element based on a first distance and a second distance corresponding to the virtual array element; In the target focus area, integrating the original sound path difference corresponding to the virtual array element to determine the target sound path difference corresponding to the virtual array element; Based on the target sound path differences corresponding to all the virtual array elements, a phase correction table corresponding to this transmission is determined.

2. The phase correction table optimization method according to claim 1, characterized in that: The determining of the effective focus area corresponding to the current emission based on the preset focus includes: Based on the focal area extension range corresponding to the previous beamforming and the theoretical focal area extension range, the focal area extension range corresponding to the current launch is determined; Based on the preset focus and the focal zone extension range corresponding to this emission, the effective focal area corresponding to this emission is determined.

3. The phase correction table optimization method according to claim 1, characterized in that: After determining the phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point, the phase correction table optimization method further includes: Determine a neighboring area corresponding to a preset focus; Performing time-of-flight sampling based on the original echo signal corresponding to the adjacent area to determine the full-channel echo signal corresponding to the adjacent area; Performing coherence analysis on the full-channel echo signal corresponding to the adjacent area to determine a target coherence factor function corresponding to the adjacent area; Based on the target coherence factor function corresponding to the adjacent area, the focal area extension range and the coherence area function corresponding to the current beamforming are determined.

4. The phase correction table optimization method according to claim 3, characterized in that: The performing coherence analysis on the full-channel echo signal corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area includes: Calculating the coherence factor of the full-channel echo signal corresponding to the adjacent area to determine the original coherence factor function corresponding to the adjacent area; The original coherence factor function corresponding to the adjacent area is smoothed and filtered to determine the target coherence factor function corresponding to the adjacent area.

5. The phase correction table optimization method according to claim 3, characterized in that: The determining, based on the target coherence factor function corresponding to the adjacent area, a focal area extension range and a coherent area function corresponding to the current beamforming includes: Determining a neighborhood focal area corresponding to the current beamforming based on a target coherence factor function corresponding to the neighborhood area and a preset proportional coefficient; Determine the focal zone extension range corresponding to this beamforming based on the real coordinates and theoretical center coordinates corresponding to the neighborhood focal area; Function fitting is performed on the target coherence factor function corresponding to the neighborhood focal area to determine the coherence area function of this beamforming.

6. The phase correction table optimization method according to claim 1, characterized in that: After determining the phase correction table corresponding to the current transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point, the phase correction table optimization method further includes: Obtaining the original echo signal; Using the phase correction table corresponding to this transmission, the original echo signal is corrected for the flight time, and the target echo signal is determined based on the corrected flight time; Performing channel compounding and multi-beam compounding on the target echo signal to determine a target composite signal; A target ultrasound image is determined based on the target composite signal.

7. A phase correction table optimization device, characterized in that: include: An effective focus area determination module, used to determine the effective focus area corresponding to this emission based on a preset focus; A target focal area determination module is used to optimize the effective focal area corresponding to this transmission using the coherent area function corresponding to the last beam synthesis, and determine the target focal area corresponding to this transmission. The coherent area function is a function related to the coherent area. The coherent area function is characterized by a value greater than 0 and less than 1, representing the coherence characteristics of the focal point in the area. The coherent area is the area with the strongest coherence. A phase correction table determination module, used to take any focal point in the target focal area as a virtual array element, and determine the phase correction table corresponding to this transmission based on the geometric relationship between the virtual array element, the imaging point and the correction point; The target focal area determination module is used to weight the coherent area function and the effective focal area corresponding to this transmission; The phase correction table determination module comprises: A first distance determining unit, used to determine a first distance between the virtual array element and the imaging point; A second distance determining unit, used to determine a second distance between the virtual array element and the calibration point; an original sound path difference determining unit, configured to determine the original sound path difference corresponding to the virtual array element based on the first distance and the second distance corresponding to the virtual array element; A target sound path difference determination unit, configured to perform integration processing on the original sound path difference corresponding to the virtual array element in the target focus area to determine the target sound path difference corresponding to the virtual array element; The phase correction table determining unit is used to determine the phase correction table corresponding to the current transmission based on the target sound path differences corresponding to all the virtual array elements.

8. The phase correction table optimization device according to claim 7, characterized in that: The effective focus area determination module comprises: A first zone extension range determination unit, configured to determine a focal zone extension range corresponding to this transmission based on a focal zone extension range corresponding to the previous beamforming and a theoretical focal zone extension range; The effective focus area determination unit is used to determine the effective focus area corresponding to this emission based on the preset focus and the focal area extension range corresponding to this emission.

9. The phase correction table optimization device according to claim 7, characterized in that: The phase correction table optimization device also includes: A neighboring area determination module, used to determine a neighboring area corresponding to a preset focus; A full-channel echo signal determination module, configured to perform time-of-flight sampling based on the original echo signal corresponding to the adjacent area to determine the full-channel echo signal corresponding to the adjacent area; A target coherence factor function determination module is used to perform coherence analysis on the full-channel echo signal corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area; The current beamforming information determination module is used to determine the focal area extension range and the coherent area function corresponding to the current beamforming based on the target coherence factor function corresponding to the adjacent area.

10. The phase correction table optimization device according to claim 9, characterized in that: The target coherence factor function determination module comprises: An original coherence factor function determining unit, configured to calculate the coherence factor of the full-channel echo signal corresponding to the adjacent area, and determine the original coherence factor function corresponding to the adjacent area; The target coherence factor function determining unit is used to perform smoothing filtering on the original coherence factor function corresponding to the adjacent area to determine the target coherence factor function corresponding to the adjacent area.

11. The phase correction table optimization device according to claim 9, characterized in that: The determining, based on the target coherence factor function corresponding to the adjacent area, a focal area extension range and a coherent area function corresponding to the current beamforming includes: A neighborhood focus area determination unit, configured to determine a neighborhood focus area corresponding to the current beamforming based on a target coherence factor function corresponding to the neighborhood area and a preset proportional coefficient; A second focal zone extension range determination unit is used to determine the focal zone extension range corresponding to the current beamforming based on the real coordinates and theoretical center coordinates corresponding to the neighborhood focal area; The coherent region function determination unit is used to perform function fitting on the target coherence factor function corresponding to the neighborhood focal region to determine the coherent region function of this beamforming.

12. The phase correction table optimization device according to claim 7, characterized in that: The phase correction table optimization device also includes: An original echo signal module, used for obtaining an original echo signal; A target echo signal determination module, used to perform flight time correction on the original echo signal using a phase correction table corresponding to this transmission, and determine the target echo signal based on the corrected flight time; A target composite signal determination module, used for performing channel composite and multi-beam composite on the target echo signal to determine the target composite signal; The target ultrasonic image determination module is used to determine the target ultrasonic image based on the target composite signal.

13. An ultrasound device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the phase correction table optimization method according to any one of claims 1 to 6 is implemented.

14. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the phase correction table optimization method according to any one of claims 1 to 6 is implemented.

Images