Data visualization enhancement through removal of dominating structures

• US 5,832,134 A
• Filed: 11/27/1996
• Issued: 11/03/1998
• Est. Priority Date: 11/27/1996
• Status: Expired due to Fees

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Avinash, Gopal Biligeri, Alyassin, Abdalmajeid Musa
• Primary Examiner(s): ROGERS, SCOTT A
• Application Number:
  US08/757,181
• Time in Patent Office: 706 Days
• Field of Search: 382/257, 382/256, 382/260, 382/262, 382/274, 382/300, 382/299, 382/131, 382/132, 382/128
• US Class Current: 382/257
• CPC Class Codes: A61B 6/5258 involving detection or reduction of artifacts or noise
• CPC Class Codes: G06T 2207/10072 Tomographic images
• CPC Class Codes: G06T 2207/30101 Blood vessel; Artery; Vein; Vascular
• CPC Class Codes: G06T 5/70 Denoising; Smoothing
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G06T 7/11 Region-based segmentation
• CPC Class Codes: G06T 7/155 involving morphological operators

Abstract

The present invention provides a system for performing fast segmentation and image processing for enhanced three dimensional (3-D) visualization of a subject. The present invention automatically extracts unwanted dominant features from images while preserving the desired information. The present invention works especially well with Computed Tomography Angiograms (CTA) for viewing vasculature of a subject'"'"'s head. Segmenting and removing dominant structures from image data permits visualization techniques such as maximum intensity projection (MIP), surface rendering and volume rendering, to provide enhanced vessel visualization.

Claims

1. A method of enhancing image data G1(x,y,z) comprising the steps of:

a) subsampling original image data G1 into lower resolution data g1(x,y,z);

b) thresholding g1 to create a threshold mask m1 with values of m1 being `1` where values of g1 are above the threshold, and zero when they are not;

c) determining connected regions of mask m1 and saving them as mask m2;

d) dilating mask m2 by a predetermined amount to create a dilated mask m3(x,y,z);

e) determining a transition region m4(x,y,z) being the difference between dilated masks m3 and threshold m2;

f) applying a filter to the region defined by m4 to create transition region replacement data g3(x,y,z);

g) interpolating m2, m4, g3 into M2, M4, G3 of original resolution, respectively;

h) determining a mean pixel intensity μ

from the image data G1;

i) assembling image data G4(x,y,z) according to;

##EQU7## j) displaying image data G4 as an enhanced image.

2. The method of enhancing images of claim 1 wherein the step of filtering comprises the step of applying an adaptive median filter to the region defined by m4 to create transition region replacement data g3(x,y,z).

3. The system for enhancing image data G1(x,y,z) comprising:

a) a subsampling device for subsampling original image data G1 into lower resolution data g1(x,y,z);

b) a thresholding device coupled to the subsampling device, for thresholding g1 to create a threshold mask m1 with values of m1 being `1` where values of g1 are above the threshold, and zero when they are not;

c) a connectivity device coupled to the thresholding device for receiving mask m1 and for determining connected regions of mask m1 and saving them as mask m2;

d) a dilation device coupled to the connectivity device for dilating mask m2 by a predetermined amount to create a dilated mask m3(x,y,z);

e) a transition region device coupled to the dilation device and the thresholding device for determining a transition region m4(x,y,z) being the difference between dilated masks m3 and threshold m2;

f) a replacement device, coupled to the subsampling device, for creating an average pixel intensity μ

, and creating a replaced data set g2(x,y,z) having value μ

replacing values of g1 where the dominant feature was (where m2=1), with all other locations being unchanged;

g) a filter coupled to the transition region device and the replacement device for replacing data values in the region defined by m4 with filtered values of g2 to create transition region replacement data g3(x,y,z) according to;

##EQU8## h) an interpolation device coupled to the filter, dilation, thresholding devices, for interpolating m2, m4, g3 into M2, M4, G3 of original resolution, respectively;

i) an image generator coupled to the interpolation device, for assembling image data G4(x,y,z) according to;

##EQU9## j) a display device coupled to the image generation device, for displaying image data G4 as an enhanced image.

4. The system for enhancing images of claim 3 wherein the filter comprises an adaptive median filter which filters data in the region defined by m4 to create transition region replacement data g3(x,y,z).