Helical interpolation for an asymmetric multi-slice scanner
Gespeichert in:
| Titel: | Helical interpolation for an asymmetric multi-slice scanner |
|---|---|
| Patent Number: | 7,062,009 |
| Publikationsdatum: | June 13, 2006 |
| Appl. No: | 10/659152 |
| Application Filed: | September 10, 2003 |
| Abstract: | A helical interpolation filter is provided for a multi-slice scanner having an asymmetric detector array. For projection data generated in the symmetric region of the scanner, in which complementary data are available, complementary interpolation is performed. For projection data generated in the asymmetric region of the scanner, in which no complementary data are available, direct interpolation is performed. A blending region is defined near each boundary between the symmetric and the asymmetric regions. In the blending region, a combination of direct and complementary interpolation is performed. |
| Inventors: | Karimi, Seemeen S. (Brookline, MA, US); Crawford, Carl R. (Brookline, MA, US) |
| Assignees: | Analogic Corporation (Peabody, MA, US) |
| Claim: | 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. |
| Claim: | 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. |
| Claim: | 3. A CT system in accordance with claim 1 , wherein the number of said rows of detector elements is about four. |
| Claim: | 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. |
| Claim: | 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: [mathematical expression included] 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; w i c (β, γ)(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 [mathematical expression included] |
| Claim: | 6. A CT system in accordance with claim 5 , wherein the number of rows of detector elements is four, and wherein said direct interpolation weights are defined by the relationship: [mathematical expression included] |
| Claim: | 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. |
| Claim: | 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. |
| Claim: | 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), w i (β,γ)=(1−α f (x)) w i c (β,γ)+α f (x) w i d (β,γ); within the range (−γ s ,−γ b), w i (β,γ)=(1−α f (x)) w i c (β,γ); where α f (x)=3x 3 −2x 2 and 1≦i≦4, and [mathematical expression included] |
| Claim: | 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. |
| Claim: | 11. A CT system in accordance with claim 1 , wherein said x-ray beam comprises a plurality of fan-shaped beams. |
| Claim: | 12. A CT system in accordance with claim 1 , wherein said translation of said object along said axis occurs at a constant speed. |
| Claim: | 13. A CT system in accordance with claim 1 , wherein the angular range of the weighted projection data is greater than 2π. |
| Claim: | 14. A CT system in accordance with claim 1 , wherein said helical interpolator includes an interpolation weight generator for generating said helical interpolation weights. |
| Claim: | 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. |
| Claim: | 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. |
| Claim: | 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. |
| Claim: | 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. |
| Claim: | 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. |
| Claim: | 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: [mathematical expression included] 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; w i c (β, γ)(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 [mathematical expression included] |
| Claim: | 21. A method in accordance with claim 20 , wherein the number of rows of detector elements is four, and wherein said direct interpolation weights are defined by the relationship: [mathematical expression included] |
| Claim: | 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. |
| Claim: | 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. |
| Claim: | 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), w i (β,γ)=(1−α f (x)) w i c (β,γ)+α f (x) w i d (β,γ); within the range (−γ s ,−γ b), w i (β,γ)=(1−α f (x)) w i c (β,γ); where α f (x)=3x 3 −2x 2 and 1≦i≦4, and [mathematical expression included] |
| Current U.S. Class: | 378/19 |
| Patent References Cited: | 6108575 August 2000 Besson 6351514 February 2002 Besson 6415012 July 2002 Taguchi et al. 6600802 July 2003 Hsieh |
| Assistant Examiner: | Thomas, Courtney |
| Primary Examiner: | Bruce, David V. |
| Attorney, Agent or Firm: | McDermott Will & Emery |
| Dokumentencode: | edspgr.07062009 |
| Datenbank: | USPTO Patent Grants |
Schreiben Sie den ersten Kommentar!