Ultrasonic beam imaging method, ultrasonic equipment and device

Abstract

The present application discloses an imaging method, ultrasonic equipment and device for ultrasonic beam images. The present application obtains ultrasonic array echo data received by ultrasonic array elements; then uses the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as basic curves to perform interpolation processing on the ultrasonic array echo data to obtain three reference array apodization data; then uses the three reference array apodization data to perform apodization processing on the ultrasonic array echo data to obtain three reference array radio frequency identification RF data; performs cross-correlation processing on the three reference array RF data to obtain target array RF data, and performs beam synthesis and imaging based on the target array RF data to obtain a target ultrasonic beam image. Therefore, the present application improves the contrast resolution of the RF image after ultrasonic beam synthesis without affecting the main lobe, and improves the contrast resolution of the image within the entire image range.

Description

Technical Field

[0001] The present application relates to the technical field of ultrasonic imaging, and in particular to an imaging method, ultrasonic equipment and device for ultrasonic beam images. Background Art

[0002] Medical ultrasound examination is widely used to examine diseases in a series of parts such as the carotid artery, abdomen, and heart. These examinations rely on the quality of ultrasound images, and the quality of ultrasound images relies on the resolution of ultrasound images. The resolution of ultrasound images is mainly divided into spatial resolution and contrast resolution. However, the RF (Radio Frequency) image synthesized by the ultrasound beam is greatly affected by the side lobes, resulting in artifacts in the image, making the contrast resolution low. Therefore, reducing the influence of side lobes is one of the effective ways to improve the contrast resolution.

[0003] The existing ultrasound image sidelobe suppression is mainly achieved through the apodization window function, and the more commonly used ones are Hanning window, Hamming window, etc. However, this kind of method will increase the width of the main lobe while reducing the side lobe, resulting in a decrease in the lateral spatial resolution of the image. Therefore, how to improve the contrast resolution of the RF image after ultrasound beam synthesis without affecting the main lobe is an urgent problem to be solved. Summary of the invention

[0004] The purpose of the present application is to provide an imaging method, ultrasound equipment and device for ultrasonic beam images, so as to solve the problem of how to improve the contrast resolution of the RF image after ultrasonic beam synthesis without affecting the main lobe.

[0005] In a first aspect, the present application provides an imaging method of an ultrasonic beam image, the method comprising:

[0006] Acquiring ultrasonic array echo data received by ultrasonic array elements;

[0007] Using a Hanning window function curve, a first alternating apodization curve, and a second alternating apodization curve as basic curves, respectively, interpolation processing is performed on the ultrasonic array echo data to obtain three reference array apodization data; wherein the activated ultrasonic array element represented by the first alternating apodization curve is the inactivated ultrasonic array element represented by the second alternating apodization curve, and the inactivated ultrasonic array element represented by the first alternating apodization curve is the activated ultrasonic array element represented by the second alternating apodization curve;

[0008] Using the three reference array apodization data respectively, the ultrasound array echo data is apodized to obtain three reference array radio frequency identification (RF) data;

[0009] The three reference array RF data are cross-correlated to obtain target array RF data, and beam synthesis and imaging are performed based on the target array RF data to obtain a target ultrasonic beam image.

[0010] In a possible implementation, the data volume of the ultrasound array echo data is the product of the number of system channels in the ultrasound system and the number of sampling depths in the depth direction;

[0011] The method uses the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as basic curves to perform interpolation processing on the ultrasound array echo data to obtain three reference array apodization data, including:

[0012] Obtaining the number of transmit apertures corresponding to each sampling depth in the depth direction;

[0013] obtaining first reference array apodization data using the first alternating apodization curve as a base curve based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the first alternating apodization curve; and

[0014] obtaining second reference array apodization data using the second alternating apodization curve as a base curve based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the second alternating apodization curve; and

[0015] Based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the Hanning window function, third reference array apodization data with the Hanning window function curve as a basic curve is obtained.

[0016] In a possible implementation, the using the three reference array apodization data respectively to perform apodization processing on the ultrasound array echo data to obtain three reference array radio frequency identification RF data includes:

[0017] Determining position information of a transmission line corresponding to the ultrasound array echo data;

[0018] Based on the position information of the transmission line, the three reference array apodization data are respectively translated to obtain three target array apodization data;

[0019] The three reference array RF data are obtained based on the three target array apodization data and the ultrasound array echo data.

[0020] In a possible implementation manner, performing cross-correlation processing on the three reference array RF data to obtain target array RF data includes:

[0021] performing cross-correlation processing on first reference array RF data using the first alternating apodization curve as a basic curve and second reference array RF data using the second alternating apodization curve as a basic curve to obtain fourth reference array RF data;

[0022] The third reference array RF data using the Hanning window function curve as a basic curve and the fourth reference array RF data are cross-correlated to obtain the target array RF data.

[0023] In a possible implementation manner, performing cross-correlation processing on first reference array RF data using the first alternating apodization curve as a base curve and second reference array RF data using the second alternating apodization curve as a base curve to obtain fourth reference array RF data includes:

[0024] Adding the RF data at the same position in the first reference array RF data and the second reference array RF data to obtain fifth reference array RF data, and determining a first correlation array coefficient based on the first reference array RF data and the second reference array RF data;

[0025] The fourth reference array RF data is obtained based on the fifth reference array RF data and the first correlation array coefficient.

[0026] In a possible implementation manner, determining a first correlation array coefficient based on the first reference array RF data and the second reference array RF data includes:

[0027] For each sampling depth in the depth direction: determining a reference coefficient factor corresponding to the sampling depth based on a preset number of RF data centered on the sampling depth in the first reference array RF data and a preset number of RF data centered on the sampling depth in the second reference array RF data; determining a target coefficient factor corresponding to the sampling depth based on the reference coefficient factor corresponding to the sampling depth and a preset threshold;

[0028] The first correlation array coefficient is determined based on a target coefficient factor corresponding to each sampling depth in the depth direction.

[0029] In a possible implementation manner, performing cross-correlation processing on the third reference array RF data using the Hanning window function curve as a base curve and the fourth reference array RF data to obtain the target array RF data includes:

[0030] performing smaller value processing on the RF data at the same position in the third reference array RF data and the fourth reference array RF data to obtain sixth reference array RF data, and determining a second correlation array coefficient based on the third reference array RF data and the fourth reference array RF data;

[0031] The target array RF data is obtained based on the sixth reference array RF data and the second correlation array coefficient.

[0032] In a possible implementation manner, determining the second correlation array coefficient based on the third reference array RF data and the fourth reference array RF data includes:

[0033] For each sampling depth in the depth direction: based on a preset number of RF data centered at the sampling depth in the third reference array RF data and the preset number of RF data centered at the sampling depth in the fourth reference array RF data, a reference coefficient factor corresponding to the sampling depth is obtained; based on the reference coefficient factor corresponding to the sampling depth and a preset threshold, a target coefficient factor corresponding to the sampling depth is determined;

[0034] The second correlation array coefficient is determined based on the target coefficient factor corresponding to each sampling depth in the depth direction.

[0035] In a second aspect, an embodiment of the present application provides an ultrasound device, the ultrasound device comprising at least one processor and at least one memory; wherein the memory stores a program code, and when the program code is executed by the processor, the processor performs the following process:

[0036] Acquiring ultrasonic array echo data received by ultrasonic array elements;

[0037] Using a Hanning window function curve, a first alternating apodization curve, and a second alternating apodization curve as basic curves, respectively, interpolation processing is performed on the ultrasonic array echo data to obtain three reference array apodization data; wherein the activated ultrasonic array element represented by the first alternating apodization curve is the inactivated ultrasonic array element represented by the second alternating apodization curve, and the inactivated ultrasonic array element represented by the first alternating apodization curve is the activated ultrasonic array element represented by the second alternating apodization curve;

[0038] Using the three reference array apodization data respectively, the ultrasound array echo data is apodized to obtain three reference array radio frequency identification (RF) data;

[0039] The three reference array RF data are cross-correlated to obtain target array RF data, and beam synthesis and imaging are performed based on the target array RF data to obtain a target ultrasonic beam image.

[0040] In a third aspect, an embodiment of the present application provides an imaging device for ultrasonic beam images, the device comprising:

[0041] An acquisition module, used for acquiring ultrasound array echo data received by ultrasound array elements;

[0042] an interpolation module, configured to perform interpolation processing on the ultrasonic array echo data using a Hanning window function curve, a first alternating apodization curve, and a second alternating apodization curve as basic curves, respectively, to obtain three reference array apodization data; wherein the activated ultrasonic array element represented by the first alternating apodization curve is the inactivated ultrasonic array element represented by the second alternating apodization curve, and the inactivated ultrasonic array element represented by the first alternating apodization curve is the activated ultrasonic array element represented by the second alternating apodization curve;

[0043] an apodization module, configured to perform apodization processing on the ultrasound array echo data using the three reference array apodization data respectively, to obtain three reference array radio frequency identification RF data;

[0044] The synthesis module is used to perform cross-correlation processing on the three reference array RF data to obtain target array RF data, and perform beam synthesis and imaging based on the target array RF data to obtain a target ultrasonic beam image.

[0045] In a fourth aspect, the present application provides a computer-readable storage medium. When the instructions in the computer-readable storage medium are executed by an ultrasound device, the ultrasound device is enabled to perform the imaging method of ultrasound beam images as described in any one of the first aspects above.

[0046] In a fifth aspect, the present application provides a computer program product, comprising a computer program: when the computer program is executed by a processor, the imaging method of an ultrasonic beam image as described in any one of the above-mentioned first aspects is implemented.

[0047] The technical solution provided by the embodiments of the present application brings at least the following beneficial effects:

[0048] The embodiment of the present application obtains the ultrasonic array echo data received by the ultrasonic array element; then uses the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as the basic curves to perform interpolation processing on the ultrasonic array echo data to obtain three reference array apodization data; then uses the three reference array apodization data to perform apodization processing on the ultrasonic array echo data to obtain three reference array radio frequency identification RF data; performs cross-correlation processing on the three reference array RF data to obtain target array RF data, and performs beam synthesis and imaging based on the target array RF data to obtain a target ultrasonic beam image. Therefore, the present application improves the contrast resolution of the RF image after ultrasonic beam synthesis without affecting the main lobe, improves the contrast resolution of the image within the entire image range, has fewer algorithm parameters, lower computational complexity, and has good real-time performance.

[0049] Other features and advantages of the present application will be described in the following description, and partly become apparent from the description, or be understood by practicing the present application. The purpose and other advantages of the present application can be realized and obtained by the structures specifically pointed out in the written description, claims, and drawings. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0051] Figure 1 This is a schematic diagram of the structure of an ultrasonic device according to an embodiment of the present application;

[0052] Figure 2 A schematic diagram of the application principle of an ultrasonic device according to an embodiment of the present application;

[0053] Figure 3 An application scenario diagram of an ultrasonic beam imaging method provided in an embodiment of the present application;

[0054] Figure 4 A schematic diagram of a flow chart of an imaging method of an ultrasonic beam image provided in an embodiment of the present application;

[0055] Figure 5 A schematic diagram of a first alternating apodization curve provided in an embodiment of the present application;

[0056] Figure 6 A schematic diagram of a second alternating apodization curve provided in an embodiment of the present application;

[0057] Figure 7 A schematic flow chart of a method for determining reference array apodization data provided in an embodiment of the present application;

[0058] Figure 8 A schematic diagram of a flow chart of a method for determining reference array RF data provided in an embodiment of the present application;

[0059] Fig. 9 A schematic diagram of the transmission line provided in the embodiment of the present application at position 1;

[0060] Fig.10 A schematic diagram of the transmission line provided in the embodiment of the present application at position 2;

[0061] Fig.11 A flowchart of a method for determining RF data of a target array provided in an embodiment of the present application;

[0062] Fig.12 A schematic flow chart of a method for determining fourth reference array RF data provided in an embodiment of the present application;

[0063] Fig.13 A schematic flow chart of a method for determining a first correlation array coefficient provided in an embodiment of the present application;

[0064] Fig.14 A flowchart of a method for determining RF data of a target array provided in an embodiment of the present application;

[0065] Fig.15 A schematic flow chart of a method for determining a second correlation array coefficient provided in an embodiment of the present application;

[0066] Fig.16 A schematic diagram of a flow chart of an imaging method of an ultrasonic beam image provided in an embodiment of the present application;

[0067] Fig.17 A schematic diagram of the results of a simulation practice of an imaging method for ultrasonic beam images provided in an embodiment of the present application;

[0068] Fig.18 A schematic diagram of the simulation results of the rectangular window apodization imaging method provided in an embodiment of the present application;

[0069] Fig.19 A schematic diagram of the structure of an imaging device for ultrasonic beam images provided in an embodiment of the present application;

[0070] Fig. 20 A schematic diagram of the structure of an imaging device for ultrasonic beam images provided in an embodiment of the present application. DETAILED DESCRIPTION

[0071] In order to make the purpose, technical scheme and advantages of the embodiments of the present application clearer, the technical scheme in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Among them, the described embodiments are part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of this application.

[0072] Furthermore, in the description of the embodiments of the present application, unless otherwise specified, “ / ” means or. For example, A / B can mean A or B. The “and / or” in the text is merely a description of the association relationship of associated objects, indicating that three relationships may exist. For example, A and / or B can mean: A exists alone, A and B exist at the same time, and B exists alone. In addition, in the description of the embodiments of the present application, “multiple” refers to two or more than two.

[0073] In the following, the terms "first" and "second" are used for descriptive purposes only and should not be understood as suggesting or implying relative importance or implicitly indicating the number of technical features indicated. Thus, features defined as "first" and "second" may explicitly or implicitly include one or more of the features. In the description of the embodiments of the present application, unless otherwise specified, "multiple" means two or more.

[0074] Medical ultrasound examination is widely used to examine diseases in a series of parts such as the carotid artery, abdomen, and heart. These examinations rely on the quality of ultrasound images, and the quality of ultrasound images relies on the resolution of ultrasound images. The resolution of ultrasound images is mainly divided into spatial resolution and contrast resolution. However, the RF (Radio Frequency) image synthesized by the ultrasound beam is greatly affected by the side lobes, resulting in artifacts in the image, making the contrast resolution low. Therefore, reducing the influence of side lobes is one of the effective ways to improve the contrast resolution.

[0075] The existing ultrasound image sidelobe suppression is mainly achieved through the apodization window function, and the more commonly used ones are the Hanning window and the Hamming window. However, usually a single amplitude apodization function can only act on a local area of ​​the entire image, so that the image resolution in this area is better, and the image resolution outside this area is significantly reduced. This type of method will increase the width of the main lobe while reducing the sidelobe, resulting in a decrease in the lateral spatial resolution of the image. In order to balance the spatial resolution and contrast resolution, some scholars have proposed a method of using a more complex nonlinear apodization function to process RF data. This method can effectively solve the former problem, but it will pay for extremely high computational complexity. Therefore, how to use a lower complexity algorithm to improve the contrast resolution of the RF image after ultrasonic beam synthesis without affecting the main lobe is an urgent problem to be solved.

[0076] In view of this, an embodiment of the present application provides an imaging method, an ultrasound beam device and an apparatus to solve the problem of how to improve the contrast resolution of the RF image after ultrasound beam synthesis without affecting the main lobe.

[0077] The embodiment of the present application obtains the ultrasonic array echo data received by the ultrasonic array element; then uses the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as the basic curves to perform interpolation processing on the ultrasonic array echo data to obtain three reference array apodization data; then uses the three reference array apodization data to perform apodization processing on the ultrasonic array echo data to obtain three reference array radio frequency identification RF data; performs cross-correlation processing on the three reference array RF data to obtain target array RF data, and performs beam synthesis and imaging based on the target array RF data to obtain a target ultrasonic beam image. Therefore, the present application improves the contrast resolution of the RF image after ultrasonic beam synthesis without affecting the main lobe, improves the contrast resolution of the image within the entire image range, has fewer algorithm parameters, lower computational complexity, and has good real-time performance.

[0078] After introducing the main inventive concept of the embodiment of the present application, the ultrasonic device applied to the imaging method of an ultrasonic beam image provided by the embodiment of the present application is introduced in conjunction with the accompanying drawings. Figure 1 As shown, it is a structural block diagram of an ultrasonic device provided in an embodiment of the present application.

[0079] It should be understood that Figure 1 The ultrasound device shown is only an example and the ultrasound device may have more Figure 1The more or less components shown in the figure can be combined with two or more components, or can have different component configurations. The various components shown in the figure can be implemented in hardware, software, or a combination of hardware and software including one or more signal processing and / or application specific integrated circuits.

[0080] Figure 1 Schematic diagram of the hardware configuration of an ultrasound device is shown in FIG.

[0081] like Figure 1 As shown, the ultrasound device may include, for example: a processor 110, a memory 120, a display unit 130 and a probe 140; wherein:

[0082] A probe 140, used for transmitting an ultrasonic signal;

[0083] A display unit 130, used to display a target ultrasonic beam image;

[0084] The memory 120 is configured to store data required for generating a target ultrasonic beam image, which may include software programs, application interface data, etc.;

[0085] The processor 110 is connected to the probe 140, the display unit 130 and the memory 120 respectively, and is configured to execute: acquiring ultrasound array echo data received by the ultrasound array element; performing interpolation processing on the ultrasound array echo data using the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as basic curves respectively, to obtain three reference array apodization data; wherein the activated ultrasound array element represented by the first alternating apodization curve is the inactivated ultrasound array element represented by the second alternating apodization curve, and the inactivated ultrasound array element represented by the first alternating apodization curve is the activated ultrasound array element represented by the second alternating apodization curve; performing apodization processing on the ultrasound array echo data using the three reference array apodization data respectively, to obtain three reference array radio frequency identification RF data; performing cross-correlation processing on the three reference array RF data, to obtain target array RF data, and performing beam synthesis and imaging based on the target array RF data, to obtain a target ultrasonic beam image.

[0086] The target ultrasound beam image is displayed through the display unit 130 .

[0087] Figure 2 Schematic diagram of the application principle according to an embodiment of the present application. Figure 1 The partial modules or functional components of the ultrasonic device shown are implemented, and only the main components will be described below, while other components, such as memory, controller, control circuit, etc., will not be described here in detail.

[0088] like Figure 2As shown, the application environment may include a user interface 210 , a display unit 220 for displaying the user interface, and a processor 230 .

[0089] The display unit 220 may include a display panel 221 and a backlight assembly 222. The display panel 221 is configured to display an ultrasound image, the backlight assembly 222 is located behind the display panel 221, and the backlight assembly 222 may include a plurality of backlight partitions (not shown in the figure), each of which may emit light to illuminate the display panel 221.

[0090] The processor 230 may be configured to control the brightness of the backlight source of each backlight partition in the backlight assembly 222 , and to control the probe to transmit an ultrasonic signal, and to receive and analyze an ultrasonic echo signal to obtain an ultrasonic image.

[0091] The processor 230 may process the ultrasonic array echo data received by the ultrasonic array element to obtain a target ultrasonic beam image.

[0092] like Figure 3 , which is an application scenario diagram of an ultrasound beam imaging method provided by an embodiment of the present application. The diagram includes: an ultrasound device 10, a patient 20, and a memory 30;

[0093] The ultrasound device 10 is used to generate an echo signal by using acoustic radiation force or low-frequency vibration to excite a target part of the patient 20, obtain ultrasound array echo data received by the ultrasound array element, and store it in the memory 30;

[0094] Then, the ultrasound device 10 uses the ultrasound array echo data received by the ultrasound array elements stored in the memory 30, and uses the Hanning window function curve, the first alternating apodization curve, and the second alternating apodization curve as basic curves to perform interpolation processing on the ultrasound array echo data to obtain three reference array apodization data; wherein, the activated ultrasound array element represented by the first alternating apodization curve is the inactivated ultrasound array element represented by the second alternating apodization curve, and the inactivated ultrasound array element represented by the first alternating apodization curve is the activated ultrasound array element represented by the second alternating apodization curve; the ultrasound array echo data is apodized using the three reference array apodization data to obtain three reference array radio frequency identification RF data; the three reference array RF data are cross-correlated to obtain target array RF data, and beam synthesis and imaging are performed based on the target array RF data to obtain a target ultrasonic beam image, which is then displayed.

[0095] Of course, the method provided in the embodiment of the present application is not limited to Figure 3 The application scenarios shown can also be used in other possible application scenarios, and the embodiments of the present application are not limited thereto. Figure 3The functions that can be implemented by each device in the application scenario shown will be described in the subsequent method embodiments, and will not be described in detail here.

[0096] To further illustrate the technical solution provided by the embodiment of the present application, this is described in detail below in conjunction with the accompanying drawings and specific implementation methods. Although the embodiment of the present application provides the method operation steps shown in the following embodiments or drawings, more or fewer operation steps may be included in the method based on routine or no creative labor. In the steps where there is no necessary causal relationship logically, the execution order of these steps is not limited to the execution order provided in the embodiment of the present application.

[0097] See also Figure 4 , is a flow chart of an imaging method of an ultrasonic beam image provided in an embodiment of the present application. Figure 4 As shown, the method comprises the following steps:

[0098] In step 401, ultrasonic array echo data received by an ultrasonic array element is acquired.

[0099] The echo signals received by the ultrasonic array elements are amplified, AD (analog-to-digital) converted, and delayed focusing processed to obtain ultrasonic array echo data.

[0100] In a possible implementation, if there is only one transmission line, the data volume of the ultrasound array echo data in the embodiment of the present application is the product of the number of system channels in the ultrasound system and the number of sampling depths in the depth direction; if there are multiple transmission lines, the data volume of the ultrasound array echo data in the embodiment of the present application is the product of the number of system channels in the ultrasound system, the number of sampling depths in the depth direction, and the number of transmission lines.

[0101] Exemplarily, if the number of system channels in the ultrasound system is m and the number of sampling depths in the depth direction is n, the ultrasound array echo data is an m*n matrix.

[0102] In step 402, the ultrasound array echo data are interpolated using the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as basic curves to obtain three reference array apodization data.

[0103] The activated ultrasonic array element represented by the first alternating apodization curve is the inactivated ultrasonic array element represented by the second alternating apodization curve, and the inactivated ultrasonic array element represented by the first alternating apodization curve is the activated ultrasonic array element represented by the second alternating apodization curve. Assuming that the number of true elements activated each time is set to 2, the specific expression of the first alternating apodization curve is as follows: Figure 5 The specific expression of the second alternating trace curve is as follows: Figure 6From the curves shown, it can be seen that the positions of activated ultrasonic array elements of the first alternating apodization curve and the second alternating apodization curve are completely opposite, and the positions of unactivated ultrasonic array elements are also completely opposite.

[0104] In order to ensure the sampling accuracy and loudness of the three basic curves, namely the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve, the three basic curves may all be represented by the same number of sampling points, for example, 1024 sampling points.

[0105] In a possible implementation manner, the amount of ultrasound array echo data in the embodiment of the present application is the product of the number of system channels in the ultrasound system and the number of sampling depths in the depth direction; therefore, the Hanning window function curve, the first alternating apodization curve, and the second alternating apodization curve are used as basic curves to perform interpolation processing on the ultrasound array echo data to obtain three reference array apodization data, which can be performed as follows: Figure 7 Steps shown:

[0106] In step 701, the number of transmit apertures corresponding to each sampling depth in the depth direction is obtained.

[0107] In step 702, based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the first alternating apodization curve, first reference array apodization data with the first alternating apodization curve as a base curve is obtained.

[0108] In step 703, based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the second alternating apodization curve, second reference array apodization data using the second alternating apodization curve as a base curve is obtained.

[0109] In step 704, based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the Hanning window function, third reference array apodization data with the Hanning window function curve as a basic curve is obtained.

[0110] Among them, steps 701, 702, and 703 are executed in no particular order. Figure 7 The execution sequence shown may also be performed simultaneously, or 702 may be performed first, and then 701 and 703. Other execution sequences may also be implemented, and this is not limited in the implementation of this application.

[0111] Exemplarily, in order to ensure the symmetry of the Hanning window apodization and optimize the value of the alternating apodization window, in the embodiment of the present application, when setting the receiving aperture, it should be ensured that the aperture value of each depth is an even number. However, the number of transmitting apertures corresponding to each sampling depth in the depth direction may be the same or different. When calculating the apodization value at different sampling depths, for the first alternating apodization curve and the second alternating apodization curve, it is only necessary to intercept the number of transmitting apertures corresponding to the sampling depth from the center of the apodization curve to both sides to obtain the marginal value corresponding to each system channel in the ultrasound system corresponding to the sampling depth. In the data with the Hanning window function curve as the basic curve, the marginal value corresponding to each system channel in the ultrasound system corresponding to the sampling depth can be calculated using the Hanning window function, and the formula (1) of the Hanning window function is:

[0112]

[0113] Among them, Y is the marginal value corresponding to each system channel; n is the number of receiving apertures, which is an even number; k is the number of receiving apertures.

[0114] By calculating in the above manner, the marginal value corresponding to each system channel in the ultrasound system corresponding to each sampling depth in the first reference array apodization data using the first alternating apodization curve as the basic curve, the marginal value corresponding to each system channel in the ultrasound system corresponding to each sampling depth in the second reference array apodization data using the second alternating apodization curve as the basic curve, and the marginal value corresponding to each system channel in the ultrasound system corresponding to each sampling depth in the third reference array apodization data using the Hanning window function curve as the basic curve are obtained, thereby determining the first reference array apodization data, the second reference array apodization data, and the third reference array apodization data.

[0115] The first reference array apodization data, the second reference array apodization data, and the third reference array apodization data are all matrices formed by marginal values ​​corresponding to each system channel in the ultrasound system corresponding to each sampling depth.

[0116] In step 403, the ultrasound array echo data are apodized using the three reference array apodization data to obtain three reference array RF data.

[0117] In a possible implementation, the present application embodiment uses three reference array apodization data to perform apodization processing on the ultrasound array echo data to obtain three reference array RF data, which can be performed as follows: Figure 8 Steps shown:

[0118] In step 801, position information of a transmission line corresponding to ultrasound array echo data is determined.

[0119] In step 802, based on the position information of the transmission line, three reference array apodization data are respectively translated to obtain three target array apodization data.

[0120] In step 803, three reference array RF data are obtained based on the three target array apodization data and the ultrasound array echo data.

[0121] In a specific implementation, during beam synthesis, due to the presence of multiple transmission lines, it is necessary to align the apodization center, that is, a column of data in the center of the first reference array apodization data, the second reference array apodization data, and the third reference array apodization data with the position of the transmission line, and then translate the first reference array apodization data, the second reference array apodization data, and the third reference array apodization data to obtain three target array apodization data.

[0122] For example, taking the third reference array apodization data with the Hanning window function curve as the base curve as an example, if the position of the transmission line is at position 1, the generated target array apodization data and the base curve are as follows: Fig. 9 As shown; if the position of the transmission line moves from position 1 to position 2, the target array apodization data and basic curve generated are as follows Fig.10 shown.

[0123] After obtaining the three target array apodization data, the three target array apodization data and the ultrasonic array echo data are used to obtain the RF data corresponding to each sampling depth in the corresponding three reference array RF data through the following formula (2):

[0124]

[0125] Among them, RF(i) represents the RF data corresponding to the sampling depth i; ch is the chth system channel; N is the number of system channels; pod(ch,i) is the data of the chth system channel corresponding to the sampling depth i in the target array apodization data; AD(ch,i) is the data of the chth system channel corresponding to the sampling depth i in the ultrasonic array echo data.

[0126] After obtaining the RF data corresponding to each sampling depth in the corresponding three reference array RF data, the RF data corresponding to all sampling depths can form the reference array RF data. If the number of system channels in the ultrasound system is m, and the number of sampling depths in the depth direction is n, the ultrasound array echo data is an m*n matrix, and the target array apodization data is also an m*n matrix, and the reference array RF data obtained by calculation is a 1*n matrix.

[0127] In step 404, cross-correlation processing is performed on the three reference array RF data to obtain target array RF data, and beam synthesis and imaging are performed based on the target array RF data to obtain a target ultrasound beam image.

[0128] In a possible implementation manner, the three reference array RF data are cross-correlated to obtain the target array RF data, which can be performed as follows: Fig.11 Steps shown:

[0129] In step 1101, cross-correlation processing is performed on first reference array RF data using the first alternating apodization curve as a base curve and second reference array RF data using the second alternating apodization curve as a base curve to obtain fourth reference array RF data.

[0130] In a possible implementation, in the embodiment of the present application, cross-correlation processing is performed on the first reference array RF data using the first alternating apodization curve as the base curve and the second reference array RF data using the second alternating apodization curve as the base curve to obtain fourth reference array RF data, which can be performed as follows: Fig.12 Steps shown:

[0131] In step 1201, the RF data at the same position in the first reference array RF data and the second reference array RF data are added to obtain fifth reference array RF data, and a first correlation array coefficient is determined based on the first reference array RF data and the second reference array RF data;

[0132] In step 1202, fourth reference array RF data is obtained based on the fifth reference array RF data and the first correlation array coefficient.

[0133] In a possible implementation manner, the first correlation array coefficient is determined based on the first reference array RF data and the second reference array RF data in the embodiment of the present application, which can be performed as follows: Fig.13 Steps shown:

[0134] In step 1301, for each sampling depth in the depth direction: based on a preset number of RF data centered on the sampling depth in the first reference array RF data and a preset number of RF data centered on the sampling depth in the second reference array RF data, a reference coefficient factor corresponding to the sampling depth is determined; based on the reference coefficient factor corresponding to the sampling depth and a preset threshold, a target coefficient factor corresponding to the sampling depth is determined;

[0135] In step 1302, a first correlation array coefficient is determined based on a target coefficient factor corresponding to each sampling depth in the depth direction.

[0136] In specific implementation, assuming that the number of emission lines in the embodiment of the present application is 1, the preset number of longitudinal convolutions centered on the sampling depth is 2A+1, and the preset number corresponding to each sampling depth can be the same or different, and generally increases first and then decreases with the increase of the sampling depth. For example, it is 3 at the beginning, then gradually changes to 5, 7, 9, and then gradually decreases back to 3. Then the reference coefficient factor corresponding to each sampling depth can be calculated using formula (3):

[0137]

[0138] Wherein, Th(i,l) is the corresponding reference coefficient factor of sampling depth i and transmission line l; RF1(k,l) is the corresponding RF data of sampling depth k and transmission line l in the first reference array RF data; RF2(k,l) is the corresponding RF data of sampling depth k and transmission line l in the second reference array RF data; It is the sum of the data corresponding to the sampling depth iA to the sampling depth i+A.

[0139] After calculating the reference coefficient factor corresponding to the sampling depth, if the reference coefficient factor corresponding to the sampling depth is greater than the preset threshold, the reference coefficient factor corresponding to the sampling depth is used as the target coefficient factor corresponding to the sampling depth; if the reference coefficient factor corresponding to the sampling depth is not greater than the preset threshold, the preset threshold is used as the target coefficient factor corresponding to the sampling depth. After calculating the target coefficient factor corresponding to each sampling depth, the target coefficient factors corresponding to all sampling depths can form the first correlation array coefficient. If each reference array RF data is a 1*n matrix, the first correlation array coefficient is also a 1*n matrix. Then the fifth reference array RF data and the RF data at the same position in the first correlation array coefficient are superimposed to obtain the fourth reference array RF data.

[0140] In step 1102, cross-correlation processing is performed on the third reference array RF data using the Hanning window function curve as a basic curve and the fourth reference array RF data to obtain target array RF data.

[0141] In a possible implementation manner, in the embodiment of the present application, the third reference array RF data using the Hanning window function curve as the base curve and the fourth reference array RF data are cross-correlated to obtain the target array RF data, which can be performed as follows: Fig.14 Steps shown:

[0142] In step 1401, a smaller value is taken for the RF data at the same position in the third reference array RF data and the fourth reference array RF data to obtain sixth reference array RF data, and a second correlation array coefficient is determined based on the third reference array RF data and the fourth reference array RF data;

[0143] In step 1402, target array RF data is obtained based on the sixth reference array RF data and the second correlation array coefficient.

[0144] In a possible implementation manner, the second correlation array coefficient is determined based on the third reference array RF data and the fourth reference array RF data in the embodiment of the present application, which can be performed as follows: Fig.15 Steps shown:

[0145] In step 1501, for each sampling depth in the depth direction: based on a preset number of RF data centered at the sampling depth in the third reference array RF data and a preset number of RF data centered at the sampling depth in the fourth reference array RF data, a reference coefficient factor corresponding to the sampling depth is obtained; based on the reference coefficient factor corresponding to the sampling depth and a preset threshold, a target coefficient factor corresponding to the sampling depth is determined;

[0146] In step 1502, a second correlation array coefficient is determined based on a target coefficient factor corresponding to each sampling depth in the depth direction.

[0147] In specific implementation, assuming that the number of emission lines in the embodiment of the present application is 1, the preset number of longitudinal convolutions centered on the sampling depth is 2A+1, and the preset number corresponding to each sampling depth can be the same or different, and generally increases first and then decreases with the increase of the sampling depth. For example, it is 3 at the beginning, then gradually changes to 5, 7, 9, and then gradually decreases back to 3. Then the reference coefficient factor corresponding to each sampling depth can be calculated using formula (3):

[0148]

[0149] Wherein, Th(i,l) is the corresponding reference coefficient factor of the sampling depth i and the transmission line l; RF1(k,l) is the corresponding RF data of the sampling depth k and the transmission line l in the third reference array RF data; RF2(k,l) is the corresponding RF data of the sampling depth k and the transmission line l in the fourth reference array RF data; It is the sum of the data corresponding to the sampling depth iA to the sampling depth i+A.

[0150] After calculating the reference coefficient factor corresponding to the sampling depth, if the reference coefficient factor corresponding to the sampling depth is greater than the preset threshold, the reference coefficient factor corresponding to the sampling depth is used as the target coefficient factor corresponding to the sampling depth; if the reference coefficient factor corresponding to the sampling depth is not greater than the preset threshold, the preset threshold is used as the target coefficient factor corresponding to the sampling depth. After calculating the target coefficient factor corresponding to each sampling depth, the target coefficient factors corresponding to all sampling depths can form the second correlation array coefficient. If each reference array RF data is a 1*n matrix, the second correlation array coefficient is also a 1*n matrix. Then the sixth reference array RF data and the RF data at the same position in the second correlation array coefficient are superimposed to obtain the target array RF data.

[0151] For ease of understanding, the following Fig.16 The complete process of an ultrasonic beam imaging method provided in an embodiment of the present application is described. Fig.16 As shown, the method comprises the following steps:

[0152] In step 1601, the echo signal received by the ultrasonic array element is amplified, AD (analog-to-digital) converted, and delayed focusing processed to obtain ultrasonic array echo data.

[0153] In step 1602, the ultrasound array echo data are interpolated using the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as basic curves to obtain three reference array apodization data.

[0154] In step 1603, the position information of the transmission line corresponding to the ultrasound array echo data is determined, and based on the position information of the transmission line, three reference array apodization data are respectively translated to obtain three target array apodization data.

[0155] In step 1604, three reference array RF data are obtained based on the three target array apodization data and the ultrasound array echo data.

[0156] In step 1605, cross-correlation processing is performed on the first reference array RF data using the first alternating apodization curve as a base curve and the second reference array RF data using the second alternating apodization curve as a base curve to obtain fourth reference array RF data.

[0157] In step 1606, cross-correlation processing is performed on the third reference array RF data using the Hanning window function curve as a basic curve and the fourth reference array RF data to obtain target array RF data.

[0158] During the specific implementation, the ultrasonic array echo data is generated by the Vantage (Advantage) ultrasonic research system, which is a set of 5 target point data arranged longitudinally and spaced 10mm apart. The Vantage ultrasonic research system is an ultrasonic experimental simulation platform developed by the American Verasonics based on acoustic principles. In order to verify the effectiveness of the proposed algorithm, a set of ultrasonic array echo data is generated using the Vantage ultrasonic research system. In the point target simulation experiment, a row of 5 point targets with a longitudinal interval of 10mm are set, evenly distributed between the depths of 10mm and 50mm. The emission fixed-point focusing and receiving dynamic focusing methods are adopted, and the emission focus is fixed at 30mm. The center frequency of the array element used is 7.6MHz, the number of array elements is 128, the spacing is 0.3mm, the sampling frequency is 31.25MHz, and the sound speed is 1450m / s. The imaging method of the ultrasonic beam image provided in the embodiment of the present application is used for simulation practice, and the following is obtained: Fig.17 The results are shown in the figure. At the same time, the rectangular window apodization imaging method is used for simulation practice, and the following results are obtained: Fig.18 The results shown clearly show that without affecting the main lobe, Fig.17 The contrast resolution of the RF image after ultrasound beam synthesis is significantly higher than that of Fig.18 The RF images shown have higher contrast resolution.

[0159] Based on the foregoing description, the embodiment of the present application obtains ultrasonic array echo data received by an ultrasonic array element; then uses the Hanning window function curve, the first alternating apodization curve, and the second alternating apodization curve as basic curves to perform interpolation processing on the ultrasonic array echo data to obtain three reference array apodization data; then uses the three reference array apodization data to perform apodization processing on the ultrasonic array echo data to obtain three reference array radio frequency identification RF data; performs cross-correlation processing on the three reference array RF data to obtain target array RF data, and performs beam synthesis and imaging based on the target array RF data to obtain a target ultrasonic beam image.

[0160] Therefore, the present application improves the contrast resolution of the RF image after ultrasonic beam synthesis without affecting the main lobe, improves the contrast resolution of the image within the entire image range, has fewer algorithm parameters, lower computational complexity, and better real-time performance.

[0161] According to the same inventive concept, an embodiment of the present application provides an imaging device for ultrasonic beam images. Fig.19 FIG. 1 is a schematic diagram of an imaging device for ultrasonic beam images provided in an embodiment of the present application, comprising:

[0162] An acquisition module 1901 is used to acquire ultrasound array echo data received by ultrasound array elements;

[0163] an interpolation module 1902, configured to perform interpolation processing on the ultrasound array echo data using a Hanning window function curve, a first alternating apodization curve, and a second alternating apodization curve as basic curves, respectively, to obtain three reference array apodization data; wherein the activated ultrasound array element represented by the first alternating apodization curve is the inactivated ultrasound array element represented by the second alternating apodization curve, and the inactivated ultrasound array element represented by the first alternating apodization curve is the activated ultrasound array element represented by the second alternating apodization curve;

[0164] An apodization module 1903 is used to perform apodization processing on the ultrasound array echo data using the three reference array apodization data respectively, so as to obtain three reference array radio frequency identification RF data;

[0165] The synthesis module 1904 is used to perform cross-correlation processing on the three reference array RF data to obtain target array RF data, and perform beam synthesis and imaging based on the target array RF data to obtain a target ultrasonic beam image.

[0166] In a possible implementation, the data volume of the ultrasound array echo data is the product of the number of system channels in the ultrasound system and the number of sampling depths in the depth direction;

[0167] When performing the interpolation processing on the ultrasound array echo data by using the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as the basic curves to obtain three reference array apodization data, the interpolation module 1902 is specifically used for:

[0168] Obtaining the number of transmit apertures corresponding to each sampling depth in the depth direction;

[0169] obtaining first reference array apodization data using the first alternating apodization curve as a base curve based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the first alternating apodization curve; and

[0170] obtaining second reference array apodization data using the second alternating apodization curve as a base curve based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the second alternating apodization curve; and

[0171] Based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the Hanning window function, third reference array apodization data with the Hanning window function curve as a basic curve is obtained.

[0172] In a possible implementation manner, when performing apodization processing on the ultrasound array echo data by respectively using the three reference array apodization data to obtain three reference array RF data, as follows: Fig. 20 As shown, the apodization module 1903 further includes:

[0173] A determination submodule 2001 is used to determine the position information of the transmission line corresponding to the ultrasound array echo data;

[0174] A translation submodule 2002 is used to translate the three reference array apodization data respectively based on the position information of the transmission line to obtain three target array apodization data;

[0175] The data determination submodule 2003 is used to obtain the three reference array RF data based on the three target array apodization data and the ultrasound array echo data.

[0176] In a possible implementation manner, when performing the cross-correlation processing on the three reference array RF data to obtain the target array RF data, the synthesis module 1904 is specifically configured to:

[0177] performing cross-correlation processing on first reference array RF data using the first alternating apodization curve as a basic curve and second reference array RF data using the second alternating apodization curve as a basic curve to obtain fourth reference array RF data;

[0178] The third reference array RF data using the Hanning window function curve as a basic curve and the fourth reference array RF data are cross-correlated to obtain the target array RF data.

[0179] In a possible implementation manner, when performing cross-correlation processing on the first reference array RF data using the first alternating apodization curve as a base curve and the second reference array RF data using the second alternating apodization curve as a base curve to obtain fourth reference array RF data, the synthesis module 1904 is specifically configured to:

[0180] Adding the RF data at the same position in the first reference array RF data and the second reference array RF data to obtain fifth reference array RF data, and determining a first correlation array coefficient based on the first reference array RF data and the second reference array RF data;

[0181] The fourth reference array RF data is obtained based on the fifth reference array RF data and the first correlation array coefficient.

[0182] In a possible implementation manner, when determining the first correlation array coefficient based on the first reference array RF data and the second reference array RF data, the synthesis module 1904 is specifically configured to:

[0183] For each sampling depth in the depth direction: determining a reference coefficient factor corresponding to the sampling depth based on a preset number of RF data centered on the sampling depth in the first reference array RF data and a preset number of RF data centered on the sampling depth in the second reference array RF data; determining a target coefficient factor corresponding to the sampling depth based on the reference coefficient factor corresponding to the sampling depth and a preset threshold;

[0184] The first correlation array coefficient is determined based on a target coefficient factor corresponding to each sampling depth in the depth direction.

[0185] In a possible implementation, when performing cross-correlation processing on the third reference array RF data using the Hanning window function curve as a base curve and the fourth reference array RF data to obtain the target array RF data, the synthesis module 1904 is specifically used to:

[0186] performing smaller value processing on the RF data at the same position in the third reference array RF data and the fourth reference array RF data to obtain sixth reference array RF data, and determining a second correlation array coefficient based on the third reference array RF data and the fourth reference array RF data;

[0187] The target array RF data is obtained based on the sixth reference array RF data and the second correlation array coefficient.

[0188] In a possible implementation, when determining the second correlation array coefficient based on the third reference array RF data and the fourth reference array RF data, the synthesis module 1904 is specifically configured to:

[0189] For each sampling depth in the depth direction: based on a preset number of RF data centered at the sampling depth in the third reference array RF data and the preset number of RF data centered at the sampling depth in the fourth reference array RF data, a reference coefficient factor corresponding to the sampling depth is obtained; based on the reference coefficient factor corresponding to the sampling depth and a preset threshold, a target coefficient factor corresponding to the sampling depth is determined;

[0190] The second correlation array coefficient is determined based on the target coefficient factor corresponding to each sampling depth in the depth direction.

[0191] In an exemplary embodiment, a computer-readable storage medium including instructions is also provided, such as a memory including instructions, and the instructions can be executed by a processor to complete the above-mentioned ultrasound beam image imaging method. Optionally, the storage medium can be a non-transitory computer-readable storage medium, for example, the non-transitory computer-readable storage medium can be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, etc.

[0192] In an exemplary embodiment, a computer program product is also provided, comprising a computer program, wherein when the computer program is executed by a processor, any one of the imaging methods of ultrasound beam images provided in the present application is implemented.

[0193] In an exemplary embodiment, various aspects of an imaging method for ultrasonic beam images provided by the present application may also be implemented in the form of a program product, which includes program code. When the program product is run on a computer device, the program code is used to enable the computer device to execute the steps of the imaging method for ultrasonic beam images according to various exemplary embodiments of the present application described above in this specification.

[0194] The program product may use any combination of one or more readable media. The readable medium may be a readable signal medium or a readable storage medium. The readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination of the above. More specific examples of readable storage media (a non-exhaustive list) include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0195] The program product of the imaging method for ultrasonic beam images of the embodiment of the present application can adopt a portable compact disk read-only memory (CD-ROM) and include program code, and can be run on an electronic device. However, the program product of the present application is not limited to this. In this document, a readable storage medium can be any tangible medium containing or storing a program, which can be used by or in combination with an instruction execution system, an apparatus or a device.

[0196] The readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, wherein the readable program code is carried. Such propagated data signals may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. The readable signal medium may also be any readable medium other than a readable storage medium, which may send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0197] The program code included in the readable medium can be transmitted with any appropriate medium, including but not limited to wireless, wired, optical cable, RF, etc., or any suitable combination of the above. The program code for performing the operation of the present application can be written in any combination of one or more programming languages, and the programming language includes object-oriented programming languages ​​such as Java, C++, etc., and also includes conventional procedural programming languages ​​such as "like" language or similar programming languages. The program code can be executed completely on the user electronic device, partially on the user device, as an independent software package, partially on the user electronic device, partially on the remote electronic device, or completely on the remote electronic device or service end. In the case of remote electronic devices, the remote electronic device can be connected to the user electronic device by any type of network including a local area network (LAN) or a wide area network (WAN), or can be connected to an external electronic device (for example, using an Internet service provider to connect to the Internet).

[0198] It should be noted that, although several units or subunits of the device are mentioned in the above detailed description, this division is merely exemplary and not mandatory. In fact, according to the embodiments of the present application, the features and functions of two or more units described above can be embodied in one unit. Conversely, the features and functions of one unit described above can be further divided into multiple units to be embodied.

[0199] In addition, although the operations of the method of the present application are described in a specific order in the drawings, this does not require or imply that the operations must be performed in this specific order, or that all the operations shown must be performed to achieve the desired results. Additionally or alternatively, some steps may be omitted, multiple steps may be combined into one step, and / or one step may be decomposed into multiple steps.

[0200] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0201] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0202] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0203] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0204] Obviously, those skilled in the art can make various changes and modifications to the present application without departing from the spirit and scope of the present application. Thus, if these modifications and variations of the present application fall within the scope of the claims of the present application and their equivalents, the present application is also intended to include these modifications and variations.

Claims

1. An imaging method of an ultrasonic beam image, characterized in that: The method comprises: Acquiring ultrasonic array echo data received by ultrasonic array elements; Using a Hanning window function curve, a first alternating apodization curve, and a second alternating apodization curve as basic curves, respectively, interpolation processing is performed on the ultrasonic array echo data to obtain three reference array apodization data; wherein the activated ultrasonic array element represented by the first alternating apodization curve is the inactivated ultrasonic array element represented by the second alternating apodization curve, and the inactivated ultrasonic array element represented by the first alternating apodization curve is the activated ultrasonic array element represented by the second alternating apodization curve; Using the three reference array apodization data respectively, the ultrasound array echo data is apodized to obtain three reference array radio frequency identification (RF) data; The three reference array radio frequency identification (RF) data are cross-correlated to obtain target array RF data, and beam synthesis and imaging are performed based on the target array RF data to obtain a target ultrasonic beam image.

2. The method according to claim 1, characterized in that The data volume of the ultrasound array echo data is the product of the number of system channels in the ultrasound system and the number of sampling depths in the depth direction; The method uses the Hanning window function curve, the first alternating apodization curve and the second alternating apodization curve as basic curves to perform interpolation processing on the ultrasound array echo data to obtain three reference array apodization data, including: Obtaining the number of transmit apertures corresponding to each sampling depth in the depth direction; obtaining first reference array apodization data using the first alternating apodization curve as a base curve based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the first alternating apodization curve; and obtaining second reference array apodization data using the second alternating apodization curve as a base curve based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the second alternating apodization curve; and Based on the number of transmit apertures corresponding to each sampling depth in the depth direction and the Hanning window function, third reference array apodization data with the Hanning window function curve as a basic curve is obtained.

3. The method according to claim 1, characterized in that The method of using the three reference array apodization data to perform apodization processing on the ultrasound array echo data to obtain three reference array radio frequency identification (RF) data comprises: Determining position information of a transmission line corresponding to the ultrasound array echo data; Based on the position information of the transmission line, the three reference array apodization data are respectively translated to obtain three target array apodization data; The three reference array radio frequency identification (RF) data are obtained based on the three target array apodization data and the ultrasound array echo data.

4. The method according to claim 1, characterized in that: The cross-correlation processing of the three reference array RF radio frequency identification data to obtain the target array RF data includes: performing cross-correlation processing on first reference array RF data using the first alternating apodization curve as a basic curve and second reference array RF data using the second alternating apodization curve as a basic curve to obtain fourth reference array RF data; The third reference array RF data using the Hanning window function curve as a basic curve and the fourth reference array RF data are cross-correlated to obtain the target array RF data.

5. The method according to claim 4, characterized in that The method of performing cross-correlation processing on first reference array RF data using the first alternating apodization curve as a basic curve and second reference array RF data using the second alternating apodization curve as a basic curve to obtain fourth reference array RF data includes: Adding the RF data at the same position in the first reference array RF data and the second reference array RF data to obtain fifth reference array RF data, and determining a first correlation array coefficient based on the first reference array RF data and the second reference array RF data; The fourth reference array RF data is obtained based on the fifth reference array RF data and the first correlation array coefficient.

6. The method according to claim 5, characterized in that The determining of a first correlation array coefficient based on the first reference array RF data and the second reference array RF data comprises: For each sampling depth in the depth direction: determining a reference coefficient factor corresponding to the sampling depth based on a preset number of RF data centered on the sampling depth in the first reference array RF data and a preset number of RF data centered on the sampling depth in the second reference array RF data; determining a target coefficient factor corresponding to the sampling depth based on the reference coefficient factor corresponding to the sampling depth and a preset threshold; The first correlation array coefficient is determined based on a target coefficient factor corresponding to each sampling depth in the depth direction.

7. The method according to claim 4, characterized in that The cross-correlation processing is performed on the third reference array RF data using the Hanning window function curve as a basic curve and the fourth reference array RF data to obtain the target array RF data, including: performing smaller value processing on the RF data at the same position in the third reference array RF data and the fourth reference array RF data to obtain sixth reference array RF data, and determining a second correlation array coefficient based on the third reference array RF data and the fourth reference array RF data; The target array RF data is obtained based on the sixth reference array RF data and the second correlation array coefficient.

8. The method according to claim 7, characterized in that The determining of second correlation array coefficients based on the third reference array RF data and the fourth reference array RF data comprises: For each sampling depth in the depth direction: based on a preset number of RF data centered at the sampling depth in the third reference array RF data and the preset number of RF data centered at the sampling depth in the fourth reference array RF data, a reference coefficient factor corresponding to the sampling depth is obtained; based on the reference coefficient factor corresponding to the sampling depth and a preset threshold, a target coefficient factor corresponding to the sampling depth is determined; The second correlation array coefficient is determined based on the target coefficient factor corresponding to each sampling depth in the depth direction.

9. An ultrasonic device, characterized in that: The ultrasound device includes at least one processor and at least one memory; wherein the memory stores a program code, and when the program code is executed by the processor, the processor performs the following process: Acquiring ultrasonic array echo data received by ultrasonic array elements; Using a Hanning window function curve, a first alternating apodization curve, and a second alternating apodization curve as basic curves, respectively, interpolation processing is performed on the ultrasonic array echo data to obtain three reference array apodization data; wherein the activated ultrasonic array element represented by the first alternating apodization curve is the inactivated ultrasonic array element represented by the second alternating apodization curve, and the inactivated ultrasonic array element represented by the first alternating apodization curve is the activated ultrasonic array element represented by the second alternating apodization curve; Using the three reference array apodization data respectively, the ultrasound array echo data is apodized to obtain three reference array radio frequency identification (RF) data; The three reference array radio frequency identification (RF) data are cross-correlated to obtain target array RF data, and beam synthesis and imaging are performed based on the target array RF data to obtain a target ultrasonic beam image.

10. An imaging device for ultrasonic beam images, characterized in that: The device comprises: An acquisition module, used for acquiring ultrasound array echo data received by ultrasound array elements; an interpolation module, configured to perform interpolation processing on the ultrasonic array echo data using a Hanning window function curve, a first alternating apodization curve, and a second alternating apodization curve as basic curves, respectively, to obtain three reference array apodization data; wherein the activated ultrasonic array element represented by the first alternating apodization curve is the inactivated ultrasonic array element represented by the second alternating apodization curve, and the inactivated ultrasonic array element represented by the first alternating apodization curve is the activated ultrasonic array element represented by the second alternating apodization curve; an apodization module, configured to perform apodization processing on the ultrasound array echo data using the three reference array apodization data respectively, to obtain three reference array radio frequency identification RF data; The synthesis module is used to perform cross-correlation processing on the three reference array radio frequency identification RF data to obtain target array RF data, and perform beam synthesis and imaging based on the target array RF data to obtain a target ultrasonic beam image.

Images