CORRECTING AND RECONSTRUCTING X-RAY IMAGES USING PATIENT MOTION VECTORS EXTRACTED FROM MARKER POSITIONS IN X-RAY IMAGES

• US 20110123081A1
• Filed: 02/04/2010
• Published: 05/26/2011
• Est. Priority Date: 11/25/2009
• Status: Active Grant

• Current Assignee: Dental Imaging Technologies Corporation (Envista Holdings Corporation)
• Original Assignee: Dental Imaging Technologies Corporation (Envista Holdings Corporation)
• Inventors: Bruning, Horst F., Gouzeev, Serguei, Charamisinau, Ivan, Sebok, David
• Granted Patent: US 9,826,942 B2
• Time in Patent Office: Days
• US Class Current: 382/131
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: A61B 6/4085 Cone-beams
• CPC Class Codes: A61B 6/501 for diagnosis of the head, e.g. neuroimaging or craniography
• CPC Class Codes: A61B 6/5264 due to motion
• 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 2207/10116 X-ray image
• CPC Class Codes: G06T 2207/30036 Dental; Teeth
• CPC Class Codes: G06T 2211/412 Dynamic
• CPC Class Codes: G06T 7/33 using feature-based methods

Abstract

A method and a system for correcting and reconstruing a plurality of projection images based on detected patient motion. The method includes obtaining a plurality of projection images generated by a scanner. Each of the plurality of projection images includes at least three markers. Each of the at least three markers has a measured three-dimensional position and measured positions on a detector panel of the scanner in a first dimension and a second dimension. The method further includes determining a position error vector for each of the plurality of projection images based on the at least three markers in each of the plurality of projection images, the position error vector defining patient motion in the projection image. The method finally includes combining each position error vector for each of the plurality of projection images with geometric parameters associated with the scanner to derive a projection matrix for each of the plurality of projection images, and generating reconstructed images corrected for patient motion from the plurality of projection images and the projection matrix for each of the plurality of projection images.

Claims

1. A method for correcting and reconstructing a plurality of projection images based on detected patient motion, the method executed by an imaging system including a computer with an electronic processing unit and a memory module storing a projection matrix calculation module and a reconstruction module, both executable by the electronic processing unit, the method comprising:

obtaining, at the computer, a plurality of projection images generated by a scanner, each of the plurality of projection images including at least three markers, each of the at least three markers having a measured three-dimensional position and measured positions on a detector panel of the scanner in a first dimension and a second dimension;

determining, with the electronic processing unit, a position error vector for each of the plurality of projection images based on the at least three markers in each of the plurality of projection images, the position error vector defining patient motion in the projection image;

combining each position error vector for each of the plurality of projection images with geometric parameters associated with the scanner to derive a projection matrix for each of the plurality of projection images; and

generating reconstructed images corrected for patient motion from the plurality of projection images and the projection matrix for each of the plurality of projection images.

2. The method of claim 1, wherein generating reconstructed images corrected for patient motion includes performing a back-projection reconstructing process.

3. The method of claim 1, wherein generating reconstructed images corrected for patient motion includes performing an iterative reconstructing process.

4. The method of claim 1, wherein each of the at least three markers has a measured three-dimensional position X_i, Y_i, Z_i and wherein determining the position error vector for each projection image includes calculating expected positions of each of the at least three markers on the detector panel in a first dimension U_e and a second dimension V_e and determining motion parameters Δ

X, Δ

Y, Δ

Z, α

, β

, and γ

defining patient motion occurring in the image based on variations between the excepted positions and the measured positions of each of the at least three markers.

5. The method of claim 4, wherein determining the motion parameters patient motion vector includes

(a) assuming a value of Δ

X, Δ

Y, and γ

,

(b) calculating a value for Δ

Z, α

, and β

using the assumed values of Δ

X, Δ

Y, and γ and

at least one equation defining positions in the second dimension, and

(c) calculating a value for Δ

X, Δ

Y, and γ

using the calculated values of Δ

X, Δ

Y, and γ and

at least one equation defining positions in the first dimension.

6. The method of claim 5, wherein determining the motion parameters further includes repeating steps (b) through (c) at least once.

7. The method of claim 4, wherein determining the motion parameters further includes approximating sines and cosines in rotational transformations using a Taylor expansion consisting of a linear first term and an estimated residual.

8. The method of claim 7, wherein determining the motion parameters further includes updating the residuals of the approximated sines and cosines during repetition of steps (b) through (c).

9. The method of claim 7, wherein approximating the sines and cosines in the rotational transformations using the Taylor expansion includes using an estimated residual that is constant.

10. The method of claim 1, wherein obtaining the image including the at least three markers includes obtaining the image including the at least three markers, each of the at least three markers having measured positions on a detector panel of the scanner in a first dimension U_m tangential to a circle of rotation of a gantry of the scanner and a second dimension V_m parallel to an axis of rotation of the gantry.

11. A system for correcting and reconstructing a plurality of projection images based on detected patient motion, the system comprising:

a scanner; and

a computer including an electronic processing unit, a memory module, and an input/output interface, where the memory is configured to store a projection matrix calculation module and a reconstruction module, both executed by the electronic processing unit;

wherein the scanner is configured to obtain and send to the computer a plurality of projection images generated by the scanner, each of the plurality of projection images including at least three markers, each of the at least three markers having a measured three-dimensional position and measured positions on a detector panel of the scanner in a first dimension and a second dimension,

wherein the electronic processing unit is configured to determine a position error vector for each of the plurality of projection images based on the at least three markers in each of the plurality of projection images, the position error vector defining patient motion in the projection image,

wherein the electronic processing unit is configured to combine each position error vector for each of the plurality of projection images with geometric parameters associated with the scanner to derive a projection matrix for each of the plurality of projection images,

wherein the electronic processing unit is configured to generate reconstructed images corrected for patient motion from the plurality of projection images and the projection matrix for each of the plurality of projection images.

12. The system according to claim 11, wherein the electronic processing unit is configured to generate reconstructed images corrected for patient motion by performing a back-projection reconstructing process.

13. The system according to claim 11, wherein the electronic processing unit is configured to generate reconstructed images corrected for patient motion by performing an iterative reconstructing process.

14. The system according to claim 11, wherein each of the at least three markers has a measured three-dimensional position X_i, Y_i, Z_i and wherein the electronic processing unit is configured to determine the position error vector for each projection image by calculating expected positions of each of the at least three markers on the detector panel in a first dimension U_e and a second dimension V_e and determining motion parameters Δ

X, Δ

Y, Δ

Z, α

, β

, and γ

defining patient motion occurring in the image based on variations between the excepted positions and the measured positions of each of the at least three markers.

15. The system of claim 14, wherein the electronic processing unit is configured to determine the motion parameters patient motion vector by:

(a) assuming a value of Δ

X, Δ

Y, and γ

,

(b) calculating a value for Δ

Z, α

, and β

using the assumed values of Δ

X, Δ

Y, and γ and

at least one equation defining positions in the second dimension, and

(c) calculating a value for Δ

X, Δ

Y, and γ

using the calculated values of Δ

X, Δ

Y, and γ and

at least one equation defining positions in the first dimension.

16. The system of claim 15, wherein the electronic processing unit is further configured to determine the motion parameters by repeating steps (b) through (c) at least once.

17. The system of claim 14, wherein the electronic processing unit is further configured to determine the motion parameters by approximating sines and cosines in rotational transformations using a Taylor expansion consisting of a linear first term and an estimated residual.

18. The system of claim 17, wherein the electronic processing unit is further configured to determine the motion parameters by updating the residuals of the approximated sines and cosines during repetition of steps (b) through (c).

19. The system of claim 17, wherein the electronic processing unit is configured to approximate the sines and cosines in the rotational transformations using the Taylor expansion by using an estimated residual that is constant.

20. The system of claim 11, wherein the scanner is configured to obtain the image including the at least three markers by obtaining the image including the at least three markers, each of the at least three markers having measured positions on a detector panel of the scanner in a first dimension U_m tangential to a circle of rotation of a gantry of the scanner and a second dimension V_m parallel to an axis of rotation of the gantry.