3D image reconstruction for helical partial cone beam scanners using wedge beam transform

• US 6,104,775 A
• Filed: 10/29/1998
• Issued: 08/15/2000
• Est. Priority Date: 10/29/1998
• Status: Expired due to Term

• Current Assignee: Picker International, Inc. (Koninklijke Philips N.V.)
• Original Assignee: Picker International, Inc. (Koninklijke Philips N.V.)
• Inventors: Tuy, Heang K.
• Primary Examiner(s): BRUCE, DAVID VERNON
• Application Number:
  US09/181,980
• Time in Patent Office: 656 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: 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

A method of image reconstruction from wedge beam data includes collecting cone beam data in two-dimensional arrays. The collected data corresponds to rays of penetrating radiation diverging in two dimensions from a common vertex which travels along a curve. Each element of the data is related to line integrals of an object being reconstructed taken along each ray. A local coordinate system is defined having three mutually orthogonal axes and an origin at the vertex. The third axis of the local coordinate system extends in a direction tangential to the curve at the vertex. The collected data is rebinned into a wedge beam format wherein sets of parallel rays are grouped to define planes of radiation that angularly diverge from a common axis. A one-dimensional convolution of the rebinned data is computed in the local coordinate system along a direction parallel to the third axis. Finally, the convolved data is weighted and three-dimensionally backprojected. In a preferred embodiment, the curve defines a helical path, and the kernel for the one-dimensional convolution of the rebinned data is a ramp convolver.

Claims

1. A method of image reconstruction from wedge beam data comprising:

(a) collecting cone beam data in two-dimensional arrays, said collected data corresponding to rays of penetrating radiation diverging in two dimensions from a common vertex as the vertex travels along a curve with each element of the data related to line integrals of an object being reconstructed taken along each ray;

(b) defining a local coordinate system having three mutually orthogonal axes and an origin at the vertex, wherein a third axis thereof extends in a direction tangential to the curve at the vertex;

(c) rebinning the collected data into a wedge beam format, wherein sets of parallel rays are grouped to define planes of radiation that angularly diverge from a common axis;

(d) computing a one-dimensional convolution of the rebinned data in the local coordinate system along a direction parallel to the third axis;

(e) weighting the convolved data; and

,

(f) three-dimensionally backprojecting the weighted convolution.

2. The method according to claim 1, wherein the common axis from which the defined planes angularly diverge is coincident with the third axis of the local coordinate system.

3. The method according to claim 2, wherein for each position of the vertex as it travels along the curve, the data collected is partial cone beam data which is not truncated in a direction parallel to the tangent of the curve at the position of the vertex and is truncated in at least one other orthogonal direction.

4. The method according to claim 3, wherein the curve defines a helical path relative to the object.

5. The method according to claim 4, wherein a kernel for the one-dimensional convolution of the rebinned data is a ramp convolver.

6. The method according to claim 5, wherein the weighting of the convolved data is proportional to an average of values of a multiplicity function, said values corresponding to planes containing a line extending from the vertex to a projection of a point of reconstruction on a plane defined by first and second axes of the local coordinate system.

7. The method according to claim 6, wherein values of the multiplicity function are determined based on a number of intersections their corresponding plane has with a portion of the helical path traveled by the vertex to collect enough data for reconstruction.

8. The method according to claim 7, wherein the values of the multiplicity function are equal to the inverse of the number of intersections.

9. The method according to claim 7, wherein the portion of the helical path traveled by the vertex to collect enough data for rebinning into the wedge beam format extends from a point where the vertex is -π

/2 revolutions away from an axial plane containing the point of reconstruction to a point where the vertex is +π

/2 revolutions away from the axial plane containing the point of reconstruction.

10. The method according to claim 1, wherein full cone beam data is collected in step (a).

11. A CT scanner comprising:

a first gantry which defines an examination region;

a rotating gantry mounted on the first gantry for rotation about the examination region;

a source of penetrating radiation arranged on the rotating gantry for rotation therewith, said source of penetrating radiation emitting a cone-shaped beam of radiation that passes through the examination region as the rotating gantry rotates;

a subject support which holds a subject being examined at least partially within the examination region, wherein at least one of the first gantry and the subject support is translated such that the subject passes through the examination region while the rotating gantry is rotated and the source of penetrating radiation follows a helical path relative to the subject;

a two-dimensional array of radiation detectors arranged to receive the radiation emitted from the source of penetrating radiation after it has traversed the examination region;

a reconstruction processor which reconstructs images of the subject from data collected by the two-dimensional array of radiation detectors, said reconstruction processor including;

a rebinning processor which rebins the data into a wedge beam format such that sets of parallel rays of radiation are grouped to define planes of radiation that angularly diverge from a common axis which is tangential to the helical path;

a convolver which takes the data from the rebinning processor and performs a one-dimensional convolution on it in a local coordinate system along a direction which is tangential to the helical path;

a data processor which takes the data from the convolver and weights it; and

,

a backprojector which takes the data from the data processor and three-dimensionally backprojects it into an image memory; and

,

a human viewable display which accesses the image memory to display reconstructed images of the subject.

12. The CT scanner according to claim 11, wherein for each position of the source of penetrating radiation along the helical path, the data collected is not truncated in a direction parallel to the tangent of the helical path at the position of the source of penetrating radiation and is truncated in at least one other orthogonal direction.

13. The CT scanner according to claim 12, wherein the convolver employs a ramp function as a kernel for the one-dimensional convolution.

14. The CT scanner according to claim 13, wherein the data processor multiplies the data by a weight that is proportional to an average of values of a multiplicity function, said values corresponding to planes containing a ray extending from the position of the source of penetrating radiation to a projection of a point of reconstruction on a plane which contains the source of penetrating radiation and is normal to the common axis.

15. The CT scanner according to claim 14, wherein each value of the multiplicity function is determined based on a number of intersections its corresponding plane has with an amount of the helical path followed by the source of penetrating radiation, said amount allowing enough data collection for reconstruction.

16. The CT scanner according to claim 15, wherein the values of the multiplicity function are equal to one divided by the number of intersections.

17. The CT scanner according to claim 15, wherein the amount of the helical path followed by the source of penetrating radiation which allows enough data collection for rebinning into the wedge beam format extends from a point where the source of penetrating radiation is -π

/2 revolutions away from an axial plane containing the point of reconstruction to a point where the source of penetrating radiation is +π

/2 revolutions away from the axial plane containing the point of reconstruction.