CN107510474B - Shear wave elastic imaging method and system - Google Patents
Shear wave elastic imaging method and system Download PDFInfo
- Publication number
- CN107510474B CN107510474B CN201710859560.6A CN201710859560A CN107510474B CN 107510474 B CN107510474 B CN 107510474B CN 201710859560 A CN201710859560 A CN 201710859560A CN 107510474 B CN107510474 B CN 107510474B
- Authority
- CN
- China
- Prior art keywords
- shear wave
- information
- shear
- data
- amplitude
- 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.)
- Active
Links
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/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (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)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention discloses a shear wave elastic imaging method and system capable of reflecting shear wave speed and attenuation information, which are used for measuring and displaying speed information and attenuation information in a shear wave propagation process, and solving the technical problem that in the prior art, the shear wave elastic imaging system has a single measuring object, so that the tissue information is not completely reflected, and are beneficial to improving the reliability and the comprehensiveness of shear wave detection. The invention can be widely applied to various shear wave elastic imaging systems.
Description
Technical Field
The invention relates to the technical field of ultrasonic imaging, in particular to a shear wave elastography method and a shear wave elastography system.
Background
The two-dimensional shear wave elastography technology is mainly characterized in that the propagation velocity of shear waves in tissues is measured, then the shear wave velocity is directly mapped into different colors to be displayed on an image according to a specific color table, and the different velocity values are correspondingly mapped; or according to a formula, combining other parameters such as tissue density and the like, converting the shear wave speed into Young modulus, and then mapping corresponding different Young modulus values into different colors through a specific color table to be displayed on an image. A doctor or an operator can intuitively judge the hardness information of the tissue of the region of interest according to the image color so as to determine the abnormal region.
However, the shear wave velocity and the young's modulus calculated from the shear wave velocity are only a part of the information of the tissue reflected in the shear wave propagation process, and only the information of the hardness of the tissue is reflected, and the reliability of shear wave detection is low.
Disclosure of Invention
In order to better develop the application of shear waves in the field of tissue detection, the invention aims to provide a shear wave elastic imaging method and a shear wave elastic imaging system, which can reflect the speed and attenuation information of shear waves and are beneficial to improving the reliability and comprehensiveness of shear wave detection.
The technical scheme adopted by the invention is as follows:
a method of shear wave elastography comprising the steps of:
s1, exciting to generate shear waves at a preset position by using acoustic radiation force;
s2, transmitting ultrasonic detection wave beams for multiple times on the shear wave propagation path to detect shear waves so as to obtain shear wave detection data;
s3, analyzing the shear wave detection data to obtain the velocity information and attenuation information of the shear wave;
and S4, processing the shear wave velocity information and the attenuation information by using a preset image mapping algorithm, and outputting corresponding image information.
Preferably, the step S3 specifically includes: the method comprises the steps of generating shear wave waveform data by utilizing shear wave detection data, and calculating speed information and attenuation information of shear waves according to the shear wave waveform data of adjacent shear wave detection scanning lines, wherein the attenuation information comprises shear wave amplitude attenuation information and/or shear wave frequency attenuation information and/or shear wave time domain waveform change information.
Preferably, the step S3 includes the sub-step of obtaining shear wave velocity information: and calculating the shear wave velocity information of the position according to the propagation arrival time of the shear wave of the adjacent shear wave detection scanning line position.
Preferably, the step S3 includes the sub-step of obtaining shear wave amplitude attenuation information: determining the peak value of the axial displacement data or the axial velocity data of the shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave amplitude attenuation information.
Preferably, the step S3 includes the sub-step of obtaining shear wave frequency attenuation information: the method comprises the steps of processing shear wave waveform data by utilizing Fourier transform, determining the peak value of axial displacement data or axial velocity data after Fourier transform, setting an amplitude cut-off region, obtaining the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to obtain shear wave frequency attenuation information.
Preferably, the step S2 includes the sub-step of obtaining the change information of the shear wave time domain waveform: determining the peak value of axial displacement data or axial velocity data of shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the squares of the amplitudes in the amplitude cut-off region, and comparing the average values of the squares of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave time domain waveform change information.
Preferably, the step S1 specifically includes: in the region of interest of shear wave elastography, shear waves are generated by utilizing acoustic radiation force excitation at a preset position, and the times of ultrasonic acoustic radiation force excitation are determined according to the size of the region of interest and the propagation distance of the shear waves generated by single excitation.
Preferably, the step S2 specifically includes: the method includes the steps of transmitting ultrasonic detection beams multiple times on a shear wave propagation path to obtain echo signals, processing the echo signals with a beam forming algorithm to generate a plurality of shear wave detection scanlines, and generating shear wave detection data at multiple times on each shear wave detection scanline.
A shear wave elastography system, comprising: the shear wave excitation module is used for exciting at a preset position by using acoustic radiation force to generate shear waves; the detection module is used for transmitting ultrasonic detection wave beams on the shear wave propagation path to detect shear waves so as to acquire shear wave detection data; the calculation module is used for analyzing the shear wave detection data to acquire the speed information and the attenuation information of the shear wave; and the image processing module is used for processing the shear wave velocity information and the attenuation information by utilizing a preset image mapping algorithm and outputting corresponding image information.
Preferably, the attenuation information includes shear wave amplitude attenuation information and/or shear wave frequency attenuation information and/or shear wave time domain waveform change information, and the calculation module specifically includes: the shear wave velocity information calculation submodule is used for detecting the propagation arrival time of the shear wave at the scanning line position according to the adjacent shear waves and calculating the shear wave velocity information of the position; the shear wave amplitude attenuation information submodule is used for determining the peak value of axial displacement data or axial speed data of shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave amplitude attenuation information; the shear wave frequency attenuation information submodule is used for determining the peak value of axial displacement data or axial velocity data after Fourier transform, setting an amplitude cut-off region, acquiring the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave frequency attenuation information; and the shear wave time domain waveform change information submodule is used for determining the peak value of the axial displacement data or the axial speed data of the shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the square of the amplitude in the amplitude cut-off region, and comparing the average values of the squares of the amplitudes of adjacent shear wave detection scanning lines to acquire the shear wave time domain waveform change information.
The invention has the beneficial effects that:
the shear wave elastic imaging method and the shear wave elastic imaging system can be used for measuring the speed information and the attenuation information in the shear wave propagation process and displaying the speed information and the attenuation information, so that the technical problem that the tissue information is not completely reflected due to the fact that a shear wave elastic imaging system is single in measuring object in the prior art is solved, the shear wave elastic imaging method and the shear wave elastic imaging system can be used for reflecting the speed and the attenuation information of the shear wave, and the reliability and the comprehensiveness of shear wave detection are improved.
In addition, the shear wave velocity information and the shear wave attenuation information are calculated by adopting the axial velocity or axial displacement data of the shear wave obtained by carrying out algorithms such as autocorrelation, cross correlation and the like on the shear wave detection data after beam synthesis, so that the method can be conveniently implemented in a system.
The invention can be widely applied to various shear wave elastic imaging systems.
Drawings
The following further describes embodiments of the present invention with reference to the accompanying drawings:
FIG. 1 is a schematic flow chart diagram of one embodiment of a method of shear wave elastography of the present invention;
FIG. 2 is a schematic flow chart diagram of another embodiment of a method of shear wave elastography according to the present invention;
FIG. 3 is a schematic diagram of one embodiment of a method of shear wave elastography of the present invention;
FIG. 4 is a schematic time domain waveform representation of one embodiment of a method of shear wave elastography of the present invention;
FIG. 5 is a schematic frequency domain waveform of an embodiment of a method of shear wave elastography according to the present invention.
Detailed Description
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict.
Biological histological studies showed that: viscoelasticity is the best index to reflect human tissue mechanics information. That is, viscosity needs to be taken into account in addition to elasticity, and it is not sufficient to evaluate the mechanical properties of the tissue as a whole with only a single shear wave velocity information. Meanwhile, most tissues in the human body have anisotropy and dynamics, such as muscles, liver and the like. Other parameters are therefore required to assess the mechanical properties of the tissue. The amplitude, frequency and waveform of the shear wave are changed during the propagation process, and the changes can be used for measuring the reliability and accuracy of the shear wave velocity estimation and as an auxiliary method for evaluating the mechanical properties of the tissues.
Based on the principle, the basic concept of the invention is as follows: the shear wave elastography method and the system thereof are provided, the velocity information and the attenuation information of the shear wave are mapped into different colors or gray scales through a specific image mapping algorithm and are displayed on a display medium or a device. By using the method and the system, not only can the shear wave velocity information of the tissue, namely the shear wave velocity or Young modulus of the tissue be displayed, but also the attenuation information of the shear wave propagation of the tissue can be displayed. The implementation steps are shown in fig. 1 and comprise the following steps:
s1, exciting to generate shear waves at a preset position by using acoustic radiation force;
s2, transmitting ultrasonic detection wave beams for multiple times on the shear wave propagation path to detect shear waves so as to obtain shear wave detection data;
s3, analyzing the shear wave detection data to obtain the velocity information and attenuation information of the shear wave;
and S4, processing the shear wave velocity information and the attenuation information by using a preset image mapping algorithm, and outputting corresponding image information.
It should be noted that the shear wave elastic imaging method and system provided by the present invention are suitable for single-angle or multi-angle plane wave imaging of shear wave detection, and are also suitable for most ultrasonic imaging systems with multi-beam synthesis capability.
Specifically, as shown in fig. 2 to 4, in the present embodiment, the shear wave elastography method includes the following steps:
a1, generating shear waves by exciting with acoustic radiation force at specific preset positions in a Region Of Interest (ROI box) Of shear wave elastography.
As shown in fig. 2, the present implementation calculates a shear wave information matrix from axial displacement or axial velocity data. The shear wave information matrix comprises shear wave speed or Young modulus, shear wave amplitude attenuation, shear wave frequency attenuation, shear wave time domain waveform change and the like.
In particular, the number of ultrasonic acoustic radiation force excitations may be determined depending on the size of the ROI-box and the approximate propagation distance of the shear wave generated by a single excitation. If the length and width of the ROI are twice the approximate distance traveled by a shear wave generated by a single excitation, at least four excitations are required to fill the entire ROI-box with shear waves. Among the methods of generating acoustic radiation is to emit high energy (typically hundreds of microseconds) ultrasonic pulses and then emit tracking pulses to detect the shear waves generated. The emission and detection can be realized on the same probe, and can also be realized by different probes respectively.
A2, transmitting ultrasonic detection wave beams for multiple times on the propagation path of the shear wave to detect the shear wave, generating shear wave detection scanning lines by using a wave beam synthesis algorithm on echo signals, and generating shear wave detection data of multiple time instants on each shear wave detection scanning line.
Specifically, the ultrasonic detection beam with a large coverage can be repeatedly transmitted for multiple times on the propagation path of the shear wave, and this type of ultrasonic beam can be a non-focused plane wave ultrasonic beam (single angle or multi-angle), and can also be a focused ultrasonic beam. Beam synthesis can be roughly classified into software beam synthesis (ultra-fast imaging technology) and hardware beam synthesis (multi-beam synthesis technology). Software beamforming can reconstruct data of all channels by single transmission, while hardware beamforming can only form a few scan lines by single transmission. By using a multi-beam synthesis technology, echo information of a single-emission shear wave detection ultrasonic pulse beam is obtained, and data of a plurality of (usually 4-8) scanning lines are reconstructed. The emitting times and the repetition frequency of the ultrasonic pulse wave beams for detecting the shear waves are determined according to the length width of the ROI box and the shear wave velocity value. For example: the larger the length or width of the ROI box is, the more the emission times are; the larger the shear wave velocity value, the corresponding repetition frequency and the number of transmissions of the detection pulse are correspondingly increased. As shown in fig. 3, wherein an ultrasonic transducer 5 generates an acoustic radiation force excitation ultrasonic pulse beam 1, the acoustic radiation force excitation ultrasonic pulse beam 1 is focused at a preset focal point 2, thereby generating a shear wave within an ROI box 6, 4 is a propagation wavefront of the shear wave, the ultrasonic transducer 5 generates an ultrasonic pulse beam 3 for detecting the shear wave, and 4 shear wave detection scan lines A, B, C and D are generated using a beam forming algorithm on an echo signal of the ultrasonic pulse beam 3. The ultrasound pulse beam for detecting the shear wave is repeatedly transmitted a plurality of times to form time series data of any point on each scan line.
And A3, utilizing the shear wave detection data on the shear wave detection scanning lines to calculate and obtain the shear wave waveform data on each shear wave detection scanning line.
The calculation of the axial displacement or the axial velocity of the tissue caused by the propagation of the shear wave is performed based on the time-series data of the detected point and the scanning line formed at the same position and different time. The axial displacement or axial velocity calculation can be obtained by calculation by utilizing an autocorrelation or cross-correlation algorithm relative to scan line data which is stored in advance before the ultrasonic acoustic radiation force excitation; or the data of the scanning lines at the same position and different moments can be obtained by performing autocorrelation or cross-correlation algorithm calculation. The axial velocity and the axial displacement can be mutually converted through a relation formula: axial velocity is the differential of axial displacement in the time dimension, and axial displacement is the integral of axial velocity in the time dimension. The axial displacement or axial velocity data calculated for multiple echoes at each scanline location are combined in time sequence to form a shear wave time waveform at that location. As shown in fig. 4, 7 in fig. 4 is the shear wave time waveform at scan line a in fig. 2, 8 in fig. 4 is the shear wave time waveform at scan line B in fig. 3, 9 in fig. 4 is the shear wave time waveform at scan line C in fig. 3, and 10 in fig. 4 is the shear wave time waveform at scan line D in fig. 3.
And A4, calculating the propagation arrival time of the shear wave at each shear wave detection scanning line position according to the shear wave waveform data on each shear wave detection scanning line. There are many methods for determining the propagation arrival time of the shear wave, and the most common method is to determine the time corresponding to the peak point of the axial velocity data or the axial displacement data, which is the arrival time of the shear wave at the position.
A5, detecting the propagation arrival time of the shear wave at the scanning line position according to several adjacent shear waves, and calculating the shear wave speed at the position.
A6, detecting the shear wave waveform data of the scanning line position according to several adjacent shear waves, and calculating the attenuation information of the shear wave, wherein the attenuation information comprises the amplitude attenuation of the shear wave, the frequency attenuation of the shear wave, the waveform change of the shear wave in the time domain, and the like.
The ultrasonic sound radiation force excites the tissue to generate shear waves, the amplitude and frequency of signals and waveforms of the shear waves change along with the increase of the propagation distance, and therefore, the situation that most of the existing machines only display the single parameter of the shear wave speed can be supplemented by evaluating the changes.
Specifically, for shear wave amplitude attenuation, taking FIG. 4 as an example, the peak of the axial displacement data or axial velocity data is determined, and then the-3 dB closest to the peak (of the peak) is found
Double) or-6The cut-off position in dB (1/2 times the peak) is averaged over all positions greater than or equal to the cut-off position. Where 11, 12, 13 and 14 in fig. 3 are the cut-off positions of 7, 8, 9 and 10, respectively. Taking 7 in fig. 4 as an example, the calculation process is as follows:
num ═ length (AmpA ≧ CutoffA) (equation 1)
Where cutoff a is the cut-off position, AmpA represents the amplitude of 7 in fig. 4, num represents the number of points (calculated using the length function) where the amplitude of 7 in fig. 4 is greater than the cut-off position, and meanAmpA represents the average (calculated using the mean average function) of all amplitudes greater than the cut-off position. Likewise, the average values of the amplitudes of 8, 9 and 10 in fig. 4 can be obtained. Therefore, the rate of change of the amplitude is evaluated by differentiating the calculation results of the adjacent shear wave scanning lines. The calculation process is as follows:
meanapdiff (i) ═ meanAmp (i) — meanAmp (i +1) (equation 3)
Wherein i represents the number of shear wave detection scan lines in the transverse direction, meanAmp (i) and meanAmp (i +1) represent the average value of the amplitudes of adjacent positions at a certain depth, and meanAmp diff (i) represents the differential result of the average values of the amplitudes of the adjacent positions at the depth, respectively, thereby evaluating the amplitude change rate.
For shear wave frequency attenuation, we can also evaluate it in a similar way. As shown in fig. 5, where 15, 16, 17 and 18 are the results of the fourier transforms of 7, 8, 9 and 10, respectively, in fig. 4. By fourier transform, the frequency components of the signal can be extracted. The peak value of the Fourier transform result of the axial displacement data or the axial velocity data is firstly determined, and then-3 dB (of the peak value) closest to the peak value is found
Multiple) or-6 dB (1/2 times the peak value), all greater than or equal to the cutoff position are averaged. Where 19, 20, 21 and 22 in figure 5 are the cut-off positions of 15, 16, 17 and 18 respectively. With the figure15 in 5 is taken as an example, and the calculation process is as follows:
fftAmpA ═ abs (fft (ampa)) (equation 4)
fftNum ═ length (fftAmpA ≧ fftCutoffa) (equation 5)
Where AmpA is the amplitude of 7 in fig. 4, fft (AmpA) is the result of fourier transform of AmpA, fftCutoffA is the cut-off position, fftAmpA represents the amplitude of 15 in fig. 5 (calculated using abs absolute value function), fftNum represents the number of points in fig. 5 where the amplitude of 15 is greater than or equal to the cut-off position (calculated using length function), and fftmean AmpA represents the average of all amplitudes greater than or equal to the cut-off position (calculated using mean function). Similarly, the average values of the amplitudes of 16, 17 and 18 can be obtained. Thus, the rate of change of the shear wave frequency is evaluated by differentiating the results of the calculations of adjacent shear wave scan lines.
And the waveform change of the shear wave time domain can be measured by the energy of the signal. It is assumed that an axial displacement or axial velocity signal at a certain axial position and lateral position is taken, since the length of the time series is finite, that is, the energy of the signal is finite. The physical meaning of the energy of a signal refers to the magnitude of the signal transfer capability. If the frequency of the signal is higher, the attenuation is quicker, and the propagation distance is smaller; the lower the frequency of the signal, the greater the energy it has, and the further the propagation distance. Taking 7 in fig. 4 as an example, the calculation process is as follows:
where AmpA is the amplitude of 7 in FIG. 4 and N is the length of the shear wave time domain waveform. | · | represents taking the absolute value of AmpA, and EnergyA represents the energy of the axial displacement or axial velocity signal of time series length N. Similarly, 16, 17 and 18 energies can be obtained. Therefore, the rate of change of the signal energy is evaluated by differentiating the results of the calculations of adjacent shear wave scan lines.
Attenuation information of the shear wave is calculated, wherein the attenuation information comprises information such as amplitude attenuation of the shear wave, frequency attenuation of the shear wave, and waveform change of a shear wave in a time domain.
A7, repeating A2-A6 until obtaining a shear wave velocity information matrix and a shear wave attenuation information matrix within the ROI box capable of covering the whole shear wave elastography.
Specifically, repeating a2-a6, i.e., calculating along the transverse direction and the axial direction, respectively, obtains information matrices such as shear wave velocity, shear wave amplitude attenuation, shear wave frequency attenuation, and shear wave time-domain waveform change in the ROI frame that can cover the entire shear wave elastography. As shown in fig. 3, since the propagation of the shear wave has an attenuation characteristic, the shear wave gradually attenuates after propagating a distance in the lateral direction; furthermore, shear waves do not exist only at a certain depth, but have a certain axial extent. Therefore, when shear wave speed information and shear wave attenuation information are calculated, different detection points are firstly calculated at the same depth along the transverse direction; then, the next depth (axial direction) is calculated for different detection points along the transverse direction, and so on. Therefore, the results of the ROI frame are all calculated, and a shear wave speed information matrix and a shear wave attenuation information matrix which can cover the whole shear wave elastography in the ROI frame are obtained.
And A8, mapping the shear wave velocity information matrix and the shear wave attenuation information matrix obtained in the step A7 into different colors or different gray scales through a special image mapping algorithm, and displaying the different colors or the different gray scales.
The image mapping algorithm may be a direct table look-up mapping algorithm implemented by using a color table or a gray scale table (the color table/gray scale table records the corresponding relationship between color/gray scale and speed information and attenuation information), or an image mapping algorithm implemented by using other calculation forms (for example, the speed information and the attenuation information are subjected to numerical transformation of a specific formula and then are correspondingly matched with color/gray scale codes). The two-dimensional shear wave amplitude attenuation image is displayed and can be used for reflecting the accuracy and reliability of shear wave velocity estimation laterally. The shear wave frequency attenuation image is displayed in a two-dimensional form by means of grayscale or color mapping by performing calculations along the transverse direction and the axial direction, respectively.
Specifically, the result can be mapped to different gray levels or different colors by a special image mapping algorithm and displayed in a two-dimensional form. The shear wave velocity information matrix reflects the elastic information of the tissue to be detected, and the shear wave amplitude attenuation, the shear wave frequency attenuation and the shear wave time domain waveform change laterally reflect the accuracy and reliability of shear wave velocity estimation, and can be used as another method for evaluating the tissue to be detected.
The invention also provides a shear wave elastography system, comprising: the shear wave excitation module is used for exciting at a preset position by using acoustic radiation force to generate shear waves; the detection module is used for transmitting ultrasonic detection wave beams on the shear wave propagation path to detect shear waves so as to acquire shear wave detection data; the calculation module is used for analyzing the shear wave detection data to acquire the speed information and the attenuation information of the shear wave; and the image processing module is used for processing the shear wave velocity information and the attenuation information by utilizing a preset image mapping algorithm and outputting corresponding image information.
In this embodiment, the attenuation information includes shear wave amplitude attenuation information and/or shear wave frequency attenuation information and/or shear wave time domain waveform change information, and the calculation module specifically includes: the shear wave velocity information calculation submodule is used for detecting the propagation arrival time of the shear wave at the scanning line position according to the adjacent shear waves and calculating the shear wave velocity information of the position; the shear wave amplitude attenuation information submodule is used for determining the peak value of axial displacement data or axial speed data of shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave amplitude attenuation information; the shear wave frequency attenuation information submodule is used for determining the peak value of axial displacement data or axial velocity data after Fourier transform, setting an amplitude cut-off region, acquiring the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave frequency attenuation information; and the shear wave time domain waveform change information submodule is used for determining the peak value of the axial displacement data or the axial speed data of the shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the square of the amplitude in the amplitude cut-off region, and comparing the average values of the squares of the amplitudes of adjacent shear wave detection scanning lines to acquire the shear wave time domain waveform change information.
In this embodiment, a working process of a shear wave elastography system corresponds to the shear wave elastography method, which is not described herein again.
The shear wave elastic imaging method and the shear wave elastic imaging system can be used for measuring the speed information and the attenuation information in the shear wave propagation process and displaying the speed information and the attenuation information, so that the technical problem that the tissue information is not completely reflected due to the fact that a shear wave elastic imaging system is single in measuring object in the prior art is solved, the shear wave elastic imaging method and the shear wave elastic imaging system can be used for reflecting the speed and the attenuation information of the shear wave, and the reliability and the comprehensiveness of shear wave detection are improved.
In addition, the shear wave velocity information and the shear wave attenuation information are calculated by adopting the axial velocity or axial displacement data of the shear wave obtained by carrying out algorithms such as autocorrelation, cross correlation and the like on the shear wave detection data after beam synthesis, so that the method can be conveniently implemented in a system.
The invention can be widely applied to various shear wave elastic imaging systems.
While the preferred embodiments of the present invention have been illustrated and described, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (9)
1. A method of shear wave elastography, comprising the steps of:
s1, exciting to generate shear waves at a preset position by using acoustic radiation force;
s2, transmitting ultrasonic detection wave beams for multiple times on the shear wave propagation path to detect shear waves so as to obtain shear wave detection data;
s3, analyzing the shear wave detection data to obtain the speed information and attenuation information of the shear wave;
s4, processing the speed information and attenuation information of the shear wave by using a preset image mapping algorithm, and outputting corresponding image information;
the step S3 specifically includes:
the method comprises the steps of generating shear wave waveform data by utilizing shear wave detection data, and calculating speed information and attenuation information of shear waves according to the shear wave waveform data of adjacent shear wave detection scanning lines, wherein the attenuation information comprises shear wave amplitude attenuation information and/or shear wave frequency attenuation information and/or shear wave time domain waveform change information.
2. A method of shear wave elastography according to claim 1, wherein said step S3 includes the sub-step of obtaining shear wave velocity information:
and calculating the shear wave velocity information of the position according to the propagation arrival time of the shear wave of the adjacent shear wave detection scanning line position.
3. A method of shear wave elastography according to claim 1, wherein said step S3 includes the sub-step of obtaining shear wave amplitude attenuation information:
determining the peak value of the axial displacement data or the axial velocity data of the shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave amplitude attenuation information.
4. A method of shear wave elastography according to claim 1, wherein said step S3 includes the sub-step of obtaining shear wave frequency attenuation information:
the method comprises the steps of processing shear wave waveform data by utilizing Fourier transform, determining the peak value of axial displacement data or axial velocity data after Fourier transform, setting an amplitude cut-off region, obtaining the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to obtain shear wave frequency attenuation information.
5. A method for shear wave elastography according to claim 1, wherein said step S3 includes a sub-step of obtaining shear wave time domain waveform change information:
determining the peak value of axial displacement data or axial velocity data of shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the squares of the amplitudes in the amplitude cut-off region, and comparing the average values of the squares of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave time domain waveform change information.
6. A method according to any of claims 1 to 5, wherein said step S1 is specifically: in the region of interest of shear wave elastography, shear waves are generated by utilizing acoustic radiation force excitation at a preset position, and the times of ultrasonic acoustic radiation force excitation are determined according to the size of the region of interest and the propagation distance of the shear waves generated by single excitation.
7. The method according to claim 6, wherein the step S2 is specifically: the method includes the steps of transmitting ultrasonic detection beams multiple times on a shear wave propagation path to obtain echo signals, processing the echo signals with a beam synthesis algorithm to generate a plurality of shear wave detection scanlines, and generating shear wave detection data at multiple times on each shear wave detection scanline.
8. A shear wave elastography system, comprising:
the shear wave excitation module is used for exciting at a preset position by using acoustic radiation force to generate shear waves;
the detection module is used for transmitting ultrasonic detection wave beams on the shear wave propagation path to detect shear waves so as to acquire shear wave detection data;
the calculation module is used for analyzing the shear wave detection data to acquire the speed information and the attenuation information of the shear wave, wherein the attenuation information comprises shear wave amplitude attenuation information and/or shear wave frequency attenuation information and/or shear wave time domain waveform change information;
and the image processing module is used for processing the shear wave velocity information and the attenuation information by utilizing a preset image mapping algorithm and outputting corresponding image information.
9. A shear wave elastography system according to claim 8, wherein said computing means particularly comprises: the shear wave velocity information calculation submodule is used for detecting the propagation arrival time of the shear wave at the scanning line position according to the adjacent shear waves and calculating the shear wave velocity information of the position;
the shear wave amplitude attenuation information submodule is used for determining the peak value of axial displacement data or axial speed data of shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave amplitude attenuation information;
the shear wave frequency attenuation information submodule is used for determining the peak value of axial displacement data or axial velocity data after Fourier transform, setting an amplitude cut-off region, acquiring the average value of the amplitudes in the amplitude cut-off region, and comparing the average values of the amplitudes of adjacent shear wave detection scanning lines to acquire shear wave frequency attenuation information;
and the shear wave time domain waveform change information submodule is used for determining the peak value of the axial displacement data or the axial speed data of the shear wave waveform data, setting an amplitude cut-off region, acquiring the average value of the square of the amplitude in the amplitude cut-off region, and comparing the average values of the squares of the amplitudes of adjacent shear wave detection scanning lines to acquire the shear wave time domain waveform change information.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710859560.6A CN107510474B (en) | 2017-09-21 | 2017-09-21 | Shear wave elastic imaging method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710859560.6A CN107510474B (en) | 2017-09-21 | 2017-09-21 | Shear wave elastic imaging method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107510474A CN107510474A (en) | 2017-12-26 |
| CN107510474B true CN107510474B (en) | 2020-07-10 |
Family
ID=60726300
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710859560.6A Active CN107510474B (en) | 2017-09-21 | 2017-09-21 | Shear wave elastic imaging method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107510474B (en) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107970043B (en) * | 2017-12-28 | 2021-01-19 | 深圳开立生物医疗科技股份有限公司 | Shear wave detection method and device |
| CN110494082B (en) * | 2018-04-11 | 2023-05-05 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasound elastography method and system |
| WO2020042020A1 (en) * | 2018-08-29 | 2020-03-05 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic elasticity detection device and shear wave elasticity imaging method and apparatus |
| CN111655154B (en) * | 2018-12-04 | 2021-07-20 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method and ultrasonic imaging system |
| CN110604598B (en) * | 2018-12-04 | 2021-01-15 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method and ultrasonic imaging system |
| CN111388012B (en) * | 2019-01-02 | 2021-10-08 | 无锡海斯凯尔医学技术有限公司 | Method, device and system for detecting tissue hardness |
| CN110243950B (en) * | 2019-05-15 | 2021-10-08 | 天津大学 | Acoustic radiation force pulse excitation method based on amplitude modulation |
| CN110477947B (en) * | 2019-08-14 | 2022-04-15 | 中国科学院苏州生物医学工程技术研究所 | Plane wave beam synthesis method, system, storage medium and equipment based on deep learning |
| CN114287968A (en) * | 2019-09-27 | 2022-04-08 | 深圳迈瑞生物医疗电子股份有限公司 | Elasticity imaging method, system and computer readable storage medium |
| CN110742644B (en) * | 2019-09-29 | 2022-07-08 | 深圳大学 | Elastography system, elastography method and storage medium |
| CN111449681B (en) * | 2020-04-08 | 2023-09-08 | 深圳开立生物医疗科技股份有限公司 | Shear wave imaging method, device, equipment and readable storage medium |
| CN112244889A (en) * | 2020-10-15 | 2021-01-22 | 中国科学院苏州生物医学工程技术研究所 | Method for determining vibration element array, puncture needle imaging method and ultrasonic equipment |
| CN112799078B (en) * | 2021-04-15 | 2021-07-02 | 深圳中科乐普医疗技术有限公司 | Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment |
| CN114403920A (en) * | 2022-01-24 | 2022-04-29 | 陕西师范大学 | Time-division multipoint excitation shear wave fusion imaging method and system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5606971A (en) * | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
| CN103054552A (en) * | 2012-12-24 | 2013-04-24 | 深圳先进技术研究院 | Method and system for measuring biological tissue viscoelasticity |
| CN103462643A (en) * | 2013-09-29 | 2013-12-25 | 深圳市开立科技有限公司 | Shear wave speed measurement method, device and system |
| CN104825195A (en) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | Shear wave viscoelasticity imaging method and system |
| CN106175831A (en) * | 2016-09-28 | 2016-12-07 | 重庆理工大学 | A kind of tissue viscoelasticity based on shearing wave amplitude and phase-detection measuring method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10136876B2 (en) * | 2012-09-27 | 2018-11-27 | Mayo Foundation For Medical Education And Research | System and method for local estimation of nonlinear tissue elasticity with acoustic radiation force |
-
2017
- 2017-09-21 CN CN201710859560.6A patent/CN107510474B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5606971A (en) * | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
| CN103054552A (en) * | 2012-12-24 | 2013-04-24 | 深圳先进技术研究院 | Method and system for measuring biological tissue viscoelasticity |
| CN103462643A (en) * | 2013-09-29 | 2013-12-25 | 深圳市开立科技有限公司 | Shear wave speed measurement method, device and system |
| CN104825195A (en) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | Shear wave viscoelasticity imaging method and system |
| CN106175831A (en) * | 2016-09-28 | 2016-12-07 | 重庆理工大学 | A kind of tissue viscoelasticity based on shearing wave amplitude and phase-detection measuring method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107510474A (en) | 2017-12-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107510474B (en) | Shear wave elastic imaging method and system | |
| US10159466B2 (en) | Sparse tracking in acoustic radiation force impulse imaging | |
| JP6053860B2 (en) | Computer program for determining shear wave information in medical ultrasound imaging systems | |
| US9468421B2 (en) | Visualization of associated information in ultrasound shear wave imaging | |
| CN104510499B (en) | Shear wave detection in medical ultrasound imaging | |
| CN106419961B (en) | Adaptive motion estimation in acoustic radiation force imaging | |
| US11154277B2 (en) | Tissue viscoelastic estimation from shear velocity in ultrasound medical imaging | |
| CN114391868A (en) | Quantitative viscoelastic ultrasound imaging | |
| JP2012081269A5 (en) | ||
| CN107049361B (en) | Sound velocity imaging using shear waves | |
| KR20130057435A (en) | Imaging method and apparatus using shear waves | |
| CN106955125B (en) | Motion independence in acoustic radiation force pulse imaging | |
| CN110893103A (en) | Angle for ultrasound-based shear wave imaging | |
| US11963824B2 (en) | Shear wave imaging based on ultrasound with increased pulse repetition interval | |
| CN108338808B (en) | Shear velocity imaging using coherence | |
| JP6492230B2 (en) | SPECTRUM ANALYZER, SPECTRUM ANALYSIS METHOD, AND ULTRASONIC IMAGING DEVICE | |
| WO2023278445A1 (en) | System and method for non-invasive determination of pressure in a biological compartment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |