Exact region of interest cone beam imaging using 3D backprojection

• US 5,926,521 A
• Filed: 03/31/1998
• Issued: 07/20/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: Tam, Kwok
• Primary Examiner(s): Bruce, David Vernon
• Application Number:
  US09/052,281
• Time in Patent Office: 476 Days
• Field of Search: 378/4, 378/15, 378/901
• 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: 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 scanning and data acquisition method and apparatus for three dimensional (3D) computerized tomography (CT) imaging of an object, wherein a plurality of line integral derivatives for a respective plurality of line segments L formed in cone beam projection data acquired on a detector at a plurality of source positions, are calculated. The extent of line segments L in the data acquired at each of the source positions is determined by a mask formed by cone beam projections onto the plane of the detector of portions of the source scan path that are above and below the source position that acquired the cone beam data in which the line integral derivatives for line segments L are being calculated. The line integral derivatives calculated for a plurality of the line segments L are then backprojected onto a 2D space corresponding to the plane of the detector. Finally, the results of the 2D backprojecting step are 3D backprojected into a 3D space, thereby reconstructing a 3D image of the object.

Claims

1. A method for performing three dimensional computerized tomography imaging of a region-of-interest (ROI) in an object using a cone beam source of radiation energy, comprising the steps of:

defining a source scanning trajectory as a scan path traversed by the source;

using the cone beam source, fixed relative to an area detector with both source and detector movably positioned relative to the object, to apply radiation energy towards the object from a plurality of source positions along the scan path as said source traverses the scan path, and cause said area detector to acquire cone beam projection data corresponding to respective portions of the object at each of said source positions;

calculating a line integral derivative for each of a plurality of line segments L formed in the cone beam projection data acquired at each of said source positions, the extent of said line segments L in the cone beam data acquired at a given one of said source positions being determined by a mask formed by cone beam projections onto the plane of the detector of portions of the source scan path that are above and below the given one of said source positions;

2D backprojecting the line integral derivatives calculated for a plurality of said line segments L onto a 2D space corresponding to the plane of the detector; and

3D backprojecting the results of the 2D backprojecting step into a 3D space, thereby reconstructing a 3D image of the ROI in the object.

2. The method of claim 1, wherein the calculating step calculates line integral derivatives for a plurality of line segments L spaced along each of a plurality of lines (r, θ

), each line having a different angular position, θ

, with respect to the coordinate system of the detector.

3. The method of claim 1, wherein the calculating step uses a mask formed for each source position that has a top curve and a bottom curve, with the extent of said line segments L for which the line integral derivatives are calculated in the cone beam projection data being bound between the top curve and the bottom curve of the mask.

4. The method of claim 3, wherein the source scanning trajectory is defined as a spiral scan path that surrounds the ROI in the object, and is connected to a first scan path circle at a top portion of the ROI in the object, and a second scan path circle at a bottom portion of the ROI in the object.

5. The method according to claim 4, wherein at source positions of the scan path at the top of the ROI in the object, the top curve of the mask includes a horizontal line portion formed from a projection of the first scan path circle, and the bottom curve of the mask is formed from a projection of the spiral scan path.

6. The method according to claim 4, wherein at source positions of the scan path at the bottom of the ROI in the object, the top curve of the mask is formed from a projection of the spiral scan path, and the bottom curve of the mask includes a horizontal line portion formed from a projection of the second scan path circle.

7. The method according to claim 5, wherein the top curve is formed from an intersection of a projection of the first scan path circle and a projection of the spiral scan path.

8. The method according to claim 6, wherein the bottom curve is formed from an intersection of a projection of the second scan path circle and a projection of the spiral scan path.

9. The method of claim 1, wherein after completion of the 2D backprojecting step, a 1D d/dt filtering step is performed on the results in said 2D space along lines in the direction of the scan path.

10. The method of claim 9, wherein the 3D backprojecting step comprises performing a weighted 3D backprojection of the 1D d/dt filtering results at discrete points along said lines onto a plurality of sample points in a 3D object volume.

11. The method of claim 10, wherein the weight assigned to each sample point in the 3D object volume is the inverse of the square of the distance between each sample point and a virtual location for the radiation source in the 3D volume.

12. The method of claim 1, wherein the 2D backprojecting step comprises 2D backprojecting the calculated line integral derivatives onto a virtual 2D space having a height that is greater than the height of said area detector.

13. The method of claim 12, wherein the height of said virtual 2D space is sufficient to cover the entire field-of-view of the ROI in the object.

14. Apparatus for performing three dimensional computerized tomography imaging of a region-of-interest (ROI) in an object using a cone beam source of radiation energy, comprising:

a source of cone beam radiation energy;

a manipulator for providing a source scanning trajectory as a scan path that encircles the ROI in the object and causes the source and detector to traverse the scan path;

means for causing the source to apply radiation energy towards the object from a plurality of source positions along the scan path as said source traverses the scan path, said area detector acquiring cone beam projection data corresponding to respective portions of the object at each of said source positions; and

an image reconstruction processor for,

calculating a line integral derivative for each of a plurality of line segments L formed in the cone beam projection data acquired at each of said source positions, the extent of said line segments L in the cone beam data acquired at a given one of said source positions being determined by a mask calculated by said image reconstruction processor that corresponds to the shape of cone beam projections onto the plane of the detector of portions of the source scan path that are above and below the given one of said source positions,

2D backprojecting the line integral derivatives calculated for a plurality of said line segments L onto a 2D space corresponding to the plane of the detector, and

3D backprojecting the 2D backprojection results into a 3D space, thereby reconstructing a 3D image of the ROI in the object.

15. The apparatus of claim 14, wherein the mask calculated for each source position that has a top curve and a bottom curve, with the extent of said line segments L for which the line integral derivatives are calculated in the cone beam projection data being bound between the top curve and the bottom curve of the mask.

16. The apparatus of claim 15, wherein the manipulator defines the source scanning trajectory as a spiral scan path the surrounds the ROI in the object, and is connected to a first scan path circle at a top portion of the ROI in the object, and a second scan path circle at a bottom portion of the ROI in the object.

17. The apparatus of claim 16, wherein the mask calculated at source positions of the scan path at the top of the ROI in the object has a top curve that includes a horizontal line portion formed from a projection of the first scan path circle, and the bottom curve of the mask is formed from a projection of the spiral scan path.

18. The apparatus of claim 16, wherein the mask calculated at source positions of the scan path at the bottom of the ROI in the object has a top curve that is formed from a projection of the spiral scan path, and a bottom curve that includes a horizontal line portion formed from a projection of the second scan path circle.

19. The apparatus of claim 17, wherein the top curve is formed from an intersection of a projection of the first scan path circle and a projection of the spiral scan path.

20. The apparatus of claim 18, wherein the bottom curve is formed from an intersection of a projection of the second scan path circle and a projection of the spiral scan path.

21. The apparatus of claim 14, wherein said image reconstruction processor performs a 1D d/dt filtering on the 2D backprojection results in said 2D space along lines extending in the direction of the scan path, before beginning the 3D backprojecting.

22. The apparatus of claim 21, wherein said image reconstruction processor performs the 3D backprojecting by performing a weighted 3D backprojection of the 1D d/dt filtering results at discrete points along said lines onto a plurality of sample points in a 3D object volume.

23. The apparatus of claim 14, wherein said image reconstruction processor performs the 2D backprojecting by 2D backprojecting the calculated line integral derivatives onto a virtual 2D space having a height that is greater than the height of said area detector.

24. The apparatus of claim 23, wherein the height of said virtual 2D space covers the entire field-of-view of the ROI in the object.