CN112799078B - Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment - Google Patents
Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment Download PDFInfo
- Publication number
- CN112799078B CN112799078B CN202110404171.0A CN202110404171A CN112799078B CN 112799078 B CN112799078 B CN 112799078B CN 202110404171 A CN202110404171 A CN 202110404171A CN 112799078 B CN112799078 B CN 112799078B
- Authority
- CN
- China
- Prior art keywords
- deformation
- curves
- time
- principal component
- shear wave
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Clinical applications
- A61B8/0833—Clinical applications involving detecting or locating foreign bodies or organic structures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Signal Processing (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
According to the shear wave propagation velocity detection method and system and the ultrasonic imaging device, provided by the invention, spectrum analysis is carried out by acquiring the deformation-time curves of the first position and the second position in the elastic observation area of the biological tissue, and whether the range difference of the principal component frequency ranges of the two spectrums exceeds a preset threshold value is judged. If the frequency range exceeds the preset frequency range, band-pass filtering is carried out on the deformation-time curves corresponding to the two frequency spectrums according to the main component frequency range of the frequency spectrum of the detection position farther away from the shear wave source, so that the approach of the main component frequency ranges of the shear waves of the two detection positions can be ensured. And finally, estimating the shear wave velocity according to the two deformation-time curves, thereby improving the accuracy of detecting the shear wave propagation velocity.
Description
Technical Field
The invention relates to the field of medical instruments, in particular to a method and a system for detecting the propagation velocity of shear waves and ultrasonic imaging equipment.
Background
The shear wave elastography technology can realize real-time hardness quantitative detection of biological tissues and provide a basis for clinically judging the pathological changes of the tissues. Its basic principle is as follows: the acoustic radiation force focusing impact energy can generate shear waves in tissues, and due to the fact that the propagation speeds of the shear waves in the tissues with different hardness are different, the hardness and softness of the positions can be indirectly reflected by detecting the propagation speeds of the shear waves in the different positions. And performing pseudo-color mapping according to the wave velocity of the shear wave to obtain shear wave elastic imaging. Therefore, the accuracy of shear wave velocity estimation is crucial to the determination of tissue elasticity.
The shear wave can cause the displacement deformation of the tissue in the tissue propagation process, the displacement deformation of a target area is observed and detected by using an ultrafast frame frequency imaging technology, and a shear wave vibration curve of an observed position is reconstructed through the displacement deformation. The conventional shear wave velocity estimation method is a Time-of-flight (tof) method, in which a peak matching is performed on shear wave vibration curves reconstructed at different positions in an observation area to find a peak difference Time interval, which can be understood as the Time taken for a shear wave to propagate from the two detection positions. And the distance between the two sensing locations is known, and the average velocity of the shear wave passing between the two sensing locations is found by dividing the distance by the time.
The shear waves impacted by the acoustic radiation force decay very rapidly as the tissue propagates, and the farther away from the source of the impact, the lower the signal-to-noise ratio of the displacement estimate. These factors directly affect the accuracy of shear wave velocity estimation.
Disclosure of Invention
The invention provides a method and a system for detecting the propagation velocity of a shear wave and an ultrasonic imaging device, which are used for improving the accuracy of the detection of the propagation velocity of the shear wave.
An embodiment provides a method for detecting a propagation velocity of a shear wave, including:
acquiring deformation-time curves of a first position and a second position in an elastic observation area of a biological tissue, wherein the deformation of the first position and the second position is caused by shear waves, and the shear waves pass through the first position and then pass through the second position;
respectively carrying out Fourier transform on the deformation-time curves of the first position and the second position to obtain corresponding frequency spectrum curves;
determining principal component frequency ranges of the frequency spectrum curves of the first position and the second position, and comparing to obtain a range difference of the two principal component frequency ranges;
judging whether the range difference exceeds a preset threshold value, if so, filtering deformation-time curves of the first position and the second position respectively by using the principal component frequency range of the second position to obtain the deformation-time curves of the first position and the second position after filtering;
obtaining a time interval corresponding to the peak value of the two curves according to the two filtered deformation-time curves;
and obtaining the speed value of the shear wave propagating between the first position and the second position according to the time interval and the distance between the first position and the second position.
An embodiment provides the method, wherein the determining the principal component frequency ranges of the spectral curves at the first and second positions includes:
determining the spectrum width corresponding to the preset amplitude in the spectrum curve of the first position as a principal component frequency range;
and determining the spectrum width corresponding to the preset amplitude value in the spectrum curve of the second position as a principal component frequency range.
In one embodiment, the predetermined amplitude is one of 30% to 70% of a peak value in a spectrum curve.
An embodiment provides the method, wherein the obtaining a "deformation-time" curve of a first location and a second location in an elastic observation region of a biological tissue comprises:
acquiring an elastic observation area of biological tissue;
exciting the biological tissue to excite shear waves to pass through the elastic observation area, and acquiring ultrasonic echo signals of a first position and a second position preset in the elastic observation area by adopting an ultrahigh frame frequency imaging technology;
and carrying out displacement estimation operation according to the ultrasonic echo signals to obtain deformation-time curves of the first position and the second position.
An embodiment provides the method, wherein the acquiring an elastic observation region of a biological tissue includes:
transmitting ultrasonic waves to biological tissues and receiving ultrasonic echo signals;
generating an ultrasonic image according to the ultrasonic echo signal and displaying the ultrasonic image;
an elastic viewing zone of the biological tissue is obtained in response to user input instructions for determining the elastic viewing zone on the ultrasound image.
An embodiment provides the method, wherein the first position and the second position have the same depth.
An embodiment provides a system for detecting a propagation velocity of a shear wave, comprising:
an acquisition module, configured to acquire "deformation-time" curves of a first position and a second position in an elastic observation region of a biological tissue, where the deformation of the first position and the deformation of the second position are caused by a shear wave, and the shear wave passes through the first position and then passes through the second position;
the transformation module is used for respectively carrying out Fourier transformation on the deformation-time curves of the first position and the second position to obtain corresponding frequency spectrum curves;
the calculation and judgment module is used for determining the principal component frequency ranges of the frequency spectrum curves of the first position and the second position and comparing the principal component frequency ranges to obtain the range difference of the two principal component frequency ranges; judging whether the range difference exceeds a preset threshold value, if so, filtering deformation-time curves of the first position and the second position respectively by using the principal component frequency range of the second position to obtain the deformation-time curves of the first position and the second position after filtering; obtaining a time interval corresponding to the peak value of the two curves according to the two filtered deformation-time curves; and obtaining the speed value of the shear wave propagating between the first position and the second position according to the time interval and the distance between the first position and the second position.
An embodiment provides the measurement system, further comprising a display module, where the display module is configured to display the speed value.
An embodiment provides an ultrasound imaging apparatus comprising:
the ultrasonic probe is used for transmitting ultrasonic waves and receiving echoes of the ultrasonic waves;
a display;
a processor for obtaining 'deformation-time' curves of a first position and a second position in an elastic observation region of a biological tissue, the deformations of the first position and the second position being caused by shear waves, the shear waves passing through the first position first and then the second position; respectively carrying out Fourier transform on the deformation-time curves of the first position and the second position to obtain corresponding frequency spectrum curves; determining principal component frequency ranges of the frequency spectrum curves of the first position and the second position, and comparing to obtain a range difference of the two principal component frequency ranges; judging whether the range difference exceeds a preset threshold value, if so, filtering deformation-time curves of the first position and the second position respectively by using the principal component frequency range of the second position to obtain the deformation-time curves of the first position and the second position after filtering; obtaining a time interval corresponding to the peak value of the two curves according to the two filtered deformation-time curves; and obtaining the speed value of the shear wave propagating between the first position and the second position according to the time interval and the distance between the first position and the second position.
An embodiment provides a computer readable storage medium having a program stored thereon, the program being executable by a processor to implement a method as described above.
According to the shear wave propagation velocity detection method, the shear wave propagation velocity detection system and the ultrasonic imaging device in the embodiments, spectral analysis is performed by acquiring the deformation-time curves of the first position and the second position in the elastic observation region of the biological tissue, and whether the range difference of the principal component frequency ranges of the two spectra exceeds a preset threshold value is judged. If the frequency range exceeds the preset frequency range, band-pass filtering is carried out on the deformation-time curves corresponding to the two frequency spectrums according to the main component frequency range of the frequency spectrum of the detection position farther away from the shear wave source, so that the approach of the main component frequency ranges of the shear waves of the two detection positions can be ensured. And finally, estimating the shear wave velocity according to the two deformation-time curves, thereby improving the accuracy of detecting the shear wave propagation velocity.
Drawings
FIG. 1 is a block diagram of an embodiment of a system for detecting the propagation velocity of a shear wave according to the present invention;
FIG. 2 is a flow chart of an embodiment of a method for detecting a propagation velocity of a shear wave according to the present invention;
FIG. 3 is a schematic diagram of an ultrasound image in the system for detecting the propagation velocity of a shear wave according to the present invention;
FIG. 4 is a "deformation-time" graph of a first position and a second position in a shear wave propagation velocity detection system provided by the present invention;
FIG. 5 is a graph of frequency spectra of a first location and a second location in a shear wave propagation velocity detection system provided by the present invention;
FIG. 6 is a schematic diagram of two principal component frequency ranges of the spectral curves of FIG. 5;
FIG. 7 is a graph of "deformation versus time" after filtering for a first location and a second location;
fig. 8 is a block diagram of an ultrasound imaging apparatus according to an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the following detailed description and accompanying drawings. Wherein like elements in different embodiments are numbered with like associated elements. In the following description, numerous details are set forth in order to provide a better understanding of the present application. However, those skilled in the art will readily recognize that some of the features may be omitted or replaced with other elements, materials, methods in different instances. In some instances, certain operations related to the present application have not been shown or described in detail in order to avoid obscuring the core of the present application from excessive description, and it is not necessary for those skilled in the art to describe these operations in detail, so that they may be fully understood from the description in the specification and the general knowledge in the art.
Furthermore, the features, operations, or characteristics described in the specification may be combined in any suitable manner to form various embodiments. Also, the various steps or actions in the method descriptions may be transposed or transposed in order, as will be apparent to one of ordinary skill in the art. Thus, the various sequences in the specification and drawings are for the purpose of describing certain embodiments only and are not intended to imply a required sequence unless otherwise indicated where such sequence must be followed.
The numbering of the components as such, e.g., "first", "second", etc., is used herein only to distinguish the objects as described, and does not have any sequential or technical meaning. The term "connected" and "coupled" when used in this application, unless otherwise indicated, includes both direct and indirect connections (couplings).
The invention judges the deviation range of the central frequency positions of two frequency spectrums by carrying out frequency spectrum analysis on the deformation-time curves of two different detection positions; if the frequency range exceeds the range threshold value, the shear wave principal component frequency range is selected according to the frequency spectrum of the detection position farther away from the shear wave source, and then the displacement estimation result is subjected to band-pass filtering by the frequency range, so that the shear wave frequency principal component frequency ranges of the two detection positions are close to each other. The method improves the influence of the principal component frequency deviation on the wave velocity estimation when the shear wave is transmitted in the tissue, so that the shear wave velocity estimation result is more accurate, and the effectiveness of the result is improved. The following is a detailed description of specific examples.
As shown in fig. 1, the system for detecting the propagation velocity of a shear wave according to the present invention includes an obtaining module 10, a transforming module 20, a calculating and judging module 30, and a displaying module 40.
The acquisition module 10 is configured to acquire a "deformation-time" curve of a first location and a second location in an elastic viewing area of a biological tissue. The deformations of the first and second locations are caused by the shear wave and the shear wave passes the first location before the second location, in other words the first location is before the second location as seen in the propagation path of the shear wave.
The transform module 20 is configured to perform fourier transform on the "deformation-time" curves of the first location and the second location, respectively, to obtain corresponding spectrum curves.
The calculation and judgment module 30 is configured to determine principal component frequency ranges of the spectrum curves at the first position and the second position, and compare the principal component frequency ranges to obtain a range difference between the two principal component frequency ranges; judging whether the range difference exceeds a preset threshold value, if so, filtering the deformation-time curves of the first position and the second position respectively by using the principal component frequency range of the second position to obtain the deformation-time curves of the first position and the second position after filtering; obtaining a time interval corresponding to the peak value of the two curves according to the two filtered deformation-time curves; and obtaining the speed value of the shear wave propagating between the first position and the second position according to the time interval and the distance between the first position and the second position.
The display module 40 is used for displaying the velocity value of the propagation velocity of the shear wave.
Therefore, when the frequency range deviation of the frequency spectrum curves at the two positions is larger, the deformation-time curves at the two positions are respectively filtered by adopting the frequency range of the second position farther away from the shear wave source, so that the deformation-time curves at the two positions are drawn closer, the subsequent shear wave propagation speed is calculated according to the deformation-time curves, and the accuracy of the propagation speed is improved.
As shown in fig. 2, the detection system detects the propagation velocity of the shear wave, and the process may include the following steps:
In this embodiment, the detection system may further include an ultrasonic probe and an input device, where the input device is used to receive input from a user, and may be, for example, one or more of a mouse, a keyboard, a control panel, a trackball, a touch pad, a touch screen, and the like. As shown in fig. 3, the acquiring module 10 acquires an elastic observation region a of a biological tissue B, for example, may control an ultrasonic probe to transmit ultrasonic waves to the biological tissue B and receive an ultrasonic echo signal; an ultrasound image P is generated according to the ultrasound echo signal and displayed on a display interface of the display module 40. The user determines the elastic observation region a on the ultrasound image P through an input device, for example, by selecting an area through a mouse or adopting a default area of the system as the elastic observation region a. The acquisition module 10 obtains the elastic observation region a of the biological tissue B in response to an instruction input by the user to determine the elastic observation region a on the ultrasound image P.
The acquisition module 10 controls the ultrasound probe to excite the biological tissue to excite shear waves to propagate through the elastic observation region a, for example by impinging acoustic radiation forces to excite shear waves to propagate through the region a. The acquisition module 10 controls the ultrasound probe to acquire ultrasound echo signals of a first position and a second position preset in the elastic observation region a by using an ultrahigh frame frequency imaging technology, for example, transmitting ultrasound waves at a preset frame frequency (frequency) and receiving the echo signals of the ultrasound waves, where the preset frame frequency is greater than or equal to 1000 frames/second. The obtaining module 10 performs a displacement estimation operation according to the ultrasound echo signal to obtain a "deformation-time" curve of the first position and the second position, for example, the ultrasound echo signal is subjected to an image processing to obtain an ultrasound image of a period of time, and the change of the first position in each ultrasound image frame in the ultrasound image is compared to obtain the deformation of the first position, so as to obtain the "deformation-time" curve of the first position, as shown in m1 of fig. 4, and in the same manner, the "deformation-time" curve of the second position is obtained, as shown in n1 of fig. 4.
The first position and the second position are two preset positions of the elastic observation area a, which is determined, and the first position and the second position are also determined, in this embodiment, the first position and the second position are the center m on the left side and the center n on the right side of the elastic observation area a, respectively. The first position m and the second position n are determined and the distance L between them is known. In this embodiment, the first position m and the second position n have the same depth, that is, the first position m and the second position n are two different positions on the same horizontal line, in other words, the first position m and the second position n are two different lateral detection positions. The example will be described with the shear wave source closer to the first position m.
And 3, determining the principal component frequency range of the spectrum curve at the first position and the principal component frequency range of the spectrum curve at the second position by the calculation and judgment module 30. The principal component frequency range is, as the name implies, the principal frequency range of the spectral curve. For example, determining a spectral width corresponding to a preset amplitude value in a spectral curve of the first position as a principal component frequency range; and determining the spectral width corresponding to the preset amplitude value in the spectral curve of the second position as the principal component frequency range. The preset amplitude may be set as desired, for example, one of 30% -70% of the peak in the spectral curve. In other words, in this embodiment, a spectrum width corresponding to 50% of a peak in a spectrum curve at a first position is used as a principal component frequency range (a width indicated by a dotted arrow in fig. 6); the spectral width corresponding to 50% of the peak in the spectral curve at the second position is taken as the principal component frequency range (the width indicated by the solid arrow in fig. 6).
The calculation and judgment module 30 compares the range difference between the two principal component frequency ranges, for example, the absolute value of the difference between the two principal component frequency ranges is taken as the range difference.
And 4, judging whether the range difference exceeds a preset threshold value by the calculation judgment module 30, if so, entering the step 5, and otherwise, entering the step 6.
And step 6, the calculation and judgment module 30 obtains the time interval corresponding to the peak value of the two curves according to the two deformation-time curves. That is, if the range difference does not exceed the preset threshold, which indicates that the difference between the "deformation-time" curves of the two positions is acceptable, the calculation and determination module 30 obtains the time interval t corresponding to the peak value of the two "deformation-time" curves according to the "deformation-time" curve of the first position and the "deformation-time" curve of the second position, as shown in fig. 4. If the two deformation-time curves are filtered, the calculation and judgment module 30 obtains the time interval t corresponding to the peak value of the two filtered deformation-time curves according to the deformation-time curve m2 after the first position filtering and the deformation-time curve n2 after the second position filtering, as shown in fig. 7. Specifically, the calculation and judgment module 30 searches for the maximum value in the "deformation-time" curve, which is the peak value, so as to obtain the time corresponding to the peak value, and the time difference between the times corresponding to the two peak values is the time interval t.
Step 7, the calculation and judgment module 30 obtains a velocity value v of the shear wave propagating between the first position m and the second position n according to the time interval t and the distance L between the first position and the second position, for example, v = L/t. The velocity value reflects the elasticity of the elastic viewing area. The calculation and judgment module 30 may also perform pseudo-color mapping (rendering) on the elastic observation region based on the velocity value v on the ultrasound image of the biological tissue, so as to facilitate observation by the user, and in the flow shown in fig. 2, a plurality of first positions and a plurality of second positions may be provided, and the first positions and the second positions appear in pairs, and one pair corresponds to one elastic observation region, so that elasticity (shear wave propagation velocity) of a plurality of elastic observation regions can be obtained, which is beneficial for the user to grasp soft and hard conditions of different positions of the same tissue.
And 8, displaying the velocity value v of the propagation velocity of the shear wave and a rendered ultrasonic image of the biological tissue by the display module 40. The display module 40 may be various types of displays.
In conclusion, the method and the system provided by the invention can effectively improve the accuracy of measuring the propagation speed of the shear wave and provide reliable basis for a doctor to evaluate the hardness and softness of tissues.
As shown in fig. 8, the present invention also provides an ultrasonic imaging apparatus including: an ultrasound probe 81, a transmit and receive module 82, a processor 83, a display 84 and an input device 85.
An ultrasound probe 81, which includes at least one transducer, transmits ultrasound waves and receives echoes of the ultrasound waves.
And a transmitting and receiving module 82, which is respectively connected to the ultrasonic probe 81 and the processor 83, and is configured to transmit the transmitting sequence and the receiving sequence to the ultrasonic probe 81, and transmit the echo of the ultrasonic wave received by the ultrasonic probe 81 to the processor 83.
The input device 85 is used for receiving input of a user, and may be one or more of a mouse, a keyboard, a control panel, a track ball, a touch pad, a touch screen, and the like.
The processor 83 is configured to generate a transmission sequence for controlling the ultrasound probe 81 to transmit the ultrasound waves, generate a reception sequence for controlling the ultrasound probe 81 to receive the echoes of the ultrasound waves, and process the echoes received by the ultrasound probe to generate the ultrasound images.
In this embodiment, the processor 83 is further configured to obtain a "deformation-time" curve for the first location and the second location in the elastic viewing area of the biological tissue; respectively carrying out Fourier transform on the deformation-time curves of the first position and the second position to obtain corresponding frequency spectrum curves; determining principal component frequency ranges of the frequency spectrum curves of the first position and the second position, and comparing to obtain a range difference of the two principal component frequency ranges; judging whether the range difference exceeds a preset threshold value, if so, filtering the deformation-time curves of the first position and the second position respectively by using the principal component frequency range of the second position to obtain the deformation-time curves of the first position and the second position after filtering; obtaining a time interval corresponding to the peak value of the two curves according to the two filtered deformation-time curves; obtaining the speed value of the shear wave propagating between the first position and the second position according to the time interval and the distance between the first position and the second position; the speed values are displayed by the display 84.
It can be seen that the functions of the above system can be performed by the processor 83, in other words, the processor 83 includes the following components: the device comprises an acquisition module 10, a transformation module 20 and a calculation judgment module 30. Since the specific functions of these functional modules included in the processor 83 and the detailed process of detecting the propagation velocity of the shear wave are described in the above embodiments, they are not described in detail herein.
Those skilled in the art will appreciate that all or part of the functions of the various methods in the above embodiments may be implemented by hardware, or may be implemented by computer programs. When all or part of the functions of the above embodiments are implemented by a computer program, the program may be stored in a computer-readable storage medium, and the storage medium may include: a read only memory, a random access memory, a magnetic disk, an optical disk, a hard disk, etc., and the program is executed by a computer to realize the above functions. For example, the program may be stored in a memory of the device, and when the program in the memory is executed by the processor, all or part of the functions described above may be implemented. In addition, when all or part of the functions in the above embodiments are implemented by a computer program, the program may be stored in a storage medium such as a server, another computer, a magnetic disk, an optical disk, a flash disk, or a removable hard disk, and may be downloaded or copied to a memory of a local device, or may be version-updated in a system of the local device, and when the program in the memory is executed by a processor, all or part of the functions in the above embodiments may be implemented.
The present invention has been described in terms of specific examples, which are provided to aid understanding of the invention and are not intended to be limiting. For a person skilled in the art to which the invention pertains, several simple deductions, modifications or substitutions may be made according to the idea of the invention.
Claims (10)
1. A method of detecting the propagation velocity of a shear wave, comprising:
acquiring deformation-time curves of a first position and a second position in an elastic observation area of a biological tissue, wherein the deformation of the first position and the second position is caused by shear waves, and the shear waves pass through the first position and then pass through the second position;
respectively carrying out Fourier transform on the deformation-time curves of the first position and the second position to obtain corresponding frequency spectrum curves;
determining principal component frequency ranges of the frequency spectrum curves of the first position and the second position, and comparing to obtain a range difference of the two principal component frequency ranges;
judging whether the range difference exceeds a preset threshold value, if so, filtering deformation-time curves of the first position and the second position respectively by using the principal component frequency range of the second position to obtain the deformation-time curves of the first position and the second position after filtering;
obtaining a time interval corresponding to the peak value of the two curves according to the two filtered deformation-time curves;
and obtaining the speed value of the shear wave propagating between the first position and the second position according to the time interval and the distance between the first position and the second position.
2. The method for detecting according to claim 1, wherein said determining the principal component frequency range of the spectral curves of the first and second locations comprises:
determining the spectrum width corresponding to the preset amplitude in the spectrum curve of the first position as a principal component frequency range;
and determining the spectrum width corresponding to the preset amplitude value in the spectrum curve of the second position as a principal component frequency range.
3. The detection method of claim 2, wherein the predetermined amplitude is one of 30% to 70% of a peak in the spectral curve.
4. The method of claim 1, wherein obtaining a "deformation-time" curve for a first location and a second location in an elastic viewing area of a biological tissue comprises:
acquiring an elastic observation area of biological tissue;
exciting the biological tissue to excite shear waves to pass through the elastic observation area, and acquiring ultrasonic echo signals of a first position and a second position preset in the elastic observation area by adopting an ultrahigh frame frequency imaging technology;
and carrying out displacement estimation operation according to the ultrasonic echo signals to obtain deformation-time curves of the first position and the second position.
5. The method for detecting according to claim 4, wherein said obtaining an elastic observation region of a biological tissue comprises:
transmitting ultrasonic waves to biological tissues and receiving ultrasonic echo signals;
generating an ultrasonic image according to the ultrasonic echo signal and displaying the ultrasonic image;
an elastic viewing zone of the biological tissue is obtained in response to user input instructions for determining the elastic viewing zone on the ultrasound image.
6. The detection method of claim 5, wherein the first and second locations are the same depth.
7. A system for detecting the propagation velocity of a shear wave, comprising:
an acquisition module, configured to acquire "deformation-time" curves of a first position and a second position in an elastic observation region of a biological tissue, where the deformation of the first position and the deformation of the second position are caused by a shear wave, and the shear wave passes through the first position and then passes through the second position;
the transformation module is used for respectively carrying out Fourier transformation on the deformation-time curves of the first position and the second position to obtain corresponding frequency spectrum curves;
the calculation and judgment module is used for determining the principal component frequency ranges of the frequency spectrum curves of the first position and the second position and comparing the principal component frequency ranges to obtain the range difference of the two principal component frequency ranges; judging whether the range difference exceeds a preset threshold value, if so, filtering deformation-time curves of the first position and the second position respectively by using the principal component frequency range of the second position to obtain the deformation-time curves of the first position and the second position after filtering; obtaining a time interval corresponding to the peak value of the two curves according to the two filtered deformation-time curves; and obtaining the speed value of the shear wave propagating between the first position and the second position according to the time interval and the distance between the first position and the second position.
8. The detection system of claim 7, further comprising a display module for displaying the velocity values.
9. An ultrasound imaging apparatus, comprising:
the ultrasonic probe is used for transmitting ultrasonic waves and receiving echoes of the ultrasonic waves;
a display;
a processor for obtaining 'deformation-time' curves of a first position and a second position in an elastic observation region of a biological tissue, the deformations of the first position and the second position being caused by shear waves, the shear waves passing through the first position first and then the second position; respectively carrying out Fourier transform on the deformation-time curves of the first position and the second position to obtain corresponding frequency spectrum curves; determining principal component frequency ranges of the frequency spectrum curves of the first position and the second position, and comparing to obtain a range difference of the two principal component frequency ranges; judging whether the range difference exceeds a preset threshold value, if so, filtering deformation-time curves of the first position and the second position respectively by using the principal component frequency range of the second position to obtain the deformation-time curves of the first position and the second position after filtering; obtaining a time interval corresponding to the peak value of the two curves according to the two filtered deformation-time curves; and obtaining the speed value of the shear wave propagating between the first position and the second position according to the time interval and the distance between the first position and the second position.
10. A computer-readable storage medium, characterized in that the medium has stored thereon a program which is executable by a processor to implement the method according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110404171.0A CN112799078B (en) | 2021-04-15 | 2021-04-15 | Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110404171.0A CN112799078B (en) | 2021-04-15 | 2021-04-15 | Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112799078A CN112799078A (en) | 2021-05-14 |
| CN112799078B true CN112799078B (en) | 2021-07-02 |
Family
ID=75811420
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110404171.0A Active CN112799078B (en) | 2021-04-15 | 2021-04-15 | Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112799078B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114136249B (en) * | 2021-11-30 | 2023-08-22 | 国网上海市电力公司 | Transformer winding deformation ultrasonic detection signal denoising method |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20200163649A1 (en) * | 2017-06-05 | 2020-05-28 | Mayo Foundation For Medical Education And Research | Estimating Phase Velocity Dispersion in Ultrasound Elastography Using a Multiple Signal Classification |
| CN107616814A (en) * | 2017-08-25 | 2018-01-23 | 深圳中科乐普医疗技术有限公司 | A kind of biological tissue's shear-wave velocity measuring method and medical supersonic wave device |
| CN107510474B (en) * | 2017-09-21 | 2020-07-10 | 深圳开立生物医疗科技股份有限公司 | Shear wave elastic imaging method and system |
| CN107753058B (en) * | 2017-11-22 | 2021-03-02 | 深圳中科乐普医疗技术有限公司 | Shear wave dynamic filtering method |
| CN108852416B (en) * | 2018-05-17 | 2021-03-30 | 深圳开立生物医疗科技股份有限公司 | Method and device for determining propagation speed of shear wave |
| US11406361B2 (en) * | 2018-07-06 | 2022-08-09 | Siemens Healthcare Gmbh | Mapping shear wave velocity and shear modulus in biological tissues |
| CN109589138A (en) * | 2018-11-26 | 2019-04-09 | 深圳中科乐普医疗技术有限公司 | A kind of shear-wave velocity calculation method and elastogram equipment |
| CN111735526B (en) * | 2020-07-20 | 2020-11-17 | 深圳中科乐普医疗技术有限公司 | Ultrasonic elastography device and shear wave velocity measurement method for elasticity measurement |
-
2021
- 2021-04-15 CN CN202110404171.0A patent/CN112799078B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112799078A (en) | 2021-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2688299C1 (en) | Method and device for determining viscoelastic medium parameter | |
| WO2020044769A1 (en) | Ultrasound diagnosis device and ultrasound diagnosis device control method | |
| KR101398948B1 (en) | Viscoelasticity measurement using amplitude-phase modulated ultrasound wave | |
| KR102332464B1 (en) | Shear wave estimation from analytic data | |
| CN104510499B (en) | Shear wave detection in medical ultrasound imaging | |
| JP6216736B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic diagnostic method | |
| CN104622509A (en) | Ultrasonic diagnostic apparatus and elastic evaluation method | |
| RU2641067C1 (en) | Ultrasonic system and method for measurement using shear wave | |
| US20170311929A1 (en) | Method and device for elasticity detection | |
| RU2723753C1 (en) | Method and device for ultrasonic measurement and visualization of biological tissues elasticity in real time | |
| CN106037796A (en) | Quantitative viscoelastic ultrasound imaging | |
| EP2903530B1 (en) | Shear wave attenuation from k-space analysis system | |
| CN110840488B (en) | Imaging method, system and device based on shear wave | |
| KR20170094521A (en) | Sound speed imaging using shear waves | |
| JP4297699B2 (en) | Method and apparatus for rendering spectral distortion | |
| KR101914706B1 (en) | Motion independence in acoustic radiation force impulse imaging | |
| US20190328363A1 (en) | Ultrasound diagnostic apparatus and ultrasound signal processing method | |
| US11540809B2 (en) | Ultrasonic diagnostic apparatus and method for propagation speed analysis of shear wave and elastic modulus measurement of a tissue | |
| CN107427281B (en) | Ultrasonic diagnostic apparatus | |
| CN112799078B (en) | Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment | |
| US9918698B2 (en) | System and method for gradient-based k-space search for shear wave velocity dispersion estimation | |
| KR100998443B1 (en) | Ultrasound System and Method for Processing Ultrasonic Data Considering Scan Conversion | |
| EP2853918B1 (en) | Shear wave detection in medical ultrasound imaging | |
| JP7302651B2 (en) | Ultrasonic signal processing device, ultrasonic diagnostic device, ultrasonic signal processing method, and program | |
| JP2003325520A (en) | Ultrasonic diagnostic apparatus |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |