CN113040816A - Ultrasonic elastography method, device, electronic equipment and storage medium - Google Patents
Ultrasonic elastography method, device, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN113040816A CN113040816A CN202110369497.4A CN202110369497A CN113040816A CN 113040816 A CN113040816 A CN 113040816A CN 202110369497 A CN202110369497 A CN 202110369497A CN 113040816 A CN113040816 A CN 113040816A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- detection image
- tissue
- ultrasound
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000002091 elastography Methods 0.000 title abstract description 27
- 238000001514 detection method Methods 0.000 claims abstract description 249
- 238000012360 testing method Methods 0.000 claims abstract description 56
- 238000003384 imaging method Methods 0.000 claims abstract description 38
- 230000005284 excitation Effects 0.000 claims abstract description 34
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 22
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 22
- 210000001519 tissue Anatomy 0.000 claims description 112
- 238000002604 ultrasonography Methods 0.000 claims description 48
- 238000002113 ultrasound elastography Methods 0.000 claims description 25
- 238000007689 inspection Methods 0.000 claims description 24
- 238000001914 filtration Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 15
- 238000012285 ultrasound imaging Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 13
- 230000003902 lesion Effects 0.000 claims description 13
- 210000005228 liver tissue Anatomy 0.000 claims description 11
- 239000003550 marker Substances 0.000 claims description 11
- 239000000523 sample Substances 0.000 claims description 10
- 230000002146 bilateral effect Effects 0.000 claims description 6
- 230000001427 coherent effect Effects 0.000 claims description 6
- 239000003086 colorant Substances 0.000 claims description 6
- 230000002708 enhancing effect Effects 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 6
- 239000002131 composite material Substances 0.000 claims 1
- 230000008859 change Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 12
- 238000003908 quality control method Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000029058 respiratory gaseous exchange Effects 0.000 description 5
- 239000006185 dispersion Substances 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000009278 visceral effect Effects 0.000 description 2
- 208000004930 Fatty Liver Diseases 0.000 description 1
- 206010019668 Hepatic fibrosis Diseases 0.000 description 1
- 206010019708 Hepatic steatosis Diseases 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 208000010706 fatty liver disease Diseases 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 238000001310 location test Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 231100000240 steatosis hepatitis Toxicity 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The embodiment of the application provides an ultrasonic elastography method, an ultrasonic elastography device, an electronic device and a storage medium, wherein an ultrasonic echo signal is obtained by applying vibration excitation to a tissue to be detected and performing ultrasonic detection; performing beam synthesis according to the ultrasonic echo signals to generate synthesized ultrasonic imaging data; generating an ultrasonic detection image and an elasticity detection image by using the synthesized ultrasonic imaging data, wherein the elasticity detection image is used for representing the tissue elasticity detectability information and the signal quality of the ultrasonic detection image at the corresponding moment; according to the position relation 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 displayed in an overlapped mode.
Description
Technical Field
The present application relates to the field of medical equipment technologies, and in particular, to an ultrasound elastography method, an ultrasound elastography device, an electronic device, and a storage medium.
Background
Currently, the ultrasound imaging technology is an examination technology widely used in clinical examination, and detects a lesion of a biological tissue by detecting a physiological structure of the biological tissue using ultrasound. However, in the early lesion stage of the biological tissue, the structure of the biological tissue does not change obviously, so that the traditional ultrasonic imaging method such as B ultrasonic and the like is not sensitive to the early lesion of the biological tissue. On the basis of the ultrasonic imaging technology, the elastography technology which is made up in recent years is a novel imaging technology, the vibration of biological tissues is excited by applying shear waves to the biological tissues, the form change of the detected tissues before and after the vibration stress is generated is tracked by ultrasonic waves, so that the mechanical characteristics of the detected tissues are obtained, the state of the biological tissues is detected according to the mechanical characteristics of the biological tissues, and whether early lesion characteristics exist is determined.
In the prior art, the result of elasticity calculation is usually combined with the ultrasound image in the form of a parameter value to form an elasticity imaging result, and the elasticity imaging result is displayed to an operator of the device, and the operator needs to observe the ultrasound image and combine the quantitative parameter value of elasticity measurement to realize the state detection of the detected tissue.
However, in the process of performing elasticity measurement, the optimal test position needs to be determined after multiple times of observation and adjustment, so the ultrasonic elasticity imaging method in the prior art has the problems of low detection efficiency and poor detection accuracy.
Disclosure of Invention
The application provides an ultrasonic elastography method, an ultrasonic elastography device, electronic equipment and a storage medium, which are used for solving the problems of low detection efficiency and poor detection accuracy in the prior art.
According to a first aspect of embodiments herein, there is provided an ultrasound elastography method, the method comprising: applying vibration excitation to the tissue to be detected and carrying out ultrasonic detection to obtain an ultrasonic echo signal; performing beam synthesis according to the ultrasonic echo signals to generate synthesized ultrasonic imaging data; generating an ultrasonic detection image and an elasticity detection image by using the synthesized ultrasonic imaging data, wherein the elasticity detection image is used for representing the tissue elasticity detectability information and the signal quality of the ultrasonic detection image at the corresponding moment; and displaying the ultrasonic detection image and the elastic detection image in an overlapping manner according to the position relation of the test region corresponding to the ultrasonic detection image and the test region corresponding to the elastic detection image.
In one possible implementation, applying a vibration excitation to a tissue to be detected and performing an ultrasonic detection to obtain an ultrasonic echo signal includes: applying shear waves to the tissue to be detected, and acquiring ultrasonic echo signals for tracking shear wave propagation; performing beam synthesis according to the ultrasonic echo signal to generate synthesized ultrasonic imaging data, including: and performing beam synthesis on the ultrasonic echo signals for tracking the shear wave propagation to generate multi-frame scanning lines, wherein the multi-frame scanning lines are used for representing ultrasonic detection results at different positions.
In one possible implementation, the generating an ultrasound inspection image and an elasticity inspection image using the synthesized ultrasound imaging data includes: processing the synthesized ultrasound imaging data by a graphics processor, generating the ultrasound detection image and the elasticity detection image in parallel, and obtaining an elastic modulus E according to a shear wave propagation velocity, wherein E is 3 rho V2ρ is the tissue density and V is the shear wave velocity.
In one possible implementation, before displaying after the generation of the ultrasound inspection image, the following processing is further performed on the ultrasound inspection image: carrying out median filtering to remove electronic noise; performing Gaussian smoothing and non-local mean filtering to smooth the structural information; carrying out bilateral filtering and enhancing edge information; and taking two thresholds of the signal as display thresholds, wherein the lower threshold is the amplitude of the signal from 1% to 20%, and the upper threshold is the amplitude of the signal from 60% to 200%, and performing histogram equalization processing according to the display thresholds.
In a possible implementation manner, the elastic detection image includes at least one position mark, where the position mark is used to represent a propagation position of the shear wave generated by the vibration excitation in the tissue to be detected, and the ultrasound detection image and the elastic detection image are displayed in an overlapping manner according to a positional relationship between a test area corresponding to the ultrasound detection image and a test area corresponding to the elastic detection image, including: determining a positioning coordinate system by taking the test area of the ultrasonic detection image as a reference; determining the position coordinates of the position marks in the positioning coordinate system according to the position relation between the position marks and the test area; and displaying the position mark on the ultrasonic detection image in an overlapping manner according to the position coordinates.
In a possible implementation manner, the elastic detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color, and the pseudo-color is used for representing corresponding time when the shear wave arrives at different positions when being transmitted inside the tissue to be detected.
In one possible implementation, each of the position markers has a different smoothness, which is used to characterize the propagation stability of the shear wave.
In a possible implementation manner, after obtaining the ultrasonic echo signal, the method further includes: calculating an ultrasonic attenuation parameter UAP according to the formula of (I1-I2)/(D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2; calculating an ultrasonic attenuation coefficient alpha according to the formula of alpha-20 log (I1/I2) -6)/2 (D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2; calculating a scattering son distribution coefficient k according to the formula k which is s/sigma; where s is the coherent signal energy; σ is the spread signal energy; k ranges from [0, 1 ]; after the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient alpha and the scatterer distribution coefficient k are normalized, calculating a liver tissue lesion comprehensive Score according to a formula of Score (a) E + b UAP + c a + d) k according to an ultrasonic echo signal and corresponding histopathological typing; wherein, 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 ].
According to a second aspect of embodiments of the present application, there is provided an ultrasound elastography device comprising:
the detection module is used for applying vibration excitation to the tissue to be detected and carrying out ultrasonic detection to obtain an ultrasonic echo signal;
the synthesis module is used for carrying out beam synthesis according to the ultrasonic echo signals to generate synthesized ultrasonic imaging data;
a generating module, configured to generate an ultrasound detection image and an elasticity detection image by using the synthesized ultrasound imaging data, where the elasticity detection image is used to represent tissue elasticity detectability information and signal quality of the ultrasound detection image at a corresponding time;
and the display module is used for displaying the ultrasonic detection image and the elastic detection image in an overlapping manner according to the position relation between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image.
In a possible implementation manner, the detection module is specifically configured to: applying shear waves to the tissue to be detected, and acquiring ultrasonic echo signals for tracking shear wave propagation; the synthesis module is specifically configured to: and performing beam synthesis on the ultrasonic echo signals for tracking the shear wave propagation to generate multi-frame scanning lines, wherein the multi-frame scanning lines are used for representing ultrasonic detection results at different positions.
In one possible implementation, the generating module has a function of: processing the synthesized ultrasound imaging data by a graphics processor, generating the ultrasound detection image and the elasticity detection image in parallel, and obtaining an elastic modulus E according to a shear wave propagation velocity, wherein E is 3 rho V2ρ is the tissue density and V is the shear wave velocity.
In a possible implementation manner, after generating the ultrasound inspection image, the generating module further performs the following processing on the ultrasound inspection image: carrying out median filtering to remove electronic noise; performing Gaussian smoothing and non-local mean filtering to smooth the structural information; carrying out bilateral filtering and enhancing edge information; and taking two thresholds of the signal as display thresholds, wherein the lower threshold is the amplitude of the signal from 1% to 20%, and the upper threshold is the amplitude of the signal from 60% to 200%, and performing histogram equalization processing according to the display thresholds.
In a possible implementation manner, the elasticity detection image includes at least one position mark, the position mark is used for characterizing a propagation position of the shear wave generated by the vibration excitation inside the tissue to be detected, and the display module is specifically configured to: determining a positioning coordinate system by taking the test area of the ultrasonic detection image as a reference; determining the position coordinates of the position marks in the positioning coordinate system according to the position relation between the position marks and the test area; and displaying the position mark on the ultrasonic detection image in an overlapping manner according to the position coordinates.
In a possible implementation manner, the elastic detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color, and the pseudo-color is used for representing corresponding time when the shear wave arrives at different positions when being transmitted inside the tissue to be detected.
In one possible implementation, each of the position markers has a different smoothness, which is used to characterize the propagation stability of the shear wave.
In a possible implementation manner, after obtaining the ultrasound echo signal, the generating module is further configured to: calculating an ultrasonic attenuation parameter UAP according to the formula of (I1-I2)/(D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2; calculating an ultrasonic attenuation coefficient alpha according to the formula of alpha-20 log (I1/I2) -6)/2 (D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2; calculating a scattering son distribution coefficient k according to the formula k which is s/sigma; where s is the coherent signal energy; σ is the spread signal energy; k ranges from [0, 1 ]; after the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient alpha and the scatterer distribution coefficient k are normalized, calculating a liver tissue lesion comprehensive Score according to a formula of Score (a) E + b UAP + c a + d) k according to an ultrasonic echo signal and corresponding histopathological typing; wherein, 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 ].
According to a third aspect of embodiments of the present application, there is provided an electronic device, comprising: a memory, a processor, and a computer program;
wherein the computer program is stored in the memory and configured to be executed by the processor for performing the ultrasound elastography method as defined in any one of the first aspects of the embodiments of the present application.
According to a fourth aspect of embodiments of the present application, there is provided a computer-readable storage medium having stored therein computer-executable instructions for implementing the ultrasound elastography method according to any one of the first aspect of embodiments of the present application when executed by a processor.
According to the ultrasonic elastography method, the ultrasonic elastography device, the electronic equipment and the storage medium, the ultrasonic echo signal is obtained by applying vibration excitation to the tissue to be detected and performing ultrasonic detection; performing beam synthesis according to the ultrasonic echo signals to generate synthesized ultrasonic imaging data; generating an ultrasonic detection image and an elasticity detection image by using the synthesized ultrasonic imaging data, wherein the elasticity detection image is used for representing the tissue elasticity detectability information and the signal quality of the ultrasonic detection image at the corresponding moment; and after the ultrasonic detection image is generated, the generated elastic detection image is superposed and displayed on the ultrasonic detection image, and the signal quality of the ultrasonic detection image is represented by utilizing the elastic detection image, so that an operator can quickly control the quality of the detection result, the detection efficiency is improved, and the detection accuracy is improved.
Drawings
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 application.
Fig. 1 is a scene diagram of an application of the ultrasound elastography method provided in the embodiment of the present application;
FIG. 2 is a flow chart of a method of ultrasound elastography provided in an embodiment of the present application;
FIG. 3 is a schematic flow chart illustrating a process of generating an ultrasonic inspection image and an elastic inspection image according to an embodiment of the present application;
FIG. 4 is a schematic diagram of a superimposed display of an elasticity inspection image and an ultrasonic inspection image according to an embodiment of the present application;
FIG. 5 is a flow chart of an elastography method provided in another embodiment of the present application;
FIG. 6 is a schematic diagram of a process for quality control according to smoothness of position marks according to an embodiment of the present disclosure;
FIG. 7 is a schematic structural diagram of an ultrasound elastography device provided in an embodiment of the present application;
fig. 8 is a schematic diagram of an electronic device according to an embodiment of the present application.
With the above figures, there are shown specific embodiments of the present application, which will be described in more detail below. These drawings and written description are not intended to limit the scope of the inventive concepts in any manner, but rather to illustrate the inventive concepts to those skilled in the art by reference to specific embodiments.
Detailed Description
Reference will now be made in detail to the exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, like numbers in different drawings represent the same or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the present application, as detailed in the appended claims.
The following explains an application scenario of the embodiment of the present application:
fig. 1 is a view of an application scenario of an ultrasound elastography method provided in an embodiment of the present application, and as shown in fig. 1, the ultrasound elastography method provided in the embodiment of the present application is applied to an electronic device, in particular, an elastography device 11 that can be applied clinically. The elastography device 11 performs ultrasonic elastography on a human body by the elastography method provided by the embodiment of the application to obtain mechanical characteristics of the structure of the human visceral tissue 12, and further judges whether the human visceral tissue 12 is in a healthy state, thereby providing important data support for screening early lesions.
In the prior art, the result of the ultrasound elastography is usually displayed to an operator of the apparatus in the form of a quantitative parameter value, and the operator observes the ultrasound detection image and combines the quantitative parameter value of the elastography to realize the state detection of the biological tissue. However, due to the non-uniform characteristics of human tissues, the test position is not good, and human breathing and other reasons can interfere with the ultrasonic detection image, affect the quality of ultrasonic imaging, and affect the judgment of an operator on the tissue state during the process of performing elasticity measurement. Meanwhile, the elastic detection result can represent the propagation process of the shear wave, and the elastic detection result in the interfered state is different from the elastic detection result in the undisturbed state, so that the quality of the ultrasonic detection image can be controlled through the elastic detection result, the signal quality of the ultrasonic detection image is improved, and the judgment accuracy of an operator on the tissue state is improved. Therefore, how to realize quality control on the ultrasonic detection image through the elastic detection result to improve the signal quality of the ultrasonic detection image is a problem which needs to be solved at present.
The following describes the technical solutions of the present application and how to solve the above technical problems with specific embodiments. The following several specific embodiments may be combined with each other, and details of the same or similar concepts or processes may not be repeated in some embodiments. Embodiments of the present application will be described below with reference to the accompanying drawings.
Fig. 2 is a flowchart of an ultrasound elastography method according to an embodiment of the present application, and as shown in fig. 2, the ultrasound elastography method according to the embodiment is applied to an elastography device, and includes the following steps:
step S101, applying vibration excitation to the tissue to be detected and carrying out ultrasonic detection to obtain an ultrasonic echo signal.
Illustratively, an operator applies vibration excitation to the tissue to be detected by operating a detection probe of the elastography device, so as to generate shear waves transmitted from the surface of the tissue to be detected to the inside of the tissue to be detected, the shear waves can generate micro displacement and deformation in the tissue to be detected, and meanwhile, ultrasonic detection is correspondingly applied to the tissue to be detected, so as to acquire data of the micro displacement and deformation generated by the tissue to be detected. Specifically, for example, an ultrasonic wave is transmitted to the tissue to be measured by the detection probe and an echo is received, so as to obtain an ultrasonic echo signal, and the ultrasonic echo signal can be used to observe the micro displacement and deformation of the tissue to be measured.
The execution sequence for performing vibration excitation and ultrasonic detection on the tissue to be detected may include multiple types, for example, first applying vibration excitation and then performing ultrasonic detection, or first performing ultrasonic detection and then applying vibration excitation, or then applying vibration excitation and performing ultrasonic detection simultaneously, which is not specifically limited herein and may be set as required.
Further, when applying vibration excitation to the tissue to be detected and performing ultrasonic detection, the vibration frequencies may be the same or different, and in one possible implementation, the vibration excitation of different vibration frequencies is sequentially applied to the tissue to be detected and the ultrasonic detection is performed. Specifically, for example, the vibration excitation is cyclically applied to the tissue to be detected, and during each cycle, the vibration frequency of the vibration excitation is changed until a preset stop condition is met, and the tissue to be detected is simultaneously subjected to ultrasonic detection, wherein the ultrasonic detection may be performed at the time after the vibration excitation with a specific frequency is applied each time, or after all the vibration excitations are applied, and this time is not particularly limited. The tissue to be detected is subjected to vibration excitation with different frequencies, so that shear waves with different frequencies are generated inside the tissue to be detected, and ultrasonic detection is correspondingly applied to the tissue to be detected after each vibration excitation, so that response detection on the tissue to be detected under different excitation frequencies is realized.
And step S102, performing beam synthesis according to the ultrasonic echo signals to generate synthesized ultrasonic imaging data.
Exemplarily, after transmitting ultrasonic waves to the tissue to be measured and receiving ultrasonic echo signals formed after the ultrasonic waves are blocked by the tissue to be measured, beam forming is performed on the ultrasonic echo signals, and the result of the beam forming is ultrasonic imaging data. In an exemplary embodiment, the vibration excitation applied to the tissue to be measured by the elastography device is a shear wave, the ultrasonic echo signal is an ultrasonic echo signal acquired by the elastography device and used for tracking shear wave transmission, and then the ultrasonic echo signal used for tracking shear wave propagation is subjected to beam forming to generate synthesized ultrasonic imaging data.
And step S103, generating an ultrasonic detection image and an elasticity detection image by using the synthesized ultrasonic imaging data, wherein the elasticity detection image is used for representing the tissue elasticity detectability information and the signal quality of the ultrasonic detection image at the corresponding moment.
Fig. 3 is a schematic flowchart of a process for generating an ultrasound detection image and an elasticity detection image according to an embodiment of the present disclosure, and referring to fig. 3, for example, the synthesized ultrasound imaging data is data carrying deformation information of a tissue to be detected, and the synthesized ultrasound imaging data may be used for performing ultrasound imaging and elasticity imaging simultaneously, that is, generating the ultrasound detection image and the elasticity detection image. Specifically, the step of applying ultrasonic detection to the tissue to be detected includes transmitting ultrasonic waves to the tissue to be detected and receiving echo waves formed after the ultrasonic waves are blocked by the tissue to be detected, and in the cyclic process of vibration excitation-ultrasonic detection, shear waves formed by mechanical vibration excite the tissue to be detected and the tissue to be detected generates micro deformation, so that deformation information of the tissue to be detected is carried in ultrasonic echo signals. And performing beam forming and imaging calculation on the ultrasonic echo signals to obtain corresponding ultrasonic detection images. The specific implementation of generating the ultrasonic detection image according to the ultrasonic echo signal is the prior art in the field, and is not described herein again. Further, by processing and imaging calculation of the ultrasonic echo signals, a corresponding elasticity detection image can be obtained, specifically, for example, first, filtering processing is performed on the ultrasonic echo signals, then, tissue displacement or strain amount is calculated on the filtered ultrasonic echo signals, and a motion propagation mode map is obtained. And determining a dispersion curve according to the motion propagation mode diagram, and determining the phase velocity according to the dispersion curve. And performing function fitting on the phase velocity to obtain the viscoelastic information of the tissue to be detected. The method for obtaining the viscosity parameter by fitting the phase velocities corresponding to different propagation rates is known to those skilled in the art and will not be described herein again. The elasticity detection image is an image representation of viscoelastic information, and may be implemented in the form of a pseudo-color image. Since the generated ultrasonic detection image and the elastic detection image are generated by using the same source data, the detection regions corresponding to the generated ultrasonic detection image and the elastic detection image are exactly the same position, and there is no time error. Compared with the scheme that the elastic imaging data are calculated through the original ultrasonic echo signals and through calculation of the original ultrasonic echo signals, the ultrasonic detection image and the elastic detection image are generated, the time error problem caused by the fact that the original ultrasonic echo signals are inconsistent with the calculated elastic imaging data can be avoided, interference is further avoided, and the reliability of elastic measurement evaluation is improved.
Optionally, in this embodiment, the method further includes: according to the propagation velocity of the shear waveTo modulus of elasticity E, where E ═ 3 ρ V2ρ is the tissue density and V is the shear wave velocity. The shear wave propagation speed is determined by the phase velocity, and the tissue density is obtained by recording preset configuration information of the tissue density to be measured.
Optionally, in this embodiment, after obtaining the ultrasound echo signal, the method further includes:
calculating an ultrasonic attenuation parameter UAP according to the formula of (I1-I2)/(D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2; calculating an ultrasonic attenuation coefficient alpha according to the formula of alpha-20 log (I1/I2) -6)/2 (D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2; calculating a scattering son distribution coefficient k according to the formula k which is s/sigma; where s is the coherent signal energy; σ is the spread signal energy; k ranges from [0, 1 ]; after the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient alpha and the scatterer distribution coefficient k are normalized, calculating a liver tissue lesion comprehensive Score according to a formula of Score (a) E + b UAP + c a + d) k according to an ultrasonic echo signal and corresponding histopathological typing; wherein, 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 ". x" represents a dot-by-dot operation.
Wherein, the elasticity modulus E is used for evaluating the elasticity of the liver and the degree of hepatic fibrosis; evaluating the degree of liver tissue fat change by using an ultrasonic attenuation parameter UAP; evaluating the degree of inflammation of the liver tissue by using an ultrasonic attenuation coefficient alpha; evaluating the degree of the liver tissue bubble sample change by using a scatterer distribution coefficient k; UAP/alpha/k comprehensive evaluation evaluates the degree of the fatty liver; score comprehensive assessment of liver tissue lesion extent; wherein, after normalizing each parameter, the Score is calculated according to the ultrasonic signal and the corresponding tissue pathological type.
In this embodiment, a quantitative ultrasound parameter is extracted according to the ultrasound echo signal, and multi-dimensional evaluation information is provided. Except the elastic detection image, quantitative ultrasonic parameters are simultaneously output and displayed, and the method assists a user in better judging the change of the tissue structure through the optimized ultrasonic detection image, better judging the elastic change of the tissue through the elastic detection image, better judging the change of the microstructure through the quantitative ultrasonic parameters and comprehensively evaluating the tissue state, thereby improving the detection efficiency and the detection accuracy.
The process of generating the ultrasound inspection image and the elasticity inspection image by using the synthesized ultrasound imaging data includes: the synthetic ultrasonic imaging data is processed through the GPU, an ultrasonic detection image and an elastic detection image are generated in parallel, the image processing speed is improved, the detection image can be displayed in real time, and the detection accuracy is improved.
Optionally, after the generation of the ultrasound inspection image and before the display of the ultrasound inspection image, the following processing is further performed on the ultrasound inspection image:
carrying out median filtering to remove electronic noise; performing Gaussian smoothing and non-local mean filtering to smooth the structural information; carrying out bilateral filtering and enhancing edge information; and taking two thresholds of the signal as display thresholds, wherein the lower threshold is the amplitude of the signal from 1% to 20%, and the upper threshold is the amplitude of the signal from 60% to 200%, and performing histogram equalization processing according to the display thresholds.
In the step of this embodiment, before the ultrasonic detection image is displayed, the ultrasonic detection image is optimized, so that the processed ultrasonic detection image is clearer, which is helpful for a user to evaluate the state of the tissue to be detected by observing the ultrasonic detection image, and improves the detection efficiency and the detection accuracy.
And step S104, overlapping and displaying the ultrasonic detection image and the elastic detection image according to the position relation of the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image.
Illustratively, the ultrasound detection image is image information representing a tissue form of a test region inside a tissue to be detected, where the image information includes description information of a spatial position inside the tissue to be detected, the elasticity detection image is image information representing a propagation process of a shear wave in the test region, fig. 4 is a schematic diagram of an overlapped display of the elasticity detection image and the ultrasound detection image provided in an embodiment of the present application, as shown in fig. 4, the elasticity detection image includes a plurality of position markers 31, the position markers 31 are used for representing propagation positions of the shear wave generated by vibration excitation inside the tissue to be detected, and the positions of the shear wave generated by vibration excitation propagating inside the tissue to be detected have a real correspondence with the spatial position of the tissue to be detected displayed by the ultrasound detection image 32. Each position mark has different smoothness, the smoothness is used for representing the propagation stability of the shear wave, and when the smoothness of the position mark 31 is better, the propagation stability of the shear wave is good, namely the propagation of the shear wave is not influenced by factors such as the respiration of a tested human body, the poor testing position and the like, so that the signal quality of the ultrasonic detection image 32 is better, and the reliability is higher; on the contrary, when the smoothness of the position mark 31 is poor, the propagation stability of the characteristic shear wave is poor, that is, the propagation of the shear wave is affected by the respiration of the tested human body, the testing position is not good, and the like, so the signal quality of the ultrasonic detection image 32 is poor, and the reliability is low.
Further, the elastography device includes a display unit, such as a display screen, for example, when the operator performs the elastography on the tissue to be tested, the position of the detection probe needs to be continuously adjusted to locate a better test area, so that the internal tissue state of the tissue to be tested can be better represented in the ultrasonic image. Through this embodiment, detect the image and elasticity with the supersound and detect the image stack and show on the display screen, at this moment, because elasticity detects the image and has played the quality control effect to the supersound and detect the image, consequently, the operator is in the position of constantly adjusting test probe, the test area's of location preferred in-process, can utilize this supersound test image to detect the quality control effect of image to elasticity, whether the current position at place of better definite test probe corresponds the test area of preferred, thereby improve the efficiency of location test area, improve the signal quality of supersound detection image, and improve the testing result accuracy.
Applying vibration excitation to the tissue to be detected and carrying out ultrasonic detection to obtain an ultrasonic echo signal; generating an ultrasonic detection image and an elastic detection image according to the ultrasonic echo signal, wherein the elastic detection image is used for representing 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 displayed in an overlapped mode according to the position relation of 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 utilized and displayed on the ultrasonic detection image in an overlapped mode, and the elastic detection image is utilized to represent the signal quality of the ultrasonic detection image, so that an operator can quickly control the quality of a detection result, the influence caused by the reasons of poor test position, the breathing of a detected human body and the like is eliminated, the signal quality of the ultrasonic detection image is improved, and the accuracy of the detection result is improved.
Fig. 5 is a flowchart of an elastography method provided in another embodiment of the present application, and as shown in fig. 5, the elastography method provided in this embodiment further refines step S104 on the basis of the elastography method provided in the embodiment shown in fig. 2, and adds a step of acquiring configuration information, so that the elastography method provided in this embodiment includes the following steps:
step S201, applying shear waves to the tissue to be measured, and acquiring an ultrasonic echo signal for tracking shear wave propagation.
Step S202, performing beam forming on the ultrasonic echo signals for tracking the shear wave propagation to generate synthetic ultrasonic imaging data, wherein the ultrasonic imaging data comprises multiple frames of scanning lines, and the multiple frames of scanning lines are used for representing ultrasonic detection results at different positions.
Step S203, generating an ultrasonic detection image and an elastic detection image according to the multi-frame scanning lines in the synthetic ultrasonic imaging data, wherein the elastic detection image comprises at least one position mark, and the position mark is used for representing the propagation position of the shear wave generated by vibration excitation in the tissue to be detected.
Illustratively, the ultrasound inspection image is image information generated by beam-forming and imaging calculations for a plurality of frames of scan lines. The elasticity test image is image information representing viscoelastic information of the tissue to be tested, and the process of determining the viscoelastic information of the tissue to be tested includes, for example: filtering the multiframe scanning lines, determining a motion propagation pattern diagram according to echo data corresponding to the filtered multiframe scanning lines, performing time-frequency analysis on the motion propagation pattern diagram, determining a dispersion curve, determining phase velocities corresponding to different frequencies according to the dispersion curve, and performing function fitting on the phase velocities corresponding to the different frequencies to obtain the viscoelastic information of the tissue to be measured. The ultrasonic detection image and the elastic detection image generated by the scanning lines of a plurality of frames have the same source data, so that the detection areas corresponding to the ultrasonic detection image and the elastic detection image are completely the same position without any time error. The time error problem caused by the inconsistency of the original ultrasonic echo signal and the calculated elastic imaging data can be avoided, so that the interference is further avoided, and the reliability of the elastic measurement evaluation is improved.
And step S204, determining a positioning coordinate system by taking the test area of the ultrasonic detection image as a reference.
Furthermore, the 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, a coordinate system for describing the spatial position can be established by taking the test area corresponding to the synthesized ultrasonic imaging data as a reference, and the position coordinates in the coordinate system are used for representing the spatial position of the coordinate point. Further, the position mark in the elastic detection image is also generated by the ultrasonic echo signal, so that the position coordinate of the position mark in the coordinate system taking the test area as the reference can be determined according to the synthetic ultrasonic imaging data.
Step S205, according to the position relation between the position mark and the test area, the position coordinate of the position mark in the positioning coordinate system is determined.
In a possible implementation manner, the range of the test area may be adjusted and set according to a user requirement, for example, obtaining configuration information input by the user, and determining the range of the test area according to the configuration information. The range of the test area is adjusted through the configuration information, only the interested area or the preset range is displayed, the influence of surrounding non-interested areas or other interference information on the judgment of the result is avoided, and the flexibility and the application range of the algorithm can be further improved.
And step S206, displaying the position mark on the ultrasonic detection image in an overlapping mode according to the position coordinates.
Specifically, according to the position coordinates, the position mark is displayed at the position coordinates, and the position mark and the ultrasonic detection image are displayed in a superposition mode. The position markers represent the spatial positions of arrival of the shear waves over time. Optionally, the position marker is a pseudo-color map, different position markers have different colors, and the different colors represent the time when the shear wave reaches the position, e.g., position marker a is red, indicating that the time when the shear wave reaches position marker a is 50 milliseconds, position marker B is yellow, indicating that the time when the shear wave reaches position marker B is 100 milliseconds; position marker C is blue indicating that the time of arrival of the shear wave at position marker C is 150 milliseconds. The position markers have corresponding smoothness according to color variations between the position markers, wherein the smoothness is used for characterizing the propagation stability of the shear wave. For example, the color change between the position markers can be changed according to the color temperature, so that the smoothness of the color change between the position markers can be more conveniently identified by an operator. The following describes how to perform quality control according to the smoothness of the position mark in a more specific embodiment.
Fig. 6 is a schematic diagram of a process of quality control according to smoothness of position marks according to an embodiment of the present application, and as shown in fig. 6, the position marks A, B, C, D, E are respectively displayed in an ultrasound detection image 32 in an overlapping manner, and each position mark 31 represents a position where a shear wave reaches a corresponding human tissue in the ultrasound detection image. When the operator performs the elastic detection by operating the detection probe 111, the position mark 31 is displayed in a different color superimposed on the ultrasonic detection image 32, that is, the position marks A, B, C, D, E have different pseudo-color colors, for example, the position mark a is red, the position mark B is orange, the position mark C is yellow, the position mark D is light blue, and the position mark E is dark blue. Different colors, corresponding to different times of arrival of the shear wave at the location.
When the operator moves the detection probe 111 to the first position, the detection environment is not good, and the shear wave is affected during the propagation process in the human tissue, so that the propagation speed of the shear wave in the human body is not uniform, which may result in uneven color change from the mark position a to the mark position E. The operator can judge that the detection environment is not good at this time by observing the color change of the position mark 31, so that the detection probe 111 is moved to the second position, the detection environment is good at this time, that is, the shear wave is not influenced by human respiration, uneven tissue blockage and the like in the propagation process in the human tissue, and 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 smoothly changed. An operator can determine whether the interference is caused in the test area by observing the smoothness of the pseudo-color image type position mark 31, so that the quality control purpose of the ultrasonic detection image is realized.
In the present embodiment, the signal quality of the ultrasonic detection image is expressed by displaying the marker positions in the form of pseudo-color images and by the smoothness of each position marker. The elastic detection image superposed in the ultrasonic detection image can improve the quality control of the ultrasonic detection image. Compared with the technical scheme that the elastic detection result is displayed through the character identification in the prior art, the display mode of the pseudo-color image enables an operator to observe the signal quality of the ultrasonic detection image more intuitively and in real time, so that the test area is adjusted, the interference on the elastic detection is avoided, and the accuracy of the detection of the tissue state to be detected is improved.
In this embodiment, the implementation manners of step S201 to step S202 are already described in the embodiment shown in fig. 2 of the present application, and are not described in detail here.
Fig. 7 is a schematic structural diagram of an ultrasound elastography device provided in an embodiment of the present application, and is applied to an elastography apparatus, as shown in fig. 7, an ultrasound elastography device 3 provided in this embodiment includes:
the detection module 31 is configured to apply vibration excitation to a tissue to be detected and perform ultrasonic detection to obtain an ultrasonic echo signal;
a synthesis module 32, configured to perform beam synthesis according to the ultrasound echo signal, and generate synthesized ultrasound imaging data;
a generating module 33, configured to generate an ultrasound detection image and an elasticity detection image by using the synthesized ultrasound imaging data, where the elasticity detection image is used to represent tissue elasticity detectability information and signal quality of the ultrasound detection image at a corresponding time;
and the display module 34 is configured to display the ultrasonic detection image and the elastic detection image in an overlapping manner according to the position relationship between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image.
In a possible implementation manner, the detection module 31 is specifically configured to: applying shear waves to the tissue to be detected, and acquiring ultrasonic echo signals for tracking shear wave propagation; the synthesis module 32 is specifically configured to: and performing beam synthesis on the ultrasonic echo signals for tracking the shear wave propagation to generate multi-frame scanning lines, wherein the multi-frame scanning lines are used for representing ultrasonic detection results at different positions.
In one possible implementation, the generating module 33 has a module for: processing the synthesized ultrasonic imaging data through a graphic processor, generating an ultrasonic detection image and an elastic detection image in parallel, and obtaining an elastic modulus E according to the shear wave propagation speed, wherein E is 3 rho V2ρ is the tissue density and V is the shear wave velocity.
In a possible implementation manner, after generating the ultrasound inspection image, the generating module 33 further performs the following processing on the ultrasound inspection image: carrying out median filtering to remove electronic noise; performing Gaussian smoothing and non-local mean filtering to smooth the structural information; carrying out bilateral filtering and enhancing edge information; and taking two thresholds of the signal as display thresholds, wherein the lower threshold is the amplitude of the signal from 1% to 20%, and the upper threshold is the amplitude of the signal from 60% to 200%, and performing histogram equalization processing according to the display thresholds.
In a possible implementation manner, the elasticity detection image includes at least one position mark, the position mark is used for characterizing a propagation position of the shear wave generated by the vibration excitation in the tissue to be measured, and the display module 34 is specifically configured to: determining a positioning coordinate system by taking a test area of the ultrasonic detection image as a reference; determining the position coordinates of the position marks in a positioning coordinate system according to the position relation between the position marks and the test area; and displaying the position mark on the ultrasonic detection image in an overlapping manner according to the position coordinates.
In a possible implementation manner, the elastic detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color, and the pseudo-color is used for representing corresponding time when the shear wave reaches different positions when being transmitted inside the tissue to be detected.
In one possible implementation, each position marker has a different smoothness, which is used to characterize the propagation stability of the shear wave.
In a possible implementation manner, after obtaining the ultrasound echo signal, the generating module 33 is further configured to: calculating an ultrasonic attenuation parameter UAP according to the formula of (I1-I2)/(D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2; calculating an ultrasonic attenuation coefficient alpha according to the formula of alpha-20 log (I1/I2) -6)/2 (D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2; calculating a scattering son distribution coefficient k according to the formula k which is s/sigma; where s is the coherent signal energy; σ is the spread signal energy; k ranges from [0, 1 ]; after the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient alpha and the scatterer distribution coefficient k are normalized, calculating a liver tissue lesion comprehensive Score according to a formula of Score (a) E + b UAP + c a + d) k according to an ultrasonic echo signal and corresponding histopathological typing; wherein, 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 detection module 31, the synthesis module 32, the generation module 33 and the display module 34 are connected in sequence. The ultrasound elastography device 3 provided in this embodiment may implement the technical solution of the method embodiment as shown in any one of fig. 2 to 6, and the implementation principle and the technical effect are similar, and are not described herein again.
Fig. 8 is a schematic diagram of an electronic device according to an embodiment of the present application, and as shown in fig. 8, an electronic device 4 according to the embodiment includes: a memory 41, a processor 42 and a computer program.
Wherein the computer program is stored in the memory 41 and configured to be executed by the processor 42 to implement the ultrasound elastography method provided by any of the embodiments corresponding to fig. 2-5 of the present application.
The memory 41 and the processor 42 are connected by a bus 43.
The relevant descriptions and effects corresponding to the steps in the embodiments corresponding to fig. 2 to fig. 6 can be understood, and are not described in detail herein.
In one possible implementation, the electrons are elastography devices.
One embodiment of the present application provides a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to implement the ultrasound elastography method provided in any one of the embodiments corresponding to fig. 2 to fig. 6 of the present application.
The computer readable storage medium may be, among others, ROM, Random Access Memory (RAM), CD-ROM, magnetic tape, floppy disk, optical data storage device, and the like.
An embodiment of the present application provides a computer program product, which includes a computer program, and the computer program is used for implementing the ultrasound elastography method provided in any one of the embodiments corresponding to fig. 2 to fig. 6 of the present application when being executed by a processor.
In the several embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, a division of modules is merely a division of logical functions, and an actual implementation may have another division, for example, a plurality of modules or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or modules, and may be in an electrical, mechanical or other form.
Other embodiments of the present application will be apparent to those skilled in the art from consideration of the specification and practice of the application disclosed herein. This application is intended to cover any variations, uses, or adaptations of the invention following, in general, the principles of the application and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the application being indicated by the following claims.
It will be understood that the present application is not limited to the precise arrangements 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 application is limited only by the appended claims.
Claims (19)
1. A method of ultrasound elastography, the method comprising:
applying vibration excitation to the tissue to be detected and carrying out ultrasonic detection to obtain an ultrasonic echo signal;
performing beam synthesis according to the ultrasonic echo signals to generate synthesized ultrasonic imaging data;
generating an ultrasonic detection image and an elasticity detection image by using the synthesized ultrasonic imaging data, wherein the elasticity detection image is used for representing the tissue elasticity detectability information and the signal quality of an ultrasonic detection signal at a corresponding moment;
and displaying the ultrasonic detection image and the elastic detection image in an overlapping manner according to the position relation of the test region corresponding to the ultrasonic detection image and the test region corresponding to the elastic detection image.
2. The method of claim 1, wherein applying a vibrational excitation to the tissue under test and performing an ultrasonic probe to obtain an ultrasonic echo signal comprises:
applying shear waves to the tissue to be detected, and acquiring ultrasonic echo signals for tracking shear wave propagation;
performing beam synthesis according to the ultrasonic echo signal to generate synthesized ultrasonic imaging data, including:
and performing beam synthesis on the ultrasonic echo signals for tracking the shear wave propagation to generate multi-frame scanning lines, wherein the multi-frame scanning lines are used for representing ultrasonic detection results at different positions.
3. The method of claim 1, wherein generating an ultrasound inspection image and an elasticity inspection image using the composite ultrasound imaging data comprises:
processing the synthesized ultrasound imaging data by a graphics processor, generating the ultrasound detection image and the elasticity detection image in parallel, and obtaining an elastic modulus E according to a shear wave propagation velocity, wherein E is 3 rho V2ρ is the tissue density and V is the shear wave velocity.
4. The method of claim 3, further comprising, prior to displaying after the ultrasound inspection image is generated, processing the ultrasound inspection image by:
carrying out median filtering to remove electronic noise;
performing Gaussian smoothing and non-local mean filtering to smooth the structural information;
carrying out bilateral filtering and enhancing edge information;
and taking two thresholds of the signal as display thresholds, wherein the lower threshold is the amplitude of the signal from 1% to 20%, and the upper threshold is the amplitude of the signal from 60% to 200%, and performing histogram equalization processing according to the display thresholds.
5. The method according to claim 1, wherein the elastic detection image includes at least one position mark, the position mark is used for representing a propagation position of the shear wave generated by the vibration excitation in the tissue to be detected, and the ultrasonic detection image and the elastic detection image are displayed in an overlapping manner according to a position relationship between a test area corresponding to the ultrasonic detection image and a test area corresponding to the elastic detection image, including:
determining a positioning coordinate system by taking the test area of the ultrasonic detection image as a reference;
determining the position coordinates of the position marks in the positioning coordinate system according to the position relation between the position marks and the test area;
and displaying the position mark on the ultrasonic detection image in an overlapping manner according to the position coordinates.
6. The method according to claim 5, wherein the elastic detection image is a pseudo-color image, each position mark corresponds to a different pseudo-color, and the pseudo-color colors are used for representing corresponding time when the shear wave arrives at different positions when the shear wave is transmitted inside the tissue to be detected.
7. The method of claim 6, wherein each of the position markers has a different smoothness, the smoothness being used to characterize the propagation stability of the shear wave.
8. The method of claim 3, further comprising, after obtaining the ultrasound echo signal:
calculating an ultrasonic attenuation parameter UAP according to the formula of (I1-I2)/(D1-D2);
calculating an ultrasonic attenuation coefficient alpha according to the formula of alpha-20 log (I1/I2) -6)/2 (D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2;
calculating a scattering son distribution coefficient k according to the formula k which is s/sigma; where s is the coherent signal energy; σ is the spread signal energy; k ranges from [0, 1 ];
after the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient alpha and the scatterer distribution coefficient k are normalized, calculating a liver tissue lesion comprehensive Score according to a formula of Score (a) E + b UAP + c a + d) k according to an ultrasonic echo signal and corresponding histopathological typing; wherein, 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 ].
9. An ultrasound elastography device, characterized in that the device comprises:
the detection module is used for applying vibration excitation to the tissue to be detected and carrying out ultrasonic detection to obtain an ultrasonic echo signal;
the synthesis module is used for carrying out beam synthesis according to the ultrasonic echo signals to generate synthesized ultrasonic imaging data;
a generating module, configured to generate an ultrasound detection image and an elasticity detection image by using the synthesized ultrasound imaging data, where the elasticity detection image is used to represent tissue elasticity detectability information and signal quality of the ultrasound detection image at a corresponding time;
and the display module is used for displaying the ultrasonic detection image and the elastic detection image in an overlapping manner according to the position relation between the test area corresponding to the ultrasonic detection image and the test area corresponding to the elastic detection image.
10. The apparatus according to claim 9, wherein the detection module is specifically configured to:
applying shear waves to the tissue to be detected, and acquiring ultrasonic echo signals for tracking shear wave propagation;
the synthesis module is specifically configured to: and performing beam synthesis on the ultrasonic echo signals for tracking the shear wave propagation to generate multi-frame scanning lines, wherein the multi-frame scanning lines are used for representing ultrasonic detection results at different positions.
11. The device of claim 9, whichCharacterized in that the generating module has means for: processing the synthesized ultrasound imaging data by a graphics processor, generating the ultrasound detection image and the elasticity detection image in parallel, and obtaining an elastic modulus E according to a shear wave propagation velocity, wherein E is 3 rho V2ρ is the tissue density and V is the shear wave velocity.
12. The apparatus of claim 11, wherein the generating module, after generating the ultrasound inspection image, further processes the ultrasound inspection image by:
carrying out median filtering to remove electronic noise;
performing Gaussian smoothing and non-local mean filtering to smooth the structural information;
carrying out bilateral filtering and enhancing edge information;
and taking two thresholds of the signal as display thresholds, wherein the lower threshold is the amplitude of the signal from 1% to 20%, and the upper threshold is the amplitude of the signal from 60% to 200%, and performing histogram equalization processing according to the display thresholds.
13. The device according to claim 9, wherein the elasticity detection image comprises at least one position marker for characterizing a propagation position of the shear wave generated by the vibration excitation inside the tissue to be measured, and the display module is specifically configured to:
determining a positioning coordinate system by taking the test area of the ultrasonic detection image as a reference;
determining the position coordinates of the position marks in the positioning coordinate system according to the position relation between the position marks and the test area;
and displaying the position mark on the ultrasonic detection image in an overlapping manner according to the position coordinates.
14. The device according to claim 13, wherein the elastic detection image is a pseudo-color image, each of the position marks corresponds to a different pseudo-color, and the pseudo-color colors are used for representing corresponding times when the shear waves arrive at different positions while being transmitted inside the tissue to be detected.
15. The apparatus of claim 14, wherein each of the position markers has a different smoothness, the smoothness being used to characterize the propagation stability of the shear wave.
16. The apparatus of claim 9, wherein the generating module, after obtaining the ultrasound echo signal, is further configured to:
calculating an ultrasonic attenuation parameter UAP according to the formula of (I1-I2)/(D1-D2);
calculating an ultrasonic attenuation coefficient alpha according to the formula of alpha-20 log (I1/I2) -6)/2 (D1-D2); wherein I1 is the ultrasonic signal intensity of the first interested position D1, and D1 is less than or equal to 5 cm; i2 is divided into the ultrasonic signal intensity of a second interested position D2, D2 is less than or equal to 30 cm; d1< D2;
calculating a scattering son distribution coefficient k according to the formula k which is s/sigma; where s is the coherent signal energy; σ is the spread signal energy; k ranges from [0, 1 ];
after the elastic modulus E, the ultrasonic attenuation parameter UAP, the ultrasonic attenuation coefficient alpha and the scatterer distribution coefficient k are normalized, calculating a liver tissue lesion comprehensive Score according to a formula of Score (a) E + b UAP + c a + d) k according to an ultrasonic echo signal and corresponding histopathological typing; wherein, 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 for displaying the obtained liver tissue lesion comprehensive Score.
17. An electronic device, comprising: a memory, a processor, and a computer program;
wherein the computer program is stored in the memory and configured to be executed by the processor to implement the ultrasound elastography method as claimed in any of claims 1 to 8.
18. A computer-readable storage medium having computer-executable instructions stored thereon which, when executed by a processor, are configured to implement the ultrasound elastography method of any of claims 1 to 8.
19. A computer program product comprising a computer program which, when executed by a processor, implements the ultrasound elastography method of any of claims 1 to 8.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110369497.4A CN113040816A (en) | 2021-04-06 | 2021-04-06 | Ultrasonic elastography method, device, electronic equipment and storage medium |
| PCT/CN2022/085128 WO2022213938A1 (en) | 2021-04-06 | 2022-04-02 | Ultrasound elastography method and apparatus, and electronic device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110369497.4A CN113040816A (en) | 2021-04-06 | 2021-04-06 | Ultrasonic elastography method, device, electronic equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113040816A true CN113040816A (en) | 2021-06-29 |
Family
ID=76517683
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110369497.4A Pending CN113040816A (en) | 2021-04-06 | 2021-04-06 | Ultrasonic elastography method, device, electronic equipment and storage medium |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN113040816A (en) |
| WO (1) | WO2022213938A1 (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022213938A1 (en) * | 2021-04-06 | 2022-10-13 | 无锡海斯凯尔医学技术有限公司 | Ultrasound elastography method and apparatus, and electronic device and storage medium |
| CN116549008A (en) * | 2023-05-30 | 2023-08-08 | 无锡海斯凯尔医学技术有限公司 | Elastic imaging system based on ultrasonic signals and application method |
| CN116671999A (en) * | 2023-05-30 | 2023-09-01 | 无锡海斯凯尔医学技术有限公司 | Tissue evaluation method, device, electronic equipment and storage medium |
| CN116862782A (en) * | 2023-05-30 | 2023-10-10 | 无锡海斯凯尔医学技术有限公司 | Image optimization method, device, electronic equipment and storage medium |
| WO2024001142A1 (en) * | 2022-06-30 | 2024-01-04 | 无锡海斯凯尔医学技术有限公司 | Computing device for assessing liver lesion conditions, liver elasticity measurement device, remote workstation, and medium |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150094579A1 (en) * | 2013-09-30 | 2015-04-02 | Siemens Medical Solutions Usa, Inc. | Shear Wave Detection in Medical Ultrasound Imaging |
| CN104825195A (en) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | Shear wave viscoelasticity imaging method and system |
| DE102017211895A1 (en) * | 2016-07-13 | 2018-01-18 | Siemens Medical Solutions Usa, Inc. | Tissue characterization in medical diagnostic ultrasound |
| CN109875608A (en) * | 2019-01-30 | 2019-06-14 | 清华大学 | Elastograph imaging method |
| CN110477954A (en) * | 2019-07-08 | 2019-11-22 | 无锡海斯凯尔医学技术有限公司 | Detection device based on elastogram |
| CN110613484A (en) * | 2019-09-26 | 2019-12-27 | 无锡海斯凯尔医学技术有限公司 | Tissue elasticity detection method and equipment |
| US20200205786A1 (en) * | 2017-04-06 | 2020-07-02 | Siemens Medical Solutions Usa, Inc. | Liver Disease Activity Estimation with Ultrasound Medical Imaging |
| CN111772677A (en) * | 2020-07-07 | 2020-10-16 | 意领科技有限公司 | Biological tissue elasticity detection method with dimension, detection system and storage medium |
| US20200390421A1 (en) * | 2019-06-14 | 2020-12-17 | Echosens | Method and device for measuring an ultrasound parameter of a viscoelastic medium |
| CN113081038A (en) * | 2021-04-06 | 2021-07-09 | 无锡海斯凯尔医学技术有限公司 | Elasticity imaging method and device, electronic equipment and storage medium |
| WO2022213938A1 (en) * | 2021-04-06 | 2022-10-13 | 无锡海斯凯尔医学技术有限公司 | Ultrasound elastography method and apparatus, and electronic device and storage medium |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103845074B (en) * | 2012-11-28 | 2017-12-12 | 深圳迈瑞生物医疗电子股份有限公司 | A kind of ultrasonic elastograph imaging system and method |
| CN112534468A (en) * | 2018-08-24 | 2021-03-19 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic elastography device and method for processing elastic image |
-
2021
- 2021-04-06 CN CN202110369497.4A patent/CN113040816A/en active Pending
-
2022
- 2022-04-02 WO PCT/CN2022/085128 patent/WO2022213938A1/en active Application Filing
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150094579A1 (en) * | 2013-09-30 | 2015-04-02 | Siemens Medical Solutions Usa, Inc. | Shear Wave Detection in Medical Ultrasound Imaging |
| CN104825195A (en) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | Shear wave viscoelasticity imaging method and system |
| DE102017211895A1 (en) * | 2016-07-13 | 2018-01-18 | Siemens Medical Solutions Usa, Inc. | Tissue characterization in medical diagnostic ultrasound |
| US20200205786A1 (en) * | 2017-04-06 | 2020-07-02 | Siemens Medical Solutions Usa, Inc. | Liver Disease Activity Estimation with Ultrasound Medical Imaging |
| CN109875608A (en) * | 2019-01-30 | 2019-06-14 | 清华大学 | Elastograph imaging method |
| US20200390421A1 (en) * | 2019-06-14 | 2020-12-17 | Echosens | Method and device for measuring an ultrasound parameter of a viscoelastic medium |
| CN110477954A (en) * | 2019-07-08 | 2019-11-22 | 无锡海斯凯尔医学技术有限公司 | Detection device based on elastogram |
| CN110613484A (en) * | 2019-09-26 | 2019-12-27 | 无锡海斯凯尔医学技术有限公司 | Tissue elasticity detection method and equipment |
| CN111772677A (en) * | 2020-07-07 | 2020-10-16 | 意领科技有限公司 | Biological tissue elasticity detection method with dimension, detection system and storage medium |
| CN113081038A (en) * | 2021-04-06 | 2021-07-09 | 无锡海斯凯尔医学技术有限公司 | Elasticity imaging method and device, electronic equipment and storage medium |
| WO2022213938A1 (en) * | 2021-04-06 | 2022-10-13 | 无锡海斯凯尔医学技术有限公司 | Ultrasound elastography method and apparatus, and electronic device and storage medium |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022213938A1 (en) * | 2021-04-06 | 2022-10-13 | 无锡海斯凯尔医学技术有限公司 | Ultrasound elastography method and apparatus, and electronic device and storage medium |
| WO2024001142A1 (en) * | 2022-06-30 | 2024-01-04 | 无锡海斯凯尔医学技术有限公司 | Computing device for assessing liver lesion conditions, liver elasticity measurement device, remote workstation, and medium |
| CN116549008A (en) * | 2023-05-30 | 2023-08-08 | 无锡海斯凯尔医学技术有限公司 | Elastic imaging system based on ultrasonic signals and application method |
| CN116671999A (en) * | 2023-05-30 | 2023-09-01 | 无锡海斯凯尔医学技术有限公司 | Tissue evaluation method, device, electronic equipment and storage medium |
| CN116862782A (en) * | 2023-05-30 | 2023-10-10 | 无锡海斯凯尔医学技术有限公司 | Image optimization method, device, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022213938A1 (en) | 2022-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2022213938A1 (en) | Ultrasound elastography method and apparatus, and electronic device and storage medium | |
| US20230243966A1 (en) | Imaging methods and apparatuses for performing shear wave elastography imaging | |
| US11344278B2 (en) | Ovarian follicle count and size determination using transvaginal ultrasound scans | |
| US6884216B2 (en) | Ultrasound diagnosis apparatus and ultrasound image display method and apparatus | |
| US10743845B2 (en) | Ultrasound diagnostic apparatus and method for distinguishing a low signal/noise area in an ultrasound image | |
| CN113081038B (en) | Elastography method, elastography device, electronic equipment and storage medium | |
| EP2030572B1 (en) | Ultrasonographic device | |
| WO2018023336A1 (en) | Method and system for displaying ultrasonic elastic measurement | |
| EP3563769A1 (en) | Method and ultrasound system for shear wave elasticity imaging | |
| CN105246415B (en) | The method of operating of ultrasound observation apparatus and ultrasound observation apparatus | |
| US20130046168A1 (en) | Method and system of characterization of carotid plaque | |
| US20090024037A1 (en) | Ultrasonic imaging apparatus and a method of acquiring ultrasonic images | |
| US20060058674A1 (en) | Optimizing ultrasound acquisition based on ultrasound-located landmarks | |
| JP6303912B2 (en) | Ultrasonic diagnostic apparatus, ultrasonic diagnostic method, and computer-readable non-transitory recording medium storing program | |
| EP2208465A1 (en) | Ultrasonographic device, operation method thereof, and ultrasonogram diagnosis program | |
| JP2006507883A (en) | A segmentation tool to identify flow areas within an imaging system | |
| CN107427278A (en) | Pass through the elastic method of ultrasound quantification material | |
| KR20140020486A (en) | Method and apparatus for analyzing elastography of tissue using ultrasound | |
| US9855025B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus | |
| US20130182120A1 (en) | Ultrasonic diagnostic apparatus, image processing device, and stain image generating method | |
| JP2020103883A (en) | Ultrasound imaging system and method for displaying target object quality level | |
| CN107833218B (en) | Biological tissue shear wave elastic image quality evaluation display method | |
| Lu et al. | A real time displacement estimation algorithm for ultrasound elastography | |
| US11540809B2 (en) | Ultrasonic diagnostic apparatus and method for propagation speed analysis of shear wave and elastic modulus measurement of a tissue | |
| CN116671996A (en) | Tissue detection method and electronic equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |