Shear wave image generation method and device and ultrasonic equipment

Abstract

The embodiment of the invention provides a shear wave image generation method and device and ultrasonic equipment, and relates to the technical field of shear wave elastic imaging. The ultrasonic device comprises: at least one probe for transmitting an ultrasonic signal; the processor is used for detecting an echo signal of the ultrasonic signal; determining the background echo intensity of each detection position point of the target detection part and the displacement peak value of each detection position point based on the detection result of the echo signal; based on the displacement peak value of each detection position point, determining the shear wave speed of the shear wave when the shear wave is propagated at each detection position point; for any detection position point, determining display parameters of the detection position point according to the background echo intensity of the detection position point and the displacement peak value of the detection position point; generating a shear wave image based on the shear wave speed corresponding to each detection position point and the display parameter of each detection position point; and the display screen is used for displaying the shear wave image.

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, device and ultrasonic equipment. Background Art

[0002] Shear wave elastography can realize real-time quantitative hardness detection of biological tissues, providing a basis for judging the pathological conditions of biological tissues. The basic principle of this technology is as follows: focusing the impact energy through acoustic radiation force can generate shear waves inside the tissue. Since the propagation speed of shear waves in tissues of different hardness is different, the propagation speed of shear waves at different locations can indirectly reflect the hardness and softness of the location. The image obtained by pseudo-color mapping according to the shear wave velocity is shear wave elastography.

[0003] However, in practical applications, due to the complex structure of human tissue, shear waves will split and refract when passing through different tissue interfaces, but these situations cannot be distinguished only by the shear wave velocity. Therefore, how to provide a shear wave image generation method to make the shear wave image include more abundant shear wave characteristic information is urgently needed. Summary of the invention

[0004] The embodiments of the present application provide a shear wave image generation method, an apparatus and an ultrasonic device, so as to provide a shear wave image generation method so that the shear wave image includes more abundant shear wave characteristic information.

[0005] In the first aspect, an embodiment of the present application provides an ultrasound device, comprising a display screen, a processor and at least one probe; the at least one probe is used to emit an ultrasound signal during an ultrasound inspection of a target detection part of an object; the processor is used to detect the echo signal of the ultrasound signal; based on the detection result of the echo signal, the background echo intensity of each detection point of the target detection part and the displacement peak of each detection point are determined; based on the displacement peak of each detection point, the shear wave velocity when the shear wave propagates at each detection point is determined; for any detection point, the display parameters of the detection point are determined according to the background echo intensity of the detection point and the displacement peak of the detection point, the display parameters are used to characterize the strength of the vibration amplitude caused at the detection point when the shear wave propagates; based on the shear wave velocity corresponding to each detection point and the display parameters of each detection point, a shear wave image is generated; the display screen is used to display the shear wave image.

[0006] In one possible implementation, when the processor determines, for any detection location point, the display parameters of the detection location point based on the background echo intensity of the detection location point and the displacement peak value of the detection location point, it is used to: for any detection location point, use the background echo intensity of the detection location point as the ambient light impact parameter in the Lambert illumination model, and correct the normal vector height corresponding to the detection location point according to the displacement peak value; determine the diffuse reflection value of the detection location point based on the Lambert illumination model and the corrected normal vector height corresponding to the detection location point; determine the display parameters of the detection location point according to the diffuse reflection value of the detection location point and the mapping relationship between the diffuse reflection value and the display parameters.

[0007] In one possible implementation, the processor determines the displacement peak value of each detection position point according to the following method: for any detection position point, the following operations are performed respectively: based on the detection result corresponding to the detection position point in the detection result of the echo signal, the axial displacement of the detection position point within a preset time length corresponding to each moment is determined by a two-dimensional autocorrelation algorithm; the maximum value of each axial displacement of the detection position point is determined by a maximum value search method, and the maximum value is used as the displacement peak value of the detection position point.

[0008] In one possible implementation, the processor determines the shear wave velocity when the shear wave propagates at each detection position point based on the displacement peak value of each detection position point, and is used to: for any detection position point, perform the following operations respectively: determine the N detection position points included in the sliding window corresponding to the detection position point; N is a positive integer; fit the displacement peak value of each detection position point in the N detection position points and the acquisition time corresponding to the displacement peak value through the least squares method to obtain the average shear wave velocity of the N detection position points; and use the average shear wave velocity as the shear wave velocity of the detection position point.

[0009] In one possible implementation, the detection result of the echo signal includes the in-phase data and the orthogonal data of each detection position point of the target detection part; when the processor determines the background echo intensity of each detection position point of the target detection part based on the detection result of the echo signal, it is used to: for any detection position point of the target detection part, perform a modulo operation on the in-phase data and the orthogonal data of the detection position point to determine the background echo intensity of the detection position point.

[0010] In one possible implementation, when the processor generates a shear wave image based on the shear wave velocities corresponding to each detection position point and the display parameters of each detection position point, it is used to: determine the color value corresponding to each detection position point according to a preset mapping relationship between shear wave velocity and color value and the shear wave velocities corresponding to each detection position point; and generate a shear wave image based on the color value corresponding to each detection position point and the display parameters of each detection position point.

[0011] In a second aspect, an embodiment of the present application provides a shear wave image generation method, comprising:

[0012] During ultrasonic inspection of a target detection part of an object, an ultrasonic signal is emitted and an echo signal of the ultrasonic signal is detected; based on the detection result of the echo signal, the background echo intensity of each detection point of the target detection part and the displacement peak of each detection point are determined; based on the displacement peak of each detection point, the shear wave velocity when the shear wave propagates at each detection point is determined; for any detection point, the display parameters of the detection point are determined according to the background echo intensity of the detection point and the displacement peak of the detection point, and the display parameters are used to characterize the strength of the vibration amplitude caused at the detection point when the shear wave propagates; based on the shear wave velocity corresponding to each detection point and the display parameters of each detection point, a shear wave image is generated.

[0013] In a possible implementation, for any detection position point, the display parameters of the detection position point are determined according to the background echo intensity of the detection position point and the displacement peak value of the detection position point, including: for any detection position point, the background echo intensity of the detection position point is used as the ambient light impact parameter in the Lambert illumination model, and the normal vector height corresponding to the detection position point is corrected according to the displacement peak value; based on the Lambert illumination model and the corrected normal vector height corresponding to the detection position point, the diffuse reflection value of the detection position point is determined; according to the diffuse reflection value of the detection position point, and the mapping relationship between the diffuse reflection value and the display parameter, the display parameters of the detection position point are determined.

[0014] In one possible implementation, the displacement peak value of each detection position point is determined according to the following method: for any detection position point, the following operations are performed respectively: based on the detection result corresponding to the detection position point in the detection result of the echo signal, the axial displacement of the detection position point within the preset time length corresponding to each moment is determined by a two-dimensional autocorrelation algorithm; the maximum value of each axial displacement of the detection position point is determined by a maximum value search method, and the maximum value is used as the displacement peak value of the detection position point.

[0015] In one possible implementation, based on the displacement peak value of each detection position point, the shear wave velocity when the shear wave propagates at each detection position point is determined, including: for any detection position point, the following operations are performed respectively: N detection position points included in the sliding window corresponding to the detection position point are determined; N is a positive integer; the displacement peak value of each detection position point in the N detection position points and the acquisition time corresponding to the displacement peak value are fitted by the least squares method to obtain the average shear wave velocity of the N detection position points; and the average shear wave velocity is used as the shear wave velocity of the detection position point.

[0016] In one possible implementation, the detection result of the echo signal includes the in-phase data and the orthogonal data of each detection position point of the target detection part; and determining the background echo intensity of each detection position point of the target detection part based on the detection result of the echo signal includes: for any detection position point of the target detection part, performing a modulo operation on the in-phase data and the orthogonal data of the detection position point to determine the background echo intensity of the detection position point.

[0017] In one possible implementation, the shear wave image is generated based on the shear wave velocities corresponding to each detection position point and the display parameters of each detection position point, including: determining the color value corresponding to each detection position point according to a preset mapping relationship between shear wave velocity and color value, and the shear wave velocities corresponding to each detection position point; and generating the shear wave image based on the color value corresponding to each detection position point and the display parameters of each detection position point.

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

[0019] A transmitting unit, for transmitting an ultrasonic signal during an ultrasonic inspection of a target detection part of an object;

[0020] A receiving unit, configured to receive an echo signal of the ultrasonic signal;

[0021] A processing unit, used to perform the following processing: detecting the echo signal of the ultrasonic signal; determining the background echo intensity of each detection position point of the target detection part and the displacement peak value of each detection position point based on the detection result of the echo signal; determining the shear wave velocity when the shear wave propagates at each detection position point based on the displacement peak value of each detection position point; for any detection position point, determining the display parameters of the detection position point according to the background echo intensity of the detection position point and the displacement peak value of the detection position point, the display parameters being used to characterize the strength of the vibration amplitude caused at the detection position point when the shear wave propagates; generating a shear wave image based on the shear wave velocity corresponding to each detection position point and the display parameters of each detection position point;

[0022] A display unit is used to display the shear wave image.

[0023] In one possible implementation, when the processing unit determines, for any detection position point, the display parameters of the detection position point based on the background echo intensity of the detection position point and the displacement peak value of the detection position point, it is used to: for any detection position point, use the background echo intensity of the detection position point as the ambient light impact parameter in the Lambert illumination model, and correct the normal vector height corresponding to the detection position point according to the displacement peak value; determine the diffuse reflection value of the detection position point based on the Lambert illumination model and the corrected normal vector height corresponding to the detection position point; determine the display parameters of the detection position point according to the diffuse reflection value of the detection position point and the mapping relationship between the diffuse reflection value and the display parameters.

[0024] In one possible implementation, the processing unit determines the displacement peak value of each detection position point according to the following method: for any detection position point, the following operations are performed respectively: based on the detection result corresponding to the detection position point in the detection result of the echo signal, the axial displacement of the detection position point within a preset time length corresponding to each moment is determined by a two-dimensional autocorrelation algorithm; the maximum value of each axial displacement of the detection position point is determined by a maximum value search method, and the maximum value is used as the displacement peak value of the detection position point.

[0025] In one possible implementation, when the processing unit determines the shear wave velocity when the shear wave propagates at each detection position point based on the displacement peak value of each detection position point, it is used to: for any detection position point, perform the following operations respectively: determine the N detection position points included in the sliding window corresponding to the detection position point; N is a positive integer; fit the displacement peak value of each detection position point in the N detection position points and the acquisition time corresponding to the displacement peak value through the least squares method to obtain the average shear wave velocity of the N detection position points; and use the average shear wave velocity as the shear wave velocity of the detection position point.

[0026] In one possible implementation, the detection result of the echo signal includes the in-phase data and the orthogonal data of each detection position point of the target detection part; when the processing unit determines the background echo intensity of each detection position point of the target detection part based on the detection result of the echo signal, it is used to: for any detection position point of the target detection part, perform a modulo operation on the in-phase data and the orthogonal data of the detection position point to determine the background echo intensity of the detection position point.

[0027] In one possible implementation, when the processing unit generates a shear wave image based on the shear wave velocities corresponding to each detection position point and the display parameters of each detection position point, it is used to: determine the color value corresponding to each detection position point according to a preset mapping relationship between shear wave velocity and color value and the shear wave velocities corresponding to each detection position point; and generate a shear wave image based on the color value corresponding to each detection position point and the display parameters of each detection position point.

[0028] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium, including:

[0029] The computer-readable storage medium stores computer instructions, and when the computer instructions are executed on a computer, the computer is enabled to perform the method as described in any one of the first aspects.

[0030] 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 above-mentioned first aspects is implemented.

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

[0032] The embodiments of the present application provide a shear wave image generation method, device and ultrasonic equipment, which can emit ultrasonic signals and detect echo signals of ultrasonic signals during ultrasonic inspection of a target detection part of an object; determine the background echo intensity of each detection point of the target detection part and the displacement peak of each detection point based on the detection result of the echo signal; determine the shear wave velocity when the shear wave propagates at each detection point based on the displacement peak of each detection point; for any detection point, determine the display parameters of the detection point based on the background echo intensity of the detection point and the displacement peak of the detection point, the display parameters are used to characterize the strength of the vibration amplitude caused at the detection point when the shear wave propagates; generate a shear wave image based on the shear wave velocity corresponding to each detection point and the display parameters of each detection point. Since the shear wave image is generated based on the display parameters of each detection position point and the shear wave velocity corresponding to each detection position point, the shear wave image can simultaneously reflect the shear wave velocity at each position of the target detection part and the vibration amplitude caused by the shear wave propagation at each position of the target detection part. In other words, the shear wave image can include richer shear wave characteristic information, thereby improving the effectiveness of the shear wave image and greatly improving the work efficiency of relevant staff.

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

[0034] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for describing the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application.

[0035] Figure 1 A schematic diagram of an application scenario of a shear wave image generation method provided in an embodiment of the present application;

[0036] Figure 2 A schematic diagram of an application scenario of another shear wave image generation method provided in an embodiment of the present application;

[0037] Figure 3 One of the exemplary flow charts of a shear wave image generation method provided in an embodiment of the present application;

[0038] Figure 4 A schematic diagram of a process for determining a displacement peak value of a detection position point provided in an embodiment of the present application;

[0039] Figure 5 A schematic diagram of the vibration displacement curves of various detection positions provided in the embodiment of the present application;

[0040] Figure 6 A schematic diagram of a vibration displacement curve of a detection position point provided in an embodiment of the present application;

[0041] Figure 7 One of the schematic diagrams of the shear wave velocity calculation process at the detection position point provided in the embodiment of the present application;

[0042] Figure 8 One of the schematic diagrams of the shear wave velocity calculation process at the detection position point provided in the embodiment of the present application;

[0043] Fig. 9 One of the schematic diagrams of the shear wave velocity calculation process at the detection position point provided in the embodiment of the present application;

[0044] Fig.10 One of the schematic diagrams of the shear wave velocity calculation process at the detection position point provided in the embodiment of the present application;

[0045] Fig.11 One of the schematic diagrams of the shear wave velocity calculation process at the detection position point provided in the embodiment of the present application;

[0046] Fig.12 One of the schematic diagrams of the shear wave velocity calculation process at the detection position point provided in the embodiment of the present application;

[0047] Fig.13A A schematic diagram of the positions of the detection position point X and the detection position point Y provided in the embodiment of the present application at the target detection part;

[0048] Fig. 13B A schematic diagram of a vibration displacement curve of a detection position point X and a detection position point Y provided in an embodiment of the present application;

[0049] Fig.14 A schematic diagram of a shear wave velocity image provided in an embodiment of the present application;

[0050] Fig.15 A schematic diagram of a shear wave image provided in an embodiment of the present application;

[0051] Fig.16 A schematic diagram of the system structure of a method for generating a shear wave image using an embodiment of the present application;

[0052] Fig.17 This is one of the exemplary flow charts of a shear wave image generation method provided in an embodiment of the present application. DETAILED DESCRIPTION

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

[0054] The terms "first" and "second" in the embodiments of the present application are used to distinguish different objects, rather than to describe a specific order. In addition, the term "comprising" and any of their variations are intended to cover non-exclusive protection. For example, a process, method, system, product or device comprising a series of steps or units is not limited to the listed steps or units, but optionally also includes steps or units that are not listed, or optionally also includes other steps or units inherent to these processes, methods, products or devices. "Multiple" in the present application can mean at least two, for example, two, three or more, and the embodiments of the present application are not limited.

[0055] In addition, the term "and / or" in this article is only a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in this article, unless otherwise specified, generally indicates that the associated objects before and after are in an "or" relationship.

[0056] At present, shear wave elastic imaging technology can realize real-time quantitative hardness detection of biological tissues, providing a basis for judging the pathological conditions of biological tissues. However, in practical applications, due to the complex structure of human tissues, shear waves will split and refract when passing through different tissue interfaces, but these situations cannot be distinguished only by the shear wave velocity. Therefore, how to provide a shear wave image generation method to make the shear wave image include richer shear wave characteristic information is urgently needed.

[0057] In view of this, an embodiment of the present application provides a shear wave image generation method, device and ultrasonic equipment, in which an ultrasonic signal can be emitted and the echo signal of the ultrasonic signal can be detected during the ultrasonic inspection of the target detection part of the object; based on the detection result of the echo signal, the background echo intensity of each detection point of the target detection part and the displacement peak of each detection point are determined; based on the displacement peak of each detection point, the shear wave velocity when the shear wave propagates at each detection point is determined; for any detection point, the display parameters of the detection point are determined according to the background echo intensity of the detection point and the displacement peak of the detection point, and the display parameters are used to characterize the strength of the vibration amplitude caused at the detection point when the shear wave propagates; based on the shear wave velocity corresponding to each detection point and the display parameters of each detection point, a shear wave image is generated. Since the shear wave image is generated based on the display parameters of each detection position point and the shear wave velocity corresponding to each detection position point, the shear wave image can simultaneously reflect the shear wave velocity at each position of the target detection part and the vibration amplitude caused by the shear wave propagation at each position of the target detection part. In other words, the shear wave image can include richer shear wave characteristic information, thereby improving the effectiveness of the shear wave image and greatly improving the work efficiency of relevant staff.

[0058] The structure of the ultrasonic equipment is introduced below with reference to the accompanying drawings. Figure 1 Schematic diagram of an application scenario of a shear wave image generation method provided in an embodiment of the present application. The method can be applied to Figure 1 The ultrasound device 10 shown includes a display screen 101, a processor 102 and at least one probe 103. Among them:

[0059] Processor 102 is used to detect the echo signal of the ultrasonic signal; based on the detection result of the echo signal, determine the background echo intensity of each detection point of the target detection part and the displacement peak of each detection point; based on the displacement peak of each detection point, determine the shear wave velocity when the shear wave propagates at each detection point; for any detection point, determine the display parameters of the detection point according to the background echo intensity of the detection point and the displacement peak of the detection point, the display parameters are used to characterize the strength of the vibration amplitude caused at the detection point when the shear wave propagates; based on the shear wave velocity corresponding to each detection point and the display parameters of each detection point, generate a shear wave image;

[0060] The display screen 101 is used to display the shear wave image.

[0061] See also Figure 2 A schematic diagram of an application scenario of another shear wave image generation method provided in an embodiment of the present application, such as Figure 2As shown, it includes an ultrasound device 10 and a server 21;

[0062] The ultrasound device 10 includes a display screen 101 and at least one probe 103;

[0063] The ultrasound device 10 is used to transmit an ultrasound signal, receive an echo signal of the ultrasound signal, and send a shear wave image generation instruction to the server 21 after receiving the echo signal during an ultrasound examination of a target detection part of an object;

[0064] The server 21 detects the echo signal of the ultrasonic signal according to the shear wave image generation instruction; determines the background echo intensity of each detection position point of the target detection part and the displacement peak value of each detection position point based on the detection result of the echo signal; determines the shear wave velocity when the shear wave propagates at each detection position point based on the displacement peak value of each detection position point; for any detection position point, determines the display parameters of the detection position point according to the background echo intensity of the detection position point and the displacement peak value of the detection position point, and the display parameters are used to characterize the strength of the vibration amplitude caused at the detection position point when the shear wave propagates; generates a shear wave image based on the shear wave velocity corresponding to each detection position point and the display parameters of each detection position point;

[0065] The ultrasound apparatus 10 displays the shear wave image generated by the server 21 on the display screen 102 .

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

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

[0068] See also Figure 3 This is one of the exemplary flow charts of a shear wave image generation method provided in an embodiment of the present application. The method can be applied to Figure 1 or Figure 2 The ultrasound equipment shown may include the following procedures:

[0069] S301, during an ultrasonic inspection of a target detection part of an object, an ultrasonic signal is emitted and an echo signal of the ultrasonic signal is detected.

[0070] During the ultrasonic inspection of the target detection part of the object, the ultrasonic device can use the acoustic radiation force of the ultrasonic signal to excite the target detection part of the object to generate shear waves, and detect the target detection part affected by the shear wave by emitting ultrasonic signals multiple times, obtain echo signals, and detect the echo signals through orthogonal demodulation and downsampling processing to obtain the detection results of the echo signals. The detection results may include the in-phase data (I data) and orthogonal data (Q data) of each detection position point of the target detection part.

[0071] S302, based on the detection result of the echo signal, determine the background echo intensity of each detection position point of the target detection part and the displacement peak value of each detection position point.

[0072] In one possible implementation, when the ultrasonic device determines the background echo intensity of each detection position point of the target detection part based on the detection result of the echo signal, it can perform a modulo operation on the I data and Q data of any detection position point to obtain an envelope signal, thereby determining the background echo intensity of the detection position point, satisfying formula (1).

[0073]

[0074] Where Y(j) represents the background echo intensity of the j-th detection position, I(j) represents the data of the j-th detection position in the I data, and Q(j) represents the data of the j-th detection position in the Q data.

[0075] Optionally, j may also refer to an index (or serial number) in a data sequence. For example, if I data includes N data, then I(j) represents the jth data in the N data, where 1≤j≤N.

[0076] In some embodiments, when determining the displacement peak value of each detection position point of the target detection part, the ultrasound device can perform the following steps for each detection position point: Figure 4 The process shown:

[0077] S401, based on the detection result corresponding to the detection position point in the detection result of the echo signal, using a two-dimensional autocorrelation algorithm to determine the axial displacement of the detection position point within a preset time period corresponding to each moment.

[0078] Here, the axial direction refers to the longitudinal direction, that is, the depth direction.

[0079] In the embodiment of the present application, when the ultrasonic device determines the displacement peak value of each detection position point of the target detection part, it is necessary to first estimate the vibration displacement of each detection position point. Specifically, the ultrasonic device can use a two-dimensional autocorrelation algorithm to determine the axial displacement of the detection position point within a preset time period corresponding to each moment based on the detection result corresponding to the detection position point in the detection result of the echo signal, and satisfy formula (2).

[0080]

[0081] In the formula, is the axial displacement of the detection position point within the preset time length corresponding to a moment; M is the sampling window size in the depth direction, N is the sampling window size in the time direction, m is the depth, n is the time, and c is the propagation speed of the sound in the target detection part; f c is the center frequency of the echo signal.

[0082] S402, determining the maximum value of each axial displacement of the detection position point by a maximum value search method, and taking the maximum value as the displacement peak value of the detection position point.

[0083] In a possible implementation, after the displacement corresponding to the detection position point at each moment is determined by formula (2), a method of separating shear waves propagating in different directions can be used for each detection position point to draw a vibration displacement curve of each detection position point along the time direction. Figure 5 A schematic diagram of the vibration displacement curves of each detection position point provided in an embodiment of the present application. Figure 5 In the diagram, the horizontal axis represents time, the vertical axis represents vibration amplitude, and the magnitude of displacement can be represented by the magnitude of vibration amplitude. Different curves can be used to represent the vibration displacement curves of the detection position points at the same depth direction and different time directions. For example, when m represents depth and n represents time, Figure 5 Different curves in the figure can represent detection position points with the same m value but different n values.

[0084] by Figure 5 As an example, see Figure 6 A schematic diagram of a vibration displacement curve of a detection position point provided in an embodiment of the present application is shown in FIG. Figure 6 As shown, the point with the largest amplitude in the curve is A, then the amplitude corresponding to A is the displacement peak value of the detection position point, and the time corresponding to A is the acquisition time corresponding to the displacement peak value.

[0085] It should be noted that the maximum value search method can use commonly used maximum value search methods such as the sorting method to find the maximum value, and this application does not limit this.

[0086] Based on the above scheme, since the shear wave will cause displacement deformation of the tissue during the process of tissue propagation, the ultrasonic signal is used to obtain the detection results of the shear wave passing through the target detection part, and the shear wave displacement curve of the target detection part is reconstructed through the displacement deformation. The displacement peak value of each detection position point can reflect the strength of the vibration amplitude caused by the shear wave at the detection position point when the shear wave propagates in the target detection part, thereby reflecting richer shear wave characteristic information when generating the shear wave image.

[0087] S303, determining the shear wave velocity when the shear wave propagates at each detection position point based on the displacement peak value of each detection position point.

[0088] The ultrasonic device can determine the shear wave velocity when the shear wave propagates at each detection position point based on the position of each detection position point at the target detection part and the acquisition time when each detection position point reaches the displacement peak.

[0089] In some embodiments, when the ultrasonic device calculates the shear wave velocity corresponding to each detection position point, the following operations can be performed for any detection position point: Determine the N detection position points included in the sliding window corresponding to a detection position point. Through the least squares method, the displacement peak value of each detection position point in the N detection position points and the acquisition time corresponding to the displacement peak value are fitted to obtain the average shear wave velocity of the N detection position points. The average shear wave velocity is used as the shear wave velocity of the detection position point. Among them, the size of the sliding window can be set according to actual conditions or experience but is less than or equal to the size of the ultrasonic image of the target detection part collected by the ultrasonic device. N is a positive integer. When the sliding window size is different, N can be different.

[0090] Specifically, the ultrasonic device can determine the vibration displacement curve of each detection position point among the N detection position points, and determine the N displacement peak points corresponding to when the displacement of the detection position point in the curve reaches the displacement peak. Then, the N displacement peak points are subjected to straight line fitting processing by the least squares method, and the slope of the fitted straight line is used as the average shear wave velocity of the N detection position points.

[0091] The following is a detailed description of a method for ultrasonic equipment to calculate the shear wave velocity corresponding to each detection position point using a specific embodiment.

[0092] See also Figure 7-Figure 12 A schematic diagram of the shear wave velocity calculation process for the detection position point provided in the embodiment of the present application, Figure 7 The target detection part shown includes detection position point A-detection position point E. The preset sliding window size is 3*1, and the step size is 1. Then, in the first sliding, the sliding window may include detection position point A, detection position point B, and detection position point C, and determine Figure 8The vibration displacement curves of the detection position A, detection position B, and detection position C are shown. Figure 8 It can be seen that in the vibration displacement curve of the detection position point A, the displacement peak point is O, in the vibration displacement curve of the detection position point B, the displacement peak point is P, and in the vibration displacement curve of the detection position point C, the displacement peak point is Q. Then, O, P, and Q can be subjected to straight line fitting processing by the least squares method, and the slope of the fitted straight line is used as the average shear wave velocity of the detection position points A, B, and C. And the average shear wave velocity is used as the shear wave velocity of the detection position point B.

[0093] like Fig. 9 As shown, during the second sliding, the sliding window may include detection position point B, detection position point C, and detection position point D, and determine as follows Fig.10 The vibration displacement curves of the detection position B, detection position C, and detection position D are shown. Fig.10 It can be seen that in the vibration displacement curve of the detection position point B, the displacement peak point is P, in the vibration displacement curve of the detection position point C, the displacement peak point is Q, and in the vibration displacement curve of the detection position point D, the displacement peak point is R. Then, P, Q, and R can be subjected to straight line fitting processing by the least squares method, and the slope of the fitted straight line is used as the average shear wave velocity of the detection position points B, C, and D. And the average shear wave velocity is used as the shear wave velocity of the detection position point C.

[0094] like Fig.11 As shown, during the second sliding, the sliding window may include detection position point C, detection position point D, and detection position point E, and determine as follows Fig.12 The vibration displacement curves of the detection position points C, D, and E are shown. Fig.12 It can be seen that in the vibration displacement curve of the detection position point C, the displacement peak point is Q, in the vibration displacement curve of the detection position point D, the displacement peak point is R, and in the vibration displacement curve of the detection position point E, the displacement peak point is S. Then, Q, R, and S can be linearly fitted by the least squares method, and the slope of the fitted straight line is used as the average shear wave velocity of the detection position points C, D, and E. And the average shear wave velocity is used as the shear wave velocity of the detection position point D.

[0095] Optionally, for detection position point A and detection position point E located at the edge of the target detection part, the shear wave velocity of detection position point B can be used as the shear wave velocity of detection position point A, and the shear wave velocity of detection position point D can be used as the shear wave velocity of detection position point E. It should be understood that the above-mentioned method for determining the shear wave velocity of detection position point A and detection position point E at the edge of the target detection part is only exemplary, and other methods can also be used for calculation, which is not limited in the present application.

[0096] In one example, the ultrasonic device can also use the time of flight (TOF) method to calculate the shear wave velocity corresponding to each detection position point. Specifically, for any two detection position points, the vibration displacement curves of the two detection position points can be determined, and the peak moment when the displacement of each detection position point reaches the displacement peak can be determined through the displacement peak point of the curve. Based on the difference in peak moments and the positions of the two detection position points in the target detection area, the average shear wave velocity of the two detection position points is calculated, and the average shear wave velocity is used as the shear wave velocity of each of the two detection position points.

[0097] For example, see Fig.13A A schematic diagram of the positions of the detection position point X and the detection position point Y at the target detection part provided in the embodiment of the present application, such as Fig.13A As shown, the detection position point X and the detection position point Y are two detection position points for calculating the shear wave velocity. According to the positions of the two detection position points in the target detection part, the distance between the detection position point X and the detection position point Y in the target detection part can be determined as d. Fig. 13B As shown in the vibration displacement curves of the detection position point X and the detection position point Y, it can be seen that the peak time of the detection position point X is t x , the peak time of the detection position point Y is t y Then the average shear wave velocity of the detection position point X and the detection position point Y is shown in formula (3). And the average shear wave velocity is used as the shear wave velocity corresponding to the detection position point X and the shear wave velocity corresponding to the detection position point Y respectively.

[0098]

[0099] S304: for any detection position point, determine the display parameters of the detection position point according to the background echo intensity of the detection position point and the displacement peak value of the detection position point.

[0100] Among them, the display parameter is used to characterize the strength of the vibration amplitude caused by the shear wave propagation at the detection location.

[0101] In a possible implementation, when the ultrasonic device determines the display parameters corresponding to the detection position point based on the background echo intensity of the detection position point and the displacement peak value of the detection position point, the background echo intensity of the detection position point can be used as the ambient light impact parameter in the Lambert illumination model, and the normal vector height corresponding to the detection position point can be corrected according to the displacement peak value. Then, based on the Lambert illumination model and the normal vector height corresponding to the corrected detection position point, the diffuse reflection value of the detection position point is determined. Then, based on the mapping relationship between the diffuse reflection value and the display parameter, the display parameter corresponding to the diffuse reflection value of the detection position point is determined.

[0102] Since the Lambert illumination model is independent of the viewing angle, the surface illumination distribution is not affected by the change of viewing angle, so the Lambert illumination model can be used to perform stereoscopic display processing on the shear wave velocity image. The Lambert illumination model satisfies the following formula: diffuse reflection value = reflection coefficient × ambient light + reflection coefficient × (vertex normal vector * incident light reverse direction vector).

[0103] Since there is usually still some light in the world (such as the moon, distant lights, etc.) even in darkness, objects are almost never completely dark, so the influence of ambient light usually needs to be considered in diffuse reflection. In the embodiment of the present application, the background echo intensity of each detection position point can be normalized and used as the ambient light influence parameter in the Lambert illumination model.

[0104] (vertex normal vector * incident light reverse direction vector) represents diffuse reflection illumination, which is used to simulate the directional influence of light source on objects, so that the deflection on the object that is closer to the light direction can obtain more brightness from the light source. In this application, the normal vector height corresponding to the detection position point is corrected according to the displacement peak value, that is, the displacement peak value is used as the perturbation amount in the same direction of the diffuse reflection vertex normal vector in the illumination model, so that the plane has a concave and convex effect.

[0105] In addition, in the embodiment of the present application, the focal position of the acoustic radiation force of the ultrasonic signal is used as the light source position. The Lambert illumination model is adjusted, and the adjusted Lambert illumination model satisfies the following formula: diffuse reflection value = reflection coefficient × background echo intensity + reflection coefficient × (vertex normal vector * shear wave displacement peak * incident light reverse direction vector).

[0106] The diffuse reflection value of each detection position point can be determined by the adjusted Lambert illumination model, and then the display parameter corresponding to the diffuse reflection value of the detection position point can be determined according to the mapping relationship between the diffuse reflection value and the display parameter. Among them, the display parameter can be represented by a pixel value between 0 and 255. That is, the brightness of the detection position point represents the strength of the vibration amplitude caused by the shear wave propagation at the detection position point, 0 represents the brightest, and 255 represents the darkest. The mapping relationship between the diffuse reflection value and the display parameter can be set according to actual conditions or experience. For example, the diffuse reflection value can be divided into 256 parts from small to large, and the diffuse reflection value in each part corresponds to a display parameter, etc. This application does not limit this.

[0107] In addition, the larger the diffuse reflection value, the larger the display parameter, or the smaller the diffuse reflection value, the larger the display parameter, and this application does not limit this.

[0108] S305: Generate a shear wave image based on the shear wave velocity corresponding to each detection position point and the display parameters of each detection position point.

[0109] In one possible implementation, the ultrasonic device performs pseudo-color mapping processing based on the shear wave velocity corresponding to each detection position point to determine the color value corresponding to each detection position point. Specifically, the ultrasonic device can determine the color value corresponding to each detection position point according to the preset mapping relationship between the shear wave velocity and the color value and the shear wave velocity corresponding to each detection position point. Among them, the color value can refer to the RGB channel value corresponding to a certain color, and the color value corresponds to the color one by one.

[0110] For example, the preset mapping relationship between shear wave velocity and color value may include that when the shear wave velocity is a, the color value is R=255, G=0, B=0. When the shear wave velocity is b, the color value is R=255, G=255, B=0. When the shear wave velocity is c, the color value is R=0, G=0, B=255.

[0111] Optionally, in the preset mapping relationship between shear wave velocity and color value, the color value can be from blue to yellow to red corresponding to the shear wave velocity from small to large, that is, red is displayed when the shear wave velocity is large, yellow is displayed when the shear wave velocity is medium, and blue is displayed when the shear wave velocity is small.

[0112] Optionally, the ultrasound device may also generate a shear wave velocity image based on the color values ​​corresponding to each detection position point, and may also display the shear wave velocity image on a display screen of the ultrasound device.

[0113] In some embodiments, when the ultrasound device generates a shear wave velocity image based on the color values ​​corresponding to each detection position point, it can generate a pseudo-color mapping image based on the corresponding colors of the color values ​​corresponding to each detection position point, and superimpose the pseudo-color mapping image on the ultrasound grayscale image to generate a shear wave velocity image.

[0114] See also Fig.14 A schematic diagram of a shear wave velocity image provided in an embodiment of the present application. Fig.14 The shear wave velocity of the detection position point in region 1 is greater than the shear wave velocity of the detection position point in region 2, and the shear wave velocity of the detection position point in region 2 is greater than the shear wave velocity of the detection position point in region 3. Therefore, the color value corresponding to region 1 has a larger R channel value, and smaller G and B channel values, and can be displayed as red or a color close to red. The color value corresponding to region 2 has a larger R and G channel value, and a smaller B channel value, and can be displayed as yellow or a color close to yellow. The color value corresponding to region 3 has a smaller R and G channel value, and a larger B channel value, and can be displayed as blue or a color close to blue.

[0115] It should be understood that the mapping relationship between the shear wave velocity and the color value may be preset based on experience or actual conditions, and the present application does not limit this.

[0116] In some embodiments, the ultrasonic device can use the display parameters corresponding to each detection position point, the color value corresponding to each detection position point, and the position of the detection position point in the target detection part as parameters to perform stereoscopic display processing on the shear wave velocity image through a stereoscopic effect display function to generate a shear wave image. Among them, the stereoscopic effect display function can be selected according to actual conditions or experience, for example, it can be Surfl function, Surf function, etc., and this application does not limit this.

[0117] In some embodiments, the detection position points located at the edge of the target detection part can be deleted after the shear wave velocity is calculated, and the color values ​​corresponding to the remaining detection position points are determined according to the shear wave velocities of the remaining detection position points. Then, the shear wave image is generated according to the display parameters corresponding to the remaining detection position points. Figure 7 Taking the detection position points A to E included in the target detection part as an example, after calculating the shear wave velocity, the shear wave velocity of detection position points B, C, and D can be obtained. Then, the detection position points A and E located at the edge of the target detection part can be deleted, and the color values ​​corresponding to the detection position points B, C, and D can be determined. Then, according to the display parameters corresponding to the detection position points B, C, and D, a shear wave image is generated.

[0118] In a possible implementation, the ultrasound device may also generate a shear wave displacement image according to the display parameters, and superimpose the shear wave displacement image on the shear wave velocity image to generate a shear wave image.

[0119] See also Fig.15 A schematic diagram of a shear wave image provided in an embodiment of the present application. In the embodiment of the present application, it is taken as an example that the larger the displacement peak value, the smaller the display parameter. Fig.15 The smaller the display parameter, the brighter the brightness. Fig.15 The brighter area shown is the location of the detection point with a larger displacement peak, and the darker area is the location of the detection point with a lower displacement peak. Through the distribution of the brighter and darker areas, and the principle of shear wave intensity attenuation from high to low, it can be seen that Fig.15 The shear wave propagation path is from left to right. It can be seen that the shear wave image generated by the shear wave image generation method provided by the present application can not only display the shear wave velocity by color, but also display the vibration displacement of the detection position point by displaying brightness, that is, the strength of the vibration amplitude of the detection position point caused by the shear wave propagation. More shear wave characteristics can be displayed through the shear wave image, thereby providing richer information reference.

[0120] Below, in order to more clearly understand the solution proposed in the embodiment of the present application, a shear wave image generation method provided by the present application will be introduced in combination with a specific embodiment.

[0121] See also Fig.16 , is a schematic diagram of a system structure of a shear wave image generation method provided in an embodiment of the present application. In the system 1600, an ultrasound transmitting module 1601, an ultrasound receiving module 1602, a processing module 1603 and a display module 1604 are included. Fig.16 The system shown can be executed Fig.17 The process shown.

[0122] S1701, Acoustic Radiation Force Focused Impact.

[0123] The ultrasonic transmitting module 1601 can use the acoustic radiation force of the ultrasonic signal to focus and impact the target detection part of the object to generate shear waves, so that the shear waves propagate in the target detection part. Then the ultrasonic transmitting module 1601 transmits the ultrasonic signal.

[0124] S1702, receiving an echo signal.

[0125] The ultrasonic receiving module 1602 receives the echo signal of the ultrasonic signal. It should be noted that the ultrasonic receiving module can be a traditional ultrasonic receiving and processing mode in the related technology, and the specific processing flow in the ultrasonic receiving module will not be repeated in this application.

[0126] S1703, calculate the background echo intensity.

[0127] The processing module 1603 performs orthogonal demodulation and downsampling processing on the echo signal to obtain the detection result of the echo signal. The background echo intensity can be calculated according to the detection result. For specific methods, see Figure 3 The relevant description in the method embodiment shown will not be repeated here.

[0128] S1704, shear wave displacement estimation.

[0129] The processing module 1603 can estimate the shear wave displacement according to the detection results of each detection position point of the target detection part, and determine the axial displacement of each detection position point within the preset time length corresponding to each moment. The specific method can be found in Figure 3 The relevant description in the method embodiment shown will not be repeated here.

[0130] S1705, calculate the displacement peak value.

[0131] The processing module 1603 can calculate the displacement peak value of the corresponding detection position point according to the shear wave displacement calculated at each detection position point in S1704. For specific methods, see Figure 3 The relevant description in the method embodiment shown will not be repeated here.

[0132] S1706, Shear wave velocity estimation.

[0133] The processing module 1603 can estimate the shear wave velocity of each detection position point according to the displacement peak value of each detection position point, and determine the shear wave velocity corresponding to each detection position point. The specific method can be found in Figure 3 The relevant description in the method embodiment shown will not be repeated here.

[0134] S1707, generate a shear wave image.

[0135] The display module 1604 can generate and display a shear wave image according to the color value corresponding to each detection position point and the display parameter corresponding to each detection position point. The specific method of generating a shear wave image can be found in Figure 3 The relevant description in the method embodiment shown will not be repeated here.

[0136] The embodiment of the present application further provides a computer-readable storage medium on which computer instructions are stored. When the computer instructions are executed on a computer, the computer executes the steps of any of the above methods.

[0137] Based on the same inventive concept, the embodiment of the present application also provides a computer program product, which includes: computer program code, when the computer program code is run on a computer, the computer executes any of the shear wave image generation methods discussed above. Since the principle of solving the problem by the above computer program product is similar to that of the shear wave image generation method, the implementation of the above computer program product can refer to the implementation of the method, and the repeated parts will not be repeated.

[0138] Those skilled in the art can understand that all or part of the steps of implementing the above-mentioned method embodiments can be completed by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, the steps of the above-mentioned method embodiments are executed; and the aforementioned storage medium includes: ROM, RAM, disk or optical disk and other media that can store program codes.

[0139] Although the specific embodiments of the present application are described above, it should be understood by those skilled in the art that these are merely examples, and the scope of protection of the present application is limited by the attached claims. Those skilled in the art may make various changes or modifications to these embodiments without departing from the principles and essence of the present application, but these changes and modifications all fall within the scope of protection of the present application. Although the preferred embodiments of the present application have been described, once those skilled in the art are aware of the basic creative concepts, additional changes and modifications may be made to these embodiments. Therefore, the attached claims are intended to be interpreted as including the preferred embodiments and all changes and modifications falling within the scope of the present application.

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

Claims

1. An ultrasonic device, characterized in that: comprising a display screen, a processor and at least one probe; The at least one probe is used to transmit an ultrasonic signal during an ultrasonic inspection of a target inspection part of the object; The processor is used to detect the echo signal of the ultrasonic signal; based on the detection result of the echo signal, determine the background echo intensity of each detection point of the target detection part and the displacement peak of each detection point; based on the displacement peak of each detection point, determine the shear wave velocity when the shear wave propagates at each detection point; for any detection point, determine the display parameters of the detection point according to the background echo intensity of the detection point and the displacement peak of the detection point, the display parameters are used to characterize the strength of the vibration amplitude caused by the shear wave at the detection point when the shear wave propagates; Generate a shear wave image based on the shear wave velocity corresponding to each detection position point and the display parameters of each detection position point; The display screen is used to display the shear wave image.

2. The ultrasonic device according to claim 1, characterized in that When the processor determines, for any detection position point, the display parameter of the detection position point according to the background echo intensity of the detection position point and the displacement peak value of the detection position point, the processor is specifically used to: For any detection position point, the background echo intensity of the detection position point is used as the ambient light impact parameter in the Lambert illumination model, and the normal vector height corresponding to the detection position point is corrected according to the displacement peak value; Determine the diffuse reflection value of the detection position point based on the Lambert illumination model and the corrected normal vector height corresponding to the detection position point; The display parameters of the detection position point are determined according to the diffuse reflection value of the detection position point and the mapping relationship between the diffuse reflection value and the display parameter.

3. The ultrasonic device according to claim 1, characterized in that The processor determines the displacement peak value of each detection position point according to the following method: For any detection location point, perform the following operations: Based on the detection result corresponding to the detection position point in the detection result of the echo signal, using a two-dimensional autocorrelation algorithm to determine the axial displacement of the detection position point within a preset time length corresponding to each moment; The maximum value among the axial displacements of the detection position point is determined by the maximum value search method, and the maximum value is used as the displacement peak value of the detection position point.

4. The ultrasonic device according to any one of claims 1 to 3, characterized in that: The processor determines the shear wave velocity when the shear wave propagates at each detection position point based on the displacement peak value of each detection position point, and is specifically used to: For any detection location point, perform the following operations: Determine N detection position points included in the sliding window corresponding to the detection position point; N is a positive integer; The displacement peak value of each of the N detection positions and the acquisition time corresponding to the displacement peak value are fitted by the least squares method to obtain the average shear wave velocity of the N detection positions; and the average shear wave velocity is used as the shear wave velocity of the detection position.

5. The ultrasonic device according to any one of claims 1 to 3, characterized in that: The detection result of the echo signal includes in-phase data and orthogonal data of each detection position point of the target detection part; When the processor determines the background echo intensity of each detection position point of the target detection part based on the detection result of the echo signal, it is specifically used to: For any detection position point of the target detection part, a modulo operation is performed on the in-phase data and the orthogonal data of the detection position point to determine the background echo intensity of the detection position point.

6. The ultrasonic device according to any one of claims 1 to 3, characterized in that: When the processor generates a shear wave image based on the shear wave velocity corresponding to each detection position point and the display parameters of each detection position point, it is specifically used to: According to a preset mapping relationship between shear wave velocity and color value and the shear wave velocity corresponding to each detection position point, respectively determine the color value corresponding to each detection position point; A shear wave image is generated based on the color values ​​corresponding to the respective detection position points and the display parameters of the respective detection position points.

7. A shear wave image generation method, characterized in that: include: During an ultrasonic inspection of a target detection part of an object, an ultrasonic signal is emitted and an echo signal of the ultrasonic signal is detected; Based on the detection result of the echo signal, determine the background echo intensity of each detection position point of the target detection part and the displacement peak value of each detection position point; Based on the displacement peak value of each detection position point, the shear wave velocity when the shear wave propagates at each detection position point is determined; For any detection position point, a display parameter of the detection position point is determined according to the background echo intensity of the detection position point and the displacement peak value of the detection position point, wherein the display parameter is used to characterize the vibration amplitude strength caused by the shear wave propagation at the detection position point; A shear wave image is generated based on the shear wave velocity corresponding to each detection position point and the display parameters of each detection position point.

8. The method according to claim 7, characterized in that The step of determining, for any detection position point, the display parameter of the detection position point according to the background echo intensity of the detection position point and the displacement peak value of the detection position point specifically includes: For any detection position point, the background echo intensity of the detection position point is used as the ambient light impact parameter in the Lambert illumination model, and the normal vector height corresponding to the detection position point is corrected according to the displacement peak value; Determine the diffuse reflection value of the detection position point based on the Lambert illumination model and the corrected normal vector height corresponding to the detection position point; The display parameters of the detection position point are determined according to the diffuse reflection value of the detection position point and the mapping relationship between the diffuse reflection value and the display parameter.

9. The method according to claim 7, characterized in that: The peak displacement of each detection position is determined in the following way: For any detection location point, perform the following operations: Based on the detection result corresponding to the detection position point in the detection result of the echo signal, using a two-dimensional autocorrelation algorithm to determine the axial displacement of the detection position point within a preset time length corresponding to each moment; The maximum value among the axial displacements of the detection position point is determined by the maximum value search method, and the maximum value is used as the displacement peak value of the detection position point.

10. A shear wave image generating device, characterized in that: include: A transmitting unit, for transmitting an ultrasonic signal during an ultrasonic inspection of a target detection part of an object; A receiving unit, configured to receive an echo signal of the ultrasonic signal; A processing unit is used to perform the following processing: detecting an echo signal of the ultrasonic signal; Based on the detection result of the echo signal, determine the background echo intensity of each detection position point of the target detection part and the displacement peak value of each detection position point; Determine the shear wave velocity when the shear wave propagates at each detection position point based on the displacement peak value of each detection position point; For any detection position point, a display parameter of the detection position point is determined according to the background echo intensity of the detection position point and the displacement peak value of the detection position point, wherein the display parameter is used to characterize the vibration amplitude strength caused by the shear wave propagation at the detection position point; Generate a shear wave image based on the shear wave velocity corresponding to each detection position point and the display parameters of each detection position point; A display unit is used to display the shear wave image.

Images