A shear wave image generation method, ultrasonic equipment and device

Abstract

The present application discloses a shear wave image generation method, ultrasonic equipment and device. The present application utilizes acoustic radiation force or low-frequency vibration to excite the target part to generate shear waves, detects the shear waves on the shear wave propagation path, and obtains the detection data of the shear waves; determines the particle vibration velocities corresponding to each particle of the shear wave at different times; for the particle vibration velocities corresponding to any particle at different times, performs filtering processing on the particle vibration velocities corresponding to the particles at different times in multiple directions, obtains the reference particle vibration velocities corresponding to each particle vibration velocity in each direction, and obtains the shear wave velocity corresponding to each particle in any direction; and generates a shear wave image based on the preset ultrasonic black and white image and the shear wave velocity corresponding to each particle in each direction. In this way, the directionality and accuracy of the shear wave propagation velocity can be improved, providing a more comprehensive shear wave information.

Description

Technical Field

[0001] The present application relates to the technical field of shear wave elastic imaging, and in particular to a shear wave image generation method, ultrasonic equipment and device. Background Art

[0002] Anisotropy is a directionally dependent property found in fibrous biological tissues. When soft tissues such as skeletal muscle, cardiac muscle, and tendons are subjected to external forces, their mechanical properties exhibit a certain degree of anisotropy. In ultrasound elastography, the propagation velocity of directional shear waves is used to describe the mechanical properties of the soft tissue.

[0003] However, current shear wave imaging generally assumes that shear waves propagate isotropically in biological soft tissues. The shear wave velocity between two shear wave particles is then calculated to estimate the mechanical properties of the tissue. This method yields a single shear wave velocity, resulting in low accuracy in the resulting shear wave propagation velocity. Summary of the Invention

[0004] The purpose of this application is to provide a shear wave image generation method, ultrasonic equipment and device to solve the problem of low accuracy of shear wave propagation velocity information represented in the prior art.

[0005] In a first aspect, the present application provides a shear wave image generation method, the method comprising:

[0006] Using acoustic radiation force or low-frequency vibration to excite a target part of the object to generate a shear wave, detecting the shear wave along a shear wave propagation path, and acquiring detection data of the shear wave;

[0007] Based on the detection data, determining the particle vibration speed corresponding to each particle of the shear wave at different times;

[0008] For the particle vibration velocities corresponding to any particle at different times, filtering is performed in multiple directions on the particle vibration velocities corresponding to the particle at different times to obtain reference particle vibration velocities corresponding to each particle vibration velocity in each direction;

[0009] For any direction, determining the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction;

[0010] A shear wave image is generated based on the preset ultrasonic black and white image and the shear wave velocity corresponding to each particle in each direction.

[0011] In a possible implementation, determining the particle vibration velocities corresponding to each particle of the shear wave at different times based on the detection data includes:

[0012] For any particle, based on the data corresponding to the particle in the detection data, a two-dimensional autocorrelation algorithm is used to determine the axial displacement of the particle within a preset time period corresponding to each moment;

[0013] Based on the axial displacement within the preset time period corresponding to each moment and the preset time period, the particle vibration speed corresponding to each particle at different moments is determined.

[0014] In a possible implementation, performing filtering processing on the particle vibration velocities corresponding to the particle at different times in multiple directions to obtain the reference particle vibration velocity corresponding to each particle vibration velocity in each direction includes:

[0015] Performing Fourier transform on the particle vibration velocities corresponding to the particles at different times to obtain frequency domain and wave number domain data corresponding to the vibration velocity of each particle;

[0016] For any frequency domain and wave number domain data, the frequency domain and wave number domain data are multiplied by the preset masks corresponding to the multiple directions respectively; and each product is subjected to inverse Fourier transform to obtain the reference particle vibration velocity corresponding to the frequency domain and wave number domain data in each direction.

[0017] In a possible implementation, the preset mask is a mask obtained by truncating edges using a truncation function.

[0018] In a possible implementation, determining the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction includes:

[0019] For any particle in the direction, perform the following operations:

[0020] determining a distance between the mass point and the candidate mass points based on the position information of the mass point and the position information of the candidate mass points; the candidate mass points are one or more continuous mass points adjacent to the mass point in the direction;

[0021] The shear wave velocity corresponding to the particle is determined based on the reference particle vibration velocities of the particle and the candidate particle in the direction and the distance between the particle and the candidate particle.

[0022] In a possible implementation, determining the shear wave velocity corresponding to the particle based on the reference particle vibration velocities of the particle and the candidate particle in the direction, and the distance between the particle and the candidate particle, includes:

[0023] determining a time difference between the first time and the second time based on a first time corresponding to a peak point of a reference mass point vibration velocity of the mass point in the direction and a second time corresponding to a peak point of a reference mass point vibration velocity of the candidate mass point in the direction;

[0024] The shear wave velocity corresponding to the particle point is obtained based on the distance between the particle point and the candidate particle point, and the time difference between the first time and the second time.

[0025] In one possible implementation, the generating of shear waves by exciting a target portion of an object with acoustic radiation force or low-frequency vibration, detecting the shear waves along a shear wave propagation path, and acquiring detection data of the shear waves includes:

[0026] Using acoustic radiation force or low-frequency vibration to excite the target part of the object to generate shear waves, and transmitting ultrasonic signals multiple times along the shear wave propagation path to detect the shear waves and obtain echo data;

[0027] The echo data is subjected to orthogonal demodulation and frequency reduction processing to obtain detection data of the shear wave.

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

[0029] Using acoustic radiation force or low-frequency vibration to excite a target part of the object to generate a shear wave, detecting the shear wave along a shear wave propagation path, and acquiring detection data of the shear wave;

[0030] Based on the detection data, determining the particle vibration speed corresponding to each particle of the shear wave at different times;

[0031] For the particle vibration velocities corresponding to any particle at different times, filtering is performed in multiple directions on the particle vibration velocities corresponding to the particle at different times to obtain reference particle vibration velocities corresponding to each particle vibration velocity in each direction;

[0032] For any direction, determining the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction;

[0033] A shear wave image is generated based on the preset ultrasonic black and white image and the shear wave velocity corresponding to each particle in each direction.

[0034] In a possible implementation, when performing filtering processing on the particle vibration velocities corresponding to the particle at different times in multiple directions to obtain a reference particle vibration velocity corresponding to each particle vibration velocity in each direction, the processor is specifically configured to:

[0035] Performing Fourier transform on the particle vibration velocities corresponding to the particles at different times to obtain frequency domain and wave number domain data corresponding to the vibration velocity of each particle;

[0036] For any frequency domain and wave number domain data, the frequency domain and wave number domain data are multiplied by the preset masks corresponding to the multiple directions respectively; and each product is subjected to inverse Fourier transform to obtain the reference particle vibration velocity corresponding to the frequency domain and wave number domain data in each direction.

[0037] In a third aspect, an embodiment of the present application provides a shear wave image generating device, the device comprising:

[0038] an acquisition module, configured to generate shear waves by exciting a target part of an object using acoustic radiation force or low-frequency vibration, detect the shear waves along a shear wave propagation path, and acquire detection data of the shear waves;

[0039] a particle vibration velocity determination module, configured to determine the particle vibration velocities corresponding to the respective particles of the shear wave at different times based on the detection data;

[0040] a reference particle vibration velocity determination module, configured to filter the particle vibration velocities corresponding to any particle at different times in multiple directions to obtain the reference particle vibration velocities corresponding to each particle vibration velocity in each direction;

[0041] A shear wave velocity determination module is used to determine, for any direction, the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction;

[0042] The image generation module is used to generate a shear wave image based on a preset ultrasonic black and white image and the shear wave velocity corresponding to each particle in each direction.

[0043] 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 shear wave image generation method as described in any one of the first aspects above.

[0044] 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 shear wave image generation method as described in any one of the first aspects above is implemented.

[0045] The technical solutions provided by the embodiments of this application bring at least the following beneficial effects:

[0046] The embodiment of the present application utilizes acoustic radiation force or low-frequency vibration to excite the target part of the object to generate shear waves, detects the shear waves on the shear wave propagation path, and obtains the detection data of the shear waves; based on the detection data, determines the particle vibration velocities corresponding to each particle of the shear wave at different times; for the particle vibration velocities corresponding to any particle at different times, performs filtering processing in multiple directions on the particle vibration velocities corresponding to the particle at different times to obtain the reference particle vibration velocities corresponding to each particle vibration velocity in each direction; for any direction, determines the shear wave velocity corresponding to each particle in the direction based on the reference particle vibration velocities of each particle in the direction; generates a shear wave image based on a preset ultrasonic black and white image and the shear wave velocities corresponding to each particle in each direction.

[0047] Therefore, by displaying and superimposing the information of shear wave propagation velocity in different directions to generate the final shear wave image, the accuracy of shear wave propagation velocity is improved, and the problem of low accuracy of shear wave propagation velocity in a single dimension (only size information) is solved. A more comprehensive shear wave information is provided, which can greatly improve the work efficiency of clinicians and has important clinical significance for disease diagnosis.

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

[0049] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following is a brief introduction to the drawings required for use in the embodiments of the present application. 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 creative work.

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

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

[0052] Figure 3 This is an application scenario diagram of a shear wave image generation method provided in an embodiment of the present application;

[0053] Figure 4 A schematic flow chart of a shear wave image generation method provided in an embodiment of the present application;

[0054] Figure 5 A schematic diagram of preset masks corresponding to multiple directions provided in an embodiment of the present application;

[0055] Figure 6 A schematic diagram of the vibration velocities of a reference particle corresponding to multiple directions provided in an embodiment of the present application;

[0056] Figure 7 A schematic flow chart of a method for determining the shear wave velocity corresponding to each particle provided in an embodiment of the present application;

[0057] Figure 8 A schematic flow chart of a method for determining the shear wave velocity corresponding to each particle provided in an embodiment of the present application;

[0058] Figure 9 A schematic diagram of a mass point and candidate mass points provided in an embodiment of the present application;

[0059] Figure 10 A schematic diagram of the vibration velocities of the reference particles corresponding to the particles and candidate particles at different times provided in an embodiment of the present application;

[0060] Figure 11 A schematic diagram of a shear wave image in one direction provided in an embodiment of the present application;

[0061] Figure 12 A schematic diagram of a shear wave image in another direction provided by an embodiment of the present application;

[0062] Figure 13 A schematic diagram of a shear wave image in another direction provided by an embodiment of the present application;

[0063] Figure 14 A schematic diagram of a shear wave image in another direction provided by an embodiment of the present application;

[0064] Figure 15 A schematic diagram of a shear wave image in another direction provided by an embodiment of the present application;

[0065] Figure 16 A schematic diagram of a shear wave image in another direction provided by an embodiment of the present application;

[0066] Figure 17A schematic diagram of a shear wave image in another direction provided by an embodiment of the present application;

[0067] Figure 18 A schematic diagram of a shear wave image in another direction provided by an embodiment of the present application;

[0068] Figure 19 A shear wave image including shear wave velocities in various directions provided in an embodiment of the present application;

[0069] Figure 20 A complete flowchart of a shear wave image generation method provided in an embodiment of the present application;

[0070] Figure 21 A schematic structural diagram of a shear wave image generating device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0071] To make the purpose, technical solutions, and advantages of the embodiments of the present application clearer, the technical solutions 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, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0072] Moreover, in the description of the embodiments of the present application, unless otherwise specified, “ / ” means or, for example, A / B can mean A or B; “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 to imply or suggest relative importance or implicitly indicate the number of the technical features indicated. Therefore, features defined as "first" or "second" may explicitly or implicitly include one or more of the features. In the description of the embodiments of this application, unless otherwise specified, "plurality" means two or more.

[0074] The embodiment of the present application utilizes acoustic radiation force or low-frequency vibration to excite the target part of the object to generate shear waves, detects the shear waves on the shear wave propagation path, and obtains the detection data of the shear waves; based on the detection data, the particle vibration velocities corresponding to each particle of the shear wave at different times are determined; for the particle vibration velocities corresponding to any particle at different times, the particle vibration velocities corresponding to the particle at different times are filtered in multiple directions to obtain the reference particle vibration velocities corresponding to each particle vibration velocity in each direction; for any direction, the shear wave velocity corresponding to each particle in the direction is determined based on the reference particle vibration velocities of each particle in the direction; a shear wave image is generated based on a preset ultrasonic black and white image and the shear wave velocity corresponding to each particle in each direction. Thus, by displaying and superimposing the information of the shear wave propagation velocity in different directions to generate the final shear wave image, the accuracy of the shear wave propagation velocity is improved, the problem of low accuracy of the shear wave propagation velocity of a single dimension (only size information) is solved, and a more comprehensive shear wave information is provided, which can greatly improve the work efficiency of clinicians and has important clinical significance for disease diagnosis.

[0075] After introducing the main inventive concept of the embodiment of the present application, the ultrasonic device to which the shear wave image generation method provided by the embodiment of the present application is applied is introduced below with reference to the accompanying drawings. Figure 1 , which is a structural block diagram of an ultrasound device provided in an embodiment of the present application.

[0076] It should be understood that Figure 1 The ultrasound device shown is only an example and the ultrasound device may have more Figure 1 The 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.

[0077] Figure 1 exemplarily shows a hardware configuration block diagram of an ultrasound device.

[0078] 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:

[0079] Probe 140, used for transmitting ultrasonic signals;

[0080] A display unit 130, configured to display a shear wave image;

[0081] The memory 120 is configured to store data required for shear wave image generation, which may include software programs, application interface data, etc.;

[0082] The processor 110 is respectively connected to the probe 140, the display unit 130 and the memory 120, and is configured to perform: using acoustic radiation force or low-frequency vibration to excite the target part of the object to generate shear waves, detecting the shear waves on the shear wave propagation path, and obtaining the detection data of the shear waves; based on the detection data, determining the particle vibration velocities corresponding to each particle of the shear wave at different times; for the particle vibration velocities corresponding to any particle at different times, filtering the particle vibration velocities corresponding to the particle at different times in multiple directions to obtain the reference particle vibration velocities corresponding to each particle vibration velocity in each direction; for any direction, determining the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocities of each particle in the direction; generating a shear wave image according to a preset ultrasonic black and white image and the shear wave velocities corresponding to each particle in each direction.

[0083] The display unit 130 displays shear wave images in various directions generated according to the shear wave velocity corresponding to each particle in various directions, and the final shear wave image generated according to the preset ultrasonic black and white image and the shear wave velocity corresponding to each particle in various directions.

[0084] Figure 2 Schematic diagram of the application principle according to an embodiment of the present application. Figure 1 The implementation of some modules or functional components of the ultrasound device shown will only be described below with respect to the main components, while other components, such as memory, controller, control circuit, etc., will not be described in detail here.

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

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

[0087] 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.

[0088] The processor 230 may process the echo signal of the shear wave to obtain a shear wave image.

[0089] like Figure 3FIG. 1 is an application scenario diagram of a shear wave image generation method provided by an embodiment of the present application. The diagram includes: an ultrasound device 10, a patient 20, and a memory 30;

[0090] The ultrasound device 10 is used to generate shear waves on a target part of a patient 20 by using acoustic radiation force or low-frequency vibration, detect the shear waves along the shear wave propagation path, obtain detection data of the shear waves, and store the data in a memory 30;

[0091] The ultrasonic device 10 then uses the detection data stored in the memory 30 to determine the particle vibration velocities corresponding to each particle of the shear wave at different times; for the particle vibration velocities corresponding to any particle at different times, the particle vibration velocities corresponding to the particle at different times are filtered in multiple directions to obtain the reference particle vibration velocities corresponding to each particle vibration velocity in each direction; for any direction, the shear wave velocity corresponding to each particle in the direction is determined based on the reference particle vibration velocities of each particle in the direction; a shear wave image is generated based on a preset ultrasonic black and white image and the shear wave velocities corresponding to each particle in each direction, and is displayed.

[0092] 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 3 The functions that can be implemented by each device in the application scenario shown will be described in subsequent method embodiments and will not be described in detail here.

[0093] To further illustrate the technical solutions provided by the embodiments of the present application, the following is a detailed description of the technical solutions in conjunction with the accompanying drawings and specific implementation methods. Although the embodiments of the present application provide the method operation steps as shown in the following embodiments or drawings, more or fewer operation steps may be included in the method based on routine or no creative work. In 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 embodiments of the present application.

[0094] See also Figure 4 , is a flow chart of a shear wave image generation method provided in an embodiment of the present application. Figure 4 As shown, the method includes the following steps:

[0095] In step 401, acoustic radiation force or low-frequency vibration is used to excite a target part of an object to generate shear waves, and the shear waves are detected along a shear wave propagation path to obtain detection data of the shear waves.

[0096] In one possible embodiment, acoustic radiation force or low-frequency vibration is used to excite the target part of the object to generate shear waves, the shear waves are detected on the shear wave propagation path, and the detection data of the shear waves are obtained. This can be performed as follows: acoustic radiation force or low-frequency vibration is used to excite the target part of the object to generate shear waves, and ultrasonic signals are emitted multiple times on the shear wave propagation path to detect the shear waves and obtain echo data; the echo data is orthogonally demodulated and down-converted to obtain the detection data of the shear waves.

[0097] Among them, after performing orthogonal demodulation processing and frequency reduction processing on the echo data, IQ data can be obtained, and the IQ data is the detection data of the shear wave.

[0098] In step 402, based on the detection data, the particle vibration velocities corresponding to each particle of the shear wave at different times are determined.

[0099] In one possible embodiment, based on the detection data, the particle vibration velocities corresponding to each particle of the shear wave at different moments are determined, which can be performed as follows: for any particle, based on the data corresponding to the particle in the detection data, a two-dimensional autocorrelation algorithm is used to determine the axial displacement of the particle within the preset time length corresponding to each moment; based on the axial displacement within the preset time length corresponding to each moment and the preset time length, the particle vibration velocities corresponding to the particle at different moments are determined.

[0100] Specifically, the echo data is subjected to orthogonal demodulation and frequency reduction processing to obtain IQ data and each particle on each shear wave detection line in the horizontal direction. Then, for any particle on any shear wave detection line, within the preset axial range of the particle, multiple pre-set vertical and time sampling points are used to calculate the axial displacement of the particle within the preset time period at each moment according to the following formula:

[0101]

[0102] in, is the axial displacement of the particle within the preset time length corresponding to a moment; M is the number of sample points in the vertical direction; N is the number of sample points in the time direction; m is the coordinate of each sample point in the vertical direction; n is the coordinate of each sample point in the time direction; c is the propagation speed of sound at the target part; π is the circumference of a circle; f c is the center frequency of the echo data signal; I and Q are the in-phase and quadrature components of the echo data signal, respectively, and are obtained after orthogonal demodulation and frequency reduction processing of the echo data signal.

[0103] Then, after using the above formula to calculate each particle, the axial displacement of each particle within the preset time length corresponding to each moment will be obtained. Then, based on the axial displacement within the preset time length corresponding to each moment and the preset time length, the particle vibration speed corresponding to the moment can be determined according to the following formula:

[0104]

[0105] in, It represents the vibration velocity of the particle at that moment. represents the axial displacement of the particle within the preset time period corresponding to that moment, and Δt represents the preset time period corresponding to that moment. The above formula indicates that the particle's vibration velocity at that moment can be obtained by dividing the axial displacement of the particle within the preset time period corresponding to that moment by the preset time period. It also indicates that by differentiating the axial displacement of the particle within the preset time period corresponding to that moment, the particle's vibration velocity at that moment can be obtained.

[0106] In step 403, for the particle vibration velocity corresponding to any particle at different times, the particle vibration velocity corresponding to the particle at different times is filtered in multiple directions to obtain the reference particle vibration velocity corresponding to each particle vibration velocity in each direction.

[0107] In one possible implementation, filtering processes are performed on the particle vibration velocities corresponding to the particles at different times in multiple directions to obtain the reference particle vibration velocities corresponding to each particle vibration velocity in each direction. This can be performed as follows:

[0108] The particle vibration velocities corresponding to the particles at different times are subjected to Fourier transform respectively to obtain the frequency domain and wave number domain data corresponding to each particle vibration velocity; for any frequency domain and wave number domain data, the frequency domain and wave number domain data are multiplied by the preset masks corresponding to multiple directions respectively; and the product is subjected to inverse Fourier transform to obtain the reference particle vibration velocity corresponding to the frequency domain and wave number domain data in each direction.

[0109] In a possible implementation manner, in order to prevent the Gibbs effect, the preset mask in the embodiment of the present application is a mask after the edge is truncated using a truncation function.

[0110] In specific implementation, the particle vibration velocities corresponding to each particle at different times are Fourier transformed respectively, and the time-space domain data are converted into frequency domain-wave number domain data to obtain the frequency domain-wave number domain data corresponding to each particle vibration velocity. Then, according to the direction of interest to the user, all the frequency domain-wave number domain data of each particle are multiplied by the preset mask corresponding to the direction of interest to the user to suppress the shear waves propagating in other directions. Then, the product is inverse Fourier transformed to obtain the reference particle vibration velocity corresponding to each particle at different times in the direction of interest to the user.

[0111] Among them, in the embodiment of the present application, when all frequency domain and wavenumber domain data of each particle are multiplied by the preset mask corresponding to the direction of the user's interest according to the direction of the user's interest, the mask after the edge is truncated by the truncation function can be used as the preset mask, or the mask without the edge being truncated by the truncation function can be used as the preset mask. The embodiment of the present application does not impose any restrictions on this, and one of them can be selected for use according to actual conditions.

[0112] For example, Figure 5 As shown in FIG. 1 , a schematic diagram of preset masks corresponding to multiple directions provided in an embodiment of the present application can be used to respectively compare the frequency domain wave number domain data with the frequency domain wave number domain data. Figure 5 The preset masks corresponding to 0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315° are multiplied respectively, and then the inverse Fourier transform is performed to obtain the following: Figure 6 The reference particle vibration speed corresponding to 0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315° is shown in FIG. Figure 5 The preset masks corresponding to 0° are multiplied respectively, and then the inverse Fourier transform is performed to obtain Figure 6 The reference particle vibration velocity corresponding to 0° after filtering is shown. The process of obtaining the reference particle vibration velocity in other directions is the same as that in the 0° direction, and will not be described in detail here. The multiple directions are set to 0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315°, respectively. Other directions can also be set according to the direction of user interest, and this embodiment of the application does not limit this.

[0113] In one possible implementation, in order to eliminate isolated particles, in the embodiment of the present application, before performing filtering processing in multiple different directions on the particle vibration velocities corresponding to the particles at different times, the particle vibration velocities corresponding to any particle at different times can also be low-pass filtered and median filtered.

[0114] Thus, interference between shear waves can be eliminated, so that the time delay between adjacent particles of the shear wave in each direction can be calculated more accurately in the following step 404 .

[0115] In step 404, for any direction, the shear wave velocity corresponding to each particle in the direction is determined based on the reference particle vibration velocity of each particle in the direction.

[0116] In one possible implementation, the shear wave velocity corresponding to each particle in the direction is determined based on the reference particle vibration velocity of each particle in the direction. The shear wave velocity can be performed for each particle in the direction. Figure 7 Steps shown:

[0117] In step 701, the distance between the particle and candidate particle points is determined based on the position information of the particle and the position information of the candidate particle points; the candidate particle points are one or more continuous particle points adjacent to the particle in a direction.

[0118] In step 702, the shear wave velocity corresponding to the particle is determined based on the reference particle vibration velocity in the direction of the particle and the candidate particle, and the distance between the particle and the candidate particle.

[0119] In one possible implementation, the shear wave velocity corresponding to the particle is determined based on the reference particle vibration velocity of the particle and the candidate particle in the direction, and the distance between the particle and the candidate particle, which can be performed as follows: Figure 8 Steps shown:

[0120] In step 801, based on a first time corresponding to a peak point of a reference particle vibration velocity in a particle direction and a second time corresponding to a peak point of a reference particle vibration velocity in a candidate particle direction, a time difference between the first time and the second time is determined;

[0121] In step 802, the shear wave velocity corresponding to the particle point is obtained based on the distance between the particle point and the candidate particle point, and the time difference between the first time and the second time.

[0122] For example, Figure 9 As shown in the figure, the shear wave velocity corresponding to the particle A in the 90° direction is calculated. If the candidate particle is a particle adjacent to the particle A in the 90° direction, first a particle adjacent to the particle A in the 90° direction is selected, that is, the particle B. Then, based on the position information of the particle A and the position information of the particle B, the distance d between the particle A and the particle B is determined. Then, based on the following equation: Figure 10 The reference particle vibration speeds of the particle A and the particle B in the 90° direction are shown, and the first time t corresponding to the peak point of the reference particle vibration speed of the particle A in the 90° direction is determined. A, and the second time t corresponding to the peak point of the reference particle vibration velocity of particle B in the 90° direction B , using the principle of cross-correlation to calculate the time delay of particle A in the 90° direction, that is, to determine the time difference t between the first time and the second time B -t A , and finally based on the distance d between mass point A and mass point B, and the time difference t between the first time and the second time B -t A , the shear wave velocity corresponding to particle A is obtained according to the following formula:

[0123]

[0124] Exemplarily, the shear wave velocity corresponding to particle A in the 90° direction is calculated. If the candidate particles are multiple continuous particles adjacent to particle A in the 90° direction, then multiple continuous particles adjacent to particle A in the 90° direction are selected, such as particle B, particle C, and particle D. Then, the distances d1, d2, and d3 between particle A and particle B, particle C, and particle D are determined respectively. At the same time, based on the reference particle vibration velocities of particle A, particle B, particle C, and particle D in the 90° direction, the time differences t1, t2, and t3 between particle A and particle B, particle C, and particle D are determined respectively. Then, the above formula is used to determine the multiple shear wave velocities corresponding to particle A as follows: Finally, the average of V1, V2, and V3 is used as the final shear wave velocity corresponding to particle A.

[0125] Thus, the shear wave velocity corresponding to the particle A with direction information can be calculated, so that the shear wave velocity has two characteristic information: magnitude and direction.

[0126] In step 405, a shear wave image is generated according to a preset ultrasonic black-and-white image and the shear wave velocity corresponding to each particle in each direction.

[0127] In specific implementation, after calculating the shear wave velocity corresponding to each particle in each direction, the shear wave velocity corresponding to each particle in each direction can be displayed after being superimposed on the preset ultrasonic black and white image according to the vector information (direction information and size information). Figure 11 The shear wave image in the 0° direction is shown. Figure 12 The shear wave image in the 45° direction is shown. Figure 13 The shear wave image in the 90° direction is shown. Figure 14 The shear wave image in the 135° direction is shown. Figure 15 The shear wave image in the 180° direction is shown. Figure 16 The shear wave image in the 225° direction is shown. Figure 17The shear wave image in the 270° direction is shown. Figure 18 The shear wave image in the 315° direction shown in the figure can finally be Figures 11-18 The shear wave images corresponding to each direction are fused according to the weight corresponding to each direction, and the result is as follows: Figure 19 The shear wave image shown in Figure 1 contains the shear wave velocities in all directions. Figures 11-19 The images shown are all stored in the ultrasound device, and the shear wave image can be displayed according to the user's needs. For example, if the user needs to view the shear wave image in the 90° direction, he can select Figure 13 For example, if the user needs to view the shear wave images in various directions, he can select Figure 19 The shear wave image shown is displayed.

[0128] exist Figures 11-19 In the figure, the length of each arrow represents the shear wave velocity corresponding to each particle. The longer the arrow, the greater the shear wave velocity corresponding to the particle. The direction of each arrow represents the direction of the shear wave velocity corresponding to each particle.

[0129] In another example, the length and color of the arrow can be used together to represent the shear wave velocity corresponding to each particle. The longer the arrow and the higher the RGB value of the color, the greater the shear wave velocity corresponding to the particle. The direction of the arrow represents the direction of the shear wave velocity corresponding to each particle.

[0130] For ease of understanding, the following Figure 20 The complete process of a shear wave image generation method provided in an embodiment of the present application is described. Figure 20 As shown, the method includes the following steps:

[0131] In step 2001, acoustic radiation force or low-frequency vibration is used to excite a target part of an object to generate shear waves, and the shear waves are detected along a shear wave propagation path to obtain detection data of the shear waves.

[0132] In step 2002, for any particle, based on the data corresponding to the particle in the detection data, a two-dimensional autocorrelation algorithm is used to determine the axial displacement of the particle within the preset time length corresponding to each moment; based on the axial displacement within the preset time length corresponding to each moment and the preset time length, the particle vibration speed corresponding to the particle at different moments is determined.

[0133] In step 2003, Fourier transform is performed on the particle vibration velocities corresponding to any particle at different times to obtain frequency domain and wave number domain data corresponding to each particle vibration velocity.

[0134] In step 2004, for any frequency domain and wave number domain data, the frequency domain and wave number domain data are multiplied by preset masks corresponding to multiple directions respectively; and the product is subjected to inverse Fourier transform to obtain the reference particle vibration velocity corresponding to the frequency domain and wave number domain data in each direction.

[0135] In step 2005, for any particle in any direction, the distance between the particle and the candidate particle is determined based on the position information of the particle and the position information of the candidate particle; based on the first time corresponding to the peak point of the reference particle vibration velocity of the particle in the direction, and the second time corresponding to the peak point of the reference particle vibration velocity of the candidate particle in the direction, the time difference between the first time and the second time is determined.

[0136] In step 2006, the shear wave velocity corresponding to the particle point is obtained based on the distance between the particle point and the candidate particle point, and the time difference between the first time and the second time.

[0137] In step 2007, a shear wave image is generated based on a preset ultrasonic black-and-white image and the shear wave velocity corresponding to each particle in each direction.

[0138] Based on the foregoing description, the embodiment of the present application utilizes acoustic radiation force or low-frequency vibration to excite the target part of the object to generate shear waves, detects the shear waves on the shear wave propagation path, and obtains the detection data of the shear waves; based on the detection data, determines the particle vibration velocities corresponding to each particle of the shear wave at different times; for the particle vibration velocities corresponding to any particle at different times, performs filtering processing in multiple directions on the particle vibration velocities corresponding to the particles at different times to obtain the reference particle vibration velocities corresponding to each particle vibration velocity in each direction; for any direction, determines the shear wave velocity corresponding to each particle in the direction based on the reference particle vibration velocities of each particle in the direction; generates a shear wave image based on a preset ultrasonic black and white image and the shear wave velocities corresponding to each particle in each direction.

[0139] Therefore, by displaying and superimposing the information of shear wave propagation velocity in different directions to generate the final shear wave image, the accuracy of shear wave propagation velocity is improved, and the problem of low accuracy of shear wave propagation velocity in a single dimension (only size information) is solved. A more comprehensive shear wave information is provided, which can greatly improve the work efficiency of clinicians and has important clinical significance for disease diagnosis.

[0140] Based on the same inventive concept, the embodiment of the present application provides a shear wave image generating device. Figure 21 FIG. 1 is a schematic diagram of a shear wave image generating device provided in an embodiment of the present application, comprising:

[0141] An acquisition module 2101 is configured to generate shear waves by exciting a target part of an object using acoustic radiation force or low-frequency vibration, detect the shear waves along a shear wave propagation path, and acquire detection data of the shear waves;

[0142] a particle vibration velocity determination module 2102 for determining the particle vibration velocity corresponding to each particle of the shear wave at different times based on the detection data;

[0143] The reference particle vibration velocity determination module 2103 is configured to perform filtering processing on the particle vibration velocity corresponding to any particle at different times in multiple directions to obtain the reference particle vibration velocity corresponding to each particle vibration velocity in each direction;

[0144] A shear wave velocity determination module 2104 is configured to determine, for any direction, the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction;

[0145] The image generation module 2105 is used to generate a shear wave image according to a preset ultrasonic black and white image and the shear wave velocity corresponding to each particle in each direction.

[0146] In a possible implementation, when determining the particle vibration velocities corresponding to each particle of the shear wave at different times based on the detection data, the particle vibration velocity determination module 2102 is specifically configured to:

[0147] For any particle, based on the data corresponding to the particle in the detection data, a two-dimensional autocorrelation algorithm is used to determine the axial displacement of the particle within a preset time period corresponding to each moment;

[0148] Based on the axial displacement within the preset time period corresponding to each moment and the preset time period, the particle vibration speed corresponding to each particle at different moments is determined.

[0149] In one possible implementation, when performing filtering processing on the particle vibration velocities corresponding to the particle at different times in multiple directions to obtain the reference particle vibration velocity corresponding to each particle vibration velocity in each direction, the reference particle vibration velocity determination module 2103 is specifically configured to:

[0150] Performing Fourier transform on the particle vibration velocities corresponding to the particles at different times to obtain frequency domain and wave number domain data corresponding to the vibration velocity of each particle;

[0151] For any frequency domain and wave number domain data, the frequency domain and wave number domain data are multiplied by the preset masks corresponding to the multiple directions respectively; and each product is subjected to inverse Fourier transform to obtain the reference particle vibration velocity corresponding to the frequency domain and wave number domain data in each direction.

[0152] In a possible implementation, the preset mask is a mask obtained by truncating edges using a truncation function.

[0153] In a possible implementation, when determining the shear wave velocity corresponding to each particle in the direction based on the reference particle vibration velocity of each particle in the direction, the shear wave velocity determination module 2104 is specifically configured to:

[0154] For any particle in the direction, perform the following operations:

[0155] determining a distance between the mass point and the candidate mass points based on the position information of the mass point and the position information of the candidate mass points; the candidate mass points are one or more continuous mass points adjacent to the mass point in the direction;

[0156] The shear wave velocity corresponding to the particle is determined based on the reference particle vibration velocities of the particle and the candidate particle in the direction and the distance between the particle and the candidate particle.

[0157] In one possible implementation, when determining the shear wave velocity corresponding to the particle based on the reference particle vibration velocity of the particle and the candidate particle in the direction, and the distance between the particle and the candidate particle, the shear wave velocity determination module 2104 is specifically configured to:

[0158] determining a time difference between the first time and the second time based on a first time corresponding to a peak point of a reference mass point vibration velocity of the mass point in the direction and a second time corresponding to a peak point of a reference mass point vibration velocity of the candidate mass point in the direction;

[0159] The shear wave velocity corresponding to the particle point is obtained based on the distance between the particle point and the candidate particle point, and the time difference between the first time and the second time.

[0160] In one possible implementation, when performing the step of using acoustic radiation force or low-frequency vibration to excite a target portion of an object to generate a shear wave, detecting the shear wave along a shear wave propagation path, and acquiring detection data of the shear wave, the acquisition module 2101 is specifically configured to:

[0161] Using acoustic radiation force or low-frequency vibration to excite the target part of the object to generate shear waves, and transmitting ultrasonic signals multiple times along the shear wave propagation path to detect the shear waves and obtain echo data;

[0162] The echo data is subjected to orthogonal demodulation and frequency reduction processing to obtain detection data of the shear wave.

[0163] In an exemplary embodiment, a computer-readable storage medium including instructions is further provided, such as a memory including instructions, wherein the instructions are executable by a processor to perform the shear wave image generation method. Alternatively, the storage medium may be a non-transitory computer-readable storage medium, such as a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, or the like.

[0164] In an exemplary embodiment, a computer program product is further provided, comprising a computer program, wherein when the computer program is executed by a processor, any of the shear wave image generation methods provided in the present application is implemented.

[0165] In an exemplary embodiment, various aspects of a shear wave image generation method provided by the present application can 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 shear wave image generation method according to various exemplary embodiments of the present application described above in this specification.

[0166] The program product may employ 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 component, or any combination thereof. 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 thereof.

[0167] The program product for the shear wave image generation method of the embodiment of the present application can be in 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 thereto, and in this document, a readable storage medium can be any tangible medium containing or storing a program that can be used by or in combination with an instruction execution system, apparatus, or device.

[0168] A readable signal medium may include a data signal transmitted in baseband or as part of a carrier wave, which carries readable program code. Such a transmitted data signal may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A readable signal medium may also be any readable medium other than a readable storage medium that can transmit, propagate, or transfer a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0169] 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 this 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 and partially on the remote electronic device, or completely on the remote electronic device or the server. In the case of a remote electronic device, 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 via the Internet).

[0170] It should be noted that although several units or subunits of the device are mentioned in the detailed description above, this division is merely exemplary and not mandatory. In fact, depending on the embodiment of the application, the features and functions of two or more units described above can be embodied in a single unit. Conversely, the features and functions of a single unit described above can be further divided and embodied by multiple units.

[0171] Furthermore, although the operations of the method of the present application are described in a particular order in the accompanying drawings, this does not require or imply that the operations must be performed in this particular order, or that all illustrated operations 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.

[0172] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can 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 can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0173] 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 block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks 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.

[0174] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work 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 The function specified in one or more boxes.

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

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

Claims

1. A shear wave image generation method, characterized in that: The method comprises: Using acoustic radiation force or low-frequency vibration to excite a target part of the object to generate a shear wave, detecting the shear wave along a shear wave propagation path, and acquiring detection data of the shear wave; Based on the detection data, determining the particle vibration speed corresponding to each particle of the shear wave at different times; For the particle vibration velocities corresponding to any particle at different times, Fourier transform is performed on the particle vibration velocities corresponding to the particle at different times to obtain frequency domain wave number domain data corresponding to each particle vibration velocity; for any frequency domain wave number domain data, the frequency domain wave number domain data is multiplied by preset masks corresponding to multiple directions; and each product is inverse Fourier transformed to obtain the reference particle vibration velocity corresponding to the frequency domain wave number domain data in each direction; For any direction, determining the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction; The shear wave velocity corresponding to each particle in each direction is superimposed on the preset ultrasonic black and white image according to the direction information and size information to obtain the shear wave images in each direction; the shear wave images in each direction are fused according to the weight corresponding to each direction to obtain a shear wave image containing the shear wave velocities in each direction.

2. The method according to claim 1, characterized in that The determining, based on the detection data, particle vibration velocities corresponding to the respective particles of the shear wave at different times, includes: For any particle, based on the data corresponding to the particle in the detection data, a two-dimensional autocorrelation algorithm is used to determine the axial displacement of the particle within a preset time period corresponding to each moment; Based on the axial displacement within the preset time period corresponding to each moment and the preset time period, the particle vibration speed corresponding to each particle at different moments is determined.

3. The method according to claim 1, characterized in that The preset mask is a mask after the edge is truncated using a truncation function.

4. The method according to claim 1, wherein Determining the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction includes: For any particle in the direction, perform the following operations: determining a distance between the mass point and the candidate mass points based on the position information of the mass point and the position information of the candidate mass points; the candidate mass points are one or more continuous mass points adjacent to the mass point in the direction; The shear wave velocity corresponding to the particle is determined based on the reference particle vibration velocities of the particle and the candidate particle in the direction and the distance between the particle and the candidate particle.

5. The method according to claim 4, characterized in that The determining of the shear wave velocity corresponding to the particle based on the reference particle vibration velocity of the particle and the candidate particle in the direction and the distance between the particle and the candidate particle includes: determining a time difference between the first time and the second time based on a first time corresponding to a peak point of a reference mass point vibration velocity of the mass point in the direction and a second time corresponding to a peak point of a reference mass point vibration velocity of the candidate mass point in the direction; The shear wave velocity corresponding to the particle point is obtained based on the distance between the particle point and the candidate particle point, and the time difference between the first time and the second time.

6. The method according to claim 1, characterized in that The method comprises: utilizing acoustic radiation force or low-frequency vibration to excite a target part of an object to generate a shear wave, detecting the shear wave along a shear wave propagation path, and acquiring detection data of the shear wave, including: Using acoustic radiation force or low-frequency vibration to excite the target part of the object to generate shear waves, and transmitting ultrasonic signals multiple times along the shear wave propagation path to detect the shear waves and obtain echo data; The echo data is subjected to orthogonal demodulation and frequency reduction processing to obtain detection data of the shear wave.

7. An ultrasonic device, characterized in that The ultrasound device includes at least one processor and at least one memory; wherein the memory stores program code, and when the program code is executed by the processor, the processor performs the following process: Using acoustic radiation force or low-frequency vibration to excite a target part of the object to generate a shear wave, detecting the shear wave along a shear wave propagation path, and acquiring detection data of the shear wave; Based on the detection data, determining the particle vibration speed corresponding to each particle of the shear wave at different times; For the particle vibration velocities corresponding to any particle at different times, Fourier transform is performed on the particle vibration velocities corresponding to the particle at different times to obtain frequency domain wave number domain data corresponding to each particle vibration velocity; for any frequency domain wave number domain data, the frequency domain wave number domain data is multiplied by preset masks corresponding to multiple directions; and each product is inverse Fourier transformed to obtain the reference particle vibration velocity corresponding to the frequency domain wave number domain data in each direction; For any direction, determining the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction; The shear wave velocity corresponding to each particle in each direction is superimposed on the preset ultrasonic black and white image according to the direction information and size information to obtain the shear wave images in each direction; the shear wave images in each direction are fused according to the weight corresponding to each direction to obtain a shear wave image containing the shear wave velocities in each direction.

8. A shear wave image generating device, characterized in that: The device comprises: an acquisition module, configured to generate shear waves by exciting a target part of an object using acoustic radiation force or low-frequency vibration, detect the shear waves along a shear wave propagation path, and acquire detection data of the shear waves; a particle vibration velocity determination module, configured to determine the particle vibration velocities corresponding to the respective particles of the shear wave at different times based on the detection data; A reference particle vibration velocity determination module is configured to perform Fourier transform on the particle vibration velocities corresponding to any particle at different times to obtain frequency-domain wavenumber-domain data corresponding to each particle vibration velocity; for any frequency-domain wavenumber-domain data, multiply the frequency-domain wavenumber-domain data by preset masks corresponding to multiple directions; and perform inverse Fourier transform on each product to obtain the reference particle vibration velocity corresponding to the frequency-domain wavenumber-domain data in each direction; A shear wave velocity determination module is used to determine, for any direction, the shear wave velocity corresponding to each particle in the direction according to the reference particle vibration velocity of each particle in the direction; The image generation module is used to superimpose the shear wave velocity corresponding to each particle in each direction according to the direction information and size information and the preset ultrasonic black and white image to obtain the shear wave images in each direction; and fuse the shear wave images in each direction according to the weight corresponding to each direction to obtain a shear wave image containing the shear wave velocity in each direction.

Images