Pre-calculated hitlist for reducing run-time processing of an exact cone beam reconstruction algorithm

• US 5,862,198 A
• Filed: 09/30/1997
• Issued: 01/19/1999
• Est. Priority Date: 09/30/1997
• 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): Porta, David P.
• Assistant Examiner(s): Bruce, David Vernon
• Application Number:
  US08/940,924
• Time in Patent Office: 476 Days
• Field of Search: 378/4, 378/901, 382/131, 600/425
• 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: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

A method and apparatus for reconstructing an image of an object using a three dimensional (3D) computerized tomographic (CT) imager having a cone beam radiation source and detector arrangement for acquiring measurement data. Before the acquisition of measurement data during an imaging operation of the apparatus, image reconstruction processing information dependant on the geometric parameters of the imager and required for processing of the acquired measurement data for developing contributions to the final reconstruction of the image, is pre-calculated and stored. During an imaging operation of the apparatus, the acquired measurement data is processed using the pre-calculated image reconstruction processing information for reconstructing the image.

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:

before the acquisition of measurement data during an imaging operation of the apparatus,

pre-calculating and storing image reconstruction processing information required for processing of the acquired measurement data to develop contributions to the final reconstruction of the image, which information is dependant on the geometric parameters of the imager and independent of the measurement data; and

during the imaging operation of the apparatus,

operating said source and detector arrangement at a plurality of source positions along a scan path so as to acquire a corresponding plurality of sets of said measurement data, and

processing the acquired measurement data using said pre-calculated image reconstruction processing information for developing said contributions and reconstructing the image.

2. The method of claim 1, wherein said pre-calculating and storing step comprises,

calculating as said image reconstruction processing information that information which is determined by predetermined geometric parameters of the scan path and the detector, and a desired sampling of the scan path and Radon space.

3. The method of claim 1, wherein said processing step develops numerical values for a plurality of sample points in Radon space in response to reconstruction information supplied thereto from said database for developing contributions to said image reconstruction, said database information being further sorted so as to supply to said processing step a subset of its information for causing said processing step to develop a subset of said sample points.

4. The method of claim 1, wherein said pre-calculating and storing step comprises,

organizing the pre-calculated image reconstruction processing information into a database.

5. The method of claim 4, wherein:

said operating step operates the source and detector arrangement at a plurality of spatially adjacent source positions along the scan path so as to cause acquisition of said corresponding plurality of sets of measurement data in a time sequential order, and

said processing step processes the acquired measurement data in the time sequential order that they are acquired, in conjunction with said database having said image reconstruction processing information sorted therein by the time sequential order of the acquisition of said sets of measurement data.

6. The method of claim 5, wherein said processing step develops numerical values for a plurality of sample points in Radon space in response to reconstruction information supplied thereto from said database for developing contributions to said image reconstruction, said database information being further sorted so as to supply to said processing step a subset of its information for causing said processing step to develop a subset of said sample points.

7. The method of claim 6, comprising the further step of:

controllably sorting said database by said sample points so as to only supply to said processing step reconstruction information which will reconstruct the image in a desired resolution.

8. The method of claim 6, comprising the further step of:

controllably sorting said database by said sample points so as to only supply to said processing step reconstruction information for a region of support in Radon space which is appropriate for the size of the object being imaged.

9. The method of claim 1, wherein said pre-calculating and storing step comprises,

calculating as said image reconstruction processing information that information which indicates the correspondence between sample points in Radon space and the source positions that develop measurement data that can be processed to contribute numerical values thereto.

10. The method of claim 9, wherein said pre-calculating and storing step includes,

defining a plurality of planes, each one of which will pass through a field of view of the imaging apparatus and intersect the scan path in at least one point as well as the detector, for determining the correspondence between sample points and the source positions from which the contributions thereto can be developed.

11. The method of claim 10, wherein said pre-calculating and storing step includes,

calculating as said image reconstruction processing information parameters that define line integral derivatives that need to be calculated in the acquired measurement data for contributing to the numerical value of said sample points.

12. The method of claim 10 wherein said image reconstruction processing information includes masking and line parameters.

13. 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:

before the acquisition of measurement data during an imaging operation of the apparatus,

pre-calculating and storing image reconstruction processing information required for processing of the acquired measurement data to develop contributions to the final reconstruction of the image, which information is stored in an organized format that indicates the correspondence between sample points in Radon space and the source positions that develop measurement data that can be processed to contribute numerical values thereto; and

during the imaging operation of the apparatus, operating said source and detector arrangement at a plurality of source positions along a scan path so as to acquire a corresponding plurality of sets of said measurement data, and

processing the acquired measurement data using said pre-calculated image reconstruction processing information for developing said contributions and reconstructing the image.

14. The method of claim 13, wherein said pre-calculating and storing step comprises,

organizing the pre-calculated image reconstruction processing information into a database.

15. The method of claim 14, wherein:

said operating step operates the source and detector arrangement at a plurality of spatially adjacent source positions along the scan path so as to cause acquisition of said corresponding plurality of sets of measurement data in a time sequential order, and

said processing step processes the acquired measurement data in the time sequential order that they are acquired, in conjunction with the organized format of said database information being sorted by the time sequential order of the acquisition of said sets of measurement data.

16. The method of claim 15, wherein said processing step develops numerical value contributions for a plurality of sample points in Radon space in response to reconstruction information supplied thereto, the organized format of said database information being further sorted so as to supply to said processing step a subset of its information for causing said processing step to develop a subset of said sample points.

17. The method of claim 16, comprising the further step of:

controllably sorting said database by said sample points so as to only supply to said processing step reconstruction information which will reconstruct the image in a desired resolution.

18. The method of claim 16, comprising the further step of:

controllably sorting said database by said sample points so as to only supply to said processing step reconstruction information for a region of support in Radon space which is appropriate for the size of the object being imaged.

19. A scanning and data acquisition apparatus for three dimensional (3D) computerized tomography (CT) imaging of an object in a field of view radially centered on a predetermined axis, comprising:

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

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

a manipulator (6) 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 (8) for defining the scan path as a plurality of spaced stages on a predetermined geometric surface surrounding the field of view, with each of a plurality of planes that pass through the field of view 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 sample points;

image reconstruction processing means (18) including means for calculating Radon data for each of said Radon sample points by processing multiple ones of the sets of measurement data, for developing image data; and

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

wherein said image reconstruction processing means (18) includes a database (21) of reconstruction processing information that is dependant on the geometric parameters of the imager, which information is precalculated and stored before acquisition of the measurement data, and then used by said image reconstruction processing means (18) during imaging operation of said apparatus for facilitating processing of the measurement data to calculate the Radon data.

20. Apparatus in accordance with claim 19, wherein said image reconstruction processing means (18) includes Radon derivative calculating means (20) responsive to the image reconstruction processing information stored in said database (21) for calculating Radon derivative data for each of said Radon sample points by processing line integral values from multiple ones of the sets of measurement data.

21. Apparatus in accordance with claim 19, wherein said image reconstruction processing means (18) includes Radon derivative data converting means (22) responsive to the image reconstruction processing information stored in said (database 21) and said Radon derivative data for calculating Radon data at equally spaced polar grid points.

22. Apparatus in accordance with claim 19, wherein said image reconstruction processing means (18) includes inverse transformation processing means (24) responsive to said Radon data at equally spaced polar grid points for performing a 3D inverse Radon transformation of said Radon data to and develop said image data.

23. Apparatus in accordance with claim 19, wherein said image reconstruction processing means (18) is response to said database (21) so as to develop numerical values for a plurality of sample points in Radon space in response to reconstruction information supplied thereto from said database.

24. Apparatus in accordance with claim 23, wherein said image reconstruction processing means (18) includes means for causing said database (21) to supply to said processing means (18) a subset of its information, thereby causing said processing means (18) to develop a subset of said sample points.

25. Apparatus in accordance with claim 23, wherein said manipulator positions the energy source (10) at a plurality of spatially adjacent source positions along the scan path so as to cause acquisition of said corresponding sets of measurement data in a time sequential order, and

said database (21) provides said image reconstruction processing information to said processing means in said time sequential order, thereby enabling said processing means to process the acquired measurement data in the time sequential order that they are acquired.