US8659975B2 - Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions - Google Patents
Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions Download PDFInfo
- Publication number
- US8659975B2 US8659975B2 US12/992,228 US99222809A US8659975B2 US 8659975 B2 US8659975 B2 US 8659975B2 US 99222809 A US99222809 A US 99222809A US 8659975 B2 US8659975 B2 US 8659975B2
- Authority
- US
- United States
- Prior art keywords
- motion
- ultrasonic
- subject
- recited
- echo signals
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active, expires
Links
Images
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
Definitions
- the field of the invention is coherent imaging using vibratory energy, such as ultrasound and, in particular, systems and methods for shearwave dispersion ultrasound vibrometry (“SDUV”).
- vibratory energy such as ultrasound and, in particular, systems and methods for shearwave dispersion ultrasound vibrometry (“SDUV”).
- SDUV shearwave dispersion ultrasound vibrometry
- the ultrasound transmitter may be placed on one side of the object and the sound transmitted through the object to the ultrasound receiver placed on the other side.
- transmission mode methods an image may be produced in which the brightness of each pixel is a function of the amplitude of the ultrasound that reaches the receiver (“attenuation” mode), or the brightness of each pixel is a function of the time required for the sound to reach the receiver (“time-of-flight” or “speed of sound” mode).
- the receiver may be positioned on the same side of the object as the transmitter and an image may be produced in which the brightness of each pixel is a function of the amplitude or time-of-flight of the ultrasound reflected from the object back to the receiver (“refraction,” “backscatter,” or “echo” mode).
- A-scan an ultrasound pulse is directed into the object by an ultrasound transducer and the amplitude of the reflected sound is recorded over a period of time.
- the amplitude of the echo signal is proportional to the scattering strength of the refractors in the object and the time delay is proportional to the range of the refractors from the transducer.
- B-scan the transducer transmits a series of ultrasonic pulses as it is scanned across the object along a single axis of motion.
- the resulting echo signals are recorded as with the A-scan method and their amplitude is used to modulate the brightness of pixels on a display.
- the location of the transducer and the time delay of the received echo signals locates the pixels to be illuminated. With the B-scan method, enough data are acquired from which a two-dimensional image of the refractors can be reconstructed. Rather than physically moving the transducer over the subject to perform a scan it is more common to employ an array of transducer elements and electronically move an ultrasonic beam over a region in the subject.
- the transducer typically has a number of piezoelectric elements arranged in an array and driven with separate voltages (“apodizing”). By controlling the time delay (or phase) and amplitude of the applied voltages, the ultrasonic waves produced by the piezoelectric elements (“transmission mode”) combine to produce a net ultrasonic wave focused at a selected point. By controlling the time delay and amplitude of the applied voltages, this focal point can be moved in a plane to scan the subject.
- the transducer is employed to receive the reflected sound (“receiver mode”). That is, the voltages produced at the transducer elements in the array are summed together such that the net signal is indicative of the sound reflected from a single focal point in the subject. As with the transmission mode, this focused reception of the ultrasonic energy is achieved by imparting separate time delay (and/or phase shifts) and gains to the echo signal received by each transducer array element.
- a linear array system includes a transducer having a large number of elements disposed in a line.
- a small group of elements are energized to produce an ultrasonic beam that travels away from the transducer, perpendicular to its surface.
- the group of energized elements is translated along the length of the transducer during the scan to produce a corresponding series of beams that produce echo signals from a two-dimensional region in the subject.
- the pulsing of the inner elements in each energized group is delayed with respect to the pulsing of the outer elements.
- the time delays determine the depth of focus which can be changed during scanning. The same delay factors are applied when receiving the echo signals to provide dynamic focusing during the receive mode.
- the second common form of ultrasonic imaging is referred to as phased array sector scanning (“PASS”).
- PASS phased array sector scanning
- Such a scan is comprised of a series of measurements in which all of the elements of a transducer array are used to transmit a steered ultrasonic beam.
- the system then switches to receive mode after a short time interval, and the reflected ultrasonic wave is received by all of the transducer elements.
- the transmission and reception are steered in the same direction, ⁇ , during each measurement to acquire data from a series of points along a scan line.
- the receiver is dynamically focused at a succession of ranges, R, along the scan line as the reflected ultrasonic waves are received.
- a series of measurements are made at successive steering angles, ⁇ , to scan a pie-shaped sector of the subject.
- the time required to conduct the entire scan is a function of the time required to make each measurement and the number of measurements required to cover the entire region of interest at the desired resolution and signal-to-noise ratio. For example, a total of 128 scan lines may be acquired over a sector spanning 90 degrees, with each scan line being steered in increments of 0.70 degrees.
- the same scanning methods may be used to acquire a three-dimensional image of the subject.
- the transducer in such case is a two-dimensional array of elements which steer a beam throughout a volume of interest or linearly scan a plurality of adjacent two-dimensional slices.
- a focused ultrasound beam operating within FDA safety limits, is applied to a subject to generate harmonic shear waves in a tissue of interest.
- the propagation speed of the induced shear wave is frequency dependent, or “dispersive,” and relates to the mechanical properties of the tissue of interest.
- Shear wave speeds at a number of frequencies are measured by pulse echo ultrasound and subsequently fit with a theoretical dispersion model to inversely solve for tissue elasticity and viscosity. These shear wave speeds are estimated from the phase of tissue vibration that is detected between two or more points with known distance along the shear wave propagation path.
- One feature of the SDUV method is the use of a so-called “binary pushing pulse” that allows the operation of one single array ultrasound transducer for both motion excitation and the echo signal detection.
- the transducer focuses ultrasound at one location, the “vibration origin,” to vibrate the tissue of interest and then electronically steers its focus to another location, a “motion detection point,” for echo signal vibration detection.
- the “pushing” ultrasound is turned on during a vibration time period to vibrate the tissue and turned off to provide a time window for the pulse echo motion detection.
- the pushing pulse is off, a series of short ultrasound pulses is transmitted to the motion detection location and a corresponding series of echo signals is received and processed to determine the tissue vibration.
- This intermittent pulse sequencing strategy allows both the production of a shear wave and the monitoring of its propagation at the same time with a single array transducer.
- a technical challenge for the SDUV method is that the shear wave generated by the pushing ultrasound is small and difficult to detect with pulse echo ultrasound.
- a Kalman filtering process is often applied to the vibration time-sequence, which is demodulated from the series of ultrasound echoes received from the same motion detection point, to reliably extract the phase of tissue vibration.
- gross tissue motion due to, for example, respiratory and cardiac activity is large and confounds the estimation of phases. With larger shear wave magnitudes the effects of this background motion can potentially be mitigated.
- Larger shear waves can be produced in the subject by increasing the intensity of the pushing pulse; however, current ultrasound scanners may not be able to provide an ultrasound beam intensity high enough to meet the requirement of this solution.
- the intensity of diagnostic ultrasound is limited by the FDA, typically rendering an increase in pushing pulse intensity unfeasible.
- SDUV shearwave dispersion ultrasound vibrometry
- the present invention overcomes the aforementioned drawbacks by providing a method for shearwave dispersion ultrasound vibrometry (“SDUV”) that compensates for large background motions.
- SDUV shearwave dispersion ultrasound vibrometry
- harmonic shear wave motion with a greater magnitude is produced through the application of an ultrasonic vibration pulse having a long duration.
- the ultrasonic vibration pulse is applied during an “on” time period having a duration, ⁇ T.
- the longer duration vibration pulse is accompanied by a series of ultrasonic detection pulses having an increased pulse repetition frequency (“PRF”) and correspondingly decreased pulse repetition period (“PRP”).
- PRF pulse repetition frequency
- PRP pulse repetition period
- Kalman filtering is applied to the echo signals received from the motion detection point to detect a signal indicative of the imparted harmonic shear wave motion.
- the Kalman filter specifically accounts for the non-uniform temporal sampling resulting from the “on” time period duration being larger than the PRP of the detection pulses; that is, ⁇ T>PRP.
- This provides a method for Kalman filtering that does not first require applying a band-pass filter to the echo signals. In turn, this allows for the application of the non-uniform temporal pattern where ⁇ T>PRP.
- the multiple frequency Kalman filter adaptation decreases errors that arise during the estimation of the phase and amplitude of the harmonic motion by properly accounting for higher harmonics in that motion.
- additional state variables include the displacement and velocity of the background tissue motion.
- the transition matrix in the state equation is correspondingly modified to account for these additional variables, as is the measurement vector.
- FIG. 1A is a pulse timing diagram indicating the application of ultrasonic vibration and detection pulses in a substantially uniform pattern
- FIG. 1B is a pulse timing diagram indicating the application of ultrasonic vibration and detection pulses in a non-uniform pattern
- FIG. 2 is a block diagram of a vibroacoustography system that employs the present invention
- FIG. 3 is a block diagram of an ultrasound system used in the system of FIG. 2 ;
- FIG. 4 is a block diagram of a linear ultrasound transducer that forms a part of one configuration of the ultrasound system of FIG. 3 ;
- FIG. 5 is a block diagram of a transmitter which forms part of the ultrasound system of FIG. 3 ;
- FIG. 6 is a block diagram of a receiver which forms part of the ultrasound system of FIG. 3 ;
- FIG. 7 is a flow chart setting forth the steps of a method for measuring the mechanical properties of a tissue of interest in accordance with the present invention.
- SDUV shearwave dispersion ultrasound vibrometry methods
- an ultrasonic imaging system that interrogates the tissue of interest with a pulsed ultrasound beam and examines the resulting echo signals to measure the phase and amplitude of the harmonic motion imparted to the subject tissue.
- the challenge is to extract this information from the echo signals where the amplitude of the harmonic motion is at the submicron level.
- two general adaptations are provided.
- a pushing pulse with a longer duration is applied to induce shear waves of larger intensity, making them easier to differentiate from background motion.
- the background motion can be compensated through modeled, thereby reducing its influence on the estimated phase and amplitude of the desired harmonic shear wave motion.
- D ( t ) D 0 sin( ⁇ s t+ ⁇ s ) Eqn. (1);
- A(t f ,t s ) is the amplitude of the complex envelope of the received echo signal; t f is fast time, which is representative of depth; t s is slow time, which is representative of pulse repetition; ⁇ f is the transmission center frequency; ⁇ s is again the tissue vibration frequency; ⁇ s is again the tissue vibration phase; and ⁇ is an angle between the ultrasound beam and direction of tissue motion.
- the vibration becomes a phase term of the angular modulation.
- the modulation index, ⁇ is:
- ⁇ 0 is a constant phase added during the quadrature demodulation to keep I(t f ,t s ) nonzero; that is, to be either all positive, or all negative, in slow time at a location.
- this constant phase term is removed, as shown below.
- the arctangent of the ratio between Q(t f ,t s ) and I(t f ,t s ) is calculated accordingly:
- s (t f ,t s ) is a mean value of s(t f ,t s ) in slow time. If the sampling frequency in fast time is high, I(t f ,t s ) and Q(t f ,t s ) can be averaged with a limited length in fast time to reduce noise before s(t f ,t s ) is calculated.
- a band-pass filter (“BPF”) centered at the vibration frequency can improve y(t s ) by reducing noise and distortions.
- ⁇ y is a standard deviation of y(t s ).
- the amplitude of the oscillatory Doppler shifts can also be directly measured by applying a turbulence estimation method to the received echo signals, r(t f ,t s ), to estimate the variance of motion velocity.
- a Kalman filter process is employed to recursively estimate the phase and amplitude of the harmonic shear wave motion.
- a Kalman filter is a numerical method used to track a time-varying signal in the presence of noise. If the signal can be characterized by some number of parameters that vary slowly with time, then Kalman filtering can be used to tell how incoming raw measurements should be processed to best estimate those parameters as a function of time.
- a Kalman filter extracts information about the imparted harmonic motion from random and noisy measurement data with known vibration frequency and unknown vibration amplitude and phase.
- an ultrasonic vibration pulse 100 is applied to a vibration origin 200 in a tissue of interest 202 in order to induce harmonic shear wave motion 204 in the tissue 202 .
- This shear wave motion 204 at a location in the tissue is interrogated by applying a series of ultrasonic detection pulses 102 to a motion detection point 206 .
- Tissue vibration obtained from the k th ultrasound echo has the form:
- y k [ sin ⁇ ( ⁇ s ⁇ kT ) cos ⁇ ( ⁇ s ⁇ kT ) ] ⁇ [ ⁇ ⁇ ⁇ cos ⁇ ( ⁇ s ) ⁇ ⁇ ⁇ sin ⁇ ( ⁇ s ) ] + n k ; Eqn . ⁇ ( 14 )
- ⁇ s is the frequency of the tissue vibration
- T is the pulse repetition period (“PRP”), which is equal to the inverse of the pulse repetition frequency (“PRF”), or 1/PRF
- PRP pulse repetition period
- PRF pulse repetition frequency
- ⁇ and ⁇ s are the vibration amplitude and phase to be estimated, respectively
- x k is a state variable having the form:
- w k is a white driving sequence vector that allows some variations in the vibration amplitude and phase
- ⁇ is a transition matrix, which has the form:
- Eqns. (15)-(18) are updated every T seconds in the Kalman filter.
- the Kalman filter is recursively applied to extract, from the noisy echo signals, the amplitude and phase of the underlying signal indicative of the harmonic shear wave motion.
- the amplitude and phase parameters have the form:
- the pushing pulse duration is increased while maintaining its intensity. This provides an increase in the induced shear wave magnitude while maintaining ultrasound beam intensities within safe limits.
- a lower pulse repetition frequency (“PRF”) is required to maintain a uniform time interval between detection pulse periods. This can result in inaccurate measurements in situations where large background tissue motion occurs, such as, when examining the heart.
- PRF is increased during the detection pulse period. However, this results in non-uniform temporal sampling illustrated in the pulse timing sequence shown in FIG. 1B .
- the ultrasonic vibration pulse 100 is applied during a time period of duration ⁇ T.
- the PRF of the ultrasonic detection pulses 102 is increased, the pulse repetition period (“PRP”) is correspondingly decreased.
- PRP pulse repetition period
- a non-uniform temporal sampling is generally defined for those situations when ⁇ T>PRP, as illustrated in FIG. 1B .
- a non-uniform temporal sampling also exists when ⁇ T ⁇ PRP; however, this will result in aliasing and is a generally undesirable configuration.
- the pulse timing sequence shown in FIG. 1A where a uniform and lower PRF is employed, is very simple and may work well under certain conditions.
- interpolation can be used on the vibration sequence demodulated from ultrasound echoes.
- the interpolated sequence will have uniform temporal sampling and, therefore, can be processed as if there are no missing data samples. While this approach works well for in vitro tissue experiments, the results will degrade if large background motion is present. Thus, this approach is likely limited for in vivo applications, which are likely to see large background motions. Also, as stated above, if the PRF is too low, the higher harmonics of the pushing pulses will exceed the Nyquist sampling limit of the motion detection ultrasound, which will result in aliasing artifacts.
- the non-uniform temporal sampling approach shown in FIG. 1B provides a degree of flexibility in pulse sequence design.
- the pushing pulse 100 can be significantly longer than in uniform approaches, such as the one in FIG. 1A .
- a higher overall PRF can also be maintained, which provides some compensation when the signal-to-noise ratio (“SNR”) is low for the detected signal indicative of the harmonic shear wave motion.
- SNR signal-to-noise ratio
- the non-uniform temporal sampling shown in FIG. 1B presents a processing challenge that is addresses during the Kalman filtering process.
- Kalman filtering can be applied directly to non-uniform sampled data, if the measurement equation properly accommodates the differences between the pushing pulse period and the pulse repetition period (“PRP”) of the detection pulses. Therefore, the state equation given by Eqn. (18) is applicable if modified, even when tissue motion is sampled non-uniformly in time. To achieve this, the measurement vector, H k , is modified to properly account for the exact timing of each iteration, thereby providing a means by which the Kalman filter can accommodate non-uniform sampled sequence. With reference to the sampling pattern provided in FIG.
- H k [sin( ⁇ s ( kT+n ( ⁇ T ⁇ T )))cos( ⁇ s ( kT+n ( ⁇ T ⁇ T ))] Eqn. (23).
- the binary pushing pulses when performing SDUV the binary pushing pulses generate harmonics in addition to the fundamental pushing frequency, and these harmonics can introduce errors during the Kalman filtering process if they are not included in the vibration model utilized by the Kalman filter. This is because the Kalman filter assumes that the signal indicative of the harmonic motion contains only one frequency and that all other out-of-band motions are white noise. While this assumption is correct under certain circumstances, such as where one transducer is dedicated to maintain continuous monochromatic vibration within the studied tissue, it is not true when the binary pushing pulses contain a fundamental frequency (e.g., 100 Hz) and higher harmonics (e.g., 200 Hz, 300 Hz, and 400 Hz).
- a fundamental frequency e.g. 100 Hz
- harmonics e.g. 200 Hz, 300 Hz, and 400 Hz
- H k [s 1 c 1 s 2 c 2 . . . s N c N ] Eqn. (24);
- s m sin(m ⁇ s kT)
- c m cos(m ⁇ s kT)
- m 1, 2, 3, . . . , N.
- s m sin( m ⁇ s ( kT+n ( ⁇ T ⁇ T ))
- c m cos( m ⁇ s ( kT+n ( ⁇ T ⁇ T )) Eqn. (26).
- the state variable vector, x k , measurement vector, H k , and the Kalman filter gain vector all have 2N elements.
- the transition matrix, ⁇ , covariance matrix of the estimation errors, and the covariance matrix of the driving sequences all have a dimension of 2N ⁇ 2N, which are extensions of the 2 ⁇ 2 matrices described above. While only the first few harmonics are generally useful when estimating the harmonic shear wave motion, all of the harmonics are included in the Kalman filtering process for more accurate estimations.
- small background motions can be accommodated by the white driving sequence vector, w k , in Eqn. (18), which allows limited variations in the vibration amplitude and phase over time.
- the variances of, w k are given by the Q k covariance matrix of the Kalman filter.
- a random process can be used to model this motion, such as, for example, a random walk model.
- two additional state variables are added to the state equation, Eqn. (18).
- D B (kT) is added to model the background displacement and another, ⁇ B (kT), to model the background velocity.
- x k [ 1 0 0 0 0 1 0 0 0 0 1 T 0 0 0 1 ] ⁇ x k - 1 + w k ; Eqn . ⁇ ( 28 )
- H k [sin( ⁇ s kT )cos( ⁇ s kT )1 0] Eqn. (30);
- H k [sin( ⁇ s ( kT+n ( ⁇ T ⁇ T )))cos( ⁇ s ( kT+n ( ⁇ T ⁇ T )))1 0] Eqn. (31).
- a vibroacoustography system which employs the present invention employs an ultrasonic transducer 311 that is operable to produce focused ultrasound beams.
- the transducer 311 such as a phased array transducer, intermittently transmits a beam of ultrasonic vibration pulses 100 to a vibration origin 200 in the tissue of interest 202 to vibrate, or oscillate, the tissue 202 at a prescribed frequency.
- the focus of the transducer is electronically steered to a motion detection point 206 at a distance, ⁇ r, from the vibration origin 200 and harmonic motion 204 at that point is detected.
- a vibration mode is multiplexed with a detection mode. This enables the detection of the harmonic motion 206 by the same transducer 311 as that transmitting the vibration pulses 100 and both vibration and detection can be achieved without mechanically moving the transducer 311 .
- a transducer array 311 includes a plurality of separately driven elements 312 which each produce a burst of ultrasonic energy when energized by a pulse produced by a transmitter 313 .
- the ultrasonic energy reflected back to the transducer array 311 from the subject under study is converted to an electrical signal by each transducer element 312 and applied separately to a receiver 314 through a set of switches 315 .
- the transmitter 313 , receiver 314 , and the switches 315 are operated under the control of a digital controller 316 responsive to the commands input by the human operator.
- a complete scan is performed by acquiring a series of echoes in which the switches 315 are set to their transmit position, the transmitter 313 is gated on momentarily to energize each transducer element 312 , the switches 315 are then set to their receive position, and the subsequent echo signals produced by each transducer element 312 are applied to the receiver 314 .
- the separate echo signals from each transducer element 312 are combined in the receiver 314 to produce a single echo signal which is employed to produce a line in an image on a display system 317 .
- the transmitter drives the transducer array 311 such that the ultrasonic energy produced is directed, or steered, in a beam.
- a B-scan can therefore be performed by moving this beam through a set of angles from point-to-point rather than physically moving the transducer array 311 .
- a sector scan is performed by progressively changing the time delays T i in successive excitations.
- the angle ⁇ is thus changed in increments to steer the transmitted beam in a succession of directions.
- the timing of the pulses 320 is reversed.
- a linear transducer array may also be employed in the ultrasonic imaging system of FIG. 3 .
- the transmitter 313 drives the transducer array 311 such that an ultrasonic beam is produced which is directed substantially perpendicular to its front surface.
- a subgroup of the elements 312 are energized to produce the beam and the pulsing of the inner elements 312 in this subgroup are delayed relative to the outer elements 312 as shown at 320 .
- a beam focused at point P results from the interference of the small separate wavelets produced by the subgroup elements.
- the time delays determine the depth of focus, or range R, and this is typically changed during a scan when a two-dimensional image is to be produced.
- the same time delay pattern is used when receiving the echo signals resulting in dynamic focusing of the echo signals received by the subgroup of elements 312 . In this manner a single scan line in the image is formed.
- the subgroup of elements to be energized are shifted one element position along the transducer length and another scan line is required.
- the focal point, P, of the ultrasonic beam is thus shifted along the length of the transducer 311 by repeatedly shifting the location of the energized subgroup of elements 312 .
- the echo signals produced by each burst of ultrasonic energy emanate from reflecting objects located at successive positions, R, along the ultrasonic beam. These are sensed separately by each segment 312 of the transducer array 311 and a sample of the magnitude of the echo signal at a particular point in time represents the amount of reflection occurring at a specific range, R. Due to the differences in the propagation paths between a focal point, P, and each transducer element 312 , however, these echo signals will not occur simultaneously and their amplitudes will not be equal.
- the function of the receiver 314 is to amplify and demodulate these separate echo signals, impart the proper time delay to each and sum them together to provide a single echo signal which accurately indicates the total ultrasonic energy reflected from each focal point, P, located at range, R, along the ultrasonic beam oriented at the angle, ⁇ .
- time delays are introduced into each separate transducer element channel of the receiver 314 .
- the delay introduced in each channel may be divided into two components, one component is referred to as a beam steering time delay, and the other component is referred to as a beam focusing time delay.
- the beam steering and beam focusing time delays for reception are precisely the same delays, T i , as the transmission delays described above.
- the focusing time delay component introduced into each receiver channel is continuously changing during reception of the echo to provide dynamic focusing of the received beam at the range, R, from which the echo signal emanates.
- the receiver 314 Under the direction of the digital controller 316 , the receiver 314 provides delays during the scan such that the steering of the receiver 314 tracks with the direction of the beam steered by the transmitter 311 and it samples the echo signals at a succession of ranges and provides the proper delays to dynamically focus at points, P, along the beam.
- each emission of an ultrasonic pulse results in the acquisition of a series of data points which represent the amount of reflected sound from a corresponding series of points, P, located along the ultrasonic beam.
- echoes from multiple focused locations can be received to measure vibration information from several points of the tissue.
- the limitation of the lateral resolution of the transducer for two closely located points can be improved by assigning different transmitting codes for different locations.
- the display system 317 receives the series of data points produced by the receiver 314 and converts the data to a form producing the desired image. For example, if an A-scan is desired, the magnitude of the series of data points is merely graphed as a function of time. If a B-scan is desired, each data point in the series is used to control the brightness of a pixel in the image, and a scan comprised of a series of measurements at successive locations along the length of the transducer 311 (linear array mode) or steering angles (PASS mode) is performed to provide the data necessary for display of an image.
- linear array mode linear array mode
- PASS mode steering angles
- the transmitter 313 includes a set of channel pulse code memories which are indicated collectively at 550 .
- Each pulse code memory 550 stores a bit pattern 551 that determines the frequency of the ultrasonic pulse 552 that is to be produced. This bit pattern is read out of each pulse code memory 550 by a master clock and applied to a driver 553 which amplifies the signal to a power level suitable for driving the transducer 311 .
- the bit pattern is a sequence of four “1” bits alternated with four “0” bits to produce a 5 megahertz (“MHz”) ultrasonic pulse 552 .
- the transducer elements 311 to which these ultrasonic pulses 552 are applied respond by producing ultrasonic energy.
- the pulses 552 for each of the N channels must be produced and delayed by the proper amount. These delays are provided by a transmit control 554 which receives control signals from the digital controller 316 . When the control signal is received, the transmit control 554 gates a clock signal through to the first transmit channel 550 . At each successive delay time interval thereafter, the clock signal is gated through to the next channel pulse code memory 550 until all the channels to be energized are producing their ultrasonic pulses 552 . Each transmit channel 550 is reset after its entire bit pattern 551 has been transmitted and the transmitter 313 then waits for the next control signal from the digital controller 316 .
- the receiver 314 is comprised of three sections: a time-gain control (“TGC”) section 600 , a beam forming section 601 , and a mid processor 602 .
- the time-gain control section 600 includes an amplifier 605 for each of the N receiver channels and a time-gain control circuit 606 .
- the input of each amplifier 605 is connected to a respective one of the transducer elements 312 to receive and amplify the echo signal which it receives.
- the amount of amplification provided by the amplifiers 605 is controlled through a control line 607 that is driven by the time-gain control circuit 606 .
- As is well known in the art as the range of the echo signal increases, its amplitude is diminished.
- the brightness of the image diminishes rapidly as a function of range, R.
- This amplification is controlled by the operator who manually sets TGC linear potentiometers 608 to values which provide a relatively uniform brightness over the entire range of the scan.
- the time interval over which the echo signal is acquired determines the range from which it emanates, and this time interval is divided into segments by the TGC control circuit 606 .
- the settings of the potentiometers are employed to set the gain of the amplifiers 605 during each of the respective time intervals so that the echo signal is amplified in ever increasing amounts over the acquisition time interval.
- the beam forming section 601 of the receiver 314 includes N separate receiver channels 610 .
- Each receiver channel 610 receives the analog echo signal from one of the TGC amplifiers 605 at an input 611 , and it produces a stream of digitized output values on an I bus 612 and a Q bus 613 .
- Each of these I and Q values represents a sample of the echo signal envelope at a specific range, R.
- These samples have been delayed in the manner described above such that when they are summed at summing points 614 and 615 with the I and Q samples from each of the other receiver channels 610 , they indicate the magnitude and phase of the echo signal reflected from a point, P, located at range, R, on the ultrasonic beam.
- the mid processor section 602 receives the beam samples from the summing points 614 and 615 .
- the I and Q values of each beam sample is a digital number which represents the in-phase and quadrature components of the magnitude of the reflected sound from a point, P.
- the resulting magnitude values output at 621 to the display system 317 result in an image in which the magnitude of the reflected echo at each image pixel is indicated by pixel brightness.
- the present invention is implemented by a mechanical property processor 622 which forms part of the mid-processor 602 .
- this processor 602 receives the I and Q beam samples acquired during a sequence of measurements of the subject tissue 202 and calculates a mechanical property of the tissue 202 .
- a method for measuring mechanical properties, such as shear elasticity and shear viscosity, of a tissue of interest using a shear wave speed dispersion technique is performed with an ultrasound system 300 , such as the one described above.
- the mechanical property processor 622 controls the measurements made by the ultrasound system 300 , the vibration pulse output of the transducer elements 312 , processes the resulting echo signals I and Q to satisfy the above-described Eqns. (5) and (6), and calculates a mechanical property of the target tissues.
- target tissues may be, for example, an artery or myocardial tissue, and the mechanical property may be stiffness.
- the ultrasound system 300 is operated first using B-mode scanning to acquire an anatomical ultrasound image of a region of interest, such as the heart or liver, as indicated at step 700 .
- a target of interest is determined such as by selected in the ultrasound image a plurality of motion detection points at which it is desired to obtain mechanical properties, as indicated at step 702 .
- a vibration origin is selected.
- the vibration origin and the plurality of motion detection points are, for example, selected to be co-linear and spaced as discussed below.
- the shear wave speed dispersion technique makes use of an estimation of wave speed using equation (22), which includes a phase difference, such that the phase ⁇ s , of the echo signals is determined at a plurality of points. This is achieved, for example, using a single vibration origin and two different motion detection points each located a different distance from the vibration origin; however, it is also possible to use a single detection point and change the vibration origin.
- the determination of a desired distance, R, between a vibration origin and a motion detection point and the determination of a desired spacing, ⁇ r, between the motion detection points can be based on a consideration of the type of tissue that is under examination. For example, an appropriate distance, R, in the liver may be on the order of one centimeter, while an appropriate distance, R, for a breast lesion may be on the order of five millimeters.
- the outgoing shear wave generated at the vibration origin using amplitude modulated pulses can be approximated as a cylindrical shear wave. Its amplitude decreases as the wave propagates outwards from the excitation center due to both geometry effects and attenuation resulting from the medium in which the wave propagates.
- the gating frequency and the modulation frequency of the ultrasonic pulses are determined, as indicated at step 706 .
- the gating frequency is 2.5 kHz and a prescribed modulation frequency for a first shear wave speed is 100 Hz.
- the transmission of intermittent ultrasonic vibration pulses to a vibration origin, and the transmission of ultrasonic detection pulses and the receipt of echo signals from the motion detection point occurs under the control of the digital controller 316 of the ultrasound system 300 , as indicated at step 708 .
- the ultrasonic vibration pulses are amplitude modulated (“AM”) and are applied to the subject in an on-off time sequence with a detection mode occurring during the off intervals of the vibration pulses.
- AM amplitude modulated
- the non-uniform temporal pulse sequencing illustrated in FIG. 1B is employed, wherein the duration of the “on” period when the pushing pulse 100 is applied, ⁇ T, is larger than the pulse repetition period (“PRP”) of the ultrasonic detection pulses 102 .
- the “on” phase of this sequence essentially comprises an intermittent or gated AM signal composed of ultrasonic vibration pulses that are modulated at a modulation frequency.
- the application of a continuous wave AM signal to a subject in this manner generates a radiation force having a frequency equal to the modulation frequency of the beam.
- the radiation force imparts a harmonic shear wave 204 to the tissue if interest 202 .
- the application of a gated rather than a continuous AM signal generates a radiation force including various frequency components and a shear wave 204 including these frequency components is imparted to the subject tissue 202 .
- two focused ultrasound beams having a beat frequency can also be employed for vibrating the tissue 202 .
- the two beams can be achieved by dividing the elements 312 of the transducer array 311 into two groups and using the first group to transmit ultrasonic pulses at a first frequency, ⁇ 1 , and the second group to transmit ultrasonic pulses at a second frequency, ⁇ 2 .
- the first frequency may be 1 MHz+100 Hz and the second frequency may be 1 MHz ⁇ 100 Hz.
- Such a dual beam will vibrate the tissue 202 at the difference frequency, which in this example is 200 Hz.
- the same concept of intermittent vibration and detection can be used to measure a characteristic of an imparted shear wave in order to characterize a mechanical property of a subject.
- the ultrasound system 300 is operated to acquire echo signals from the subject tissues at a series of motion detection points.
- 100 echoes sampled at a 40 MHz sample rate are acquired at each of 11 motion detection points that are spread evenly along 10 to 20 mm of the length of the artery are measured. Eight echo samples at the desired location are used to obtain average I and Q values.
- the ultrasonic detection pulses can be applied to the motion detection points in a number of ways. For example, each motion detection point can be fully sampled before the detection pulses are steered to the next point. However, in the alternative the plurality of motion detection points can be sampled substantially contemporaneously.
- parallel beamforming is employed it both the transmission of the detection pulses and the reception of the resulting echo signals.
- a substantially plane wave is produced in the region of the tissue of interest undergoing harmonic shear wave motion by properly phasing the transducer elements 312 . This results in a rather broad beam instead of a focused ultrasound beam.
- the backscattered ultrasound is then formed into focused beams in a beamformer after being received by the transducer 311 .
- the transducer elements 312 are energized in subgroups such that a plurality of focused ultrasound beams are directed to the plurality of motion detection points.
- the motion detection points can be sampled in a so-called interleaved manner.
- an ultrasonic detection pulse is applied to the j th motion detection point either before or after an echo is received from the (j ⁇ 1) th motion detection point.
- the time difference between the transmission of an ultrasonic detection pulse to the (j ⁇ 1) th and j th detection point is thereby defined as T d .
- This interleaved sampling pattern can be readily employed using a phased array transducer with proper beamforming, as is known in the art.
- phase estimate (for example, as estimated by Kalman filtering) of the harmonic shear wave motion needs to be corrected accordingly.
- the I values are checked and if a zero crossing occurs, all the I and Q values are reprocessed to add a constant phase, ⁇ 0 , as indicated above in equations (5) and (6). Such a constant phase is added until no zero crossings are detected.
- the amplitude and phase of the tissue motion at each point is then estimated from the acquired I and Q echo samples.
- the arctangent of the ratio of the Q and I beam samples are calculated and the mean value is removed to obtain the harmonic motion in slow time, as indicated above in Equations (7) and (8).
- the harmonic motion is modeled by a second order differential equation with random amplitude and phase and the known beat frequency.
- x k [ 1 0 0 0 0 1 0 0 0 0 1 T 0 0 0 1 ] ⁇ x k - 1 + w k ;
- H k [sin( ⁇ s ( kT+n ( ⁇ T ⁇ T )))cos( ⁇ s ( kT+n ( ⁇ T ⁇ T ))1 0].
- the change in tissue oscillation phase as a function of distance is then calculated for this beat frequency using the calculated phase values at the 11 points along the artery.
- the digital controller determines whether the last frequency has been measured. If not, at process block 718 , another frequency is selected and process blocks 708 - 714 are repeated at each desired prescribed frequency. For example, shear wave speeds are calculated using a set of modulation frequencies including 100 Hz, 200 Hz, 300 Hz, 400 Hz and 500 Hz.
- the next step is to calculate the shear wave speeds in the subject tissue 202 at the different beat frequencies.
- Linear regression is applied to the 11 phase changed measurements to yield a phase change over 10 mm distance along the artery. From this phase change over distance information, the shear wave speed at each beat frequency is estimated as described with reference to equation (22).
- the final step is to calculate a mechanical property of the tissue 202 from the shear wave speed information.
- the shear elasticity and viscosity of the tissue 202 is estimated from the set of shear wave speeds.
- These mechanical properties indicate the stiffness of the artery, which is a valuable clinical measurement, but also, for example, the stiffness of myocardial and liver tissue.
- This calculation is based on shear wave speed dispersion as described, for example, by S. Chen, et al., in “Complex Stiffness Quantification Using Ultrasound Stimulated Vibrometry,” IEEE Ultrasonics Symposium, 2003; 941-944.
- the shear wave speeds at multiple frequencies are fit with appropriate theoretical models to solve for the shear elasticity and viscosity. For example, one appropriate equation is the so-called Voigt model:
- ⁇ 1 is the shear modulus
- ⁇ 2 is the shear viscosity
- ⁇ is frequency
- ⁇ is the density of the tissue, which can be assumed to be 1000 kilograms per cubic meter (kg/m 3 ).
- harmonic motion is detected at a motion detection point by transmitting detection pulses to the motion detection point and receiving echo pulses therefrom. These signals are then analyzed as described above in the mid-processor 602 of the receiver 314 .
- a signal indicative of the induced harmonic shear wave motion is detected at the prescribed frequency in the received ultrasonic echo signals and a characteristic of the detected signal, such as amplitude and/or phase, is determined.
- the mechanical property is then calculated using the measured characteristic.
- an amplitude could be measured and the change in amplitude over distance could be determined and used in conjunction with an appropriate model to determine one or more other mechanical properties of the subject target tissue.
- tone bursts of ultrasonic pulses are not modulated but have a specific amplitude, duration and period such that they impart a force having desired frequency components.
- tone bursts of ultrasonic pulses having a duration of 1 ms and repeated every 10 ms will generate a radiation force that includes frequency components at a fundamental frequency and multiples thereof, for example at 100 Hz, 200 Hz, 300 Hz, 400 Hz and 500 Hz.
- Using such a waveform for vibration provides advantages in faster data acquisition as described below as well as in lower tissue and transducer heating.
- the measurements need not be repeated over different frequencies, since the tone bursts include frequency components at multiple frequencies. These tone bursts are also easier to implement using a conventional ultrasound system and the average intensity and power would be lower as compared to for example a similar gated AM modulation method. Further, because the off interval is generally longer when utilizing tone bursts, the detection pulses can be steered to more than one motion detection point during the off intervals.
- an ultrasound system can be used for excitation in a manner such as described herein, and other known means, such as MRI or optical methods, can be used to detect the ultrasound motion.
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Biophysics (AREA)
- Nonlinear Science (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Radar, Positioning & Navigation (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
D(t)=D 0 sin(ωs t+φ s) Eqn. (1);
r(t f ,t s)=A(t f ,t s)cos(ωf ,t f+φf+β sin(ωf t s+φs)) Eqn. (3);
I(t f ,t s)=A(t f ,t s)cos(β sin(ωs t s+φs)+φf+φ0) Eqn. (5);
and
Q(t f ,t s)=A(t f ,t s)sin(β sin(ωs t s+φs)+φf+φ0) Eqn. (6);
and
y(t s)=s(t f ,t s)−
β=√{square root over (2)}σy Eqn. (9);
I(t s)=β cos(φs) Eqn. (10);
and
Q(t s)=β sin(φs) Eqn. (11).
β(t s)=√{square root over (I(t s)2 +Q(t s)2)}{square root over (I(t s)2 +Q(t s)2)} Eqn. (12);
y k =H k x k +n k Eqn. (15);
H k=[sin(ωs kT)cos(ωs kT)] Eqn. (16);
x k =Φx k-1 +w k Eqn. (18);
H k=[sin(ωs(kT+n(ΔT−T)))cos(ωs(kT+n(ΔT−T)))] Eqn. (23).
Hk=[s1 c1 s2 c2 . . . sN cN] Eqn. (24);
s m=sin(mω s(kT+n(ΔT−T))) Eqn. (25);
and,
c m=cos(mω s(kT+n(ΔT−T))) Eqn. (26).
H k=[sin(ωs kT)cos(ωs kT)1 0] Eqn. (30);
H k=[sin(ωs(kT+n(ΔT−T)))cos(ωs(kT+n(ΔT−T)))1 0] Eqn. (31).
M=√{square root over (I 2 +Q 2)} Eqn. (33).
y k(kT s ,z j)=β sin(ωs kT s+ωs(j−1)T d+φs(z j)) Eqn. (34);
{tilde over (φ)}s(z j)=ωs(j−1)T d+φs(z j) Eqn. (35).
φs(z j)={tilde over (φ)}s(z j)−ωs(j−1)T d Eqn. (36).
P k =E[(x k −{circumflex over (x)} k)(x k −{circumflex over (x)} k)T] Eqn. (37);
H k=[sin(ωs(kT+n(ΔT−T)))cos(ωs(kT+n(ΔT−T)))1 0].
Claims (15)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/992,228 US8659975B2 (en) | 1997-07-21 | 2009-05-15 | Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US5326497P | 1997-07-21 | 1997-07-21 | |
PCT/US2009/044163 WO2009140607A1 (en) | 2008-05-15 | 2009-05-15 | Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions |
US12/992,228 US8659975B2 (en) | 1997-07-21 | 2009-05-15 | Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions |
Publications (2)
Publication Number | Publication Date |
---|---|
US20110063950A1 US20110063950A1 (en) | 2011-03-17 |
US8659975B2 true US8659975B2 (en) | 2014-02-25 |
Family
ID=21983010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/992,228 Active 2030-12-31 US8659975B2 (en) | 1997-07-21 | 2009-05-15 | Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions |
Country Status (2)
Country | Link |
---|---|
US (1) | US8659975B2 (en) |
ZA (1) | ZA985834B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10024757B1 (en) | 2017-05-04 | 2018-07-17 | United Launch Alliance, L.L.C. | Non-uniform sampling in bandwidth constrained data acquisition systems |
US11857373B2 (en) | 2016-10-04 | 2024-01-02 | The Regents Of The University Of California | Multi-frequency harmonic acoustography for target identification and border detection |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8187187B2 (en) * | 2008-07-16 | 2012-05-29 | Siemens Medical Solutions Usa, Inc. | Shear wave imaging |
US8913816B2 (en) * | 2009-04-06 | 2014-12-16 | Hitachi Medical Corporation | Medical image dianostic device, region-of-interest setting method, and medical image processing device |
WO2011007278A2 (en) * | 2009-07-17 | 2011-01-20 | Koninklijke Philips Electronics, N.V. | Spatially-fine shear wave dispersion ultrasound vibrometry sampling |
CA2826770C (en) | 2011-02-15 | 2020-06-23 | Hemosonics, Llc | Characterization of blood hemostasis and oxygen transport parameters |
US8532430B2 (en) * | 2011-07-28 | 2013-09-10 | General Electric Company | Methods for reducing motion artifacts in shear wave images |
CN103300890B (en) | 2012-03-16 | 2016-06-08 | 通用电气公司 | For measuring the system and method for tissue mechanical properties |
US10624609B2 (en) | 2012-10-07 | 2020-04-21 | Mayo Foundation For Medical Education And Research | System and method for shear wave elastography by transmitting ultrasound with subgroups of ultrasound transducer elements |
WO2014201020A1 (en) * | 2013-06-10 | 2014-12-18 | Mayo Foundation For Medical Education And Research | System and method for acoustic radiation force creep-recovery and shear wave propagation for elasticity imaging |
JP5730978B2 (en) | 2013-11-08 | 2015-06-10 | 日立アロカメディカル株式会社 | Ultrasonic diagnostic apparatus and method |
JP6361528B2 (en) * | 2015-03-03 | 2018-07-25 | コニカミノルタ株式会社 | Ultrasonic treatment equipment |
US9726647B2 (en) | 2015-03-17 | 2017-08-08 | Hemosonics, Llc | Determining mechanical properties via ultrasound-induced resonance |
US10444203B2 (en) * | 2016-09-15 | 2019-10-15 | Texas Instruments Incorporated | Ultrasonic vibration sensing |
US11357478B2 (en) | 2017-04-11 | 2022-06-14 | Massachusetts Institute Of Technology | Methods and systems for shear wave elastography |
JP7253540B2 (en) * | 2017-10-12 | 2023-04-06 | コーニンクレッカ フィリップス エヌ ヴェ | Ultrasound Shear Wave Imaging with Patient-Adaptive Shear Wave Generation |
KR102665895B1 (en) * | 2017-10-16 | 2024-05-16 | 삼성메디슨 주식회사 | Ultrasound diagnosis apparatus and operating method for the same |
US11154277B2 (en) * | 2017-10-31 | 2021-10-26 | Siemens Medical Solutions Usa, Inc. | Tissue viscoelastic estimation from shear velocity in ultrasound medical imaging |
US11607199B2 (en) * | 2018-11-20 | 2023-03-21 | Siemens Medical Solutions Usa, Inc. | Switched capacitor for elasticity mode imaging with ultrasound |
WO2020112586A1 (en) * | 2018-11-29 | 2020-06-04 | The Penn State Research Foundation | Harmonic shear wave imaging |
US10611497B1 (en) | 2019-02-18 | 2020-04-07 | Amazon Technologies, Inc. | Determining vehicle integrity using vibrometric signatures |
US10861164B1 (en) * | 2019-03-28 | 2020-12-08 | Amazon Technologies, Inc. | Visually determining vibrometric behavior of aerial vehicles |
US11079303B1 (en) | 2019-06-11 | 2021-08-03 | Amazon Technologies, Inc. | Evaluating joints using vibrometric signatures |
US11562483B2 (en) * | 2019-06-27 | 2023-01-24 | University Of Rochester | 2D shear wave dispersion imaging using a reverberant shear wave field |
JP7077433B1 (en) | 2021-02-10 | 2022-05-30 | ジーイー・プレシジョン・ヘルスケア・エルエルシー | Control program for ultrasonic diagnostic equipment, ultrasonic diagnostic system and ultrasonic diagnostic equipment |
CN113633313B (en) * | 2021-08-23 | 2024-07-09 | 深圳欢影医疗科技有限公司 | Acoustic radiation force elastic imaging method and system based on multi-frequency ultrasonic probe |
Citations (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3023609A (en) | 1958-06-04 | 1962-03-06 | Gen Motors Corp | Flaw detection apparatus |
US4061017A (en) | 1975-11-17 | 1977-12-06 | Time/Data Corporation | Structural analysis system |
US4947851A (en) | 1988-02-19 | 1990-08-14 | Institute for Physical Chemistry | Method and device for acoustic testing of elasticity of biological tissues |
US5086775A (en) | 1990-11-02 | 1992-02-11 | University Of Rochester | Method and apparatus for using Doppler modulation parameters for estimation of vibration amplitude |
US5487306A (en) | 1994-09-26 | 1996-01-30 | General Electric Company | Phase aberration correction in phased-array imaging systems |
US5505090A (en) | 1993-11-24 | 1996-04-09 | Holographics Inc. | Method and apparatus for non-destructive inspection of composite materials and semi-monocoque structures |
US5606971A (en) | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
US5672830A (en) | 1994-10-04 | 1997-09-30 | Massachusetts Institute Of Technology | Measuring anisotropic mechanical properties of thin films |
US5810731A (en) | 1995-11-13 | 1998-09-22 | Artann Laboratories | Method and apparatus for elasticity imaging using remotely induced shear wave |
US5903516A (en) | 1996-05-08 | 1999-05-11 | Mayo Foundation For Medical Education And Research | Acoustic force generator for detection, imaging and information transmission using the beat signal of multiple intersecting sonic beams |
US5921928A (en) | 1996-12-05 | 1999-07-13 | Mayo Foundation For Medical Education And Research | Acoustic force generation by amplitude modulating a sonic beam |
US5974881A (en) | 1997-07-16 | 1999-11-02 | The Trustees Of The Stevens Institute Of Technology | Method and apparatus for acoustic detection of mines and other buried man-made objects |
US5974889A (en) | 1998-01-02 | 1999-11-02 | General Electric Company | Ultrasonic multi-transducer rotatable scanning apparatus and method of use thereof |
US5991239A (en) | 1996-05-08 | 1999-11-23 | Mayo Foundation For Medical Education And Research | Confocal acoustic force generator |
US5995620A (en) * | 1995-02-15 | 1999-11-30 | Telefonaktiebolaget Lm Ericsson | Echo canceller having Kalman filter for optimal adaptation |
WO2000010638A2 (en) * | 1998-08-24 | 2000-03-02 | Baskent University | An asynchronous oversampling beamformer |
US6045504A (en) | 1998-06-09 | 2000-04-04 | Muzilla; David John | Method and apparatus for polynomial approximation of nonlinear operations in medical ultrasound imaging |
US6068597A (en) | 1999-04-13 | 2000-05-30 | Lin; Gregory Sharat | Vibrational resonance ultrasonic Doppler spectrometer and imager |
US6245016B1 (en) | 1999-03-12 | 2001-06-12 | General Electric Company | Ultrasound imaging system having post-beamformer signal processing using deconvolution algorithm |
US20020007118A1 (en) | 2000-03-15 | 2002-01-17 | Hideo Adachi | Ultrasonic wave transducer system and ultrasonic wave transducer |
US20020095087A1 (en) | 2000-11-28 | 2002-07-18 | Mourad Pierre D. | Systems and methods for making noninvasive physiological assessments |
US6711954B2 (en) | 2001-01-19 | 2004-03-30 | Lockheed Martin Corporation | Method and apparatus for improving the dynamic range of laser detected ultrasound in attenuative materials |
US6715354B2 (en) | 1998-02-24 | 2004-04-06 | Massachusetts Institute Of Technology | Flaw detection system using acoustic doppler effect |
US6764448B2 (en) | 2002-10-07 | 2004-07-20 | Duke University | Methods, systems, and computer program products for imaging using virtual extended shear wave sources |
US20040167403A1 (en) | 2000-04-05 | 2004-08-26 | Nightingale Kathryn R. | Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing |
US20040260180A1 (en) | 2001-08-20 | 2004-12-23 | Hiroshi Kanai | Tissue identifying method in ultrasonography and ultrasonograph |
US20050075565A1 (en) | 2003-09-19 | 2005-04-07 | Fuji Photo Film Co., Ltd. | Ultrasonic transmission/reception apparatus |
US6899680B2 (en) | 2000-10-19 | 2005-05-31 | Odetect As | Ultrasound measurement techniques for bone analysis |
US20050119568A1 (en) | 2003-10-14 | 2005-06-02 | Salcudean Septimiu E. | Method for imaging the mechanical properties of tissue |
US20050165306A1 (en) | 2003-10-03 | 2005-07-28 | Yi Zheng | Detection of motion in vibro-acoustography |
US20050167596A1 (en) | 2002-08-28 | 2005-08-04 | Siemens Westinghouse Power Corporation | System and method for multiple mode flexible excitation in sonic infrared imaging |
US6951544B2 (en) | 2000-04-05 | 2005-10-04 | Duke University | Method and apparatus for the identification and characterization of regions of altered stiffness |
US20050252295A1 (en) | 2002-09-02 | 2005-11-17 | Centre National De La Recherche Scientifique- Cnrs | Imaging method and device employing sherar waves |
US20060000281A1 (en) | 2000-07-21 | 2006-01-05 | Harris Paul D | Method and apparatus for assessing or predicting characteristics of wood or other wooden materials |
US6984209B2 (en) | 2003-07-02 | 2006-01-10 | The Brigham And Women's Hospital, Inc. | Harmonic motion imaging |
US20060036174A1 (en) * | 2004-07-26 | 2006-02-16 | Siemens Medical Solutions Usa, Inc. | Contrast agent imaging with agent specific ultrasound detection |
US20060169029A1 (en) | 2004-12-20 | 2006-08-03 | Heyman Joseph S | Acoustic concealed item detector |
US20070038095A1 (en) * | 2003-10-03 | 2007-02-15 | Greenleaf James F | Ultrasound vibrometry |
US20080302187A1 (en) | 2004-12-20 | 2008-12-11 | Mayo Foundation For Medical Education And Research | Vibroacoustic System for Vibration Testing |
US20100010346A1 (en) | 2008-04-07 | 2010-01-14 | Greenleaf James F | System and Method For Non-Invasively Measuring Tissue Viscoelasticity Using Surface Waves |
-
1998
- 1998-07-02 ZA ZA985834A patent/ZA985834B/en unknown
-
2009
- 2009-05-15 US US12/992,228 patent/US8659975B2/en active Active
Patent Citations (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3023609A (en) | 1958-06-04 | 1962-03-06 | Gen Motors Corp | Flaw detection apparatus |
US4061017A (en) | 1975-11-17 | 1977-12-06 | Time/Data Corporation | Structural analysis system |
US4947851A (en) | 1988-02-19 | 1990-08-14 | Institute for Physical Chemistry | Method and device for acoustic testing of elasticity of biological tissues |
US5086775A (en) | 1990-11-02 | 1992-02-11 | University Of Rochester | Method and apparatus for using Doppler modulation parameters for estimation of vibration amplitude |
US5505090A (en) | 1993-11-24 | 1996-04-09 | Holographics Inc. | Method and apparatus for non-destructive inspection of composite materials and semi-monocoque structures |
US5487306A (en) | 1994-09-26 | 1996-01-30 | General Electric Company | Phase aberration correction in phased-array imaging systems |
US5672830A (en) | 1994-10-04 | 1997-09-30 | Massachusetts Institute Of Technology | Measuring anisotropic mechanical properties of thin films |
US5995620A (en) * | 1995-02-15 | 1999-11-30 | Telefonaktiebolaget Lm Ericsson | Echo canceller having Kalman filter for optimal adaptation |
US5606971A (en) | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
US5810731A (en) | 1995-11-13 | 1998-09-22 | Artann Laboratories | Method and apparatus for elasticity imaging using remotely induced shear wave |
US5903516A (en) | 1996-05-08 | 1999-05-11 | Mayo Foundation For Medical Education And Research | Acoustic force generator for detection, imaging and information transmission using the beat signal of multiple intersecting sonic beams |
US5991239A (en) | 1996-05-08 | 1999-11-23 | Mayo Foundation For Medical Education And Research | Confocal acoustic force generator |
US5921928A (en) | 1996-12-05 | 1999-07-13 | Mayo Foundation For Medical Education And Research | Acoustic force generation by amplitude modulating a sonic beam |
US5974881A (en) | 1997-07-16 | 1999-11-02 | The Trustees Of The Stevens Institute Of Technology | Method and apparatus for acoustic detection of mines and other buried man-made objects |
US5974889A (en) | 1998-01-02 | 1999-11-02 | General Electric Company | Ultrasonic multi-transducer rotatable scanning apparatus and method of use thereof |
US6715354B2 (en) | 1998-02-24 | 2004-04-06 | Massachusetts Institute Of Technology | Flaw detection system using acoustic doppler effect |
US6045504A (en) | 1998-06-09 | 2000-04-04 | Muzilla; David John | Method and apparatus for polynomial approximation of nonlinear operations in medical ultrasound imaging |
WO2000010638A2 (en) * | 1998-08-24 | 2000-03-02 | Baskent University | An asynchronous oversampling beamformer |
US6245016B1 (en) | 1999-03-12 | 2001-06-12 | General Electric Company | Ultrasound imaging system having post-beamformer signal processing using deconvolution algorithm |
US6068597A (en) | 1999-04-13 | 2000-05-30 | Lin; Gregory Sharat | Vibrational resonance ultrasonic Doppler spectrometer and imager |
US20020007118A1 (en) | 2000-03-15 | 2002-01-17 | Hideo Adachi | Ultrasonic wave transducer system and ultrasonic wave transducer |
US6951544B2 (en) | 2000-04-05 | 2005-10-04 | Duke University | Method and apparatus for the identification and characterization of regions of altered stiffness |
US20040167403A1 (en) | 2000-04-05 | 2004-08-26 | Nightingale Kathryn R. | Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing |
US20060000281A1 (en) | 2000-07-21 | 2006-01-05 | Harris Paul D | Method and apparatus for assessing or predicting characteristics of wood or other wooden materials |
US6899680B2 (en) | 2000-10-19 | 2005-05-31 | Odetect As | Ultrasound measurement techniques for bone analysis |
US20020095087A1 (en) | 2000-11-28 | 2002-07-18 | Mourad Pierre D. | Systems and methods for making noninvasive physiological assessments |
US6711954B2 (en) | 2001-01-19 | 2004-03-30 | Lockheed Martin Corporation | Method and apparatus for improving the dynamic range of laser detected ultrasound in attenuative materials |
US20040260180A1 (en) | 2001-08-20 | 2004-12-23 | Hiroshi Kanai | Tissue identifying method in ultrasonography and ultrasonograph |
US20050167596A1 (en) | 2002-08-28 | 2005-08-04 | Siemens Westinghouse Power Corporation | System and method for multiple mode flexible excitation in sonic infrared imaging |
US20050252295A1 (en) | 2002-09-02 | 2005-11-17 | Centre National De La Recherche Scientifique- Cnrs | Imaging method and device employing sherar waves |
US6764448B2 (en) | 2002-10-07 | 2004-07-20 | Duke University | Methods, systems, and computer program products for imaging using virtual extended shear wave sources |
US6984209B2 (en) | 2003-07-02 | 2006-01-10 | The Brigham And Women's Hospital, Inc. | Harmonic motion imaging |
US20050075565A1 (en) | 2003-09-19 | 2005-04-07 | Fuji Photo Film Co., Ltd. | Ultrasonic transmission/reception apparatus |
US20050165306A1 (en) | 2003-10-03 | 2005-07-28 | Yi Zheng | Detection of motion in vibro-acoustography |
US20070038095A1 (en) * | 2003-10-03 | 2007-02-15 | Greenleaf James F | Ultrasound vibrometry |
US7753847B2 (en) | 2003-10-03 | 2010-07-13 | Mayo Foundation For Medical Education And Research | Ultrasound vibrometry |
US7785259B2 (en) | 2003-10-03 | 2010-08-31 | Mayo Foundation For Medical Education And Research | Detection of motion in vibro-acoustography |
US20050119568A1 (en) | 2003-10-14 | 2005-06-02 | Salcudean Septimiu E. | Method for imaging the mechanical properties of tissue |
US20060036174A1 (en) * | 2004-07-26 | 2006-02-16 | Siemens Medical Solutions Usa, Inc. | Contrast agent imaging with agent specific ultrasound detection |
US20060169029A1 (en) | 2004-12-20 | 2006-08-03 | Heyman Joseph S | Acoustic concealed item detector |
US20080302187A1 (en) | 2004-12-20 | 2008-12-11 | Mayo Foundation For Medical Education And Research | Vibroacoustic System for Vibration Testing |
US20100010346A1 (en) | 2008-04-07 | 2010-01-14 | Greenleaf James F | System and Method For Non-Invasively Measuring Tissue Viscoelasticity Using Surface Waves |
Non-Patent Citations (17)
Title |
---|
Amendment (Response to Aug. 18, 2008 Office Action), U.S. Appl. No. 11/536,330, filed Nov. 18, 2008. |
Amendment (Response to Nov. 29, 2007 Office Action), U.S. Appl. No. 11/536,330, filed Feb. 29, 2008. |
Amendment (Response to Oct. 8, 2009 Office Action), U.S. Appl. No. 11/536,330, filed Jan. 8, 2010. |
Appellant's Brief on Appeal, U.S. Appl. No. 11/536,330, filed Jun. 8, 2009. |
Chen, et al., Complex Stiffness Quantification Using Ultrasound Stimulated Vibrometry, 2003 IEEE Ultrasonics Symposium, pp. 941-944. |
International Search Report under date of mailing of Jul. 17, 2009 corresponding to PCT/US2009/044163. |
Michishita, et al., Ultrasonic Measurement of Minute Displacement of Object Cyclically Actuated by Acoustic Radiation Force, Jpn. J. Appl. Phys., 2003, 42:4608-4612. |
Nightingale, et al., On the Feasibility of Remote Palpation Using Acoustic Radiation Force, J. Acoust. Soc. Am., 2001, 110(1):625-634. |
PCT International Search Report and Written Opinion, PCT/US05/45964, Jun. 21, 2006. |
Ruvolo, et al., Skin Viscoelasticity Displays Site-and Age-Dependent Angular Anisotropy, Skin Pharmacol. Physiol., 2007, 20:313-321. |
United States Patent and Trademark Office, Notice of Allowance, U.S. Appl. No. 11/536,330, filed May 7, 2010. |
United States Patent and Trademark Office, Office Action Summary and Detailed Action, U.S. Appl. No. 11/536,330, filed Aug. 18, 2008. |
United States Patent and Trademark Office, Office Action Summary and Detailed Action, U.S. Appl. No. 11/536,330, filed Feb. 6, 2009. |
United States Patent and Trademark Office, Office Action Summary and Detailed Action, U.S. Appl. No. 11/536,330, filed Nov. 29, 2007. |
United States Patent and Trademark Office, Office Action Summary and Detailed Action, U.S. Appl. No. 11/536,330, filed Oct. 8, 2009. |
Vexler, et al., Evaluation of Skin Viscoelasticity and Anisotropy by Measurement of Speed of Shear Wave Propagation with Viscoelasticity Skin Analyzer, J. Invest Dermatol., 1999, 113:732-739. |
Zhang, et al., Excitation and Measurement of Flexural Waves in Arterial Vessels, 2003 IEEE Ultrasonics Symposium, pp. 1883-1886. |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11857373B2 (en) | 2016-10-04 | 2024-01-02 | The Regents Of The University Of California | Multi-frequency harmonic acoustography for target identification and border detection |
US10024757B1 (en) | 2017-05-04 | 2018-07-17 | United Launch Alliance, L.L.C. | Non-uniform sampling in bandwidth constrained data acquisition systems |
Also Published As
Publication number | Publication date |
---|---|
ZA985834B (en) | 1999-01-14 |
US20110063950A1 (en) | 2011-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8659975B2 (en) | Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions | |
US7753847B2 (en) | Ultrasound vibrometry | |
US9986973B2 (en) | Method for shear wave ultrasound vibrometry with interleaved push and detection pulses | |
US7785259B2 (en) | Detection of motion in vibro-acoustography | |
EP0008517B1 (en) | Duplex ultrasonic imaging system with repetitive excitation of common transducer in doppler modality | |
JP4150866B2 (en) | Method of operating an ultrasound imaging system | |
US8602994B2 (en) | Method for ultrasound vibrometry using orthogonal basis functions | |
US5806520A (en) | Method and device for evaluating and characterizing the properties of bone | |
Jensen | Medical ultrasound imaging | |
US8734352B2 (en) | Spatially-fine shear wave dispersion ultrasound vibrometry sampling | |
US20210251612A1 (en) | Ultrasound imaging with spectral compounding for speckle reduction | |
US8021302B2 (en) | Ultrasonic apparatus and ultrasonic diagnostic method | |
JP5715594B2 (en) | Method and apparatus for flow parameter imaging | |
US5172343A (en) | Aberration correction using beam data from a phased array ultrasonic scanner | |
WO2009140607A1 (en) | Vibration generation and detection in shear wave dispersion ultrasound vibrometry with large background motions | |
EP2287632A1 (en) | Ultrasound imaging using non-linear manipulation of forward propagation properties of a pulse | |
US20050277835A1 (en) | Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties | |
EP0146073B1 (en) | Ultrasonic diagnosing apparatus | |
JP2005536309A (en) | System and method for improved harmonic imaging | |
JP2001252276A (en) | Ultrasonic imaging device | |
US7666142B2 (en) | Ultrasound doppler diagnostic apparatus and image data generating method | |
JP2002534187A (en) | Development of Doppler angle in ultrasonic color flow / Doppler | |
US11369343B2 (en) | Ultrasound imaging by nonlinear localization | |
JPH09224938A (en) | Ultrasonic diagnostic device and method for optimizing delay time | |
US6533728B1 (en) | Method and apparatus for recovery and parametric display of contrast agents in ultrasound imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: MAYO FOUNDATON FOR MEDICAL EDUCATION AND RESEARCH, Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:GREENLEAF, JAMES F.;CHEN, SHIGAO;ZHENG, YI;AND OTHERS;REEL/FRAME:025352/0846 Effective date: 20101008 |
|
AS | Assignment |
Owner name: NATIONAL INSTITUTES OF HEALTH (NIH), U.S. DEPT. OF Free format text: CONFIRMATORY LICENSE;ASSIGNOR:MAYO FOUNDATION FOR MEDICAL EDUCATION AND RESEARCH;REEL/FRAME:025656/0457 Effective date: 20101118 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551) Year of fee payment: 4 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 8 |