CN106618635A - Shear wave elastic imaging method and device - Google Patents
Shear wave elastic imaging method and device Download PDFInfo
- Publication number
- CN106618635A CN106618635A CN201710022522.5A CN201710022522A CN106618635A CN 106618635 A CN106618635 A CN 106618635A CN 201710022522 A CN201710022522 A CN 201710022522A CN 106618635 A CN106618635 A CN 106618635A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- image data
- ultrasound image
- shearing wave
- data
- 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.)
- Granted
Links
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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
-
- 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)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Vascular Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention relates to a shear wave elastic imaging method and device. The method comprises the steps of obtaining sampling data of a subject to be detected, wherein the sampling data comprises ultrasonic image data of a rectangular coordinate system; conducting coordinate transformation on the ultrasonic image data of the rectangular coordinate system so as to transforming the ultrasonic image data of the rectangular coordinate system into the ultrasonic image data of a polar coordinate system; conducting treatment on the ultrasonic image data of the polar coordinate system through a preset processing method so as to obtain shear wave elastic imaging used for representing an elasticity modulus of the subject to be detected; conducting visual processing on the shear wave elastic imaging so as to be used for the elasticity modulus of the subject to be detected. The shear wave elastic imaging of the subject to be detected is achieved based on coordinate change, and a shear wave of a long-axis tangent plane and the shear wave of a short-axis tangent plane are subjected to treatment so as to obtain an precise elastic quantitative of the subject to be detected, and to improve pathological detection precision of the subject to be detected.
Description
Technical field
The present invention relates to ultrasonic therapy technical field, relates more specifically to a kind of shearing wave elastograph imaging method and dress
Put.
Background technology
Ultrasonic imaging is widely used for clinic due to having the advantages that real-time, cheap, Noninvasive and Non-ionizing radiation
Diagnosis.Ultrasonic elastograph imaging, is based particularly on the shearing wave imaging of acoustic radiation force, in the qualitative and quantitative measurement of tissue elasticity
Huge effect is played, for example, blood vessel shearing wave elastogram quantitatively plays vital effect to the elasticity of tube wall.But
Due to some characteristics of some tissues itself, the shearing wave elastogram in correlation technique is caused to carry out bullet to it exactly
Property is qualitative and quantitative, by taking blood vessel as an example for, blood vessel has anisotropic characteristic so that in blood vessel shearing wave elastogram
When, need simultaneously consider the shearing wave and the shearing wave of short axis view of the long axis view of blood vessel, the shearing wave bullet in correlation technique
Property imaging be merely capable of the long axis view to blood vessel and carry out shearing wave elastogram, the elasticity of blood vessel can be caused quantitative inaccurate,
Reduce the pathological examination precision of blood vessel.
Accordingly, it is desirable to provide a kind of shearing wave elastograph imaging method and device, to solve the above problems at least in part.
The content of the invention
A kind of shearing wave elastograph imaging method of the invention and device are proposed in view of the problems referred to above, can be directed to be checked
Object is surveyed, such as bodily tissue, blood vessel etc. obtain accurate elasticity quantitatively, improve the pathological examination precision of object to be detected.
According to an aspect of the present invention, there is provided a kind of shearing wave elastograph imaging method, methods described includes：
The sampled data of object to be detected is obtained, the sampled data includes the ultrasound image data of rectangular coordinate system；It is right
The ultrasound image data of the rectangular coordinate system carries out Coordinate Conversion, and the ultrasound image data of the rectangular coordinate system is changed
For the ultrasound image data of polar coordinate system；At default ultrasound image data of the processing method to the polar coordinate system
Reason, to obtain the shearing wave elastogram for characterizing the elastic modelling quantity of the object to be detected；To the shearing wave elasticity into
As carrying out visualization processing, for showing the elastic modelling quantity of the object to be detected.
Alternatively, the ultrasound image data to the rectangular coordinate system carries out Coordinate Conversion includes：To the right angle
The ultrasound image data of coordinate system carries out space interpolation, to obtain the ultrasound image data of polar coordinate system.
Alternatively, the ultrasound image data to the rectangular coordinate system carries out space interpolation includes：It is determined that described treat
The geometric center of detection object；And with the geometric center as polar origin, according to default sample rate respectively to polar coordinates
The angle direction and radial direction of system carries out data sampling, to obtain sampled data points；With the ultrasound figure of the rectangular coordinate system
As data carry out space interpolation to the employing data point, to obtain the ultrasound image data of polar coordinate system.
Alternatively, it is described process is carried out to the ultrasound image data of the polar coordinate system by default processing method to include：
By the ultrasound image data of polar coordinate system described in the first disposal methods, to obtain the ultrasonoscopy number of the polar coordinate system
According to any data point displacement and along the object to be detected interested parties to propagate shearing wave；And, based on described
The displacement of any data point, is processed the shearing wave propagated along the object ring to be detected by second processing method,
To obtain shearing wave dispersion curve.
Alternatively, described process the ultrasound image data of the polar coordinate system by default processing method is also wrapped
Include：Displacement based on any data point and the shearing wave propagated along the object ring to be detected obtain the transmission of shearing wave
Speed.
Alternatively, the ultrasound image data for passing through polar coordinate system described in the first disposal methods includes：To described
The ultrasound image data of polar coordinate system carries out speckle tracking, to obtain the radial direction position in the polar coordinate system of any data point
Move and ring displacement；And displacement time curve is set up based on the radial displacement and ring displacement；It is bent to the displacement time
Line carries out two-dimensional Fourier transform, to obtain Fourier transformation result, wherein, the different angle models of the Fourier transformation result
Enclose the movable information of the different directions for representing shearing wave；Windowing process is carried out to angle interested, after through inverse Fourier become
Change, to obtain the interested parties along the object to be detected to the shearing wave propagated.
Alternatively, the ultrasound image data to the polar coordinate system carries out speckle tracking includes：To the polar coordinates
Any data point in the ultrasound image data of system carries out windowing process, obtain to should data point data window；Connect before and after
Described in the ultrasound image data of continuous two frames during the correlation maximum of data window, determine the data point in the polar coordinate system
Radial displacement and ring displacement.
Alternatively, the displacement based on any data point, and by second processing method to along described to be checked
The shearing wave that survey object ring is propagated carries out process to be included：Based on the displacement of any data point, the shearing wave ripple is obtained
The front displacement time curve obtained in different time in-position；Two-dimensional Fourier transform is carried out to the displacement time curve,
To obtain Fourier transformation result, the Fourier transformation result includes frequency values and the corresponding all wave numbers of arbitrary frequency values
Value；For angle interested, set up based on the corresponding maximum wave number value of each described frequency values in the angle interested
Shearing wave dispersion curve.
Alternatively, it is described visualization processing is carried out to the shearing wave elastogram to include：The polar coordinate system is cut
Cutting ripple elastogram carries out space interpolation, to obtain the visualization data of the shearing wave elastogram.
According to a further aspect of the invention, there is provided shearing wave elastogram device, described device includes：Data acquisition list
Unit, for obtaining the sampled data of object to be detected, the sampled data includes the ultrasound image data of rectangular coordinate system；Coordinate
Converting unit, for carrying out Coordinate Conversion to the ultrasound image data of the rectangular coordinate system, by the rectangular coordinate system
Ultrasound image data is converted to the ultrasound image data of polar coordinate system；Processing unit, for by presetting processing method to described
The ultrasound image data of polar coordinate system is processed, to obtain the shearing wave for characterizing the elastic modelling quantity of the object to be detected
Elastogram；Visualization, it is described to be checked for showing for carrying out visualization processing to the shearing wave elastogram
Survey the elastic modelling quantity of object.
Shearing wave elastograph imaging method provided in an embodiment of the present invention and device, based on coordinate transform, (rectangular coordinate system turns
It is changed to polar coordinate system) the shearing wave elastogram of object to be detected is realized, the shearing of the long axis view of detection object can be treated
The shearing wave of ripple and short axis view is processed, quantitative with the accurate elasticity for obtaining object to be detected, is improved to be detected right
The pathological examination precision of elephant.
Description of the drawings
The embodiment of the present invention is described in more detail by combining accompanying drawing, above-mentioned and other purposes of the present invention,
Feature and advantage will be apparent from.Accompanying drawing is used for providing further understanding the embodiment of the present invention, and constitutes explanation
A part for book, is used to explain the present invention together with the embodiment of the present invention, is not construed as limiting the invention.In the accompanying drawings,
Identical reference number typically represents same parts or step.
Fig. 1 illustrates the schematic flow sheet of shearing wave elastograph imaging method according to an embodiment of the invention；
Fig. 2 illustrates that space interpolation according to an embodiment of the invention realizes the schematic diagram of process；
Fig. 3 illustrates the principle schematic of speckle tracking according to an embodiment of the invention；
Fig. 4 illustrates the structured flowchart of shearing wave elastogram device according to an embodiment of the invention.
Specific embodiment
In order that the object, technical solutions and advantages of the present invention become apparent from, root is described in detail below with reference to accompanying drawings
According to the example embodiment of the present invention.Obviously, described embodiment is only a part of embodiment of the present invention, rather than this
Bright whole embodiments, it should be appreciated that the present invention is not limited by example embodiment described herein.Described in the present invention
The embodiment of the present invention, those skilled in the art's all other embodiment resulting in the case where creative work is not paid
All should fall under the scope of the present invention.
It is a kind of Ultrasonic Elasticity Imaging of assessment tissue hardness based on the shearing wave elastography of acoustic radiation force.
Its general principle is：Launch the ultrasonic wave of high-energy from probe to the soft tissue of organism, acoustic radiation force is produced, in acoustic radiation force
In the presence of the shear stress of tissue, the soft tissue in specific region can produce the vibration propagated to surrounding, cut so as to produce
Ripple is cut, because the hardness and shear wave velocity of organism soft tissue have incidence relation, therefore can be by detecting shearing wave
Speed analyzing the hardness of organism soft tissue.
The shearing wave elastogram of correlation technique can include two key points：1st, axial displacement, i.e. shearing wave are considered
Direction of vibration；2nd, some direction interested is only considered when anisotropic filter is arranged.For the soft tissue as blood vessel,
Shearing wave elastogram in correlation technique is only capable of processing the shearing wave imaging of its long axis view, and it is short to process its
The shearing wave elastogram of axial section.But for the shearing wave elastogram of blood vessel short axis view, it is envisaged that radial direction position
Move rather than axial displacement, and need the anisotropic filter along ring filtering.Therefore, by the shearing wave elasticity of correlation technique
Imaging, it is difficult to carry out elasticity to blood vessel exactly quantitatively, and reduces the pathological examination precision of blood vessel.
In order to solve problem as described above, the embodiment of the present invention proposes a kind of shearing wave elastograph imaging method and device,
The shearing wave elastogram of object to be detected is realized based on coordinate transform, the accurate elasticity for obtaining object to be detected is quantitative, carries
Height treats the pathological examination precision of detection object.
The shearing wave elastograph imaging method and device provided the present invention below in conjunction with accompanying drawing is described in detail, so that
Those skilled in the art can clearly and accurately understand technical scheme.
Fig. 1 illustrates the schematic flow sheet of shearing wave elastograph imaging method according to an embodiment of the invention.Such as Fig. 1 institutes
Show, a kind of shearing wave elastograph imaging method provided in an embodiment of the present invention comprises the steps：
Step 110, obtains the sampled data of object to be detected, wherein, sampled data includes the ultrasound figure of rectangular coordinate system
As data.
In this step, it is possible to use acoustic radiation force encourages the area-of-interest of object to be detected, and area-of-interest is entered
Row ultrahigh speed ultrasonic imaging, to obtain the ultrasound image data of the area-of-interest of object to be detected.Exemplarily, ultrasonoscopy
Data can include at least one of RF (Radio Frequency, radio frequency) data, envelope data, B moduluses according in.
By taking blood vessel as an example, internal elastic imaging can encourage blood vessel using air bag, blood pressure or external compression,
The motion for estimating blood vessel is displacement (generally longitudinal), the Strain Distribution of blood vessel is obtained, so as to characterize the elasticity of blood vessel.At this
In bright embodiment, it is possible to use acoustic radiation force encourages blood vessel, obtains the ultrasound image data of blood vessel.Exemplarily, in this step
The ultrasound image data of middle acquisition is the ultrasound image data of rectangular coordinate system.
Step 120, to the ultrasound image data of rectangular coordinate system Coordinate Conversion is carried out, by the ultrasound figure of rectangular coordinate system
As data are converted to the ultrasound image data of polar coordinate system.
In this step, the ultrasound image data of rectangular coordinate system can be converted to by polar coordinate system by coordinate transform
Ultrasound image data.The Coordinate Conversion of rectangular coordinate system to polar coordinate system can be realized by Spatial Interpolation Method.Show as one
Example, can be defined as polar origin, according to default sample rate respectively to polar coordinate system by the geometric center of object to be detected
Angle direction and radial direction carry out data sampling, to obtain sampled data points, with the ultrasound image data pair of rectangular coordinate system
Space interpolation is carried out using data point, to obtain the ultrasound image data of polar coordinate system.Exemplarily, it is horizontal in polar coordinate system
Coordinate is angle, and ordinate is radius.Geometric center with regard to determining object to be detected, the side that for example can manually select
Formula determines, it is also possible to carry out Image Automatic Segmentation to realize the selection to its geometric center by treating the image of detection object.
Still by taking blood vessel as an example, by the geometric center (for example, the center of circle) of blood vessel as polar origin, from origin wire is launched
Sampling, wherein, it is radial direction along directions of rays, it is angle direction perpendicular to directions of rays.In embodiments of the present invention, root
The radial direction and angle direction of blood vessel are sampled according to default sample rate, to obtain data point is adopted, recycle right angle
The ultrasound image data of coordinate system realizes that rectangular co-ordinate is tied to the coordinate of polar coordinate system and turns to carrying out space interpolation using data point
Change, so as to obtain the ultrasound image data of polar coordinate system.Default sample rate for example can be that any one direction of polar coordinate system is adopted
The quantity of the data point of collection, it can set according to the actual requirements, this is not defined in embodiments of the present invention.Insert in space
Value method can not be limited this for example with any one realization such as linear interpolation, fitting of a polynomial interpolation, the present invention.
The schematic diagram of process is realized below in conjunction with the space interpolation of the embodiment of the present invention shown in Fig. 2, it is possible to understand that
It is that the interpolation embodiment shown in Fig. 2 is merely illustrative, and does not cause limitation of the present invention.It is with two-dimensional space interpolation in Fig. 2
Example, if the coordinate of any one data point in rectangular coordinate system can be expressed as (x, y), it is determined that in the geometry of object to be detected
After the heart (for example, the center of circle of blood vessel), it can be (r, θ) that any one data point is represented in polar coordinates.Point solid in figure is represented
The ultrasound image data of rectangular coordinate system, hollow point represents the data (polar coordinates for needing to obtain under polar coordinates after the center of circle determines
Ultrasound image data under system), can be obtained by pair radius direction and angle direction uniform sampling.Can be according to solid dot
Information obtains the information of hollow point by space interpolation.
Step 130, is processed the ultrasound image data of polar coordinate system, to obtain for table by default processing method
Levy the shearing wave elastogram of the elastic modelling quantity of object to be detected.
Shearing wave elastogram is that the propagation for producing shearing wave by detecting acoustic radiation force to excite tissue is imaged
Method.Can be in the hope of the modulus of shearing of tissue, to enter to histopathology by carrying out calculation process to shearing wave elastogram
Row analysis.
In this step, at can be by ultrasound image data of the processing method set in advance to polar coordinate system
Reason, so as to obtain the shearing wave elastogram of object to be detected.According to one embodiment of the invention, the first processing method can be passed through
The ultrasound image data of polar coordinate system is processed, to obtain displacement and the edge of any data point of the ultrasound image data of polar coordinate system
The interested parties of the object to be detected are to the shearing wave propagated, and are displaced through second processing side based on any data point
Method is processed the shearing wave propagated along object ring to be detected, to obtain shearing wave dispersion curve.
Further, in certain embodiments, the place using the first processing method to the ultrasound image data of polar coordinate system
Reason may include steps of：
Step A, the ultrasound image data to polar coordinate system carry out speckle tracking, with obtain any data point in polar coordinates
Radial displacement and ring displacement in system.
Displacement Estimation can be realized by speckle tracking in this step, wherein, the implementation method of speckle tracking includes
All methods that can obtain two-dimension displacement in polar coordinates of the Block matchings such as cross-correlation, auto-correlation, optical flow method are not limited.
Exemplarily, windowing process can be carried out to the point of any data in the ultrasound image data of polar coordinate system, is obtained
To should data point data window, when in the ultrasound image data of two continuous frames data window correlation (evaluate correlation ginseng
Number is many, for example, it may be auto-correlation coefficient, or cross-correlation coefficient etc.) it is maximum when, determine sitting in pole for the data point
Radial displacement and ring displacement in mark system.In embodiments of the present invention, relative coefficient can by auto-correlation, cross-correlation,
Absolute error and, the mode such as amplitude difference quadratic sum and optical flow method calculate determination, it will be understood that the calculating of relative coefficient is not limited
In this several implementation, can also realize that this is no longer going to repeat them by the motion estimation algorithm of other Block- matchings.As schemed
The principle schematic of the speckle tracking shown in 3.In a images of Fig. 3, determine an area-of-interest (such as data window position in figure
Put), because shearing wave is continuous motion, it is therefore desirable to tracking is scanned for it at least two field pictures, in Fig. 3
B and c images, (matching degree is can be understood as best when the correlation maximum of data window in the ultrasound image data of two continuous frames
Position, i.e. best match position in Fig. 3 in d images), respectively obtain the displacement along angle direction and radial direction, i.e. footpath
To displacement and ring displacement, the mark of the d images in such as Fig. 3 two at displacement.
Step B, displacement time curve is set up based on radial displacement and ring displacement.
Exemplarily, for example can be on the propagation path of the area-of-interest of object to be detected, any one moment pair
The locus for answering wavefront to reach, continuous several moment are combined, and a dimension represents the time, another dimension
Locus is represented, because front propagation is continuous, therefore the situation of front propagation can be determined in two-dimensional space, that is, built
Vertical displacement time curve.
Step C, two-dimensional Fourier transform is carried out to displacement time curve, to obtain Fourier transformation result, wherein, Fu Li
The different angular ranges of leaf transformation result represent the movable information of the different directions of shearing wave.
Step D, windowing process is carried out to angle interested, after through inverse Fourier transform, to obtain along to be detected right
The interested parties of elephant are to the shearing wave propagated.
Step C and step D realize trend pass filtering, i.e., displacement time curve is carried out Fourier transformation, inverse Fourier transform,
Direction windowing process obtains interested parties to the shearing wave propagated.
And step 130 based on any data point be displaced through second processing method to along object ring to be detected propagate
Shearing wave can include the step of processed：Based on the displacement of any data point, shearing wave wavefront is obtained in different time
The displacement time curve that in-position obtains, to displacement time curve two-dimensional Fourier transform is carried out, to obtain Fourier transformation
As a result, Fourier transformation result includes frequency values and the corresponding all wave number values of arbitrary frequency values, for angle interested,
Shearing wave dispersion curve is set up based on the corresponding maximum wave number value of each frequency values in angle interested.
Exemplarily, select angle direction interested, for example, when angle direction interested to the right when, corresponding Fu
In information in leaf transformation be located at first quartile, wherein, the information in Fourier transformation includes that respectively frequency information and wave number are believed
Breath (frequency values and wave number value), takes wave number value maximum under any one frequency values, and the shearing wave frequency dispersion for obtaining object to be detected is bent
The shearing phase velocity of wave curve of line, i.e. different frequency.Dispersion curve based on shearing wave evaluates the modulus of shearing of object to be detected,
Accuracy rate is higher, and for object to be detected as blood vessel, based under polar coordinate system to the selection of angle interested, can
To process long axis view and short axis view simultaneously, the drawbacks of make up correlation technique, the more accurately elasticity of object to be detected is obtained
Quantitatively.
According to another embodiment of the present invention, further, this step can also include the displacement based on any data point and
The shearing wave propagated along object ring to be detected obtains the transmission speed of shearing wave.Exemplarily, in polar coordinate system, according to cutting
Cut ripple wavefront to measure shear-wave velocity in different time in-position, so as to obtain the spread speed of shearing wave.Using
The embodiment, with reference to the spread speed of shearing wave and the dispersion curve of shearing wave the modulus of shearing of object to be detected is evaluated, both
Can complement each other, further improve the elastic quantitative precision of object to be detected.
Step 140, to shearing wave elastogram visualization processing is carried out, for showing the springform of object to be detected
Amount.
Because Displacement Estimation, trend pass filtering etc. are carried out under polar coordinates, the propagation of the shearing wave after trend pass filtering is obtained
Afterwards, coordinate can be carried out to the data (for example, the sub- vibration displacement of scattering, speed, acceleration etc. can be included) that shearing wave is propagated
Conversion, i.e., be converted into rectangular co-ordinate by polar coordinates, and transfer principle and process are identical with step 120, will not be described here, by sitting
The visualization that shearing wave is propagated is realized in mark conversion, that is, realize the visualization processing of shearing wave elastogram.
Using shearing wave elastograph imaging method provided in an embodiment of the present invention, based on coordinate transform, (rectangular coordinate system is changed
For polar coordinate system) the shearing wave elastogram of object to be detected is realized, the shearing wave of the long axis view of detection object can be treated
Processed with the shearing wave of short axis view, quantitative with the accurate elasticity for obtaining object to be detected, detection object is treated in raising
Pathological examination precision.
Additionally, by trend pass filtering process, fortune of the shearing wave in object to be detected (for example, vascular wall) radial direction can be obtained
It is dynamic, it is possible to effectively to remove the back wave of object to be detected (for example, along vascular wall) ring, it is that the elasticity of object to be detected is determined
Quantity research provides basis.
In addition, being imaged by the shearing wave of object to be detected (such as blood vessel) cross section, tube wall can be further obtained
Shear wave velocity, dispersion curve, the elasticity for treating detection object with reference to shearing velocity of wave propagation and dispersion curve is quantitatively surveyed
Amount, can improve certainty of measurement and the degree of accuracy.
The present invention also provides a kind of shearing wave elastogram device using above-mentioned shearing wave elastograph imaging method.Fig. 4 shows
Go out the structured flowchart of shearing wave elastogram device according to an embodiment of the invention.As shown in figure 4, shearing wave elastogram
Device can include data capture unit 410, coordinate transformation unit 420, processing unit 430 and visualization 440.It is exemplary
Ground, data capture unit 410, coordinate transformation unit 420, processing unit 430 and visualization 440 can be arranged on calculating
In machine, realized by the operation processing unit of computer.
Data capture unit 410 can be used for obtaining the sampled data of object to be detected, wherein, sampled data includes right angle
The ultrasound image data of coordinate system.
Coordinate transformation unit 420 can be used for the ultrasonoscopy number of the rectangular coordinate system obtained to data capture unit 410
According to Coordinate Conversion is carried out, the ultrasound image data of super rectangular coordinate system is converted to into the ultrasound image data of polar coordinate system.
Processing unit 430 can be used for by presetting the polar coordinate system that processing method is obtained to coordinate transformation unit 420
Ultrasound image data is processed, to obtain the shearing wave elastogram for characterizing the elastic modelling quantity of object to be detected.
Visualization 440 can be used for carrying out visualization processing to the shearing wave elastogram that processing unit 430 is obtained,
For showing the elastic modelling quantity of the object to be detected.
Using shearing wave elastogram device provided in an embodiment of the present invention, coordinate transform (rectangular coordinate system can be based on
Be converted to polar coordinate system) realize the shearing wave elastogram of object to be detected, the long axis view that can treat detection object is cut
The shearing wave for cutting ripple and short axis view is processed, quantitative with the accurate elasticity for obtaining object to be detected, is improved to be detected
The pathological examination precision of object.
Although the herein example embodiment by reference to Description of Drawings, it should be understood that above-mentioned example embodiment is merely exemplary
, and be not intended to limit the scope of the invention to this.Those of ordinary skill in the art can wherein carry out various changes
And modification, it is made without departing from the scope of the present invention and spirit.All such changes and modifications are intended to be included in claims
Within required the scope of the present invention.
Those of ordinary skill in the art are it is to be appreciated that the list of each example with reference to the embodiments described herein description
Unit and algorithm steps, being capable of being implemented in combination in electronic hardware or computer software and electronic hardware.These functions are actually
Performed with hardware or software mode, depending on the application-specific and design constraint of technical scheme.Professional and technical personnel
Each specific application can be used different methods to realize described function, but this realization it is not considered that exceeding
The scope of the present invention.
In several embodiments provided herein, it should be understood that disclosed apparatus and method, it can be passed through
Its mode is realized.For example, apparatus embodiments described above are only schematic, and for example, the division of unit is only
A kind of division of logic function, can there is an other dividing mode when actually realizing, such as multiple units or component can with reference to or
Person is desirably integrated into another equipment, or some features can be ignored, or does not perform.
In specification mentioned herein, a large amount of details are illustrated.It is to be appreciated, however, that the enforcement of the present invention
Example can be put into practice in the case of without these details.In some instances, known method, structure is not been shown in detail
And technology, so as not to obscure the understanding of this description.
Similarly, it will be appreciated that in order to simplify the present invention and help understand one or more in each inventive aspect, exist
To the present invention exemplary embodiment description in, the present invention each feature be grouped together into sometimes single embodiment, figure,
Or in descriptions thereof.However, the method for the present invention should be construed to reflect following intention：It is i.e. required for protection
The more features of feature that application claims ratio is expressly recited in each claim.More precisely, such as corresponding power
As sharp claim reflects, its inventive point is can be with the spy of all features less than certain disclosed single embodiment
Levy to solve corresponding technical problem.Therefore, it then follows it is concrete that thus claims of specific embodiment are expressly incorporated in this
Separate embodiments of the embodiment, wherein each claim as the present invention itself.
It will be understood to those skilled in the art that in addition to mutually exclusive between feature, any combinations pair can be adopted
All features and so disclosed any method disclosed in this specification (including adjoint claim, summary and accompanying drawing)
Or all processes or unit of equipment are combined.Unless expressly stated otherwise, this specification (will including adjoint right
Ask, make a summary and accompanying drawing) disclosed in each feature can, equivalent identical by offer or similar purpose alternative features replacing.
Although additionally, it will be appreciated by those of skill in the art that some embodiments described herein include other embodiments
In included some features rather than further feature, but the combination of the feature of different embodiments means in of the invention
Within the scope of and form different embodiments.For example, in detail in the claims, embodiment required for protection one of arbitrarily
Can in any combination mode using.
The present invention all parts embodiment can be realized with hardware, or with one or more processor operation
Software module realize, or with combinations thereof realization.It will be understood by those of skill in the art that can use in practice
During microprocessor or digital signal processor (DSP) are to realize wire suspension location device according to embodiments of the present invention
The some or all functions of some modules.The present invention is also implemented as the part for performing method as described herein
Or the program of device (for example, computer program and computer program) of whole.Such program for realizing the present invention can
To store on a computer-readable medium, or there can be the form of one or more signal.Such signal can be from
Download on internet website and obtain, or provide on carrier signal, or provide in any other form.
It should be noted that above-described embodiment the present invention will be described rather than limits the invention, and ability
Field technique personnel can design without departing from the scope of the appended claims alternative embodiment.In the claims,
Any reference symbol between bracket should not be configured to limitations on claims.Word "comprising" is not excluded the presence of not
Component listed in the claims or step.Word "a" or "an" before component does not exclude the presence of multiple such
Component.The present invention can come real by means of the hardware for including some different components and by means of properly programmed computer
It is existing.If in the unit claim for listing equipment for drying, several in these devices can be by same hardware branch
To embody.The use of word first, second, and third does not indicate that any order.These words can be explained and be run after fame
Claim.
The above, the only specific embodiment of the present invention or the explanation to specific embodiment, the protection of the present invention
Scope is not limited thereto, any those familiar with the art the invention discloses technical scope in, can be easily
Expect change or replacement, all should be included within the scope of the present invention.Protection scope of the present invention should be with claim
Protection domain is defined.
Claims (10)
1. a kind of shearing wave elastograph imaging method, it is characterised in that methods described includes：
The sampled data of object to be detected is obtained, the sampled data includes the ultrasound image data of rectangular coordinate system；
Coordinate Conversion is carried out to the ultrasound image data of the rectangular coordinate system, by the ultrasonoscopy number of the rectangular coordinate system
According to the ultrasound image data for being converted to polar coordinate system；
The ultrasound image data of the polar coordinate system is processed by default processing method, to obtain described treating is characterized
The shearing wave elastogram of the elastic modelling quantity of detection object；
Visualization processing is carried out to the shearing wave elastogram, for showing the elastic modelling quantity of the object to be detected.
2. method according to claim 1, it is characterised in that the ultrasound image data to the rectangular coordinate system enters
Row Coordinate Conversion includes：
Space interpolation is carried out to the ultrasound image data of the rectangular coordinate system, to obtain the ultrasound image data of polar coordinate system.
3. method according to claim 2, it is characterised in that the ultrasound image data to the rectangular coordinate system enters
Row space interpolation includes：
Determine the geometric center of the object to be detected；And
With the geometric center as polar origin, according to default sample rate angle direction respectively to polar coordinate system and radius side
To data sampling is carried out, to obtain sampled data points；
Space interpolation is carried out to the employing data point with the ultrasound image data of the rectangular coordinate system, to obtain polar coordinate system
Ultrasound image data.
4. method according to claim 1, it is characterised in that it is described by default processing method to the polar coordinate system
Ultrasound image data carries out process to be included：
By the ultrasound image data of polar coordinate system described in the first disposal methods, to obtain the ultrasound figure of the polar coordinate system
Displacement as any data point of data and the interested parties along the object to be detected are to the shearing wave propagated；
And,
Based on the displacement of any data point, by shearing of the second processing method to propagating along the object ring to be detected
Ripple is processed, to obtain shearing wave dispersion curve.
5. method according to claim 4, it is characterised in that it is described by default processing method to the polar coordinate system
Ultrasound image data is processed also to be included：
Displacement based on any data point and the shearing wave propagated along the object ring to be detected obtain the biography of shearing wave
Defeated speed.
6. method according to claim 4, it is characterised in that described to pass through polar coordinate system described in the first disposal methods
Ultrasound image data include：
Speckle tracking is carried out to the ultrasound image data of the polar coordinate system, with obtain any data point in the polar coordinate system
In radial displacement and ring displacement；And
Displacement time curve is set up based on the radial displacement and ring displacement；
Two-dimensional Fourier transform is carried out to the displacement time curve, to obtain Fourier transformation result, wherein, the Fourier
The different angular ranges of transformation results represent the movable information of the different directions of shearing wave；
Windowing process is carried out to angle interested, after through inverse Fourier transform, to obtain the sense along the object to be detected
The shearing wave that interest direction is propagated.
7. method according to claim 6, it is characterised in that the ultrasound image data to the polar coordinate system is carried out
Speckle tracking includes：
Windowing process is carried out to any data point in the ultrasound image data of the polar coordinate system, is obtained to should data point
Data window；
When the relative coefficient of the data window described in the ultrasound image data of two continuous frames is maximum, determine the data point in institute
State the radial displacement and ring displacement in polar coordinate system.
8. method according to claim 4, it is characterised in that the displacement based on any data point, Yi Jitong
Crossing second processing method and carrying out process to the shearing wave propagated along the object ring to be detected includes：
Based on the displacement of any data point, when obtaining the displacement that the shearing wave wavefront is obtained in different time in-position
Half interval contour；
Two-dimensional Fourier transform is carried out to the displacement time curve, to obtain Fourier transformation result, the Fourier transformation
As a result frequency values and the corresponding all wave number values of arbitrary frequency values are included；
For angle interested, set up based on the corresponding maximum wave number value of each described frequency values in the angle interested
Shearing wave dispersion curve.
9. method according to claim 1, it is characterised in that described that the shearing wave elastogram is carried out at visualization
Reason includes：
Space interpolation is carried out to the shearing wave elastogram of the polar coordinate system, to obtain the visual of the shearing wave elastogram
Change data.
10. a kind of shearing wave elastogram device, it is characterised in that described device includes：
Data capture unit, for obtaining the sampled data of object to be detected, the sampled data includes the super of rectangular coordinate system
Acoustic image data；
Coordinate transformation unit, for carrying out Coordinate Conversion to the ultrasound image data of the rectangular coordinate system, will be described super straight
The ultrasound image data of angular coordinate system is converted to the ultrasound image data of polar coordinate system；
Processing unit, for being processed the ultrasound image data of the polar coordinate system by default processing method, to obtain
For characterizing the shearing wave elastogram of the elastic modelling quantity of the object to be detected；
Visualization is described to be detected right for showing for carrying out visualization processing to the shearing wave elastogram
The elastic modelling quantity of elephant.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710022522.5A CN106618635B (en) | 2017-01-12 | 2017-01-12 | Shearing wave elastograph imaging method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710022522.5A CN106618635B (en) | 2017-01-12 | 2017-01-12 | Shearing wave elastograph imaging method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106618635A true CN106618635A (en) | 2017-05-10 |
| CN106618635B CN106618635B (en) | 2019-11-08 |
Family
ID=58843177
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710022522.5A Active CN106618635B (en) | 2017-01-12 | 2017-01-12 | Shearing wave elastograph imaging method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106618635B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107577642A (en) * | 2017-08-25 | 2018-01-12 | 深圳中科乐普医疗技术有限公司 | A kind of shear-wave direction filter achieving method and medical supersonic wave device |
| CN107970043A (en) * | 2017-12-28 | 2018-05-01 | 深圳开立生物医疗科技股份有限公司 | A kind of detection method and device of shearing wave |
| CN109567872A (en) * | 2018-11-05 | 2019-04-05 | 清华大学 | Blood vessel guided wave elastograph imaging method and system based on machine learning |
| CN110420036A (en) * | 2019-08-01 | 2019-11-08 | 无锡海斯凯尔医学技术有限公司 | Polymorphic type probe compatible data processing method, device and readable storage medium storing program for executing |
| KR20200030580A (en) * | 2017-07-21 | 2020-03-20 | 우시 히스키 메디칼 테크놀로지스 컴퍼니., 리미티드. | Method and apparatus for measuring viscoelasticity of medium |
| CN110897656A (en) * | 2019-12-12 | 2020-03-24 | 首都医科大学附属北京友谊医院 | Ultrasonic shear wave elasticity measuring method and device |
| CN111449681A (en) * | 2020-04-08 | 2020-07-28 | 深圳开立生物医疗科技股份有限公司 | Shear wave imaging method, device and equipment and readable storage medium |
| CN111671464A (en) * | 2020-06-19 | 2020-09-18 | 深圳开立生物医疗科技股份有限公司 | Image marking method, device, equipment and readable storage medium |
| CN114173672A (en) * | 2020-04-24 | 2022-03-11 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic elastography method, ultrasonic frequency dispersion measurement device and storage medium |
| CN114402197A (en) * | 2019-05-21 | 2022-04-26 | P & L管道测量有限责任公司 | Apparatus and method for pipe inspection using shear waves generated by EMAT |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100010346A1 (en) * | 2008-04-07 | 2010-01-14 | Greenleaf James F | System and Method For Non-Invasively Measuring Tissue Viscoelasticity Using Surface Waves |
| CN102283679A (en) * | 2011-08-04 | 2011-12-21 | 中国科学院深圳先进技术研究院 | Ultrasonic imaging system for elasticity measurement and method for measuring elasticity of biological tissue |
| CN102805616A (en) * | 2011-06-02 | 2012-12-05 | 曹铁生 | Novel index for testing partial elasticity of artery |
| CN102961166A (en) * | 2011-08-31 | 2013-03-13 | 通用电气公司 | Method for detecting and tracing needle |
| CN104398271A (en) * | 2014-11-14 | 2015-03-11 | 西安交通大学 | Method using three-dimensional mechanics and tissue specific imaging of blood vessels and plaques for detection |
| US20160157828A1 (en) * | 2014-06-05 | 2016-06-09 | Chikayoshi Sumi | Beamforming method, measurement and imaging instruments, and communication instruments |
-
2017
- 2017-01-12 CN CN201710022522.5A patent/CN106618635B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100010346A1 (en) * | 2008-04-07 | 2010-01-14 | Greenleaf James F | System and Method For Non-Invasively Measuring Tissue Viscoelasticity Using Surface Waves |
| CN102805616A (en) * | 2011-06-02 | 2012-12-05 | 曹铁生 | Novel index for testing partial elasticity of artery |
| CN102283679A (en) * | 2011-08-04 | 2011-12-21 | 中国科学院深圳先进技术研究院 | Ultrasonic imaging system for elasticity measurement and method for measuring elasticity of biological tissue |
| CN102961166A (en) * | 2011-08-31 | 2013-03-13 | 通用电气公司 | Method for detecting and tracing needle |
| US20160157828A1 (en) * | 2014-06-05 | 2016-06-09 | Chikayoshi Sumi | Beamforming method, measurement and imaging instruments, and communication instruments |
| CN104398271A (en) * | 2014-11-14 | 2015-03-11 | 西安交通大学 | Method using three-dimensional mechanics and tissue specific imaging of blood vessels and plaques for detection |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2020527723A (en) * | 2017-07-21 | 2020-09-10 | 无錫海斯凱尓医学技術有限公司Ｗｕｘｉ Ｈｉｓｋｙ Ｍｅｄｉｃａｌ Ｔｅｃｈｎｏｌｏｇｉｅｓ Ｃｏ．，Ｌｔｄ． | Method and device for quantifying viscoelasticity of medium |
| KR20200030580A (en) * | 2017-07-21 | 2020-03-20 | 우시 히스키 메디칼 테크놀로지스 컴퍼니., 리미티드. | Method and apparatus for measuring viscoelasticity of medium |
| US11719613B2 (en) | 2017-07-21 | 2023-08-08 | Wuxi Hisky Medical Technologies Co., Ltd. | Method and device for quantifying viscoelasticity of a medium |
| JP7006888B2 (en) | 2017-07-21 | 2022-01-24 | 无錫海斯凱尓医学技術有限公司 | Method and device for quantifying viscoelasticity of medium |
| KR102317869B1 (en) | 2017-07-21 | 2021-10-25 | 우시 히스키 메디칼 테크놀로지스 컴퍼니., 리미티드. | Method and apparatus for quantifying the viscoelasticity of a medium |
| CN107577642A (en) * | 2017-08-25 | 2018-01-12 | 深圳中科乐普医疗技术有限公司 | A kind of shear-wave direction filter achieving method and medical supersonic wave device |
| CN107577642B (en) * | 2017-08-25 | 2020-06-23 | 深圳中科乐普医疗技术有限公司 | Shear wave direction filtering implementation method and medical ultrasonic equipment |
| CN107970043A (en) * | 2017-12-28 | 2018-05-01 | 深圳开立生物医疗科技股份有限公司 | A kind of detection method and device of shearing wave |
| CN109567872A (en) * | 2018-11-05 | 2019-04-05 | 清华大学 | Blood vessel guided wave elastograph imaging method and system based on machine learning |
| CN114402197A (en) * | 2019-05-21 | 2022-04-26 | P & L管道测量有限责任公司 | Apparatus and method for pipe inspection using shear waves generated by EMAT |
| CN110420036A (en) * | 2019-08-01 | 2019-11-08 | 无锡海斯凯尔医学技术有限公司 | Polymorphic type probe compatible data processing method, device and readable storage medium storing program for executing |
| CN110897656A (en) * | 2019-12-12 | 2020-03-24 | 首都医科大学附属北京友谊医院 | Ultrasonic shear wave elasticity measuring method and device |
| CN110897656B (en) * | 2019-12-12 | 2022-06-21 | 首都医科大学附属北京友谊医院 | Ultrasonic shear wave elasticity measuring method and device |
| CN111449681A (en) * | 2020-04-08 | 2020-07-28 | 深圳开立生物医疗科技股份有限公司 | Shear wave imaging method, device and equipment and readable storage medium |
| CN111449681B (en) * | 2020-04-08 | 2023-09-08 | 深圳开立生物医疗科技股份有限公司 | Shear wave imaging method, device, equipment and readable storage medium |
| CN114173672A (en) * | 2020-04-24 | 2022-03-11 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic elastography method, ultrasonic frequency dispersion measurement device and storage medium |
| CN111671464A (en) * | 2020-06-19 | 2020-09-18 | 深圳开立生物医疗科技股份有限公司 | Image marking method, device, equipment and readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106618635B (en) | 2019-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106618635A (en) | Shear wave elastic imaging method and device | |
| CN107049360B (en) | Shear wave elastography method and device | |
| KR101984824B1 (en) | Method and apparatus for analyzing elastography of tissue using ultrasound | |
| US11517289B2 (en) | Method and apparatus for acquiring motion information | |
| CN114972567B (en) | Medical ultrasonic CT multi-parameter image reconstruction method based on wave equation | |
| CN106983524A (en) | A kind of parameter and its measuring method for reflecting that biological tissue is abnormal | |
| TWI703962B (en) | Medium viscoelasticity quantitative method and device | |
| Engel et al. | A new method for shear wave speed estimation in shear wave elastography | |
| Kijanka et al. | Fast local phase velocity-based imaging: Shear wave particle velocity and displacement motion study | |
| CN111789632B (en) | Ultrasonic diagnostic apparatus, signal processing apparatus, and storage medium | |
| JP2022542356A (en) | Ultrasonic signal processing method, device, equipment and storage medium | |
| Tran et al. | 2D shear wave imaging in Gaussian noise and reflection media | |
| WO2018164181A1 (en) | Tissue elasticity measurement device and measurement method | |
| Lasaygues et al. | Circular antenna for breast ultrasonic diffraction tomography | |
| Montagnon et al. | Viscoelastic characterization of elliptical mechanical heterogeneities using a semi-analytical shear-wave scattering model for elastometry measures | |
| Carbente et al. | Image reconstruction utilizing median filtering applied to elastography | |
| Grube et al. | Ultrasound shear wave velocity estimation in a small field of view via spatio-temporal deep learning | |
| RU2795207C1 (en) | Method and device for processing ultrasonic signals, device and information media | |
| Quang-Huy et al. | Two-dimensional viscoelastic imaging using an enhanced FDTD-AHI approach | |
| Kauati et al. | Mean scatterer space estimation from ultrasound signals combining singular spectral analysis and entropy | |
| CN108427655B (en) | Local frequency estimation method based on lognormal filter | |
| Curry | Prostate Tissue Motion Tracking and Ultrasound Elastography Using Tissue Mechanics Derived Constraints | |
| Zwiebel | Limited view sound speed imaging for breast cancer detection | |
| Hewener et al. | Deconvolution of Medical Ultrasound Data with Consideration of the Pressure Field, the Excitation Pulse and Focussing | |
| KR20190060740A (en) | Method and apparatus for analyzing elastography of tissue using ultrasound |
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 |