Method and apparatus for eliminating boundary errors in cone beam imaging

• US 5,748,697 A
• Filed: 12/20/1996
• Issued: 05/05/1998
• Est. Priority Date: 12/20/1996
• Status: Expired due to Term

• Current Assignee: Siemens AG
• Original Assignee: Siemens Corporate Research Incorporated (Siemens AG)
• Inventors: Tam, Kwok C.
• Primary Examiner(s): Porta, David P.
• Assistant Examiner(s): Bruce, David Vernon
• Application Number:
  US08/771,401
• Time in Patent Office: 501 Days
• Field of Search: 378/4, 378/8, 378/15, 378/19, 378/901
• US Class Current: 378/19
• 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/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 accordance with the principles of the present invention, image reconstruction errors which occur when using an area detector having a plurality of rows of detector elements defining a height for the detector that is smaller than the cone beam image, can be avoided by correlating an amount of orthogonal translation applied to line segments L extending across the detector towards top and bottom edges thereof that are used for calculating Radon derivative data, with the spacing between adjacent rows of the detector elements that are at the top and bottom edges of the detector. In a preferred embodiment the detector comprises a plurality of M rows of detector elements centered between the top and bottom edges of the detector, and an additional N rows of detector elements adjacent the M rows, at both of the top and bottom edges of the detector (where N may equal 1).

Claims

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

applying cone beam energy from a cone beam source to at least a portion of the object;

defining a source scanning trajectory as a 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 scan about the object;

specifying the source scanning trajectory as a spiral path defining a plurality of spaced stages on a predetermined geometric surface surrounding the field of view such that each plane passing through the field of view intersects the scanning trajectory in at least one point and intersects the area detector along a line segment L that extends towards top and bottom edges of the area detector, the area detector comprising a plurality of detector elements arranged in an array of rows and columns and having a predetermined number M of rows of detector elements between the top and bottom edges that extend sufficiently along a direction generally parallel to the predetermined axis so as to span at least two consecutive stages of the spiral path having the largest spacing therebetween, plus an extra N rows of detector elements adjacent both of the top and bottom edges of the detector, where N is greater than or equal to 1;

scanning at a plurality of positions along the source scanning trajectory to cause the detector elements of said area detector to acquire cone beam projection data corresponding to respective portions of the object;

calculating Radon derivative data by processing line integral values from cone beam projection data in adjacent detector line segments L₁ and L₂ that are parallel to line segment L and translated in an orthogonal direction therefrom by an amount which is less than or equal to N times the spacing between said extra N rows of detector elements contained in said area detector, where N is greater than of equal to 1; and

reconstructing an image of the object using said Radon derivative data.

2. The method of claim 1, wherein N=1.

3. The method of claim 1, wherein said reconstructing step includes performing an inverse Radon transformation of said Radon derivative data.

4. The method of claim 1, wherein said calculating step is repeated so as to develop Radon derivative data for all points in the Radon space of said object in the field of view.

5. The method of claim 1, wherein said calculating step is repeated for a plurality of closely spaced line segments L developed as a result of the intersection of said detector with a respective plurality of rotated planes that pass through the field of view and also intersect the scanning trajectory in at least one point.

6. The method of claim 1, wherein said calculating step is repeated at each of said plurality of positions along the source scanning trajectory provided by said scanning step.

7. The method of claim 4, wherein said calculating step is repeated at each of said plurality of positions along the source scanning trajectory provided by said scanning step.

8. The method of claim 5, wherein said calculating step is repeated at each of said plurality of positions along the source scanning trajectory provided by said scanning step.

9. The method of claim 1, wherein said source provides x-ray energy.

10. A method in accordance with claim 1, wherein the geometric surface surrounding the field of view is cylindrical.

11. A method in accordance with claim 1, wherein said using step comprises moving the source along the scanning trajectory in a step-wise scan.

12. A method in accordance with claim 1, wherein said using step comprises moving the source along the scanning trajectory in a continuous scan.

13. A method in accordance with claim 1, wherein the step of using the cone beam source comprises maintaining the source and detector in stationary positions while translating and rotating the object so as to effect relative movement of the source with respect to the object.

14. A method in accordance with claim 1, wherein the step of using the cone beam source comprises maintaining the object in a stationary position while moving the source and detector to achieve relative movement of the source with respect to the object.

15. A method in accordance with claim 1, wherein the step of using the cone beam source comprises translating the object while rotating the source and detector to achieve spiral relative movement of the source with respect to the object.

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

a cone beam source for applying cone beam energy to at least a portion of the object;

a two-dimensional area detector positioned to receive cone beam energy from the source, the area detector comprising a plurality of detector elements arranged in an array of rows and columns and having a predetermined number M of rows of detector elements between top and bottom edges thereof that extend along a direction generally parallel to the predetermined axis, plus an extra number N rows of detector elements adjacent both of the top and bottom edges thereof, where N is greater than or equal to 1;

a scanning device operatively coupled to one of the object or the source and area detector for causing relative motion between the source and the object such that the source moves along a scanning trajectory relative to the object;

trajectory defining apparatus operatively coupled to the scanning device to cause the scanning device to provide as said scanning trajectory a spiral path defining a plurality of spaced stages on a predetermined geometric surface surrounding the field of view such that each plane of a plurality of planes passing through the field of view intersects the scanning trajectory in at least one point and intersects the area detector along a line segment L that extends towards top and bottom edges of the area detector, said predetermined number M of rows of detector elements of the area detector spanning at least two consecutive stages of the spiral path having the largest spacing therebetween;

beam energy detecting apparatus for acquiring cone beam projection data in response to said source being at a plurality of positions along the source scanning trajectory, said cone beam projection data corresponding to respective portions of the object;

a processor for calculating Radon derivative data from cone beam projection data in adjacent detector line segments L₁ and L₂ that are parallel to line segment L and translated in an orthogonal direction therefrom by an amount which is less than or equal to N times the spacing between said extra N rows of detector elements contained at the top and bottom edges, respectively, of said area detector, where N is greater than or equal to 1; and

means for reconstructing an image of the object using said Radon derivative data.

17. The system in accordance with claim 16, wherein said M rows of detector elements span at least two stages, but less than four stages defined by said spiral path.

18. The system in accordance with claim 16, wherein N=1.

19. The system in accordance with claim 16, wherein said means for reconstructing an image includes a display for displaying an image of the object based upon said Radon Derivitive data.

20. The system in accordance with claim 16, wherein said means for reconstructing an image includes a processor for performing an inverse Radon transformation of said Radon derivative data.

21. The system in accordance with claim 16, wherein said processor repeats said calculation of Radon derivative data for a plurality of closely spaced line segments L developed as a result of the intersection of said detector with a respective plurality of rotated planes that pass through the field of view and also intersect the scanning trajectory in at least one point.

22. The system in accordance with claim 16, wherein said processor repeats said calculation of Radon derivative data as to develop Radon derivative data for all points in the Radon space of said object in the field of view.

23. The system in accordance with claim 21, wherein said processor repeats said calculation of Radon derivative data as to develop Radon derivative data for all points in the Radon space of said object in the field of view.

24. The system in accordance with claim 16, wherein said cone beam source comprises a source of x-ray energy.

25. The system in accordance with claim 16, wherein said scanning device causes the source to move along the scanning trajectory in a continuous scan.

26. The system in accordance with claim 16, wherein said scanning device causes the source to move along the scanning trajectory in a step-wise scan.

27. The system in accordance with claim 16, wherein said scanning device maintains the source and detector in stationary positions while translating and rotating the object so as to effect relative movement of the source with respect to the object.

28. The system in accordance with claim 16, wherein said scanning device maintains the object in a stationary position while moving the source and detector to achieve relative movement of the source with respect to the object.

29. The system in accordance with claim 16, wherein said scanning device translates the object while rotating the source and detector to achieve spiral relative movement of the source with respect to the object.