Ultrasound image texture control using adaptive speckle control algorithm

• US 5,841,889 A
• Filed: 12/29/1995
• Issued: 11/24/1998
• Est. Priority Date: 12/29/1995
• Status: Expired due to Fees

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Seyed-Bolorforosh, Mir Said
• Primary Examiner(s): Couso, Jose L.
• Assistant Examiner(s): Bella, Matthew C.
• Application Number:
  US08/581,668
• Time in Patent Office: 1,061 Days
• Field of Search: 382/128, 382/254, 382/274, 128/660.07, 128/661.08, 128/661.01, 600/443, 600/448, 600/463, 600/441
• US Class Current: 382/128
• CPC Class Codes: G01N 2291/0289 Internal structure, e.g. defects, grain size, texture
• CPC Class Codes: G01S 7/52046 Techniques for image enhancement involving transmitter or receiver image enhancement by image data processing G06T5/00

Abstract

A method for controlling the contrast resolution and the tissue texture of ultrasonic images employs a signal processing algorithm to adjust an image by changing the tissue texture (speckle grain size). The algorithm increases the image contrast by allocating a different speckle pattern to different grey scale levels. This signal processing algorithm is based on the Hilbert transform. It uses the property of the Hilbert transform to generate a quadrature component of the given image which is modulated by the speckle pattern as well as all other image information. The quadrature component is then squared before being selectively added to or subtracted from the square of the in-phase component of the given image to control the information in the image. The added or subtracted signals have different textures (or spatial frequency components). The arithmetic operation between the in-phase and quadrature phase components determines the image texture. If this operation is done selectively, then the contrast and the texture of the image can be selectively controlled.

Claims

1. A method for adjusting speckle in an image, comprising the steps of:

acquiring an input waveform S_in (x, y) representing an input image;

transforming said input waveform to form an in-phase component (Re(S_in (x, y)) and a quadrature phase component (Im(S_in (x, y)) of said input waveform; and

computing an output waveform S_out (x, y) representing a speckle-modulated region of said input image in accordance with the following equation;

space="preserve" listing-type="equation">S.sub.out (x,y)=|{(Re(S.sub.in (x,y)).sup.2 +M(Im(S.sub.in (x,y)).sup.2 }|.sup.0.5

wherein M is an adaptive speckle control index having a value in a range -1≦

M<

1.

2. The method as defined in claim 1, wherein the quadrature phase component is computed using the Hilbert transform.

3. The method as defined in claim 1, further comprising the step of concurrently displaying said speckle-modulated region of said input image in conjunction with a non-speckle-adjusted region of said input image.

4. The method as defined in claim 3, further comprising the steps of scan converting and interpolating said output waveform prior to displaying.

5. The method as defined in claim 1, wherein said computing step comprises the steps of:

computing a quantity (Re(S_in (x, y))² +(Im(S_in (x, y))² !^1/2 for each sample of said input waveform;

determining a size of said kernel in dependence on said computed quantity;

integrating said computed quantity over said kernel;

dividing the integrated quantity by the area of the kernel, resulting in a mean integrated value; and

selecting a value for said adaptive speckle control index M by inputting said mean integrated value into a mapping of M values.

6. The method as defined in claim 5, wherein the product of kernel area and said input sample data is constant.

7. The method as defined in claim 5, wherein the kernel size is proportional to the localized spatial frequency contents of the input image waveform.

8. The method as defined in claim 1, wherein said acquiring step comprises the step of beamforming.

9. A method for ultrasound imaging of an object, comprising the steps of:

transmitting an ultrasound beams along a scan line into an object to be imaged;

forming an acoustic data signal representing the amplitude of a net ultrasound beam reflected by said object from a range point along said scan line, said acoustic data signal beamforming step being performed for each of a multiplicity of range points along said scan line;

transforming each acoustic data signal to form an in-phase component and a quadrature phase component thereof;

for at least one of said transformed acoustic data signals, generating adaptive acoustic data signals which represent the square root of the absolute value of a sum of the square of said in-phase component plus the product of an adaptive speckle control index times the square of said quadrature phase components, wherein said adaptive speckle control index has a value ≧

-1<

+1; and

displaying a specular image of said adaptive acoustic data signals.

10. The method as defined in claim 9, wherein the quadrature phase component is computed using the Hilbert transform.

11. The method as defined in claim 9, further comprising the steps of:

computing an amplitude of an envelope of said in-phase component and said quadrature component for each of said acoustic data signals;

determining a size of a kernel in dependence on said computed envelope amplitude;

integrating said computed envelope amplitude over said kernel;

adjusting the kernel size to have a constant integrated value, but different mean integrated value; and

selecting a value for said adaptive speckle control index by inputting the result of said integrating step into a mapping of adaptive speckle control index values.

12. An apparatus for adjusting speckle in an image, comprising:

a beamformer for acquiring an input waveform S_in (x, y) representing an input image;

means for transforming said input waveform to form an in-phase component (Re(S_in (x, y)) and a quadrature phase component (Im(S_in (x, y)) of said input waveform;

processing means for computing an output waveform S_out (x, y) representing a speckle-modulated region of said input image in accordance with the following equation;

space="preserve" listing-type="equation">S.sub.out (x,y)=|{(Re(S.sub.in (x,y)).sup.2 +M(Im(S.sub.in (x,y)).sup.2 }|.sup.0.5

wherein M is an adaptive speckle control index having a value in a range -1≦

M<

1; and

means for displaying said speckle-modulated region of said input image in conjunction with a non-speckle-adjusted region of said input image.

13. The apparatus as defined in claim 12, wherein the quadrature phase component is computed using the Hilbert transform.

14. The apparatus as defined in claim 12, further comprising a scan converter and an interpolator for scan converting and interpolating said output waveform prior to display.

15. The apparatus as defined in claim 12, wherein said processing means comprise:

means for computing a quantity .(Re(S_in (x, y))² +(Im(S_in (x, y))² !^1/2 for each sample of said input waveform;

means for determining a size of said kernel in dependence on said computed quantity;

means for integrating said computed quantity over said kernel;

means for storing a mapping of M values; and

means for selecting a value for said adaptive speckle control index M by inputting the result of said integrating step into said mapping of M values.

16. The apparatus as defined in claim 15, wherein the product of kernel area and said input sample data is constant.

17. The apparatus as defined in claim 15, wherein the kernel size is proportional to the localized spatial frequency contents of the input image waveform.