CN106373166A - Multi-shape plane shear wave composite imaging method based on Mach cone effects - Google Patents
Multi-shape plane shear wave composite imaging method based on Mach cone effects Download PDFInfo
- Publication number
- CN106373166A CN106373166A CN201610893496.9A CN201610893496A CN106373166A CN 106373166 A CN106373166 A CN 106373166A CN 201610893496 A CN201610893496 A CN 201610893496A CN 106373166 A CN106373166 A CN 106373166A
- Authority
- CN
- China
- Prior art keywords
- shearing
- shear
- tissue
- shearing wave
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 30
- 230000000694 effects Effects 0.000 title claims abstract description 28
- 239000002131 composite material Substances 0.000 title claims abstract description 19
- 238000006073 displacement reaction Methods 0.000 claims abstract description 36
- 230000005855 radiation Effects 0.000 claims abstract description 15
- 239000000523 sample Substances 0.000 claims abstract description 12
- 230000001427 coherent effect Effects 0.000 claims abstract description 7
- 238000010008 shearing Methods 0.000 claims description 141
- 238000000034 method Methods 0.000 claims description 14
- 238000012544 monitoring process Methods 0.000 claims description 12
- 230000033228 biological regulation Effects 0.000 claims description 11
- 230000006698 induction Effects 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 9
- 230000006854 communication Effects 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 7
- 230000005284 excitation Effects 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 6
- 230000001737 promoting effect Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 4
- 238000002366 time-of-flight method Methods 0.000 claims description 4
- 230000021615 conjugation Effects 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 3
- 230000008520 organization Effects 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 241001274660 Modulus Species 0.000 claims 1
- 230000000149 penetrating effect Effects 0.000 claims 1
- 230000001105 regulatory effect Effects 0.000 abstract description 4
- 210000001519 tissue Anatomy 0.000 description 42
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 230000001276 controlling effect Effects 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 208000009453 Thyroid Nodule Diseases 0.000 description 1
- 201000008275 breast carcinoma Diseases 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 201000001514 prostate carcinoma Diseases 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention provides a multi-shape plane shear wave composite imaging method based on Mach cone effects. Through regulating an acoustic radiation force to excite a probe to be quickly focused on different longitudinal depths of a tissue, a quickly-moving shear wave source is formed, coherent interference happens to shear waves generated by multiple wave sources due to Mach cone effects, a large-range shear wave plane moving to two sides is formed, the shear wave propagation condition is monitored in an ultra high-speed and ultrasonic mode, the shear wave speed is predicted, a tissue shear modulus is inverted, and a tissue elastic graph is generated; speeds or directions of shear wave sources are regulated, different shapes of shear wave planes are formed, a shear wave plane with a different shape generates a different sub elastic graph, and the sub elastic graphs are weighted and composited to form a final elastic graph. The tissue can be induced to generate large displacement, the imaging area is expanded, the shear wave elastic imaging robustness is enhanced, different shapes of shear waves are formed, the sub elastic graphs formed by the different shapes of shear waves are weighted, the imaging noise can be reduced, and the shear wave elastic imaging quality can be significantly improved.
Description
Technical field
The present invention relates to Medical ultrasound image technology field is and in particular to a kind of Multiple Shape plane based on Mach coning effect
Shearing wave composite imaging method.
Background technology
The elasticity (hardness) of tissue is closely related with its pathological state, and therefore elastogram can provide weight for medical diagnosis on disease
Will foundation.When occurring as diseases such as breast carcinoma, hepatocarcinoma, carcinoma of prostate, thyroid nodules, tissue elasticity would generally change.
Different with response detection method according to tissue excitation mode, in succession propose at present including quasistatic compression elasticity
Imaging harmony radiation force pulses elastogram etc. is in interior multiple elastograph imaging methods.Wherein, quasistatic compression elastogram makes
With probe static pressure tissue, organize relative resilient figure by detecting that the deformation organized produces, but quasistatic compression elastogram
Artifact noise substantially and is appropriate only for detecting superficial tissue；Acoustic radiation force pulse elastogram utilizes acoustic radiation force excitation group
Knit, then knit change in displacement situation using the burnt district's groups of detection pulse monitoring, because change in displacement is related to tissue elasticity, therefore may be used
Produce tissue elasticity figure to analyze by change in displacement.
Shearing wave elastogram is a kind of new imaging mode proposing in recent years, and it uses acoustic radiation force excitation tissue to produce
Raw shearing wave, then monitors shearing wave communication process using detection beam, by estimating shear wave velocity to echo analysis.?
In tissue linear elasticity simplified model, there is following relation with tissue modulus of shearing in shear wave velocity:
Here, c is shear wave velocity, and μ is tissue modulus of shearing, and ρ is tissue density's (constant).Therefore, can according to formula (1)
With by shear wave velocity direct inversion organize modulus of shearing μ it may be assumed that
μ=ρ c2Formula (2)
The present inventor finds although existing shearing wave elastogram overcomes quasistatic compression elasticity through research
Imaging is only limitted to superficial organ and acoustic radiation force pulse elastogram only detects that burnt district's groups knit the shortcoming of relative resilient, and is provided that
The quantitative elastic parameter value of human internal organs, but yet suffer from following deficiency: and 1, for security consideration, acoustic radiation force encourages
It is impossible to produce significantly shearing wave limited by intensity, the displacement of tissue of therefore induction is less, and detection difficulty is big；2nd, shearing wave exists
Decay in communication process quickly, imaging region is limited, out of focus area area image noise farther out is big；3rd, single shearing wave source is leaned on to produce
Raw shearing wave comes the elasticity of image tissue, poor robustness, is difficult to produce reliable high-quality elastic graph.
Content of the invention
The technical problem existing for art methods, the present invention provides a kind of Multiple Shape based on Mach coning effect to put down
Face shearing wave composite imaging method, the method is passed through to regulate and control acoustic radiation force in the longitudinal dynamic sequential focusing of different depth of tissue, produces
The shearing wave source that life is quickly moved, using Mach coning effect, the shearing wave that multiple wave sources produce is formed significantly due to coherent interference
Spending the shearing wave plane propagated to both sides, thus expanding detection zone, reducing picture noise；Change shearing wave source translational speed or
Direction, forms shearing wave plane of different shapes, is individually imaged the different son of i.e. generation to the shearing wave plane of each shape
Each elastic image finally is weighted being compounded to form final high-quality quantitation elastic graph, improves and cut by elastic image
Cut the robustness of ripple elastograph imaging method.
In order to solve above-mentioned technical problem, present invention employs following technical scheme:
A kind of Multiple Shape plane shear ripple composite imaging method based on Mach coning effect, the method comprises the following steps:
S1, the collection of reference frame echo data: collection tissue echo radiofrequency signal s0As reference frame, use it for calculating and cut
Cut the displacement of tissue of ripple induction；
S2, tissue excitation and plane shear ripple produce: regulation and control probe produces the acoustic radiation force promoting tissue so as to focus exists
The longitudinal different depth of tissue quickly focuses on, and forming speed is v1Mobile shearing wave source, due to Mach coning effect, many wave sources inductions
Shearing wave produce coherent interference, formed to both sides propagate shearing wave plane 1；Wherein, translational speed v1More than shear wave velocity
Degree c；
S3, shearing wave propagate Monitoring Data collection: in shearing wave communication process, to gather tissue more than 1000hz frame per second
Echo radiofrequency signal s1, s2, s3..., sm；
S4, displacement of tissue estimate with when m- Displacement Sequence image generate: by reference frame s0With Monitoring Data frame s1, s2,
s3..., smCarry out cross-correlation respectively, for estimating the tissue axial displacement of shearing wave propagation induction, using two-dimensional cross correlation displacement
Algorithm for estimating, according to following formula iteration generation axial displacement estimated value:
Wherein,Represent that in the i-th width displacement diagram picture, m row line n tracking window carries out the Displacement Estimation after l iteration
Value；X is lateral window width, and y is that axial window is long；δ x and δ y is laterally and axially to be spaced between adjacent estimating window respectively；M, n represent
The position subscript of estimating window；s0And siRepresent reference frame echo and the baseband signal of monitoring frame echo respectively；ω0Represent central angle
Frequency, * represents conjugation, and j is complex signal, and arg is tangent operation of negating；And when i takes 1,2,3 ..., during m, permissible using formula (3)
Produce m not displacement diagram picture in the same time；
S5, shear wave velocity estimate: according to when m- Displacement Sequence image, using time Peak Intensity Method estimate shearing wave propagate
The shearing wave time of advent of each lateral test position of direction, then estimate each detection of detection zone using time-of-flight method
The shear wave velocity c of position；
S6, sub- shearing wave elastic image g (v1) generate: according to known formula μ=ρ c2Computation organization's modulus of shearing μ, by group
Knit modulus of shearing value and be mapped to the gray level image shearing wave elastic image g (v of generation1)；
S7, regulation and control probe change shearing wave source translational speed viSize or direction, formed n-1 kind Mach of different shapes
Cone, thus produces n-1 kind shearing wave of different shapes plane 2, shearing wave plane 3, shearing wave plane 4 ..., shearing wave plane n；
Wherein, viIn i=2,3 ..., n, and n >=2；
S8, to each shearing wave plane of different shapes, reuse step s3, s4, s5 and s6 produces sub- shearing wave
Elastic image g (v2), g (v3), g (v4) ..., g (vn)；
S9, the final shearing wave composite elastic image of generation: by n width shearing wave elastic image g (v1), g (v2), g (v3), g
(v4) ..., g (vn) it is weighted being combined according to following formula, the final elastic image e of generation:
Wherein, wiIt is weights, n is compound subimage number, g (vi) it is to shear wave source translational speed for viWhen the son that produces
Shearing wave image, e (n) is the image of n sub- shearing wave elastic image composition generation.
Compared with prior art, what the present invention provided is combined into image space based on the Multiple Shape plane shear ripple of Mach coning effect
Method, by regulating and controlling acoustic radiation force incentive probe in the longitudinal different depth sequential focusing of tissue, forms the shearing wave source quickly moving,
The shearing wave coherent interference being produced due to Mach coning effect, many wave sources, forms the significantly shearing wave plane to two side shiftings；Profit
With ultrahigh speed ultrasonic monitoring shearing wave propagation condition, estimate shear wave velocity, modulus of shearing is organized in inverting, produce the elasticity of tissue
Figure；Regulation and control shearing wave source speed or direction, form shearing wave plane of different shapes, produce for each shape shearing wave plane
The weighting of elastic image is finally compounded to form final elastic graph by different elastic images.Therefore, the present invention can strengthen shearing wave
The robustness of elastogram, produces high-quality shearing wave elastic image；Using Mach coning effect, many wave sources shearing wave is made to be concerned with
Interfere and produce significantly shearing wave plane, expand imaging region, decrease picture noise；Elastic image weighting is compound can to enter
One step reduces picture noise, enhances imaging robustness, improves image quality.
Further, in described step s2, regulation and control are popped one's head in and are produced the acoustic radiation force promoting tissue so as to focus is longitudinal in tissue
15 millimeters, 20 millimeters, 25 millimeters and 30 millimeters of 4 different depths quickly focus on.
Further, in described step s3, in shearing wave communication process, with 3000hz frame per second collection tissue echo radio frequency letter
Number s1, s2, s3..., sm.
Further, in described step s6, described tissue modulus of shearing value is mapped to gray level image and includes: by minimum shearing
Modulus value is mapped to gray scale 0, and maximum shear modulus value is mapped to gray scale 255, and other modulus of shearing values are linearly reflected using following formula
It is mapped between 0-255:
I=255 (μ-μmin)/(μmax-μmin) formula (5)
Wherein, μmaxIt is maximum shear modulus value, μminIt is minimum modulus of shearing value, μ is mapped modulus of shearing value, i
It is the gray value after mapping.
Further, in described step s9, n width shearing wave elastic image is carried out averagely, carrying out averagely adding according to following formula
Power is compound, generation average weighted elastic image e:
Wherein, n is compound subimage number, g (vi) it is to shear wave source translational speed for viWhen the sub- shearing wave figure that produces
Picture, e (n) is the image of n sub- shearing wave elastic image composition generation.
Brief description
Fig. 1 is that the Multiple Shape plane shear ripple composite imaging method flow process based on Mach coning effect that the present invention provides is illustrated
Figure.
Fig. 2 is that positive 3 Mach of plane shear ripples provided in an embodiment of the present invention produce schematic diagram.
Fig. 3 is that positive 5 Mach of plane shear ripples provided in an embodiment of the present invention produce schematic diagram.
Fig. 4 is that reverse 3 Mach of plane shear ripples provided in an embodiment of the present invention produce schematic diagram.
Specific embodiment
In order that technological means, creation characteristic, reached purpose and effect that the present invention realizes are easy to understand, tie below
Conjunction is specifically illustrating, and the present invention is expanded on further.
In describing the invention it is to be understood that term " longitudinal ", " radially ", " length ", " width ", " thickness ",
" on ", D score, "front", "rear", "left", "right", " vertical ", " level ", " top ", " bottom ", " interior ", the orientation of instruction such as " outward " or
Position relationship is based on orientation shown in the drawings or position relationship, is for only for ease of the description present invention and simplifies description, and not
It is instruction or the hint device of indication or element must have specific orientation, with specific azimuth configuration and operation, therefore not
It is understood that as limitation of the present invention.In describing the invention, unless otherwise stated, " multiple " are meant that two or two
More than.
In describing the invention, it should be noted that unless otherwise clearly defined and limited, term " installation ", " phase
Even ", " connection " should be interpreted broadly, for example, it may be being fixedly connected or being detachably connected, or is integrally connected；Can
To be to be mechanically connected or electrical connection；Can be to be joined directly together it is also possible to be indirectly connected to by intermediary, Ke Yishi
The connection of two element internals.For the ordinary skill in the art, above-mentioned term can be understood at this with concrete condition
Concrete meaning in invention.
Refer to shown in Fig. 1-Fig. 4, the present invention is with positive 3 Mach, positive 5 Mach and reverse 3 Mach of three kinds of plane shears
As a example ripple, shearing wave composite elastic imaging method is described.Assume to be organized as even resilient tissue, shear wave velocity is 2 meter per seconds,
Then shearing wave source translational speed is respectively 6 meter per seconds, 10 meter per seconds and -6 meter per seconds, n=3.Specifically, the present invention provide based on
The Multiple Shape plane shear ripple composite imaging method of Mach coning effect comprises the following steps:
S1, the collection of reference frame echo data: collection tissue echo radiofrequency signal s0As reference frame, use it for calculating and cut
Cut the displacement of tissue of ripple induction；
S2, tissue excitation and plane shear ripple produce: regulation and control probe produces the acoustic radiation force promoting tissue so as to focus exists
Longitudinally 15 millimeters, 20 millimeters, 25 millimeters and 30 millimeters of 4 different depths quickly focus on, such as with shearing wave tissue from shallow to deep
The 2-10 speed of speed c focuses on, and forms the mobile shearing wave source that speed is 6 meter per seconds, specifically please join shown in Fig. 2, due to Mach
Coning effect, the shearing wave of 4 wave source inductions produces coherent interference, forms the shearing wave plane 1 propagated to both sides；Certainly, this area
Technical staff on the basis of the aforementioned embodiment focusing on number, can also be according to practical situation to longitudinally different deep in tissue
The quick number focusing on of degree is increased and decreased or is changed, and is actually needed with disclosure satisfy that；
S3, shearing wave propagate Monitoring Data collection: in shearing wave communication process, to gather tissue more than 1000hz frame per second
Echo radiofrequency signal s1, s2, s3..., sm, for example, tissue echo radiofrequency signal s is gathered with 3000hz frame per second1, s2, s3..., sm；
S4, displacement of tissue estimate with when m- Displacement Sequence image generate: by reference frame s0With Monitoring Data frame s1, s2,
s3..., smCarry out cross-correlation respectively, for estimating the tissue axial displacement of shearing wave propagation induction, using two-dimensional cross correlation displacement
Algorithm for estimating, according to following formula iteration generation axial displacement estimated value:
Wherein,Represent that in the i-th width displacement diagram picture, m row line n tracking window carries out the Displacement Estimation after l iteration
Value；X is lateral window width, and y is that axial window is long；δ x and δ y is laterally and axially to be spaced between adjacent estimating window respectively；M, n represent
The position subscript of estimating window；s0And siRepresent reference frame echo and the baseband signal of monitoring frame echo respectively；ω0Represent central angle
Frequency, * represents conjugation, and j is complex signal, and arg is tangent operation of negating；And when i takes 1,2,3 ..., during m, permissible using formula (3)
Produce m not displacement diagram picture in the same time；Certainly, those skilled in the art, on the basis of present embodiment, can also adopt
Calculated with the Displacement Estimation algorithm known to other skilled in the art；
S5, shear wave velocity estimate: according to when m- Displacement Sequence image, using time Peak Intensity Method (time-to-peak,
Ttp) estimate the shearing wave time of advent of shearing each lateral test position of direction of wave travel, then using time-of-flight method
(time-of-flight, tof) estimates the shear wave velocity c of each test position of detection zone；Wherein, using time Peak Intensity Method
Estimate the shearing wave time of advent and estimate the concrete grammar of shear wave velocity for those skilled in the art using time-of-flight method
Known, will not be described here；
S6, sub- shearing wave elastic image g (v1) generate: according to known formula μ=ρ c2I.e. formula (2) computation organization modulus of shearing
Then modulus of shearing value is mapped to gray level image and produces sub- shearing wave elastic image g (v by μ1)；As a kind of specific embodiment party
Minimum modulus of shearing value can be mapped to gray scale 0, maximum shear modulus value is mapped to gray scale 255 by formula, other modulus of shearing values
Between following formula Linear Mapping to 0-255:
I=255 (μ-μmin)/(μmax-μmin) formula (5)
Wherein, μmaxIt is maximum shear modulus value, μminIt is minimum modulus of shearing value, μ is mapped modulus of shearing value, i
It is the gray value after mapping；Certainly, those skilled in the art, on the basis of present embodiment, can also adopt other abilities
Method known to field technique personnel is mapped；
S7, again regulation and control probe are so as to focus is organizing longitudinal 15 millimeters, 20 millimeters, 25 millimeters and 30 millimeters 4 not
Quickly focus on from shallow to deep with depth, for example, focused on the 2-10 speed of shear wave velocity c, forming speed is 10 meter per seconds
Mobile shearing wave source, due to Mach coning effect, produces another shearing wave plane 2 of different shapes, specifically please join shown in Fig. 3；
For shearing wave plane 2, reuse step s3, s4, s5 and s6 produces sub- shearing wave elastic image g (v2)；
S8, again regulation and control probe, change focus direction so as to focus is organizing longitudinally 30 millimeters, 25 millimeters, 20 millimeters
Quickly focus on from deep to shallow with 15 millimeters of 4 different depths, for example, focused on the 2-10 speed of shear wave velocity c, form speed
Spend the mobile shearing wave source for -6 meter per seconds, due to Mach coning effect, produce another shearing wave plane 3 of different shapes, specifically
Please join shown in Fig. 4；For shearing wave plane 3, reuse step s3, s4, s5 and s6 produces sub- shearing wave elastic image g
(v3)；
S9, the final shearing wave composite elastic image of generation: by 3 width shearing wave elastic image g (v1), g (v2), g (v3) root
It is weighted being combined according to following formula, final elastic graph e of generation:
Wherein, wiIt is weights, g (vi) it is to shear wave source translational speed for viWhen the sub- shearing wave image that produces, e (3) is 3
The final elastic image of individual sub- shearing wave elastic image composition generation；As a kind of embodiment, weights take w here1=1, w2=
2, w3=1.
As the presently preferred embodiments, in described step s9, described weight w1、w2And w3Value equal, will 3 width shearing
Ripple elastic image carries out averagely, thus can carry out average weighted according to following formula and be combined, generation average weighted elastic image e:
Wherein, g (vi) it is to shear wave source translational speed for viWhen the sub- shearing wave image that produces, e (3) is 3 sub- shearing waves
The image of elastic image composition generation.
Certainly, those skilled in the art, on the basis of aforementioned enforcement, with other multiple plane shear ripples can be
Example, to realize the Multiple Shape plane shear ripple composite imaging method based on Mach coning effect of present invention offer, it is specifically cut
Cut ripple composite elastic imaging method similar with aforementioned three kinds of plane shear ripples, will not be described here.
Compared with prior art, what the present invention provided is combined into image space based on the Multiple Shape plane shear ripple of Mach coning effect
Method, by regulating and controlling acoustic radiation force incentive probe in the longitudinal different depth sequential focusing of tissue, forms the shearing wave source quickly moving,
The shearing wave coherent interference being produced due to Mach coning effect, many wave sources, forms the significantly shearing wave plane to two side shiftings；Profit
With ultrahigh speed ultrasonic monitoring shearing wave propagation condition, estimate shear wave velocity, modulus of shearing is organized in inverting, produce the elasticity of tissue
Figure；Regulation and control shearing wave source speed or direction, form shearing wave plane of different shapes, produce for each shape shearing wave plane
The weighting of elastic image is finally compounded to form final elastic graph by different elastic images.Therefore, the present invention can strengthen shearing wave
The robustness of elastogram, produces high-quality shearing wave elastic image；Using Mach coning effect, many wave sources shearing wave is made to be concerned with
Interfere and produce significantly shearing wave plane, expand imaging region, decrease picture noise；Elastic image weighting is compound can to enter
One step reduces picture noise, enhances imaging robustness, improves image quality.
Finally illustrate, above example only in order to technical scheme to be described and unrestricted, although with reference to relatively
Good embodiment has been described in detail to the present invention, it will be understood by those within the art that, can be to the skill of the present invention
Art scheme is modified or equivalent, the objective without deviating from technical solution of the present invention and scope, and it all should be covered at this
In the middle of the right of invention.
Claims (5)
1. a kind of Multiple Shape plane shear ripple composite imaging method based on Mach coning effect is it is characterised in that the method includes
Following steps:
S1, the collection of reference frame echo data: collection tissue echo radiofrequency signal s0As reference frame, use it for calculating shearing wave
The displacement of tissue of induction；
S2, tissue excitation and plane shear ripple produce: regulation and control probe produces the acoustic radiation force promoting tissue so as to focus is in tissue
Longitudinal different depth quickly focuses on, and forming speed is v1Mobile shearing wave source, due to Mach coning effect, the cutting of many wave sources inductions
Cut ripple and produce coherent interference, form the shearing wave plane 1 propagated to both sides；Wherein, translational speed v1More than shear wave velocity c；
S3, shearing wave propagate Monitoring Data collection: in shearing wave communication process, to gather tissue echo more than 1000hz frame per second
Radiofrequency signal s1, s2, s3..., sm；
S4, displacement of tissue estimate with when m- Displacement Sequence image generate: by reference frame s0With Monitoring Data frame s1, s2, s3..., sm
Carry out cross-correlation respectively, for estimating the tissue axial displacement of shearing wave propagation induction, calculated using two-dimensional cross correlation Displacement Estimation
Method, according to following formula iteration generation axial displacement estimated value:
Wherein,Represent that in the i-th width displacement diagram picture, m row line n follows the trail of the Displacement Estimation value after window carries out l iteration；X is
Lateral window width, y is that axial window is long；δ x and δ y is laterally and axially to be spaced between adjacent estimating window respectively；M, n represent estimating window
Position subscript；s0And siRepresent reference frame echo and the baseband signal of monitoring frame echo respectively；ω0Represent center angular frequency, *
Represent conjugation, j is complex signal, arg is tangent operation of negating；And when i takes 1,2,3 ..., during m, m can be produced using formula (3)
Not displacement diagram picture in the same time；
S5, shear wave velocity estimate: according to when m- Displacement Sequence image, using time Peak Intensity Method estimate shearing direction of wave travel
The shearing wave time of advent of each lateral test position, then estimate each test position of detection zone using time-of-flight method
Shear wave velocity c；
S6, sub- shearing wave elastic image g (v1) generate: according to known formula μ=ρ c2Computation organization's modulus of shearing μ, by tissue shear
Shear modulu value is mapped to gray level image and produces sub- shearing wave elastic image g (v1)；
S7, regulation and control probe change shearing wave source translational speed viSize or direction, formed n-1 kind Mach cone of different shapes, by
This produces n-1 kind shearing wave of different shapes plane 2, shearing wave plane 3, shearing wave plane 4 ..., shearing wave plane n；Wherein,
viIn i=2,3 ..., n, and n >=2；
S8, to each shearing wave plane of different shapes, reuse step s3, s4, s5 and s6 produces sub- shearing wave elasticity
Image g (v2), g (v3), g (v4) ..., g (vn)；
S9, the final shearing wave composite elastic image of generation: by n width shearing wave elastic image g (v1), g (v2), g (v3), g
(v4) ..., g (vn) it is weighted being combined according to following formula, the final elastic image e of generation:
Wherein, wiIt is weights, n is compound subimage number, g (vi) it is to shear wave source translational speed for viWhen produce son shearing
Ripple image, e (n) is the image of n sub- shearing wave elastic image composition generation.
2. the Multiple Shape plane shear ripple composite imaging method based on Mach coning effect according to claim 1, its feature
Be, in described step s2, regulation and control probe produce the acoustic radiation force promoting tissue so as to focus longitudinally 15 millimeters of tissue, 20
Millimeter, 25 millimeters and 30 millimeters of 4 different depths quickly focus on.
3. the Multiple Shape plane shear ripple composite imaging method based on Mach coning effect according to claim 1, its feature
It is, in described step s3, in shearing wave communication process, tissue echo radiofrequency signal s is gathered with 3000hz frame per second1, s2,
s3..., sm.
4. the Multiple Shape plane shear ripple composite imaging method based on Mach coning effect according to claim 1, its feature
It is, in described step s6, described tissue modulus of shearing value is mapped to gray level image and includes: minimum modulus of shearing value is mapped
Become gray scale 0, maximum shear modulus value is mapped to gray scale 255, other modulus of shearing values use following formula Linear Mapping to 0-255 it
Between:
I=255 (μ-μmin)/(μmax-μmin) formula (5)
Wherein, μmaxIt is maximum shear modulus value, μminIt is minimum modulus of shearing value, μ is mapped modulus of shearing value, and i is to reflect
Gray value after penetrating.
5. the Multiple Shape plane shear ripple composite imaging method based on Mach coning effect according to claim 1, its feature
It is, in described step s9, n width shearing wave elastic image is carried out averagely, carrying out average weighted according to following formula and being combined, produce
Average weighted elastic image e:
Wherein, n is compound subimage number, g (vi) it is to shear wave source translational speed for viWhen the sub- shearing wave image that produces, e
N () is the image of n sub- shearing wave elastic image composition generation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610893496.9A CN106373166B (en) | 2016-10-13 | 2016-10-13 | A kind of Multiple Shape plane shear wave composite imaging method based on Mach coning effect |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610893496.9A CN106373166B (en) | 2016-10-13 | 2016-10-13 | A kind of Multiple Shape plane shear wave composite imaging method based on Mach coning effect |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106373166A true CN106373166A (en) | 2017-02-01 |
| CN106373166B CN106373166B (en) | 2019-06-14 |
Family
ID=57895093
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610893496.9A Expired - Fee Related CN106373166B (en) | 2016-10-13 | 2016-10-13 | A kind of Multiple Shape plane shear wave composite imaging method based on Mach coning effect |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106373166B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107518918A (en) * | 2017-10-13 | 2017-12-29 | 无锡祥生医疗科技股份有限公司 | Ultrasonic elastograph imaging method and its system |
| CN107970043A (en) * | 2017-12-28 | 2018-05-01 | 深圳开立生物医疗科技股份有限公司 | A kind of detection method and device of shearing wave |
| CN109512465A (en) * | 2018-12-19 | 2019-03-26 | 深圳中科乐普医疗技术有限公司 | A kind of acoustic radiation force double direction shear wave composite imaging method and device |
| CN116269506A (en) * | 2023-05-17 | 2023-06-23 | 汕头市超声仪器研究所股份有限公司 | Shear wave and quasi-static composite elastic imaging method |
-
2016
- 2016-10-13 CN CN201610893496.9A patent/CN106373166B/en not_active Expired - Fee Related
Non-Patent Citations (3)
| Title |
|---|
| ANDREAS PESAVENTO 等: "《A Time-Efficient and Accurate Strain Estimation Concept for Ultrasonic Elastography Using Iterative Phase Zero Estimation》", 《IEEE TRANSACTIONS ON ULTRASCONICS，FERROELECTRICS，AND FREQUENCY CONTROL》 * |
| NED C ROUZE 等: "《Robust Estimation of Time-of-Flight Shear Wave》", 《2010 IEEE INTERNATIONAL ULTRASONICS SYMPOSIUM PROCEEDINGS》 * |
| 崔少国 等: "《弹性成像中不同频率子图像的相关性与频率复合》", 《中国生物医学工程学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107518918A (en) * | 2017-10-13 | 2017-12-29 | 无锡祥生医疗科技股份有限公司 | Ultrasonic elastograph imaging method and its system |
| CN107518918B (en) * | 2017-10-13 | 2020-06-26 | 无锡祥生医疗科技股份有限公司 | Ultrasonic elastography method and system thereof |
| CN107970043A (en) * | 2017-12-28 | 2018-05-01 | 深圳开立生物医疗科技股份有限公司 | A kind of detection method and device of shearing wave |
| CN109512465A (en) * | 2018-12-19 | 2019-03-26 | 深圳中科乐普医疗技术有限公司 | A kind of acoustic radiation force double direction shear wave composite imaging method and device |
| CN109512465B (en) * | 2018-12-19 | 2021-04-06 | 深圳中科乐普医疗技术有限公司 | Acoustic radiation force bidirectional shear wave composite imaging method and device |
| CN116269506A (en) * | 2023-05-17 | 2023-06-23 | 汕头市超声仪器研究所股份有限公司 | Shear wave and quasi-static composite elastic imaging method |
| CN116269506B (en) * | 2023-05-17 | 2023-07-28 | 汕头市超声仪器研究所股份有限公司 | Shear wave and quasi-static composite elastic imaging method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106373166B (en) | 2019-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104622502B (en) | Diagnostic ultrasound equipment and method | |
| US8961418B2 (en) | Solving for shear wave information in medical ultrasound imaging | |
| Jensen | Medical ultrasound imaging | |
| CN106373166A (en) | Multi-shape plane shear wave composite imaging method based on Mach cone effects | |
| CN103110432B (en) | Adaptive image optimization in induced wave ultrasound imaging | |
| CN104042247B (en) | It is imaged using the ultrasonic ARFI displacements in auto-adaptive time stage | |
| US10338203B2 (en) | Classification preprocessing in medical ultrasound shear wave imaging | |
| CN103239258B (en) | Hyperacoustic coaxial shear wave is adopted to characterize | |
| US9370338B2 (en) | Ultrasound diagnosis device | |
| US11154277B2 (en) | Tissue viscoelastic estimation from shear velocity in ultrasound medical imaging | |
| CN106618638B (en) | A kind of quantitative shearing wave elastogram system | |
| CN103732152A (en) | Diagnostic ultrasound apparatus and image processing method | |
| US20160018364A1 (en) | Methods, systems and computer program products for estimating shear wave speed using statistical inference | |
| CN109730722A (en) | Based on the elastograph imaging method for focusing ultrasonic sound and vibration signal | |
| CN106618639B (en) | A kind of quantitative shearing wave elastograph imaging method | |
| CN103330576B (en) | Micro-elasticity imaging method based on tissue microbubble dynamics model | |
| CN109259801B (en) | Shear wave elastic imaging method and device | |
| CN110507360B (en) | Shear wave imaging method and system | |
| CN110243950A (en) | Acoustic radiation force impulse Excitation Method for Frequency based on amplitude modulation | |
| KR101117545B1 (en) | Ultrasound Diagnostic System and Method For Forming Elasticity Image Using Overlapped Plane Wave | |
| CN106214182A (en) | HIFU damaging shear ripple elastic characteristic method of estimation based on LK optical flow method | |
| Nillesen et al. | Cardiac motion estimation using ultrafast ultrasound imaging tested in a finite element model of cardiac mechanics | |
| CN113827278B (en) | Method and device for determining propagation speed of shear wave | |
| CN106232016A (en) | Diagnostic ultrasound equipment | |
| Nightingale et al. | Comparison of qualitative and quantitative acoustic radiation force based elasticity imaging methods |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190614 Termination date: 20201013 |