Practical cone beam image reconstruction using local regions-of-interest

• US 6,009,142 A
• Filed: 03/31/1998
• Issued: 12/28/1999
• Est. Priority Date: 03/31/1998
• Status: Expired due to Term

• Current Assignee: Siemens Medical Solutions USA Incorporated (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:
  US09/052,415
• Time in Patent Office: 637 Days
• Field of Search: 378/4, 378/15, 378/901
• US Class Current: 378/15
• CPC Class Codes: A61B 6/027 characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• 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 method and apparatus of operating a CT imaging apparatus having a cone beam radiation source and a 2D detector arrangement, for exactly reconstructing an image of a 3D region of interest (ROI) in an object. Measurement data acquired by the 2D detector is processed to develop sub-sets of Radon data. Each of the sub-sets of Radon data is targeted for reconstructing a corresponding local 2D ROI in a 2D parallel projection of the object. Accordingly, after a sub-set of Radon data is completely developed, it is subjected to a first inversion processing step for developing the corresponding local 2D ROI. Multiple ones of the local ROIs are then grouped together and subjected to a second inversion processing step to develop an image reconstruction of a part of the 3D ROI in the object.

Claims

1. A method of operating a CT imaging apparatus having a cone beam radiation source and a 2D detector arrangement, for reconstructing an image of a 3D region of interest (ROI) of an object, comprising the following steps:

operating the source and detector arrangement at a plurality of source positions along a scan path that encircles the region of interest in the object, so as to acquire a set of measurement data at each of the source positions that is representative of radiation attenuation caused by the object as the radiation source and the detector arrangement are operated along the scan path;

conversion processing the acquired measurement data to develop multiple sub-sets of Radon data, each of the sub-sets being targeted for reconstructing a corresponding local 2D ROI in a 2D parallel projection of the object; and

Radon inversion processing of multiple ones of the sub-sets that are grouped together, to develop an image reconstruction of a part of the 3D ROI in the object.

2. The method of claim 1, wherein said conversion processing develops said sub-sets of Radon data in a sampled spherical Radon space partitioned by a plurality of coaxial ψ

-planes, and each of said targeted local 2D ROI'"'"'s comprise an area on a given one of the ψ

-planes that is defined by a 2D parallel projection of the object onto the given ψ

-plane, and having upper and lower boundaries determined by a 2D parallel projection of the scan path that cuts through the object and onto the given ψ

-plane.

3. The method of claim 1, wherein said Radon inversion processing step comprises first and second inversion processing steps, said first step developing said local 2D ROI'"'"'s, and said second step developing said image reconstruction of said part of the object in a piecewise manner.

4. The method of claim 3, wherein the first inversion processing step comprises processing a portion of the Radon data so as to develop a 2D image on an image grid (z,r'"'"') of an intermediate space (r'"'"',ψ

, z), and after a number of local 2D ROI'"'"'s on a plurality of ψ

-planes have been subjected to the first inversion processing step so as to be sufficient to fill up at least one z-plane of image data in the intermediate space, subjecting the at least one z-plane of image data to the second inversion processing step.

5. The method of claim 3, wherein the conversion processing step comprises using a database of pre-calculated image processing information that identifies which of the acquired measurement data should be processed to develop Radon data for the local 2D ROI'"'"'s.

6. The method of claim 5, wherein the conversion processing step comprises using the database information to also identify how the acquired measurement data should be processed to develop Radon data for developing the local 2D ROI'"'"'s.

7. The method of claim 5, wherein the conversion processing step comprises using the database information as a list of detector weight factors, for directly developing Radon data in a region of support for the local 2D ROI'"'"'s in each of the plurality of ψ

-planes, by multiplication of the acquired measurement data with the detector weight factors.

8. The method of claim 5, wherein the conversion processing step comprises using the database information in a relative manner so that it is re-usable for calculating sub-sets of Radon data for the local 2D ROI'"'"'s in a plurality of the ψ

-planes.

9. The method of claim 5, wherein the conversion processing step comprises using the database information for calculating the Radon data sub-set for developing a local 2D ROI in one of the ψ

-planes, and then re-using the database information for calculating multiple sub-sets of the Radon data for developing the local 2D ROI'"'"'s in multiple ones of the plurality of ψ

-planes.

10. The method of claim 9, wherein the conversion processing step comprises using the database information useful for developing a local 2D ROI in one of the ψ

-planes, for calculating Radon data for successively adjacent local 2D ROI'"'"'s in the one ψ

-plane.

11. The method of claim 4, wherein said first and second inversion steps each comprise a filtered backprojection step.

12. The method of claim 1, wherein said Radon inversion processing includes a first step that develops a local 2D ROI on each of the ψ

-planes that includes a sub-area therein that spans a common portion of a plurality of the ψ

-planes.

13. The method of claim 12, wherein the sub-areas of the local 2D ROI'"'"'s that span a common portion of the ψ

-planes, are grouped together and subjected to a second inversion processing step for developing a image reconstruction of the object in a piecewise planar manner.

14. The method of claim 12, wherein said first inversion processing step develops a plurality of 2D local ROI'"'"'s which are adjacent one another on a given ψ

-plane.

15. The method of claim 14, wherein said first inversion processing step develops said adjacent 2D local ROI'"'"'s so that they are positioned on a given ψ

-plane in an overlapped manner, and only a central portion of each 2D local ROI is subjected to said second inversion processing step, so as to reduce the occurrence of artifacts that may result from subjecting portions of each 2D local ROI that is near it'"'"'s upper and lower boundaries to said second inversion processing step.

16. A scanning and data acquisition apparatus for three dimensional (3D) computerized tomography (CT) imaging of a region of interest (ROI) in an object, comprising:

an energy source for emitting energy in the form of a cone beam;

an area detector for detecting as measurement data cone beam energy after attenuation by passage through an object to be imaged;

a manipulator for movably positioning said source and detector along a scan path relative to the object, to cause irradiation of said object by said source at multiple source positions along said scan path about said object, said detector acquiring a corresponding set of measurement data at each said source position;

control means for defining the scan path as a plurality of positions on a predetermined geometric surface surrounding the object, with each of a plurality of planes that pass through the object and intersect the scan path in at least one point as well as the area detector, being useful for calculating Radon data at a given one of a plurality of Radon samples;

image reconstruction processing means including means for calculating Radon data for each of said Radon samples by processing multiple ones of the sets of measurement data, for developing image data; and

display means responsive to said image data for displaying a reconstructed image of the object;

wherein said image reconstruction processing means comprises;

a conversion processor for processing the acquired measurement data to develop multiple sub-sets of Radon data, each of the sub-sets targeted for reconstructing a corresponding local 2D ROI in a 2D parallel projection of the object, and

a Radon inversion processor for processing multiple ones of the sub-sets that are grouped together, to develop an image reconstruction of a part of the ROI in the object.

17. The scanning and data acquisition apparatus of claim 16, wherein said conversion processor develops said sub-sets of Radon data in a sampled spherical Radon space partitioned by a plurality of coaxial ψ

-planes, with each of said targeted local 2D ROI'"'"'s comprising an area on a given one of the ψ

-planes that is defined by a 2D parallel projection of the object onto the given ψ

-plane, and having upper and lower boundaries determined by a 2D parallel projection of the scan path that cuts through the object and onto the given ψ

-plane.

18. The scanning and data acquisition apparatus of claim 16, wherein said Radon inversion processor includes a first inversion processor for developing said local 2D ROI'"'"'s, and a second inversion processor for processing groups of said local 2D ROI'"'"'s to develop a piecewise planar image reconstruction of the object.

19. The scanning and data acquisition apparatus of claim 18, wherein the first inversion processor processes a portion of the Radon data so as to develop a 2D image on an image grid (z,r'"'"') of an intermediate space (r'"'"',ψ

, z), and after a number of local 2D ROI'"'"'s on a plurality of ψ

-planes have been developed so as to be sufficient to fill up at least one z-plane of image data in the intermediate space, said second inversion processor processes the at least one z-plane of image data to develop a planar image reconstruction of the object.

20. The scanning and data acquisition apparatus of claim 18, wherein the conversion processor includes a database of pre-calculated image processing information that identifies which of the acquired measurement data should be processed to develop Radon data for the local 2D ROI'"'"'s.

21. The scanning and data acquisition apparatus of claim 20, wherein said database information identifies how the acquired measurement data should be processed to develop Radon data for developing the local 2D ROI'"'"'s.

22. The scanning and data acquisition apparatus of claim 20, wherein said database information is a list of detector weight factors, which is used by said conversion processor for directly developing Radon data in a region of support for the local 2D ROI'"'"'s in each of the plurality of ψ

-planes, by multiplication of the acquired measurement data with the detector weight factors.

23. The scanning and data acquisition apparatus of claim 20, wherein said database information is used by said conversion processor in a relative manner for calculating sub-sets of Radon data for the local 2D ROI'"'"'s in a plurality of the ψ

-planes.

24. The scanning and data acquisition apparatus of claim 20, wherein said conversion processor uses the database information for calculating the Radon data sub-set for developing a local 2D ROI in one of the ψ

-planes, and then re-uses the database information for calculating multiple sub-sets of the Radon data for developing the local 2D ROI'"'"'s in multiple ones of the plurality of ψ

-planes.

25. The scanning and data acquisition apparatus of claim 24, wherein said conversion processor uses the database information for developing a local 2D ROI in one of the ψ

-planes, and then re-uses it for calculating Radon data for successively adjacent local 2D ROI'"'"'s in the one ψ

-plane.

26. The scanning and data acquisition apparatus of claim 18, wherein said first inversion processor develops a local 2D ROI on each of the ψ

-planes that includes a sub-area therein that spans a common portion of a plurality of the ψ

-planes.

27. The scanning and data acquisition apparatus of claim 26, wherein the sub-areas of the local 2D ROI'"'"'s that span a common portion of the ψ

-planes, are grouped together by a second inversion processor for developing a image reconstruction of the object in a piecewise planar manner.