Method of reducing noise in volume imaging

• US 6,862,334 B2
• Filed: 07/17/2002
• Issued: 03/01/2005
• Est. Priority Date: 07/19/2001
• Status: Expired due to Term

• Current Assignee: Koninklijke Philips Electronics N.V. (Koninklijke Philips N.V.)
• Original Assignee: Koninklijke Philips Electronics N.V. (Koninklijke Philips N.V.)
• Inventors: Thelissen, Guillaume Rudolf Petrus, Van Liere, Filips
• Primary Examiner(s): BRUCE, DAVID VERNON
• Application Number:
  US10/484,303
• Publication Number:
  US 20040199064A1
• Time in Patent Office: 958 Days
• Field of Search: 378/4, 378/8, 378/15, 378/19, 378/901
• US Class Current: 378/4
• CPC Class Codes: G06T 11/005 Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

In a method of visualizing a three-dimensional volume image data set (1), a two-dimensional image (15) is formed by projecting volume image data on a projection plane in that for each pixel (16) there is determined an intensity value which corresponds to the maximum value or minimum value of the mean intensity along a projection path (17) through the imaged volume and associated with the relevant pixel. In order to determine the mean intensity, averaging is performed over a plurality of intensity values which neighbor one another along the projection path (17). In order to achieve effective noise suppression as well as adequate along image contrast , averaging is performed over at least partly overlapping interval (18 to 22) along the projection path (17). The method offer the advantage that the user can choose the width of the averaging interval as an additional free parameter, so that image contrast and background noise can be influenced independently.

Claims

1. A method of visualizing a three-dimensional volume image data set (1) which contains a plurality of intensity values on a discrete spatial grid, in which method a two-dimensional image (15) is formed by projecting the volume image data on a projection plane in that for each pixel (16) there is determined an intensity value which corresponds to the maximum value or minimum value of the mean intensity along a projection path (17) which is associated with the relevant pixel and extends through the imaged volume, the mean intensity being calculated by averaging over a plurality of intensity values of the volume image data set (1) which neighbor one another along the projection path, characterized in that the averaging is carried out over at least partly overlapping intervals (18 to 22) along the projection path (17).

2. A method as claimed in claim 1, characterized in that the averaging is carried out by convolution of the intensity variation along the projection path with a weighting function which has a selectable width which is larger than the distance in space between two neighboring data points of the volume image data set.

3. A method as claimed in claim 1, characterized in that the intervals over which the averaging is performed overlap each time by one half.

4. An imaging diagnostic apparatus, notably a CT apparatus or an MR apparatus, for carrying out the method claimed in claim 1, which apparatus includes an imaging unit (27, 28, 31, 32) for the acquisition of coarse data of an object (30) to be examined and also includes a program-controlled reconstruction unit (33) which reconstructs volume image data from the coarse data, said volume image data consisting of a plurality of intensity values on a discrete spatial grid, and also forms a two-dimensional image (15) therefrom by projecting the volume image data on a projection plane, characterized in that projection of the volume image data is carried out by the reconstruction unit (33) in such a manner that for each pixel (16) there is determined an intensity value which corresponds to the maximum value or minimum value of the mean intensity along a projection path (17) which extends through the imaged volume and is associated with the relevant pixel, the mean intensity being calculated by averaging over at least partly overlapping intervals (18 to 22) along the projection path.

5. A computer program for carrying out the method claimed in claim 1, characterized in that it forms a two-dimensional image (15) from a three-dimensional volume image data set (1), consisting of a plurality of intensity values on a discrete spatial grid, by projecting the volume image data on a projection plane in that for each pixel (16) it determines an intensity value which corresponds to the maximum value or minimum value of the mean intensity along a projection path (17) which extends through the imaged volume and is associated with the relevant pixel, the mean intensity being calculated by averaging over at least partly overlapping intervals (18 to 22) along the projection path (17).

6. A computer program as claimed in claim 5, characterized in that it performs the averaging by convolution of the intensity variation along the projection path with a weighting function which has a width which can be selected by the user and is larger than the distance in space between two neighboring data points of the volume image data set.