Two-step radon inversion processing for .phi.-planes having local radon origins

• US 5,901,195 A
• Filed: 09/30/1997
• Issued: 05/04/1999
• Est. Priority Date: 09/30/1997
• Status: Expired due to Term

• Current Assignee: Siemens Corp. (Siemens AG)
• Original Assignee: Siemens Corporate Research Incorporated (Siemens AG)
• Inventors: Samarasekera, Supun, Sauer, Frank, Tam, Kwok
• Primary Examiner(s): Bruce, David Vernon
• Application Number:
  US08/940,324
• Time in Patent Office: 581 Days
• Field of Search: 378/4, 378/901, 382/131
• US Class Current: 378/4
• CPC Class Codes: A61B 6/027 characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
• CPC Class Codes: G06T 11/006 Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
• CPC Class Codes: G06T 2211/421 Filtered back projection [FBP]
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

A two-step 3D Radon inversion processor for processing Radon data on a set of vertically oriented co-axial φ-planes that partition Radon space, wherein each of the φ-planes is sampled during the processing so as to have its own independent local Radon origin. In accordance with further principles of the present invention, a 3D Radon data generator generates Radon data on a set of vertically oriented co-axial φ-planes that partition Radon space, wherein each of the φ-planes is sampled so as to have its own independent local Radon origin, and a two-step 3D Radon inversion processor independently processes each of the φ-planes. Any shift of the local Radon origins with respect to a global coordinate system is compensated for during the 3D Radon inversion processing. Compensation for the origin shift may be made during or after a first step of the Radon inversion processing.

Claims

1. A method of operating a 3D CT imaging apparatus having a cone beam radiation source and detector arrangement, for reconstructing an image of an object that is scanned by said source and detector arrangement, comprising the following steps:

operating said source and detector arrangement at a plurality of source positions along a scan path that traverses the object so as to acquire a corresponding plurality of sets of measurement data representative of radiation attenuation caused by said object;

processing the acquired measurement data so as to develop Radon data on a set of vertically oriented co-axial φ

-planes that partition a Radon space; and

reconstructing an image of the object by performing two 2D radon inversions of the Radon data, wherein each of the Radon space φ

-planes is sampled during a first of said inversion processing'"'"'s so as to have its own independent local Radon origin.

2. The method of claim 1, wherein said two-step Radon inversion processing includes a processing step that compensates for a shift of the local Radon origin of each φ

-plane with respect to a global coordinate system.

3. A method of operating a 3D CT imaging apparatus having a cone beam radiation source and detector arrangement, for reconstructing an image of an object that is scanned by said source and detector arrangement, comprising the following steps:

operating said source and detector arrangement at a plurality of source positions along a scan path that traverses the object so as to acquire a corresponding plurality of sets of measurement data representative of radiation attenuation caused by said object;

processing the acquired measurement data so as to develop Radon data on a set of vertically oriented co-axial φ

-planes that partition a Radon space, wherein each of the φ

-planes is sampled so as to have its own independent local Radon origin; and

reconstructing an image of the object by performing two-step 3D Radon inversion processing of the Radon data wherein each of the φ

-planes is independently processed.

4. The method of claim 3, wherein said two-step Radon inversion processing includes a first step that compensates for a shift of the local Radon origins with respect to a global coordinate system by backprojecting the Radon data from the φ

-planes onto φ

-plane sample grids (z,r'"'"') which are not shifted, i.e., are accordingly offset, from the local Radon origins.

5. The method of claim 3, wherein said two-step Radon inversion processing includes a first step that compensates for a shift of the local Radon origins with respect to a global coordinate system by backprojecting the Radon data onto φ

-plane sample grids (z,r'"'"') which are shifted along with the local Radon origins, and then, before performing a second step of the two-step inversion processing, shifting the backprojection results a corresponding amount toward the global origin to compensate for the previous shift in the local Radon origins.

6. The method of claim 4, wherein the backprojecting of the first step includes introducing a corresponding offset into coordinate variables that define the backprojecting.

7. The method of claim 6, wherein the shifting of the backprojection results is accomplished by interpolating the results onto a new grid which is aligned with the global coordinate system.