Search Images Maps Play YouTube News Gmail Drive More »
Sign in
Screen reader users: click this link for accessible mode. Accessible mode has the same essential features but works better with your reader.

Patents

  1. Advanced Patent Search
Publication numberUSRE43048 E1
Publication typeGrant
Application numberUS 09/481,814
Publication dateDec 27, 2011
Filing dateJan 11, 2000
Priority dateOct 10, 1995
Also published asUS5706819, US5951478
Publication number09481814, 481814, US RE43048 E1, US RE43048E1, US-E1-RE43048, USRE43048 E1, USRE43048E1
InventorsJuin-Jet Hwang, David Hope Simpson
Original AssigneeAdvanced Technology Laboratories, Inc.
Export CitationBiBTeX, EndNote, RefMan
External Links: USPTO, USPTO Assignment, Espacenet
Ultrasonic diagnostic imaging with harmonic contrast agents
US RE43048 E1
Abstract
Apparatus and methods are disclosed for the detection and imaging of ultrasonic harmonic contrast agents. The harmonic echo effect is detected through alternate polarity acquisition of harmonic contrast agent effects, which provides the benefits of suppressing the harmonic components of the transmitted signal while eliminating clutter.
Images(6)
Previous page
Next page
Claims(25)
1. A method of ultrasonically detecting the ultrasonic response of an ultrasonic contrast agent comprising the steps of:
transmitting a first ultrasonic pulse to said ultrasonic contrast agent to cause a first harmonic response;
transmitting a second ultrasonic pulse of a different polarity than said first ultrasonic pulse to said harmonic contrast agent to cause a second harmonic response;
detecting said first and second harmonic responses; and
combining said first and second harmonic responses.
2. The method of claim 1, wherein said step of combining comprises summing said first and second harmonic responses.
3. The method of claim 1, wherein said step of combining comprises integrating said first and second harmonic responses.
4. The method of claim 1, wherein said transmitting step comprises transmitting pulses which exhibit a pulse energy which is within a range which causes microbubbles of said ultrasonic contrast agent to oscillate without substantial microbubble destruction.
5. A method of ultrasonically detecting the nonlinear response of a substance within the body comprising the steps of:
transmitting at least three ultrasonic pulses into the body which exhibit first and second characteristics that cause a reduction in the linear echo response when echoes received in response to such pulses are combined;
receiving echoes in response to said ultrasonic pulses; and
combining said echoes to produce a nonlinear response.
6. The method of claim 5, wherein said step of receiving echoes comprises receiving echoes from a given location in the body.
7. The method of claim 5, wherein said ultrasonic pulses are transmitted in a sequence in which said first and second characteristics are alternated from pulse to pulse.
8. The method of claim 5, wherein said step of combining comprises summing pairs of echoes.
9. The method of claim 5, wherein said ultrasonic pulses are transmitted in a sequence in which said first and second characteristics are alternated from pulse to pulse; and
wherein said step of combining comprises summing pairs of echoes from successive pulses.
10. The method of claim 5, wherein said first and second characteristics comprise first and second polarities.
11. The method of claim 10, wherein said transmitted ultrasonic pulses are of the form {p −p p . . . }.
12. The method of claim 5, wherein said step of combining produces a sum result S which is substantially equal to
S = j = 1 n - 1 ( E j + E j + 1 ) ,
where Ej and Ej+1 are pulse echoes.
13. The method of claim 12, wherein the number of ultrasonic pulses which is transmitted is three.
14. A method of ultrasonically detecting the nonlinear ultrasonic response of a medium inside the body comprising the steps of:
transmitting a first ultrasonic pulse to said medium to cause a first echo response;
transmitting a second ultrasonic pulse to said medium to cause a second echo response;
transmitting a third ultrasonic pulse to said medium to cause a third echo response which is substantially the same as said first echo response; and
combining said first, second and third echo responses to produce a nonlinear response.
15. The method of claim 14, wherein said transmitted ultrasonic pulses are of the form {p −p p}.
16. The method of claim 14, wherein said step of combining produces a sum result S which is substantially equal to
S = j = 1 n - 1 ( E j + E j + 1 ) ,
where Ej and Ej+1 are pulse echoes.
17. A method of ultrasonically detecting the nonlinear response of a substance within the body comprising the steps of:
transmitting at least three ultrasonic pulses into the body in a sequence which is of the form {p −p p −p . . . −p p};
receiving echoes in response to said ultrasonic pulses which comprise a sequence of the form {E1 E2 E3 E4 . . . En-1 En}; and
accumulating said echoes to produce a nonlinear response.
18. The method of claim 17, wherein said step of accumulating comprises accumulating pairs of consecutive echoes.
19. The method of claim 17, wherein said step of accumulating produces a sum result S which is substantially equal to
S = j = 1 n - 1 ( E j + E j + 1 ) ,
where Ej and Ej+1 are pulse echoes.
20. A method of ultrasonically detecting the nonlinear response of a substance within the body comprising the steps of:
transmitting a sequence of at least three ultrasonic pulses into the body which exhibit a transmit characteristic which alternates from pulse to pulse;
receiving echoes in response to said ultrasonic pulses; and
combining said echoes to produce a nonlinear response.
21. The method of claim 20, wherein said pulses are transmitted to a given location in the body; and
wherein said step of combining reduces the primary component of said echoes and produces a harmonic response.
22. The method of claim 20, wherein said step of transmitting produces a sequence of echoes relating to a given location in the body in which the phase of the primary component of echoes produced by one transmit characteristic is out of phase with the phase of the primary component of echoes produced by the alternate transmit characteristic.
23. The method of claim 22, wherein said step of combining reduces the primary component of the combined echoes and produces a harmonic response.
24. The method of claim 23, wherein said transmit characteristic is a polarity differential from pulse to pulse.
25. The method of claim 23, wherein said transmit characteristic is a phase differential from pulse to pulse.
Description

This application claims the benefit of U.S. Provisional Application No. 60/005,009, filed Oct. 10, 1995.

This invention relates to ultrasonic diagnosis and imaging of the body with ultrasonic contrast agents and, in particular, to new methods and apparatus for ultrasonically detecting and imaging with contrast agents.

Ultrasonic diagnostic imaging systems are capable of imaging and measuring the physiology within the body in a completely noninvasive manner. Ultrasonic waves are transmitted into the body from the surface of the skin and are reflected from tissue and cells within the body. The reflected echoes are received by an ultrasonic transducer and processed to produce an image or measurement of blood flow. Diagnosis is thereby possible with no intervention into the body of the patient.

However materials known as ultrasonic contrast agents can be introduced into the body to enhance ultrasonic diagnosis. Contrast agents are substances which will strongly interact with ultrasonic waves, returning echoes which may be clearly distinguished from those returned by blood and tissue. One class of substances which has been found to be especially useful as an ultrasonic contrast agent is gases, in the form of tiny bubbles called microbubbles. Microbubbles present a significant acoustic impedance mismatch in comparison to tissue and fluids, and nonlinear behavior in certain acoustic fields which is readily detectable through special ultrasonic processing. In order to infuse bubbles into the body so that they will survive passage through the pulmonary system and circulate throughout the vascular system, gases have been stabilized in solutions in the form of tiny microbubbles. Microbubble contrast agents are useful for imaging the body's vascular system, for instance, as the contrast agent can be injected into the bloodstream and will pass through the veins and arteries of the body with the blood supply until filtered from the blood stream in the lungs, kidneys and liver.

One property of microbubble contrast agents currently under investigation is harmonic response. These harmonic contrast agents exhibit significant, detectable responses at frequencies which are harmonics of the transmitted ultrasonic frequency. This property is useful for clutter rejection of the received signals. When the transmitted frequency band is used as the received frequency band, echoes will be returned from the microbubbles, but also from surrounding tissue, the latter comprising clutter in the received echo signals. But with harmonic contrast agents, reception occurs at harmonic frequencies, where fundamental band clutter from tissue is ignored. Since tissue generally reflects very minimal harmonic components, the received harmonic band enables the microbubble echoes to be received with a high signal to noise ratio.

In accordance with the principles of present invention, a technique is provided for the detection and imaging of harmonic ultrasonic contrast agents. The harmonic contrast agent is insonified by alternate polarity transmitted pulses, and the echo signals received from the transmitted pulses are combined. The result is a suppression of harmonic components of the transmitted ultrasonic waves and the elimination of clutter.

In the drawings:

FIG. 1 illustrates in block diagram form apparatus constructed in accordance with the present invention which provides performance advantages for harmonic contrast agent detection;

FIGS. 2 and 3 illustrate passband characteristics used to explain the performance of the embodiment of FIG. 1;

FIGS. 4a and 4b illustrates the alternate polarity pulsing of harmonic contrast agents; and

FIGS. 5a5c illustrate nonlinear response waveforms produced by alternate polarity acquisition of contrast agent echoes.

Referring to FIG. 1, an ultrasonic diagnostic system for use with harmonic contrast agents in accordance with the present invention is shown in block diagram form. In this system an array transducer 112 of a probe 110 transmits ultrasonic energy and receives echoes returned in response to this transmission. The response characteristic of the transducer can exhibit two passbands, one around the central transmit frequency and another about the center of the received passband. For imaging harmonic contrast agents, a broadband transducer having a passband encompassing both the transmit and receive passbands is preferred. The transducer may be manufactured and tuned to exhibit a response characteristic as shown in FIG. 2, in which the lower hump 60 of the response characteristic is centered about the center transmit frequency ft, and the upper hump 62 is centered about the center frequency fr of the response passband. The transducer response characteristic of FIG. 3 is preferred, however, as the single dominant characteristic 64 allows the probe to be suitable for both harmonic contrast imaging and imaging without harmonic contrast agents. The characteristic 64 encompasses the central transmit frequency ft, and also the harmonic receive passband bounded between frequencies fL and fc, and centered about frequency fr. A typical harmonic contrast agent can have a response such that transmission about a central transmit frequency of 1.7 MHz will result in harmonic returning echo signals about a frequency of 3.4 MHz. A bandwidth characteristic 64 of approximately 2 MHz would be suitable for these harmonic frequencies.

In FIG. 1 a central controller 120 provides a control signal ftr to a transmit frequency control circuit or pulser 121 to control the center frequency and time of transmission of the transmitted ultrasonic energy. The transmit frequency control circuit pulses the elements of the transducer array 112 by means of a transmit/receive switch 114.

Echoes received by the transducer array 112 are coupled through the T/R switch 114 and digitized by analog to digital converters 115. The sampling frequency fs of the A/D converters 115 is controlled by the central controller. The desired sampling rate dictated by sampling theory is at least twice the highest frequency fc of the received passband and, for the preceding exemplary frequencies, might be on the order of at least 8 MHz. Sampling rates higher than the minimum requirement are also desirable.

The echo signal samples from the individual transducer elements are delayed and summed by a beamformer 116 to form coherent echo signals. The digital coherent echo signals are then filtered by a digital filter 118. In this embodiment, the transmit frequency ftr is not tied to the receiver, and hence the receiver is free to receive a band of frequencies which is separate from the transmitted band. The digital filter 118 bandpass filters the signals in the passband bounded by frequencies fL and fc in FIG. 3, and can also shift the frequency band to a lower or baseband frequency range. The digital filter could be a filter with a 1 MHz passband and a center frequency of 3.4 MHz in the above example. A preferred digital filter is a parallel arrangement of serially coupled multipliers and accumulators. This arrangement is controlled by the central controller 120, which provides multiplier weights and decimation control which control the characteristics of the digital filter. Preferably the arrangement is controlled to operate as a finite impulse response (FIR) filter, and performs both filtering and decimation.

Filtered echo signals from tissue, generally filtered by a passband centered about or demodulated from the transmit frequency, are coupled to a B mode processor 37 for conventional B mode processing. Filtered echo signals of the harmonic contrast agent passband are coupled to a contrast signal detector 128 which performs pulse to pulse summation or integration of temporally discrete echoes from a given spatial location, amplitude or envelope detects the combined signals. Simple two pulse summation of the form P1+P2 may be employed where P1 represents the echoes received following one pulse and P2 represents the echoes received following another pulse. The combination of echoes from consecutive pulses may, if desired, be performed before the digital filter 118 rather than after, the decision being a matter of choice of system design.

The filtered echo signals from the digital filter 118 are also coupled to a Doppler processor 130 for conventional Doppler processing to produce velocity and power Doppler signals. The outputs of these processors are coupled to a 3D image rendering processor 132 for the rendering of three dimensional images, which are stored in a 3D image memory 134. Three dimensional rendering may be performed as described in U.S. patent application Ser. No. 08/638,710, and in U.S. Pat. Nos. 5,474,073 and 5,485,842, the latter two patents illustrating three dimensional power Doppler ultrasonic imaging techniques. The signals from the contrast signal detector 128, the processors 37 and 130, and the three dimensional image signals are coupled to a video processor 140 where they may be selected for display on an image display 50 as dictated by user selection. The video processor preferably includes persistence processing, whereby momentary intensity peaks of detected contrast agents can be sustained in the image. One technique for providing persistence is through frame averaging, whereby new image frames are combined with previous frame information on a spatial basis. The combination can be done by weighting the contributions of the old and new frame information and the frame information can be combined in a recursive manner; that is, old frame information is fed back for combining with new frame information. A preferred persistence technique is the fast attack, slow decay technique described in U.S. Pat. No. 5,215,094, which can be applied to both Doppler and contrast agent images.

The apparatus of FIG. 1 performs alternate polarity pulse transmission as illustrated in FIGS. 4a and 4b. In the first transmission of FIG. 4a, the central controller 120 provides a first polarity control signal ftr1 to the pulser 121, which drives the transducer elements 112 to transmit a first polarity pulse 20. For the second transmission of FIG. 4b, the central controller 120 provides a second polarity control signal ftr2 to the pulser 121, which drives the transducer elements 112 to transmit a second polarity pulse 22.

The echoes received from microbubbles in response to these alternate polarity transmissions are shown in FIGS. 5a and 5b. FIG. 5a illustrates an echo waveform 310 received from the first pulsing of a microbubble contrast agent. The nonuniform amplitudes on either side of the zero reference level illustrate nonlinear reflexive action of microbubbles in the presence of acoustic waves, as the microbubbles nonlinearly compress and expand. The echo waveform of 310 FIG. 5a results from transmission of an ultrasonic pulse exhibiting a first polarity.

Following transmission of the ultrasonic pulse exhibiting the opposite polarity, the echo waveform 312 of FIG. 5b results. This waveform is similarly nonlinear, but out of phase with the first waveform due to the change in pulse polarity. When the two waveforms are combined, a harmonic response is obtained, as shown in FIG. 5c. The highly nonlinear waveform of FIG. 5c is readily detected, causing the system to become highly sensitive to the contrast agent which produced the nonlinear echo responses.

A mathematical analysis of this effect and response is as follows. To detect the harmonic response of microbubbles, the harmonic component in the incident pressure wave must be suppressed. Based on the analytical solution of the dynamic motion of microbubbles, the primary component of the backscattering pressure magnitude is linearly proportional to the incident pressure and the harmonic component is quadratically proportional to the incident pressure pi or ps(ω)αpi and ps(2ω)αpi 2. Thus, neglecting the higher order terms, one may write the backscattering pressure magnitude pB(ω) from a microbubble in a generic form
PB(ω)=k1(ω)p+k2(ω)p2  (1)
where k1 and k2 are parametrically related to the acoustic properties of the microbubble such as size, viscosity, surface tension, ambient pressure, etc.

Now assume that the microbubble is excited by two narrow band signals at different times but with the same magnitude p and at the same frequency ω, but with opposite polarity: pi1=p cos ωt and pi2=−p cos ωt. Then the back-scattered pressure wave from pi1=p cos ωt is
pB1(ω,t)=k1(ω,t)p+k2(ω,t)p2  (2)
and from pi1=−p cos ωt is
pB2(ω,t+δt)=k1(ω,t+δt)p+k2(ω,t+δt)p2  (3)
Then the total backscattered pressure magnitude may be obtained by summing Equations (2) and (3),
S=pB1+pB2=(k1(ω,t)−k1(ω,t+δt))p+(k2(ω,t)+k2(ω,t+δt))p2  (4)
≈2k2(ω)p2
Equation (4) shows that the primary component is eliminated if k1(ω) and k2(ω) do not change substantially in the time duration δt, where δt is small.

Assume the backscattering from microbubbles is quasi-stationary over T, where T is the pulse repetition interval. Therefore, the average nonlinear acoustic properties are not changed over time T, or
E{k1(ω,t)}≅E{k1(ω,t+T)}
and
E{k2(ωt)}≅E{k2(ω,t+T)}.
The relationship of Equation (4) will hold by summing the pulse echoes from two pulses which are time-diverse in T. The quasi-stationary assumption is valid for slow perfused flow, such as myocardial perfusion.

When the bandwidth of the incident pressure wave is wide, the wideband excitation wave P(t) may be represented by a Fourier series

P ( t ) = i A ( ω i ) cos ω i t
Thus the backscattered pressure magnitude of the microbubbles from P(t) may be written as

P B 1 = i k 1 ( ω i ) A ( ω i ) + i k 2 ( ω i ) A 2 ( ω i ) ( 5 )
and the backscattered pressure magnitude of the microbubbles from −P(t) may be written as

P B 2 = - i k 1 ( ω i ) A ( ω i ) + i k 2 ( ω i ) A 2 ( ω i ) ( 6 )
Summing Equations (5) and (6), one may obtain

S = p B 1 + p B 2 = 2 i k 2 ( ω i ) A 2 ( ω i ) (7)
Again, the harmonic component is extracted and the primary component is eliminated.

Let us assume the nonlinearity in tissue is negligible. Since the backscattered pressure in a linear medium is linearly proportional to the incident pressure wave, the polarity of the backscattered wave will be changed as the polarity of the incident pressure wave is changed. Assuming the tissue is relatively stationary during the period of two consecutive pulses, summing the pulse echoes from consecutive pulses with opposite polarity will cancel the echo response from tissue. Thus, tissue clutter will be suppressed.

The concept of summing the pulse echoes from two pulses of opposite polarity may be generalized into processing echoes from multiple pulses with alternate polarity to maximize the sensitivity and minimize the variance, assuming the tissue is stationary during the pulsing interval. Let the pulse sequence be
P={p −p p− p p− p * *−p p}
and the pulse echoes be
E={E1 E2 E3 E4 E5 E6 En}
Accumulating the partial sum of consecutive pairs of echoes results in

S = j = 1 n - 1 E j + E j + 1 = 2 ( n - 1 ) i k 2 ( ω i ) A 2 ( ω i )

Patent Citations
Cited PatentFiling datePublication dateApplicantTitle
US4112411Dec 9, 1976Sep 5, 1978U.S. Phillips CorporationDevice for echography by means of focussed ultrasonic beams
US4119938Nov 6, 1975Oct 10, 1978Agence Nationale De Valorisation De La Rechere (Anvar)Methods and devices for ultrasonic imaging
US4282452Mar 15, 1979Aug 4, 1981Siemens AktiengesellschaftMethod and circuit arrangement for energizing ultrasonic transducers which are utilized in impulse echo technology
US4572203Jan 27, 1983Feb 25, 1986Feinstein Steven BContact agents for ultrasonic imaging
US4844082Feb 8, 1988Jul 4, 1989Matsushita Electric Industrial Co., Ltd.Ultrasonic examination apparatus
US4865042Aug 8, 1986Sep 12, 1989Hitachi, Ltd.Ultrasonic irradiation system
US5086775Nov 2, 1990Feb 11, 1992University Of RochesterMethod and apparatus for using Doppler modulation parameters for estimation of vibration amplitude
US5135000Sep 17, 1991Aug 4, 1992Raizot University Authority For Applied Research & Industrial Development Ltd.Method of measuring regional tissue blood flow
US5224481Sep 6, 1991Jul 6, 1993Ken IshiharaImage displaying method and device for realizing same in an ultrasonic diagnostic apparatus
US5233993Nov 19, 1991Aug 10, 1993Hitachi Medical CorporationUltrasonic diagnostic apparatus and image forming method therefor
US5241473Oct 11, 1991Aug 31, 1993Ken IshiharaUltrasonic diagnostic apparatus for displaying motion of moving portion by superposing a plurality of differential images
US5255683Dec 30, 1991Oct 26, 1993Sound Science Limited PartnershipMethods of and systems for examining tissue perfusion using ultrasonic contrast agents
US5302372Jul 27, 1992Apr 12, 1994National Science CouncilSonicating aqueous albumin-glucose solution with glycerine; forms gas-in-liquid microbubbles; injecting into patient
US5410516 *Jun 14, 1993Apr 25, 1995Schering AktiengesellschaftUltrasonic processes and circuits for performing them
US5453575 *Apr 28, 1994Sep 26, 1995Endosonics CorporationApparatus and method for detecting blood flow in intravascular ultrasonic imaging
US5456257Nov 23, 1994Oct 10, 1995Advanced Technology Laboratories, Inc.Ultrasonic detection of contrast agents
US5469849Jun 14, 1994Nov 28, 1995Kabushiki Kaisha ToshibaUltrasound diagnosis apparatus
US5482044Jan 14, 1992Jan 9, 1996Diasonics Ultrasound, Inc.Direct demodulation in ultrasound instruments
US5577505Feb 6, 1996Nov 26, 1996Hewlett-Packard CompanyMeans for increasing sensitivity in non-linear ultrasound imaging systems
US5632277Jun 28, 1996May 27, 1997Siemens Medical Systems, Inc.Ultrasound imaging system employing phase inversion subtraction to enhance the image
US5706819Oct 9, 1996Jan 13, 1998Advanced Technology Laboratories, Inc.Ultrasonic diagnostic imaging with harmonic contrast agents
US5902243Apr 15, 1998May 11, 1999Acuson CorporationUltrasonic imaging method with multiple pulse cancellation
US6506158Apr 9, 2001Jan 14, 2003Kabushiki Kaisha ToshibaUltrasonic diagnosis apparatus
WO1999030617A1Dec 18, 1998Jun 24, 1999AcusonDiagnostic ultrasound imaging method and system with improved frame rate
Non-Patent Citations
Reference
1Abiru, I. et al., "Nonlinear Propagation of Ultrasonic Pulses," (Technical Report of IEICE, US89-23, p. 53).
2Charles C. Church, "The effect of an elastic solid surface layer on the radial pulsation of gas bubbles," JASA v. 97, No. 3, Mar. 1995, 1510-21.
3D.L. Miller, "Ultrasonic detection of resonant cavitation bubbles in a flow tube by their second-harmonic emissions," Ultrasonics, Sep. 1981, pp. 217-224.
4Rielly, M.R., "A Theoretical and Experimental Investigation of Nonlinear Propagation of Ultrasound Through Tissue Mimicking Media," (Ph. D. Thesis, University of Bath Library, 2000).
Classifications
U.S. Classification600/458, 367/7, 424/9.51, 600/437, 424/9.5, 424/9.52
International ClassificationA61B8/00, G01S15/89, G01S7/52, A61B8/14, A61B8/08
Cooperative ClassificationG01S15/8963, G01S15/8959, A61B8/488, A61B8/463, G01S7/52026, G01S7/52038, A61B8/08, G01S7/52041, G01S7/52039, A61B8/06, A61B8/481
European ClassificationA61B8/48B, G01S15/89D4, G01S7/52S2F1, G01S7/52S2B1, A61B8/08T, A61B8/06, A61B8/46B4, G01S7/52S2F1B, G01S7/52S2F1C, G01S15/89D4D