Cone beam reconstruction using helical data collection paths

• US 5,404,293 A
• Filed: 10/05/1992
• Issued: 04/04/1995
• Est. Priority Date: 06/11/1991
• Status: Expired due to Term

• Current Assignee: University of Utah (State of Utah)
• Original Assignee: University of Utah (State of Utah)
• Inventors: Weng, Yi, Gullberg, Grant T., Zeng, Gengsheng L.
• Primary Examiner(s): Hayes, Gail O.
• Application Number:
  US07/956,805
• Time in Patent Office: 911 Days
• Field of Search: 364/413.19, 364/413.13, 364/413.15, 364/413.2, 364/413.16, 128/660.01, 128/660.06
• 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/037 Emission tomography
• 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/416 Exact reconstruction
• CPC Class Codes: G06T 2211/421 Filtered back projection [FBP]

Abstract

Radiation passing through a cone beam collimator is received by a radiation detector (10) such as a gamma camera head, as the gamma camera head is moved in a helical orbit. Data g(n,u,v) collected during the helical orbit is scaled (42) to scaled data G(n,u,v). A first partial derivative ∂G (n,u,v)/∂u is taken (46u) with respect to a horizontal direction and a second partial derivative ∂G (n,u,v)/∂v is taken (46v) with respect to a vertical direction. The partial derivatives are linearly combined (48) by being multiplied by sine and cosine values of an angle α between the horizontal direction u and an arbitrary direction p in the detector plane to form partial derivatives ∂G(n,u,v)/∂p. The coordinate system of the derivatives is converted (60) from the (n,u,v) coordinate system to an (n,α,p) coordinate system. The first derivatives are projected (62) i.e. summed row by row, onto a q axis which is perpendicular to the p arbitrary direction. The one-dimensional projection arrays are rebinned (70) to form first derivative Radon domain data R'"'"'(θ,ρ). A second derivative R" (θ,ρ) is taken (74) of the Radon domain data. The second derivative Radon domain data is backprojected (80) into an image memory (82) and displayed on a video monitor (84).

Claims

1. An apparatus for generating an image representation of an interior portion of a subject, the apparatus comprising:

a radiation detection means for receiving radiation travelling along a cone of rays which converge at a focal point and for generating electrical data indicative thereof;

a means for moving the radiation detecting means such that the focal point moves in a helical orbit and the radiation detection means generates helical orbit data;

a partial derivative means for generating a partial derivative of the helical orbit data with respect to an arbitrary direction in a plane of the detector means;

an integrating means for projecting the derivative along a direction perpendicular to the arbitrary direction in the detector plane;

a rebinning means for rebinning the projections into first derivative Radon domain data;

a second derivative means for generating a second derivative of the Radon domain data;

a backprojecting means for reconstructing the second derivative Radon domain data into a three-dimensional image representation;

an image memory means for storing the three-dimensional image representation.

2. The apparatus as set forth in claim 1 wherein the rebinning means includes a means for mapping from preselected points in a cone beam domain to corresponding points in the first derivative Radon domain.

3. The apparatus as set forth in claim 1 wherein the rebinning means includes a means for mapping from preselected points in the first derivative Radon domain to corresponding points in a cone beam domain.

4. The apparatus as set forth in claim 1 including a monitor means for converting a portion of the image representation into a human readable display.

5. The apparatus as set forth in claim 1 further including a scaling means for scaling the helical data, the scaling means being disposed between radiation detection means and the partial derivative means.

6. The apparatus as set forth in claim 1 wherein the partial derivative means include:

a first partial derivative means for taking a partial derivative along a horizontal direction;

a second partial derivative means for taking a partial derivative along a vertical direction; and

a linear combining means for linearly combining the partial derivatives along the horizontal and vertical directions to produce the partial derivative along the arbitrary direction.

7. The apparatus as set forth in claim 6 wherein the linear combining means includes a means for multiplying the partial derivatives along the vertical and horizontal directions with sine and cosine values of an angle between one of the horizontal and vertical directions and the arbitrary direction in the detector plane.

8. The apparatus as set forth in claim 7 wherein for each angle, the integrating means sums row by row to produce a one-dimensional array.

9. A method for generating an image representation of an interior portion of a subject, the method comprising:

generating cone beam electrical data from detected radiation travelling along a cone of rays which converge at a focal point, the radiation being detected by a radiation detection means along a detection plane;

moving the radiation detecting means such that the focal point moves in a helical orbit;

taking a partial derivative of the electrical data with respect to an arbitrary direction along the detection plane;

projecting the partial derivative along a direction perpendicular to the arbitrary direction along the detection plane;

rebinning the projections into first derivative Radon domain data;

taking a second derivative of the Radon domain data;

backprojecting the second derivative Radon domain data into a three-dimensional image representation;

converting at least a portion of the image representation into a human readable display.

10. The method as set forth in claim 9 further including storing the three-dimensional image representation.

11. The method as set forth in claim 9 further including scaling the electrical data prior to taking the partial derivative along the arbitrary direction in the detection plane.

12. The method as set forth in claim 9 wherein the rebinning includes mapping from preselected points of the projected partial derivative electrical data along the detection plane to corresponding points in the first derivative Radon domain.

13. The method as set forth in claim 9 wherein the rebinning includes mapping from preselected points in the first derivative Radon domain to corresponding points on the detection plane.

14. The method as set forth in claim 9 wherein taking the partial derivative includes:

taking a partial derivative along a horizontal direction;

taking a partial derivative along a vertical direction;

linearly combining the partial derivatives along the horizontal and vertical directions to produce the partial derivative along the arbitrary direction.

15. The method as set forth in claim 14 wherein the linear combining includes multiplying the partial derivatives along the vertical and horizontal directions with sine and cosine values of an angle between one of the horizontal and vertical directions and the arbitrary direction in the detection plane.

16. The method as set forth in claim 15 wherein for each angle, the projecting step includes summing row by row to produce one-dimensional arrays which are rebinned.

17. A method for generating an image representation of an interior portion of a subject, the method comprising:

generating cone beam electrical data from detected radiation travelling along a cone of rays which converge at a focal point, the radiation being detected along a detection plane such that the focal point moves along a helical orbit;

taking partial derivatives of the electrical data with respect to an arbitrary direction along the detection plane;

projecting the partial derivative along a direction perpendicular to the arbitrary direction along the detection plane;

determining points in a Radon domain which correspond to preselected points in the detection plane;

mapping the projected partial derivatives into first derivative Radion domain data;

generating second derivative Radon domain data from the first derivative Radon domain data;

backprojecting the second derivative Radon domain data into a three-dimensional image representation;

converting at least a portion of the image representation into a human readable display.