US4800891A - Doppler velocity processing method and apparatus - Google Patents

Doppler velocity processing method and apparatus Download PDF

Info

Publication number
US4800891A
US4800891A US07/120,577 US12057787A US4800891A US 4800891 A US4800891 A US 4800891A US 12057787 A US12057787 A US 12057787A US 4800891 A US4800891 A US 4800891A
Authority
US
United States
Prior art keywords
velocity
echo signals
estimation
groups
signals
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
US07/120,577
Other versions
US4830389A (en
Inventor
Jin H. Kim
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
Siemens Medical Laboratories Inc
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 Siemens Medical Laboratories Inc filed Critical Siemens Medical Laboratories Inc
Priority to US07/120,577 priority Critical patent/US4800891A/en
Assigned to SIEMENS MEDICAL LABORATORIES, INC. reassignment SIEMENS MEDICAL LABORATORIES, INC. ASSIGNMENT OF ASSIGNORS INTEREST. Assignors: KIM, JIN H.
Application granted granted Critical
Publication of US4800891A publication Critical patent/US4800891A/en
Assigned to SIEMENS QUANTUM, INC. reassignment SIEMENS QUANTUM, INC. ASSIGNMENT OF ASSIGNORS INTEREST. Assignors: SIEMENS MEDICAL LABORATORIES, INC.
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow

Definitions

  • This invention relates to velocity estimation using Doppler techniques, and more particularly to a spatial vector averaging method and apparatus for estimating blood flow velocities in a medical diagnostic apparatus using ultrasonic pulse signals.
  • Medical diagnostic ultrasound apparatus generate images of anatomical structures within a body of a patient by transmitting ultrahigh-frequency sound waves (typically on the order of 3.0 MHz) into the patient and then analyzing the echoes, i.e., ultrasonic signals reflected from the structure of the body being scanned.
  • ultrahigh-frequency sound waves typically on the order of 3.0 MHz
  • the most widely used ultrasound diagnostic apparatus displays in real-time anatomical information in the form of a two-dimensional image of a selected cross-section of the structure.
  • the ultrasound signals are swept across the structure in the form of a sector scan.
  • the sector scan is performed in real-time so that the image is available during the examination of the patient. In such cases, motion of the structure produces a corresponding moving image (i.e., a B-mode image).
  • anatomical defects can be relatively small and beyond the resolution capabilities of conventional anatomical ultrasound imaging.
  • a small anatomical defect may manifest itself by a pronounced change in blood flow velocity, for example aortic stenosis, mitral or aortic insufficiency or a congenital defect
  • a display of blood flow velocity would allow these abnormalities to be more easily detected.
  • One known method for velocity display is to use FFT techniques to process the echo signals reflected from a selected small volume so as to generate a numerical display. This method is severely limited by the fact that velocity is only determined for the small sample volume and is also not a two dimensional real-time image. Real-time imaging of velocity in a larger area is highly desireable.
  • velocity blood flow imaging has become an increasingly important portion of the ultrasound imaging device used in the medical diagnostic field, wherein a real-time blood flow image is superimposed over a real-time anatomical image.
  • European Patent Publication No. 0 100 094 by Namekawa of the Aloka Company describes an ultrasound blood flow imaging apparatus for two-dimensional display i.e., mapping, of blood flow velocities within a body.
  • the flow mapping is superimposed over a B-mode scan and displays blood flow in a color which is representative of the direction of the blood flow with respect to the ultrasonic transducer while variations in the intensity of the colors are representative of the blood flow velocities.
  • the blood flow velocity signal is generated using a pulsed-Doppler method.
  • Dopper velocity estimators use time-domain processing techniques wherein a series of N repetitive pulses of RF-ultrasonic signal separated by time periods of T are transmitted towards a moving target along a given scan direction. A relatively large number of echo signals (from a minimum of 10 to a maximum of 256) are processed using Doppler techniques for determining the velocity and turbulence of the moving object.
  • Doppler techniques for determining the velocity and turbulence of the moving object.
  • One such technique for determining Doppler frequency shift is shown, for example, by U.S. Pat. No. 4,542,657 issued to Barber et al. which uses I and Q sampled signals of a plurality (i.e., from 16 to 256) of the demodulated echo signals. FFT and zero-crossing velocity estimators are described in U.S. Pat. No.
  • time-domain processing technique uses autocorrelaton (pulse-pair) algorithms, such as known, for example, from an article by Kasai et al. published in the IEEE Transactions on Sonics and Ultrasonics, Volume SU-32, No. 3, May 1985.
  • pulse-pair pulse-pair
  • Each of the above time-domain processing techniques evaluate a relatively large number of return echoes (from 10 to 256) for estimating blood flow velocity and turbulence.
  • FIG. 1a shows groups (or pulses) of ultrasonic signal 1, 2, 3, . . . n which are excited and groups of return echo signals e 1 , e 2 , e 3 . . . e n which are received in response thereto.
  • the echoes are reflections from a target (i.e., blood) moving along a scan line A in a direction toward the ultrasonic transducer.
  • FIG. 1b illustrates a rearrangement of six groups of the return echoes, wherein the t axis represents the axial depth or spatial direction and the tau axis represents the temporal direction.
  • T represents the time delay, typically 200 microseconds, between the start of successive ultrasonic pulse transmissions. Note that the time shift between successive ones of the echoes is substantially uniform due to the relatively short time period between pulses as compared with the velocity of the moving target.
  • the autocorrelation method which is generally recognized to provide superior performance for real-time blood flow imaging than the other known techniques, will be briefly described in conjunction with FIG. 1b.
  • the autocorrelation type of time domain processing can be represented by the following equations:
  • a(t) is mostly determined by the impulse response of the ultrasound transducer. For multiple targets having the same velocity, Eq. (3) is valid. However, the phase response of a(t) will include interference from multiple targets. If the multiple targets have different velocities, Eq. (3) will not be valid. However, we can approximate this situation by considering a(t) as a broadband signal which has several frequency components. Consequently, our aim in flow imaging is to estimate mean and variance of the frequency spectrum of the demodulated signal given by Eq. (3).
  • each echo signal vector e n is multiplied in the axial direction by the complex conjugate of the adjacent echo signal vector e n-1 , which results in a plurality of pulse-pair vector signals e 1 *e 2 , e 2 *e 3 , e 3 *e 4 , . . . e n e n-1 *.
  • the amplitude of the resulting N-1 pulse-pair signals, at a certain axial depth represented by the Doppler time axis (the dashed line in FIG. 1b) are then averaged in the temporal (tau) direction.
  • the autocorrelation can be represented by two steps of pulse-pair vector calculation and averaging.
  • N is the number of temporal averagings.
  • the phase of the pulse-pair vector in Eq. (4) represents instantaneous frequency of the input Doppler signal which is changing from pulse to pulse.
  • the temporal averaging of Eq. (5) provides averaged vectors from which we can find a mean frequency.
  • Known autocorrelation velocity processing techniques suffer from the following two problems. Firstly, it is desirable that blood flow mapping have information updated on the order of 24-30 frames/second for best diagnostic effectiveness. Since existing autocorrelation techniques use temporal averaging, as exemplified by equations 1-7, averaging of a plurality (n) of received echo vectors is required. Thus, it is required to wait for the receipt of a plurality of echoes before an accurate analysis of the Doppler signal can be made. This results in a relatively low frame rate i.e., in the order of 15 frames/second. If less echoes were used the frame rate could be increased, however the accuracy of the velocity estimation would be severely degraded.
  • the Doppler signal being the signal amplitude variations along a Doppler axis, such as shown in FIG. 1c
  • the time shift between adjacent echoes is reduced.
  • this causes the period of the Doppler signal to become longer. Consequently, a greater number of echoes are required to be processed for accurate estimation of a slower velocity.
  • the number of received echoes which are averaged is decided by the maximum flow velocity (and is at least seven echo signals) so as to obtain accurate blood flow velocity mapping at about 15 frames/second.
  • the maximum flow velocity and is at least seven echo signals
  • even more echo signals are required to be processed to develope an accurate blood flow map which includes slower velocities, thereby reducing the frame rate even more.
  • a target i.e., blood
  • the sampling point is at a fixed distance from the transducer (i.e., depth in the body)
  • the result of conventional temporal averaging at that fixed depth is approximately the same as spatial averaging (i.e., averaging in the axial direction) of a single pulse-pair signal.
  • estimation by temporal averaging of a plurality of pulse-pair signals is equivalent to estimation by spatial averaging of a single pulse-pair signal.
  • a single pair of echo signals are spatially averaged in order to develope a velocity estimation.
  • the spatial averaging of the echo signals is accomplished in the vector domain.
  • each individual echo signal is digitized in accordance with the Nyquist rate for the pulse-pair signal. That is, each echo signal is digitized so as to provide signal samples at approximately twice the Nyquist rate of the individual echo signals.
  • the spatially averaged pulse-pair signal is also used to develop a turbulence estimation.
  • the spatial vector averaging technique is combined with the temporal averaging technique under noisy signal conditions or when it is desired to make performance trade-offs between signal-to-noise and resolution. This results in improved performance for weak signals as compared with the use of spatial averaging only, yet results in improved speed of calculation if only temporal averaging were used.
  • FIGS. 1a and 1b illustrate a series of pulse echos and a rearrangement of these echos, respectively, for illustrating the temporal averaging technique.
  • FIG. 2 illustrates two adjacent echo signals, which, when compared with FIG. 1(b) pictorially illustrates the approximate equivalence of the temporal and spatial averaging techniques.
  • FIG. 3 illustrates, in block diagram form, an ultrasound imaging apparatus including pulse-pair spatial vector averaging in accordance with the principles of the present invention.
  • FIG. 4 illustrates, in functional block diagram form, the spatial vector averaging technique used in the velocity and turbulence estimation processing portion of the apparatus of FIG. 3.
  • FIG. 5 illustrates, in block diagram form, a hardware implementation of the functional block diagram of FIG. 4.
  • FIG. 6 illustrates, in functional block diagram form, a modified version of the spatial vector averager of FIG. 4 which is particularly useful for weak signals.
  • FIG. 2 when viewed in conjunction with FIG. 1b, is useful for pictorially illustrating an equivalency between the temporal and spatial vector averaging techniques.
  • the amplitude of the detected baseband signal of the return echoes is periodically sampled along the Doppler axis, generating amplitude samples A 1 , A 2 . . . A i (however, only A 1 -A 6 are shown).
  • Processing in accordance with temporal averaging techniques for velocity estimation can be expressed in a digital form for equation (5) as follows: ##EQU4## where ⁇ is the sampling interval in the spatial direction.
  • FIG. 2 shows only two of the sequentially received echoes.
  • spatial sampling in the t direction results in paired amplitude samples A 1 A 2 , A 2 A 3 , A 3 A 4 . . . being generated which correspond with the same amplitude sample pairs A 1 A 2 , A 2 A 3 , A 3 A 4 . . . which were obtained by sampling along the Doppler axis in the temporal averaging technique.
  • Spatial vector averaging (averaging along the t axis) of the amplitude sample pairs A 1 A 2 , A 2 A 3 . . . of the single pulse-pair signal for velocity estimation can be expressed by ##EQU5## where M is the number of spatial averagings.
  • FIG. 3 illustrates a system block diagram of the inventive ultrasound medical imaging apparatus, and includes a reference oscillator 302 which provides a reference signal selected by the operator of the apparatus to be the operating frequency of an ultrasound transducer array 304.
  • the reference signal is applied to controllable delay transmission circuitry 306, which includes a signal generator for developing a series of bursts of the reference signals, each series of bursts lasting no more than 8 microseconds and being repeated at 200 microseconds intervals.
  • a plurality of controllable delay circuits is responsive to each series of bursts for causing transducer array 304 to produce a series of beams of ultrasound B 1 through B n which are steered over an angle alpha so as to provide a standard sector scan image of the interior of an object 308, as well known.
  • Transducer array 304 is also connected to a receiver circuit 310 which includes a plurality of controllable delay circuits, each delay circuit having a digitally selectable delay for properly combining the echo signals received from the individual transducers of transducer array 304, in a conventional and well known manner. Additionally, as also conventionally known, controllable gain amplifiers for providing a time-gain compression characteristic are also provided in receiver circuitry 310. During the approximately 190 microsecond intervals between ultrasound transmissions, transducer array 304 serves as an ultrasound receiver and converts the reflected ultrasound signals (echoes) into electrical signals. The electrical signals are combined in accordance with the controllable delay circuitry, as well known, so that the signals provided by each transducer of the phase transducer array 304 are processed simultaneously as reflections from a central point.
  • the output signal of receiver circuitry 310 is applied to an envelope detector 312 which produces a positive output signal proportional to the amplitude of the signal provided from receiver circuitry 10.
  • the amplitude signal provided from envelope detector 312 is then digitized by an analog-to-digital (A/D) converter 314.
  • A/D analog-to-digital
  • this digitized echo signal is applied to a digital scan converter for converting the sector scan echo into a raster scan format for subsequent display and/or recording of the ultrasound generated image signals.
  • the output signal of receiver circuitry 310 is also coupled to the input of signal multipliers 316 and 318.
  • Multipliers 316 and 318 comprise phase detectors having reference signal inputs S1 and S2, respectively.
  • S1 is the reference signal provided from reference oscillator 302 and S2 is of the same frequency as S1 but has a phase which is shifted 90° with respect to the phase of S1 due to phase shifter circuit 320.
  • multipliers 316 and 318 serve as in-phase and quadrature-phase synchronous detectors, and thereby provide in-phase (I) and quadrature-phase (Q) baseband vector components of the received echo signals.
  • A/D converters 322 and 324 generate digitized versions of the I and Q signals, respectively.
  • a spatial vector processor 326 produces digital blood flow estimation signals which, over time are representative of a color blood flow image, and which are applied to a digital scan converter 328.
  • Sptial vector processor 326 will be described in detail later on.
  • the digitized echo signal is applied to digital scan converter 330.
  • a sector controller 332 controls the reading of the digital signals into digital scan converters 328 and 330 in accordance with the sector format and a raster controller 334 controls read-out of the digital signals from digital scan converters 328 and 330 in accordance with the raster format.
  • a combiner 336 combines the digital signals from scan converters 328 and 330 so as to generate a signal representative of the anatomical image having superimposed thereover a color blood flow image representative of blood flow velocities.
  • Combiner 336 has an additional input connected to the output of a central processing unit (CPU) 338 which, among other things, provides textual information including, for example, the name of a patient, medical history and operating instructions which may be selected from an operating menu for the ultrasound apparatus, in the conventional manner.
  • CPU central processing unit
  • Combiner 336 provides red, green and blue output signals which are processed in a conventional manner by a video processor 340 and then applied to display and/or record devices 342.
  • the display devices may comprise a color and a black and white CRT's and the record devices may comprise a strip recorder and a video tape recorder.
  • FIG. 4 shows a functional block diagram of the spatial vector averager shown in FIG. 3.
  • the demodulated I and Q vector signal components from quadrature multipliers 316 and 318 of FIG. 3 are applied to low pass filters (LPF) 402 and 404 for removing noise signals which are received with the incoming echo signals, as well as unwanted high frequency signal components supplied from multipliers 316 and 318.
  • LPF low pass filters
  • A/D converters 406 and 408 digitize the I and Q components, respectively, in order that the remainder of the signal processing can be done in digital form.
  • each echo signal must provide signal samples to spatial processor 326 in accordance with the Nyquist rate of the pulse-pair signal which is processed therein. For example, if each echo signal has a bandwidth of 2.25 MHz, then echo signal samples must be provided to processor 326 at the 4.5 MHz rate.
  • A/D converters 406 and 408 would digitize the echo signals at this rate, however, it is presently more cost-effective to use slower speed, and hence lower cost A/D converter at 2.25 MHz and to use digital interpolators 410 and 412 (of conventional and well known construction) to arrive at the 4.5 MHz sample rate.
  • Third order digital high pass filters (HPF's) 414 and 416 filter the I and Q signals, respectively, for the purpose of wall motion rejection, i.e., for rejection of wall motions which are relatively slow as compared with the blood flow velocity.
  • HPF's 414 and 416 are also adjustable, for example, from 66 Hz to 526 Hz when using a 5 MHz sampling rate for the transmitted ultrasound signals in order to properly reject varying amounts of wall motion.
  • HPF's 414 and 416 are conventionally implemented as digital delay line cancellers (DLC) and are well known to those skilled in the art.
  • DLC digital delay line cancellers
  • spatial averaging velocity processing includes a covariance algorithm wherein:
  • Pulse-pair vectors are calculated by complex conjugate multiplication of successive echoes
  • the mean frequency (velocity) of the blood flow is determined by calculating the phase of the averaged vectors.
  • a complex conjugate multiplier 418 receives the I and Q echo signal components and calculates pulse-pair vector components using successive echo signals by complex conjugate multiplications in accordance with equation (4).
  • the real and imaginary parts of the pulse-pair vectors are spatially averaged over their coherence interval (pulse width) by spatial averagers 420 and 422, respectively.
  • the mean frequency (velocity) is determined by computing the phase of the echo vector using a tan -1 circuit 424. That is, by calculating the tan -1 of the ratio of imaginary and real portions of the pulse-pair vector, a velocity estimation is produced.
  • a turbulence estimation is also produced. As shown in equation 7, the turbulence estimation is proportional to the ratio of the absolute value of R(T) and R(0).
  • R(0) corresponds to the power of the received echo signal, and is calculated by summing the squares of the I and Q echo signal components. This is accomplished by squaring circuit 426.
  • the calculated power is spatially averaged by averager 428 to the same degree of spatial averaging provided by averagers 420 and 422, and applied to one input of a turbulence estimation calculator 430.
  • An absolute value circuit 432 provides the absolute value of R(T) to another input of turbulence estimation calculator 430.
  • the turbulence estimation is provided at the output of calculator 430.
  • FIG. 5 illustrates in block diagram form a hardware implementation for the FIG. 4 arrangement. Complex conjugate multiplication of the echo vector signals is accomplished as shown below:
  • the output of mixer 506 is I n I n-1
  • mixer 508 is I n-1 Q n
  • mixer 510 is Q n Q n-1
  • mixer 512 is I n Q n-1
  • the output of adder 514 is R i
  • the output of adder 516 is R j .
  • Adders 518 and 520 provide recursive additions of the R i and R j vector components, respectively, each addition providing a spatial averaging of the respective vector components due to time delay circuits 522 and 524 coupled in a feedback manner to adders 518 and 520, respectively.
  • the delay ( ⁇ ) provided by time delay circuits 522 and 524 is significantly shorter than the delay provided by circuits 502 and 504 so as to effectuate spatial averaging, i.e., the delays correspond with the time periods between successive samples A 1 , A 2 , A 3 . . . A i of each echo signal of FIG. 2.
  • the number of spatial averagings is set by periodically resetting the output of averagers 420 and 422 to zero after e.g., 5-20 successive spatial averagings of the pulse-pair signal samples. It is noted that although averaging beyond the coherence interval increases the signal-to-noise of the estimation, it reduces its resolution.
  • the tan -1 calculation is provided by a look-up table (memory) 525 which receives the spatially averaged digital R i and R j signals as an address and is programmed to provide as its output the tan -1 of R j /R i , which is representative of the mean velocity of the blood flow in a selected sample volume of the ultrasound sector scan.
  • a look-up table (memory) 525 which receives the spatially averaged digital R i and R j signals as an address and is programmed to provide as its output the tan -1 of R j /R i , which is representative of the mean velocity of the blood flow in a selected sample volume of the ultrasound sector scan.
  • circuits 426 and 432 comprise squaring circuits 526, 528 and adder 530, and squaring circuits 532, 534 and adder 536, respectively.
  • the turbulence calculation is provided by a look-up table 538 which receives the processed digital R(0) and R(T) signals as an address and is programmed to provide as its output a signal corresponding to ##EQU6## which is proportional to the turbulence estimation for the blood flow in the selected sample volume of the ultrasound sector scan.
  • a square root circuit should be provided after adder 536, however, in practice the effect of the square root circuit can be pre-programmed into look-up table 538, thereby eliminating the need for a square root circuit.
  • the previously described spatial averaging technique can be modified so as to include temporal signal processing in conjunction with spatial averaging, when desired.
  • This type of processing is referred to herein as two-dimensional vector averaging because the echo vector signals are processed and averaged in a matrix of two dimensions, one axis of the matrix is the axial (spatial) direction and an orthogonal axis is the t (temporal) direction.
  • spatial averagers 420 and 422 of FIG. 4 would each be replaced by a spatial/temporal averager 602, such as shown in FIG. 6.
  • Averager 602 comprises a spatial averager such as illustrated in FIGS. 4 and 5, followed by a temporal averager of conventional design.
  • the temporal averager would temporally process corresponding signal samples from pairs of successive echo signals, for example 4 pairs (using 5 sequential echo signals), so as to generate in combination with, for example, 5 spatial averagings for each pulse-pair, 20 averagings.
  • performance trade-offs are also possible, e.g., for a higher frame-rate mode, 20 spatial averagings only could be used, i.e., no temporal averagings.

Abstract

Method and apparatus for developing an estimation of the velocity of a moving target, comprising, transmitting groups of pulse signals towards the target, receiving groups of echo signals caused by reflections of the transmitted groups of pulse signals from the target and spatially processing the echo signals for developing the estimation of the velocity of the target. In a preferred embodiment, the spatial processing is for only a single pair of adjacent ones of ultrasonic echo signals and the estimation developes a two dimensional image of blood flow velocities in the body of a patient.

Description

BACKGROUND OF THE INVENTION
1. Field of the Invention
This invention relates to velocity estimation using Doppler techniques, and more particularly to a spatial vector averaging method and apparatus for estimating blood flow velocities in a medical diagnostic apparatus using ultrasonic pulse signals.
2. Background of the Invention
Medical diagnostic ultrasound apparatus generate images of anatomical structures within a body of a patient by transmitting ultrahigh-frequency sound waves (typically on the order of 3.0 MHz) into the patient and then analyzing the echoes, i.e., ultrasonic signals reflected from the structure of the body being scanned. Perhaps the most widely used ultrasound diagnostic apparatus displays in real-time anatomical information in the form of a two-dimensional image of a selected cross-section of the structure. The ultrasound signals are swept across the structure in the form of a sector scan. The sector scan is performed in real-time so that the image is available during the examination of the patient. In such cases, motion of the structure produces a corresponding moving image (i.e., a B-mode image).
In some clinical applications, e.g., cardiac imaging, anatomical defects can be relatively small and beyond the resolution capabilities of conventional anatomical ultrasound imaging. However, since a small anatomical defect may manifest itself by a pronounced change in blood flow velocity, for example aortic stenosis, mitral or aortic insufficiency or a congenital defect, a display of blood flow velocity would allow these abnormalities to be more easily detected. One known method for velocity display is to use FFT techniques to process the echo signals reflected from a selected small volume so as to generate a numerical display. This method is severely limited by the fact that velocity is only determined for the small sample volume and is also not a two dimensional real-time image. Real-time imaging of velocity in a larger area is highly desireable. Thus, velocity blood flow imaging has become an increasingly important portion of the ultrasound imaging device used in the medical diagnostic field, wherein a real-time blood flow image is superimposed over a real-time anatomical image. However, it is difficult to acquire sufficient ultrasound data to develop an accurate and high resolution blood flow image which can be displayed in real-time at a sufficiently high rate. This is because many echos need to be processed in a short time while the physical reality is that ultrasound signals have a relatively slow propagation speed in a human body, thereby limiting the number of echoes which can be received in a short time.
European Patent Publication No. 0 100 094 by Namekawa of the Aloka Company describes an ultrasound blood flow imaging apparatus for two-dimensional display i.e., mapping, of blood flow velocities within a body. The flow mapping is superimposed over a B-mode scan and displays blood flow in a color which is representative of the direction of the blood flow with respect to the ultrasonic transducer while variations in the intensity of the colors are representative of the blood flow velocities. It is only noted therein that the blood flow velocity signal is generated using a pulsed-Doppler method.
Present Dopper velocity estimators use time-domain processing techniques wherein a series of N repetitive pulses of RF-ultrasonic signal separated by time periods of T are transmitted towards a moving target along a given scan direction. A relatively large number of echo signals (from a minimum of 10 to a maximum of 256) are processed using Doppler techniques for determining the velocity and turbulence of the moving object. One such technique for determining Doppler frequency shift is shown, for example, by U.S. Pat. No. 4,542,657 issued to Barber et al. which uses I and Q sampled signals of a plurality (i.e., from 16 to 256) of the demodulated echo signals. FFT and zero-crossing velocity estimators are described in U.S. Pat. No. 4,318,413 for processing the Doppler signal. Another time-domain processing technique uses autocorrelaton (pulse-pair) algorithms, such as known, for example, from an article by Kasai et al. published in the IEEE Transactions on Sonics and Ultrasonics, Volume SU-32, No. 3, May 1985. Each of the above time-domain processing techniques evaluate a relatively large number of return echoes (from 10 to 256) for estimating blood flow velocity and turbulence.
For a better understanding of the time-domain processing technique and its relation to the present invention, reference is now made to FIG. 1a herein, which shows groups (or pulses) of ultrasonic signal 1, 2, 3, . . . n which are excited and groups of return echo signals e1, e2, e3 . . . en which are received in response thereto. The echoes are reflections from a target (i.e., blood) moving along a scan line A in a direction toward the ultrasonic transducer. FIG. 1b illustrates a rearrangement of six groups of the return echoes, wherein the t axis represents the axial depth or spatial direction and the tau axis represents the temporal direction. T represents the time delay, typically 200 microseconds, between the start of successive ultrasonic pulse transmissions. Note that the time shift between successive ones of the echoes is substantially uniform due to the relatively short time period between pulses as compared with the velocity of the moving target.
The autocorrelation method, which is generally recognized to provide superior performance for real-time blood flow imaging than the other known techniques, will be briefly described in conjunction with FIG. 1b. The autocorrelation type of time domain processing can be represented by the following equations:
For simplicity of formulations, we consider reflected echoes from a single target to be represented by
Z(t) =a(t) cos [w.sub.0 t+φ(t)]                        (1)
where a(t) is the ultrasonic signal pulse envelope, w0 is the carrier frequency, and φ(t) is the phase response. When the target moves by a time shift αT during a pulse repetition period T, Eq. 1 for the n-th echo becomes
Z.sub.n (t)=a(t-αnT) cos [w.sub.0 (t-αnT)+φ(t-αnT)](2)
where α is the Doppler ratio given by α=wd /w0 and wd is the Doppler frequency.
After quadrature demodulation, the demodulated signal en (t) of Eq. (2) can be written as ##EQU1##
In Eq. (3), a(t) is mostly determined by the impulse response of the ultrasound transducer. For multiple targets having the same velocity, Eq. (3) is valid. However, the phase response of a(t) will include interference from multiple targets. If the multiple targets have different velocities, Eq. (3) will not be valid. However, we can approximate this situation by considering a(t) as a broadband signal which has several frequency components. Consequently, our aim in flow imaging is to estimate mean and variance of the frequency spectrum of the demodulated signal given by Eq. (3).
In the known autocorrelation processing method, each echo signal vector en is multiplied in the axial direction by the complex conjugate of the adjacent echo signal vector en-1, which results in a plurality of pulse-pair vector signals e1 *e2, e2 *e3, e3 *e4, . . . en en-1 *. The amplitude of the resulting N-1 pulse-pair signals, at a certain axial depth represented by the Doppler time axis (the dashed line in FIG. 1b) are then averaged in the temporal (tau) direction.
The autocorrelation can be represented by two steps of pulse-pair vector calculation and averaging. ##EQU2## where N is the number of temporal averagings. The phase of the pulse-pair vector in Eq. (4) represents instantaneous frequency of the input Doppler signal which is changing from pulse to pulse. The temporal averaging of Eq. (5) provides averaged vectors from which we can find a mean frequency.
From the amplitude and phase of the resultant averaged vectors, the mean frequency, which corresponds with the velocity of the blood, can be obtained. Furthermore, the variance (sigma squared) of the velocity can be obtained, which corresponds to the turbulence of the blood flow. Turbulence estimations also provide useful diagnostic information when displayed. The phase (velocity) and variance are calculated as follows: ##EQU3## wherein R(T) is an abbreviated notation of R(T; nT, t) and R(T)=Rr (T)+jRi (T)
For further details concerning this pulse-pair autocorrelation technique, the reader is referred to the forenoted article by Kasai et al.
Known autocorrelation velocity processing techniques suffer from the following two problems. Firstly, it is desirable that blood flow mapping have information updated on the order of 24-30 frames/second for best diagnostic effectiveness. Since existing autocorrelation techniques use temporal averaging, as exemplified by equations 1-7, averaging of a plurality (n) of received echo vectors is required. Thus, it is required to wait for the receipt of a plurality of echoes before an accurate analysis of the Doppler signal can be made. This results in a relatively low frame rate i.e., in the order of 15 frames/second. If less echoes were used the frame rate could be increased, however the accuracy of the velocity estimation would be severely degraded. Furthermore, since turbulence is in effect a differentiation of velocity with respect to time, due to amplitude fluctuations of the velocity signal, the receipt of even more echoes is required in order to arrive at an accurate turbulence estimation. The second problem results from the fact that the period of the amplitude and phase variations of the Doppler signal (the Doppler signal being the signal amplitude variations along a Doppler axis, such as shown in FIG. 1c) depends upon the blood flow velocity. In slower flows, the time shift between adjacent echoes is reduced. As can be easily visualized in FIG. 1a and 1b, this causes the period of the Doppler signal to become longer. Consequently, a greater number of echoes are required to be processed for accurate estimation of a slower velocity. In the known ultrasound systems employing velocity estimation, the number of received echoes which are averaged is decided by the maximum flow velocity (and is at least seven echo signals) so as to obtain accurate blood flow velocity mapping at about 15 frames/second. However, even more echo signals are required to be processed to develope an accurate blood flow map which includes slower velocities, thereby reducing the frame rate even more.
It is an object of the present invention to provide a Doppler blood flow velocity and turbulence estimator which minimizes the number of received echoes required for accurate estimation of blood flow velocity and turbulence in order that an accurate color blood flow mapping can be provided at a relatively high frame rate, for example 24 or 30 frames/second.
SUMMARY OF THE INVENTION
In accordance with the principles of the present invention, advantage is taken by the inventor of recognition of the fact that under some conditions, for example when a target (i.e., blood) is moving and the sampling point is at a fixed distance from the transducer (i.e., depth in the body), the result of conventional temporal averaging at that fixed depth is approximately the same as spatial averaging (i.e., averaging in the axial direction) of a single pulse-pair signal. Thus, in this case, estimation by temporal averaging of a plurality of pulse-pair signals is equivalent to estimation by spatial averaging of a single pulse-pair signal. Consequently, in accordance with the present invention, a single pair of echo signals are spatially averaged in order to develope a velocity estimation. In the preferred embodiment the spatial averaging of the echo signals is accomplished in the vector domain.
In accordance with a further aspect of the invention, since the processed pulse-pair signal has approximately twice the bandwidth of an individual echo signal, in order to ensure proper characterization of the echo signals for accurate digital processing, each individual echo signal is digitized in accordance with the Nyquist rate for the pulse-pair signal. That is, each echo signal is digitized so as to provide signal samples at approximately twice the Nyquist rate of the individual echo signals.
In accordance with a still further aspect of the invention the spatially averaged pulse-pair signal is also used to develop a turbulence estimation.
In accordance with an even further aspect of the invention the spatial vector averaging technique is combined with the temporal averaging technique under noisy signal conditions or when it is desired to make performance trade-offs between signal-to-noise and resolution. This results in improved performance for weak signals as compared with the use of spatial averaging only, yet results in improved speed of calculation if only temporal averaging were used.
Other features and advantages of the invention will be apparent from the following description of the preferred embodiments and from the claims.
For a fuller understanding of the present invention, reference should now be made to the following detailed description of the preferred embodiments of the invention and the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
FIGS. 1a and 1b illustrate a series of pulse echos and a rearrangement of these echos, respectively, for illustrating the temporal averaging technique.
FIG. 2 illustrates two adjacent echo signals, which, when compared with FIG. 1(b) pictorially illustrates the approximate equivalence of the temporal and spatial averaging techniques.
FIG. 3 illustrates, in block diagram form, an ultrasound imaging apparatus including pulse-pair spatial vector averaging in accordance with the principles of the present invention.
FIG. 4 illustrates, in functional block diagram form, the spatial vector averaging technique used in the velocity and turbulence estimation processing portion of the apparatus of FIG. 3.
FIG. 5 illustrates, in block diagram form, a hardware implementation of the functional block diagram of FIG. 4.
FIG. 6 illustrates, in functional block diagram form, a modified version of the spatial vector averager of FIG. 4 which is particularly useful for weak signals.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
FIG. 2, when viewed in conjunction with FIG. 1b, is useful for pictorially illustrating an equivalency between the temporal and spatial vector averaging techniques. As shown in FIG. 1b, the amplitude of the detected baseband signal of the return echoes is periodically sampled along the Doppler axis, generating amplitude samples A1, A2 . . . Ai (however, only A1 -A6 are shown). Processing in accordance with temporal averaging techniques for velocity estimation, can be expressed in a digital form for equation (5) as follows: ##EQU4## where Δ is the sampling interval in the spatial direction.
FIG. 2 shows only two of the sequentially received echoes. As visually illustrated therein, spatial sampling (in the t direction) results in paired amplitude samples A1 A2, A2 A3, A3 A4 . . . being generated which correspond with the same amplitude sample pairs A1 A2, A2 A3, A3 A4 . . . which were obtained by sampling along the Doppler axis in the temporal averaging technique. Spatial vector averaging (averaging along the t axis) of the amplitude sample pairs A1 A2, A2 A3 . . . of the single pulse-pair signal for velocity estimation can be expressed by ##EQU5## where M is the number of spatial averagings.
Note that in equation 8, the summation variable k modifies T, the time period between successive ones of the echo signals, while in equation 9, the summation variable 1 modifies Δ, the time period between adjacent spatial samples. By comparing equations 8 and 9 and FIGS. 1b and 2 it is intuitively seen that with these given circumstances, temporal vector averaging of paired samples of a plurality of pulse-pair signals for obtaining a velocity estimation is approximately equivalent to spatial vector averaging of paired samples of a single pulse-pair signal for obtaining a velocity estimation. The equivalence is only approximate since in practice many times more paired samples can be spatially averaged as can be temporally averaged in a given time period (since only two echoes have to be received for this spatial averaging technique), thus resulting in a more accurate velocity estimation when using spatial vector averaging as compared with temporal averaging.
FIG. 3 illustrates a system block diagram of the inventive ultrasound medical imaging apparatus, and includes a reference oscillator 302 which provides a reference signal selected by the operator of the apparatus to be the operating frequency of an ultrasound transducer array 304. The reference signal is applied to controllable delay transmission circuitry 306, which includes a signal generator for developing a series of bursts of the reference signals, each series of bursts lasting no more than 8 microseconds and being repeated at 200 microseconds intervals. A plurality of controllable delay circuits, each having a digitally selectable delay, is responsive to each series of bursts for causing transducer array 304 to produce a series of beams of ultrasound B1 through Bn which are steered over an angle alpha so as to provide a standard sector scan image of the interior of an object 308, as well known.
Transducer array 304 is also connected to a receiver circuit 310 which includes a plurality of controllable delay circuits, each delay circuit having a digitally selectable delay for properly combining the echo signals received from the individual transducers of transducer array 304, in a conventional and well known manner. Additionally, as also conventionally known, controllable gain amplifiers for providing a time-gain compression characteristic are also provided in receiver circuitry 310. During the approximately 190 microsecond intervals between ultrasound transmissions, transducer array 304 serves as an ultrasound receiver and converts the reflected ultrasound signals (echoes) into electrical signals. The electrical signals are combined in accordance with the controllable delay circuitry, as well known, so that the signals provided by each transducer of the phase transducer array 304 are processed simultaneously as reflections from a central point.
The output signal of receiver circuitry 310 is applied to an envelope detector 312 which produces a positive output signal proportional to the amplitude of the signal provided from receiver circuitry 10. The amplitude signal provided from envelope detector 312 is then digitized by an analog-to-digital (A/D) converter 314. As will be subsequently described, this digitized echo signal is applied to a digital scan converter for converting the sector scan echo into a raster scan format for subsequent display and/or recording of the ultrasound generated image signals.
The output signal of receiver circuitry 310 is also coupled to the input of signal multipliers 316 and 318. Multipliers 316 and 318 comprise phase detectors having reference signal inputs S1 and S2, respectively. S1 is the reference signal provided from reference oscillator 302 and S2 is of the same frequency as S1 but has a phase which is shifted 90° with respect to the phase of S1 due to phase shifter circuit 320. As well known, with this arrangement multipliers 316 and 318 serve as in-phase and quadrature-phase synchronous detectors, and thereby provide in-phase (I) and quadrature-phase (Q) baseband vector components of the received echo signals. A/D converters 322 and 324 generate digitized versions of the I and Q signals, respectively.
A spatial vector processor 326 produces digital blood flow estimation signals which, over time are representative of a color blood flow image, and which are applied to a digital scan converter 328. Sptial vector processor 326 will be described in detail later on. As previously noted, the digitized echo signal is applied to digital scan converter 330. As is conventional in the art, a sector controller 332 controls the reading of the digital signals into digital scan converters 328 and 330 in accordance with the sector format and a raster controller 334 controls read-out of the digital signals from digital scan converters 328 and 330 in accordance with the raster format. A combiner 336 combines the digital signals from scan converters 328 and 330 so as to generate a signal representative of the anatomical image having superimposed thereover a color blood flow image representative of blood flow velocities. Combiner 336 has an additional input connected to the output of a central processing unit (CPU) 338 which, among other things, provides textual information including, for example, the name of a patient, medical history and operating instructions which may be selected from an operating menu for the ultrasound apparatus, in the conventional manner.
Combiner 336 provides red, green and blue output signals which are processed in a conventional manner by a video processor 340 and then applied to display and/or record devices 342. The display devices may comprise a color and a black and white CRT's and the record devices may comprise a strip recorder and a video tape recorder.
For further details concerning the structure and operation of FIG. 3, except for spatial processor 326, the reader is referred to U.S. Pat. No. 4,612,937 of Lawrence R. Miller, assigned to the same assignee as the present invention and which is incorporated herein by reference.
Reference is now made to FIG. 4 which shows a functional block diagram of the spatial vector averager shown in FIG. 3. The demodulated I and Q vector signal components from quadrature multipliers 316 and 318 of FIG. 3 are applied to low pass filters (LPF) 402 and 404 for removing noise signals which are received with the incoming echo signals, as well as unwanted high frequency signal components supplied from multipliers 316 and 318. Next, A/ D converters 406 and 408 digitize the I and Q components, respectively, in order that the remainder of the signal processing can be done in digital form. As previously noted, since spatial vector averaging uses significantly more information (bandwidth) of each received echo signal as compared with temporal processing, it is important that the digital characterization of each echo signal have enough samples so that the subsequent digital processing is accurate. Thus, each echo signal must provide signal samples to spatial processor 326 in accordance with the Nyquist rate of the pulse-pair signal which is processed therein. For example, if each echo signal has a bandwidth of 2.25 MHz, then echo signal samples must be provided to processor 326 at the 4.5 MHz rate. Ideally, A/ D converters 406 and 408 would digitize the echo signals at this rate, however, it is presently more cost-effective to use slower speed, and hence lower cost A/D converter at 2.25 MHz and to use digital interpolators 410 and 412 (of conventional and well known construction) to arrive at the 4.5 MHz sample rate. Third order digital high pass filters (HPF's) 414 and 416 filter the I and Q signals, respectively, for the purpose of wall motion rejection, i.e., for rejection of wall motions which are relatively slow as compared with the blood flow velocity. The cut-off frequency of HPF's 414 and 416 are also adjustable, for example, from 66 Hz to 526 Hz when using a 5 MHz sampling rate for the transmitted ultrasound signals in order to properly reject varying amounts of wall motion. HPF's 414 and 416 are conventionally implemented as digital delay line cancellers (DLC) and are well known to those skilled in the art.
As in the case with the temporal type of velocity processing, spatial averaging velocity processing includes a covariance algorithm wherein:
(1) Pulse-pair vectors are calculated by complex conjugate multiplication of successive echoes,
(2) The real and imaginary parts of the resultant vectors are averaged, and
(3) The mean frequency (velocity) of the blood flow is determined by calculating the phase of the averaged vectors.
Thus, a complex conjugate multiplier 418 receives the I and Q echo signal components and calculates pulse-pair vector components using successive echo signals by complex conjugate multiplications in accordance with equation (4).
Next, the real and imaginary parts of the pulse-pair vectors are spatially averaged over their coherence interval (pulse width) by spatial averagers 420 and 422, respectively. The mean frequency (velocity) is determined by computing the phase of the echo vector using a tan-1 circuit 424. That is, by calculating the tan-1 of the ratio of imaginary and real portions of the pulse-pair vector, a velocity estimation is produced.
A turbulence estimation is also produced. As shown in equation 7, the turbulence estimation is proportional to the ratio of the absolute value of R(T) and R(0). R(0) corresponds to the power of the received echo signal, and is calculated by summing the squares of the I and Q echo signal components. This is accomplished by squaring circuit 426. For accuracy, the calculated power is spatially averaged by averager 428 to the same degree of spatial averaging provided by averagers 420 and 422, and applied to one input of a turbulence estimation calculator 430. An absolute value circuit 432 provides the absolute value of R(T) to another input of turbulence estimation calculator 430. The turbulence estimation is provided at the output of calculator 430.
FIG. 5 illustrates in block diagram form a hardware implementation for the FIG. 4 arrangement. Complex conjugate multiplication of the echo vector signals is accomplished as shown below:
e.sub.n =I.sub.n +jQ.sub.n
e.sub.n-1 =I.sub.n-1 +jQ.sub.n-1
e.sub.n e.sub.n-1 *=R.sub.i +jR.sub.j
wherein
R.sub.i =I.sub.n I.sub.n-1 +Q.sub.n Q.sub.n-1
and
R.sub.j =I.sub.n Q.sub.n-1 -I.sub.n-1 Q.sub.n
Thus, with the delay provided by delay circuits 502 and 504 equal to the time delay T between successive echo signals, the output of mixer 506 is In In-1, mixer 508 is In-1 Qn, mixer 510 is Qn Qn-1 and mixer 512 is In Qn-1. Thus, the output of adder 514 is Ri and the output of adder 516 is Rj. Adders 518 and 520 provide recursive additions of the Ri and Rj vector components, respectively, each addition providing a spatial averaging of the respective vector components due to time delay circuits 522 and 524 coupled in a feedback manner to adders 518 and 520, respectively. The delay (Δ) provided by time delay circuits 522 and 524 is significantly shorter than the delay provided by circuits 502 and 504 so as to effectuate spatial averaging, i.e., the delays correspond with the time periods between successive samples A1, A2, A3 . . . Ai of each echo signal of FIG. 2. The number of spatial averagings is set by periodically resetting the output of averagers 420 and 422 to zero after e.g., 5-20 successive spatial averagings of the pulse-pair signal samples. It is noted that although averaging beyond the coherence interval increases the signal-to-noise of the estimation, it reduces its resolution. The tan-1 calculation is provided by a look-up table (memory) 525 which receives the spatially averaged digital Ri and Rj signals as an address and is programmed to provide as its output the tan-1 of Rj /Ri, which is representative of the mean velocity of the blood flow in a selected sample volume of the ultrasound sector scan.
For turbulence estimation, circuits 426 and 432 comprise squaring circuits 526, 528 and adder 530, and squaring circuits 532, 534 and adder 536, respectively. The turbulence calculation is provided by a look-up table 538 which receives the processed digital R(0) and R(T) signals as an address and is programmed to provide as its output a signal corresponding to ##EQU6## which is proportional to the turbulence estimation for the blood flow in the selected sample volume of the ultrasound sector scan. Note that in theory, a square root circuit should be provided after adder 536, however, in practice the effect of the square root circuit can be pre-programmed into look-up table 538, thereby eliminating the need for a square root circuit.
It is noted that good turbulence estimation generally requires averaging of ten or more echoes. Thus, the above-described spatial vector averaging which uses only a single pulse-pair (i.e., two echoes) will provide a lower resolution turbulence estimation than conventional temporal averaging. Furthermore, it is known that for weak Doppler signals, more averagings are generally required for accurate estimations. Therefore, in accordance with a further aspect of the present invention, the previously described spatial averaging technique can be modified so as to include temporal signal processing in conjunction with spatial averaging, when desired. This type of processing is referred to herein as two-dimensional vector averaging because the echo vector signals are processed and averaged in a matrix of two dimensions, one axis of the matrix is the axial (spatial) direction and an orthogonal axis is the t (temporal) direction. In this case, spatial averagers 420 and 422 of FIG. 4 would each be replaced by a spatial/temporal averager 602, such as shown in FIG. 6. Averager 602 comprises a spatial averager such as illustrated in FIGS. 4 and 5, followed by a temporal averager of conventional design. The temporal averager would temporally process corresponding signal samples from pairs of successive echo signals, for example 4 pairs (using 5 sequential echo signals), so as to generate in combination with, for example, 5 spatial averagings for each pulse-pair, 20 averagings. This results in improved resolution turbulence estimations as compared with the use of spatial averagings only, and improved resolution velocity estimations as compared with the use of spatial averaging only during the reception of weak signals. Of course, performance trade-offs are also possible, e.g., for a higher frame-rate mode, 20 spatial averagings only could be used, i.e., no temporal averagings.
Thus, there has been shown and described novel apparatus for processing ultrasound echo signals for developing blood flow estimation signals which fulfull all the objects and advantages sought therefore. Many changes, modifications, variations and other uses and applications of the subject invention will, however, become apparent to those skilled in the art after considering this specification and the accompanying drawings which disclose only preferred embodiments thereof. For example, the described spatial processing estimation techniques could be applied in fields other than ultrasound, for example, in radar or other types of Doppler systems. Additionally, portions of the signal processing shown as being implemented in hardware in FIG. 5, could be implemented in software. All such changes, modifications, variations and other uses and applications which do not depart from the spirit and scope of the invention are deemed to be covered by the invention as defined only by the claims which follow.

Claims (37)

I claim:
1. A method for developing an estimation of the velocity of a moving target, comprising:
transmitting groups of pulse signals towards said target;
receiving groups of echo signals caused by reflections of said transmitted groups of pulse signals from said target;
converting said received echo signals into real and imaginary vector signal components;
vector processing of two temporally successive groups of said echo signals in their vector component format so as to develope real and imaginary processed vector signal components;
spatial averaging of said real and imaginary processed vector components; and calculating said velocity estimation from said spatially averaged real and imaginary processed vector signal components.
2. The method of claim 1, wherein said vector processing step comprises:
vector processing of only a single pair of temporally successive groups of said echo signals for use in developing said estimation of the velocity of said target.
3. The method of claim 1, wherein said vector processing and spatial averaging steps comprise:
correlation processing of said two temporally successive groups of echo signals.
4. The method of claim 3, wherein said correlation processing step comprises:
complex conjugate multiplication of said two temporally successive groups of echo signals.
5. The method of claim 3, wherein said calculating step comprises:
calculating the arctan of a ratio of the spatially averaged imaginary and real processed vector components for developing a signal representative of said velocity estimation.
6. The method of claim 1, wherein said calculating step comprises:
calculating the arctan of a ratio of the spatially averaged imaginary and real processed vector components for developing a signal representative of said velocity estimation.
7. The method of claim 1, further comprising:
digitizing said groups of echo signals before said vector processing step, said digitizing resulting in echo signals samples being generated at a rate of approximately twice the Nyquist rate for each of said groups of echo signals.
8. The method of claim 1 wherein said spatial averaging step comprises:
spatial averaging of a plurality of pairs of temporally successive ones of said echo signals in combination with temporal averaging of the results of said spatial averaging for each of said plurality of pairs of temporally successive echo signals for developing said velocity estimation.
9. The method of claim 8, comprising the further step of:
calculating a turbulence estimation using said spatially and temporally averaged echo signals.
10. The method of claim 1, including the further steps of:
repeatedly calculating velocity estimations in accordance with said prior steps so as to develope a two dimensional array of said velocity estimations over a given scan area; and
displaying said two dimensional array of velocity estimations so as to develope an image of velocity estimations within said given scan area.
11. A method according to claim 1, comprising the further step of:
calculating a turbulence estimation using said spatially averaged processed vector components.
12. A method for developing an estimation of the velocity of a moving target, comprising:
transmitting groups of signals toward said target;
receiving groups of echo signals caused by reflections of said transmitted groups of pulse signals from said target; and
digital autocorrelation processing with spatial averaging of said echo signals for developing said estimation of the velocity of said moving target, wherein a digitizing step precedes said autocorrelation processing and spatial averaging step and generates digitized samples of said echo signals at a minimum rate of approximately twice the Nyquist rate for an individual one of said echo signals.
13. The method according to claim 12 wherein said autocorrelation processing with spatial averaging step comprises:
complex conjugate multiplication of two adjacent groups of said echo signals for developing vector component signals and spatial averaging of said vector component signals.
14. The method according to claim 13 further including:
calculating the arctan of the ratio of the spatially averaged vector components for developing a signal representative of said velocity estimation.
15. The method according to claim 14 further including:
calculating a turbulence estimation using said spatially averaged vector components.
16. The method of claim 12, including the further steps of:
repeatedly calculating velocity estimations in accordance with said prior steps so as to develope a two dimensional array of said velocity estimations over a given scan area; and
displaying said two dimensional array of velocity estimations so as to develope an image of velocity estimations within said given scan area.
17. A method for developing an estimation of the velocity of a moving target, comprising:
transmitting groups of pulse signals towards said target;
receiving groups of echo signals caused by reflections of said transmitted groups of pulse signals from said target;
complex conjugate multiplication of a pair of temporally successive groups of said echo signals for developing vector component signals;
spatial averaging of said vector component signals; and
repetition of said above two steps for a successive pair of temporally successive groups of received echo signals for developing said estimation of the velocity of said target.
18. The method of claim 17, including the further steps of:
repeatedly calculating velocity estimations in accordance with said prior steps so as to develope a two dimensional array of said velocity estimations over a given scan area; and
displaying said two dimensional array of velocity estimations so as to develope an image of velocity estimations within said given scan area.
19. Apparatus for developing an estimation of the velocity of a target, comprising:
means for transmitting groups of pulse signals towards said target;
means for receiving groups of echo signals caused by reflections of said transmitted groups of pulse signals from said target;
spatial averaging means for spatially averaging said echo signals for developing said estimation of the velocity of said target, said spatial averaging means including delay means having a time delay equal to the time delay between adjacent groups of said transmitted pulse signals; and
an adder having said delay means coupled thereto in a feedback manner so as to repeatedly perform spatial averagings of at least one pair of temporally successive ones of said echo signals.
20. Apparatus according to claim 19 further including:
means for calculating a turbulence estimation from said spatially averaged echo signals.
21. The apparatus of claim 19, wherein:
said spatial averaging means averages only a single pair of temporally successive ones of said echo signals for developing said estimation of the velocity of said target.
22. The apparatus of claim 19, wherein said spatial averaging means includes:
means for converting said received echo signals into real and imaginary vector signal components;
means for vector processing of two temporally successive groups of said echo signals in their vector component format so as to develop real and imaginary processed vector signal components;
means for spatial averaging of said real and imaginary processed vector components; and
means for calculating said velocity estimation from said spatially averaged real and imaginary processed vector signal components.
23. The apparatus of claim 22, wherein said vector processing means comprises:
means for correlation processing of said two temporally successive groups of said echo signals.
24. The apparatus of claim 23, wherein said correlation processing means comprises:
a complex conjugate multiplier for multiplying said two temporally successive groups of echo signals together in a complex conjugate manner.
25. The apparatus of claim 24, wherein said calculating means comprises:
means for calculating the arctan of a ratio of the spatially averaged imaginary and real processed vector components for developing a signal representative of said velocity estimation.
26. The apparatus of claim 25, wherein said means for calculating the arctan comprises:
a look-up table memory which is addressed by said spatially averaged imaginary and real vector components.
27. The apparatus of claim 19, further including:
means for digitizing said groups of echo signals so as to supply to said spatial averaging means echo signal samples at a rate of approximately twice the Nyquist rate for each of said groups of echo signals.
28. Apparatus according to claim 19, further including:
means for temporally averaging the results of said spatial averaging means, for developing said estimation of the velocity of said target.
29. Apparatus for developing an estimation of the velocity of a moving target, comprising:
means for transmitting groups of signals toward said target;
means for receiving groups of echo signals caused by reflections of said transmitted groups of pulse signals from said target;
means for digitizing said echo signals; and
means for digital autocorrelation processing with spatial averaging of said digitized echo signals for developing said estimation of the velocity of said moving target, wherein said digitizing means generates digitized samples of said echo signals at a minimum rate of approximately twice the Nyquist rate for an individual one of said echo signals.
30. Apparatus according to claim 29 wherein said means for autocorrelation processing with spatial averaging comprises:
a complex conjugate multiplier for multiplying said two adjacent groups of said echo signals together in a complex conjugate manner for developing vector component signals and spatial averaging of said vector component signals.
31. Apparatus according to claim 30 further including:
means for calculating the arctan of a ratio of the spatially averaged vector components for developing a signal representative of said velocity estimation.
32. Apparatus according to claim 31 further including:
means for calculating a turbulence estimation using said spatially averaged vector components.
33. Apparatus according to claim 29, further including:
means for repeatedly calculating velocity estimations so as to develope a two dimensional array of said velocity estimations over a given scan area; and
means for displaying said two dimensional array of velocity estimations so as to develope an image of velocity estimations within said given scan area.
34. Apparatus according to claim 29, further including:
means for temporally averaging the results of said spatial averaging means, for developing an estimation of the turbulence of said target.
35. Apparatus for developing an estimation of the velocity of a moving target, comprising:
means for transmitting groups of pulse signals towards said target;
means for receiving groups of echo signals caused by reflections of said transmitted groups of pulse signals from said target;
means for complex conjugate multiplication of a pair of temporally successive groups of said echo signals for developing vector component signals;
means for spatial averaging of said vector component signals; and
means for temporal averaging of the results of said spatial averagings for a successive pair of temporally successive groups of said echo signal for developing said estimation of the velocity of said target.
36. Apparatus according to claim 35 further including:
means for calculating a turbulence estimation from said spatially and temporally averaged vector components.
37. Apparatus according to claim 35, further including:
means for repeatedly calculating velocity estimations so as to develope a two dimensional array of said velocity estimations over a given scan area; and
means for displaying said two dimensional array of velocity estimations so as to develope an image of velocity estimations within said given scan area.
US07/120,577 1987-11-13 1987-11-13 Doppler velocity processing method and apparatus Expired - Lifetime US4800891A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US07/120,577 US4800891A (en) 1987-11-13 1987-11-13 Doppler velocity processing method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US07/120,577 US4800891A (en) 1987-11-13 1987-11-13 Doppler velocity processing method and apparatus

Publications (1)

Publication Number Publication Date
US4800891A true US4800891A (en) 1989-01-31

Family

ID=22391199

Family Applications (1)

Application Number Title Priority Date Filing Date
US07/120,577 Expired - Lifetime US4800891A (en) 1987-11-13 1987-11-13 Doppler velocity processing method and apparatus

Country Status (1)

Country Link
US (1) US4800891A (en)

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1990004793A1 (en) * 1988-10-19 1990-05-03 Urs Moser Process and device for measuring a two-dimensional reflective structure
US4979513A (en) * 1987-10-14 1990-12-25 Matsushita Electric Industrial Co., Ltd. Ultrasonic diagnostic apparatus
US5022400A (en) * 1988-02-08 1991-06-11 Walters Glenn A Large time bandwidth echographic signal processor
EP0447597A1 (en) * 1990-03-23 1991-09-25 Siemens Aktiengesellschaft Method and apparatus to avoid ambiguities in a pulse-doppler apparatus
US5058594A (en) * 1990-08-29 1991-10-22 Quantum Medical Systems, Incorporated Direct velocity estimator for ultrasound blood flow imaging
EP0475100A1 (en) * 1990-08-20 1992-03-18 Matsushita Electric Industrial Co., Ltd. Ultrasonic doppler blood flowmeter
US5215094A (en) * 1990-11-14 1993-06-01 Advanced Technology Laboratories, Inc. Ultrasonic flow velocity imaging systems with velocity image presistence
US5313946A (en) * 1991-01-09 1994-05-24 Hewlett-Packard Company Method and apparatus for the characterization of tissue or other structure
WO1994023421A1 (en) * 1993-03-26 1994-10-13 Siemens Medical Systems, Inc. Digital beamformer having multi-phase parallel processing
WO1994023422A1 (en) * 1993-03-26 1994-10-13 Siemens Medical Systems, Inc. Partial beamforming
DE4424506A1 (en) * 1993-07-29 1995-02-02 Siemens Ag Method of determining blood flow parameter values
US5386830A (en) * 1993-10-25 1995-02-07 Advanced Technology Laboratories, Inc. Ultrasonic pulsed doppler flow measurement system with two dimensional autocorrelation processing
US5494037A (en) * 1994-12-29 1996-02-27 Siemens Medical Systems, Inc. Ultrasonic doppler imager having a spatially smoothed control signal for an adaptive tissue rejection filter
US5513640A (en) * 1994-06-24 1996-05-07 Kabushiki Kaisha Toshiba Diagnostic ultrasound system
US5544128A (en) * 1994-07-05 1996-08-06 Siemens Medical Systems, Inc. Multi-beam digital beamforming method and apparatus
US5544658A (en) * 1995-09-18 1996-08-13 Siemens Medical Systems, Inc. Doppler ultrasound velocity estimation
US5555886A (en) * 1995-09-28 1996-09-17 Siemens Medical Systems, Inc. Apparatus and method for detecting blood vessel size and direction for doppler flow measurement system
US5562098A (en) * 1995-03-20 1996-10-08 Lifesigns Corporation Ultrasonic measurement of blood flow velocity independent of probe angle
US5623929A (en) * 1995-06-30 1997-04-29 Siemens Medical Systems, Inc. Ultrasonic doppler flow imaging method for eliminating motion artifacts
US5653234A (en) * 1995-09-29 1997-08-05 Siemens Medical Systems, Inc. Method and apparatus for adaptive spatial image filtering
US5724974A (en) * 1996-03-22 1998-03-10 Acuson Corporation Two-dimensional ultrasound display system
US6030344A (en) * 1996-12-04 2000-02-29 Acuson Corporation Methods and apparatus for ultrasound image quantification
US6093149A (en) * 1996-12-04 2000-07-25 Acuson Corporation Method and apparatus for setting the integration interval for time integrated surface integral in an ultrasound imaging system
US6177923B1 (en) 1994-12-30 2001-01-23 Acuson Corporation Imaging modality showing energy and velocity
WO2001097704A3 (en) * 2000-06-23 2002-03-07 Koninkl Philips Electronics Nv Automatic flow angle correction by ultrasonic vector processing
US20050033174A1 (en) * 2003-07-10 2005-02-10 Moehring Mark A. Doppler ultrasound method and apparatus for monitoring blood flow and hemodynamics
US20050075568A1 (en) * 1998-11-11 2005-04-07 Moehring Mark A. Doppler ultrasound method and apparatus for monitoring blood flow
US20060264759A1 (en) * 2005-05-20 2006-11-23 Moehring Mark A System and method for grading microemboli monitored by a multi-gate doppler ultrasound system
US20070016050A1 (en) * 2005-06-13 2007-01-18 Moehring Mark A Medical Doppler ultrasound system for locating and tracking blood flow
EP2609870A1 (en) * 2011-12-29 2013-07-03 Samsung Medison Co., Ltd. Providing turbulent flow information based on vector doppler in ultrasound system
EP2628449A3 (en) * 2012-02-17 2013-09-04 Samsung Electronics Co., Ltd. Ultrasound apparatus and method of generating ultrasound image
US20160273957A1 (en) * 2015-03-19 2016-09-22 National Instruments Corporation Machine Condition Monitoring Using Phase Adjusted Frequency Referenced Vector Averaging
US20170192088A1 (en) * 2015-12-31 2017-07-06 Herbert U Fluhler Least squares fit classifier for improved sensor performance
US20180088212A1 (en) * 2016-09-28 2018-03-29 Robert Bosch Gmbh Method for evaluating radar radiation, and radar apparatus

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4391148A (en) * 1980-02-29 1983-07-05 National Research Development Corporation Methods and apparatus for measuring the cross-sectional area of a duct and the volume flow rate of fluid in the duct
EP0214654A2 (en) * 1985-09-12 1987-03-18 Aloka Co. Ltd. Ultrasonic device for measuring acceleration of moving reflective member
US4671294A (en) * 1985-08-14 1987-06-09 Hewlett-Packard Company Pulsed doppler flow mapping apparatus
US4697594A (en) * 1985-08-21 1987-10-06 North American Philips Corporation Displaying a single parameter image

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4391148A (en) * 1980-02-29 1983-07-05 National Research Development Corporation Methods and apparatus for measuring the cross-sectional area of a duct and the volume flow rate of fluid in the duct
US4671294A (en) * 1985-08-14 1987-06-09 Hewlett-Packard Company Pulsed doppler flow mapping apparatus
US4697594A (en) * 1985-08-21 1987-10-06 North American Philips Corporation Displaying a single parameter image
EP0214654A2 (en) * 1985-09-12 1987-03-18 Aloka Co. Ltd. Ultrasonic device for measuring acceleration of moving reflective member

Cited By (54)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4979513A (en) * 1987-10-14 1990-12-25 Matsushita Electric Industrial Co., Ltd. Ultrasonic diagnostic apparatus
US5022400A (en) * 1988-02-08 1991-06-11 Walters Glenn A Large time bandwidth echographic signal processor
WO1990004793A1 (en) * 1988-10-19 1990-05-03 Urs Moser Process and device for measuring a two-dimensional reflective structure
EP0616231A2 (en) * 1988-10-19 1994-09-21 Institut für Biomedizinische Technik und Medizinische Informatik der Universität und ETH Zürich Process and device for measuring the current volume in a layer with reflecting structure
EP0616231A3 (en) * 1988-10-19 1994-10-12 Inst Biomedizinische Technik Process and device for measuring the current volume in a layer with reflecting structure.
EP0447597A1 (en) * 1990-03-23 1991-09-25 Siemens Aktiengesellschaft Method and apparatus to avoid ambiguities in a pulse-doppler apparatus
US5081994A (en) * 1990-03-23 1992-01-21 Siemens Aktiengesellschaft Method and apparatus for avoiding ambiguities in a pulse doppler apparatus
EP0475100A1 (en) * 1990-08-20 1992-03-18 Matsushita Electric Industrial Co., Ltd. Ultrasonic doppler blood flowmeter
US5220923A (en) * 1990-08-20 1993-06-22 Matsushita Electric Industrial Co., Ltd. Ultrasonic doppler blood flowmeter
US5058594A (en) * 1990-08-29 1991-10-22 Quantum Medical Systems, Incorporated Direct velocity estimator for ultrasound blood flow imaging
US5215094A (en) * 1990-11-14 1993-06-01 Advanced Technology Laboratories, Inc. Ultrasonic flow velocity imaging systems with velocity image presistence
US5313946A (en) * 1991-01-09 1994-05-24 Hewlett-Packard Company Method and apparatus for the characterization of tissue or other structure
US5369624A (en) * 1993-03-26 1994-11-29 Siemens Medical Systems, Inc. Digital beamformer having multi-phase parallel processing
WO1994023422A1 (en) * 1993-03-26 1994-10-13 Siemens Medical Systems, Inc. Partial beamforming
WO1994023421A1 (en) * 1993-03-26 1994-10-13 Siemens Medical Systems, Inc. Digital beamformer having multi-phase parallel processing
US5388079A (en) * 1993-03-26 1995-02-07 Siemens Medical Systems, Inc. Partial beamforming
DE4424506C2 (en) * 1993-07-29 2001-02-22 Siemens Ag Method for determining blood flow parameter values
DE4424506A1 (en) * 1993-07-29 1995-02-02 Siemens Ag Method of determining blood flow parameter values
US5386830A (en) * 1993-10-25 1995-02-07 Advanced Technology Laboratories, Inc. Ultrasonic pulsed doppler flow measurement system with two dimensional autocorrelation processing
US5513640A (en) * 1994-06-24 1996-05-07 Kabushiki Kaisha Toshiba Diagnostic ultrasound system
USRE38209E1 (en) * 1994-06-24 2003-08-05 Kabushiki Kaisha Toshiba Diagnostic ultrasound system
US5544128A (en) * 1994-07-05 1996-08-06 Siemens Medical Systems, Inc. Multi-beam digital beamforming method and apparatus
US5494037A (en) * 1994-12-29 1996-02-27 Siemens Medical Systems, Inc. Ultrasonic doppler imager having a spatially smoothed control signal for an adaptive tissue rejection filter
US6177923B1 (en) 1994-12-30 2001-01-23 Acuson Corporation Imaging modality showing energy and velocity
US5562098A (en) * 1995-03-20 1996-10-08 Lifesigns Corporation Ultrasonic measurement of blood flow velocity independent of probe angle
US5623929A (en) * 1995-06-30 1997-04-29 Siemens Medical Systems, Inc. Ultrasonic doppler flow imaging method for eliminating motion artifacts
US5544658A (en) * 1995-09-18 1996-08-13 Siemens Medical Systems, Inc. Doppler ultrasound velocity estimation
US5555886A (en) * 1995-09-28 1996-09-17 Siemens Medical Systems, Inc. Apparatus and method for detecting blood vessel size and direction for doppler flow measurement system
US5653234A (en) * 1995-09-29 1997-08-05 Siemens Medical Systems, Inc. Method and apparatus for adaptive spatial image filtering
US5724974A (en) * 1996-03-22 1998-03-10 Acuson Corporation Two-dimensional ultrasound display system
US5785655A (en) * 1996-03-22 1998-07-28 Acuson Corporation Two-dimensional ultrasound display system
US6093149A (en) * 1996-12-04 2000-07-25 Acuson Corporation Method and apparatus for setting the integration interval for time integrated surface integral in an ultrasound imaging system
US6030344A (en) * 1996-12-04 2000-02-29 Acuson Corporation Methods and apparatus for ultrasound image quantification
US6193664B1 (en) 1996-12-04 2001-02-27 Acuson Corporation Method and apparatus for ultrasound image quantification
US6241677B1 (en) 1996-12-04 2001-06-05 Acuson Corporation Method and apparatus for ultrasound image quantification
US6322511B1 (en) 1996-12-04 2001-11-27 Acuson Corporation Methods and apparatus for ultrasound image quantification
US6464640B1 (en) 1996-12-04 2002-10-15 Acuson Corporation Methods and apparatus for ultrasound imaging with automatic color image positioning
US6110118A (en) * 1996-12-04 2000-08-29 Acuson Corporation Method and apparatus for ultrasound image quantification
US20050075568A1 (en) * 1998-11-11 2005-04-07 Moehring Mark A. Doppler ultrasound method and apparatus for monitoring blood flow
US7537568B2 (en) * 1998-11-11 2009-05-26 Spentech, Inc. Doppler ultrasound method and apparatus for monitoring blood flow
WO2001097704A3 (en) * 2000-06-23 2002-03-07 Koninkl Philips Electronics Nv Automatic flow angle correction by ultrasonic vector processing
US7128713B2 (en) 2003-07-10 2006-10-31 Spentech, Inc. Doppler ultrasound method and apparatus for monitoring blood flow and hemodynamics
US20050033174A1 (en) * 2003-07-10 2005-02-10 Moehring Mark A. Doppler ultrasound method and apparatus for monitoring blood flow and hemodynamics
US20060264759A1 (en) * 2005-05-20 2006-11-23 Moehring Mark A System and method for grading microemboli monitored by a multi-gate doppler ultrasound system
US7771358B2 (en) 2005-05-20 2010-08-10 Spentech, Inc. System and method for grading microemboli monitored by a multi-gate doppler ultrasound system
US20070016050A1 (en) * 2005-06-13 2007-01-18 Moehring Mark A Medical Doppler ultrasound system for locating and tracking blood flow
US8162837B2 (en) 2005-06-13 2012-04-24 Spentech, Inc. Medical doppler ultrasound system for locating and tracking blood flow
EP2609870A1 (en) * 2011-12-29 2013-07-03 Samsung Medison Co., Ltd. Providing turbulent flow information based on vector doppler in ultrasound system
EP2628449A3 (en) * 2012-02-17 2013-09-04 Samsung Electronics Co., Ltd. Ultrasound apparatus and method of generating ultrasound image
US20160273957A1 (en) * 2015-03-19 2016-09-22 National Instruments Corporation Machine Condition Monitoring Using Phase Adjusted Frequency Referenced Vector Averaging
US20170192088A1 (en) * 2015-12-31 2017-07-06 Herbert U Fluhler Least squares fit classifier for improved sensor performance
US10871551B2 (en) * 2015-12-31 2020-12-22 Herbert U Fluhler Least squares fit classifier for improved sensor performance
US20180088212A1 (en) * 2016-09-28 2018-03-29 Robert Bosch Gmbh Method for evaluating radar radiation, and radar apparatus
US10514443B2 (en) * 2016-09-28 2019-12-24 Robert Bosch Gmbh Method for evaluating radar radiation, and radar apparatus

Similar Documents

Publication Publication Date Title
US4800891A (en) Doppler velocity processing method and apparatus
US5170792A (en) Adaptive tissue velocity compensation for ultrasonic Doppler imaging
US5662115A (en) Method for determining the velocity-time spectrum of blood flow
US5860930A (en) Energy weighted parameter spatial/temporal filter
US5971927A (en) Ultrasonic diagnostic apparatus for obtaining blood data
US6390980B1 (en) Spatial compounding with ultrasonic doppler signal information
US4993417A (en) Method and system for controlling ultrasound scanning sequence
KR100742466B1 (en) Method and apparatus for visualization of motion in ultrasound flow imaging using continuous data acquisition
US5876341A (en) Removing beam interleave effect on doppler spectrum in ultrasound imaging
US5445156A (en) Method for adaptively filtering doppler signals using a complex time domain filter
JP4627366B2 (en) Method and apparatus for motion visualization in ultrasonic flow imaging using packet data acquisition
US5188112A (en) Ultrasonic Doppler imaging systems with improved flow sensitivity
JP2003501177A (en) Simultaneous tissue and motion ultrasound diagnostic imaging
US4759375A (en) Ultrasonic doppler diagnostic apparatus
US4918605A (en) Method and system for detecting and processing ultrasonic doppler signals
US5058594A (en) Direct velocity estimator for ultrasound blood flow imaging
US4884448A (en) Ultrasonic doppler meter
JP2763126B2 (en) Color ultrasonic diagnostic equipment
US6544184B1 (en) Imaging with reduced artifacts for medical diagnostic ultrasound
US20070112269A1 (en) Clutter filtering with small ensemble lengths in ultrasound imaging
JP3281435B2 (en) Ultrasound Doppler diagnostic equipment
JP4155615B2 (en) Ultrasound imaging system
JP2788926B2 (en) Velocity evaluation method and device by Doppler method
JP2678124B2 (en) Ultrasound Doppler diagnostic device
JP2563656B2 (en) Ultrasonic Doppler imaging device

Legal Events

Date Code Title Description
AS Assignment

Owner name: SIEMENS MEDICAL LABORATORIES, INC., 2404 NORTH MAI

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST.;ASSIGNOR:KIM, JIN H.;REEL/FRAME:004809/0239

Effective date: 19871210

Owner name: SIEMENS MEDICAL LABORATORIES, INC.,CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:KIM, JIN H.;REEL/FRAME:004809/0239

Effective date: 19871210

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

AS Assignment

Owner name: SIEMENS QUANTUM, INC., WASHINGTON

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST.;ASSIGNOR:SIEMENS MEDICAL LABORATORIES, INC.;REEL/FRAME:006096/0687

Effective date: 19920414

FPAY Fee payment

Year of fee payment: 4

FPAY Fee payment

Year of fee payment: 8

FPAY Fee payment

Year of fee payment: 12