CN106572838B - Elasticity measurement detection method and system - Google Patents
Elasticity measurement detection method and system Download PDFInfo
- Publication number
- CN106572838B CN106572838B CN201480081082.5A CN201480081082A CN106572838B CN 106572838 B CN106572838 B CN 106572838B CN 201480081082 A CN201480081082 A CN 201480081082A CN 106572838 B CN106572838 B CN 106572838B
- Authority
- CN
- China
- Prior art keywords
- data
- target area
- time
- wave
- shearing wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52019—Details of transmitters
- G01S7/5202—Details of transmitters for pulse systems
- G01S7/52022—Details of transmitters for pulse systems using a sequence of pulses, at least one pulse manipulating the transmissivity or reflexivity of the medium
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
A kind of elasticity measurement detection method and a kind of elasticity measurement detection system, the described method comprises the following steps, generate shearing wave (S101) in target area；Emit tracking pulse to target area and obtains the echo data (S102) of the tracking pulse；Data are carried out according to the echo data to calculate to obtain survey calculation result (S103)；Show the survey calculation result (S104).The dynamics that this method does not need operator's control pressing probe can carry out calculating acquisition to the relevant parameter of reflection tissue elasticity, so that the repeated and stability of measurement is guaranteed.
Description
Technical field
The present invention relates to medical detection field more particularly to a kind of elasticity measurement detection methods and system.
Background technique
Ultrasonic elastograph imaging is used to reflect the elastic or soft or hard degree of tissue (Tissue), in the auxiliary of tissue cancer lesion
Detection, good pernicious differentiation, prognosis restoration evaluation etc. obtain more and more applications.Conventional Ultrasound elastograph imaging method mainly leads to
It crosses probe pressing tissue and generates certain deformation, then the parameter relevant to tissue elasticity such as dependent variable, strain rate is calculated
And be imaged, reflect the elastic difference in area-of-interest between different tissues indirectly.But due to strain parameter and elastic parameter
Between corresponding relationship be under pressure the influence of size, even identical tissue or elastic parameter having the same are not
With tissue, if the pressure applied is different, the strain of generation is also different, may cause the erroneous judgement to elastic difference.And operator
It needs artificially to control the pressure operated every time when use and all be consistent, and remain uniform and stable, this makes the repeatability of image
It is challenged with stability.
Summary of the invention
It provides one kind not limited by operation requirement, elasticity measurement detection method and system with stably measured result.
A kind of elasticity measurement detection method, includes the following steps,
Shearing wave is generated in target area；
Emit tracking pulse to the target area and receives the echo data of the tracking pulse；
Data are carried out according to the echo data to calculate to obtain survey calculation result；
Show the survey calculation result.
Further, the elasticity measurement detection method is further comprising the steps of before target area generates shearing wave,
Emit reference pulse to target area, receives the reference pulse in the echo information in the target area, obtain
Take the echo data of the reference pulse.
Further, when carrying out data calculating according to the echo data to obtain survey calculation result, pass through data
It calculates and obtains displacement data, include the following steps when obtaining displacement data,
Take the echo data of any moment of any position in target area as Basic wave data；
Basic wave data are divided into multistage Basic wave Nuclear Data, the Basic wave of the Basic wave Nuclear Data is set
Nuclear Data Center；
Obtain on the echo data at each moment the maximum position of cross correlation between the Basic wave Nuclear Data Center
It sets, taking the difference of the maximum position of the cross correlation and the center of the Basic wave Nuclear Data is the displacement number at each moment
According to.
Further, when take the echo data of any moment of any position in target area as Basic wave data
When, the Basic wave data are using the echo data of the tracking pulse of any moment or the echo of the reference pulse of any moment
Data.
Further, when carrying out data calculating according to the echo data to obtain survey calculation result, pass through data
The particle velocity for obtaining different moments is calculated, is included the following steps when obtaining particle velocity,
After the displacement data for obtaining the tissue of target area, the displacement data of the same position different moments is taken；
Gradient is sought along time orientation to obtain the particle velocity curve of different moments at the position.
Further, when carrying out data calculating according to the echo data to obtain survey calculation result, pass through data
The propagation distance for obtaining the shearing wave in certain time is calculated, includes following step when obtaining the propagation distance and being calculated
Suddenly,
The displacement data at any position corresponding each moment in target area is taken, displacement-time curve is formed, finds out this
At the time of curve upward peak corresponds to；
Obtain displacement-lateral position curve at the time of the peak value in target area corresponds in the certain time of front and back
Perhaps displacement-lateral position curve or particle speed of the particle velocity-lateral position curve to each moment in the period
Degree-lateral position curve carries out cross-correlation judgement, obtains the propagation distance of the shearing wave in the period.
Further, it when carrying out data calculating according to the echo data to obtain survey calculation result, obtains different
Young's modulus ratio between target area meets following formula when obtaining the Young's modulus ratio,
E1/E2≈(c1*Tlong)2/(c2*Tlong)2≈d1 2/d2 2
In formula, d1、d1Indicate different target area in a period of time TlongThe propagation distance of the interior shearing wave, c1
And c2Indicate the spread speed of the shearing wave of different target areas.
Further, it when showing the survey calculation result, using wavefront film figure, propagation distance distribution map, propagates
Time distribution map, propagation distance curve graph, propagation distance schematic diagram, elastic ratio schematic diagram, spread speed duplicate ratio schematic diagram,
At least one of propagation distance duplicate ratio schematic diagram.
Further, when showing the survey calculation result, using grayscale or coloud coding mode.
A kind of elasticity measurement detection system, including ultrasonic probe, control module, signal processing module, computing module and aobvious
Showing device, the ultrasonic probe are equipped with transceiver module, the transceiver module of the ultrasonic probe, signal processing module, computing module,
Display device is sequentially connected, and the control module is connected to the transceiver module,
The transceiver module for emitting driving pulse, tracking pulse and reference pulse, and receive the tracking pulse and
The echo data of the reference pulse；
The control module is for controlling the transceiver module transmitting driving pulse, tracking pulse and reference pulse；
The signal processing module is used to carry out Signal Pretreatment to echo data；
The computing module is for carrying out processing calculating to the signal that Beam synthesis exports；
The display device is for showing the elastic image that the computing module generates.
Further, the computing module includes:
Particle velocity computing unit, for calculating the particle velocity of destination organization；
Flight distance calculation unit, for calculating the propagation distance of shearing wave；
Young's modulus ratio calculation unit, for calculating the Young's modulus ratio of destination organization.
The dynamics that method of the invention does not need operator's control pressing probe can be to the related ginseng of reflection tissue elasticity
Number carries out calculating acquisition, so that the repeatability and stability of measurement are guaranteed.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 is the elasticity measurement detection method flow diagram that the first better embodiment provides；
Fig. 2 to Fig. 4 is that pushing away for shearing wave propagation is generated in the elasticity measurement detection method of the first better embodiment offer
Dynamic wave emits schematic diagram；
Fig. 5 is that the displacement in the elasticity measurement detection method that the first better embodiment provides calculates schematic diagram；
Fig. 6 be propagation in elasticity measurement detection method that the first better embodiment provides in different target region away from
From-time graph schematic diagram；
Fig. 7 is the elasticity measurement detection method flow diagram that the second better embodiment provides；
Fig. 8 is the composition schematic diagram of elasticity measurement detection system provided by the invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Referring to Fig. 1, the first better embodiment of the invention provides a kind of elasticity measurement detection method, including following step
It is rapid:
Step S101 generates shearing wave in target area.In the present embodiment, target area can be according to the demand of elasticity measurement
It determines, to target area transmitting driving pulse (push pulse) to generate shearing wave behind target area.The target area can lead to
It crosses any suitable mode to select to determine, such as can be used, conventional two-dimensional B mode, which is imaged, common elastic imaging E mode is all kinds of is applicable in
It is determined after image checking mode Preliminary detection, it can also be selected according to detection demand.
The quantity of the target area can be one or more.When the quantity of target area is multiple, multiple mesh
The respective longitudinal depth or lateral position for marking region can be different.It is understood that being carried out when simultaneously to multiple target areas
When measurement, the average distance ratio between multiple target areas can be obtained by subsequent step, to reflect the multiple target
Elastic difference between region.
After determining target area, certain pulses, i.e. driving pulse can be emitted to close-proximity target zone by probe.It is described
Driving pulse transmitting time length ratio conventional Ultrasound transmitting pulse is long, and about tens arrive several hundred musec orders.The driving pulse can be with
Shearing wave source is generated in transmitting position, self cleavage wave source emits shearing wave and along different from the driving pulse direction of the launch
It propagates in direction.
As shown in Figures 2 to 4, in the present embodiment, emitting the driving pulse can be using focus emission or without focusing
Transmitting, the driving pulse can also be used transmitting and once or continuously emit repeatedly, can also adjust gathering for continuous several times transmitting
Depth of focus degree or/and lateral position, to form specific shearing direction of wave travel and widen the spread scope of shearing wave.
When confirming and multiple target areas being arranged, each target area can emit poly- with identical or different transmitting
Burnt mode, continuous emitting times, focal position multiple groups driving pulse, every group of driving pulse emit simultaneously for respective target area
Shearing wave is generated in the close-proximity target zone, so that the echo data to relevant shear wave carries out analytical calculation.Due to pushing arteries and veins
Punching transmitting needs biggish sound field energy to improve shearing intensity of wave, may need to be spaced between the transmitting of multiple groups driving pulse certain
Time to guarantee energy safety, the time interval of the transmitting of the intensity and each group driving pulse of the transmitting of every group of driving pulse by
System predefines.
Multiple target areas can also share same group of driving pulse transmitting, i.e. after one group of driving pulse transmitting, generation is cut
Wave propagation is cut, shearing wave, which is propagated, can pass through multiple target areas, then obtain the number of echoes of multiple target areas of above-mentioned process
According to progress analytical calculation.System can decide whether to share promotion according to target area depth and the distance of lateral position
Impulse ejection.
When emitting multiple groups pulse into target area, every group of transmitting may include repeatedly continuous transmitting, the arteries and veins emitted every time
Punching can have different focal positions, longitudinal depth, lateral position.It can be by between each position emitted of adjusting and/or transmitting
Every the time, to control the direction of propagation for the shearing wave that it is generated, propagate width.The shearing of multiple groups different directions or different location
Wave propagate after calculated result, may be coextruded, weight etc. and handle, to increase result accuracy.
Step S102 emits tracking pulse (tracking wave) to target area and receives its echo information, and obtains back
Wave number evidence.
In this step, the echo data refers to the echo data of the tracking pulse of tracking pulse.Emit to target area
Driving pulse and after generating shearing wave, transmitting tracking pulse simultaneously receives its echo information, thus when obtaining one section in target area
The echo data of tracking pulse in interior one section of spread scope.The emission duration of tracking pulse can predefine.It is described
The echo data of tracking pulse has recorded in shearing wave communication process, the tissue letter in above-mentioned one section of spread scope in each position
Breath.
Step S103 carries out data calculating according to the echo data, obtains survey calculation result.In the present embodiment,
It can be by tracking the echo data progress data calculating of pulse and obtaining displacement, particle velocity, each moment shearing in a period of time
The calculated results such as propagation time of shearing wave in the propagation distance of wave, a distance, and shown.
As shown in figure 5, the corresponding position of the tissue of target area can generate longitudinal position when shearing wave is by target area
It moves, i.e. change in location.The displacement calculates i.e. finger and is calculated the displacement of tissue of target area and obtain displacement data.Displacement
Calculating can further include steps of
Step S1031 takes the echo data of any moment of any position in target area as Basic wave data.
Step S1032, by Basic wave data sectional, every section is set as Basic wave Nuclear Data (kernel), described in setting
The Basic wave Nuclear Data Center of Basic wave Nuclear Data.It is understood that the length of every section of Basic wave Nuclear Data
Interval between the Basic wave Nuclear Data Center can be predefined by system.
Step S1033 obtains on the echo data at each moment the cross-correlation between the Basic wave Nuclear Data Center
Property maximum position, when to take the difference of the maximum position of the cross correlation and the center of the Basic wave Nuclear Data be each
The displacement data at quarter.
It is maximum with the Basic wave Nuclear Data Center cross correlation on the echo data at each moment in this step
Position and the difference of center of the Basic wave Nuclear Data be the corresponding displacement in the Basic wave Nuclear Data position.This
In step, can use the echo data of the tracking pulse of different moments on a certain lateral position in target area, respectively with the benchmark
Echo data carries out cross-correlation comparison, and the group on available each depth location is woven in the displacement data of different moments.The position
Shifting, which reflects, organizes interior relative longitudinal displacement of the particle relative to the above-mentioned propagation moment caused by shearing wave is propagated.
In step S1033, Block- matching is can be used in the echo data of the driving pulse of Basic wave data and different moments
(block-matching) mode carries out cross-correlation comparison, to obtain the shift value of change in location.It is understood that block
The matching criteria for algorithms selection and the Block- matching displacement matched can also sets itself.And since the displacement that shearing wave propagates generation is non-
It is often small, the segment data of same position on echo data can also be subjected to phase shift (phase shift) and calculated, directly pass through phase
Shifting calculates shift value, which is displacement data.
To one or more positions in target area, above-mentioned displacement calculating process can be carried out.Number is displaced to gained
According to certain filtering processing can be carried out to reduce noise in space or time dimension.
When carrying out data calculating according to the echo data in step S103, it may include particle velocity calculates.The particle
Speedometer the following steps are included:
Step S1034 after the displacement data for obtaining the tissue of target area, takes the displacement number of the same position different moments
According to seeking gradient along time orientation, the particle velocity curve of different moments at the position can be obtained.
The particle velocity curve reflects in shearing wave communication process movement wink of the particle in different moments on each position
Shi Sudu.To one or more positions in target area, above-mentioned particle velocity calculating process can be carried out.To gained particle
Speed data can carry out certain filtering processing in space or time dimension to reduce noise.
When carrying out data calculating according to the echo data in step S103, it may include flight distance calculation.The propagation
Distance calculate the following steps are included:
Step S1035 takes the displacement data at any position corresponding each moment in target area, forms displacement versus time
Curve, at the time of finding out the curve upward peak and correspond to.At the time of the peak value reflects shearing wave by the position at the time of correspondence.
It is understood that at the time of shearing wave is also that particle velocity-time graph peak value corresponds at the time of passing through the position, it can also
At the time of to find corresponding shearing wave using particle velocity-time graph of the position by the position.It, can in this step
To directly take out displacement or speed data near multiple particular moments in a period of time Δ T, above-mentioned flight distance calculation is carried out.
Difference is only that the distance calculated is difficult to which position that accurate judgement is located in target area, but can judge displacement shearing
After wave propagation how long.
Step S1036 obtains the position at the time of peak value in target area corresponds in the certain time Δ T of front and back
Shifting-lateral position curve or particle velocity-lateral position curve, to the displacement-transverse direction at each moment in the period Δ T
Position curve or particle velocity-lateral position curve carry out cross-correlation judgement, and the shearing wave obtained in the period is lateral
The difference of position, the difference reflect propagation distance of the shearing wave in above-mentioned a period of time Δ T.To the institute in target area
There is position to carry out above-mentioned calculation processing, in a period of time Δ when available shearing wave is by near each position in target area
Propagation distance in T.In this step, the length Δ T of the period can be determined voluntarily, which can be the above-mentioned moment
Before to later, before the above-mentioned moment or after the above-mentioned moment.
It is understood that can be directly taken out near multiple particular moments in a period of time Δ T in step S1035
Displacement or speed data, carry out above-mentioned flight distance calculation.It averages to above-mentioned propagation distance, in available target area
Average propagation distance in average propagation distance or shearing wave communication process in a period of time Δ T in a period of time Δ T.
Propagation distance in corresponding Δ T of multiple continuous moment is taken out, as shown in fig. 6, it is available propagation distance-when
Half interval contour can also be accumulated and obtain longer period of time TlongInterior total propagation distance.Two different positions can also be taken out
The displacement data at corresponding each moment is set, cross-correlation comparison is carried out, finds out the time difference of the two, then available shearing wave
Propagation time between these two positions.
When carrying out data calculating according to the echo data in step S103, it may include elastic ratio calculation, the elasticity
When ratio refers mainly to system and there are multiple target areas, the Young's modulus ratio of each target area can be calculated, it is each to reflect
Nonhomogeneous hardness degree between target area.
Under certain condition, the spread speed with tissue hardness of shearing wave have approximate fixed relationship:
E=3 ρ c2
In above formula, ρ indicates tissue density, and E indicates the Young's modulus value of tissue, and c indicates the spread speed of shearing wave.Centainly
Under the conditions of, Young's modulus is bigger, it is meant that tissue hardness is bigger.Assuming that the Young's modulus between two tissues is different, respectively E1
And E2, then its ratio meets:
E1/E2=3 ρ1c1 2/(3ρ2c2 2)
Assuming that tissue density is close, i.e. ρ1≈ρ2, then:
E1/E2≈c1 2/c2 2
A period of time T can be usedlongSquare ratio of interior shearing wave propagation distance calculates, it may be assumed that
E1/E2≈(c1*Tlong)2/(c2*Tlong)2≈d1 2/d2 2
In formula, d1、d1Different target areas is represented in a period of time TlongInterior shearing wave propagation distance, c1And c2It indicates
The spread speed of the shearing wave of different target areas.
It is of course also possible to directly calculate the spread speed in each target area, then the ratio between calculating speed square is come
Obtain elastic ratio.
When carrying out data calculating according to the echo data in step S103, it may include spread speed calculates.
To any depth in each target area, particle velocity-time graph or the displacement-of its each lateral position are taken out
Time graph relatively finds the time shift between each curve using cross-correlation, and it is mutual by above-mentioned participation which represents shearing wave
Close the time difference between two lateral positions compared.Lateral position difference and the ratio of time shift are the propagation for representing corresponding position
Speed.Also such as can use the propagation equation direct derivation of wave, to go out approximate formula as follows:
In formula, c indicates spread speed, uzLength travel data or longitudinal velocity data can be considered as, x represents laterally seat
Mark, z represent longitudinal coordinate.
Step S104 shows survey calculation result.After the calculating of above-mentioned elasticity measurement, result is shown there are many mode
Show.
Wavefront film figure is such as used, i.e., is taken out the displacement data at certain moment or particle velocity data, current time is formed
The distribution map of each position in target area, the distribution map can reflect the propagation position of shearing wave under current time.When multiple
When the displacement data or particle velocity data distribution map at quarter continuously play, shearing wave front propagation film figure can be formed, intuitively
Show the communication process of shearing wave.
Such as use propagation distance distribution map, i.e., by the propagation distance in the Δ T time at position each in target area into
Row imaging display.
Propagation time distribution map is such as used, i.e., by the propagation time between certain fixed range near each position in target area
Imaging is carried out to show.It can certainly the curve form display propagation time.
Propagation distance curve graph is such as used, i.e., by the propagation in target area since sometime in a period of time backward
Distance-time curve is drawn out, and whens multiple target areas can show a plurality of curve.
The position average propagation in target area is such as directly displayed apart from schematic diagram, or directly displays shearing wave and exists
Average propagation when each different moments in short time is apart from schematic diagram.
Such as when there are multiple target areas, show that the elastic ratio schematic diagram or spread speed in multiple target areas are put down
Side is than schematic diagram or propagation distance duplicate ratio schematic diagram.
During above-mentioned display, can be used grayscale or coloud coding, or with the image superposition of other modes, the side such as merge
Formula reinforces display effect.
As shown in fig. 7, the second preferred embodiment of the invention provides a kind of elasticity measurement detection method, the bullet of the present embodiment
Property measurement detection method it is roughly the same with the first preferred embodiment, comprising the following steps:
Step S201 emits reference pulse (reference wave) to target area, and receives the reference pulse in described
Echo information in target area obtains the echo data of reference pulse, and taking the echo data of reference pulse is benchmark number of echoes
According to.In the present embodiment, emit reference pulse before every group of driving pulse transmitting and receive its echo letter in target area
Breath, obtains the echo data of reference pulse and sets it as benchmark echo data.The echo information of the reference pulse has recorded
It generates shearing wave and starts the organizational information of the target area before propagating.
Step S202 generates shearing wave in target area.
Step S203 emits tracking pulse to target area and receives its echo information, and obtains echo data.In this reality
In this step for applying example, the echo data refers to the echo data of the tracking pulse of tracking pulse.
Step S204 carries out data calculating according to the echo data, obtains survey calculation result.In the present embodiment,
The echo data of reference pulse is taken as Basic wave data, by tracking the echo data of pulse and as Basic wave data
The echo data of reference pulse carry out data calculating and obtaining displacement, particle velocity, each moment shearing wave in a period of time
The calculated results such as propagation time of shearing wave in propagation distance, a distance.
Take in target area the echo data of the tracking pulse of different moments on a certain lateral position, respectively with the transverse direction position
The Basic wave data set carry out cross-correlation comparison, and the group on available each depth location is woven in the displacement of different moments
Data.The displacement mainly reflects particle in the tissue before organizing interior particle to propagate caused by shearing wave is propagated relative to shearing wave
Length travel.Particle velocity, in a period of time in the propagation distance, a distance of each moment shearing wave when the propagation of shearing wave
Between etc. the calculating process such as calculated results preferably implement in first roughly the same, details are not described herein.
Step S205 shows survey calculation result.
In the present embodiment, it is to emit reference pulse to target area in advance before target area generates shearing wave, and obtain
It takes the echo data of the reference pulse as Basic wave data, organizes interior particle opposite caused by shearing wave is propagated to obtain
The length travel of particle in tissue before shearing wave propagation, corresponding particle velocity, in a period of time each moment shearing wave biography
Broadcast distance, the calculated results such as calculated results such as propagation time of shearing wave in a distance.Other each steps in the present embodiment
Suddenly roughly the same with the first preferred embodiment, details are not described herein.
As shown in figure 8, the present invention also provides a kind of Ultrasonic elasticity measuring system, including ultrasonic probe 11, control module 12,
Signal processing module 13, computing module 15 and display device 17, the ultrasonic probe 11 are equipped with transceiver module 110, the ultrasound
The transceiver module 110 of probe 11, signal processing module 13, computing module 15, display device 17 are sequentially connected, the control module
12 are connected to the transceiver module 110.Wherein:
The transceiver module 110 receives the tracking pulse for emitting driving pulse, tracking pulse and reference pulse
With the echo data of the reference pulse.
The control module 12 emits driving pulse, tracking pulse and reference pulse for controlling the transceiver module 110.
In actual use, the control module 12 emits according to prefixed time interval to be made of driving pulse, reference pulse, tracking pulse
Specific ultrasonic sequence, consequently facilitating generating shearing wave in target area and for the transceiver module 110 of the ultrasonic probe 11
Receive corresponding echo data.
The signal processing module 13 is used to carry out Signal Pretreatment to echo data, consequently facilitating computing module 15 carries out
Subsequent calculating, the Signal Pretreatment may include Beam synthesis processing, may also include such as signal amplification, analog-to-digital conversion, orthogonal point
Solution etc..
Computing module 15 is for carrying out processing calculating to the signal that Beam synthesis exports.The computing module 15 can be used for counting
All kinds of propagation parameters for calculating shearing wave, specifically, the computing module 15 includes:
Particle velocity computing unit 151, for calculating the particle velocity of destination organization；
Flight distance calculation unit 152, for calculating the propagation distance of shearing wave；
Young's modulus ratio calculation unit 153, for calculating the Young's modulus ratio of destination organization.
Computing module 15 is started in the shearing wave propagation distance or certain time at each moment by calculating some position
Average distance or different location start the average distance ratio in certain time or the average propagation speed in several target areas
Ratio etc. is spent, to reflect the physical characteristic of target area, and generates corresponding elastic image.
Display device 17 is for showing the elastic image that the computing module 15 generates.
Elasticity measurement detection method of the invention and system, by emitting one into the target area of organization internal to be detected
Serial certain pulses generate the propagation of shearing wave, then constantly emit receives echo-signal within the scope of shearing propagation path,
The propagation distance curve of shearing wave in a period of time is extracted from the echo-signal, and calculates average distance, it is final to show.
If under identical launching condition, within the same time, shearing wave is propagated remoter, then illustrates that the hardness of tissue is higher, otherwise say
Bright tissue is softer.
The dynamics that method of the invention does not need operator's control pressing probe can be to the related ginseng of reflection tissue elasticity
Number carries out calculating acquisition, so that the repeatability and stability of measurement are guaranteed.On the other hand, the present invention can obtain quantitative
Parameter value, especially suitable for the elasticity measurement to uniform formation.Meanwhile the elastic difference between different tissues can also be provided in the present invention
Parameter, to the elastic difference degree reflected between different tissues.
Above disclosed is only a preferred embodiment of the present invention, cannot limit the power of the present invention with this certainly
Sharp range, those skilled in the art can understand all or part of the processes for realizing the above embodiment, and weighs according to the present invention
Benefit requires made equivalent variations, still belongs to the scope covered by the invention.
Claims (11)
1. a kind of elasticity measurement detection method, it is characterised in that: include the following steps,
Shearing wave is generated in target area；
Emit tracking pulse to the target area and receives the echo data of the tracking pulse；
Data are carried out according to the echo data to calculate to obtain survey calculation result；
Show the survey calculation result；
It is wherein, described that data calculating is carried out to obtain survey calculation result according to the echo data, comprising the following steps:
The displacement data that data calculate acquisition target area is carried out according to the echo data；
Take the displacement data at any position corresponding each moment in target area, so obtain shearing wave by the position when
It carves；
Obtain the shearing wave propagation distance at the time of shearing wave in target area passes through the position in a period of time of front and back.
2. elasticity measurement detection method as described in claim 1, it is characterised in that:
The elasticity measurement detection method is further comprising the steps of before target area generates shearing wave,
Emit reference pulse to target area, receives the reference pulse in the echo information in the target area, obtain institute
State the echo data of reference pulse.
3. elasticity measurement detection method as claimed in claim 2, it is characterised in that: described to be counted according to the echo data
According to calculate obtain target area displacement data the following steps are included:
Take the echo data of any moment of any position in target area as Basic wave data；
Basic wave data are divided into multistage Basic wave Nuclear Data, the Basic wave nucleus number of the Basic wave Nuclear Data is set
According to center；
The maximum position of cross correlation between the Basic wave Nuclear Data Center is obtained on the echo data at each moment, is taken
The difference of the maximum position of the cross correlation and the center of the Basic wave Nuclear Data is the displacement data at each moment.
4. elasticity measurement detection method as claimed in claim 3, it is characterised in that: as times for taking any position in target area
When the echo data at one moment is as Basic wave data, the Basic wave data track returning for pulse using any moment
The echo data of the reference pulse of wave number evidence or any moment.
5. elasticity measurement detection method according to any one of claims 1-4, it is characterised in that: described to obtain shearing wave process
At the time of the position, comprising the following steps:
The displacement data for taking any position corresponding each moment in target area seeks gradient along time orientation to obtain the position
At the time of locating the particle velocity data at each moment, and then obtain shearing wave by the position；Or
The displacement data at any position corresponding each moment in target area is taken, displacement-time curve is formed, finds out the curve
At the time of upward peak corresponds to, at the time of as shearing wave passes through the position.
6. elasticity measurement detection method as claimed in claim 5, it is characterised in that: obtaining the shearing wave propagation distance includes
Following steps,
Displacement-lateral position curve before and after the peak value in target area is obtained at the time of correspond in certain time or
Particle velocity-lateral position curve, to the displacement-lateral position curve or particle velocity-at each moment in the period
Lateral position curve carries out cross-correlation judgement, obtains the propagation distance of the shearing wave in the period.
7. elasticity measurement detection method as claimed in claim 6, it is characterised in that: when according to echo data progress data
When calculating to obtain survey calculation result, the Young's modulus ratio between different target region is obtained, when the acquisition Young mould
When measuring ratio, meet following formula,
E1/E2≈(c1*Tlong)2/(c2*Tlong)2≈d1 2/d2 2
In formula, d1、d1Indicate different target area in a period of time TlongThe propagation distance of the interior shearing wave, c1And c2Table
Show the spread speed of the shearing wave of different target areas.
8. elasticity measurement detection method as described in claim 1, it is characterised in that: when showing the survey calculation result,
Using wavefront film figure, propagation distance distribution map, propagation time distribution map, propagation distance curve graph, propagation distance schematic diagram, bullet
At least one of property ratio schematic diagram, spread speed duplicate ratio schematic diagram, propagation distance duplicate ratio schematic diagram.
9. elasticity measurement detection method as claimed in claim 8, it is characterised in that: when showing the survey calculation result,
Using grayscale or coloud coding mode.
10. a kind of elasticity measurement detection system, it is characterised in that: including ultrasonic probe, control module, signal processing module, meter
Module and display device are calculated, the ultrasonic probe is equipped with transceiver module, transceiver module, the signal processing mould of the ultrasonic probe
Block, computing module, display device are sequentially connected, and the control module is connected to the transceiver module,
The transceiver module for emitting driving pulse, tracking pulse and reference pulse, and receive the tracking pulse with it is described
The echo data of reference pulse；
The control module is for controlling the transceiver module transmitting driving pulse, tracking pulse and reference pulse；
The signal processing module is used to carry out Signal Pretreatment to echo data；
The computing module is for carrying out processing calculating to the signal that Beam synthesis exports；
The display device is for showing the elastic image that the computing module generates；
Wherein, the computing module is for carrying out processing calculating to the signal that Beam synthesis exports, comprising the following steps:
The displacement data that data calculate acquisition target area is carried out according to the signal of Beam synthesis output；
Take the displacement data at any position corresponding each moment in target area, so obtain shearing wave by the position when
It carves；
Obtain the shearing wave propagation distance at the time of shearing wave in target area passes through the position in a period of time of front and back.
11. elasticity measurement detection system as claimed in claim 10, it is characterised in that: the computing module includes:
Particle velocity computing unit, for calculating the particle velocity of destination organization；
Flight distance calculation unit, for calculating the propagation distance of shearing wave；
Young's modulus ratio calculation unit, for calculating the Young's modulus ratio of destination organization.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910741685.8A CN110368031B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
| CN201910741683.9A CN110432926B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2014/085821 WO2016033752A1 (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement and detection method and system |
Related Child Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910741685.8A Division CN110368031B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
| CN201910741683.9A Division CN110432926B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106572838A CN106572838A (en) | 2017-04-19 |
| CN106572838B true CN106572838B (en) | 2019-09-06 |
Family
ID=55439001
Family Applications (3)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910741683.9A Active CN110432926B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
| CN201910741685.8A Active CN110368031B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
| CN201480081082.5A Active CN106572838B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
Family Applications Before (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910741683.9A Active CN110432926B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
| CN201910741685.8A Active CN110368031B (en) | 2014-09-03 | 2014-09-03 | Elasticity measurement detection method and system |
Country Status (2)
| Country | Link |
|---|---|
| CN (3) | CN110432926B (en) |
| WO (1) | WO2016033752A1 (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6996035B2 (en) * | 2017-11-02 | 2022-01-17 | 富士フイルムヘルスケア株式会社 | Ultrasonic diagnostic device and method for evaluating physical properties of living tissue |
| CN108158610B (en) * | 2018-01-16 | 2024-02-06 | 苏州国科昂卓医疗科技有限公司 | Elastic imaging method, device, equipment and ultrasonic imaging probe |
| CN110494082B (en) * | 2018-04-11 | 2023-05-05 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasound elastography method and system |
| CN111401108B (en) * | 2019-01-02 | 2021-04-06 | 无锡海斯凯尔医学技术有限公司 | Slope parameter extraction method and device and computer readable storage medium |
| CN111399039B (en) * | 2019-01-02 | 2023-04-07 | 无锡海斯凯尔医学技术有限公司 | Slope parameter extraction method, device and storage medium |
| CN111399038B (en) * | 2019-01-02 | 2023-04-07 | 无锡海斯凯尔医学技术有限公司 | Slope parameter extraction method and device and computer readable storage medium |
| CN109893172B (en) * | 2019-02-22 | 2020-06-19 | 清华大学 | Method and device for determining mechanical parameters based on elastography and computer equipment |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101699280A (en) * | 2009-10-15 | 2010-04-28 | 北京索瑞特医学技术有限公司 | Method and device for ultrasonic and nondestructive detection ofelasticity of viscoelastic medium |
| CN102469980A (en) * | 2009-07-17 | 2012-05-23 | 皇家飞利浦电子股份有限公司 | Spatially-fine shear wave dispersion ultrasound vibrometry sampling |
| CN102525568A (en) * | 2012-01-17 | 2012-07-04 | 北京索瑞特医学技术有限公司 | Subtraction elastography method |
| CN102824194A (en) * | 2011-06-14 | 2012-12-19 | 深圳迈瑞生物医疗电子股份有限公司 | Displacement detecting method and device thereof in elasticity imaging |
| CN103269639A (en) * | 2010-12-22 | 2013-08-28 | 皇家飞利浦电子股份有限公司 | Shear wave velocity estimation using center of mass |
| CN103300890A (en) * | 2012-03-16 | 2013-09-18 | 通用电气公司 | System and method for measuring tissue mechanical characteristics |
| CN103908289A (en) * | 2012-12-31 | 2014-07-09 | Ge医疗系统环球技术有限公司 | Method for removing background noise in shear wave and corresponding ultrasonic imaging system |
Family Cites Families (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7374538B2 (en) * | 2000-04-05 | 2008-05-20 | Duke University | Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing |
| US6702744B2 (en) * | 2001-06-20 | 2004-03-09 | Advanced Cardiovascular Systems, Inc. | Agents that stimulate therapeutic angiogenesis and techniques and devices that enable their delivery |
| WO2006013916A1 (en) * | 2004-08-05 | 2006-02-09 | Hitachi Medical Corporation | Method for displaying elastic image and ultrasonograph |
| US8118744B2 (en) * | 2007-02-09 | 2012-02-21 | Duke University | Methods, systems and computer program products for ultrasound shear wave velocity estimation and shear modulus reconstruction |
| WO2010104863A1 (en) * | 2009-03-09 | 2010-09-16 | Mayo Foundation For Medical Education And Research | Method for ultrasound vibrometry using orthogonal basis functions |
| BRPI1010214B8 (en) * | 2009-06-30 | 2021-06-22 | Koninklijke Philips Electronics Nv | shear wave dispersion (sduv) vibrometry method, non-transient computer readable media, and ultrasound device configured for shear wave dispersion (sduv) vibrometry |
| JP5559788B2 (en) * | 2009-07-07 | 2014-07-23 | 株式会社日立メディコ | Ultrasonic diagnostic equipment |
| US8500639B2 (en) * | 2009-09-11 | 2013-08-06 | Mr Holdings (Hk) Limited | Systems and methods for shear wave field formation |
| CN102667522B (en) * | 2009-11-25 | 2014-10-08 | 皇家飞利浦电子股份有限公司 | Ultrasonic shear wave imaging with focused scanline beamforming |
| BR112012027443A2 (en) * | 2010-04-27 | 2016-07-19 | Halliburton Energy Services Inc | apparatus, system, processing method implemented, and article. |
| CN101869485B (en) * | 2010-06-23 | 2012-07-04 | 深圳大学 | Ultrasonic imaging method and device |
| WO2012065151A1 (en) * | 2010-11-14 | 2012-05-18 | University Of Utah Research Foundation | Ultrasonic method and system for determining tissue pathology |
| WO2012116364A1 (en) * | 2011-02-25 | 2012-08-30 | Mayo Foundation For Medical Education And Research | Ultrasound vibrometry with unfocused ultrasound |
| US20120253194A1 (en) * | 2011-03-30 | 2012-10-04 | Tadashi Tamura | Methods and apparatus for ultrasound imaging |
| US10357226B2 (en) * | 2012-03-12 | 2019-07-23 | Mayo Foundation For Medical Education And Research | System and method for model-independent quantification of tissue viscoelastic properties using ultrasound |
| US9220479B2 (en) * | 2012-03-30 | 2015-12-29 | Hitachi Aloka Medical, Ltd. | Methods and apparatus for ultrasound imaging |
| JP5733835B2 (en) * | 2012-08-27 | 2015-06-10 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Measuring device and ultrasonic diagnostic device |
| CN103845074B (en) * | 2012-11-28 | 2017-12-12 | 深圳迈瑞生物医疗电子股份有限公司 | A kind of ultrasonic elastograph imaging system and method |
| CN103462643B (en) * | 2013-09-29 | 2015-02-11 | 深圳市开立科技有限公司 | Shear wave speed measurement method, device and system |
| CN103720489B (en) * | 2013-12-30 | 2015-10-28 | 中国科学院深圳先进技术研究院 | Pathological tissues growth monitoring method and system |
-
2014
- 2014-09-03 CN CN201910741683.9A patent/CN110432926B/en active Active
- 2014-09-03 WO PCT/CN2014/085821 patent/WO2016033752A1/en active Application Filing
- 2014-09-03 CN CN201910741685.8A patent/CN110368031B/en active Active
- 2014-09-03 CN CN201480081082.5A patent/CN106572838B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102469980A (en) * | 2009-07-17 | 2012-05-23 | 皇家飞利浦电子股份有限公司 | Spatially-fine shear wave dispersion ultrasound vibrometry sampling |
| CN101699280A (en) * | 2009-10-15 | 2010-04-28 | 北京索瑞特医学技术有限公司 | Method and device for ultrasonic and nondestructive detection ofelasticity of viscoelastic medium |
| CN103269639A (en) * | 2010-12-22 | 2013-08-28 | 皇家飞利浦电子股份有限公司 | Shear wave velocity estimation using center of mass |
| CN102824194A (en) * | 2011-06-14 | 2012-12-19 | 深圳迈瑞生物医疗电子股份有限公司 | Displacement detecting method and device thereof in elasticity imaging |
| CN102525568A (en) * | 2012-01-17 | 2012-07-04 | 北京索瑞特医学技术有限公司 | Subtraction elastography method |
| CN103300890A (en) * | 2012-03-16 | 2013-09-18 | 通用电气公司 | System and method for measuring tissue mechanical characteristics |
| CN103908289A (en) * | 2012-12-31 | 2014-07-09 | Ge医疗系统环球技术有限公司 | Method for removing background noise in shear wave and corresponding ultrasonic imaging system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110368031B (en) | 2022-10-04 |
| CN106572838A (en) | 2017-04-19 |
| CN110432926A (en) | 2019-11-12 |
| CN110432926B (en) | 2022-06-07 |
| CN110368031A (en) | 2019-10-25 |
| WO2016033752A1 (en) | 2016-03-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106572838B (en) | Elasticity measurement detection method and system | |
| CN103988092B (en) | Position-measurement device | |
| CN104203112B (en) | Diagnostic ultrasound equipment and photoelastic evaluation method | |
| CN105816205B (en) | Sparse tracking in acoustic radiation force Pulse Imageing | |
| CN103462643B (en) | Shear wave speed measurement method, device and system | |
| CN101852773B (en) | Ultrasonic imaging apparatus with a virtual apex | |
| CN107613878A (en) | Ultrasonic image-forming system and method for detection object motion | |
| EP3153884A1 (en) | Detection apparatus, fish finder, and radar | |
| EP3190976B1 (en) | Ultrasound system and method for measurement using shear waves | |
| US20210356573A1 (en) | Two dimension and three dimension imaging based on speed changes of sound/ultrasound | |
| US20160363561A1 (en) | Methods and systems for measuring properties with ultrasound | |
| US20140276003A1 (en) | Systems and Methods to Detect and Present Interventional Devices via Ultrasound Imaging | |
| CN105877783B (en) | Two-dimensional shear wave elastic imaging method and device | |
| CN103239258A (en) | On-axis shear wave characterization with ultrasound | |
| CN101828929B (en) | Vector measurement method of Doppler blood flow velocity by utilizing apparent displacement | |
| CN106955125A (en) | Motion independence in acoustic radiation power Pulse Imageing | |
| EP2901939A1 (en) | Ultrasonic measurement apparatus and ultrasonic measurement method | |
| CN105911537B (en) | A method of reducing active sonar detection blind area | |
| RU2338220C1 (en) | Method of measurement of shell speed | |
| CN106443646B (en) | A kind of ultrasonic ranging system, echo processing techniques and device | |
| US20190350558A1 (en) | Shear wave imaging based on ultrasound with increased pulse repetition frequency | |
| CN111050661B (en) | Shear wave elasticity measurement method and shear wave elasticity imaging system | |
| RU167164U1 (en) | HYDROACOUSTIC MULTI-BEAM SYSTEM FOR DETERMINING NUMBER AND PARAMETERS OF FISH | |
| RU2798390C1 (en) | Method for passive determination of coordinates of a noise-generating object | |
| RU2496117C1 (en) | Method of measuring displacement parameters of probing signal source |
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 |