US6132377A - Medical diagnostic ultrasonic imaging system and method using differential sub-band detection techniques - Google Patents

Medical diagnostic ultrasonic imaging system and method using differential sub-band detection techniques Download PDF

Info

Publication number
US6132377A
US6132377A US09/282,402 US28240299A US6132377A US 6132377 A US6132377 A US 6132377A US 28240299 A US28240299 A US 28240299A US 6132377 A US6132377 A US 6132377A
Authority
US
United States
Prior art keywords
signal
receive
signals
receive signals
ratio
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.)
Expired - Lifetime
Application number
US09/282,402
Inventor
Mirsaid S. Bolorforosh
Edward A. Gardner
Gregory L. Holley
Ting Lan Ji
Sriram Krishnan
Bhaskar S. Ramamurthy
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Medical Solutions USA Inc
Original Assignee
Acuson Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Acuson Corp filed Critical Acuson Corp
Priority to US09/282,402 priority Critical patent/US6132377A/en
Assigned to ACUSON CORPORATION reassignment ACUSON CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: BOLORFOROSH, MIRSAID S., KRISHNAN, SRIRAM, RAMAMURTHY, BHASKAR, HOLLEY, GREGORY L., JI, TING LAN, GARDNER, EDWARD A.
Application granted granted Critical
Publication of US6132377A publication Critical patent/US6132377A/en
Assigned to SIEMENS MEDICAL SOLUTIONS USA, INC. reassignment SIEMENS MEDICAL SOLUTIONS USA, INC. CHANGE OF NAME (SEE DOCUMENT FOR DETAILS). Assignors: SIEMENS MEDICAL SYSTEMS, INC.
Assigned to SIEMENS MEDICAL SOLUTIONS USA, INC. reassignment SIEMENS MEDICAL SOLUTIONS USA, INC. RE-RECORD TO CORRECT CONVEYING PARTY NAME PREVIOUSLY RECORDED AT REEL 024563 FRAME 0051 Assignors: ACUSON LLC, ACUSON CORPORATION
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/481Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • G01S15/8988Colour Doppler imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52038Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/024Mixtures
    • G01N2291/02491Materials with nonlinear acoustic properties

Definitions

  • This invention relates to medical diagnostic ultrasound imaging, and in particular to methods and systems for distinguishing contrast agent from surrounding tissue.
  • Contrast agent imaging is an important medical diagnostic ultrasound imaging mode.
  • One limitation in many contrast imaging systems is the difficulty of distinguishing echo signals from contrast agents from echo signals from surrounding tissue. This is because both contrast agents and tissue generate nonlinear return signals at frequencies other than the frequency of the insonifying signals. It is well recognized that nonlinear signals from tissue are generated by nonlinear propagation of the insonifying ultrasound wave.
  • Prior-art harmonic imaging methods that exploit non-linear behaviour of contrast agent include B-mode harmonic imaging, B-mode harmonic pulse inversion imaging, harmonic power Doppler imaging, and color harmonic pulse inversion imaging.
  • the signal is transmitted at a fundamental frequency f, and the receive signal is filtered to emphasize frequency components near the second harmonic, 2f.
  • Contrast agents are known to have a stronger second harmonic response than tissue, and for this reason the receive signal from contrast agent is enhanced over that from tissue. This method is currently used by most major manufacturers of ultrasound imaging equipment. Specific examples are described in Mine U.S. Pat. No. 5,724,976, Uhlendorf U.S. Pat. No. 5,410,516, Schutt U.S. Pat. No. 5,733,527, and "Simulated Capillary Blood Measurement Using a Nonlinear Ultrasonic Contrast Agent," Schrope et al.; Ultrasonic Imaging, Vol. 14, pp. 134-158, 1992.
  • B-mode harmonic pulse inversion imaging two pulses are transmitted along the same ultrasound line, where one pulse is shifted by 180° with respect to the other. Receive signals from the two pulses are then summed, and the resultant signal is displayed.
  • the pulse inversion technique cancels stationary fundamental frequency signals and retains second harmonic signals as well as some nonstationary fundamental signals. Examples of such methods include Chapman U.S. Pat. No. 5,632,277 and Hwang U.S. Pat. No. 5,706,819.
  • harmonic power Doppler imaging multiple pulses are transmitted along the same ultrasound line, and the receive signal is filtered about the second harmonic frequency. The resultant signals are then filtered with a high-pass filter to remove stationary signals. Examples of such methods are described in Averkiou U.S. Pat. No. 5,833,613, Johnson U.S. Pat. No. 5,456,257, and "Harmonic Power Mode Doppler Using Microbubble Contrast Agents: An Improved Method for Small Vessel Flow Imaging," Burns, et al., 1994 IEEE Ultrasonic Symposium, pp. 1547-1550, 1994.
  • B-mode harmonic imaging does not differentiate between second harmonic signals generated by nonlinear propagation through tissue and second harmonic signals generated by contrast agents.
  • Pulse inversion imaging further enhances second harmonic signals, but it still does not differentiate between second harmonic signals generated by nonlinear propagation through tissue and second harmonic signals generated by contrast agents.
  • Harmonic power Doppler imaging attempts to differentiate contrast agent from tissue by looking for a loss of correlation between successive pulses due to agent destruction, agent motion or other methods. However, tissue motion will also result in a loss of correlation and may appear as a displayed signal. Tissue motion is reduced currently by a combination of increasing the pulse repetition frequency and/or the use of more aggressive clutter filters. However, increasing the pulse repetition frequency may reduce the signal from destroyed contrast agent, and more aggressive clutter filters may filter out contrast signal as well as tissue signal.
  • the present inventors have discovered specific spectral properties of the backscattered fundamental and harmonic receive signals that can be used to distinguish contrast agent from tissue in a medical diagnostic ultrasound imaging system.
  • One preferred embodiment described below uses the ratio of the energy in the harmonic sub-band of the receive signal to the energy in the fundamental sub-band of the receive signal to detect the presence or absence of contrast agent. This ratio is then used to modulate the image signal to enhance the display of contrast agent.
  • FIG. 1 is a frequency versus magnitude plot of receive signals associated with tissue and contrast agent.
  • FIG. 2 is a block diagram showing a medical diagnostic ultrasound imaging system that incorporates a presently preferred embodiment to this invention.
  • FIG. 3 is a graph showing one example of the manner in which the gain signal can vary as a function of the ratio signal in the system of FIG. 3.
  • FIG. 4 is a flow chart of a preferred embodiment of the method of this invention, as practiced by the system of FIG. 1.
  • FIG. 5 in a flow chart is a second preferred embodiment of this invention.
  • contrast agent and tissue produce harmonic signals in response to fundamental transmit signals
  • the level of the harmonic signal with respect to the fundamental signal in the backscattered receive signal is significantly different for contrast agent as compared to tissue.
  • this property is used to discriminate between contrast agent and tissue.
  • FIG. 1 shows a receive signal R T associated with backscattered signals from tissue insonified at a transmit center frequency f.
  • the fundamental frequency f is 1.7 MHz and the harmonic frequency 2f is 3.4 MHz.
  • the receive signal R T peaks at the fundamental frequency f and has a magnitude at the harmonic frequency 2f that is suppressed by about 15 dB as compared to the magnitude at the fundamental frequency f.
  • FIG. 1 also shows a frequency versus magnitude plot of a receive signal R C associated with backscattered echo signals from contrast agent.
  • the receive signal R C is generated with a transmit signal centered at a fundamental frequency f of 1.7 MHz.
  • the receive signal R C also shows a peak magnitude at about the fundamental frequency f, and this peak magnitude is approximately 10 dB greater than the magnitude at the harmonic frequency 2f.
  • FIG. 1 shows that the ratio of amplitude near the fundamental frequency f to magnitude near the harmonic frequency 2f is substantially different for echo signals from tissue and contrast agent. As explained below, this ratio can be used to discriminate echo signals from tissue versus echo signals from contrast agent.
  • an ultrasonic imaging system 10 includes a transmitter 12 that generates transmit signals T.
  • the transmit signals T are applied via a transmit/receive switch 14 to a transducer probe 16.
  • the transducer probe 16 responds to the transmit signals T by generating ultrasonic energy at the frequency of the transmit signal T and directing this ultrasonic energy into an imaged region that includes tissue T and contrast agent C.
  • the tissue T and the contrast agent C backscatter ultrasonic energy from the imaged region to the transducer probe 16. This backscattered ultrasonic energy causes the transducer probe 16 to generate receive signals R C and R T that are applied via the transmit/receive switch 14 to a receiver 18.
  • the receiver 18 includes a receive beamformer 20 and two filters 22, 24.
  • the first filter 22 is a bandpass filter that passes a sub-band of the receive signals centered around the fundamental frequency f.
  • the transmit signals T are also centered around the fundamental frequency f.
  • the second filter 24 is another bandpass filter having a bandpass centered around the harmonic frequency 2f.
  • the output signals generated by the filters 22, 24 are fundamental receive signals F and harmonic receive signals H, respectively.
  • the fundamental and harmonic receive signals F, H are applied to a ratio former 26 that in this embodiment includes first and second detectors 28, 30.
  • the detectors 28, 30 preferably square the respective input signals and take the log of the squared result.
  • the detected output signals generated by the detectors 28, 30 are applied to integrators 32, 34 that integrate the detected output signals over frequency.
  • the output signals generated by the integrators 32, 34 are then subtracted from one another in a subtraction circuit 36 that in this application performs a division function.
  • the output signals of the integrators 32, 34 vary as the logarithm of the fundamental and harmonic receive signals, after detection and integration.
  • the difference between the output signals of the integrators 32, 34 creates a ratio signal R that varies as a function of the ratio of the detected, integrated fundamental and harmonic receive signals F, H.
  • the ratio signal R is applied to an image signal generator 38 that in this embodiment includes a lookup table 40.
  • the lookup table 40 generates a gain signal G in response to the ratio signal R, and this gain signal G is applied to a multiplier 42.
  • the multiplier 42 is also responsive to one of the fundamental and harmonic receive signals F, H, as determined by the state of a switch 44.
  • the output of the multiplier 42 is an image signal I that is applied to an image processor 46.
  • the output of the image processor 46 is applied to a display 48.
  • the lookup table 40 allows the gain applied to the fundamental or harmonic receive signal F, H to be varied as a linear or nonlinear function of the ratio signal R.
  • FIG. 3 shows one example of the manner in which the gain signal G can vary from a value equal to zero for low values of the ratio signal R to a value equal to one for high values of the ratio signal R.
  • Other functions can be used in the lookup table.
  • FIG. 4 provides a flowchart for a method practiced by the system 10 of FIG. 1.
  • tissue and contrast agent are insonified with ultrasound transmit signals centered at the fundamental frequency f.
  • fundamental receive signals centered at the fundamental frequency f and harmonic receive signals centered at the harmonic frequency 2f are acquired.
  • the ratio signal R is formed as a function of the fundamental and harmonic receive signals.
  • the image signal I is generated as a function of the ratio signal R.
  • regions with high concentration of contrast agent will have a relatively smaller difference between the integrated fundamental and harmonic receive signals F, H, while the regions with a small concentration of contrast agent will have a larger difference between the integrated backscattered signal levels in these two bandpass.
  • the image signal I can be made more specific for contrast agent or for tissue harmonics, as desired.
  • FIG. 5 relates to another embodiment of this invention, which can be implemented using a modified version of the system 10 of FIG. 2.
  • the transmitter is of the type capable of transmitting pulses at two different transmit center frequencies f 1 and f 2 on successive firing events.
  • f 1 can be equal to the fundamental frequency f described above
  • f 2 can be equal to the harmonic frequency 2f described above.
  • the other modification to the system 10 is that the ratio former 26 is designed to take a ratio of receive signals from two separate firing events filtered by the filter 24 at 2f. This can be accomplished in any suitable way, as for example by providing a buffer in the ratio former.
  • step 70 tissue and contrast agent are insonified with ultrasonic transmit signals centered at frequency f 1 .
  • step 72 harmonic receive signals R1 are acquired in response to the transmit event of step 70. These harmonic receive signals R1 emphasize harmonic components in a bandpass centered at 2f 1 , as passed by the filter 24 of FIG. 2.
  • step 74 the tissue and contrast agent are insonified with ultrasound transmit signals in a second transmit event centered at ultrasonic frequency 2f 1 .
  • step 76 fundamental receive signals R2 are acquired from the second transmit event of step 74. These fundamental receive signals R2 emphasize the same frequency band as that of the receive signals R1, namely 2f 1 .
  • step 78 a ratio signal is formed as a function of R1 and R2, and in step 80 an image signal is generated as a function of the ratio signal of step 78.
  • the formation of the ratio signal in step 78 and the generation of the image signal in step 80 can employ any of the techniques described elsewhere in this specification.
  • the method of FIG. 5 uses a ratio of receive signals in the same bandpass acquired in response to two separate transmit events. These two separate transmit events are centered at the fundamental transmit frequency f 1 and the second harmonic of this frequency 2f 1 , respectively.
  • the embodiment of FIG. 5 is based on the recognition that the second harmonic response of tissue differs in strength from the second harmonic response of contrast agent, and that the ratio signal of FIG. 5 thereby provides an indication of the presence or absence of contrast agent.
  • any suitable technology can be used for the transmitter 12 and the receiver 18, including both analog and digital systems.
  • any suitable transducer probe 16 can be used, including phased arrays, steered transducers, 1-, 1.5- and 2-dimensional arrays, and flat as well as curved arrays.
  • the filters 22, 24 can be implemented as desired.
  • the filter 22 can be a lowpass filter and the filter 24 can be a highpass filter, both of which can be considered examples of bandpass filters.
  • Both digital and analog filtering techniques can be used.
  • Demodulation techniques can be used to achieve desired filtering, and if desired the filters 22, 24 can be placed upstream of the receive beamformer in the signal path.
  • the ratio former 26 can be implemented with a dividing circuit if the detectors 28, 30 are linear detectors as opposed to log detectors.
  • Various elements, including the filters and the integrator for example, can be time-shared between the two processing channels.
  • the image signal generator 38 does not in all embodiments include a lookup table. Rather, the image signal can be modulated as a linear function of the ratio signal and one of the receive signals, or the ratio signal can be thresholded or applied to an analog or digital function generating circuit. As another alternative, the image signal generator can generate the image signal as a function only of the ratio signal, and not use the ratio signal to modulate another receive signal. Also, the ratio signal can be used as an input to a color image, and the resulting color signal can be displayed either alone or overlaid on a gray-scale image. As before, either analog or digital techniques can be used.
  • the image processor 46 and the display 48 can take any suitable form, and many variants are known to those skilled in the art.
  • the image processor can include a B-mode processor and/or an M-mode processor.
  • the ratio signal may be displayed as an M-mode strip, either in an overlay mode or by itself. In the overlay mode the ratio signal can be mixed with underlying image data.
  • the detectors are of the type that perform a squaring operation.
  • the integration frequency range is determined by the signal bandwidth in the fundamental and harmonic frequency bandpass.
  • the filters can be placed before the beamformer when signal processing is performed in the frequency domain.
  • the system and method described above provide the advantage of improved differentiation of contrast agent from tissue. Because of nonlinear propagation, the presence of signal at the second harmonic frequency does not uniquely identify the presence of contrast and therefore tissue with strong linear scattering can appear brighter in a harmonic image than contrast agent. By modifying the gain based on additional frequency components, the brightness of strong linear scatterers can be reduced and the differentiation of contrast agent from tissue can be improved.
  • differential sub-band techniques can be applied to other techniques developed for contrast agent imaging.
  • the method described above can be used to improve color Doppler imaging, destruction-based imaging techniques, and pulse inversion imaging techniques. With all these techniques, processing information from several bands simultaneously can provide additional information about the character of the scatterers involved that improves the ability to discriminate between contrast agent and tissue.
  • the methods described above can be used in conjunction with two-pulse techniques, Doppler techniques, or other harmonic image processing techniques to enhance these methods.
  • the image formed using the methods described above can be constructed over multiple frames including one or more frames designed to destroy contrast agent.
  • the rate at which the ratio of the harmonic over the fundamental signal varies is measured in order to establish a perfusion rate.
  • harmonic is intended broadly to encompass both harmonics and subharmonics.
  • signal is intended broadly to encompass both analog and digital signals.
  • one signal is said to vary as a function of another whether or not the first signal also varies as a function of additional, unnamed parameters.
  • an image signal is said to vary as a function of a gain signal, whether or not the image signal also varies as a function of the fundamental or harmonic receive signals F, H.
  • the term "function of” is intended broadly to encompass systems that include integrators, as well as systems that do not.
  • a ratio signal is said to vary as a function of a ratio of fundamental and harmonic receive signal, whether or not integrators are placed in the signal path.
  • fundamental receive signal is intended to encompass a receive signal that has been filtered or demodulated to emphasize the fundamental component of the receive signal, recognizing that a fundamental receive signal will have a bandwidth that extends to some extent beyond the fundamental frequency itself.
  • harmonic receive signal is intended to refer to a receive signal that has been filtered or demodulated to emphasize harmonic spectral components, and a harmonic receive signal will also have a bandwidth that extends around the specific harmonic frequency.
  • ratio is intended broadly both to encompass the quotient of two linear signals and the difference of two logarithmic signals.
  • ratio is also intended to encompass regularized ratios. For example, one type of regularized ratio that is commonly used to reduce artifacts adds a non-zero constant to the denominator of the ratio in order to prevent the denominator from ever being set equal to zero.
  • the ratio F/H can be set equal to F/(H+X), where X is a non-zero constant.

Abstract

A medical diagnostic ultrasound imaging system and method modulate the image signal as a function of the ratio of the harmonic receive signal to the fundamental receive signal. Tissue harmonic backscatter differs substantially in spectral shape as compared to contrast agent backscatter, and this method allows improved discrimination between contrast agent and tissue backscatter.

Description

BACKGROUND
This invention relates to medical diagnostic ultrasound imaging, and in particular to methods and systems for distinguishing contrast agent from surrounding tissue.
Contrast agent imaging is an important medical diagnostic ultrasound imaging mode. One limitation in many contrast imaging systems is the difficulty of distinguishing echo signals from contrast agents from echo signals from surrounding tissue. This is because both contrast agents and tissue generate nonlinear return signals at frequencies other than the frequency of the insonifying signals. It is well recognized that nonlinear signals from tissue are generated by nonlinear propagation of the insonifying ultrasound wave.
Prior-art harmonic imaging methods that exploit non-linear behaviour of contrast agent include B-mode harmonic imaging, B-mode harmonic pulse inversion imaging, harmonic power Doppler imaging, and color harmonic pulse inversion imaging.
In B-mode harmonic imaging, the signal is transmitted at a fundamental frequency f, and the receive signal is filtered to emphasize frequency components near the second harmonic, 2f. Contrast agents are known to have a stronger second harmonic response than tissue, and for this reason the receive signal from contrast agent is enhanced over that from tissue. This method is currently used by most major manufacturers of ultrasound imaging equipment. Specific examples are described in Mine U.S. Pat. No. 5,724,976, Uhlendorf U.S. Pat. No. 5,410,516, Schutt U.S. Pat. No. 5,733,527, and "Simulated Capillary Blood Measurement Using a Nonlinear Ultrasonic Contrast Agent," Schrope et al.; Ultrasonic Imaging, Vol. 14, pp. 134-158, 1992.
In B-mode harmonic pulse inversion imaging, two pulses are transmitted along the same ultrasound line, where one pulse is shifted by 180° with respect to the other. Receive signals from the two pulses are then summed, and the resultant signal is displayed. The pulse inversion technique cancels stationary fundamental frequency signals and retains second harmonic signals as well as some nonstationary fundamental signals. Examples of such methods include Chapman U.S. Pat. No. 5,632,277 and Hwang U.S. Pat. No. 5,706,819.
In harmonic power Doppler imaging, multiple pulses are transmitted along the same ultrasound line, and the receive signal is filtered about the second harmonic frequency. The resultant signals are then filtered with a high-pass filter to remove stationary signals. Examples of such methods are described in Averkiou U.S. Pat. No. 5,833,613, Johnson U.S. Pat. No. 5,456,257, and "Harmonic Power Mode Doppler Using Microbubble Contrast Agents: An Improved Method for Small Vessel Flow Imaging," Burns, et al., 1994 IEEE Ultrasonic Symposium, pp. 1547-1550, 1994.
In color harmonic pulse inversion imaging, the approach used is similar to that used in harmonic power Doppler imaging described above, but alternate pulses are shifted in phase by 180°. Such methods are described in "Pulse Inversion Doppler: A New Method for Detecting Nonlinear Echoes from Microbubble Contrast Agents," Simpson and Burns, 1997 IEEE Ultrasonic Symposium, 1997.
Though these techniques succeed in presenting contrast agent enhanced images, they do not entirely meet the needs of clinicians in the field. B-mode harmonic imaging does not differentiate between second harmonic signals generated by nonlinear propagation through tissue and second harmonic signals generated by contrast agents. Pulse inversion imaging further enhances second harmonic signals, but it still does not differentiate between second harmonic signals generated by nonlinear propagation through tissue and second harmonic signals generated by contrast agents. Harmonic power Doppler imaging attempts to differentiate contrast agent from tissue by looking for a loss of correlation between successive pulses due to agent destruction, agent motion or other methods. However, tissue motion will also result in a loss of correlation and may appear as a displayed signal. Tissue motion is reduced currently by a combination of increasing the pulse repetition frequency and/or the use of more aggressive clutter filters. However, increasing the pulse repetition frequency may reduce the signal from destroyed contrast agent, and more aggressive clutter filters may filter out contrast signal as well as tissue signal. These problems and drawbacks apply also to color harmonic pulse inversion techniques.
It would be advantageous to improve the specificity for imaging contrast agent as opposed to tissue. Increased specificity would allow contrast agent to be used for sensitive qualitative and quantitative measurements of blood flow in tissue.
SUMMARY
The present inventors have discovered specific spectral properties of the backscattered fundamental and harmonic receive signals that can be used to distinguish contrast agent from tissue in a medical diagnostic ultrasound imaging system. One preferred embodiment described below uses the ratio of the energy in the harmonic sub-band of the receive signal to the energy in the fundamental sub-band of the receive signal to detect the presence or absence of contrast agent. This ratio is then used to modulate the image signal to enhance the display of contrast agent.
The foregoing introductory paragraph is intended to present one aspect of the preferred embodiment described below by way of introduction. Nothing in this paragraph should be taken as a limitation on the following claims.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a frequency versus magnitude plot of receive signals associated with tissue and contrast agent.
FIG. 2 is a block diagram showing a medical diagnostic ultrasound imaging system that incorporates a presently preferred embodiment to this invention.
FIG. 3 is a graph showing one example of the manner in which the gain signal can vary as a function of the ratio signal in the system of FIG. 3.
FIG. 4 is a flow chart of a preferred embodiment of the method of this invention, as practiced by the system of FIG. 1.
FIG. 5 in a flow chart is a second preferred embodiment of this invention.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
General Discussion
Though both contrast agent and tissue produce harmonic signals in response to fundamental transmit signals, the level of the harmonic signal with respect to the fundamental signal in the backscattered receive signal is significantly different for contrast agent as compared to tissue. In the preferred embodiments that follow, this property is used to discriminate between contrast agent and tissue.
FIG. 1 shows a receive signal RT associated with backscattered signals from tissue insonified at a transmit center frequency f. In RT the fundamental frequency f is 1.7 MHz and the harmonic frequency 2f is 3.4 MHz. Note that the receive signal RT peaks at the fundamental frequency f and has a magnitude at the harmonic frequency 2f that is suppressed by about 15 dB as compared to the magnitude at the fundamental frequency f.
FIG. 1 also shows a frequency versus magnitude plot of a receive signal RC associated with backscattered echo signals from contrast agent. As before, the receive signal RC is generated with a transmit signal centered at a fundamental frequency f of 1.7 MHz. The receive signal RC also shows a peak magnitude at about the fundamental frequency f, and this peak magnitude is approximately 10 dB greater than the magnitude at the harmonic frequency 2f.
FIG. 1 shows that the ratio of amplitude near the fundamental frequency f to magnitude near the harmonic frequency 2f is substantially different for echo signals from tissue and contrast agent. As explained below, this ratio can be used to discriminate echo signals from tissue versus echo signals from contrast agent.
Specific Embodiments
Turning now to FIG. 2, an ultrasonic imaging system 10 includes a transmitter 12 that generates transmit signals T. The transmit signals T are applied via a transmit/receive switch 14 to a transducer probe 16. The transducer probe 16 responds to the transmit signals T by generating ultrasonic energy at the frequency of the transmit signal T and directing this ultrasonic energy into an imaged region that includes tissue T and contrast agent C.
The tissue T and the contrast agent C backscatter ultrasonic energy from the imaged region to the transducer probe 16. This backscattered ultrasonic energy causes the transducer probe 16 to generate receive signals RC and RT that are applied via the transmit/receive switch 14 to a receiver 18.
In this embodiment, the receiver 18 includes a receive beamformer 20 and two filters 22, 24. The first filter 22 is a bandpass filter that passes a sub-band of the receive signals centered around the fundamental frequency f. As before, the transmit signals T are also centered around the fundamental frequency f. The second filter 24 is another bandpass filter having a bandpass centered around the harmonic frequency 2f. The output signals generated by the filters 22, 24 are fundamental receive signals F and harmonic receive signals H, respectively.
The fundamental and harmonic receive signals F, H are applied to a ratio former 26 that in this embodiment includes first and second detectors 28, 30. The detectors 28, 30 preferably square the respective input signals and take the log of the squared result. The detected output signals generated by the detectors 28, 30 are applied to integrators 32, 34 that integrate the detected output signals over frequency. The output signals generated by the integrators 32, 34 are then subtracted from one another in a subtraction circuit 36 that in this application performs a division function. In this embodiment, the output signals of the integrators 32, 34 vary as the logarithm of the fundamental and harmonic receive signals, after detection and integration. Thus, the difference between the output signals of the integrators 32, 34 creates a ratio signal R that varies as a function of the ratio of the detected, integrated fundamental and harmonic receive signals F, H.
The ratio signal R is applied to an image signal generator 38 that in this embodiment includes a lookup table 40. The lookup table 40 generates a gain signal G in response to the ratio signal R, and this gain signal G is applied to a multiplier 42. The multiplier 42 is also responsive to one of the fundamental and harmonic receive signals F, H, as determined by the state of a switch 44. The output of the multiplier 42 is an image signal I that is applied to an image processor 46. The output of the image processor 46 is applied to a display 48.
The lookup table 40 allows the gain applied to the fundamental or harmonic receive signal F, H to be varied as a linear or nonlinear function of the ratio signal R. FIG. 3 shows one example of the manner in which the gain signal G can vary from a value equal to zero for low values of the ratio signal R to a value equal to one for high values of the ratio signal R. Other functions can be used in the lookup table.
FIG. 4 provides a flowchart for a method practiced by the system 10 of FIG. 1. In step 60, tissue and contrast agent are insonified with ultrasound transmit signals centered at the fundamental frequency f. In step 62, fundamental receive signals centered at the fundamental frequency f and harmonic receive signals centered at the harmonic frequency 2f are acquired. In step 64, the ratio signal R is formed as a function of the fundamental and harmonic receive signals. In step 66, the image signal I is generated as a function of the ratio signal R.
As described above, regions with high concentration of contrast agent will have a relatively smaller difference between the integrated fundamental and harmonic receive signals F, H, while the regions with a small concentration of contrast agent will have a larger difference between the integrated backscattered signal levels in these two bandpass. By modulating either the fundamental or harmonic receive signal with a gain signal that varies as a function of this ratio, the image signal I can be made more specific for contrast agent or for tissue harmonics, as desired.
FIG. 5 relates to another embodiment of this invention, which can be implemented using a modified version of the system 10 of FIG. 2. In this modified version the transmitter is of the type capable of transmitting pulses at two different transmit center frequencies f1 and f2 on successive firing events. For example, f1 can be equal to the fundamental frequency f described above, and f2 can be equal to the harmonic frequency 2f described above. The other modification to the system 10 is that the ratio former 26 is designed to take a ratio of receive signals from two separate firing events filtered by the filter 24 at 2f. This can be accomplished in any suitable way, as for example by providing a buffer in the ratio former.
Returning to FIG. 5, in step 70 tissue and contrast agent are insonified with ultrasonic transmit signals centered at frequency f1. In step 72 harmonic receive signals R1 are acquired in response to the transmit event of step 70. These harmonic receive signals R1 emphasize harmonic components in a bandpass centered at 2f1, as passed by the filter 24 of FIG. 2. In step 74 the tissue and contrast agent are insonified with ultrasound transmit signals in a second transmit event centered at ultrasonic frequency 2f1. In step 76 fundamental receive signals R2 are acquired from the second transmit event of step 74. These fundamental receive signals R2 emphasize the same frequency band as that of the receive signals R1, namely 2f1. In step 78 a ratio signal is formed as a function of R1 and R2, and in step 80 an image signal is generated as a function of the ratio signal of step 78. The formation of the ratio signal in step 78 and the generation of the image signal in step 80 can employ any of the techniques described elsewhere in this specification.
The method of FIG. 5 uses a ratio of receive signals in the same bandpass acquired in response to two separate transmit events. These two separate transmit events are centered at the fundamental transmit frequency f1 and the second harmonic of this frequency 2f1, respectively. The embodiment of FIG. 5 is based on the recognition that the second harmonic response of tissue differs in strength from the second harmonic response of contrast agent, and that the ratio signal of FIG. 5 thereby provides an indication of the presence or absence of contrast agent.
Of course, many modifications and alterations can be made to the preferred embodiments described above. For example, any suitable technology can be used for the transmitter 12 and the receiver 18, including both analog and digital systems. Similarly, any suitable transducer probe 16 can be used, including phased arrays, steered transducers, 1-, 1.5- and 2-dimensional arrays, and flat as well as curved arrays.
The filters 22, 24 can be implemented as desired. For example, the filter 22 can be a lowpass filter and the filter 24 can be a highpass filter, both of which can be considered examples of bandpass filters. Both digital and analog filtering techniques can be used. Demodulation techniques can be used to achieve desired filtering, and if desired the filters 22, 24 can be placed upstream of the receive beamformer in the signal path.
Any suitable detector and integrator can be used in the ratio former. For example, the ratio former 26 can be implemented with a dividing circuit if the detectors 28, 30 are linear detectors as opposed to log detectors. Various elements, including the filters and the integrator for example, can be time-shared between the two processing channels.
The image signal generator 38 does not in all embodiments include a lookup table. Rather, the image signal can be modulated as a linear function of the ratio signal and one of the receive signals, or the ratio signal can be thresholded or applied to an analog or digital function generating circuit. As another alternative, the image signal generator can generate the image signal as a function only of the ratio signal, and not use the ratio signal to modulate another receive signal. Also, the ratio signal can be used as an input to a color image, and the resulting color signal can be displayed either alone or overlaid on a gray-scale image. As before, either analog or digital techniques can be used.
The image processor 46 and the display 48 can take any suitable form, and many variants are known to those skilled in the art. For example, the image processor can include a B-mode processor and/or an M-mode processor. The ratio signal may be displayed as an M-mode strip, either in an overlay mode or by itself. In the overlay mode the ratio signal can be mixed with underlying image data.
Using Parseval's theorem, it can be shown that signal processing can be performed in the time domain instead of the frequency domain to achieve results similar to those described above. The method described above compares the energy of the backscattered signal in two different frequency bands. Ideally, the detectors are of the type that perform a squaring operation. In the frequency domain, the integration frequency range is determined by the signal bandwidth in the fundamental and harmonic frequency bandpass. As another alternative, the filters can be placed before the beamformer when signal processing is performed in the frequency domain.
The system and method described above provide the advantage of improved differentiation of contrast agent from tissue. Because of nonlinear propagation, the presence of signal at the second harmonic frequency does not uniquely identify the presence of contrast and therefore tissue with strong linear scattering can appear brighter in a harmonic image than contrast agent. By modifying the gain based on additional frequency components, the brightness of strong linear scatterers can be reduced and the differentiation of contrast agent from tissue can be improved.
In addition to being beneficial in single pulse harmonic imaging, differential sub-band techniques can be applied to other techniques developed for contrast agent imaging. Specifically, the method described above can be used to improve color Doppler imaging, destruction-based imaging techniques, and pulse inversion imaging techniques. With all these techniques, processing information from several bands simultaneously can provide additional information about the character of the scatterers involved that improves the ability to discriminate between contrast agent and tissue.
The methods described above can be used in conjunction with two-pulse techniques, Doppler techniques, or other harmonic image processing techniques to enhance these methods.
The image formed using the methods described above can be constructed over multiple frames including one or more frames designed to destroy contrast agent. The rate at which the ratio of the harmonic over the fundamental signal varies is measured in order to establish a perfusion rate.
As used herein, the term "harmonic" is intended broadly to encompass both harmonics and subharmonics. The term "signal" is intended broadly to encompass both analog and digital signals.
As used herein, one signal is said to vary as a function of another whether or not the first signal also varies as a function of additional, unnamed parameters. Thus, an image signal is said to vary as a function of a gain signal, whether or not the image signal also varies as a function of the fundamental or harmonic receive signals F, H. Also, the term "function of" is intended broadly to encompass systems that include integrators, as well as systems that do not. Thus, a ratio signal is said to vary as a function of a ratio of fundamental and harmonic receive signal, whether or not integrators are placed in the signal path.
The term "fundamental receive signal" is intended to encompass a receive signal that has been filtered or demodulated to emphasize the fundamental component of the receive signal, recognizing that a fundamental receive signal will have a bandwidth that extends to some extent beyond the fundamental frequency itself. Similarly, the term "harmonic receive signal" is intended to refer to a receive signal that has been filtered or demodulated to emphasize harmonic spectral components, and a harmonic receive signal will also have a bandwidth that extends around the specific harmonic frequency.
The term "ratio" is intended broadly both to encompass the quotient of two linear signals and the difference of two logarithmic signals. The term "ratio" is also intended to encompass regularized ratios. For example, one type of regularized ratio that is commonly used to reduce artifacts adds a non-zero constant to the denominator of the ratio in order to prevent the denominator from ever being set equal to zero. In one example, the ratio F/H can be set equal to F/(H+X), where X is a non-zero constant.
The foregoing detailed description has described only a few of the many forms that this invention can take. For this reason, this detailed description is intended by way of illustration and not by way of limitation. It is only the following claims, including all equivalents, that are intended to define the scope of this invention.

Claims (32)

What is claimed is:
1. A medical diagnostic ultrasound imaging method comprising:
(a) insonifying a tissue containing a contrast agent with ultrasonic transmit signals;
(b) acquiring backscattered ultrasonic first receive signals in a first bandpass;
(c) acquiring backscattered ultrasonic second receive signals in a second bandpass;
(d) forming a gain signal as a function of the first and second receive signals; and
(e) applying the gain signal to backscattered ultrasonic receive signals.
2. A medical diagnostic ultrasound image processing method responsive to a first receive signal and a second receive signal, said first and second receive signals characterized by separate respective first and second bandpasses, said method comprising:
forming a gain signal as a function of the first and second receive signals; and
applying the gain signal to backscattered ultrasonic receive signals.
3. The method of claim 1 or 2, wherein the gain signal forming act comprises
forming a ratio signal as a function of the first and second receive signals; and
applying the ratio signal as an input to a lookup table that generates the gain signal in response thereto.
4. The method of claim 3, wherein the gain signal applying act comprises
applying the gain signal and at least one of the first and second receive signals to a multiplier that generates an image signal in response thereto.
5. The method of claim 3 wherein the lookup table generates the gain signal as an non-linear function of the ratio signal.
6. The method of claim 1 or 2, wherein the gain signal forming act step comprises
detecting and log compressing the first receive signal;
integrating the detected and log compressed first receive signal;
detecting and log compressing the second receive signal;
integrating the detected and log compressed second receive signal; and
applying the integrated, detected, log compressed first receive signal and the integrated, detected, log compressed second receive signal to a ratio former which generates a ratio signal.
7. The method of claim 1 wherein (b) comprises applying echo signals from a first transmit event to a first filter, and wherein (c) comprises applying echo signals from the first transmit event to a second filter.
8. The method of claim 7 wherein the transmit signals of the first transmit event are characterized by a fundamental frequency, wherein the first filter selectively passes fundamental components of the echo signals, and wherein the second filter selectively passes harmonic components of the echo signals.
9. A medical diagnostic ultrasound imaging method comprising:
(a) insonifying a tissue containing a contrast agent with ultrasonic transmit signals;
(b) acquiring backscattered ultrasonic first receive signals in a first bandpass;
(c) acquiring backscattered ultrasonic second receive signals in a second bandpass;
(d) forming a ratio signal as a function of the first and second receive signals; and
(e) generating an image signal as a function of the ratio signal;
wherein (a) comprises
(a1) insonifying the tissue containing the contrast agent with a first transmit event characterized by ultrasonic transmit signals having a first center frequency f1 ; and
(a2) insonifying the tissue containing the contrast agent with a second transmit event characterized by ultrasonic transmit signals having a second center frequency f2.
10. The method of claim 9 wherein the first and second bandpasses are characterized by first and second filter center frequencies, respectively, wherein the first and second filter center frequencies are both substantially equal to f2, and wherein f2 is about twice f1.
11. The method of claim 9 wherein the first and second bandpasses are characterized by substantially equal filter center frequencies.
12. The method of claim 1 wherein (e) generates an image signal selected from the group consisting of B-mode, M-mode, color flow image signals, and combinations thereof.
13. A medical diagnostic ultrasound imaging system comprising:
means for insonifying a tissue containing a contrast agent with ultrasonic transmit signals;
first means for acquiring backscattered ultrasonic first receive signals in a first bandpass;
second means for acquiring backscattered ultrasonic second receive signals in a second bandpass;
means for forming a gain signal as a function of the first and second receive signals; and
means for applying the gain signal to backscattered ultrasonic receive signals.
14. The invention of claim 13, wherein the gain signal forming means comprises
a ratio former responsive to the first and second receive signals; and
a lookup table responsive to the ratio former to generate the gain signal.
15. The invention of claim 14, wherein the gain signal applying means comprises a multiplier responsive to the gain signal and to at least one of the first and second receive signals to generate an image signal.
16. The invention of claim 13, wherein the gain signal forming means comprises:
a first log detector responsive to the first receive signal;
a first integrator responsive to the first detector;
a second log detector responsive to the second receive signal;
a second integrator responsive to the second detector; and
an element responsive to the first and second integrators to generate a ratio signal.
17. The invention of claim 13, wherein the first and second means comprise first and second filters, respectively, and wherein the first and second means are operative to acquire the first and second receive beams, respectively, from a common transmit event.
18. The invention of claim 17, wherein the transmit signals are characterized by a fundamental frequency, wherein the first bandpass emphasizes fundamental components of the first receive signals, and wherein the second bandpass emphasizes harmonic components of the second receive signals.
19. A medical diagnostic ultrasound imaging system comprising:
means for insonifying a tissue containing a contrast agent with ultrasonic transmit signals;
first means for acquiring backscattered ultrasonic first receive signals in a first bandpass;
second means for acquiring backscattered ultrasonic second receive signals in a second bandpass;
means for forming a ratio signal as a function of the first and second receive signals; and
means for generating an image signal as a function of the ratio signal;
wherein the insonifying means is operative to insonify the tissue containing the contrast agent with (1) first transmit signals from a first transmit event having a first center frequency f1 and (2) second transmit signals from a second transmit event having a second center frequency f2.
20. The invention of claim 19, wherein the first and second bandpasses both selectively pass receive signals at the second center frequency f2, and wherein f2 is about twice f1.
21. The invention of claim 19, wherein the first and second bandpasses are characterized by substantially equal center frequencies.
22. The invention of claim 13 wherein the gain signal applying means generates an image signal selected from the group consisting of B-mode, M-mode, color flow image signals and combination thereof.
23. A medical diagnostic ultrasound image processor responsive to a first receive signal and a second receive signal, said first and second receive signals characterized by separate respective first and second bandpasses, said image processor comprising:
a gain signal generator operative to form a gain signal as a function of the first receive signal and the second receive signal; and
an image signal generator operative to apply the gain signal to backscattered ultrasonic receive signals.
24. The invention of claim 23, wherein the gain signal generator comprises:
a ratio former operative to form a ratio signal as a function of the first and second receive signals; and
a lookup table responsive to the ratio signal to generate the gain signal.
25. The invention of claim 24, wherein the image signal generator further comprises:
a multiplier responsive to the gain signal and to at least one of the first and second receive signals to generate an image signal.
26. The invention of claim 14 or 24 wherein the lookup table generates the gain signal as an non-linear function of the ratio signal.
27. The invention of claim 23, wherein the gain signal generator comprises:
a first log detector responsive to the first receive signal;
a first integrator responsive to the first detector;
a second log detector responsive to the second receive signal;
a second integrator responsive to the second detector; and
an element responsive to the first and second integrators to generate a ratio signal.
28. The invention of claim 23, wherein the first and second bandpasses are centered at respective center frequencies f1 and f2, and wherein f2 is substantially equal to two times f1.
29. A medical diagnostic ultrasound image processor responsive to a first receive signal and a second receive signal, said first and second receive signals characterized by separate respective first and second bandpasses, said image processor comprising:
a ratio former operative to form a ratio signal as a function of the first receive signal and the second receive signal; and
an image signal generator operative to generate an image signal as a function of the ratio signal;
wherein the first and second bandpasses are centered at substantially equal center frequencies, and wherein the first and second receive signals are associated with separate respective transmit events.
30. The invention of claim 23 wherein the image signal generator generates an image signal selected from the group consisting of B-mode, M-mode, color flow image signals and combination thereof.
31. The invention of claim 13 wherein the gain signal applying means applies the gain signal to fundamental ultrasonic receive signals prior to display.
32. The method of claim 1 wherein (e) comprises applying the gain signal to fundamental ultrasonic receive signals prior to display.
US09/282,402 1999-03-31 1999-03-31 Medical diagnostic ultrasonic imaging system and method using differential sub-band detection techniques Expired - Lifetime US6132377A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US09/282,402 US6132377A (en) 1999-03-31 1999-03-31 Medical diagnostic ultrasonic imaging system and method using differential sub-band detection techniques

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US09/282,402 US6132377A (en) 1999-03-31 1999-03-31 Medical diagnostic ultrasonic imaging system and method using differential sub-band detection techniques

Publications (1)

Publication Number Publication Date
US6132377A true US6132377A (en) 2000-10-17

Family

ID=23081363

Family Applications (1)

Application Number Title Priority Date Filing Date
US09/282,402 Expired - Lifetime US6132377A (en) 1999-03-31 1999-03-31 Medical diagnostic ultrasonic imaging system and method using differential sub-band detection techniques

Country Status (1)

Country Link
US (1) US6132377A (en)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6423007B2 (en) * 1998-05-26 2002-07-23 Riverside Research Institute Ultrasonic systems and methods for contrast agent concentration measurement
US6438258B1 (en) * 1998-01-23 2002-08-20 Koninklijke Philips Electronics N.V. Ultrasound image processing employing non-linear tissue response backscatter signals
US6458083B1 (en) * 1996-11-26 2002-10-01 Koninklijke Philips Electronics N.V. Ultrasonic harmonic imaging with adaptive image formation
US6461303B2 (en) * 2000-01-19 2002-10-08 Bjorn Angelsen Method of detecting ultrasound contrast agent in soft tissue, and quantitating blood perfusion through regions of tissue
US20030004414A1 (en) * 2001-05-31 2003-01-02 Mclaughlin Glen System and method for phase inversion ultrasonic imaging
US6508767B2 (en) 2000-02-16 2003-01-21 Koninklijke Philips Electronics N.V. Ultrasonic harmonic image segmentation
US20030081820A1 (en) * 2001-11-01 2003-05-01 Avinash Gopal B. Method for contrast matching of multiple images of the same object or scene to a common reference image
US20030092992A1 (en) * 2001-11-08 2003-05-15 Kabushiki Kaisha Toshiba Ultrasound diagnosis apparatus for imaging with a contrast agent
US20030176792A1 (en) * 2002-02-08 2003-09-18 Kabushiki Kaisha Toshiba Apparatus and method for ultrasonic diagnostic imaging
US6626836B2 (en) * 2001-04-04 2003-09-30 Siemens Medical Solutions Usa, Inc. Adaptive signal processing scheme for contrast agent imaging
US6638230B2 (en) * 2001-07-31 2003-10-28 Koninklijke Philips Electronics N.V. Apparatus and method of frequency compounding to perform contrast imaging
US20040030251A1 (en) * 2002-05-10 2004-02-12 Ebbini Emad S. Ultrasound imaging system and method using non-linear post-beamforming filter
US20040059218A1 (en) * 2001-02-01 2004-03-25 Hiroshi Kanda Ultrasonic enhanced-contrast imager and its method
EP1411368A1 (en) * 2002-10-15 2004-04-21 Matsushita Electric Industrial Co., Ltd. Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
US20050054928A1 (en) * 2003-09-09 2005-03-10 Marino Cerofolini Ultrasound imaging method combined with the presence of contrast media in the body under examination
US20050087016A1 (en) * 2003-10-24 2005-04-28 Gilmore Robert S. Inspection method and apparatus for determining incipient mechanical failure
US20070016024A1 (en) * 2005-06-28 2007-01-18 Siemens Medical Solutions Usa, Inc. Ultrasound imaging system having motion adaptive gain
US20070073146A1 (en) * 2005-09-27 2007-03-29 Siemens Medical Solutions Usa, Inc. Advanced characterization of contrast agents with ultrasound
US20070125174A1 (en) * 2005-12-01 2007-06-07 General Electric Company Ultrasonic inspection system and method
EP1884195A1 (en) * 2005-05-27 2008-02-06 Hitachi Medical Corporation Ultrasonograph and ultrasonic image display method
US8167803B2 (en) * 2007-05-16 2012-05-01 Verathon Inc. System and method for bladder detection using harmonic imaging
US8679018B2 (en) 1999-08-20 2014-03-25 Zonare Medical Systems, Inc. Broad-beam imaging
US20170071569A1 (en) * 2015-09-14 2017-03-16 Toshiba Medical Systems Corporation Ultrasonic diagnostic apparatus and signal processing apparatus

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5255683A (en) * 1991-12-30 1993-10-26 Sound Science Limited Partnership Methods of and systems for examining tissue perfusion using ultrasonic contrast agents
US5410516A (en) * 1988-09-01 1995-04-25 Schering Aktiengesellschaft Ultrasonic processes and circuits for performing them
US5456257A (en) * 1994-11-23 1995-10-10 Advanced Technology Laboratories, Inc. Ultrasonic detection of contrast agents
US5526816A (en) * 1994-09-22 1996-06-18 Bracco Research S.A. Ultrasonic spectral contrast imaging
US5632277A (en) * 1996-06-28 1997-05-27 Siemens Medical Systems, Inc. Ultrasound imaging system employing phase inversion subtraction to enhance the image
US5706819A (en) * 1995-10-10 1998-01-13 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic imaging with harmonic contrast agents
US5724976A (en) * 1994-12-28 1998-03-10 Kabushiki Kaisha Toshiba Ultrasound imaging preferable to ultrasound contrast echography
US5733527A (en) * 1994-09-28 1998-03-31 Alliance Pharmaceutical Corp. Methods for harmonic imaging with ultrasound
US5740128A (en) * 1995-03-02 1998-04-14 Acuson Corporation Ultrasonic harmonic imaging system and method
US5833613A (en) * 1996-09-27 1998-11-10 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic imaging with contrast agents
US5961464A (en) * 1998-09-16 1999-10-05 Hewlett-Packard Company Ultrasound contrast agent detection using spectral analysis from acoustic scan lines

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5410516A (en) * 1988-09-01 1995-04-25 Schering Aktiengesellschaft Ultrasonic processes and circuits for performing them
US5255683A (en) * 1991-12-30 1993-10-26 Sound Science Limited Partnership Methods of and systems for examining tissue perfusion using ultrasonic contrast agents
US5526816A (en) * 1994-09-22 1996-06-18 Bracco Research S.A. Ultrasonic spectral contrast imaging
US5733527A (en) * 1994-09-28 1998-03-31 Alliance Pharmaceutical Corp. Methods for harmonic imaging with ultrasound
US5456257A (en) * 1994-11-23 1995-10-10 Advanced Technology Laboratories, Inc. Ultrasonic detection of contrast agents
US5724976A (en) * 1994-12-28 1998-03-10 Kabushiki Kaisha Toshiba Ultrasound imaging preferable to ultrasound contrast echography
US5740128A (en) * 1995-03-02 1998-04-14 Acuson Corporation Ultrasonic harmonic imaging system and method
US5706819A (en) * 1995-10-10 1998-01-13 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic imaging with harmonic contrast agents
US5632277A (en) * 1996-06-28 1997-05-27 Siemens Medical Systems, Inc. Ultrasound imaging system employing phase inversion subtraction to enhance the image
US5833613A (en) * 1996-09-27 1998-11-10 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic imaging with contrast agents
US5961464A (en) * 1998-09-16 1999-10-05 Hewlett-Packard Company Ultrasound contrast agent detection using spectral analysis from acoustic scan lines

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
"Harmonic Power Mode Doppler Using Microbubble Contrast Agents: An Improved Method for Small Vessel Flow Imaging," Burns et al., 1994 IEEE Ultrasonic Symposium, pp. 1547-1550, 1994.
"Pulse Inversion Doppler: A New Method for Detecting Nonlinear Echoes from Microbubble Contrast Agents," Simpson and Burns, 1997 IEEE Ultrasonic Symposium, 1997.
"Simulated Capillary Blood Measurement Using a Nonlinear Ultrasonic Contrast Agent," Schrope et al.; Ultrasonic Imaging, vol. 14, pp. 134-158, 1992.
Harmonic Power Mode Doppler Using Microbubble Contrast Agents: An Improved Method for Small Vessel Flow Imaging, Burns et al., 1994 IEEE Ultrasonic Symposium, pp. 1547 1550, 1994. *
Pulse Inversion Doppler: A New Method for Detecting Nonlinear Echoes from Microbubble Contrast Agents, Simpson and Burns, 1997 IEEE Ultrasonic Symposium, 1997. *
Simulated Capillary Blood Measurement Using a Nonlinear Ultrasonic Contrast Agent, Schrope et al.; Ultrasonic Imaging, vol. 14, pp. 134 158, 1992. *

Cited By (48)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6458083B1 (en) * 1996-11-26 2002-10-01 Koninklijke Philips Electronics N.V. Ultrasonic harmonic imaging with adaptive image formation
US6438258B1 (en) * 1998-01-23 2002-08-20 Koninklijke Philips Electronics N.V. Ultrasound image processing employing non-linear tissue response backscatter signals
US6423007B2 (en) * 1998-05-26 2002-07-23 Riverside Research Institute Ultrasonic systems and methods for contrast agent concentration measurement
US8679018B2 (en) 1999-08-20 2014-03-25 Zonare Medical Systems, Inc. Broad-beam imaging
US8764661B2 (en) 1999-08-20 2014-07-01 Zonare Medical Systems, Inc. Echolocation data generation
US6461303B2 (en) * 2000-01-19 2002-10-08 Bjorn Angelsen Method of detecting ultrasound contrast agent in soft tissue, and quantitating blood perfusion through regions of tissue
US6508767B2 (en) 2000-02-16 2003-01-21 Koninklijke Philips Electronics N.V. Ultrasonic harmonic image segmentation
US7198601B2 (en) * 2001-02-01 2007-04-03 Hitachi Medical Corporation Ultrasonic contrast medium imaging apparatus and method
US7850612B2 (en) 2001-02-01 2010-12-14 Hitachi Medical Corporation Ultrasonic enhanced-contrast imager and method
US20040059218A1 (en) * 2001-02-01 2004-03-25 Hiroshi Kanda Ultrasonic enhanced-contrast imager and its method
US6626836B2 (en) * 2001-04-04 2003-09-30 Siemens Medical Solutions Usa, Inc. Adaptive signal processing scheme for contrast agent imaging
US20040158149A1 (en) * 2001-05-31 2004-08-12 Mclaughlin Glen System for phase inversion ultrasonic imaging
US7699781B2 (en) 2001-05-31 2010-04-20 Zonare Medical Systems, Inc. System for phase inversion ultrasonic imaging
US10222461B2 (en) 2001-05-31 2019-03-05 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Phase inversion ultrasonic imaging
US20030004414A1 (en) * 2001-05-31 2003-01-02 Mclaughlin Glen System and method for phase inversion ultrasonic imaging
US20080103394A1 (en) * 2001-05-31 2008-05-01 Mclaughlin Glen Phase Inversion Ultrasonic Imaging
US6638230B2 (en) * 2001-07-31 2003-10-28 Koninklijke Philips Electronics N.V. Apparatus and method of frequency compounding to perform contrast imaging
US20030128892A1 (en) * 2001-11-01 2003-07-10 Ge Medical Systems Technology Company Llc Method for contrast matching of multiple images of the same object or scene to a common reference image
US7194145B2 (en) * 2001-11-01 2007-03-20 General Electric Company Method for contrast matching of multiple images of the same object or scene to a common reference image
US20030081820A1 (en) * 2001-11-01 2003-05-01 Avinash Gopal B. Method for contrast matching of multiple images of the same object or scene to a common reference image
US7206460B2 (en) * 2001-11-01 2007-04-17 General Electric Company Method for contrast matching of multiple images of the same object or scene to a common reference image
US6726630B2 (en) * 2001-11-08 2004-04-27 Kabushiki Kaisha Toshiba Ultrasound diagnosis apparatus for imaging with a contrast agent
US20030092992A1 (en) * 2001-11-08 2003-05-15 Kabushiki Kaisha Toshiba Ultrasound diagnosis apparatus for imaging with a contrast agent
US20030176792A1 (en) * 2002-02-08 2003-09-18 Kabushiki Kaisha Toshiba Apparatus and method for ultrasonic diagnostic imaging
US7044914B2 (en) * 2002-02-08 2006-05-16 Kabushiki Kaisha Toshiba Apparatus and method for ultrasonic diagnostic imaging
US20040030251A1 (en) * 2002-05-10 2004-02-12 Ebbini Emad S. Ultrasound imaging system and method using non-linear post-beamforming filter
US6705993B2 (en) 2002-05-10 2004-03-16 Regents Of The University Of Minnesota Ultrasound imaging system and method using non-linear post-beamforming filter
US6951540B2 (en) 2002-05-10 2005-10-04 Regents Of The University Of Minnesota Ultrasound imaging system and method using non-linear post-beamforming filter
US6923766B2 (en) * 2002-10-15 2005-08-02 Matsushita Electric Industrial Co., Ltd Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
EP1411368A1 (en) * 2002-10-15 2004-04-21 Matsushita Electric Industrial Co., Ltd. Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
US20040077947A1 (en) * 2002-10-15 2004-04-22 Manabu Migita Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
US8047994B2 (en) 2003-09-09 2011-11-01 Esaote, S.P.A. Ultrasound imaging method combined with the presence of contrast media in the body under examination
US20050054928A1 (en) * 2003-09-09 2005-03-10 Marino Cerofolini Ultrasound imaging method combined with the presence of contrast media in the body under examination
US20100145196A1 (en) * 2003-09-09 2010-06-10 Marino Cerofolini Ultrasound imaging method combined with the presence of contrast media in the body under examination
US7691064B2 (en) * 2003-09-09 2010-04-06 Esaote, S.P.A. Ultrasound imaging method combined with the presence of contrast media in the body under examination
US20050087016A1 (en) * 2003-10-24 2005-04-28 Gilmore Robert S. Inspection method and apparatus for determining incipient mechanical failure
US6988410B2 (en) * 2003-10-24 2006-01-24 General Electric Company Inspection method and apparatus for determining incipient mechanical failure
EP1884195A4 (en) * 2005-05-27 2011-03-30 Hitachi Medical Corp Ultrasonograph and ultrasonic image display method
US20090299182A1 (en) * 2005-05-27 2009-12-03 Katsunori Asafusa Ultrasonic diagnostic apparatus and ultrasonic image display method
EP1884195A1 (en) * 2005-05-27 2008-02-06 Hitachi Medical Corporation Ultrasonograph and ultrasonic image display method
US8055036B2 (en) 2005-05-27 2011-11-08 Hitachi Medical Corporation Ultrasonic diagnostic apparatus and ultrasonic image display method
US7645236B2 (en) 2005-06-28 2010-01-12 Siemens Medical Solutions Usa, Inc. Ultrasound imaging system having motion adaptive gain
US20070016024A1 (en) * 2005-06-28 2007-01-18 Siemens Medical Solutions Usa, Inc. Ultrasound imaging system having motion adaptive gain
US20070073146A1 (en) * 2005-09-27 2007-03-29 Siemens Medical Solutions Usa, Inc. Advanced characterization of contrast agents with ultrasound
US20070125174A1 (en) * 2005-12-01 2007-06-07 General Electric Company Ultrasonic inspection system and method
US7546769B2 (en) 2005-12-01 2009-06-16 General Electric Compnay Ultrasonic inspection system and method
US8167803B2 (en) * 2007-05-16 2012-05-01 Verathon Inc. System and method for bladder detection using harmonic imaging
US20170071569A1 (en) * 2015-09-14 2017-03-16 Toshiba Medical Systems Corporation Ultrasonic diagnostic apparatus and signal processing apparatus

Similar Documents

Publication Publication Date Title
US6132377A (en) Medical diagnostic ultrasonic imaging system and method using differential sub-band detection techniques
US6213951B1 (en) Medical diagnostic ultrasound method and system for contrast specific frequency imaging
US5628322A (en) Method of ultrasound imaging and diagnostic ultrasound system
USRE43048E1 (en) Ultrasonic diagnostic imaging with harmonic contrast agents
JP5150267B2 (en) Ultrasound diagnostic imaging system to detect liver damage
US6508767B2 (en) Ultrasonic harmonic image segmentation
US6171246B1 (en) Realtime ultrasonic imaging of perfusion using ultrasonic contrast agents
JP3432204B2 (en) Ultrasound diagnostic equipment
US8454517B2 (en) Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
US5961464A (en) Ultrasound contrast agent detection using spectral analysis from acoustic scan lines
US7985183B2 (en) Ultrasonic diagnostic equipment and ultrasonic image generation method
US6626836B2 (en) Adaptive signal processing scheme for contrast agent imaging
US6814703B2 (en) Apparatus and method for ultrasonic diagnostic imaging using a contrast medium
CA2127127A1 (en) Ultrasound contrast agent examination of tissue perfusion
JPH09505761A (en) Contrast imaging of ultrasonic spectrum
JP2519777B2 (en) Ultrasonic device and method for performing ultrasonic measurement and ultrasonic inspection
JPH11128231A (en) Ultrasonic echographic system
JPH057588A (en) Ultrasonic doppler diagnostic apparatus
JPS62117535A (en) Ultrasonic doppler apparatus
JP3029315B2 (en) Ultrasound blood flow imaging device
JP4574790B2 (en) Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
US8189427B2 (en) Clutter signal filtering for doppler signal
JPH0224140B2 (en)
US20040127795A1 (en) Method of ultrasonic speckle reduction using wide-band frequency compounding with tissue-generated harmonics
US8388537B2 (en) Ultrasonic imaging apparatus

Legal Events

Date Code Title Description
AS Assignment

Owner name: ACUSON CORPORATION, CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:BOLORFOROSH, MIRSAID S.;GARDNER, EDWARD A.;HOLLEY, GREGORY L.;AND OTHERS;REEL/FRAME:010009/0836;SIGNING DATES FROM 19990511 TO 19990527

STCF Information on status: patent grant

Free format text: PATENTED CASE

FEPP Fee payment procedure

Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

FPAY Fee payment

Year of fee payment: 4

FPAY Fee payment

Year of fee payment: 8

AS Assignment

Owner name: SIEMENS MEDICAL SOLUTIONS USA, INC.,PENNSYLVANIA

Free format text: CHANGE OF NAME;ASSIGNOR:SIEMENS MEDICAL SYSTEMS, INC.;REEL/FRAME:024563/0051

Effective date: 20010801

AS Assignment

Owner name: SIEMENS MEDICAL SOLUTIONS USA, INC., PENNSYLVANIA

Free format text: RE-RECORD TO CORRECT CONVEYING PARTY NAME PREVIOUSLY RECORDED AT REEL 024563 FRAME 0051;ASSIGNORS:ACUSON CORPORATION;ACUSON LLC;ACUSON CORPORATION;SIGNING DATES FROM 20021218 TO 20050926;REEL/FRAME:024651/0673

FPAY Fee payment

Year of fee payment: 12