AUTOMATIC SYSTEM CALIBRATION METHOD OF X-RAY CT

• US 20180068467A1
• Filed: 03/21/2016
• Published: 03/08/2018
• Est. Priority Date: 03/20/2015
• Status: Active Grant

• Current Assignee: Rensselaer Polytechnic Institute
• Original Assignee: Rensselaer Polytechnic Institute
• Inventors: Wang, Ge, Chen, Mianyi, Xi, Yan, Cong, Wenxiang
• Granted Patent: US 10,643,354 B2
• Time in Patent Office: Days
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: A61B 6/5205 involving processing of raw data to produce diagnostic data
• CPC Class Codes: A61B 6/582 Calibration
• 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: G06T 11/006 Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
• CPC Class Codes: G06T 2207/10116 X-ray image
• CPC Class Codes: G06T 2207/20 Special algorithmic details
• CPC Class Codes: G06T 2211/412 Dynamic
• CPC Class Codes: G06T 2211/424 Iterative

Abstract

Systems and methods for geometric calibration and image reconstruction in computed tomography (CT) scanning using iterative reconstruction algorithms are provided. An iterative reconstruction algorithm can be used to reconstruct an improved image, and then the improved image can be used to adjust inaccurate parameters by using a Locally Linear Embedding (LLE) method. Adjusted parameters can then be used to reconstruct new images, which can then be used to further adjust the parameters. The steps of this iterative process can be repeated until a quality threshold is met.

Claims

1. A method of reconstructing a computed tomography (CT) image, the method comprising:

i) obtaining an initial CT image;

ii) performing a reconstruction algorithm on the initial CT image to obtain a reconstructed image of the initial CT image;

iii) using the reconstructed image to adjust one or more parameters associated with the image;

iv) using the adjusted one or more parameters to perform the reconstruction algorithm on the reconstructed image to obtain an updated reconstructed image; and

v) repeating steps iii)-iv), using the most-recently updated reconstructed image in each repeated step iii), and updating the same one or more parameters in each repeated step iii) until a threshold value of a predetermined characteristic is met.

2. The method according to claim 1, wherein the reconstruction algorithm is based on Total Variation (TV).

3. The method according to claim 1, wherein using the reconstructed image to adjust one or more parameters associated with the image comprises using a Locally Linear Embedding (LLE) method on each of the one or more parameters.

4. The method according to claim 3, wherein the LLE method includes:

a) finding the K nearest neighbors for each high-dimensional point in terms of the Euclidean distance;

b) representing each point linearly with its neighbors and calculating weight coefficients for its neighbors; and

c) mapping high dimensional data with the weight coefficients to a low dimensional representation on an intrinsic manifold.

5. The method according to claim 3, wherein the LLE method includes:

finding the K nearest neighbors of a data vector b_i in its dataset vectors {tilde over (b)}_ij according to the Euclidean distance;

d_ij=∥

b_i−

{tilde over (b)}_ij∥

₂².

6. The method according claim 3, wherein the LLE method (further) includes:

with the K nearest vectors, representing the original data vector linearly with its neighboring vectors;

8. The method according to claim 3, wherein the LLE method (further) includes, calculating, with the weight coefficients W=(w_ik), the global internal coordinate Y=(y_i) by solving the equation:

9. The method according to claim 8, wherein the solution of the equation that can be solved to calculate Y=(y_i) is given by the bottom d+1 eigenvectors of the generalized eigenvalue problem:

MY=λ

Y,

where λ

is the eigenvalue, and M=(I−

W)^T(I−

W).

11. The method according to claim 3, wherein the number of nearest neighbors used in the LLE method is at least 2.

18. The method according to claim 1, wherein the predetermined characteristic, of which a threshold value must be met to stop repeating steps iii)-iv) is universal quality index (UQI).

19. The method according to claim 18, wherein the threshold value is at least 0.6.

23. The method according to claim 18, wherein the threshold value is at least 0.95.

24. The method according to claim 1, wherein the one or more parameters includes at least one of projection angle, source-to-object distance (SOD), object-to-detector distance (ODD), detector offset, detector tilt angle, object x-coordinate offset, object y-coordinate offset, and projection angular error.

30. The method according to claim 1, wherein obtaining the initial CT image comprises obtaining an initial parameter vector including initial value(s) for the one or more parameters associated with the image, and

wherein performing the reconstruction algorithm comprises updating the parameter vector with updated values of the one or more parameters based on the most recently reconstructed image.

32. The method according to claim 1, wherein the reconstruction algorithm comprises using an Ordered Subset Simultaneous Algebraic Reconstruction Technique (OS-SART).

34. The method according to claim 1, wherein each of step ii), step iii), step iv), and step v) is performed by a processor.

38. A computer readable medium having computer-executable instructions for performing the method according to claim 1.

39. A CT system, comprising:

an X-ray source;

a detector for detecting X-ray radiation from the X-ray source; and

a computer having the computer readable medium according to claim 38.

7. (canceled)

10. (canceled)

13-16. -16. (canceled)

20-22. -22. (canceled)

25-29. -29. (canceled)

31. (canceled)

33. (canceled)

35-37. -37. (canceled)

40. (canceled)