CN107510474A - Shearing wave elastograph imaging method and system - Google Patents
Shearing wave elastograph imaging method and system Download PDFInfo
- Publication number
- CN107510474A CN107510474A CN201710859560.6A CN201710859560A CN107510474A CN 107510474 A CN107510474 A CN 107510474A CN 201710859560 A CN201710859560 A CN 201710859560A CN 107510474 A CN107510474 A CN 107510474A
- Authority
- CN
- China
- Prior art keywords
- shearing wave
- shearing
- information
- amplitude
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention discloses the present invention by measuring velocity information and dampening information in shearing wave communication process, and shown, it is single to overcome shearing wave elastogram systematic survey object present in prior art, organizational information is caused to reflect incomplete technical problem, the shearing wave elastograph imaging method and system of a kind of speed for reflecting shearing wave and dampening information are realized, and then is advantageous to improve the reliability and comprehensive of shearing wave detection.It the composite can be widely applied to various shearing wave elastogram systems.
Description
Technical field
It the present invention relates to the use of ultrasound imaging, more particularly to a kind of shearing wave elastograph imaging method and be
System.
Background technology
Two dimension shearing ripple elastography is mainly by measuring the spread speed of shearing wave in the tissue, then directly
It is corresponding different velocity amplitudes is mapped to different colours to be shown on image by shear wave velocity with specific color table；Or
Shear wave velocity is scaled Young's modulus, then passes through specific face by person according to the other specification such as formula, conjunctive tissue density
Color table, corresponding different Young's modulus value are mapped as different colours and are shown on image.Doctor or operating personnel can be straight
The hardness information that region of interest tract tissue is judged according to color of image seen, and then determine abnormal area.
But shear wave velocity and the Young's modulus calculated according to shear wave velocity are shearing wave communication process
The middle part that can reflect organizational information, only reflect that the reliability of the hardness information shearing wave detection of tissue is low.
The content of the invention
In order to preferably develop purposes of the shearing wave in tissue detection field, it is an object of the invention to provide a kind of shearing wave
Elastograph imaging method and system, this method and system can reflect the speed and dampening information of shearing wave, be advantageous to improve shearing
The reliability and comprehensive of ripple detection.
The technical solution adopted in the present invention is：
A kind of shearing wave elastograph imaging method, including step：
S1, encouraged using acoustic radiation force in predeterminated position and produce shearing wave；
S2, repeatedly transmitting ultrasound detection wave beam detects shearing wave in shearing propagation path, to obtain shearing wave detection
Data；
S3, parsing shearing wave testing number obtain the velocity information and dampening information of shearing wave according to this；
S4, shear wave velocity information and dampening information, image corresponding to output are handled using default image mapping algorithm
Information.
Preferably, the step S3 is specially：Using shearing wave detection data generation shearing waveform data, according to adjacent
Shearing wave detection scan line shearing waveform data calculate shearing wave velocity information and dampening information, the dampening information
Including shearing wave amplitude fading information and/or shearing wave frequency decay information and/or shearing wave time domain waveform change information.
Preferably, the step S3 includes obtaining shear wave velocity information sub-step：Swept according to adjacent shearing wave detection
The shearing wave for retouching line position propagates arrival time, calculates the shear wave velocity information of the position.
Preferably, the step S3 includes obtaining shearing wave amplitude fading information sub-step：It is determined that shearing waveform data
Axial displacement data or axial velocity data peak value, set an amplitude cut-off region, obtain in amplitude cut-off region
Amplitude average value, more adjacent shearing wave detection scan line amplitude average value with obtain shearing wave amplitude fading letter
Breath.
Preferably, the step S3 includes obtaining shearing wave frequency decay information sub-step：Handled using Fourier transformation
Waveform data is sheared, determines the peak value of the axial displacement data after Fourier transformation or axial velocity data, a width is set
It is worth cut-off region, obtains the average value of the amplitude in amplitude cut-off region, more adjacent shearing wave detects the width of scan line
The average value of degree is to obtain shearing wave frequency decay information.
Preferably, the step S2 includes obtaining shearing wave time domain waveform change information sub-step：Determine shearing wave waveform
The axial displacement data of data or the peak value of axial velocity data, an amplitude cut-off region is set, obtained in amplitude cut-off region
Amplitude in domain square average value, the amplitude of more adjacent shearing wave detection scan line square average value to obtain
Shearing wave time domain waveform change information.
Preferably, the step S1 is specially：In the area-of-interest of shearing wave elastogram, utilized in predeterminated position
Acoustic radiation force excitation produces shearing wave, according to the propagation distance of shearing wave caused by the size and single activation of area-of-interest
To determine the number of ultrasonic acoustic radiation force excitation.
Preferably, the step S2 is specially：Repeatedly the detection of transmitting ultrasound detection wave beam is cut in shearing propagation path
Ripple is cut, to obtain echo-signal, using Beam synthesis algorithm process echo-signal to generate multiple shearing wave detection scan lines, and
The shearing wave at multiple moment, which is generated, in each shearing wave detection scan line detects data.
A kind of shearing wave elastogram system, including：Shearing wave encourages module, for utilizing acoustic radiation force in predeterminated position
Excitation produces shearing wave；Detection module, for launching ultrasound detection wave beam detection shearing wave in shearing propagation path, to obtain
Shearing wave is taken to detect data；Computing module, obtain velocity information and the decay of shearing wave according to this for parsing shearing wave testing number
Information；Image processing module, for utilizing default image mapping algorithm processing shear wave velocity information and dampening information, output
Corresponding image information.
Preferably, the dampening information include shearing wave amplitude fading information and/or shearing wave frequency decay information and/or
Shearing wave time domain waveform change information, the computing module specifically include：Shear wave velocity information calculating sub module, for basis
The shearing wave of adjacent shearing wave detection scan line position propagates arrival time, calculates the shear wave velocity information of the position；Cut
Wave amplitude dampening information submodule is cut, for determining the axial displacement data of shearing waveform data or the peak of axial velocity data
Value, an amplitude cut-off region is set, obtains the average value of the amplitude in amplitude cut-off region, the inspection of more adjacent shearing wave
The average value of the amplitude of scan line is surveyed to obtain shearing wave amplitude fading information；Shearing wave frequency decay message sub-module, is used for
The peak value of the axial displacement data after Fourier transformation or axial velocity data is determined, an amplitude cut-off region is set, is obtained
The average value of amplitude in amplitude cut-off region, the average value of the amplitude of more adjacent shearing wave detection scan line is to obtain
Shearing wave frequency decay information；Shearing wave time domain waveform change information submodule, for determining the axial direction of shearing waveform data
The peak value of displacement data or axial velocity data, an amplitude cut-off region is set, obtains the amplitude in amplitude cut-off region
Square average value, the amplitude of more adjacent shearing wave detection scan line square average value to obtain shearing wave time domain
Waveform change information.
The beneficial effects of the invention are as follows：
The present invention is shown by measuring velocity information and dampening information in shearing wave communication process, is overcome
Shearing wave elastogram systematic survey object present in prior art is single, causes organizational information to reflect that incomplete technology is asked
Topic, realize the shearing wave elastograph imaging method and system of a kind of speed for reflecting shearing wave and dampening information, Jin Eryou
Beneficial to the reliability and comprehensive for improving shearing wave detection.
In addition, calculating of the present invention for shear wave velocity information and shearing wave dampening information is all using by ripple
Shearing wave detection data after Shu Hecheng carry out auto-correlation, the axial velocity for the shearing wave that cross-correlation scheduling algorithm obtains or axial position
Data are moved, therefore can easily be implemented in systems.
It the composite can be widely applied to various shearing wave elastogram systems.
Brief description of the drawings
The embodiment of the present invention is described further below in conjunction with the accompanying drawings：
Fig. 1 is the schematic flow sheet of shearing wave elastograph imaging method one embodiment of the present invention；
Fig. 2 is the schematic flow sheet of shearing wave elastograph imaging method another kind embodiment of the present invention；
Fig. 3 is the principle schematic of shearing wave elastograph imaging method one embodiment of the present invention；
Fig. 4 is the time domain waveform schematic diagram of shearing wave elastograph imaging method one embodiment of the present invention；
Fig. 5 is the frequency-domain waveform schematic diagram of shearing wave elastograph imaging method one embodiment of the present invention.
Embodiment
It should be noted that in the case where not conflicting, the feature in embodiment and embodiment in the application can phase
Mutually combination.
The research of biological histology shows：Viscoplasticity is to reflect the optimal parameter of tissue mechanical information.That is,
In addition to elasticity, viscosity is also required to take into account, is only not enough to total evaluation tissue by single shear wave velocity information
Mechanical attribute.Meanwhile the most tissues in human body have anisotropy and dynamic, such as muscle, liver again.Therefore need
Other specification assesses the mechanical attribute of tissue.For shearing wave during propagation, amplitude, frequency and waveform can produce change,
The reliability and accuracy of shear wave velocity estimation can be used for weighing by these changes, and as a kind of method of auxiliary
To assess the mechanical attribute of tissue.
Based on above-mentioned principle, basic conception of the invention is：A kind of method and system of shearing wave elastogram are provided, will
The velocity information and dampening information of shearing wave are mapped to different color or gray scale, shown by specific image mapping algorithm
Show on display medium or device.Can not only be by the shear wave velocity information of tissue, i.e. group using this method and system
The shear wave velocity or Young's modulus knitted are shown, and can show the dampening information that the shearing wave of tissue is propagated
Come.Step is realized as shown in figure 1, including step：
S1, encouraged using acoustic radiation force in predeterminated position and produce shearing wave；
S2, repeatedly transmitting ultrasound detection wave beam detects shearing wave in shearing propagation path, to obtain shearing wave detection
Data；
S3, parsing shearing wave testing number obtain the velocity information and dampening information of shearing wave according to this；
S4, shear wave velocity information and dampening information, image corresponding to output are handled using default image mapping algorithm
Information.
It should be noted that a kind of shearing wave elastograph imaging method proposed by the present invention and system detect suitable for shearing wave
Single angle or multi-angle plane wave imaging, be also applied for mostly with multi-beam synthesis capability ultrasonic image-forming system on.
Specifically, as in Figure 2-4, in the present embodiment, shearing wave elastograph imaging method comprises the following steps：
A1, shearing wave elastogram area-of-interest (ROI frames, Region Of Interest) in, specific
Predeterminated position produces shearing wave using acoustic radiation force excitation.
As shown in Fig. 2 this implementation calculates shearing wave information matrix according to axial displacement or axial velocity data.It is described to cut
Cutting ripple information matrix includes shear wave velocity or Young's modulus, and shearing wave amplitude fading, shearing wave frequency decay and shearing
Ripple time domain waveform change etc..
Specifically, it can be determined according to the substantially propagation distance of shearing wave caused by the size and single activation of ROI frames
The number of ultrasonic acoustic radiation force excitation.If the length and width of ROI region is the substantially biography of shearing wave caused by single activation
Broadcast distance twice, then just at least need four excitations that shearing wave could be allowed to be full of whole ROI frames.Wherein, acoustic radiation is produced
The method of power is to launch the ultrasonic pulse of high-energy (generally hundreds of microseconds), then launches trace pulse and goes to cut caused by detection
Cut ripple.Transmitting and detection can be realized on same probe, and different probes can also be used to realize respectively.
A2, on the path that shearing wave is propagated, repeatedly launch ultrasound detection wave beam detection shearing wave, echo-signal used
Beam synthesis algorithm generation shearing wave detection scan line, and the shearing that multiple moment are generated in scan line is detected in each shearing wave
Ripple detects data.
Specifically, can repeatedly launch the larger ultrasound detection wave beam of coverage on the path that shearing wave is propagated,
This kind of ultrasonic beam can be with right and wrong focal plane ripple ultrasonic beam (single angle or multi-angle) or focusing ultrasonic wave
Beam.Wherein, Beam synthesis can substantially be divided into software Beam synthesis (ultra-fast imaging techniques) and hardware Beam synthesis (multi-beam
Synthetic technology).Software Beam synthesis can realize that single emission can rebuild the data of all passages, and hardware Beam synthesis list
Secondary transmitting can be only formed several scan lines.Using multi-beam synthetic technology, the shearing wave detection ultrasonic pulse of single emission is obtained
The echo information of wave beam, rebuild the data of a plurality of (being usually 4-8 bars) scan line.Detect the ultrasonic pulse wave beam transmitting of shearing wave
Number and repetition rate are required for depending on the Length x Width of ROI frames and shear wave velocity angle value.Such as：The length of ROI frames or
Width is bigger, then emitting times are more；Shear wave velocity angle value is bigger, then the repetition rate and emitting times of detection pulse accordingly
It will be correspondingly improved.As shown in figure 3, wherein, ultrasonic transducer 5 produces acoustic radiation force excitation ultrasonic pulse wave beam 1, acoustic radiation
Power excitation ultrasonic pulse wave beam 1 focuses on default focus 2, and so as to produce shearing wave in ROI frames 6,4 be the propagation of shearing wave
Wavefront, ultrasonic transducer 5 produce the ultrasonic pulse wave beam 3 for detecting shearing wave, the echo-signal of ultrasonic pulse wave beam 3 are made
4 shearing wave detection scan lines A, B, C and D are generated with Beam synthesis algorithm.The repeatedly ultrasonic arteries and veins of transmitting detection shearing wave
Wave beam is rushed, forms the time series data of any point in each scan line.
A3, using the shearing wave in shearing wave detection scan line data are detected, the detection scanning of each shearing wave is calculated
Shearing waveform data on line.
Calculate and caused tissue axial displacement or axial velocity propagated by shearing wave, be according to this be detected measuring point with same
Time series data on one position, the scan line formed at different moments calculates.This axial displacement or axial velocity
It can be the scan-line data stored in advance before being encouraged relative to ultrasonic acoustic radiation force to calculate, and utilize auto-correlation or cross correlation algorithm
It is calculated；Can also be same position, scan-line data at different moments carries out auto-correlation or cross correlation algorithm is calculated.
Can mutually it be changed by relational expression between axial velocity and axial displacement：Axial velocity is axial displacement on time dimension
Differential, axial displacement are integration of the axial velocity on time dimension.The axle calculated the multiecho of each scan line position
Combined sequentially in time to displacement or axial velocity data, form the shearing wave time waveform of the position.Such as Fig. 4 institutes
Show, 7 be scan line A locates in Fig. 2 shearing wave time waveform in Fig. 4, and 8 be the shearing wave time at scan line B in Fig. 3 in Fig. 4
Waveform, 9 be scan line C locates in Fig. 3 shearing wave time waveform in Fig. 4, in Fig. 4 10 for shearing wave at scan line D in Fig. 3 when
Between waveform.
A4, the shearing waveform data in scan line detected according to each shearing wave, the inspection of each shearing wave is calculated
The shearing wave for surveying scan line position propagates arrival time.Wherein it is determined that the method that shearing wave propagates arrival time has many kinds, most
Conventional is according to the time corresponding to the peak point of axial velocity data or axial displacement data, and it is shearing wave arrival should
The arrival time of position.
A5, the shearing wave propagation arrival time according to adjoining several shearing waves detection scan line position, calculate the position
The shear wave velocity put.
A6, the shearing waveform data according to adjoining several shearing waves detection scan line position, calculate shearing wave
Dampening information, dampening information here include shearing wave amplitude fading, shearing wave frequency decay and the change of shearing wave time domain waveform
The information such as change.
Ultrasonic acoustic radiation force excitation tissue produces shearing wave, and shearing wave is with the increase of propagation distance, the amplitude of its signal
Change can be produced with frequency and waveform, therefore, can only be shown by assessing these changes to supplement existing most of machine
The situation of this single parameter of shear wave velocity.
Specifically, by taking Fig. 4 as an example, for shearing wave amplitude fading, axial displacement data or axial velocity are first determined
The peak value of data, then find closest to peak value -3dB (peak valueTimes) or -6dB (1/2 times of peak value) cut-off
Position, take all average values more than or equal to rest position.Wherein, 11,12,13 and 14 in Fig. 3 are 7,8,9 and 10 respectively
Rest position.Exemplified by 7 in Fig. 4, calculating process is as follows：
Num=length (AmpA >=CutoffA) (formula 1)
Wherein, CutoffA is off position, and AmpA represents in Fig. 47 amplitude, and num represents that in Fig. 47 amplitude is more than and cut
The points (being calculated using length length functions) that stop bit is put, meanAmpA are represented for all width more than rest position
The average value (being calculated using mean mean value functions) of value.Similarly, 8,9 and 10 amplitude average value in Fig. 4 can be obtained.
Therefore, by the way that the result of calculation of adjacent cutouts ripple scan line is carried out into derivation operation, the rate of change of amplitude is assessed with this.Calculate
Process is as follows：
MeanAmpDiff (i)=meanAmp (i)-meanAmp (i+1) (formula 3)
Wherein, i represents the shearing wave detection number of scanning lines on horizontal direction, and meanAmp (i) and meanAmp (i+1) are respectively
The amplitude average value of adjacent position under a certain depth is represented, meanAmpDiff (i) represents that the amplitude of adjacent position under the depth is put down
The differentiation result of average, amplitude change rate is assessed with this.
And for shearing wave frequency decay, we can also be assessed with similar method.As shown in figure 5, wherein 15,
16th, 17 and 18 be respectively 7,8,9 and 10 Fourier transformation in Fig. 4 result.By Fourier transformation, signal can be extracted
Frequency content.The peak value of axial displacement data or axial velocity data Fourier transformation result is determined first, is then found most
Close to peak value -3dB (peak valueTimes) or -6dB (1/2 times of peak value) rest position, take all be more than or equal to
The average value of rest position.19,20,21 and 22 in wherein Fig. 5 be 15,16,17 and 18 rest position respectively.With in Fig. 5
15 exemplified by, calculating process is as follows：
FftAmpA=abs (fft (AmpA)) (formula 4)
FftNum=length (fftAmpA >=fftCutoffA) (formula 5)
Wherein, AmpA is 7 amplitude in Fig. 4, and fft (AmpA) is the result for seeking AmpA Fourier transformation,
FftCutoffA is off position, and fftAmpA represents in Fig. 5 15 amplitude (being calculated using abs ABS functions),
FftNum represents that in Fig. 5 15 amplitude is more than or equal to the points (being calculated using length length functions) of rest position,
FftMeanAmpA represents (to calculate using mean mean value functions for the average value of all amplitudes more than or equal to rest position
Obtain).Similarly, 16,17 and 18 amplitude average value can be obtained.Therefore, by by the calculating knot of adjacent cutouts ripple scan line
Fruit carries out derivation operation, and the rate of change of shearing wave frequency rate is assessed with this.
And the change of shearing wave time domain waveform can be weighed with the energy of signal.Assuming that take a certain axial location and horizontal position
The axial displacement or axial velocity signal put, because the length of time series is limited, that is to say, that the energy of the signal
It is limited.The physical significance of the energy of signal refers to the size of signal transmission capacity.If the frequency of signal is higher, decay is then got over
Rapidly, propagation distance is smaller；If the frequency of signal is lower, the energy having is bigger, and propagation distance is then more remote.With 7 in Fig. 4
Exemplified by, calculating process is as follows：
Wherein, AmpA is 7 amplitude in Fig. 4, and N is the length of shearing wave time domain waveform.| | expression takes absolutely to AmpA
To value, EnergyA represents the energy of the axial displacement that length of time series is N or axial velocity signal.Similarly, can obtain
16th, 17 and 18 energy.Therefore, by the way that the result of calculation of adjacent cutouts ripple scan line is carried out into derivation operation, assessed with this
The rate of change of signal energy.
According to the dampening information of shearing wave is calculated, dampening information here includes shearing wave amplitude fading, shearing wave frequency rate
The information such as decay and the change of shearing wave time domain waveform.
A7, repeat A2-A6, the shear wave velocity of the ROI inframes until obtaining covering whole shearing wave elastogram
Information matrix and shearing wave dampening information matrix.
Specifically, repeating A2-A6, i.e., by being calculated respectively along horizontal direction and axial direction, obtain covering
Cover the shear wave velocities of the ROI inframes of whole shearing wave elastogram, shearing wave amplitude fading, shearing wave frequency decay and cut
Cut the information matrixs such as ripple time domain waveform change.As shown in figure 3, because attenuation characteristic be present in the propagation of shearing wave, so shearing wave
Propagating a segment distance in a lateral direction can gradually decay afterwards；And shearing wave there is more than in a certain depth, but have
Certain axial range.So when we calculate shear wave velocity information and shearing wave dampening information, first same depth along
Horizontal direction is calculated different test points；Then again to next depth (axial direction) along horizontal direction to different detections
Point is calculated, by that analogy.So as to which the results of ROI frames all calculated, obtain covering whole shearing wave elasticity into
The shear wave velocity information matrix and shearing wave dampening information matrix of the ROI inframes of picture.
A8, by the obtained shear wave velocity information matrixs of A7 and shearing wave dampening information matrix, reflected by special image
Algorithm mapping is penetrated into different colours or different gray scales, is shown.
Wherein, image mapping algorithm can be that (color table/gray scale chart have recorded color/ash using color table or gray scale chart
The corresponding relation of degree and velocity information and dampening information) the table look-at mapping algorithm realized or utilize other calculating
Velocity information and dampening information (are such as carried out corresponding by form after the numerical transformation of specific formulation with color/gray-coded
With) realize image mapping algorithm.The shearing wave amplitude attenuation image of two dimension is shown, side reflection shearing can be used for
The accuracy and reliability of wave velocity estimation.By being calculated respectively along horizontal direction and axial direction, with gray scale or coloured silk
The method of color mapping, shearing wave frequency decay image is shown in a two-dimensional manner.
Specifically, result can be mapped to by different gray scales or different colours by special image mapping algorithm, with two
The form of dimension is shown.Shear wave velocity information matrix reflects the elastic information of test serum, and shearing wave amplitude attenuation,
Shearing wave frequency decay and the change of shearing wave time domain waveform accuracy and reliability that then reflection shear wave velocity in side is estimated, with
And can be as another method for assessing test serum.
Present invention also offers a kind of shearing wave elastogram system, including：Shearing wave encourages module, for utilizing acoustically radiating
Penetrate power and encourage generation shearing wave in predeterminated position；Detection module, for launching ultrasound detection wave beam in shearing propagation path
Shearing wave is detected, to obtain shearing wave detection data；Computing module, shearing wave is obtained according to this for parsing shearing wave testing number
Velocity information and dampening information；Image processing module, for utilizing default image mapping algorithm processing shear wave velocity information
And dampening information, image information corresponding to output.
In the present embodiment, the dampening information includes shearing wave amplitude fading information and/or shearing wave frequency decay information
And/or shearing wave time domain waveform change information, the computing module specifically include：Shear wave velocity information calculating sub module, use
Arrival time is propagated in the shearing wave that scan line position is detected according to adjacent shearing wave, calculates the shear wave velocity letter of the position
Breath；Shearing wave amplitude fading message sub-module, for determining the axial displacement data or axial velocity number of shearing waveform data
According to peak value, set an amplitude cut-off region, obtain the average value of the amplitude in amplitude cut-off region, more adjacent cuts
The average value of the amplitude of ripple detection scan line is cut to obtain shearing wave amplitude fading information；Shearing wave frequency decay information submodule
Block, for determining the peak value of the axial displacement data after Fourier transformation or axial velocity data, an amplitude cut-off region is set
Domain, obtains the average value of the amplitude in amplitude cut-off region, and the amplitude of more adjacent shearing wave detection scan line is averaged
Value is to obtain shearing wave frequency decay information；Shearing wave time domain waveform change information submodule, for determining shearing wave waveform number
According to axial displacement data or axial velocity data peak value, set an amplitude cut-off region, obtain in amplitude cut-off region
Interior amplitude square average value, more adjacent shearing wave detection scan line amplitude square average value cut with obtaining
Cut ripple time domain waveform change information.
In the present embodiment, a kind of course of work of shearing wave elastogram system and a kind of foregoing shearing wave elastogram side
Method is corresponding, will not be described here.
The present invention is shown by measuring velocity information and dampening information in shearing wave communication process, is overcome
Shearing wave elastogram systematic survey object present in prior art is single, causes organizational information to reflect that incomplete technology is asked
Topic, realize the shearing wave elastograph imaging method and system of a kind of speed for reflecting shearing wave and dampening information, Jin Eryou
Beneficial to the reliability and comprehensive for improving shearing wave detection.
In addition, calculating of the present invention for shear wave velocity information and shearing wave dampening information is all using by ripple
Shearing wave detection data after Shu Hecheng carry out auto-correlation, the axial velocity for the shearing wave that cross-correlation scheduling algorithm obtains or axial position
Data are moved, therefore can easily be implemented in systems.
It the composite can be widely applied to various shearing wave elastogram systems.
Above is the preferable implementation to the present invention is illustrated, but the invention is not limited to the implementation
Example, those skilled in the art can also make a variety of equivalent variations on the premise of without prejudice to spirit of the invention or replace
Change, these equivalent deformations or replacement are all contained in the application claim limited range.
Claims (10)
1. a kind of shearing wave elastograph imaging method, it is characterised in that including step：
S1, encouraged using acoustic radiation force in predeterminated position and produce shearing wave；
S2, repeatedly transmitting ultrasound detection wave beam detects shearing wave in shearing propagation path,
To obtain shearing wave detection data；
S3, parse velocity information and dampening information that the shearing wave testing number obtains shearing wave according to this；
S4, the velocity information and dampening information of the shearing wave are handled using default image mapping algorithm, export corresponding figure
As information.
2. a kind of shearing wave elastograph imaging method according to claim 1, it is characterised in that the step S3 is specially：
Using shearing wave detection data generation shearing waveform data, according to the shearing wave ripple of adjacent shearing wave detection scan line
Graphic data calculate shearing wave velocity information and dampening information, the dampening information include shearing wave amplitude fading information and/or
Shearing wave frequency decay information and/or shearing wave time domain waveform change information.
3. a kind of shearing wave elastograph imaging method according to claim 2, it is characterised in that the step S3 includes obtaining
Shear wave velocity information includes：
Arrival time is propagated according to the shearing wave of adjacent shearing wave detection scan line position, calculates the shear wave velocity of the position
Information.
4. a kind of shearing wave elastograph imaging method according to claim 2, it is characterised in that the step S3 includes obtaining
Shearing wave amplitude fading information includes：
It is determined that the axial displacement data of shearing waveform data or the peak value of axial velocity data,
One amplitude cut-off region is set, obtains the average value of the amplitude in amplitude cut-off region, more adjacent shearing wave
The average value of the amplitude of scan line is detected to obtain shearing wave amplitude fading information.
5. a kind of shearing wave elastograph imaging method according to claim 2, it is characterised in that the step S3 includes obtaining
Shearing wave frequency decay information attached bag includes：
Shearing waveform data is handled using Fourier transformation, determines axial displacement data or axial velocity after Fourier transformation
The peak value of data, an amplitude cut-off region is set,
The average value of the amplitude in amplitude cut-off region is obtained, the amplitude of more adjacent shearing wave detection scan line is averaged
Value is to obtain shearing wave frequency decay information.
6. a kind of shearing wave elastograph imaging method according to claim 2, it is characterised in that the step S2 includes obtaining
Shearing wave time domain waveform change information includes：
It is determined that the axial displacement data of shearing waveform data or the peak value of axial velocity data, set an amplitude cut-off region
Domain, obtain amplitude in amplitude cut-off region square average value, more adjacent shearing wave detects the amplitude of scan line
Square average value to obtain shearing wave time domain waveform change information.
A kind of 7. shearing wave elastograph imaging method according to any one of claim 1 to 6, it is characterised in that the step
S1 is specially：In the area-of-interest of shearing wave elastogram, shearing wave is produced using acoustic radiation force excitation in predeterminated position,
Ultrasonic acoustic radiation force excitation is determined according to the propagation distance of shearing wave caused by the size and single activation of area-of-interest
Number.
8. a kind of shearing wave elastograph imaging method according to claim 7, it is characterised in that the step S2 is specially：
Repeatedly transmitting ultrasound detection wave beam to obtain echo-signal, utilizes Beam synthesis algorithm process echo in shearing propagation path
Signal detects the shearing wave that multiple moment are generated in scan line in each shearing wave to generate multiple shearing wave detection scan lines
Detect data.
A kind of 9. shearing wave elastogram system, it is characterised in that including：
Shearing wave encourages module, and shearing wave is produced for being encouraged using acoustic radiation force in predeterminated position；
Detection module, for launching ultrasound detection wave beam detection shearing wave in shearing propagation path, to obtain shearing wave inspection
Survey data；
Computing module, obtain the velocity information and dampening information of shearing wave according to this for parsing shearing wave testing number；
Image processing module, for utilizing default image mapping algorithm processing shear wave velocity information and dampening information, output
Corresponding image information.
10. a kind of shearing wave elastogram system according to claim 9, it is characterised in that the dampening information includes
Shearing wave amplitude fading information and/or shearing wave frequency decay information and/or shearing wave time domain waveform change information, the calculating
Module specifically includes：
Shear wave velocity information calculating sub module, the shearing wave for detecting scan line position according to adjacent shearing wave travel to
Up to the time, the shear wave velocity information of the position is calculated；
Shearing wave amplitude fading message sub-module, for determining the axial displacement data or axial velocity number of shearing waveform data
According to peak value, set an amplitude cut-off region, obtain the average value of the amplitude in amplitude cut-off region, more adjacent cuts
The average value of the amplitude of ripple detection scan line is cut to obtain shearing wave amplitude fading information；
Shearing wave frequency decay message sub-module, for determining axial displacement data or axial velocity data after Fourier transformation
Peak value, an amplitude cut-off region is set, obtains the average value of the amplitude in amplitude cut-off region, more adjacent shearing
The average value of the amplitude of ripple detection scan line is to obtain shearing wave frequency decay information；
Shearing wave time domain waveform change information submodule, for determining the axial displacement data or axially fast of shearing waveform data
The peak value of degrees of data, an amplitude cut-off region is set, obtain amplitude in amplitude cut-off region square average value, than
The amplitude of more adjacent shearing wave detection scan line square average value to obtain shearing wave time domain waveform change information.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710859560.6A CN107510474B (en) | 2017-09-21 | 2017-09-21 | Shear wave elastic imaging method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710859560.6A CN107510474B (en) | 2017-09-21 | 2017-09-21 | Shear wave elastic imaging method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107510474A true CN107510474A (en) | 2017-12-26 |
| CN107510474B CN107510474B (en) | 2020-07-10 |
Family
ID=60726300
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710859560.6A Active CN107510474B (en) | 2017-09-21 | 2017-09-21 | Shear wave elastic imaging method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107510474B (en) |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107970043A (en) * | 2017-12-28 | 2018-05-01 | 深圳开立生物医疗科技股份有限公司 | A kind of detection method and device of shearing wave |
| CN110243950A (en) * | 2019-05-15 | 2019-09-17 | 天津大学 | Acoustic radiation force impulse Excitation Method for Frequency based on amplitude modulation |
| CN110494082A (en) * | 2018-04-11 | 2019-11-22 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic elastograph imaging method and system |
| CN110477947A (en) * | 2019-08-14 | 2019-11-22 | 中国科学院苏州生物医学工程技术研究所 | Plane wave beam synthesizing method based on deep learning, system, storage medium, equipment |
| CN110573084A (en) * | 2018-08-29 | 2019-12-13 | 深圳迈瑞生物医疗电子股份有限公司 | ultrasonic elasticity detection equipment and shear wave elasticity imaging method and device |
| CN110604598A (en) * | 2018-12-04 | 2019-12-24 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method and ultrasonic imaging system |
| CN110742644A (en) * | 2019-09-29 | 2020-02-04 | 深圳大学 | Elastography system, elastography method and storage medium |
| WO2020140917A1 (en) * | 2019-01-02 | 2020-07-09 | 无锡海斯凯尔医学技术有限公司 | Tissue stiffness measuring method, device, and system |
| CN111449681A (en) * | 2020-04-08 | 2020-07-28 | 深圳开立生物医疗科技股份有限公司 | Shear wave imaging method, device and equipment and readable storage medium |
| CN112135567A (en) * | 2018-05-03 | 2020-12-25 | 皇家飞利浦有限公司 | Shear wave amplitude reconstruction for tissue elasticity monitoring and display |
| CN112244889A (en) * | 2020-10-15 | 2021-01-22 | 中国科学院苏州生物医学工程技术研究所 | Method for determining vibration element array, puncture needle imaging method and ultrasonic equipment |
| CN112799078A (en) * | 2021-04-15 | 2021-05-14 | 深圳中科乐普医疗技术有限公司 | Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment |
| CN113081054A (en) * | 2018-12-04 | 2021-07-09 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method and ultrasonic imaging system |
| CN114072066A (en) * | 2019-09-27 | 2022-02-18 | 深圳迈瑞生物医疗电子股份有限公司 | Elasticity imaging method, system and computer readable storage medium |
| CN114403920A (en) * | 2022-01-24 | 2022-04-29 | 陕西师范大学 | Time-division multipoint excitation shear wave fusion imaging method and system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5606971A (en) * | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
| CN103054552A (en) * | 2012-12-24 | 2013-04-24 | 深圳先进技术研究院 | Method and system for measuring biological tissue viscoelasticity |
| CN103462643A (en) * | 2013-09-29 | 2013-12-25 | 深圳市开立科技有限公司 | Shear wave speed measurement method, device and system |
| US20140088421A1 (en) * | 2012-09-27 | 2014-03-27 | Mayo Foundation For Medical Education And Research | System and method for local estimation of nonlinear tissue elasticity with acoustic radiation force |
| CN104825195A (en) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | Shear wave viscoelasticity imaging method and system |
| CN106175831A (en) * | 2016-09-28 | 2016-12-07 | 重庆理工大学 | A kind of tissue viscoelasticity based on shearing wave amplitude and phase-detection measuring method |
-
2017
- 2017-09-21 CN CN201710859560.6A patent/CN107510474B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5606971A (en) * | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
| US20140088421A1 (en) * | 2012-09-27 | 2014-03-27 | Mayo Foundation For Medical Education And Research | System and method for local estimation of nonlinear tissue elasticity with acoustic radiation force |
| CN103054552A (en) * | 2012-12-24 | 2013-04-24 | 深圳先进技术研究院 | Method and system for measuring biological tissue viscoelasticity |
| CN103462643A (en) * | 2013-09-29 | 2013-12-25 | 深圳市开立科技有限公司 | Shear wave speed measurement method, device and system |
| CN104825195A (en) * | 2015-05-25 | 2015-08-12 | 无锡海斯凯尔医学技术有限公司 | Shear wave viscoelasticity imaging method and system |
| CN106175831A (en) * | 2016-09-28 | 2016-12-07 | 重庆理工大学 | A kind of tissue viscoelasticity based on shearing wave amplitude and phase-detection measuring method |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107970043A (en) * | 2017-12-28 | 2018-05-01 | 深圳开立生物医疗科技股份有限公司 | A kind of detection method and device of shearing wave |
| CN110494082A (en) * | 2018-04-11 | 2019-11-22 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic elastograph imaging method and system |
| CN112135567A (en) * | 2018-05-03 | 2020-12-25 | 皇家飞利浦有限公司 | Shear wave amplitude reconstruction for tissue elasticity monitoring and display |
| CN110573084B (en) * | 2018-08-29 | 2022-07-08 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic elasticity detection equipment and shear wave elasticity imaging method and device |
| CN110573084A (en) * | 2018-08-29 | 2019-12-13 | 深圳迈瑞生物医疗电子股份有限公司 | ultrasonic elasticity detection equipment and shear wave elasticity imaging method and device |
| WO2020042020A1 (en) * | 2018-08-29 | 2020-03-05 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic elasticity detection device and shear wave elasticity imaging method and apparatus |
| CN110604598A (en) * | 2018-12-04 | 2019-12-24 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method and ultrasonic imaging system |
| CN113081054A (en) * | 2018-12-04 | 2021-07-09 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging method and ultrasonic imaging system |
| CN111388012B (en) * | 2019-01-02 | 2021-10-08 | 无锡海斯凯尔医学技术有限公司 | Method, device and system for detecting tissue hardness |
| WO2020140917A1 (en) * | 2019-01-02 | 2020-07-09 | 无锡海斯凯尔医学技术有限公司 | Tissue stiffness measuring method, device, and system |
| CN111388012A (en) * | 2019-01-02 | 2020-07-10 | 无锡海斯凯尔医学技术有限公司 | Method, device and system for detecting tissue hardness |
| CN110243950B (en) * | 2019-05-15 | 2021-10-08 | 天津大学 | Acoustic radiation force pulse excitation method based on amplitude modulation |
| CN110243950A (en) * | 2019-05-15 | 2019-09-17 | 天津大学 | Acoustic radiation force impulse Excitation Method for Frequency based on amplitude modulation |
| CN110477947A (en) * | 2019-08-14 | 2019-11-22 | 中国科学院苏州生物医学工程技术研究所 | Plane wave beam synthesizing method based on deep learning, system, storage medium, equipment |
| CN114072066A (en) * | 2019-09-27 | 2022-02-18 | 深圳迈瑞生物医疗电子股份有限公司 | Elasticity imaging method, system and computer readable storage medium |
| CN114072066B (en) * | 2019-09-27 | 2024-06-07 | 深圳迈瑞生物医疗电子股份有限公司 | Elastography method, elastography system and computer-readable storage medium |
| CN110742644A (en) * | 2019-09-29 | 2020-02-04 | 深圳大学 | Elastography system, elastography method and storage medium |
| CN110742644B (en) * | 2019-09-29 | 2022-07-08 | 深圳大学 | Elastography system, elastography method and storage medium |
| CN111449681A (en) * | 2020-04-08 | 2020-07-28 | 深圳开立生物医疗科技股份有限公司 | Shear wave imaging method, device and equipment and readable storage medium |
| CN111449681B (en) * | 2020-04-08 | 2023-09-08 | 深圳开立生物医疗科技股份有限公司 | Shear wave imaging method, device, equipment and readable storage medium |
| CN112244889A (en) * | 2020-10-15 | 2021-01-22 | 中国科学院苏州生物医学工程技术研究所 | Method for determining vibration element array, puncture needle imaging method and ultrasonic equipment |
| CN112799078A (en) * | 2021-04-15 | 2021-05-14 | 深圳中科乐普医疗技术有限公司 | Detection method and system for shear wave propagation velocity and ultrasonic imaging equipment |
| CN114403920A (en) * | 2022-01-24 | 2022-04-29 | 陕西师范大学 | Time-division multipoint excitation shear wave fusion imaging method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107510474B (en) | 2020-07-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107510474A (en) | Shearing wave elastograph imaging method and system | |
| CN102551791B (en) | A kind of ultrasonic imaging method and device | |
| CN103462643B (en) | Shear wave speed measurement method, device and system | |
| CN105816205B (en) | Sparse tracking in acoustic radiation force Pulse Imageing | |
| Jensen | Users’ guide for the Field II program | |
| KR101922522B1 (en) | Sound speed imaging using shear waves | |
| EP2356941A1 (en) | Ultrasonographic system and ultrasonographic device | |
| CN104605890B (en) | A kind of method of shearing wave peak value waveform correction, device, system and its application | |
| CN107613878A (en) | Ultrasonic image-forming system and method for detection object motion | |
| JPH0120899B2 (en) | ||
| JP7340868B2 (en) | Ultrasonic method and equipment | |
| US20240245391A1 (en) | Shear wave imaging based on ultrasound with increased pulse repetition interval | |
| CN110301936A (en) | Frequency scanning for acoustic radiation force pulse | |
| CN102365054A (en) | Ultrasonograph and urtrasonic processing device | |
| CN108338808B (en) | Shear velocity imaging using coherence | |
| CN110893103A (en) | Angle for ultrasound-based shear wave imaging | |
| CN110243950A (en) | Acoustic radiation force impulse Excitation Method for Frequency based on amplitude modulation | |
| CN108852416A (en) | A kind of determination method and device for shearing velocity of wave propagation | |
| CN109259801A (en) | A kind of shearing wave elastograph imaging method and device | |
| WO2020197646A1 (en) | Asymmetry for acoustic radiation force impulse | |
| CN1095356C (en) | Mutual correlation-spectrum figure method in supersonic blood current measurement | |
| JP7302651B2 (en) | Ultrasonic signal processing device, ultrasonic diagnostic device, ultrasonic signal processing method, and program | |
| Coetzee | The physics of ultrasound and Doppler | |
| CN108982674A (en) | Wave property is sheared using single driving pulse and with different in width and multiple driving pulses of shape measurement tissue | |
| Li et al. | Nonlinear propagation in Doppler ultrasound |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |