Claims
- 1. A CT system for reconstructing at least one image of an object, said system comprising:
a. an x-ray source for generating x-rays, said x-ray source being mounted on a gantry for rotation about an axis; b. an x-ray detector system mounted opposite said x-ray source for providing a set of projection data with respect to said object as the object is translated along said axis, said x-ray detector system including a plurality of rows of detector elements; wherein for each detector element, a corresponding detector ray is defined by x-ray photons traveling from said x-ray source to said detector element; and wherein said detector rays define an x-ray beam that is asymmetric in a plane perpendicular to said axis, said x-ray beam including a symmetric region in which one or more complementary rays can be found for each detector ray, and an asymmetric region in which no complementary ray can be found for any detector ray; c. an interpolator for interpolating said projection data from each detector element onto a slice plane by multiplying said data with helical interpolation weights;
wherein for projection data resulting from detector rays lying within the symmetric region of said x-ray beam, the helical interpolation weights are complementary interpolation weights that weigh the data from complementary rays in proportion to the distance from each detector element row to said slice plane; and wherein for projection data resulting from detector rays that lie within the asymmetric region of said x-ray beam, the helical interpolation weights are direct interpolation weights that weigh the data at a given azimuthal angle from different rows in proportion to the distance from each detector element row to said slice plane; and d. an image reconstructor for reconstructing a tomographic image of said object using said helically interpolated projection data.
- 2. A CT system in accordance with claim 1, wherein said plurality of rows of detector elements are substantially parallel, and wherein said substantially parallel rows are disposed side-by-side along said axis, so that during a given sampling period, projection data can be acquired that are representative of a plurality of sections of said object.
- 3. A CT system in accordance with claim 1, wherein the number of said rows of detector elements is about four.
- 4. A CT system in accordance with claim 1, wherein for all detector rays lying within said symmetric region of said x-ray beam, the data for each detector ray and the data for its complementary rays used to interpolate the x-ray beam come from different ones of said plurality of rows of detector elements.
- 5. A CT system in accordance with claim 1, wherein the number of rows of detector elements is four;
wherein said rotation of said gantry defines a gantry rotation plane; and wherein said complementary interpolation weights are defined by the relationship: 8w1c(β,γ)=α(x1){β-β2c(-)β1-β2c(-),β2c(-)<β≤β1β3c(-)-ββ3c(-)-β1,β1<β<β3c(-) w2c(β,γ)={β-β3c(-)β2-β3c(-),β3c(-)<β≤β2β3c(+)-ββ1c(+)-β2,β2<β<β1c(+) w3c(β,γ)={β-β1c(+)β3-β1c(+),β1c(+)<β≤β3β2c(+)-ββ2c(+)-β3,β3<β<β2c(-) w4c(β,γ)=(1-α(x4)){β-β2c(+)β4-β2c(+),β2c(+)<β≤β4β3c(+)-ββ3c(+)-β4,β4<β<β3c(+) where β denotes the fan angle, said fan angle being defined as the angle between a line from said x-ray source to the isocenter of said CT system, and a fixed coordinate axis in said gantry rotation plane; βi(i=1, . . . 4) denotes the fan angle for the i-th row indicative of the view at which the row i is directly under the point of intersection of said slice plane with said axis of rotation; γ denotes the angle between the projection of said detector ray on said gantry rotation plane, and a center ray disposed in the center of said x-ray beam; wci(β, γ)(i=1, . . . 4) denotes the complementary interpolation weight for a detector element in the i-th row and ray angle γ with respect to said center ray, at a fan angle β; βic(i=1, . . . 4) denotes the fan angle of a ray complementary to the ray at βi, γ, (+) indicates that said complementary ray comes from a current rotation; (−) indicates that said complementary ray comes from a previous rotation; and 9α(x)=3x2-2x3,x1=β-β2c(1)β3c(-)-β2c(-) and x4=β-β2c+β3c(+)-β2c(+).
- 6. A CT system in accordance with claim 1, wherein the number of rows of detector elements is four, and
wherein said direct interpolation weights are defined by the relationship: 10w1d(β,γ)={α(x1)(1-β1-ββ2-β1), β2c(-)<β≤β1α(x1)β2-ββ2-β1, β<β≤B3c(-)β2-ββ2-β1, β3c(-)<β≤β2w2d(β,γ)={β-β1β2-β1,β1<β≤β2β3-ββ3-β2,β2<β<β3w3d(β,γ)={β-β2β3-β2,β2<β≤β3β4-ββ4-β3,β3<β<β4w4d(β,γ)={β-β3β4-β3, β2c(-)<β≤β1(1-α(x4))β-β3β4-β3, β2c(+)<β≤β4(1-α(x4))(1-β-β4β4-β3,) β4<β≤β3c(+)and where β, βi(i=1, . . . 4), γ, wci(β, γ)(i=1, . . . 0.4),
βic(i=1, . . . 4), (+), (−), α(x), x1, and X4 are defined as in claim 5.
- 7. A CT system in accordance with claim 1, wherein said x-ray beam further comprises at least one blending region within said symmetric region adjacent to the boundary between said symmetric region and said asymmetric region, and
wherein for projection data resulting from detector rays lying within said blending region, the helical interpolation weights are blending interpolation weights created by a combination of direction and complementary interpolation weights.
- 8. A CT system in accordance with claim 7, wherein said blending region occurs at γ angle ranges between (γb, γs), and (−γb, −γs), where γb represents the value of γ at the start of said blending region within said symmetric region of said x-ray beam, and γs represents the value of γ at said boundary between said symmetric region and said asymmetric region.
- 9. A CT system in accordance with claim 7, wherein said blending interpolation weights are defined by the following relationship:
within the range (γb, γs), wi(β, γ)=(1−αf(x))wci(β, γ)+αf(x)wdi(β, γ); within the range (−γs, −γb), wi(β, γ)=(1−αf(x))wci(β, γ); where αf(x)=3x3−2x2 and 1≦i≦4, and 11x=γ-γbγs-γb.
- 10. A CT system in accordance with claim 1, wherein said image reconstructor includes means for performing 2D (two-dimensional) backprojection of said helically interpolated projection data.
- 11. A CT system in accordance with claim 1, wherein said x-ray beam comprises a plurality of fan-shaped beams.
- 12. A CT system in accordance with claim 1, wherein said translation of said object along said axis occurs at a constant speed.
- 13. A CT system in accordance with claim 1, wherein the angular range of the weighted projection data is greater than 2π.
- 14. A CT system in accordance with claim 1, wherein said helical interpolator includes an interpolation weight generator for generating said helical interpolation weights.
- 15. A CT system in accordance with claim 1, wherein for all detector rays lying within said symmetric region of said x-ray beam, the data for each detector ray and the data for its complementary rays used to interpolate the x-ray beam come from every one of said plurality of rows of detector elements.
- 16. A method of reconstructing at least one image of an object, the method comprising:
a. helically scanning said object with x-rays to acquire tomographic projection data representative of said object while said object is translated along an axis, said x-rays being generated by an x-ray source mounted on a gantry for rotation about said axis, said x-rays being incident upon a multi-slice x-ray detector system having a plurality of substantially parallel rows of detector elements;
wherein for each detector element, a corresponding detector ray is defined by x-ray photons traveling from said x-ray source to said detector element; and wherein said detector rays define an x-ray beam that is asymmetric in a plane perpendicular to said axis, said x-ray beam including a symmetric region in which one or more complementary rays can be found for each detector ray, and an asymmetric region in which no complementary ray can be found for any detector ray; b. helically interpolating the projection data of the detector elements by multiplying the data with helical interpolation weights;
wherein for projection data resulting from detector rays lying within the symmetric region of said x-ray beam, the helical interpolation weights are complementary interpolation weights that weigh the data from complementary rays in proportion to the distance from each detector element row to said slice plane; and wherein for projection data resulting from detector rays that lie within the asymmetric region of said x-ray beam, the helical interpolation weights are direct interpolation weights that weigh the data at a given azimuthal angle from different rows in proportion to the distance from each detector element row to said slice plane; and C. reconstructing a tomographic image of said object, using said helically interpolated projection data.
- 17. A method in accordance with claim 16, wherein said plurality of rows of detector elements are substantially parallel, and wherein said substantially parallel rows are disposed side-by-side along said axis, so that during a given sampling period, projection data can be acquired that are representative of a plurality of sections of said object.
- 18. A method in accordance with claim 16, wherein for all detector rays lying within said symmetric region, the data for each detector ray and the data for its complementary ray used to interpolate the x-ray beam come from different ones of said plurality of rows of detector elements.
- 19. A method in accordance with claim 16, wherein for all detector rays lying within said symmetric region of said x-ray beam, the data for each detector ray and the data for its complementary rays used to interpolate the x-ray beam come from every one of said plurality of rows of detector elements.
- 20. A method in accordance with claim 16, wherein the number of rows of detector elements is four;
wherein said rotation of said gantry defines a gantry rotation plane; and wherein said complementary interpolation weights are defined by the relationship: 12w1c(β,γ)=α(x1){β-β2c(-)β1-β2c(-),β2c(-)<β≤β1β3c(-)-ββ3c(-)-β1,β1<β<β3c(-)w2c(β,γ)={β-β3c(-)β2-β3c(-),β3c(-)<β≤β2β3c(+)-ββ1c(+)-β2,β2<β<β1c(+)w3c(β,γ)={β-β1c(+)β3-β1c(+),β1c(+)<β≤β3β2c(+)-ββ2c(+)-β3,β3<β<β2c(-)w4c(β,γ)=(1-α(x4)){β-β2c(+)β4-β2c(+),β2c(+)<β≤β4β3c(+)-ββ3c(+)-β4,β4<β<β3c(+)where β denotes the fan angle, said fan angle being defined as the angle between a line from said x-ray source to the isocenter of said CT system, and a fixed coordinate axis in said gantry rotation plane; βi(i=1, . . . 4) denotes the fan angle for the i-th row indicative of the view at which the row i is directly under the point of intersection of said slice plane with said axis of rotation; γ denotes the angle between the projection of a detector ray onto said gantry rotation plane and a center ray disposed in the center of said x-ray beam; wci(β, γ)(i=1, . . . 4) denotes the complementary interpolation weight for a detector element in the i-th row and ray angle γ with respect to said center ray, at a fan angle β; βic(i=1, . . . 4) denotes the fan angle of a ray complementary to the ray at βi, γ, (+) indicates that said complementary ray comes from a current rotation; (−) indicates that said complementary ray comes from a previous rotation; and 13α(x)=3x2-2x3,x1=β-β2c(1)β3c(-)-β2c(-) and x4=β-β2c+β3c(+)-β2c(+).
- 21. A method in accordance with claim 16, wherein the number of rows of detector elements is four, and
wherein said direct interpolation weights are defined by the relationship: 14w1d(β,γ)={α(x1)(1-β1-ββ2-β1),β2c(-)<β≤β1α(x1)β2-ββ2-β1, β<β≤B3c(-)β2-ββ2-β1, β3c(-)<β≤β2w2d(β,γ)={β-β1β2-β1,β1<β≤β2β3-ββ3-β2,β2<β<β3w2d(β,γ)={β-β2β3-β2,β2<β≤β3β4-ββ4-β3,β3<β<β4w4d(β,γ)={β-β3β4-β3, β3<β≤β2c(+)(1-α(x4))β-β3β4-β3, β2c(+)<β≤β4(1-α(x4))(1-β-β4β4-β3,) β4<β≤β3c(+)and where β, βi(i=1, . . . 4), γ, wci(β, γ)(i=1, . . . 4),
βic(i=1, . . . 4), (+), (−), α(x), x1, and x4 are defined as in claim 20.
- 22. A method in accordance with claim 16, wherein the step of reconstructing said image of said object includes the step of performing a 2D (two-dimensional) backprojection of said helically interpolated projection data.
- 23. A method in accordance with claim 16, wherein said x-ray beam further comprises at least one blending region that is contained within said symmetric region and that is proximate to the boundary between said symmetric region and said asymmetric region, and
wherein for projection data resulting from detector rays lying within said blending region, the helical interpolation weights are blending interpolation weights created by a combination of direct and complementary interpolation weights.
- 24. A method in accordance with claim 23, wherein said blending region occurs at γ angle ranges between (γb, γs), and (−γb, −γs), where γb represents the value of γ at the start of said blending region within said symmetric region of said x-ray beam, and γs represents the value of γ at said boundary between said symmetric region and said asymmetric region; and
wherein said blending interpolation weights are defined by the following relationship: within the range (γb, γs), wi(β, γ)=(1−αf(x))wci(β, γ)+αf(x)wdi(β, γ); within the range (−γs, −γb), wi(β, γ)=(1−αf(x))wci(β, γ); where αf(x)=3x3−2x2 and 1≦i≦4, and 15x=γ-γbγs-γb.
CROSS REFERENCES TO RELATED APPLICATIONS
[0001] This application claims the benefit of priority under 35 U.S.C. §119(e) from co-pending, commonly owned U.S. provisional patent application serial No. 60/410,244 (entitled “Helical Interpolation For An Asymmetric Multi-Slice Scanner” and having an attorney docket number ANAK-240PR), which was filed on Sep. 12, 2002.
Provisional Applications (1)
|
Number |
Date |
Country |
|
60410244 |
Sep 2002 |
US |