CN106618639A - Quantitative shear wave elastography method - Google Patents
Quantitative shear wave elastography method Download PDFInfo
- Publication number
- CN106618639A CN106618639A CN201610973110.5A CN201610973110A CN106618639A CN 106618639 A CN106618639 A CN 106618639A CN 201610973110 A CN201610973110 A CN 201610973110A CN 106618639 A CN106618639 A CN 106618639A
- Authority
- CN
- China
- Prior art keywords
- wave
- shear
- point
- wave velocity
- quantitative
- 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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Vascular Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The invention relates to the technical field of medical ultrasonic imaging, in particular to a quantitative shear wave elastography method. According to the ultrasonic quantitative elastography method and system based on strain of sliding window linear fitting and a shear wave velocity detection algorithm using two-dimensional linear fitting, higher noise immunity can be achieved, and a more reliable result can be obtained; meanwhile, under the condition that the load of an ultrasonic front-end storage and transmission module is not increased, ultrasonic whole-domain quantitative elastography is achieved, and the design difficulty of the ultrasonic quantitative elastography system and the equipment cost are greatly lowered.
Description
Technical field
The present invention relates to Medical ultrasound image technology field, more particularly to a kind of quantitative shearing wave elastograph imaging method.
Background technology
In recent years, tumor disease has become the healthy primary killers of harm our people.The preventing and treating of tumor disease
It is critical only that early detection.Clinically, often tumor disease of paramount importance early stage is pre- for the mechanical characteristics change of tissue
Alert signal, such as mammary gland, thyroid gland neoplastic lesion.With the growth of tumour, pathological tissues compared with normal structure, send out by viscoelasticity
Raw large change, such as breast cancer about differ 90 times with normal surrounding tissue, and fibrosed tissue, non-invasive carcinoma become and wellability
Also there is larger difference in the shearing elasticity coefficient of cancerous issue, therefore, obtain the tissue parameter information related to viscoplasticity
(such as deformational displacement, shear wave velocity etc.) is early stage area of medical diagnostics is especially to breast cancer, thyroid cancer and prostate cancer
Detection is significant.
Medical ultrasonic elastogram be it is a kind of by histokinesis show the ultrasound of tissue elasticity modulus or soft or hard degree into
Image space formula, traditional squash type elastogram mode needs doctor to make its movement under force using probe pinch detection position, so as to
Tissue elasticity information is obtained, the deficiency of this mode is：Doctor is needed with appropriate operation pinch detection position, subjectivity
By force；Due to needing doctor to extrude, it is right that the result of different doctors operation or the operation of same doctor's different times is difficult directly to carry out
Than testing result is difficult to reappear, often can only be used as qualitatively testing result, it is impossible to quantitative detection information is obtained, to disease
The tracking of disease and clinical follow.
The Ultrasonic Elasticity Imaging of quantitative analysis can be carried out at present, be mainly based upon acoustic radiation force elastogram
(Acoustic Radiation Force Imaging, ARFI), ARFI is using the focusing ultrasound in medical ultrasonic power bracket
Driving pulse produces acoustic radiation force in bioviscoelasticity tissue focal zone, deforms upon tissue, then using detection pulse
Echo-signal the deformational displacement situation of tissue is detected in different time points by the time-delay calculation method of correlation, be imaged
Qualitatively reflect the Viscoelastic Properties of tissue afterwards.ARFI elastograms overcome conventional freehand elastogram cannot be from vitro to depth
Portion's tissue effectively presses and because the use habit difference of operator causes repeatable poor shortcoming.But it is disadvantageous in that
Image is shown the relative deformational displacement difference of tissue, and cannot quantitatively estimate the elastic modelling quantity of tissue completely.It is based on
Ultrasonic shear wave imaging (SWI, shear wave imaging) derived from ARFI technologies, by adjusting pulse excitation mode, root
The shear-wave velocity of horizontal transmission is estimated according to the displacement-time curve of the multiple points of horizontal direction, then propagates special using shearing wave
Property and biological tissue elasticity feature between inner link, the elastic modelling quantity of final quantitative reconstruct tissue is forming the figure of two dimension
Picture.
But, existing quantitative elastogram SWI is needed using the original radio frequency signal number before the Beam synthesis of ultrasonic front end
According to, needing to launch based on plane wave and receive, plane wave technology is per second to need N*c/2z calculating, and wherein N is parallel reception sound
Beam number, c is the velocity of sound, and z is scan depths.And the amount of calculation of traditional beam synthesizing method is then much smaller, only M/P/
One, wherein M are emitting times, and P is simultaneously line number, such as when M is equal to 100, P equal to 2, the amount of calculation of plane-wave method is tradition
At least 50 times (sometimes even up to hundreds times) of ARFI methods；In addition to the greatest differences of amount of calculation, based on plane wave skill
The receiving terminal software Beam synthesis of art also need to store a large amount of original echoed signals, data storage is greatly improved and transmission is hard
Part cost, is not easy to and popularization and application integrated with medical ultrasound system.
The content of the invention
It is an object of the invention to overcome existing quantitative elastography to need using original radio frequency signal (without ripple
Computationally intensive, necessary hardware difficult design Shu Hecheng) brought and problem with high costs, there is provided it is a kind of it is low for equipment requirements,
The little ultrasonic quantitative elastograph imaging method of amount of calculation and system.
In order to realize foregoing invention purpose, the invention provides technical scheme below：
A kind of quantitative shearing wave elastograph imaging method, comprises the steps of：
ARFI detections are carried out, each position strain stress in target area is obtained, wherein comprising specified reference position strain stressref；
Shearing wave detection is carried out, is obtained and specify in target area reference position single-point shear-wave velocity cref；
Each position shear-wave velocity c is calculated according to the single-point shear-wave velocity；
The quantitative shearing wave elastic image E in each position is obtained according to each position strain gauge shear-wave velocity.
Further, each position strain in the target areaWherein, N is gradient cunning
Dynamic calculation window size, ziThe concrete depth coordinate of each position is represented,Represent each position coordinates in gradient slip calculation window
Average, yiCurrent point displacement in gradient slip calculation window is represented,Represent each position displacement in gradient slip calculation window
Average.
Further, lateral separation, axial distance of the single-point shear-wave velocity according to specified measurement point, and each refers to
Determine the time that measurement point reached required for maximum transversal displacement is calculated using two-dimensional linear fitting process.
Further, the single-point shear-wave velocity crefDrawn by formula X=A β, wherein, In formula, xiRepresent the lateral separation of specified measurement point, tiRepresent that specified measurement point reaches most
The time required to big displacement, ziRepresent the axial distance of specified measurement point, β1As required specified measurement point single-point shear wave velocity
cref, β0And β2It is to calculate single-point shear wave velocity c using above formularefDuring the regression parameter that calculates simultaneously, in this method not
Using the two parameters.
Further, described each position shear-wave velocityWherein, in focal zone, γ values 1；
Axial region outside focal zone,Wherein, z is axial distance, and n and σ is to represent that the system of excitation sound field is normal
Number, span is the value of two parameters of 0～10, n and σ can be with identical or different.
Further, the focal zone DOF=8 (f#)2λ, wherein, λ represents wavelength, and f# is aperture control parameter, its
It is the real number between 0～5, it is seen that focal zone is generally determined that its magnitude range is only related to device parameter by device parameter,
Unrelated with target area size, target area is likely larger than focal zone, it is also possible to equal to focal zone, or less than focal zone
Domain.
Further, the quantitative shearing wave elastic image E=3 ρ c2, wherein, ρ is tested region Media density.
Above method may operate in following ultrasonic quantitative elastogram system, and the system includes,
Shearing pulse generator, for launching shearing wave driving pulse and detection pulse signal；
ARFI impulse generators, for launching ARFI driving pulses and detection pulse signal；
Probe, for receiving shearing wave echo-signal or ARFI echo-signals；
Control device, for control shearing pulse generator, the switch of ARFI impulse generators, receives shearing wave echo
Signal and ARFI echo-signals, and, the strain stress of each position of target area is calculated according to ARFI echo-signals, wherein comprising finger
Determine reference position strain stressref；Calculated according to shearing wave echo-signal and specify reference position single-point shear-wave velocity cref；According to list
Point shear-wave velocity calculates each position shear-wave velocity c；According to each position strain gauge shear-wave velocity c and strain stressrefObtain
Take the quantitative shearing wave elastic image E in each position；
Display device, for the quantitative shearing wave elastic image E to be shown.
The probe includes beam synthesizer, and it is used to close FRFI echo-signals or shearing wave echo-signal through wave beam
Into synthesizing echo radiofrequency signal.
The control device includes shearing wave signal processing module, single-point shear wave velocity computing module；
The shearing wave signal processing module receives shearing wave echo radiofrequency signal from popping one's head in, and demodulated, displacement is estimated
Calculation draws tested tissue target area each position or only obtains the shearing wave displacement versus time data for specifying reference position；
The single-point shear wave velocity computing module is calculated according to shearing wave displacement versus time data and specifies reference position shearing wave
Velocity of wave cref, it passes through formula X=A β and draws, wherein,In formula, xiRepresent
The lateral separation of specified measurement point (specified measurement point includes specified reference position), tiRepresent that specified measurement point reaches maximum displacement
Required time, ziRepresent the axial distance of specified measurement point, β1As required specified measurement point single-point shear wave velocity cref；
The control device also includes ARFI signal processing modules, strain calculation module；
The ARFI signal processing modules are used to receive ARFI echo radiofrequency signals from probe, and demodulated, displacement is estimated
Calculation draws the ARFI displacement versus time data of tested tissue target area each position；
The strain calculation module calculates the strain of target area each position according to the ARFI displacement versus times dataWherein, N be gradient slip calculation window size, ziThe concrete depth coordinate of each position is represented,
The average of each position coordinates in gradient slip calculation window is represented, yi represents current point displacement in gradient slip calculation window,The average of each position displacement in gradient slip calculation window is represented, comprising the strain stress for specifying reference position in εref。
The control device also includes elastic modelling quantity computing module；
The control device also includes velocity of wave computing module and elastic modelling quantity computing module；
The velocity of wave computing module is used to calculate each position shearing wave according to specified reference position single-point shear-wave velocity
Velocity of waveWherein, in focal zone, γ values 1；Axial region outside focal zone,Its
In, z is axial distance, and n and σ is the system constants for representing excitation sound field, and span is two parameters of 0～10, n and σ
Value can be with identical or different；
The elastic image computing module is used for according to c and strain stressrefObtain the quantitative shearing wave elastic image E in each position
=3 ρ c2, wherein, ρ is tested region Media density.
The control device control is described to shear pulse generator, ARFI impulse generator alternate emissions.
The shearing pulse generator and the ARFI impulse generators are same circuit realiration, and it is filled by the control
Put control and send shearing wave impulse or ARFI pulses；,
The shearing pulse generator and the ARFI impulse generators are the realization of respective independent circuits, and it is by the control
Device control processed starts, the shut-in time.
The shearing pulse generator and/or ARFI impulse generators have more than 24 tunnels.
The shearing pulse generator or ARFI impulse generators include driving amplifier, the pulse generation being sequentially connected in series
Device, overcurrent-overvoltage protecting circuit, Tx/Rx switches, and the limiting amplifier from Tx/Rx switch feedbacks, D/A converting circuit.
Compared with prior art, beneficial effects of the present invention：
The ultrasonic quantitative elastograph imaging method that the present invention is provided is based on the strain of sliding window linear fit and using two dimension
The shear-wave velocity detection algorithm of linear fit, its result robustness is more preferable, noise resisting ability is higher, as a result more reliable.Simultaneously
In the case where ultrasonic front end storage and transport module load is not additionally increased, the quantitative elastogram of ultrasonic universe is realized, greatly
Reduce the design difficulty and equipment cost of ultrasonic quantitative elastogram system width.
Efficient circuit switching and protection module are provided with the system of the inventive method operation, allow system to meet ARFI and cut
While cutting ripple different transmitting timing requirements, effective protection human body and circuit system, and acoustically-driven is carried out to greatest extent to reach
To more preferable Detection results.
Description of the drawings：
The ultrasonic quantitative elastograph imaging method flow chart that Fig. 1 is provided for the present invention.
Fig. 2 is the ultrasonic quantitative elastogram system block diagram that the present invention is provided.
Fig. 3 is the structural frames that pulse generator and/or ARFI impulse generators composition example are sheared in the present invention
Figure.
Fig. 4 is the transmitting pulse train schematic diagram that acoustically-driven is used in the present invention in the present invention.
Fig. 5 is that shear-wave velocity calculates the impulse ejection schematic diagram for using.
Specific embodiment
Below in conjunction with the accompanying drawings and specific embodiment the present invention is described in further detail.But this should not be interpreted as this
The scope for inventing above-mentioned theme is only limitted to below example, and all technologies realized based on present invention belong to the present invention
Scope.
Embodiment 1：As shown in figure 1, the present invention provides a kind of quantitative shearing wave elastograph imaging method, comprise the steps of：
S100：ARFI detections are carried out, each position strain stress in target area, in the present embodiment, described each position is obtained
Can be that having in a position, or target area for target area specifies the aiming spot of number (such as target area
N number of equally distributed position), or the N number of regularly arranged point position chosen of target area (regularly arranged can be with
It is equidistant arrangement, or the Unequal distance arrangement of rule, each point may be in line, it is also possible to line up other regular shapes)；No matter
Described each position refers to target area Suo Youdian positions, or N number of equally distributed aiming spot, or N number of rule row
The point position of row, wherein specify reference position comprising at least one, that is, comprising specified reference in each position strain stress for obtaining
Position strain stressref；
S200：Shearing wave detection is carried out, is obtained and specify in target area reference position single-point shear-wave velocity cref；
The execution sequence of step S100 and S200 has no particular/special requirement, can such as first carry out S100, then performs S200, also may be used
To first carry out S200, then perform S100.
S300：Each position shear-wave velocity c is calculated according to the single-point shear-wave velocity；
S400：The quantitative shearing wave elastic image E in each position is obtained according to each position strain gauge shear-wave velocity.
Specifically, each position strain in the target areaWherein, N is gradient slip
Calculation window size, ziThe concrete depth coordinate of each position is represented,Represent the equal of interior each position coordinates of gradient slip calculation window
Value, yiCurrent point displacement in gradient slip calculation window is represented,Represent the equal of interior each position displacement of gradient slip calculation window
Value.
Lateral separation, axial distance of the single-point shear-wave velocity according to specified measurement point, and each specified measurement point
Reach the time required for maximum transversal displacement is calculated using two-dimensional linear fitting process.
The single-point shear-wave velocity crefDrawn by formula X=A β, wherein, In formula, generally we preset I depth in axial depth, and J specified location (point) is preset in each depth,
Altogether, m specified location (point), m=I*J；Each position (point) is issued to the time of maximum displacement in shearing wave effect for we
Regard that shearing wave travels to the time of the position (point) as, by shearing wave displacement versus time data shearing wave transmission can be obtained
To the time of the position (point), i.e. ti；Meanwhile, ziRepresent that specified measurement position (point) (is characterized deep in the axial distance of axial depth
Degree position), xiRepresent the lateral separation (sign lateral attitude) of specified measurement position (point), i.e. when each specified location (point)
When selected, each element value and z in XiValue is just it has been determined that we only need to find out each from shearing wave displacement versus time data
Specified location (point) reaches maximum displacement required time tiCan complete to calculate；In formula, β1As required specified measurement point list
Point shear wave velocity cref；Simultaneously, it is noted that β0And β2It is to calculate single-point shear wave velocity c using above formularefDuring count simultaneously
The regression parameter for drawing, uses the two parameters in this method, its effect and implication are to the no essential meaning of the present invention.
Described each position shear-wave velocityWherein, in focal zone γ values 1；Focal zone it
Outer axial region (when target area is more than focal zone, need to use herein),Wherein, z is axial distance,
N and σ are the system constants for representing excitation sound field, and span is the value of two parameters of 0～10, n and σ can be with identical or not
Together.
The quantitative shearing wave elastic image E=3 ρ c2, wherein, ρ is tested region Media density.
The method that the present embodiment is provided may operate in ultrasonic quantitative elastogram system as shown in Figure 2 and Figure 3, should
System includes,
Shearing pulse generator, for launching shearing wave driving pulse and detection pulse signal；
ARFI impulse generators, for launching ARFI driving pulses and detection pulse signal；
Probe, for receiving shearing wave echo-signal or ARFI echo-signals；
Control device 1, for control shearing pulse generator, the switch of ARFI impulse generators, receives shearing wave and returns
Ripple signal and ARFI echo-signals, and, the strain stress of each position of target area is calculated according to ARFI echo-signals, wherein including
Specified reference position strain stressref；Calculated according to shearing wave echo-signal and specify reference position single-point shear-wave velocity cref；According to
Single-point shear-wave velocity calculates each position shear-wave velocity c；According to each position strain gauge shear-wave velocity c and strain stressref
Obtain the quantitative shearing wave elastic image E in each position.
Display device, for the quantitative shearing wave elastic image E to be shown.Specifically, before display, it is right also to include
The elastic image E data that control device is calculated are scanned a series of processing procedure such as conversion, smoothing processing, finally will
The image is presented to user by display.
The probe includes beam synthesizer 2, and it is used for ARFI echo-signals or shearing wave echo-signal through wave beam
Synthesis, synthesizes echo radiofrequency signal.
The control device 1 includes shearing wave signal processing module 13, single-point shear wave velocity computing module 14；
The shearing wave signal processing module 13 receives shearing wave echo radiofrequency signal from popping one's head in, and demodulated becomes IQ
Data, further carry out the shearing wave displacement versus time that offset estimation draws tested tissue target area each position to the I/Q data
Data；
The single-point shear wave velocity computing module 14 is calculated according to shearing wave displacement versus time data specifies reference position shearing
Ripple velocity of wave cref, it passes through formula X=A β and draws, wherein,In formula, xiTable
Show the lateral separation of specified measurement point, tiThe time required to representing that specified measurement point reaches maximum displacement, ziRepresent specified measurement point
Axial distance, β1As required specified measurement point single-point shear wave velocity cref；
The control device 1 also includes ARFI signal processing modules 11, strain calculation module 12；
The ARFI signal processing modules 11 are used to receive ARFI echo radiofrequency signals, and demodulated formation IQ from probe
Data, and further by I/Q data is carried out offset estimation draw tested tissue target area each position ARFI displacement-when
Between data；
The strain calculation module 12 calculates the strain of target area each position according to the ARFI displacement versus times dataWherein, N be gradient slip calculation window size, zi(each position herein is led to represent each position
It is also often to choose from target area) concrete depth coordinate,Represent the equal of interior each position coordinates of gradient slip calculation window
Value, yiCurrent point displacement in gradient slip calculation window is represented,Represent the equal of interior each position displacement of gradient slip calculation window
Value.
The control device 1 also includes velocity of wave computing module 15 and elastic modelling quantity computing module 16；
The velocity of wave computing module 15 is used to calculate each position shearing according to specified reference position single-point shear-wave velocity
Ripple velocity of waveWherein, in focal zone, γ values 1；Axial region outside focal zone,
Wherein, z is axial distance, and n and σ is the system constants for representing excitation sound field, and span is two parameters of 0～10, n and σ
Value can be with identical or different；
The elastic image computing module 16 is used for according to c and strain stressrefObtain the quantitative shearing wave elastic graph in each position
As E=3 ρ c2, wherein, ρ is tested region Media density.
The control device 1 controls the shearing pulse generator, ARFI impulse generator alternate emissions.
Specifically, as described in Figure 3, in some embodiments, the hardware components of control device 1 can by PC, control board and
Fpga chip group is constituted, and in other embodiment, control device can also be completely integrated integrated circuit together and constitute,
Above-mentioned shearing wave signal processing module, single-point shear wave velocity computing module, ARFI signal processing modules, strain calculation module,
Elastic modelling quantity computing module is distributed in above-mentioned hardware unit by function.
The shearing pulse generator and the ARFI impulse generators are same circuit realiration, and it is filled by the control
Put control and send shearing wave impulse or ARFI pulses；Or,
The shearing pulse generator and the ARFI impulse generators are the realization of respective independent circuits, and it is by the control
Device control processed starts, the shut-in time.
Specifically, in the shearing pulse generator and enforcement that the ARFI impulse generators are same circuit realiration
In example, as shown in figure 3, generator includes driving amplifier 100, impulse generator 101, the over-current over-voltage protection electricity being sequentially connected in series
Road 102, Tx/Rx switch 103, and the limiting amplifier 104 from Tx/Rx switch feedbacks, D/A converting circuit 105.Cross and flow through
Voltage protection circuit 102 can allow system meet ARFI it is different from shearing wave transmitting timing requirements while, effective protection human body and
Circuit system, and acoustically-driven is carried out to greatest extent to reach more preferable Detection results.
The shearing pulse generator and/or ARFI impulse generators have more than 24 tunnels, such as 48 tunnels, 64 tunnels, 128
Road, 256 tunnels etc..
When using, the control pulse of the parameter information such as voltage, pulse length and phase place that control device sets according to user is sent out
Raw device (shearing pulse generator and/or the ARFI impulse generators) produces excitation and detects pulse signal, in T/R switches
Control under the ultrasonic wave that produces of radiating circuit enter biological tissue, receiving circuit receives echo-signal.System is being launched or is being connect
During the collection of letters, by arranging aperture control parameter f# and depth of focus z, and formula is utilizedTo control activity array element (arteries and veins
Rush generator) number, in formula, D is probe pore size, i.e. user by arranging aperture control parameter f# and the depth of focus
, come the control that controls the size in movable probe aperture to realize to movable element number of array, user can according to circumstances in maximum activity for z
The activity array element quantity being actually needed is selected in array number (such as 64 tunnels), for example, if system includes No. 64 impulse generators,
The maximum probe pore size that system can be realized is the area that 64 array elements are distributed, and aperture control parameter f# is usually big
In 0 real number less than or equal to 5.Because it is the option transmitting set according to user by probe that excitation tissue produces the power of deformation
Short duration pulse ripple is formed in focal zone, unrelated with the operation technique of operator, it can be considered that in effective imaging area
The size of domain internal force be it is uniform, it is consistent.The depth (Depth of Focus) of focal zone can be described as：DOF=8
(f#)2λ, λ represent wavelength；For example, the driving pulse for 5MHz, if aperture control parameter f#=2, the sound in biological tissue
Speed is approximately equal to 1540m/s, then the close 1cm of value of DOF.ARFI elastograms use by-line mode and scan, each
Horizontal level all can launch driving pulse and cause tissue to produce miniature deformation according to identical systematic parameter, it can be considered that
Acoustic radiation force in the range of near focal point 1cm can be regarded as uniformly, consistent, and we are referred to as focal zone in this region,
As described above, in focal zone, shear-wave velocity computing formulaIn γ values be 1.
Transmitting pulse train is illustrated as shown in figure 4, according to ROI window depth locations, signal center frequency size, PRF
The state modulator such as (pulse repeat frequency, pulse recurrence frequency), driving voltage and pulse length launches pulse
Sequence.The head of sequence is one or more detections pulse (Detect beam, high voltage (such as 80V, typically with B-mode phase
Together), short pulse (such as 2 cycles)) as the reference signal of displacement of tissue information；Followed by multigroup driving pulse/detection
Pulse pair, driving pulse (Push beam, low voltage (such as 20～40V), long pulse (100～250 cycles)) is used for
The regional area excitation of near focal point in tissue produces micro-displacement.
Detection pulse (high-tension short pulse) is used to track the deformation that acoustic radiation force load is organized during loading
Situation；It is afterwards the deformation situation organized after a series of detection pulse tracking acoustic radiation force load disappears.The work of driving pulse
With being the local organization at focal position is vibrated under safety condition, long echo signal is due to its spatial resolution
Difference cannot be used for offset estimation, so probe array element launch driving pulse when shutdown signal apodization, all activity array elements send out
Penetrate that signal amplitude is identical, strengthen energy that signal carries to produce bigger acoustic radiation force under safety condition.For encouraging
Long pulse wave frequency rate, and for detection short pulse ripple centre frequency may be selected it is identical or different, it is preferable that centre frequency
Difference, so can in order to avoid driving pulse to detect pulse echo signal interference, while being easy to area in signal transacting
The echo-signal that shunt excitation is encouraged and detected.
In addition the long pulse wave number mesh in transmitting pulse train be it is variable, the wave number mesh of whole sequence be also it is variable,
Here the wave number mesh for once encouraging whole sequence is defined as sampling volume number (ensemblesize) by us, and such as 16,24,
32.Long pulse wave number mesh (field wave)+short pulse wave number mesh (reference and detection ripple)=sampling volume number, gradient is slided and is calculated
The size of window N is the arbitrary value between 0～sampling volume number, and specifically it is to be automatically selected using adaptive algorithm, herein not
Repeat again.
One group of shearing wave transmitting pulse includes several different spatials transmitting pulse train as shown in Figure 4, wherein
Driving pulse locus is fixed, and the detection pulse interval position of different transmitting sequences is different；As shown in figure 5, D is detection
Ripple, P is field wave, if there is one group of shearing wave transmitting pulse using 5 transmitting pulse trains at intervals of 1mm, this 5 sequences
Field wave P1, P2, P3, P4, P5 in the same space position, and detection ripple D1, D2, D3 in 5 transmitting pulse trains,
D4, D5 in identical or different locus, meanwhile, at least one detection pulse locus it is identical with pulse location,
D1 such as in figure.
The transmitting pulse train that ARFI impulse generators are used is identical with shearing pulse generator, but a frame ARFI is imaged
It is the same space position that field wave is kept with detection ripple in scanning process.
Claims (7)
1. a kind of quantitative shearing wave elastograph imaging method, it is characterised in that comprise the steps of：
ARFI detections are carried out, each position strain stress in target area is obtained, wherein comprising specified reference position strain stressref；
Shearing wave detection is carried out, is obtained and specify in target area reference position single-point shear-wave velocity cref；
Each position shear-wave velocity c is calculated according to the single-point shear-wave velocity；
The quantitative shearing wave elastic image E in each position is obtained according to each position strain gauge shear-wave velocity.
2. imaging method as claimed in claim 1, it is characterised in that
Each position strain in the target areaWherein, N is gradient slip calculation window size,
ziThe concrete depth coordinate of each position is represented,Represent the average of each position coordinates in gradient slip calculation window, yiRepresent ladder
Current point displacement in degree slip calculation window,Represent the average of each position displacement in gradient slip calculation window.
3. imaging method as claimed in claim 1, it is characterised in that the single-point shear-wave velocity is according to specified measurement point
Lateral separation, axial distance, and each specified measurement point reached time required for maximum transversal displacement and intended using two-dimensional linear
It is legal to calculate.
4. imaging method as claimed in claim 3, it is characterised in that the single-point shear-wave velocity crefBy formula X=A β
Draw, wherein,In formula, xiRepresent the lateral separation of specified measurement point, ti
The time required to representing that specified measurement point reaches maximum displacement, ziRepresent the axial distance of specified measurement point, β1As required specifies
Measurement point single-point shear wave velocity cref。
5. imaging method as claimed in claim 1, it is characterised in that described each position shear-wave velocity
In focal zone γ values 1；Axial region outside focal zone,Wherein, z is axial distance, and n and σ are
The system constants of excitation sound field are represented, span is the value of two parameters of 0～10, n and σ can be with identical or different.
6. imaging method as claimed in claim 5, it is characterised in that the focal zone DOF=8 (f#)2λ, wherein, λ is represented
Wavelength, f# is aperture control parameter, and it is the real number between 0～5.
7. imaging method as claimed in claim 1, it is characterised in that the quantitative shearing wave elastic image E=3 ρ c2, wherein,
ρ is tested region Media density.
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610973110.5A CN106618639B (en) | 2016-11-04 | 2016-11-04 | A kind of quantitative shearing wave elastograph imaging method |
| PCT/CN2017/107120 WO2018082458A1 (en) | 2016-11-04 | 2017-10-20 | Quantitative shear wave elasticity imaging method and system |
| US16/346,079 US11357480B2 (en) | 2016-11-04 | 2017-10-20 | Quantitative shear wave elasticity imaging method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610973110.5A CN106618639B (en) | 2016-11-04 | 2016-11-04 | A kind of quantitative shearing wave elastograph imaging method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106618639A true CN106618639A (en) | 2017-05-10 |
| CN106618639B CN106618639B (en) | 2019-01-04 |
Family
ID=58821960
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610973110.5A Active CN106618639B (en) | 2016-11-04 | 2016-11-04 | A kind of quantitative shearing wave elastograph imaging method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106618639B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018082458A1 (en) * | 2016-11-04 | 2018-05-11 | 声泰特（成都）科技有限公司 | Quantitative shear wave elasticity imaging method and system |
| CN108670303A (en) * | 2018-02-26 | 2018-10-19 | 长庚大学 | Method and system for detecting uniformity of ultrasonic image |
| CN109875607A (en) * | 2019-01-29 | 2019-06-14 | 中国科学院苏州生物医学工程技术研究所 | Infiltrate tissue testing method, apparatus and system |
| CN110927729A (en) * | 2019-11-09 | 2020-03-27 | 天津大学 | Acoustic radiation force pulse elastography method based on displacement attenuation characteristics |
| CN114245726A (en) * | 2020-04-26 | 2022-03-25 | 深圳迈瑞生物医疗电子股份有限公司 | Prostate elasticity measuring method and ultrasonic imaging system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130261452A1 (en) * | 2012-03-30 | 2013-10-03 | Tadashi Tamura | Methods and apparatus for ultrasound 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 |
| CN103720489A (en) * | 2013-12-30 | 2014-04-16 | 中国科学院深圳先进技术研究院 | Lesion tissue growth monitoring method and system |
| CN104055541A (en) * | 2014-06-26 | 2014-09-24 | 中国科学院苏州生物医学工程技术研究所 | Method for intravascular ultrasound multi-slice shear wave elastography |
| CN104306026A (en) * | 2014-11-18 | 2015-01-28 | 声泰特（成都）科技有限公司 | Echo displacement detection system and imaging system based on acoustic radiation force |
| CN104546014A (en) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Shear wave velocity estimation method for biological tissue elasticity measurement |
-
2016
- 2016-11-04 CN CN201610973110.5A patent/CN106618639B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130261452A1 (en) * | 2012-03-30 | 2013-10-03 | Tadashi Tamura | Methods and apparatus for ultrasound 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 |
| CN103720489A (en) * | 2013-12-30 | 2014-04-16 | 中国科学院深圳先进技术研究院 | Lesion tissue growth monitoring method and system |
| CN104055541A (en) * | 2014-06-26 | 2014-09-24 | 中国科学院苏州生物医学工程技术研究所 | Method for intravascular ultrasound multi-slice shear wave elastography |
| CN104306026A (en) * | 2014-11-18 | 2015-01-28 | 声泰特（成都）科技有限公司 | Echo displacement detection system and imaging system based on acoustic radiation force |
| CN104546014A (en) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Shear wave velocity estimation method for biological tissue elasticity measurement |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018082458A1 (en) * | 2016-11-04 | 2018-05-11 | 声泰特（成都）科技有限公司 | Quantitative shear wave elasticity imaging method and system |
| US11357480B2 (en) | 2016-11-04 | 2022-06-14 | Saset Chengdu Technology Ltd. | Quantitative shear wave elasticity imaging method and system |
| CN108670303A (en) * | 2018-02-26 | 2018-10-19 | 长庚大学 | Method and system for detecting uniformity of ultrasonic image |
| CN108670303B (en) * | 2018-02-26 | 2021-06-25 | 长庚大学 | Method and system for detecting uniformity of ultrasonic image |
| CN109875607A (en) * | 2019-01-29 | 2019-06-14 | 中国科学院苏州生物医学工程技术研究所 | Infiltrate tissue testing method, apparatus and system |
| CN110927729A (en) * | 2019-11-09 | 2020-03-27 | 天津大学 | Acoustic radiation force pulse elastography method based on displacement attenuation characteristics |
| CN110927729B (en) * | 2019-11-09 | 2022-04-01 | 天津大学 | Acoustic radiation force pulse elastography method based on displacement attenuation characteristics |
| CN114245726A (en) * | 2020-04-26 | 2022-03-25 | 深圳迈瑞生物医疗电子股份有限公司 | Prostate elasticity measuring method and ultrasonic imaging system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106618639B (en) | 2019-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106618638B (en) | A kind of quantitative shearing wave elastogram system | |
| AU2018203912B2 (en) | Coherent spread-spectrum coded waveforms in synthetic aperture image formation | |
| JP6000569B2 (en) | Ultrasonic diagnostic apparatus and control program | |
| Jensen | Medical ultrasound imaging | |
| CN106618639A (en) | Quantitative shear wave elastography method | |
| JP5530685B2 (en) | System and method for detecting areas of varying stiffness | |
| CN103505243B (en) | Measure hyperacoustic sound absorption or decay | |
| US20130296743A1 (en) | Ultrasound for Therapy Control or Monitoring | |
| EP2536338B1 (en) | Subject information processing apparatus using acoustic waves received from the subject | |
| EP2578162B1 (en) | Ultrasound diagnosis device | |
| CN105167802B (en) | Doppler imaging method and device | |
| WO2018082458A1 (en) | Quantitative shear wave elasticity imaging method and system | |
| JP6559808B2 (en) | Ultrasonic system and method of operating an ultrasonic system | |
| EP2962642A1 (en) | Method for detecting cavitation and ultrasonic medical apparatus therefor | |
| CN110301936B (en) | Frequency scanning for acoustic radiation force pulses | |
| WO2020113397A1 (en) | Ultrasonic imaging method and ultrasonic imaging system | |
| Hemmsen et al. | Tissue harmonic synthetic aperture ultrasound imaging | |
| JP2004181209A (en) | Ultrasonic diagnostic apparatus | |
| KR20100000881A (en) | Ultrasound system and method for forming an elastic image | |
| KR101857346B1 (en) | Portable ultrasonic diagnostic apparatus and operating method thereof | |
| JP4599208B2 (en) | Ultrasonic diagnostic equipment | |
| KR101117547B1 (en) | Ultrasound Diagnostic System and Method For Forming Elasticity Image Using Plane Wave | |
| JP2005245479A (en) | Ultrasonic diagnostic apparatus | |
| CN111281423A (en) | Ultrasonic image optimization method and ultrasonic imaging equipment | |
| JP7280713B2 (en) | ultrasound diagnostic equipment |
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 |