Ultrasonic elastic imaging method, device, electronic equipment and storage medium

Abstract

The embodiments of the present application provide an ultrasonic elastic imaging method, device, electronic device and storage medium, which obtain ultrasonic echo signals by applying vibration excitation to the tissue to be tested and performing ultrasonic detection; perform beam synthesis based on the ultrasonic echo signals to generate synthetic ultrasonic imaging data; use the synthetic ultrasonic imaging data to generate ultrasonic detection images and elastic detection images, wherein the elastic detection image is used to characterize tissue elasticity detection information and the signal quality of the ultrasonic detection image at the corresponding moment; and display the ultrasonic detection image and the elastic detection image in a superimposed manner based on the positional relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image. After the ultrasonic detection image is generated, the generated elastic detection image is superimposed on the ultrasonic detection image and displayed, and the signal quality of the ultrasonic detection image is characterized by the elastic detection image, so that the operator can quickly control the detection results and improve the detection efficiency and accuracy.

Description

Technical Field

[0001] The present application relates to the technical field of medical equipment, and in particular to an ultrasonic elastic imaging method, device, electronic equipment and storage medium. Background Art

[0002] Currently, ultrasound imaging technology is a widely used inspection technology in clinical testing. It detects the physiological structure of biological tissues through ultrasound and realizes the detection of lesions in biological tissues. However, in the early stage of biological tissue lesions, since its structure does not change significantly, traditional ultrasound imaging methods such as B-ultrasound are not sensitive to early lesions in biological tissues. Based on ultrasound imaging technology, elastic imaging technology, which has emerged in recent years, is a new imaging technology. It applies shear waves to biological tissues to stimulate the vibration of biological tissues, and then uses ultrasound to track the morphological changes of the measured tissues before and after vibration stress, thereby obtaining the mechanical characteristics of the measured tissues, and detecting the state of the biological tissues based on the mechanical characteristics of the biological tissues to determine whether there are early lesion characteristics.

[0003] In the prior art, the results of elasticity calculation are usually in the form of parameter values, combined with ultrasound images to form elasticity imaging results, and the elasticity imaging results are displayed to the operator of the equipment. The operator needs to observe the ultrasound image and combine it with the quantitative parameter value of the elasticity measurement to realize the state detection of the measured tissue.

[0004] However, since the best test position can only be determined after multiple observations and adjustments during the elasticity measurement process, the ultrasonic elasticity imaging method in the prior art has the problems of low detection efficiency and poor detection accuracy. Summary of the invention

[0005] The present application provides an ultrasonic elastic imaging method, device, electronic device and storage medium to solve the problems of low detection efficiency and poor detection accuracy in the prior art.

[0006] According to the first aspect of an embodiment of the present application, the present application provides an ultrasonic elasticity imaging method, the method comprising: applying vibration excitation to a tissue to be tested and performing ultrasonic detection to obtain an ultrasonic echo signal; performing beam synthesis based on the ultrasonic echo signal to generate synthetic ultrasonic imaging data; using the synthetic ultrasonic imaging data to generate an ultrasonic detection image and an elasticity detection image, wherein the elasticity detection image is used to characterize tissue elasticity detection information and the signal quality of the ultrasonic detection image at a corresponding moment; based on the positional relationship between a test area corresponding to the ultrasonic detection image and a test area corresponding to the elasticity detection image, the ultrasonic detection image and the elasticity detection image are superimposed and displayed.

[0007] In one possible implementation, vibration excitation is applied to the tissue to be tested and ultrasonic detection is performed to obtain an ultrasonic echo signal, including: applying a shear wave to the tissue to be tested, and collecting an ultrasonic echo signal for tracking the propagation of the shear wave; performing beam synthesis based on the ultrasonic echo signal to generate synthetic ultrasonic imaging data, including: performing beam synthesis on the ultrasonic echo signal for tracking the propagation of the shear wave to generate multiple frames of scan lines, wherein the multiple frames of scan lines are used to characterize the ultrasonic detection results at different positions.

[0008] In a possible implementation, the use of the synthetic ultrasonic imaging data to generate an ultrasonic detection image and an elastic detection image includes: processing the synthetic ultrasonic imaging data through a graphics processor to generate the ultrasonic detection image and the elastic detection image in parallel, and obtaining an elastic modulus E according to a shear wave propagation velocity, wherein E = 3ρV 2 , ρ is tissue density, and V is shear wave velocity.

[0009] In a possible implementation, after the ultrasonic detection image is generated and before it is displayed, the ultrasonic detection image is also processed as follows: median filtering is performed to remove electronic noise; Gaussian smoothing and non-local mean filtering are performed to smooth the structural information; bilateral filtering is performed to enhance the edge information; two thresholds of the signal are used as display thresholds, the lower threshold is the signal at an amplitude of 1% to 20%, and the higher threshold is the signal at an amplitude of 60% to 200%, and histogram equalization is performed according to the display threshold.

[0010] In one possible implementation, the elastic detection image includes at least one position mark, and the position mark is used to characterize the propagation position of the shear wave generated by the vibration excitation inside the tissue to be tested. According to the positional relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image, the ultrasonic detection image and the elastic detection image are superimposed and displayed, including: determining a positioning coordinate system based on the test area of ​​the ultrasonic detection image; determining the position coordinates of the position mark in the positioning coordinate system according to the positional relationship between the position mark and the test area; and superimposing and displaying the position marker on the ultrasonic detection image according to the position coordinates.

[0011] In a possible implementation, the elasticity detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color color, and the pseudo-color color is used to represent the time corresponding to when the shear wave arrives at different positions when it is transmitted inside the tissue to be tested.

[0012] In a possible implementation manner, each of the position marks has a different smoothness, and the smoothness is used to characterize the propagation stability of the shear wave.

[0013] In a possible implementation, after obtaining the ultrasonic echo signal, it further includes: calculating the ultrasonic attenuation parameter UAP according to the formula UAP = (I1 - I2) / (D1 - D2); where, I1 is the ultrasonic signal intensity at the first position of interest D1, D1 ≤ 5 cm; I2 is the ultrasonic signal intensity at the second position of interest D2, D2 ≤ 30 cm; D1 < D2; calculating the ultrasonic attenuation coefficient α according to the formula α = (20log(I1 / I2) - 6) / 2(D1 - D2); where, I1 is the ultrasonic signal intensity at the first position of interest D1, D1 ≤ 5 cm; I2 is the ultrasonic signal intensity at the second position of interest D2, D2 ≤ 30 cm; D1 < D2; calculating the scatterer distribution coefficient k according to the formula k = s / σ; where, s is the coherent signal energy; σ is the diffused signal energy; the range of k is [0, 1]; after normalizing the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient α, and the scatterer distribution coefficient k, according to the ultrasonic echo signal and the corresponding tissue pathological classification, according to the formula calculate the comprehensive score Score of liver tissue lesions; where, the value range of a is [0.5, 1], the value range of b is [0.2, 0.8], the value range of c is [0, 0.5], and the value range of d is [0.2, 1].

[0014] According to the second aspect of the embodiments of the present application, the present application provides an ultrasonic elastography device, including:

[0015] A detection module, configured to apply a vibration excitation to the tissue to be measured and perform ultrasonic detection to obtain an ultrasonic echo signal;

[0016] A synthesis module, configured to perform beam synthesis according to the ultrasonic echo signal to generate synthesized ultrasonic imaging data;

[0017] A generation module, configured to use the synthesized ultrasonic imaging data to generate an ultrasonic detection image and an elasticity detection image, where the elasticity detection image is used to characterize tissue elasticity detection information and the signal quality of the ultrasonic detection image at the corresponding moment;

[0018] A display module, configured to superimpose and display the ultrasonic detection image and the elasticity detection image according to the positional relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elasticity detection image.

[0019] In a possible implementation, the detection module is specifically configured to: apply a shear wave to the tissue to be measured and collect an ultrasonic echo signal for tracking the propagation of the shear wave; the synthesis module is specifically configured to: perform beam synthesis on the ultrasonic echo signal for tracking the propagation of the shear wave to generate multiple frames of scan lines, where the multiple frames of scan lines are used to characterize the ultrasonic detection results at different positions.

[0020] In a possible implementation, the generating module is used to: process the synthesized ultrasonic imaging data through a graphics processor, generate the ultrasonic detection image and the elastic detection image in parallel, and obtain the elastic modulus E according to the shear wave propagation velocity, where E = 3ρV 2 , ρ is tissue density, and V is shear wave velocity.

[0021] In a possible implementation, after generating the ultrasonic detection image, the generation module further processes the ultrasonic detection image as follows: performs median filtering to remove electronic noise; performs Gaussian smoothing and non-local mean filtering to smooth the structural information; performs bilateral filtering to enhance the edge information; uses two thresholds of the signal as display thresholds, the lower threshold is the signal amplitude of 1% to 20%, and the higher threshold is the signal amplitude of 60% to 200%, and performs histogram equalization processing according to the display thresholds.

[0022] In one possible implementation, the elastic detection image includes at least one position marker, which is used to characterize the propagation position of the shear wave generated by the vibration excitation inside the tissue to be tested. The display module is specifically used to: determine a positioning coordinate system based on the test area of ​​the ultrasound detection image; determine the position coordinates of the position marker in the positioning coordinate system based on the positional relationship between the position marker and the test area; and superimpose and display the position marker on the ultrasound detection image based on the position coordinates.

[0023] In a possible implementation, the elasticity detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color color, and the pseudo-color color is used to represent the time corresponding to when the shear wave arrives at different positions when it is transmitted inside the tissue to be tested.

[0024] In a possible implementation manner, each of the position marks has a different smoothness, and the smoothness is used to characterize the propagation stability of the shear wave.

[0025] In a possible implementation, after obtaining the ultrasonic echo signal, the generating module is further configured to: calculate the ultrasonic attenuation parameter UAP according to the formula UAP = (I1 - I2) / (D1 - D2); where, I1 is the ultrasonic signal intensity at the first position of interest D1, D1 ≤ 5 cm; I2 is the ultrasonic signal intensity at the second position of interest D2, D2 ≤ 30 cm; D1 < D2; calculate the ultrasonic attenuation coefficient α according to the formula α = (20log(I1 / I2) - 6) / 2(D1 - D2); where, I1 is the ultrasonic signal intensity at the first position of interest D1, D1 ≤ 5 cm; I2 is the ultrasonic signal intensity at the second position of interest D2, D2 ≤ 30 cm; D1 < D2; calculate the scatterer distribution coefficient k according to the formula k = s / σ; where, s is the coherent signal energy; σ is the diffused signal energy; the range of k is [0, 1]; after normalizing the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient α, and the scatterer distribution coefficient k, according to the ultrasonic echo signal and the corresponding tissue pathological classification, according to the formula calculate the comprehensive score Score of the liver tissue lesion; where, the value range of a is [0.5, 1], the value range of b is [0.2, 0.8], the value range of c is [0, 0.5], and the value range of d is [0.2, 1].

[0026] According to the third aspect of the embodiments of the present application, the present application provides an electronic device, including: a memory, a processor, and a computer program;

[0027] wherein, the computer program is stored in the memory and is configured to be executed by the processor to perform the ultrasonic elastography method according to any one of the first aspect of the embodiments of the present application.

[0028] According to the fourth aspect of the embodiments of the present application, the present application provides a computer-readable storage medium, in which computer execution instructions are stored, and when the computer execution instructions are executed by a processor, they are used to implement the ultrasonic elastography method according to any one of the first aspect of the embodiments of the present application.

[0029] The ultrasonic elastic imaging method, device, electronic device and storage medium provided in the present application obtain an ultrasonic echo signal by applying vibration excitation to the tissue to be tested and performing ultrasonic detection; beam synthesis is performed based on the ultrasonic echo signal to generate synthetic ultrasonic imaging data; ultrasonic detection images and elastic detection images are generated using the synthetic ultrasonic imaging data, wherein the elastic detection image is used to characterize tissue elasticity detection information and the signal quality of the ultrasonic detection image at the corresponding moment; the ultrasonic detection image and the elastic detection image are superimposed and displayed based on the positional relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image; after the ultrasonic detection image is generated, the generated elastic detection image is superimposed and displayed on the ultrasonic detection image, and the signal quality of the ultrasonic detection image is characterized by the elastic detection image, so that the operator can quickly control the detection results, improve the detection efficiency, and improve the detection accuracy. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0031] Figure 1 An application scenario diagram of the ultrasonic elastic imaging method provided in an embodiment of the present application;

[0032] Figure 2 A flowchart of an ultrasonic elastic imaging method provided in one embodiment of the present application;

[0033] Figure 3 A schematic diagram of a process for generating an ultrasonic detection image and an elastic detection image provided in an embodiment of the present application;

[0034] Figure 4 A schematic diagram of an elasticity detection image and an ultrasonic detection image superimposed and displayed in an embodiment of the present application;

[0035] Figure 5 A flowchart of an elastic imaging method provided in another embodiment of the present application;

[0036] Figure 6 A schematic diagram of a process for quality control based on the smoothness of position marks provided in an embodiment of the present application;

[0037] Figure 7 A schematic diagram of the structure of an ultrasonic elastic imaging device provided in one embodiment of the present application;

[0038] Figure 8 A schematic diagram of an electronic device provided for one embodiment of the present application.

[0039] The above drawings have shown clear embodiments of the present application, which will be described in more detail later. These drawings and text descriptions are not intended to limit the scope of the present application in any way, but to illustrate the concept of the present application to those skilled in the art by referring to specific embodiments. DETAILED DESCRIPTION

[0040] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the following exemplary embodiments do not represent all implementations consistent with the present application. Instead, they are merely examples of devices and methods consistent with some aspects of the present application as detailed in the appended claims.

[0041] The application scenarios of the embodiments of the present application are explained below:

[0042] Figure 1 An application scenario diagram of the ultrasonic elastic imaging method provided in the embodiment of the present application, such as Figure 1 As shown, the ultrasonic elastic imaging method provided in the embodiment of the present application is applied to electronic equipment, specifically, for example, an elastic imaging device 11 that can be used clinically. The elastic imaging device has a detection probe 111, which can apply vibration excitation and ultrasonic waves to the human body simultaneously or sequentially. The elastic imaging device 11 performs ultrasonic elastic detection on the human body through the elastic imaging method provided in the embodiment of the present application, obtains the mechanical characteristics of the structure of the human visceral tissue 12, and then determines whether the human visceral tissue 12 is in a healthy state, providing important data support for the screening of early lesions.

[0043] In the prior art, the results of ultrasonic elastic detection are usually displayed to the operator of the equipment in the form of quantitative parameter values. The operator observes the ultrasonic detection image and combines the quantitative parameter values ​​of elastic imaging to detect the state of biological tissue. However, due to the non-uniform characteristics of human tissue, during the elastic measurement process, poor test position, human breathing and other reasons will interfere with the ultrasonic detection image, affect the quality of ultrasonic imaging, and affect the operator's judgment of the state of the tissue. At the same time, the elastic detection results can characterize the propagation process of the shear wave. There is a difference between the elastic detection results in the disturbed state and the elastic detection results in the undisturbed state. Therefore, the elastic detection results can be used to control the quality of the ultrasonic detection image to improve the signal quality of the ultrasonic detection image and improve the accuracy of the operator's judgment of the state of the tissue. Therefore, how to achieve quality control of the ultrasonic detection image through the elastic detection results to improve the signal quality of the ultrasonic detection image is a problem that needs to be solved urgently.

[0044] The technical solution of the present application and how the technical solution of the present application solves the above-mentioned technical problems are described in detail below with specific embodiments. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments. The embodiments of the present application will be described below in conjunction with the accompanying drawings.

[0045] Figure 2 A flowchart of an ultrasonic elastic imaging method provided in one embodiment of the present application is as follows: Figure 2 As shown, the elastic imaging method provided in this embodiment is applied to an elastic imaging device, and includes the following steps:

[0046] Step S101, applying vibration excitation to the tissue to be tested and performing ultrasonic detection to obtain an ultrasonic echo signal.

[0047] For example, the operator applies vibration excitation to the tissue to be tested by operating the detection probe of the elastic imaging device, which can generate shear waves transmitted from the surface of the tissue to be tested to the inside of the tissue to be tested, and the shear waves will cause tiny displacements and deformations inside the tissue to be tested. At the same time, ultrasonic detection is applied to the tissue to be tested to obtain data on the tiny displacements and deformations of the tissue to be tested. Specifically, for example, the detection probe transmits ultrasonic waves to the tissue to be tested and receives echoes to obtain ultrasonic echo signals, and the ultrasonic echo signals can be used to observe the tiny displacements and deformations of the tissue to be tested.

[0048] Among them, the execution sequence of applying vibration excitation to the tissue to be tested and performing ultrasonic detection can include multiple types, such as applying vibration excitation first and then performing ultrasonic detection, or performing ultrasonic detection first and then applying vibration excitation, or applying vibration excitation and performing ultrasonic detection at the same time. No specific limitation is made here and it can be set as needed.

[0049] Furthermore, when applying vibration excitation to the tissue to be tested and performing ultrasonic detection, the vibration frequencies may be the same or different. In a possible implementation, vibration excitations of different frequencies are sequentially applied to the tissue to be tested and ultrasonic detection is performed. Specifically, for example, vibration excitation is cyclically applied to the tissue to be tested, and during each cycle, the vibration frequency of the vibration excitation is changed until a preset stop condition is met, and ultrasonic detection is performed on the tissue to be tested at the same time, wherein the timing of ultrasonic detection may be after each application of vibration excitation of a specific frequency, or after all vibration excitations are applied, which is not specifically limited this time. By applying vibration excitations of different frequencies to the tissue to be tested, shear waves of different frequencies are generated inside the tissue to be tested, and by applying ultrasonic detection to the tissue to be tested after each vibration excitation, the response detection of the tissue to be tested at different excitation frequencies can be achieved.

[0050] Step S102: performing beam synthesis according to the ultrasonic echo signal to generate synthetic ultrasonic imaging data.

[0051] Exemplarily, after transmitting ultrasound to the tissue to be tested and receiving the ultrasonic echo signal formed after the ultrasound is blocked by the tissue to be tested, the ultrasonic echo signal is beamformed, and the result of the beamformation is the ultrasonic imaging data. Exemplarily, the vibration excitation applied by the elastic imaging device to the tissue to be tested is a shear wave, and the ultrasonic echo signal is an ultrasonic echo signal collected by the elastic program device for tracking the transmission of the shear wave. After that, the ultrasonic echo signal for tracking the propagation of the shear wave is beamformed to generate synthesized ultrasonic imaging data. The implementation process of beamforming the ultrasonic signal involved in this process is a prior art in the art and will not be described in detail here.

[0052] Step S103, using the synthesized ultrasonic imaging data, generating an ultrasonic detection image and an elastic detection image, wherein the elastic detection image is used to characterize tissue elasticity detection information and the signal quality of the ultrasonic detection image at a corresponding moment.

[0053] Figure 3 A schematic diagram of a process for generating an ultrasonic detection image and an elastic detection image provided in an embodiment of the present application, referring to Figure 3, exemplarily, the synthetic ultrasonic imaging data is data carrying deformation information of the tissue to be tested, and the synthetic ultrasonic imaging data can be used for ultrasonic imaging and elastic imaging at the same time, that is, generating ultrasonic detection images and elastic detection images. Specifically, the step of applying ultrasonic detection to the tissue to be tested includes transmitting ultrasonic waves to the tissue to be tested and receiving echoes formed after the ultrasonic waves are blocked by the tissue to be tested. In each cycle of vibration excitation-ultrasonic detection, the shear wave formed by mechanical vibration excites the tissue to be tested, and the tissue to be tested undergoes a slight deformation. Therefore, the ultrasonic echo signal carries the deformation information of the tissue to be tested. The ultrasonic echo signal is subjected to beam synthesis and imaging calculation to obtain the corresponding ultrasonic detection image. The specific implementation method of generating the ultrasonic detection image according to the ultrasonic echo signal is the prior art in the art and will not be described here. Further, the corresponding elastic detection image can be obtained by processing and imaging calculation of the ultrasonic echo signal. Specifically, for example, the ultrasonic echo signal is first filtered, and then the tissue displacement or strain is calculated for the filtered ultrasonic echo signal, and a motion propagation mode diagram is obtained. Then, the dispersion curve is determined according to the motion propagation mode diagram, and the phase velocity is determined according to the dispersion curve. By fitting the phase velocity function, the viscoelastic information of the tissue to be tested is obtained. The method of fitting the phase velocity corresponding to different propagation rates to obtain the viscosity parameter is a prior art known to those skilled in the art and will not be repeated here. For example, the elastic detection image is an image representation of the viscoelastic information, and the elastic detection image can be realized in the form of a pseudo-color image. Since the generation of the ultrasonic detection image and the elastic detection image is generated by using the same source data, the detection area corresponding to the generation of the ultrasonic detection image and the elastic detection image is exactly the same position, and there is no time error. Compared with the scheme of generating the ultrasonic detection image and the elastic detection image by calculating the elastic imaging data calculated by the original ultrasonic echo signal and the original ultrasonic echo signal, the time error problem caused by the inconsistency between the original ultrasonic echo signal and the calculated elastic imaging data can be avoided, thereby further avoiding interference and improving the reliability of elastic measurement evaluation.

[0054] Optionally, the steps of this embodiment further include: obtaining the elastic modulus E according to the shear wave propagation velocity, where E = 3ρV 2 , ρ is tissue density, V is shear wave velocity. The shear wave propagation velocity is determined by the phase velocity, and the tissue density is obtained by the preset configuration information for recording the density of the tissue to be measured.

[0055] Optionally, in the steps of this embodiment, after obtaining the ultrasonic echo signal, the following steps are also included:

[0056] Calculate the ultrasonic attenuation parameter UAP according to the formula UAP = (I1 - I2) / (D1 - D2); where, I1 is the ultrasonic signal intensity at the first position of interest D1, D1 ≤ 5 cm; I2 is the ultrasonic signal intensity at the second position of interest D2, D2 ≤ 30 cm; D1 < D2; calculate the ultrasonic attenuation coefficient α according to the formula α = (20log(I1 / I2) - 6) / 2(D1 - D2); where, I1 is the ultrasonic signal intensity at the first position of interest D1, D1 ≤ 5 cm; I2 is the ultrasonic signal intensity at the second position of interest D2, D2 ≤ 30 cm; D1 < D2; calculate the scatterer distribution coefficient k according to the formula k = s / σ; where, s is the coherent signal energy; σ is the diffused signal energy; the range of k is [0, 1]; after normalizing the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient α, and the scatterer distribution coefficient k, according to the ultrasonic echo signal and the corresponding tissue pathological classification, according to the formula Calculate the comprehensive score Score of liver tissue lesions; where, the value range of a is [0.5, 1], the value range of b is [0.2, 0.8], the value range of c is [0, 0.5], and the value range of d is [0.2, 1]. Where, " " represents dot product operation.

[0057] Among them, the elastic modulus E evaluates the liver elasticity and the degree of liver fibrosis; the ultrasonic attenuation parameter UAP evaluates the degree of fatty change in liver tissue; the ultrasonic attenuation coefficient α evaluates the degree of inflammation in liver tissue; the scatterer distribution coefficient k evaluates the degree of bubble-like change in liver tissue; UAP / α / k comprehensively evaluates the degree of fatty liver in the liver; Score comprehensively evaluates the degree of liver tissue lesions; among them, after normalizing each parameter, calculate Score according to the ultrasonic signal and the corresponding tissue pathological classification.

[0058] In the steps of this embodiment, according to the ultrasonic echo signal, extract quantitative ultrasonic parameters and provide multi-dimensional evaluation information. In addition to the shown elastography images, the quantitative ultrasonic parameters are output and displayed simultaneously, assisting the user to better judge the changes in tissue structure through the optimized ultrasonic detection images, better judge the changes in tissue elasticity through the elastography images, better judge the changes in tissue microstructure through the quantitative ultrasonic parameters, comprehensively evaluate the tissue state, thereby improving the detection efficiency and detection accuracy.

[0059] Among them, exemplarily, a Graphics Processing Unit (GPU) is provided in the elastography device. The process of using the synthetic ultrasonic imaging data to generate ultrasonic detection images and elastography images includes: processing the synthetic ultrasonic imaging data through the GPU to generate ultrasonic detection images and elastography images in parallel, improving the image processing speed, enabling the detection images to be displayed in real time, and improving the detection accuracy.

[0060] Optionally, after the ultrasonic detection image is generated and before it is displayed, the ultrasonic detection image is further processed as follows:

[0061] Perform median filtering to remove electronic noise; perform Gaussian smoothing and non-local mean filtering to smooth the structural information; perform bilateral filtering to enhance the edge information; use the two thresholds of the signal as display thresholds, the lower threshold is the signal at 1% to 20% amplitude, and the higher threshold is the signal at 60% to 200% amplitude, and perform histogram equalization processing according to the display thresholds.

[0062] In the steps of this embodiment, the ultrasonic detection image is optimized before being displayed to make the processed ultrasonic detection image clearer, which helps the user to evaluate the state of the tissue to be tested by observing the ultrasonic detection image, thereby improving the detection efficiency and accuracy.

[0063] Step S104 , based on the positional relationship between the test area corresponding to the ultrasonic test image and the test area corresponding to the elastic test image, the ultrasonic test image and the elastic test image are superimposed and displayed.

[0064] For example, the ultrasonic detection image is image information representing the tissue morphology of the test area inside the tissue to be tested, including descriptive information about the spatial position inside the tissue to be tested. The elastic detection image is image information representing the propagation process of the shear wave in the test area. Figure 4 A schematic diagram of an elastic detection image and an ultrasonic detection image superimposed and displayed in an embodiment of the present application is shown as follows: Figure 4 As shown, the elastic detection image includes a plurality of position marks 31, and the position mark 31 is used to characterize the propagation position of the shear wave generated by vibration excitation inside the tissue to be tested. The position of the shear wave generated by vibration excitation inside the tissue to be tested has a real corresponding relationship with the spatial position of the tissue to be tested shown in the ultrasonic detection image 32. Each position mark has a different smoothness, and the smoothness is used to characterize the propagation stability of the shear wave. When the smoothness of the position mark 31 is good, it characterizes that the propagation stability of the shear wave is good, that is, at this time, the propagation of the shear wave is not affected by the breathing of the human body under test, the poor test position and other factors. Therefore, the signal quality of the ultrasonic detection image 32 is good and the credibility is high; on the contrary, when the smoothness of the position mark 31 is poor, it characterizes that the propagation stability of the shear wave is poor, that is, at this time, the propagation of the shear wave is affected by the breathing of the human body under test, the poor test position and other factors. Therefore, the signal quality of the ultrasonic detection image 32 is poor and the credibility is low.

[0065] Further, exemplarily, the elastic imaging device includes a display unit, such as a display screen. When the operator performs elasticity detection on the tissue to be tested, it is necessary to continuously adjust the position of the detection probe to locate a better test area so that the internal tissue state of the tissue to be tested can be better represented in the ultrasound image. Through this embodiment, the ultrasound detection image and the elasticity detection image are superimposed and displayed on the display screen. At this time, since the elasticity detection image plays a quality control role on the ultrasound detection image, the operator can use the quality control role of the ultrasound detection image on the elasticity detection image in the process of continuously adjusting the position of the detection probe and locating the better test area, so as to better determine whether the current position of the detection probe corresponds to the better test area, thereby improving the efficiency of locating the test area, improving the signal quality of the ultrasound detection image, and improving the accuracy of the detection results.

[0066] An ultrasonic echo signal is obtained by applying vibration excitation to the tissue to be tested and performing ultrasonic detection; an ultrasonic detection image and an elastic detection image are generated based on the ultrasonic echo signal, wherein the elastic detection image is used to characterize the tissue elastic detection information and the signal quality of the ultrasonic detection image at the corresponding moment; the ultrasonic detection image and the elastic detection image are superimposed and displayed based on the positional relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image; after the ultrasonic detection image is generated, the generated elastic detection image is simultaneously used and superimposed on the ultrasonic detection image, and the elastic detection image is used to characterize the signal quality of the ultrasonic detection image, so that the operator can quickly control the test results, eliminate the influence caused by poor test position, breathing of the tested person, etc., improve the signal quality of the ultrasonic detection image, and improve the accuracy of the test results.

[0067] Figure 5 A flowchart of an elastic imaging method provided in another embodiment of the present application is as follows: Figure 5 As shown, the elastic imaging method provided in this embodiment is Figure 2 On the basis of the elastic imaging method provided in the illustrated embodiment, step S104 is further refined, and a step of obtaining configuration information is added. The elastic imaging method provided in this embodiment includes the following steps:

[0068] Step S201 : applying shear waves to the tissue to be tested, and collecting ultrasonic echo signals for tracking the propagation of the shear waves.

[0069] Step S202 , beamforming the ultrasonic echo signal used to track the propagation of shear waves to generate synthesized ultrasonic imaging data, where the ultrasonic imaging data includes multiple frames of scan lines, and the multiple frames of scan lines are used to represent ultrasonic detection results at different positions.

[0070] Step S203, generating an ultrasonic detection image and an elastic detection image according to multiple frames of scan lines in the synthesized ultrasonic imaging data, wherein the elastic detection image includes at least one position mark, and the position mark is used to characterize the propagation position of the shear wave generated by vibration excitation inside the tissue to be tested.

[0071] Exemplarily, the ultrasonic detection image is image information generated by beam synthesis and imaging calculation of multiple frames of scanning lines. The elastic detection image is image information that characterizes the viscoelastic information of the tissue to be tested. The process of determining the viscoelastic information of the tissue to be tested includes, for example: filtering the multiple frames of scanning lines, determining the motion propagation mode diagram according to the echo data corresponding to the filtered multiple frames of scanning lines, performing time-frequency analysis on the motion propagation mode diagram, determining the dispersion curve, determining the phase velocity corresponding to different frequencies according to the dispersion curve, performing function fitting on the phase velocity corresponding to different frequencies, and obtaining the viscoelastic information of the tissue to be tested. The ultrasonic detection image and the elastic detection image generated by the multiple frames of scanning lines have the same source data, so the detection area corresponding to the generated ultrasonic detection image and the elastic detection image is exactly the same position, and there is no time error. It can avoid the time error problem caused by the inconsistency between the original ultrasonic echo signal and the calculated elastic imaging data, thereby further avoiding interference and improving the reliability of elastic measurement evaluation.

[0072] Step S204: determine a positioning coordinate system based on the test area of ​​the ultrasonic detection image.

[0073] Furthermore, a set of synthetic ultrasonic imaging data corresponds to a set of ultrasonic echo information, and the sum of the areas detected by each ultrasonic echo signal is the test area. Therefore, based on the test area corresponding to the synthetic ultrasonic imaging data, a coordinate system for describing the spatial position can be established, and the position coordinates in the coordinate system are used to characterize the spatial position of the coordinate point. Furthermore, the position mark in the elastic detection image is also generated by the ultrasonic echo signal, so the position coordinates of the position mark in the coordinate system based on the test area can be determined based on the synthetic ultrasonic imaging data.

[0074] Step S205: determining the position coordinates of the position marker in the positioning coordinate system according to the position relationship between the position marker and the test area.

[0075] In a possible implementation, the range of the test area can be adjusted and set according to user needs, for example, the configuration information input by the user is obtained, and the range of the test area is determined according to the configuration information. By adjusting the range of the test area through the configuration information, only the area of ​​interest or the preset range is displayed to avoid the influence of the surrounding non-interest area or other interference information on the judgment of the result, which can further improve the flexibility and scope of application of the algorithm.

[0076] Step S206: superimpose and display the position mark on the ultrasonic detection image according to the position coordinates.

[0077] Specifically, according to the position coordinates, the position marker is displayed at the position coordinates to realize the superposition display of the position marker and the ultrasonic detection image. The position marker represents the spatial position where the shear wave arrives over time. Optionally, the position marker is a pseudo-color image, and different position markers have different colors, and different colors represent the time when the shear wave arrives at the position. For example, position marker A is red, indicating that the time when the shear wave arrives at position marker A is 50 milliseconds, position marker B is yellow, indicating that the time when the shear wave arrives at position marker B is 100 milliseconds; position marker C is blue, indicating that the time when the shear wave arrives at position marker C is 150 milliseconds. According to the color change between the position markers, the position marker has a corresponding smoothness, wherein the smoothness is used to characterize the propagation stability of the shear wave. Exemplarily, the color change between each position marker can be changed according to the color temperature, which is more convenient for the operator to identify the smoothness of the color change between each position marker. The following is a more specific embodiment to illustrate the process of how to perform quality control based on the smoothness of the position marker.

[0078] Figure 6 A schematic diagram of a process for quality control based on the smoothness of position marks provided in an embodiment of the present application, such as Figure 6 As shown, position markers A, B, C, D, and E are superimposed and displayed in the ultrasonic detection image 32, and each position marker 31 represents the position where the shear wave arrives at the corresponding human tissue in the ultrasonic detection image. When the operator performs elasticity detection by operating the detection probe 111, the position marker 31 is superimposed and displayed on the ultrasonic detection image 32 in different colors, that is, the position markers A, B, C, D, and E have different pseudo-color colors, for example, position marker A is red, position marker B is orange-red, position marker C is yellow, position marker D is light blue, and position marker E is dark blue. Different colors correspond to different times when the shear wave arrives at the position.

[0079] When the operator moves the detection probe 111 to the first position, the detection environment is not good at this time, and the shear wave is affected during the propagation process in the human tissue, making the propagation speed of the shear wave in the human body uneven, which will cause the color change from the marked position A to the marked position E to be uneven. The operator observes the color change of the position mark 31 and determines that the detection environment is not good at this time, so he moves the detection probe 111 to the second position. At this time, the detection environment is good, that is, the shear wave is not affected by human breathing, uneven tissue obstruction, etc. during the propagation process in the human tissue, then the propagation speed of the shear wave should be uniform, that is, the color change from the position mark A to the position mark E is smooth. By observing the smoothness of the pseudo-color image-type position mark 31, the operator can determine whether there is interference in the test area, thereby achieving the purpose of quality control of the ultrasonic detection image.

[0080] In this embodiment, the marked position is displayed in the form of a pseudo-color image, and the signal quality of the ultrasonic detection image is represented by the smoothness of each position mark. The elastic detection image superimposed on the ultrasonic detection image can improve the quality control of the ultrasonic detection image. Compared with the technical solution of displaying the elasticity detection result by text mark in the prior art, the display method of the pseudo-color image can enable the operator to observe the signal quality of the ultrasonic detection image more intuitively and in real time, thereby adjusting the test area, avoiding interference with the elasticity detection, and improving the accuracy of the detection of the state of the tissue to be tested.

[0081] In this embodiment, the implementation of step S201-step S202 is described in this application. Figure 2 The above has been introduced in the illustrated embodiment and will not be described in detail here.

[0082] Figure 7 A schematic diagram of the structure of an ultrasonic elastic imaging device provided in an embodiment of the present application is applied to elastic imaging equipment, such as Figure 7 As shown, the ultrasonic elastic imaging device 3 provided in this embodiment includes:

[0083] The detection module 31 is used to apply vibration excitation to the tissue to be detected and perform ultrasonic detection to obtain an ultrasonic echo signal;

[0084] A synthesis module 32, used for performing beam synthesis according to the ultrasonic echo signal to generate synthetic ultrasonic imaging data;

[0085] A generating module 33, used to generate an ultrasonic detection image and an elastic detection image using the synthesized ultrasonic imaging data, wherein the elastic detection image is used to characterize tissue elasticity detection information and the signal quality of the ultrasonic detection image at a corresponding moment;

[0086] The display module 34 is used to overlay and display the ultrasonic detection image and the elastic detection image according to the positional relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image.

[0087] In one possible implementation, the detection module 31 is specifically used to apply shear waves to the tissue to be tested and collect ultrasonic echo signals for tracking the propagation of shear waves; the synthesis module 32 is specifically used to perform beam synthesis on the ultrasonic echo signals for tracking the propagation of shear waves to generate multiple frames of scan lines, wherein the multiple frames of scan lines are used to characterize the ultrasonic detection results at different positions.

[0088] In a possible implementation, the generating module 33 is used to: process the synthesized ultrasonic imaging data through a graphics processor, generate an ultrasonic detection image and an elastic detection image in parallel, and obtain an elastic modulus E according to the shear wave propagation velocity, where E = 3ρV 2 , ρ is tissue density, and V is shear wave velocity.

[0089] In a possible implementation, after generating the ultrasonic detection image, the generation module 33 further processes the ultrasonic detection image as follows: performs median filtering to remove electronic noise; performs Gaussian smoothing and non-local mean filtering to smooth the structural information; performs bilateral filtering to enhance the edge information; uses two thresholds of the signal as display thresholds, the lower threshold is the signal amplitude of 1% to 20%, and the higher threshold is the signal amplitude of 60% to 200%, and performs histogram equalization processing according to the display threshold.

[0090] In one possible implementation, the elastic detection image includes at least one position marker, which is used to characterize the propagation position of the shear wave generated by vibration excitation inside the tissue to be tested. The display module 34 is specifically used to: determine the positioning coordinate system based on the test area of ​​the ultrasonic detection image; determine the position coordinates of the position marker in the positioning coordinate system based on the positional relationship between the position marker and the test area; and superimpose the position marker on the ultrasonic detection image based on the position coordinates.

[0091] In a possible implementation, the elasticity detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color color, and the pseudo-color color is used to represent the time corresponding to when the shear wave arrives at different positions when it is transmitted inside the tissue to be tested.

[0092] In a possible implementation, each position mark has a different smoothness, and the smoothness is used to characterize the propagation stability of the shear wave.

[0093] In a possible implementation, after obtaining the ultrasonic echo signal, the generating module 33 is further configured to: calculate the ultrasonic attenuation parameter UAP according to the formula UAP = (I1 - I2) / (D1 - D2); where I1 is the ultrasonic signal intensity at the first position of interest D1, D1 ≤ 5 cm; I2 is the ultrasonic signal intensity at the second position of interest D2, D2 ≤ 30 cm; D1 < D2; calculate the ultrasonic attenuation coefficient α according to the formula α = (20log(I1 / I2) - 6) / 2(D1 - D2); where I1 is the ultrasonic signal intensity at the first position of interest D1, D1 ≤ 5 cm; I2 is the ultrasonic signal intensity at the second position of interest D2, D2 ≤ 30 cm; D1 < D2; calculate the scatterer distribution coefficient k according to the formula k = s / σ; where s is the coherent signal energy; σ is the diffused signal energy; the range of k is [0, 1]; after normalizing the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient α, and the scatterer distribution coefficient k, according to the ultrasonic echo signal and the corresponding tissue pathological classification, according to the formula calculate the comprehensive score Score of the liver tissue lesion; where the value range of a is [0.5, 1], the value range of b is [0.2, 0.8], the value range of c is [0, 0.5], and the value range of d is [0.2, 1].

[0094] Among them, the detection module 31, the synthesis module 32, the generation module 33, and the display module 34 are connected in sequence. The ultrasonic elastography device 3 provided in this embodiment can execute the technical solution of the method embodiment as described in Figure 2-Figure 6 any one of them. The implementation principle and technical effects are similar and will not be elaborated here.

[0095] Figure 8 It is a schematic diagram of an electronic device provided in an embodiment of the present application. As shown in Figure 8 shown, the electronic device 4 provided in this embodiment includes: a memory 41, a processor 42, and a computer program.

[0096] Among them, the computer program is stored in the memory 41 and is configured to be executed by the processor 42 to implement the ultrasonic elastography method provided in any one of the embodiments corresponding to the present application Figure 2-Figure 5 corresponding to the embodiments.

[0097] Among them, the memory 41 and the processor 42 are connected through a bus 43.

[0098] For relevant descriptions, reference can be made to the relevant descriptions and effects corresponding to the steps in the corresponding embodiments Figure 2-Figure 6 corresponding to the embodiments, and no further elaboration will be made here.

[0099] In a possible implementation, the electronic device is an elastography device.

[0100] An embodiment of the present application provides a computer-readable storage medium having a computer program stored thereon, and the computer program is executed by a processor to implement the present application. Figure 2-Figure 6 The ultrasonic elastic imaging method provided in any one of the corresponding embodiments.

[0101] Among them, the computer readable storage medium may be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, and the like.

[0102] One embodiment of the present application provides a computer program product, including a computer program, which implements the present application when executed by a processor. Figure 2-Figure 6 The ultrasonic elastic imaging method provided in any one of the corresponding embodiments.

[0103] In the several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of modules is only a logical function division. There may be other division methods in actual implementation, such as multiple modules or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or modules, which can be electrical, mechanical or other forms.

[0104] Those skilled in the art will readily appreciate other embodiments of the present application after considering the specification and practicing the application disclosed herein. The present application is intended to cover any modification, use or adaptation of the present application, which follows the general principles of the present application and includes common knowledge or customary techniques in the art that are not disclosed in the present application. The specification and examples are intended to be exemplary only, and the true scope and spirit of the present application are indicated by the following claims.

[0105] It should be understood that the present application is not limited to the precise structures that have been described above and shown in the drawings, and that various modifications and changes may be made without departing from the scope thereof. The scope of the present application is limited only by the appended claims.

Claims

1. An ultrasound elastic imaging system, the system is used to perform the following steps by a full computer: Obtaining an ultrasonic echo signal that applies vibration excitation to the tissue to be tested; Performing beam synthesis according to the ultrasonic echo signal to generate synthetic ultrasonic imaging data; Using the synthesized ultrasonic imaging data, an ultrasonic detection image and an elastic detection image are generated, wherein: The elasticity detection image is used to characterize tissue elasticity detection information and the signal quality of the ultrasonic detection signal at the corresponding moment; According to the positional relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image, the ultrasonic detection image and the elastic detection image are superimposed and displayed; After obtaining the ultrasonic echo signal, the system is further used for: The ultrasonic attenuation parameter UAP is calculated according to the formula UAP=(I1-I2) / (D1-D2); The ultrasonic attenuation coefficient α is calculated according to the formula α=(20log(I1 / I2)-6) / 2(D1-D2); where I1 is the ultrasonic signal intensity at the first location of interest D1, D1≤5cm; I2 is the ultrasonic signal intensity at the second location of interest D2, D2≤30cm; D1 <D2; The scatterer distribution coefficient k is calculated according to the formula k=s / σ; where s is the coherent signal energy; σ is the diffuse signal energy; the range of k is [0, 1]; After normalizing the elastic modulus E, ultrasonic attenuation parameter UAP, ultrasonic attenuation coefficient α, and scatterer distribution coefficient k, according to the ultrasonic echo signal and the corresponding tissue pathological typing, the formula Calculate the comprehensive score of liver tissue lesions Score; where the value range of a is [0.5, 1], the value range of b is [0.2, 0.8], the value range of c is [0, 0.5], and the value range of d is [0.2, 1].

2. The system according to claim 1, characterized in that The system performs beam synthesis according to the ultrasonic echo signal to generate synthetic ultrasonic imaging data, including: The ultrasonic echo signal used to track the propagation of shear waves is beam synthesized to generate multiple frames of scan lines, wherein the multiple frames of scan lines are used to represent ultrasonic detection results at different positions.

3. The system according to claim 1, characterized in that The system generates an ultrasound detection image and an elasticity detection image using the synthesized ultrasound imaging data, including: The synthetic ultrasonic imaging data is processed by a graphics processor to generate the ultrasonic detection image and the elastic detection image in parallel, and the elastic modulus E is obtained according to the shear wave propagation velocity, where E = 3ρV 2 , ρ is tissue density, and V is shear wave velocity.

4. The system according to claim 3, characterized in that After the ultrasonic detection image is generated but before it is displayed, the system is further used to process the ultrasonic detection image as follows: Perform median filtering to remove electronic noise; Gaussian smoothing and non-local mean filtering are performed to smooth the structural information; Perform bilateral filtering to enhance edge information; The two thresholds of the signal are used as display thresholds, the lower threshold is the signal at 1% to 20% amplitude, and the higher threshold is the signal at 60% to 200% amplitude, and histogram equalization processing is performed according to the display thresholds.

5. The system according to claim 1, characterized in that The elastic detection image includes at least one position mark, and the position mark is used to characterize the propagation position of the shear wave generated by the vibration excitation inside the tissue to be tested. The system superimposes and displays the ultrasonic detection image and the elastic detection image according to the positional relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image, including: Determine a positioning coordinate system based on the test area of ​​the ultrasonic detection image; Determining the position coordinates of the position mark in the positioning coordinate system according to the position relationship between the position mark and the test area; The position mark is superimposed and displayed on the ultrasonic detection image according to the position coordinates.

6. The system according to claim 5, characterized in that The elasticity detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color color, and the pseudo-color color is used to represent the time corresponding to when the shear wave arrives at different positions when it is transmitted inside the tissue to be tested.

7. The system according to claim 6, characterized in that Each of the position marks has a different smoothness, and the smoothness is used to characterize the propagation stability of the shear wave.

8. An ultrasonic elastic imaging device, characterized in that: The device comprises: A detection module is used to apply vibration excitation to the tissue to be tested and perform ultrasonic detection to obtain an ultrasonic echo signal; A synthesis module, used for performing beam synthesis according to the ultrasonic echo signal to generate synthesized ultrasonic imaging data; A generating module, used to generate an ultrasonic detection image and an elastic detection image using the synthesized ultrasonic imaging data, wherein the elastic detection image is used to characterize tissue elasticity detection information and the signal quality of the ultrasonic detection image at a corresponding moment; A display module, configured to display the ultrasonic detection image and the elastic detection image in a superimposed manner according to a positional relationship between a test area corresponding to the ultrasonic detection image and a test area corresponding to the elastic detection image; After obtaining the ultrasonic echo signal, the generating module is further used to: The ultrasonic attenuation parameter UAP is calculated according to the formula UAP=(I1-I2) / (D1-D2); The ultrasonic attenuation coefficient α is calculated according to the formula α=(20log(I1 / I2)-6) / 2(D1-D2); where I1 is the ultrasonic signal intensity at the first location of interest D1, D1≤5cm; I2 is the ultrasonic signal intensity at the second location of interest D2, D2≤30cm; D1 <D2; The scatterer distribution coefficient k is calculated according to the formula k = s / σ, where s is the coherent signal energy; σ is the diffuse signal energy; the range of k is [0, 1]; After normalizing the elastic modulus E, ultrasonic attenuation parameter UAP, ultrasonic attenuation coefficient α, and scatterer distribution coefficient k, according to the ultrasonic echo signal and the corresponding tissue pathological typing, the formula Calculate the comprehensive score of liver tissue lesions Score; where the value range of a is [0.5, 1], the value range of b is [0.2, 0.8], the value range of c is [0, 0.5], and the value range of d is [0.2, 1]; The display module is also used to display the obtained liver tissue lesion comprehensive score Score.

9. The device according to claim 8, characterized in that The detection module is specifically used for: Applying shear waves to the tissue to be tested and collecting ultrasonic echo signals for tracking the propagation of the shear waves; The synthesis module is specifically used to: perform beam synthesis on the ultrasonic echo signal used to track the propagation of shear waves to generate multiple frames of scan lines, wherein the multiple frames of scan lines are used to represent the ultrasonic detection results at different positions.

10. The device according to claim 8, characterized in that The generating module is used to: process the synthesized ultrasonic imaging data through a graphics processor, generate the ultrasonic detection image and the elastic detection image in parallel, and obtain the elastic modulus E according to the shear wave propagation velocity, wherein E = 3ρV 2 , ρ is tissue density, and V is shear wave velocity.

11. The device according to claim 10, characterized in that After generating the ultrasonic detection image, the generating module further processes the ultrasonic detection image as follows: Perform median filtering to remove electronic noise; Gaussian smoothing and non-local mean filtering are performed to smooth the structural information; Perform bilateral filtering to enhance edge information; The two thresholds of the signal are used as display thresholds, the lower threshold is the signal at 1% to 20% amplitude, and the higher threshold is the signal at 60% to 200% amplitude, and histogram equalization processing is performed according to the display thresholds.

12. The device according to claim 8, characterized in that The elasticity detection image includes at least one position mark, and the position mark is used to characterize the propagation position of the shear wave generated by the vibration excitation inside the tissue to be tested. The display module is specifically used to: Determine a positioning coordinate system based on the test area of ​​the ultrasonic detection image; Determining the position coordinates of the position mark in the positioning coordinate system according to the position relationship between the position mark and the test area; The position mark is superimposed and displayed on the ultrasonic detection image according to the position coordinates.

13. The device according to claim 12, characterized in that The elasticity detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color color, and the pseudo-color color is used to represent the time corresponding to when the shear wave arrives at different positions when it is transmitted inside the tissue to be tested.

14. The device according to claim 13, characterized in that Each of the position marks has a different smoothness, and the smoothness is used to characterize the propagation stability of the shear wave.

15. An electronic device, characterized in that: include: Memory, processor and computer program; The computer program is stored in the memory and is configured to be executed by the processor to implement the ultrasonic elastic imaging method according to any one of claims 1 to 7.

16. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions, and when the computer-executable instructions are executed by a processor, they are used to implement the ultrasonic elasticity imaging method according to any one of claims 1 to 7. 17 . A computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the method for ultrasonic elastography according to claim 1 is implemented.

Images