New! View global litigation for patent families

US6126605A - Ultrasound color flow display optimization by adjusting dynamic range - Google Patents

Ultrasound color flow display optimization by adjusting dynamic range Download PDF

Info

Publication number
US6126605A
US6126605A US09224342 US22434298A US6126605A US 6126605 A US6126605 A US 6126605A US 09224342 US09224342 US 09224342 US 22434298 A US22434298 A US 22434298A US 6126605 A US6126605 A US 6126605A
Authority
US
Grant status
Grant
Patent type
Prior art keywords
range
dynamic
data
memory
flow
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 - Fee Related
Application number
US09224342
Inventor
Michael J. Washburn
Gary E. MacLeod
Sean D. Lucas
David J. Muzilla
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.)
General Electric Co
Original Assignee
General Electric Co
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
Grant date

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
    • 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
    • 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
    • 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/52034Data rate converters
    • 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/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52071Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/58Testing, adjusting or calibrating the diagnostic device
    • A61B8/585Automatic set-up of the device

Abstract

An ultrasound color flow imaging system is programmed to optimize display images of power and velocity by automatically adjusting thresholds and data compression by using histograms and samplings of color flow data.

Description

CROSS-REFERENCE TO RELATED APPLICATIONS

Not applicable.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

Not applicable.

BACKGROUND OF THE INVENTION

This invention generally relates to ultrasound color flow Doppler imaging of fluid flow fields. In particular, the invention relates to a method and an apparatus for improving the display of such imaging.

Ultrasonic scanners for detecting blood flow based on the Doppler effect are well known. Such systems operate by actuating an ultrasonic transducer array to transmit ultrasonic waves into the object and receiving ultrasonic echoes backscattered from the object. In the measurement of blood flow characteristics, returning ultrasonic waves are compared to a frequency reference to determine the frequency shift imparted to the returning waves by flowing scatterers such as blood cells. This frequency, i.e., phase shift, translates into the velocity of the blood flow. The blood velocity is calculated by measuring the phase shift from firing to firing at a specific range gate.

The change or shift in backscattered frequency increases when blood flows toward the transducer and decreases when blood flows away from the transducer. Color flow images are produced by superimposing a color image of the velocity of moving material, such as blood, over a black and white anatomical B-mode image. Typically, color flow mode displays hundreds of adjacent sample volumes simultaneously, all laid over a B-mode image and color-coded to represent each sample volume's velocity.

Typically, color flow processors estimate blood flow velocity, blood flow power, and blood flow variance. Typically, color flow data is used to modify the color of a region of interest on a display screen. The user selects the type of data used for the display. The modes typically available are power only, velocity only or velocity and variance combined.

In current ultrasound scanners, various color flow display parameters are either fixed with no user selectability or are preset to some specific setting and can only be changed if action is taken by the user, one parameter at a time. This limits image quality and user productivity for any given application and scanning situation. There is a need for a scanner in which these same parameters can all be automatically adjusted at the same time to optimize image quality related to color flow display for a specific scanning situation, thus increasing user productivity.

In the color flow power mode of operation, known ultrasound scanners typically provide a color flow dynamic range based on a compression curve preselected at the factory depending on the type of scanning application. For example, one dynamic range based on one compression curve is selected for scanning of the kidney, whereas another dynamic range based on another compression curve is selected for scanning of the carotid artery. Frequently, the actual scan data has a dynamic range different from the range upon which the compression curve is based. As a result, the dynamic range of the display is less than optimal. Accordingly, there is a need for a color flow ultrasound scanner which can automatically adjust for changes in the dynamic range of the received signals.

BRIEF SUMMARY OF THE INVENTION

This invention is useful in an ultrasound imaging system generating color flow signals in response to ultrasound signals backscattered from a subject under study. This aspect of the invention displays images in response to the color flow signals by storing first memory values in response to the color flow signals, preferably in a digital memory. A dynamic range compression scheme is determined based on an analysis of the first memory values, preferably by a logic unit. Second memory values are generated based on the dynamic range compression scheme and a color flow image is displayed in response to the second memory values, preferably by a color display unit.

By using the foregoing techniques, the display of an ultrasound imaging device can be automatically adjusted for changes in the dynamic range of the received signals.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic block diagram showing the signal processing chain for a conventional color flow and B-mode ultrasound imaging system.

FIG. 2 is a schematic block diagram showing the mid processor color flow apparatus illustrated in FIG. 1.

FIG. 3 is a schematic block diagram showing the mid processor B-mode apparatus illustrated in FIG. 1.

FIG. 4 is a schematic block diagram showing additional details of portions of the system illustrated in FIG. 1.

FIG. 5 is a graph illustrating one possible form of data stored in the cine memory shown in FIG. 4.

FIG. 6 is a graph showing another possible form of data stored in the cine memory shown in FIG. 4.

FIG. 7 is a graph showing another possible form of data resulting from use of the preferred form of the invention.

DETAILED DESCRIPTION OF THE INVENTION

Referring to FIG. 1 the basic signal processing chain for a color flow and gray scale imaging system comprises an ultrasound transducer array 2, which is activated to transmit pulse sequences comprising tone bursts of length P which are fired repeatedly at a pulse repetition frequency (PRF) which typically is in the kilohertz range. The pulse sequences, including burst lengths P, are different for the color flow and B-mode processing. For color flow imaging, P may be 4 to 8 cycles, and the tone bursts are focused at the same transmit focal position with the same transmit characteristics.

A series of color flow transmit firings focused at the same transmit focal position are referred to as a "packet". Each transmit beam propagates through the object being scanned and is reflected by ultrasound scatterers in the object.

The return RF signals are detected by the transducer elements and received by the respective receive channels in the beamformer 4. The beamformer sums the delayed channel data and outputs in a beam summed signal which is demodulated into in-phase and quadrature (I/Q) signal components by a demodulator 6. The B-mode I, Q outputs from demodulator 6 are transmitted to a mid processor 8G for gray scale B-mode processing, and the color flow I, Q outputs from demodulator 6 are transmitted to a mid-processor 8C for color processing.

FIG. 2 illustrates mid-processor 8C. The I/Q signal components from demodulator 6 are stored in a corner turner memory 7, whose purpose is to buffer data from possibly interleaved firings and output the data as vectors of points across firings at a given range cell. Data is received in "fast time", or sequentially down range (along a vector) for each firing. The output of the corner turner memory is reordered into "slow time", or sequentially by firing for each range cell. The resultant "slow time" I/Q signal samples are passed through a wall filter 9 which rejects any clutter corresponding to stationary or very slow-moving tissue. The filtered outputs are then fed into a parameter estimator 11, which converts the range cell information into the intermediate autocorrelation parameters N, D and R(O). N and D are the numerator and denominator for the autocorrelation equation, as shown below: ##EQU1## where Ii and Qi are the demodulated, basebanded input data for firing i, and M is the number of firings in the packet. R(O) is approximated as a finite sum over the number of firings in a packet, as follows: ##EQU2##

A processor converts N and D into a magnitude and phase for each range cell. The equations used are as follows:

|R(T)|=√N.sup.2 +D.sup.2          (5) ##EQU3##

The parameter estimator processes the magnitude and phase values into signals having values representing estimates of power, velocity and turbulence or variance which are transmitted on conductors 11A, and 11B and 11C, respectively. The phase is used to calculate the mean Doppler frequency, which is proportional to the velocity as shown below; R(O) and |R(T)| (magnitude) are used to estimate the turbulence.

The mean Doppler frequency in hertz is obtained from the phase of N and D and the pulse repetition from T: ##EQU4##

The mean velocity is calculated using the Doppler shift equation below. Since θ, the angle between the flow direction and the sampling direction, is not known, cos θ is assumed to be 1.0. ##EQU5## Preferably, the parameter estimator does not calculate the mean Doppler frequency as an intermediate output, but calculates v directly from the phase output of the processor using a look-up table.

The turbulence may be calculated in the time domain as a second-order series expansion of the variance of the mean Doppler frequency. The time domain expression for turbulence involves calculating the zero-lag and one-lag autocorrelation functions, R(O) and R(T) respectively. The exact autocorrelation functions are approximated by finite sums over the known data in the number of firings in a packet: ##EQU6##

The mean value signal θ (R(T)) is an estimate of the mean Doppler frequency shift of the flowing reflectors, which in turn is proportional to the mean blood flow velocity. The variance signal σ2 indicates the frequency spread of the flow signal component of the baseband echo signal. This value is indicative of flow turbulence, since laminar flow has a very narrow range of velocities, while turbulent flow is a mixture of many velocities. To indicate the strength of the signal from the flowing reflectors, the signal R(O) indicates the amount of the returned power in the Doppler-shifted flow signal.

The signal power on conductor 11A is passed through a data compression module 13 which compresses the data according to families of data compression curves. A different family of curves can be provided for different scanning applications. For example, one family of curves is provided for renal scanning, while another family of curves is provided for carotid artery scanning. Typically, there are about three curves per family. The dynamic range of the signals is changed according to the curve used for the data compression. The curves in each family are arranged in order of increasing dynamic range. Controller 26 sets the default curve when a user selects the scan application. The dynamic range controls the range of intensities or lumens created on display 18.

Referring to FIG. 3, gray scale B-mode mid-processor 8G includes an envelope detector 10 for forming the envelope of the beamsummed receive signal by computing the quantity (I2 +Q2)1/2. The envelope of the signal undergoes some additional B-mode processing, such as logarithmic compression (block 12 in FIG. 3), to form display data which is output to the scan converter 14 (FIG. 1).

Referring again to FIG. 1, the color flow estimates and gray scale display data are sent to the scan converter 14, which converts the data into X-Y format for video display. The scan-converted frames are passed to a video processor 16, which basically maps the video data to a display color map and gray scale image frames for video display. The image frames are then sent to the video monitor 18 for display. Typically, for color images, either velocity or power are displayed alone or velocity is displayed in conjunction with either power or turbulence. System control is centered in a host computer (not shown), which accepts operator inputs through an operator interface (e.g., a keyboard) and in turn controls the various subsystems.

In general, for B-mode gray scale images, the display data is converted by the scan converter 14 into X-Y format for video display. The scan-converted frames are passed to the video processor 16, which maps the video data to a gray scale or mapping for video display. The gray scale image frames are then sent to the video monitor 18 for display.

The images displayed by the video monitor 18 are produced from an image frame of data in which each datum indicates the intensity or brightness of a respective pixel in the display. An image frame may, e.g., comprise a 256×256 data array in which each intensity datum is an 8-bit binary number that indicates pixel brightness. The brightness of each pixel on the display monitor 18 is continuously refreshed by reading the value of its corresponding element in the data array in a well-known manner. Each pixel has an intensity value which is a function of the backscatter cross section of a respective sample volume in response to interrogating ultrasonic pulses and the gray map employed.

Referring to FIG. 4, system control is centered in a master controller or host computer 26, which accepts operator inputs through an operator interface (not shown) and in turn controls the various subsystems. The master controller 26 also generates the system timing and control signals. The master controller 26 comprises a central processing unit (CPU) 30 and a random access memory 32. A keyboard 29 is used to enter data into CPU 30. The CPU 30 has read only memory incorporated therein for storing routines used in constructing gray and color maps based on acquired raw data.

The scan converter 14 comprises an acoustic line memory 22 and an X-Y memory 24. The B-mode and color mode intensity data stored in polar coordinate (R-θ) sector format in acoustic line memory 22 is transformed to appropriately scaled Cartesian coordinate pixel display data, which is stored in X-Y memory 24. The color data is stored in memory locations 24C, and the gray scale data is stored in memory locations 24G. The scan-converted frames are passed to video processor 16, which maps the data to a gray map for video display. The gray scale image frames are then sent to the video monitor for display.

Successive frames of acoustic sample data are stored in cine memory 28 on a first-in, first-out basis. Color frames are stored in memory locations 28C, and gray scale frames are stored in memory locations 28G. In the color region of interest, for every word of color data corresponding to a display pixel, there is a corresponding word of B-mode gray scale data corresponding to that pixel. The cine memory is like a circular image buffer that runs in the background, continually capturing acoustic sample data that is displayed in real time to the user. When the user freezes the system, the user has the capability to view acoustic sample data previously captured in cine memory.

The CPU 30 controls the XY memory 24 and the cine memory 28 via the system control bus 34. In particular, the CPU 30 controls the flow of raw data from the XY memory 24 to the video processor 16 and to the cine memory 28 and from the cine memory to the video processor 16 and to the CPU 26 itself. The CPU also loads the gray maps and color maps into the video processor.

Image frames are collected in cine memory 28 on a continuous basis. The cine memory 28 provides resident digital image storage for single image review and multiple image loop review and various control functions. The region of interest displayed during single-image cine replay is that used during the image's acquisition. The cine memory also acts as a buffer for transfer of images to digital archival devices (not shown) via the master controller 26.

The CPU 30 has random access memory for storing routines used in acquiring a raw data histogram, determining the end points of a new gray map input range, constructing a new gray map based on the end points of the new gray map input range, comparing the slope and gain of the new gray map to predetermined slope and gain limits, and if either limit is exceeded, reconstructing the new gray map to conform to the limit or limits.

In accordance with the preferred embodiments of the invention, the contrast of the ultrasound images is adjusted by the master controller 26 by creating a mapping of raw acoustic sample data into adjusted gray and color map values. First, the master controller 26 retrieves one or more image frames of raw data from the X-Y memory 24 or from the cine memory 28, storing that raw data in memory 32. The CPU 30 then compiles a histogram of the number of acoustic samples having an amplitude or value within each of a multiplicity of prescribed ranges or bins for the retrieved image frames of raw data.

According to the preferred embodiment, a color flow auto display processing mode is initiated by the user through keyboard 29 (FIG. 4) and can then be re-initiated by the user for updating of post-processing parameters or turned off altogether as the scanning situation changes.

The preferred embodiment uses the above-described B-mode gray scale and color flow scan data to optimize image quality of the color display. A composite histogram (histogram of the data over several frames) and/or a single frame histogram are constructed from the cine memory 28 data for color flow and/or B-mode by controller 26. Algorithms then are applied to the histogram results by controller 26 to determine how to properly adjust various parameters for a specific scanning situation or application.

By operating keyboard 29, the user may select the power only mode of color flow display (i.e., the PDI mode). According to the preferred embodiment, in the PDI mode, several discrete families of color flow dynamic range selections and data for corresponding compression curves are entered into memory 32 (FIG. 4). There is one family for each type of scanning application. For example, one family of dynamic range selections is used for ultrasound examination of the kidney, whereas another family of dynamic range selections is used for ultrasound examination of the carotid artery. In each family, there are three dynamic ranges available for automatic selection by the system. Each dynamic range is controlled by a different compression curve defined by digital data in memory 32.

The user enters the type of application on keyboard 29. In response, controller 26 presets the middle value dynamic range selection and corresponding compression curve which represents the typical dynamic range encountered when scanning the application selected by the user. The lower dynamic range in the family provides less dynamic range and the higher dynamic range in the family provides increased dynamic range compared to the middle dynamic range setting. Using the families of dynamic ranges, the auto dynamic range selection algorithm part of the auto color flow display processing mode automates dynamic range selection.

First, n frames of PDI color flow data are collected from cine memory 28C which represents the amplitude of the flow data present in the color flow region of interest (ROI) using the current preset dynamic range setting and compression curve. The n frames of data are required to account for flow pulsatility. Controller 26 executes the auto dynamic range selection algorithm which generates a first composite histogram of the data.

In the event the preset dynamic range is too low, a substantial percent of the data may be saturated at the maximum output value of 127. That is, the values are clustered in a range which is too high to create an optimum image on display monitor 18. Such a condition results in a first histogram of the type shown in FIG. 5. The auto dynamic range selection algorithm analyzes the first histogram, and, if at least x percent of the data is saturated at the maximum output value of 127 (i.e., condition 1 which is illustrated in FIG. 5), then the dynamic range is flagged as being too low, and the next higher dynamic range and corresponding compression curve stored in memory 32 are selected.

In the event the present dynamic range is too high, all of the data have values substantially below the maximum output value of 127. That is, all of the values are clustered in a range which is too low to create an optimum image on display monitor 18. Such a condition results in a first histogram of the type shown in FIG. 6. The auto dynamic range selection algorithm analyzes the first histogram, and, if less than t percent of the data occurs at output values between z and the maximum value of 127 (i.e., condition 2 which is illustrated in FIG. 6), then the dynamic range is flagged as being too high, and the next lower dynamic range and corresponding compression curve stored in memory 32 are selected.

Based on the newly selected dynamic range and compression curve, the auto dynamic range selection algorithm calculates a second histogram based on the same data used to calculate the first histogram. The original dynamic range and compression curve are used to decompress the original data. The original data then is applied to the newly selected dynamic range and compression curve to create new compressed data in memory 32. The second histogram calculated from the new compressed data using the new dynamic range typically shows a spread of values more appropriate for optimum viewing on display monitor 18. Such an exemplary second histogram may be of the type shown in FIG. 7.

Alternatively, if more than three dynamic ranges and three compression curves per family are provided, the process may be repeated as new dynamic ranges and compression curves are selected. If the first histogram showed that the dynamic range was too low (i.e., condition 1 shown in FIG. 5), the process of selecting dynamic ranges and compression curves, decompressing the data, compressing the data using the newly selected curve and calculating second and subsequent histograms is repeated until y percent of the data or less is determined to be saturated (the condition shown in FIG. 7) or until the highest dynamic range setting is reached, effectively spreading the data out over more of the lower output values. If the original histogram showed that the dynamic range was too high (i.e., condition 2 shown in FIG. 6), the process of selecting dynamic ranges and compression curves, decompressing the data, compressing the data using the newly selected curve and calculating second and subsequent histograms is repeated until q percent of the data or more is determined to be between z and 127 (the condition shown in FIG. 7) or until the lowest dynamic range setting is reached, effectively spreading the data out over more of the higher output values.

As a second alternative, the first histogram may be used to calculate a new compression curve not previously stored which more nearly spreads the actual data over the dynamic range available for display. The new compression curve then is used to recompress the original data.

If neither condition 1 nor condition 2 is met, then the current dynamic range and compression curve are maintained.

As an alternative to the preceding steps, controller 26 forms a histogram of the PDI data before any dynamic range compression is applied or uncompresses the data and forms a histogram. Then the histogram of uncompressed PDI data is analyzed to determine its statistics, and an optimal dynamic range compression scheme is calculated. Various dynamic range compression curves such as logarithmic, cube root, S-curve, or others are optimally applied across the data based on the statistics of the histogram. This alternative preferably is implemented by a digital signal processor (DSP).

The user can choose to re-activate the color flow auto display processing mode, causing new data to be formed into a histogram, based on the current dynamic range setting, and causing the auto dynamic range selection algorithm to be re-employed. Or the user can turn off the color flow auto display processing mode, causing the current dynamic range setting to be maintained until the user changes the selection manually.

The embodiments described here can be extended to automatically adjust other post-processing parameters such as power thresholds, wall filter cutoffs, baseline shifts, and velocity scales. The same basic idea of collecting B and/or color flow frames of cine data would be applied and associated algorithms would be employed to determine exactly how to adjust the particular post processing parameter.

The foregoing preferred embodiments have been disclosed for the purpose of illustration. Variations and modifications of the concept of the invention will be readily apparent to persons skilled in the art. All such variations and modifications are intended to be encompassed by the claims set forth hereinafter.

Claims (21)

What is claimed is:
1. In an ultrasound imaging system generating color flow signals in response to ultrasound signals backscattered from a subject under study, improved apparatus for displaying images in response to the color flow signals comprising in combination:
a memory connected to store first memory values in response to the color flow signals;
a logic unit connected to determine a dynamic range compression scheme based on an analysis of the first memory values and to generate second memory values based on the dynamic range compression scheme; and
a display connected to display a color flow image in response to the second memory values.
2. Apparatus, as claimed in claim 1, wherein the first memory values are altered from the values of the color flow signals based on a second dynamic range compression scheme.
3. Apparatus, as claimed in claim 1, wherein the logic unit is connected to create the dynamic range compression scheme based on an algorithm which analyzes the first memory values.
4. Apparatus, as claimed in claim 1, wherein the first memory values equal the values of the color flow signals.
5. Apparatus, as claimed in claim 1, wherein the memory stores compression data defining first and second dynamic range compression schemes, wherein the logic unit determines the dynamic range compression scheme by selecting one of the first or second stored dynamic range compression schemes and wherein the first memory values are altered based on the selected dynamic range compression scheme.
6. Apparatus, as claimed in claim 5, wherein the first dynamic range compression scheme comprises a first type and the second dynamic ranges comprises a second type different from the first type.
7. Apparatus, as claimed in claim 5, wherein the first memory values are altered from the values of the color flow signals based on a third dynamic range compression scheme.
8. Apparatus, as claimed in claim 7, wherein the logic unit is connected to analyze the first memory values by decompressing the first memory values, analyzing the decompressed first memory values, selecting one of the first or second stored dynamic range compression schemes, altering the decompressed first memory values based on the selected dynamic range compression scheme to generate second memory values and enabling storage of the second memory values.
9. Apparatus, as claimed in claim 1, wherein the logic unit is connected to analyze the first memory values by generating a histogram of the first memory values.
10. Apparatus, as claimed in claim 1, wherein the first memory values represent power estimates calculated in response to the backscattered signals.
11. In an ultrasound imaging system generating color flow signals in response to ultrasound signals backscattered from a subject under study, an improved method for displaying images in response to the color flow signals comprising the steps of:
storing first memory values in response to the color flow signals;
determining a dynamic range compression scheme based on an analysis of the first memory values;
generating second memory values based on the dynamic range compression scheme; and
displaying a color flow image in response to the second memory values.
12. A method, as claimed in claim 11, wherein the step of storing the first memory values comprises the step of altering the values of the color flow signals based on a second dynamic range compression scheme.
13. A method, as claimed in claim 11, wherein the step of determining a dynamic range compression scheme comprises the step of creating the dynamic range compression scheme by analyzing the first memory values.
14. A method, as claimed in claim 11, wherein the first memory values equal the values of the color flow signals.
15. A method, as claimed in claim 11, wherein the step of storing comprises the step of storing compression data defining first and second dynamic range compression schemes, wherein the step of determining a dynamic range compression scheme comprises the step of selecting one of the first or second stored dynamic range compression schemes and wherein the step of generating second memory values comprises the step of altering the values of the color flow signals based on the selected dynamic range compression scheme.
16. A method, as claimed in claim 15, wherein the first dynamic range compression scheme comprises a first type and the second dynamic ranges comprises a second type different from the first type.
17. A method, as claimed in claim 15, wherein the step of storing comprises the step of altering from the values of the color flow signals based on a third dynamic range compression scheme.
18. A method, as claimed in claim 17, wherein the step of selecting comprises the steps of:
decompressing the first memory values; and
analyzing the decompressed first memory values.
19. A method, as claimed in claim 18, wherein the step of generating comprises the steps of:
altering the decompressed first memory values based on the selected dynamic range compression scheme to generate second memory values; and
enabling storage of the second memory values.
20. A method, as claimed in claim 11, wherein the step of determining a dynamic range compression scheme comprises the step of analyzing the first memory values by generating a histogram of the first memory values.
21. A method, as claimed in claim 11, wherein the first memory values represent power estimates calculated in response to the backscattered signals.
US09224342 1998-12-31 1998-12-31 Ultrasound color flow display optimization by adjusting dynamic range Expired - Fee Related US6126605A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US09224342 US6126605A (en) 1998-12-31 1998-12-31 Ultrasound color flow display optimization by adjusting dynamic range

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US09224342 US6126605A (en) 1998-12-31 1998-12-31 Ultrasound color flow display optimization by adjusting dynamic range
PCT/US1999/028920 WO2000041004A1 (en) 1998-12-31 1999-12-06 Ultrasound color flow display optimization by adjusting dynamic range
JP2000592669A JP2002534189A (en) 1998-12-31 1999-12-06 Optimization of the ultrasonic color flow display by adjusting the dynamic range
EP19990967216 EP1060412A1 (en) 1998-12-31 1999-12-06 Ultrasound color flow display optimization by adjusting dynamic range
US09472048 US6379306B1 (en) 1998-12-31 1999-12-27 Ultrasound color flow display optimization by adjusting dynamic range including remote services over a network

Publications (1)

Publication Number Publication Date
US6126605A true US6126605A (en) 2000-10-03

Family

ID=22840246

Family Applications (2)

Application Number Title Priority Date Filing Date
US09224342 Expired - Fee Related US6126605A (en) 1998-12-31 1998-12-31 Ultrasound color flow display optimization by adjusting dynamic range
US09472048 Expired - Fee Related US6379306B1 (en) 1998-12-31 1999-12-27 Ultrasound color flow display optimization by adjusting dynamic range including remote services over a network

Family Applications After (1)

Application Number Title Priority Date Filing Date
US09472048 Expired - Fee Related US6379306B1 (en) 1998-12-31 1999-12-27 Ultrasound color flow display optimization by adjusting dynamic range including remote services over a network

Country Status (4)

Country Link
US (2) US6126605A (en)
EP (1) EP1060412A1 (en)
JP (1) JP2002534189A (en)
WO (1) WO2000041004A1 (en)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6379306B1 (en) * 1998-12-31 2002-04-30 General Electric Company Ultrasound color flow display optimization by adjusting dynamic range including remote services over a network
US20020186868A1 (en) * 2001-06-12 2002-12-12 Steinar Bjaerum Ultrasound color characteristic mapping
US6500122B1 (en) * 1999-03-30 2002-12-31 Ge Medical Systems Global Technology Company, Llc Method and apparatus for remotely positioning region of interest in image
US6592522B2 (en) * 2001-06-12 2003-07-15 Ge Medical Systems Global Technology Company, Llc Ultrasound display of displacement
US20030187354A1 (en) * 2002-04-01 2003-10-02 David Rust System and method of dynamic automatic sensing of available dynamic range
US20040102706A1 (en) * 2001-08-28 2004-05-27 Donald Christopher Automatic optimization of doppler display parameters
US20040208364A1 (en) * 2003-04-15 2004-10-21 Jamal Haque System and method for image segmentation
US20070038083A1 (en) * 2005-08-12 2007-02-15 Siemens Medical Solutions Usa, Inc. Automatic velocity scale identification for medical diagnostic ultrasound
US20080097212A1 (en) * 2006-08-16 2008-04-24 Siemens Medical Solutions Usa, Inc. Automatic identification of orientation in medical diagnostic ultrasound
US20090069675A1 (en) * 2007-09-11 2009-03-12 Siemens Medical Solutions Usa, Inc. Adaptive motion imaging in medical diagnostic ultrasound imaging
US8145833B1 (en) * 2009-10-30 2012-03-27 CSR Technology Holdings Inc. Compression solution for embedded devices and systems
US20120220877A1 (en) * 2011-02-28 2012-08-30 Kabushiki Kaisha Toshiba Data compression method and data compression apparatus
WO2013046089A1 (en) 2011-09-30 2013-04-04 Koninklijke Philips Electronics N.V. Ultrasound system with dynamically automated doppler flow settings as a sample volume is moved
WO2013046087A1 (en) 2011-09-30 2013-04-04 Koninklijke Philips Electronics N.V. Ultrasound system with automated doppler flow settings.
US20140093150A1 (en) * 2012-03-09 2014-04-03 Seno Medical Instruments, Inc Statistical mapping in an optoacoustic imaging system
US9733119B2 (en) 2011-11-02 2017-08-15 Seno Medical Instruments, Inc. Optoacoustic component utilization tracking
US9730587B2 (en) 2011-11-02 2017-08-15 Seno Medical Instruments, Inc. Diagnostic simulator
US9743839B2 (en) 2011-11-02 2017-08-29 Seno Medical Instruments, Inc. Playback mode in an optoacoustic imaging system
US9814394B2 (en) 2011-11-02 2017-11-14 Seno Medical Instruments, Inc. Noise suppression in an optoacoustic system

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6530887B1 (en) * 1996-12-24 2003-03-11 Teratech Corporation Ultrasound probe with integrated electronics
US6669633B2 (en) 1999-06-22 2003-12-30 Teratech Corporation Unitary operator control for ultrasonic imaging graphical user interface
US9402601B1 (en) 1999-06-22 2016-08-02 Teratech Corporation Methods for controlling an ultrasound imaging procedure and providing ultrasound images to an external non-ultrasound application via a network
US6969352B2 (en) * 1999-06-22 2005-11-29 Teratech Corporation Ultrasound probe with integrated electronics
US6542626B1 (en) * 1999-11-05 2003-04-01 General Electric Company Method and apparatus for adapting imaging system operation based on pixel intensity histogram
US8360973B2 (en) * 2000-11-29 2013-01-29 L'oreal Process for acquiring scanned image data relating to an external body portion and/or a product applied thereto
US8292811B2 (en) * 2003-03-20 2012-10-23 Siemens Medical Solutions Usa, Inc. Advanced application framework system and method for use with a diagnostic medical ultrasound streaming application
US6932767B2 (en) * 2003-03-20 2005-08-23 Siemens Medical Solutions Usa, Inc. Diagnostic medical ultrasound system having a pipes and filters architecture
JP4381028B2 (en) * 2003-05-08 2009-12-09 株式会社東芝 The ultrasonic diagnostic apparatus
JPWO2004107981A1 (en) * 2003-06-03 2006-07-20 株式会社日立メディコ The ultrasonic diagnostic apparatus
US20050043620A1 (en) * 2003-08-20 2005-02-24 Siemens Medical Solutions Usa, Inc. Diagnostic medical ultrasound system communication network architecture and method
JP2007505672A (en) * 2003-09-17 2007-03-15 コニンクリユケ フィリップス エレクトロニクス エヌ.ブイ. Iterative inspection report of
JP5203605B2 (en) * 2004-06-18 2013-06-05 株式会社日立メディコ The ultrasonic diagnostic apparatus
US7004906B1 (en) * 2004-07-26 2006-02-28 Siemens Medical Solutions Usa, Inc. Contrast agent imaging with agent specific ultrasound detection
US20080242988A1 (en) * 2007-03-28 2008-10-02 Tetsuya Yoshida Ultrasound diagnostic apparatus, ultrasound image processing apparatus, and ultrasound image processing method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5235984A (en) * 1992-03-30 1993-08-17 Hewlett-Packard Company On-line acoustic densitometry tool for use with an ultrasonic imaging system
US5347600A (en) * 1989-07-05 1994-09-13 Interated Systems, Inc. Method and apparatus for compression and decompression of digital image data
US5718229A (en) * 1996-05-30 1998-02-17 Advanced Technology Laboratories, Inc. Medical ultrasonic power motion imaging
US5732705A (en) * 1995-09-12 1998-03-31 Kabushiki Kaisha Toshiba Ultrasound diagnostic apparatus
US5735797A (en) * 1996-12-30 1998-04-07 General Electric Company Method and apparatus for combining topographic flow power imagery with a B-mode anatomical imagery
US5993392A (en) * 1996-11-05 1999-11-30 Atl Ultrasound, Inc. Variable compression of ultrasonic image data with depth and lateral scan dimensions
US6017309A (en) * 1998-12-31 2000-01-25 Washburn; Michael J. Ultrasound color flow display optimization by adjusting color maps

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5456255A (en) * 1993-07-12 1995-10-10 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus
US5492125A (en) * 1995-02-10 1996-02-20 University Of Washington Ultrasound signal processing apparatus
US5873830A (en) * 1997-08-22 1999-02-23 Acuson Corporation Ultrasound imaging system and method for improving resolution and operation
JP3862838B2 (en) * 1997-11-26 2006-12-27 株式会社東芝 The ultrasonic diagnostic apparatus
US5997478A (en) * 1998-02-03 1999-12-07 Acuson Corporation Ultrasound system and method for facilitating a reproducible ultrasound imaging environment
US6126605A (en) * 1998-12-31 2000-10-03 General Electric Company Ultrasound color flow display optimization by adjusting dynamic range

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5347600A (en) * 1989-07-05 1994-09-13 Interated Systems, Inc. Method and apparatus for compression and decompression of digital image data
US5235984A (en) * 1992-03-30 1993-08-17 Hewlett-Packard Company On-line acoustic densitometry tool for use with an ultrasonic imaging system
US5732705A (en) * 1995-09-12 1998-03-31 Kabushiki Kaisha Toshiba Ultrasound diagnostic apparatus
US5718229A (en) * 1996-05-30 1998-02-17 Advanced Technology Laboratories, Inc. Medical ultrasonic power motion imaging
US5993392A (en) * 1996-11-05 1999-11-30 Atl Ultrasound, Inc. Variable compression of ultrasonic image data with depth and lateral scan dimensions
US5735797A (en) * 1996-12-30 1998-04-07 General Electric Company Method and apparatus for combining topographic flow power imagery with a B-mode anatomical imagery
US6017309A (en) * 1998-12-31 2000-01-25 Washburn; Michael J. Ultrasound color flow display optimization by adjusting color maps

Cited By (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6379306B1 (en) * 1998-12-31 2002-04-30 General Electric Company Ultrasound color flow display optimization by adjusting dynamic range including remote services over a network
US6500122B1 (en) * 1999-03-30 2002-12-31 Ge Medical Systems Global Technology Company, Llc Method and apparatus for remotely positioning region of interest in image
US20020186868A1 (en) * 2001-06-12 2002-12-12 Steinar Bjaerum Ultrasound color characteristic mapping
US6592522B2 (en) * 2001-06-12 2003-07-15 Ge Medical Systems Global Technology Company, Llc Ultrasound display of displacement
US7245746B2 (en) 2001-06-12 2007-07-17 Ge Medical Systems Global Technology Company, Llc Ultrasound color characteristic mapping
US20040102706A1 (en) * 2001-08-28 2004-05-27 Donald Christopher Automatic optimization of doppler display parameters
US20030187354A1 (en) * 2002-04-01 2003-10-02 David Rust System and method of dynamic automatic sensing of available dynamic range
US6666824B2 (en) * 2002-04-01 2003-12-23 Koninklijke Philips Electronics N.V. System and method of dynamic automatic sensing of available dynamic range
US7298898B2 (en) * 2003-04-15 2007-11-20 Honeywell International Inc. System and method for image segmentation
US20040208364A1 (en) * 2003-04-15 2004-10-21 Jamal Haque System and method for image segmentation
US7887484B2 (en) 2005-08-12 2011-02-15 Siemens Medical Solutions Usa, Inc. Automatic velocity scale identification for medical diagnostic ultrasound
US20070038083A1 (en) * 2005-08-12 2007-02-15 Siemens Medical Solutions Usa, Inc. Automatic velocity scale identification for medical diagnostic ultrasound
US20080097212A1 (en) * 2006-08-16 2008-04-24 Siemens Medical Solutions Usa, Inc. Automatic identification of orientation in medical diagnostic ultrasound
US20090227871A2 (en) * 2006-08-16 2009-09-10 Siemens Medical Solutions Usa, Inc. Automatic Identification of Orientation in Medical Diagnostic Ultrasound
US8047991B2 (en) 2006-08-16 2011-11-01 Siemens Medical Solutions Usa, Inc. Automatic identification of orientation in medical diagnostic ultrasound
US8098910B2 (en) 2007-09-11 2012-01-17 Siemens Medical Solutions Usa, Inc. Adaptive motion imaging in medical diagnostic ultrasound imaging
US20090069675A1 (en) * 2007-09-11 2009-03-12 Siemens Medical Solutions Usa, Inc. Adaptive motion imaging in medical diagnostic ultrasound imaging
US8423711B1 (en) * 2009-10-30 2013-04-16 Cambridge Silicon Radio Limited Compression solution for embedded devices and systems
US8145833B1 (en) * 2009-10-30 2012-03-27 CSR Technology Holdings Inc. Compression solution for embedded devices and systems
US20120220877A1 (en) * 2011-02-28 2012-08-30 Kabushiki Kaisha Toshiba Data compression method and data compression apparatus
US8666132B2 (en) * 2011-02-28 2014-03-04 Kabushiki Kaisha Toshiba Data compression method and data compression apparatus
WO2013046087A1 (en) 2011-09-30 2013-04-04 Koninklijke Philips Electronics N.V. Ultrasound system with automated doppler flow settings.
WO2013046089A1 (en) 2011-09-30 2013-04-04 Koninklijke Philips Electronics N.V. Ultrasound system with dynamically automated doppler flow settings as a sample volume is moved
US9733119B2 (en) 2011-11-02 2017-08-15 Seno Medical Instruments, Inc. Optoacoustic component utilization tracking
US9730587B2 (en) 2011-11-02 2017-08-15 Seno Medical Instruments, Inc. Diagnostic simulator
US9743839B2 (en) 2011-11-02 2017-08-29 Seno Medical Instruments, Inc. Playback mode in an optoacoustic imaging system
US9814394B2 (en) 2011-11-02 2017-11-14 Seno Medical Instruments, Inc. Noise suppression in an optoacoustic system
US20140093150A1 (en) * 2012-03-09 2014-04-03 Seno Medical Instruments, Inc Statistical mapping in an optoacoustic imaging system
US9330452B2 (en) * 2012-03-09 2016-05-03 Seno Medical Instruments, Inc. Statistical mapping in an optoacoustic imaging system
US9836838B2 (en) * 2012-03-09 2017-12-05 Seno Medical Instruments, Inc. Statistical mapping in an optoacoustic imaging system

Also Published As

Publication number Publication date Type
WO2000041004A1 (en) 2000-07-13 application
JP2002534189A (en) 2002-10-15 application
US6379306B1 (en) 2002-04-30 grant
EP1060412A1 (en) 2000-12-20 application

Similar Documents

Publication Publication Date Title
US6102861A (en) Method and apparatus for three-dimensional ultrasound imaging using surface-enhanced volume rendering
US6176830B1 (en) Method and system for pre-determining spectral doppler user parameters
US5016641A (en) Spectral interpolation of ultrasound Doppler signal
US6755787B2 (en) Medical diagnostic ultrasound system and method for versatile processing
US6368277B1 (en) Dynamic measurement of parameters within a sequence of images
US4848354A (en) Method and apparatus for investigating a circulatory system in living biological structures
US6102864A (en) Three-dimensional ultrasound imaging of velocity and power data using average or median pixel projections
US5735797A (en) Method and apparatus for combining topographic flow power imagery with a B-mode anatomical imagery
US6733455B2 (en) System and method for adaptive clutter filtering in ultrasound color flow imaging
US6139501A (en) Coincident tissue and motion ultrasonic diagnostic imaging
US5282471A (en) Ultrasonic imaging system capable of displaying 3-dimensional angiogram in real time mode
US6210335B1 (en) Acoustic flash to increase penetration
US6296612B1 (en) Method and apparatus for adaptive wall filtering in spectral Doppler ultrasound imaging
US4790321A (en) Method of displaying stream lines of an inhomogeneous flowing medium and a device therefor
US5105817A (en) Ultrasonic bloodstream imaging apparatus
US6530885B1 (en) Spatially compounded three dimensional ultrasonic images
US6547732B2 (en) Adaptive image processing for spatial compounding
US6221020B1 (en) System and method for providing variable ultrasound analyses in a post-storage mode
US6048312A (en) Method and apparatus for three-dimensional ultrasound imaging of biopsy needle
US5107841A (en) Maximum entropy velocity estimator for ultrasonic flow imaging system
US6071240A (en) Method and apparatus for coherence imaging
US5609155A (en) Energy weighted parameter spatial/temporal filter
US4800891A (en) Doppler velocity processing method and apparatus
US5429137A (en) Acoustic scan conversion method and apparatus for velocity flow
US4930513A (en) Two dimensional processing of pulsed Doppler signals

Legal Events

Date Code Title Description
AS Assignment

Owner name: GENERAL ELECTRIC COMPANY, NEW YORK

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:WASHBURN, MICHAEL J.;MACLEOD, GARY E.;LUCAS, SEAN D.;ANDOTHERS;REEL/FRAME:009835/0704

Effective date: 19990311

FPAY Fee payment

Year of fee payment: 4

REMI Maintenance fee reminder mailed
LAPS Lapse for failure to pay maintenance fees
FP Expired due to failure to pay maintenance fee

Effective date: 20081003