Method and apparatus for X-ray computed tomography

• US 5,640,436 A
• Filed: 01/19/1996
• Issued: 06/17/1997
• Est. Priority Date: 01/26/1995
• Status: Expired due to Term

• Current Assignee: Hitachi Medical Corporation (Hitachi, Ltd.)
• Original Assignee: Hitachi Medical Corporation (Hitachi, Ltd.)
• Inventors: Kawai, Hiroyuki, Sekihara, Kensuke
• Primary Examiner(s): Porta, David P.
• Assistant Examiner(s): Bruce, David Vernon
• Application Number:
  US08/588,918
• Time in Patent Office: 515 Days
• Field of Search: 378/4, 378/901, 364/413, 364/18
• US Class Current: 378/4
• CPC Class Codes: A61B 6/5258 involving detection or reduction of artifacts or noise
• 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 for X-ray computed tomography includes a process for rotating around an object an X-ray source for applying a cone-beam X ray and a scanner having a two-dimensional X-ray detector mounted thereon for deriving projection data and a process for preparing a geometric distortion correction table for correcting an image geometric distortion of the two-dimensional X-ray detector for the projection data, for reconstructing a distribution of X-ray attenuation coefficients of the object from said corrected projection data. The method further includes a process for calculating a function for evaluating the correction table, a process for comparing the calculated value with a threshold value, a process for correcting the correction table through the interpolation based on the compared result, a process for correcting an image geometric distortion of the projection data by using the corrected table, a process for detecting a location of a view field boundary of the detector from the distortion-corrected projection data, and a process for generating the projection data of a predetermined width from the detected location of the view field boundary through the effect of the extrapolation. The method operates to reconstruct a distribution of X-ray attenuation coefficients of the object from the extrapolated projection data.

Claims

1. A method for cone-beam X-ray computed tomography in which a scanner mounting an X-ray source for applying a cone-beam X ray and a two-dimensional X-ray detector for detecting an X-ray is rotated around an object to measure omnidirectional projection data around said object and a distribution of X-ray attenuation coefficients of said object is reconstructed, comprising the steps of:

correcting an image geometric distortion of said projection data;

detecting a view field boundary of said two-dimensional X-ray detector from said projection data whose geometric distortion is corrected;

extrapolating projection data of a predetermined width from said view field boundary detected in the step of detecting a view field and boundary; and

reconstructing said distribution of X-ray attenuation coefficients of said object from said projection data extrapolated in the step of extrapolating projection data.

2. A method for cone-beam X-ray computed tomography as claimed in claim 1, wherein the step of correcting comprises a step of setting, a value indicating being out of view field to a region whose geometric distortion is corrected in said projection data before the image geometric distortion is corrected, said region located out of the view field of said two-dimensional X-ray detector.

3. A method for cone-beam X-ray computed tomography as claimed in claim 2, wherein said value to be set is zero.

4. A method for cone-beam X-ray computed tomography as claimed in claim 2, wherein the step of detecting the view field boundary comprises a step of sequentially scanning said distortion-corrected projection data in parallel to the tangent of the rotation of said scanner, and if a point has any value other than the value indicating out of the view field, detecting the point as said boundary of the view field.

5. A method for cone-beam X-ray computed tomography as claimed in any one of claims 1 to 4, wherein the step of extrapolating comprises a step of, in said projection data whose boundary of the view field is detected, extrapolating the projection data of only said predetermined width from said detected view field boundary to said region of out of the view field in parallel to the rotation of said scanner.

6. A method for cone-beam X-ray computed tomography as claimed in claim 5, wherein the step of extrapolating comprises a step of extrapolating so that said projection data in the region located out of the view field is concatenated with said projection data located within said view field and forms such a part of an arc of an ellipse as taking zero at a point spaced by a certain width from said view field boundary.

7. A method for cone-beam X-ray computed tomography as claimed in claim 5, wherein the step of extrapolating comprises a step of, if the rotating direction of said scanner is set to a u axis and the center axis of rotation of said scanner is set to a v axis on a projection plane of said two-dimensional X-ray detector, assuming that a projection value at the k-th projection location (u, v) is P(u, v, k), the location of said view field boundary detected about a certain value of v is (a(v), v) and (b(v), v) (wherein a(v)<

b(v)), the projection of the center axis of rotation is a straight line u=c, and the certain width to be extrapolated is D, in the region located out of the view field of u<

a(v) and b(v)<

u, extrapolating the projection data located out of the view field according to the expression of;

##EQU12##

8. An apparatus for cone-beam X-ray computed tomography having a scanner mounting an X-ray source for applying a cone-beam X ray and a two-dimensional X-ray detector for detecting an X-ray, means for rotating said scanner around an object, and means for measuring omnidirectional projection data around said object, said apparatus serving to reconstruct a distribution of X-ray attenuation coefficients of said object based on said projection data, comprising:

correcting means for correcting an image geometric distortion of said projection data;

detecting means for detecting a location of a view field boundary of said detector of said projection data whose geometric distortion is corrected by said correcting means;

extrapolating means for generating said projection data of a predetermined width from said view field boundary location detected by said detecting means; and

means for reconstructing the distribution of X-ray attenuation coefficients of said object from said projection data extrapolated in said extrapolating means.

9. An apparatus for cone-beam X-ray computed tomography as claimed in claim 8, wherein said means for correcting an image geometric distortion of said projection data comprises means for setting a value indicating the out of view field to the distortion-corrected region corresponding to the region located out of the view field of said two-dimensional X-ray detector.

10. An apparatus for cone-beam X-ray computed tomography as claimed in claim 9, wherein said value to be set is zero.

11. An apparatus for cone-beam X-ray computed tomography as claimed in claim 9, wherein said means for detecting a location of said view field boundary comprises means for sequentially scanning said projection data in parallel to the rotating direction of said scanner and detecting as the view field boundary points where any value other than the value indicating the out of view field appears.

12. An apparatus for cone-beam X-ray computed tomography as claimed in any one of claims 9 to 11, wherein said extrapolating means comprises means for extrapolating said projection data of said predetermined of said detected view field boundary to said region located out of the view field in parallel to the rotating direction of said scanner in said projection data whose view field boundary is detected.

13. An apparatus for cone-beam X-ray computed tomography as claimed in claim 12, wherein said extrapolating means comprises means for extrapolating the projection data so that said projection data in the region located out of the view field is concatenated with said projection data located within said view field and forms such a part of an arc of an ellipse as assuming zero at a point spaced by a certain width from said view field boundary.

14. An apparatus for cone-beam X-ray computed tomography as claimed in claim 12, wherein if the rotating direction of said scanner is set to a u axis and the center axis of rotation of said scanner is set to a v axis on a projection plane of said two-dimensional X-ray detector, assuming that a projection value at the k-th projection location (u, v) is P(u, v, k), the location of said view field boundary detected about a certain value of v is (a(v), v) and (b(v), v) (a(v)<

b(v)), the projection of the center axis of rotation is a straight line u=c, and the certain width to be extrapolated is D, in the region located out of the view field of u<

a(v) and b(v)<

u, said extrapolating means comprises means for extrapolating said projection data located out of the view field according to the expression of;

##EQU13##

15. A method for cone-beam X-ray computed tomography in which a scanner mounting an X-ray source for applying a cone-beam X ray and a two-dimensional X-ray detector for detecting an X-ray is rotated around an object to measure omnidirectional projection data around said object, an image geometric distortion of said projection data is corrected by referring to a geometric distortion correction table, and a distribution of X-ray attenuation coefficients of said object is reconstructed by using said distortion-corrected projection data comprising the steps of:

preparing said correction table and storing said table in means for saving a correction table;

calculating an evaluation function for evaluating continuity of said prepared correction table between the previous and the subsequent projections of each subject projection on which said projection data is obtained;

deleting said prepared correction table that is not continuously changed between said previous and subsequent projections from said means for saving a correction table, replacing said correction table that is not continuously changed with a new correction table interpolated from said prepared correction tables of the previous and the subsequent projections, and saving said new correction table in said means for saving a correction table; and

reconstructing the distribution of X-ray attenuations of said object.

16. A method for cone-beam X-ray computed tomography as claimed in claim 15, wherein if on the plane of said two-dimensional X-ray detector for detecting said projection data a u axis is set to the rotating direction of said scanner, a v axis is set in parallel to the center axis of rotation of said scanner, and said prepared correction table of said k-directional projection is represented by U(u, v, k) and V(u, v, k) (k=1, 2, . . . , N_k, u=1, 2, . . . , N_i, v=1, 2, . . . , N), said evaluation function is operated to calculate a correlation between said k-directional prepared correction tables in the step of calculating the evaluation function.

17. A method for cone-beam X-ray computed tomography as claimed in claim 15, wherein the step of calculating an evaluation function comprises a step of obtaining if k is any natural number, valU(u, v, k) and valV(u, v, k) (where if k+l is overflown out of the range of k=1, 2, . . . , N_k during the calculation, the tail is circulated to the head of k) as be said evaluation function according to the expression of;

##EQU14##

18. A method for cone-beam X-ray computed tomography as claimed in claim 15, wherein if on the plane of said two-dimensional X-ray detector for detecting said projection data a u axis is set to the rotating direction of said scanner, a v axis is set in parallel to the center axis of rotation of said scanner, and said prepared correction table of said k-directional projection is represented by U(u, v, k) and V(u, v, k) (k=1, 2, . . . , N_k, u=1, 2, . . . , N_i, v=1, 2, . . . , N), and said evaluation function is operated to calculate a correlation between said k-directional prepared correction tables in the step of calculating said evaluation function, and wherein said prepared correction table in which the values of said evaluation function has a value meeting the specific conditions is determined to be the correction table that is not continuously changed.

19. A method for cone-beam X-ray computed tomography as claimed in claim 15, wherein the step of calculating an evacuation function comprises a step of obtaining if k is any natural number, valU(u, v, k) and valV(u, v, k) (where if k+l is overflown out of the range of k=1, 2, . . . , N_k during the calculation, the tail is circulated to the head of k) as be said evaluation function according to the expression of:

##EQU15## and wherein said prepared correction table in which said evaluation function has a greater value than a specific threshold value is determined to be said correction table that is not continuously changed.

20. A method for cone-beam X-ray computed tomography as claimed in claim 15, wherein the step of calculating an evacuation function comprises a step of obtaining if k is any natural number, valU(u, v, k) and valV(u, v, k) (where if k+l is overflown out of the range of k=1, 2, . . . , N_k during the calculation, the tail is circulated to the head of k) as be said evaluation function according to the expression of:

##EQU16## and wherein if said evaluation function has a greater value than a specific threshold value of said projections in the serial direction (k=k₀, k₀ +1, . . . , k₀ +n), said prepared correction table of said projection in the direction of k=k₀ +1, . . . , k₀ +n-1 is deleted from said means for saving correction table, the correction table of said projection in the k=k₀ +1, . . . , k=k₀ +n-1 is derived by doing the interpolation of said prepared correction table in said k₀ and (k₀ ±

n) directions and then is saved in said means for saving correction table.

21. A method for cone-beam X-ray computed tomography as claimed in claim 20, wherein said interpolation is executed for deriving said prepared correction table U(u, v, k) and V(u, v, k) (where k=k₀ +1, . . . , k₀ +n-1, k₁ <

k₀ +1, k₀ +n-1<

k₂) deleted from said means for saving correction table from said prepared correction table U(u, v, k₁) and V(u, v, k₁) immediately before said deleted correction table in said projecting direction and said prepared correction table U(u, v, k₂) and V(u, v, k₂) immediately after said deleted correction table in said projecting direction according to the expression of;

##EQU17##

22. An apparatus for cone-beam X-ray computed tomography having a scanner mounting an x-ray source for applying a cone-beam X ray and a two-dimensional X-ray detector for detecting an X-ray, means for rotating said scanner around an object, means for measuring omnidirectional projection data around said object, means for preparing a correction table on which an image geometric distortion of said projection data is corrected, and means for correcting an image geometric distortion of said projection data by referring to said correction table, for reconstructing a distribution of X-ray attenuation coefficients of said object on said distortion-corrected projection data, comprising:

means for preparing said correction table having;

means for preparing said correction table of each projection for obtaining said projection data in advance;

evaluating means for calculating an evaluation function for evaluating continuity of said correction tables between the previous and the subsequent projections of a subject projection for obtaining said projection data, and comparing the value of the calculated evaluation function with a predetermined threshold value for evaluating said prepared correction table;

interpolating means for deriving said prepared correction table that is not continuously changed between said previous and subsequent projections from said prepared correction tables for said previous and subsequent projections through the effect of the interpolation, based on the evaluated result of said evaluating means; and

means for reconstructing the distribution of X-ray attenuation coefficients of said object.

23. An apparatus for cone-beam X-ray computed tomography as claimed in claim 22, wherein if on the plane of said two-dimensional X-ray detector for detecting said projection data a u axis is set to the rotating direction of said scanner, a v axis is set in parallel to the center axis of rotation of said scanner, and said prepared correction table of said projection in the k direction is represented by U(u, v, k) and V(u, v, k) (k=1, 2, . . . , N_k, u=1, 2, . . . , N_i, v=1, 2, . . . , N_j), said evaluating means comprises means for calculating a correlation between said prepared correction tables adjacent to each other in said k direction for evaluating continuity of said prepared correction tables adjacent to each other.

24. An apparatus for cone-beam X-ray computed tomography as claimed in claim 22, if k is any natural number, said evaluating means comprises means for calculating valU(u, v, k) and valV(u, v, K) (if k+l is overflown out of the range of k=1, 2, . . . , N_k in the calculation, the tail of k is circulated to the head of k.) according to the expression of;

##EQU18##

25. An apparatus for cone-beam X-ray computed tomography as claimed in claim 22, wherein if on the plane of said two-dimensional X-ray detector for detecting said projection data a u axis is set to the rotating direction of said scanner, a v axis is set in parallel to the center axis of rotation of said scanner, and said prepared correction table of said projection in the k direction is represented by U(u, v, k) and V(u, v, k) (k=1, 2, . . . , N_k, u=1, 2, . . . , N_i, v=1, 2, . . . , N_j), said evaluating means comprises means for calculating a correlation between said prepared correction tables adjacent to each other in said k direction for evaluating continuity of said prepared correction tables adjacent to each other, and wherein it is determined if the values calculated by said evaluating means meet a specific condition.

26. An apparatus for cone-beam X-ray computed tomography as claimed in claim 22, wherein if on the plane of said two-dimensional X-ray detector for detecting said projection data a u axis is set to the rotating direction of said scanner, a v axis is set in parallel to the center axis of rotation of said scanner, and said prepared correction table of said projection in the k direction is represented by U(u, v, k) and V(u, v, k) (k=1, 2, . . . , N_k, u=1, 2, . . . , N_i, v=1, 2, . . . , N_j), said evaluating means comprises means for calculating a correlation between said prepared correction tables adjacent to each other in said k direction for evaluating continuity of said prepared correction tables adjacent to each other, and wherein it is determined if the values calculated by said evaluating means are greater than a specific threshold value.

27. An apparatus for cone-beam X-ray computed tomography as claimed in claim 22, wherein if k is any natural number, said evaluating means comprises means for calculating valU(u, v, k) and valV(u, v, k) (if k+l is overflown out of the range of k=1, 2, . . . , N_k in the calculation, the tail is circulated to the head of k) according to the expression of:

##EQU19## and wherein if the value calculated by said evaluating means is greater than a specific threshold value of said projection in the serial k (=k₀, k₀ +1, . . . , K₀ +n) direction, said prepared correction table of said projection in the direction of k=k₀ +1, . . . , k=k₀ +n-1 is deleted from said means for saving correction table, and said interpolating means operates to derive the correction table of said projection in the direction of k=k₀ +1, . . . , k=k₀ +n-1 from said prepared correction table in said k₀ and said (k₀ +n) directions through the effect of the interpolation.

28. An apparatus for cone-beam X-ray computed tomography as claimed in claim 27, wherein said interpolating means comprises means for deriving said prepared correction table U(u, v, k) and V(u, v, k) (k=k₀ +1, . . . , k₀ +n-1, k₁ <

k₀ +1, k₀ +n-1<

k₂) deleted from said means for saving correction table from said prepared correction table U(u, v, k₁) and V(u, v, k₁) immediately before said deleted correction table in said projection direction and said prepared correction table U(u, v, k₂) and V(u, v, k₂) immediately after said deleted correction table in said projection direction according to the expression of;

##EQU20##

29. A method for cone-beam X-ray computed tomography in which a scanner mounting an X-ray source for applying a cone-beam X ray and a scanner having a two-dimensional X-ray detector for detecting an X-ray is rotated around an object to measure omnidirectional projection data around said object, an image geometric distortion of said projection data is corrected by referring a correction table, and a distribution of X-ray attenuation coefficients is reconstructed by using said distortion-corrected projection data, comprising the steps of:

preparing said correction table in advance and saving said table in means for saving correction table;

calculating an evaluation function for evaluating continuity of said prepared correction table between the previous and the subsequent projections of a subject projection for obtaining said projection data;

deleting said prepared correction table that is not continuously changed between the previous and the subsequent projections from said means for saving correction table based on the calculated value of said evaluation function, replacing said correction table that is not continuously changed with a correction table obtained by interpolating said prepared correction tables of the previous and the subsequent projections, and saving the replaced correction table in said means for saving correction table;

detecting a location of a view field boundary of said two-dimensional X-ray detector from said distortion-corrected projection data by using the correction table saved in said means for saving correction table; and

extrapolating projection data of a predetermined width from said location of the view field boundary detected in the step of detecting; and

reconstructing said X-ray attenuation coefficients distribution from the projection data extrapolated in the step of extrapolating.

30. An apparatus for X-ray computed tomography having a scanner mounting an X-ray source for applying a cone-beam X ray and a two-dimensional X-ray detector for detecting an X-ray, means for rotating said scanner around an object, means for measuring omnidirectional projection data around said object, and means for measuring omnidirectional projection data around said object, for reconstructing an X-ray attenuation distribution, comprising:

means for preparing a correction table on which an image geometric distortion of said projection data is corrected, said preparing means having means for preparing said correction table of each projection for obtaining said projection data in advance and means for saving said prepared correction table;

evaluating means for calculating an evaluation function for evaluating continuity of said correction table between the previous and the subsequent projections of each subject projection for obtaining said projection data and comparing the calculate value of said evaluation function with a predetermined threshold value for evaluating said prepared correction table;

interpolating means for replacing said prepared correction table that is not continuously changed between said previous and subsequent projections with a correction table interpolated from the prepared correction tables of said previous and subsequent projections, based on the evaluating result of said evaluation function;

means for correcting an image geometric distortion of said projection data by referring to said correction table;

means for detecting a location of a view field boundary of said two-dimensional X-ray detector about said distortion-corrected projection data;

extrapolating means for generating said projection data of a predetermined width from said location of the view field boundary detected by said detecting means; and

means for reconstructing said X-ray attenuation distribution of said object from said projection data extrapolated in said extrapolating means.