System and method for motion estimation and compensation in helical computed tomography

• US 10,339,678 B2
• Filed: 11/16/2017
• Issued: 07/02/2019
• Est. Priority Date: 11/16/2016
• Status: Active Grant

• Current Assignee: iTomography Corporation, University of Central Florida Research Foundation Inc. (State University System of Florida)
• Original Assignee: University of Central Florida Research Foundation Inc. (State University System of Florida)
• Inventors: Katsevich, Alexander, Frenkel, Michael
• Primary Examiner(s): Fitzpatrick, Atiba O
• Application Number:
  US15/815,662
• Publication Number:
  US 20180137657A1
• Time in Patent Office: 593 Days
• Field of Search: None
• 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/10081 Computed x-ray tomography [CT]
• CPC Class Codes: G06T 5/70 Denoising; Smoothing
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G06T 7/20 Analysis of motion motion estimation for coding, decoding, compressing or decompressing digital video signals H04N19/43, H04N19/51

Abstract

An improved system and method for estimating and compensating for motion by reducing motion artifacts produced during image reconstruction from helical computed tomography (CT) scan data. In a particular embodiment, the reconstruction may be based on helical partial angle reconstruction (PAR).

Claims

1. A method for estimating and compensating for motion by reducing motion artifacts in an image reconstruction from helical computed tomography (CT) scan data of an object of interest, the method comprising:

collecting helical computer tomography (CT) scan data of an object of interest, wherein the scan data is acquired using a radiation source to generate a cone beam and a radiation detector to detect the cone beam;

selecting a plurality of center-points along a trajectory of the radiation source;

identifying a plurality of pairs of sections along the trajectory of the radiation source, wherein each of the plurality of pairs of sections is associated with one of the plurality of center-points and wherein a first section of each of the pairs of sections and a second section of each of the pairs of sections are positioned on opposite sides of the center-point;

selecting a subset of the plurality of pairs of sections;

reconstructing, for each pair of the subset, a first partial image from the scan data of the first section and a second partial image from the scan data of the second section;

performing image registration of the first partial image and the second partial image for each pair of the subset to estimate a deformation that transforms the first partial image into the second partial image, wherein the deformation is representative of motion of the object of interest during the scan; and

generating a motion compensated image by reconstructing the object of interest using the scan data and the estimated deformation.

2. The method of claim 1, wherein reconstructing, for each pair of the subset, a first partial image from the scan data of the first section and a second partial image from the scan data of the second section further comprises, reconstructing the first partial image and the second partial image in a region of interest having a substantial overlap.

3. The method of claim 2, further comprising identifying at least one chord between the first section and the second section and wherein reconstructing, for each pair of the subset, a first partial image from the scan data of the first section and a second partial image from the scan data of the second section further comprises, reconstructing the first partial image and the second partial image at a set of points in a neighborhood of at least one chord.

4. The method of claim 1, wherein the first section of each of the plurality of pairs of sections and the second section of each of the plurality of pairs of sections are separated from each other by an angular distance equal to π

.

5. The method of claim 1, wherein the plurality of center-points are equally spaced along the trajectory of the helical CT radiation source.

6. The method of claim 1, further comprising acquiring a motion information signal and wherein a spacing between each of the plurality of center-points and a location of each of the plurality of center-points is based upon the motion information signal.

7. The method of claim 1, further comprising acquiring a motion information signal and wherein a length of the first section and a length of the second section of the subset is increased or decreased based upon the motion information signal.

8. The method of claim 1, wherein performing image registration of the first partial image and the second partial image for each pair of the subset to estimate a deformation that transforms the first partial image into the second partial image further comprises estimating a deformation of a point in the object of interest at a point in time, which is inside an interval of time during which the point was irradiated by the cone beam of the radiation source.

9. The method of claim 1, wherein generating a motion compensated image by reconstructing the object of interest using the scan data and the estimated deformation further comprising, estimating a deformation of the object of interest by interpolating the estimates of deformations for the subset.

10. The method of claim 1, where reconstructing, for each pair of sections of the subset, a first partial image from the scan data of the first section and a second partial image from the scan data of the second section is based upon Local Tomography (LT), exact reconstruction or quasi-exact reconstruction techniques.

11. The method of claim 2, further comprising identifying at least one chord between the first section and the second section for each pair of the subset and reconstructing the first partial image and the second partial image of each pair in a neighborhood of the identified chord to generate a two-dimensional first partial image and a two-dimensional second partial image and performing image registration using the two-dimensional first partial image and the two-dimensional second partial image.

12. The method of claim 2, further comprising identifying a chord surface between the first section and the second section for each pair of the subset and reconstructing the first partial image and the second partial image of each pair in a neighborhood of the identified chord surface to generate a three-dimensional first partial image and a three-dimensional second partial image and performing image registration using the three-dimensional first partial image and the three-dimensional second partial image.

13. A system for estimating and compensating for motion by reducing motion artifacts in an image reconstruction from helical computed tomography (CT) scan data of an object of interest, the method comprising:

a memory for storing a helical computer tomography (CT) scan data of an object of interest;

a data processor for estimating and compensating for motion by reducing motion artifacts in an image reconstruction from the helical computed tomography (CT) scan data of an object of interest, wherein the data processor is adapted for performing the following operations;

loading the helical CT scan data from the memory;

selecting a plurality of center-points along a trajectory of the radiation source;

identifying a plurality of pairs of sections along the trajectory of the radiation source, wherein each of the plurality of pairs of section is associated with one of the plurality of center-points and wherein a first section of each of the pairs of sections and a second section of each of the pairs of sections are positioned on opposite sides of the center-point;

selecting a subset of the plurality of pairs of sections;

reconstructing, for each pair of the subset, a first partial image from the scan data of the first section and a second partial image from the scan data of the second section;

performing image registration of the first partial image and the second partial image for each pair of the subset to estimate a deformation that transforms the first partial image into the second partial image, wherein the deformation is representative of motion of the object during the scan; and

generating a motion artifact compensated image by reconstructing the object of interest using the scan data and the estimated deformation.

15. The media of claim 13, wherein reconstructing, for each pair of the subset, a first partial image from the scan data of the first section and a second partial image from the scan data of the second section further comprises, reconstructing the first partial image and the second partial image in a region of interest having a substantial overlap.

17. The media of claim 13, wherein performing image registration of the first partial image and the second partial image for each pair of the subset to estimate a deformation that transforms the first partial image into the second partial image further comprises estimating a deformation of a point in the object of interest at a point in time, which is inside an interval of time during which the point was irradiated by the cone beam of the radiation source.

18. The media of claim 13, wherein generating a motion compensated image by reconstructing the object of interest using the scan data and the estimated deformation further comprising, estimating a deformation of the object of interest by interpolating the estimate of deformations for the subset.

14. One or more non-transitory computer-readable media having computer-executable instructions for performing a method of estimating and compensating for motion by reducing motion artifacts in an image reconstruction from helical computed tomography (CT) scan data of an object of interest, the method comprising:

collecting helical computer tomography (CT) scan data of an object of interest, wherein the scan data is acquired using a radiation source to generate a cone beam and a radiation detector to detect the cone beam;

selecting a plurality of center-points along a trajectory of the radiation source;

identifying a plurality of pairs of sections along the trajectory of the radiation source, wherein each of the plurality of pairs of sections is associated with one of the plurality of center-points and wherein a first section of each of the pairs of sections and a second section of each of the pairs of sections are positioned on opposite sides of the center-point;

selecting a subset of the plurality of pairs of sections;

reconstructing, for each pair of the subset, a first partial image from the scan data of the first section and a second partial image from the scan data of the second section;

performing image registration for each pair of the first partial image and the second partial image of the subset to estimate a deformation that transforms the first partial image into the second partial image, wherein the deformation is representative of motion of the object of interest during the scan; and

generating a motion artifact compensated image by reconstructing the object of interest using the scan data and the estimated deformation.

16. The media of claim 14, further comprising identifying at least one chord between the first section and the second section and wherein reconstructing, for each pair of the subset, a first partial image from the scan data of the first section and a second partial image from the scan data of the second section further comprises, reconstructing the first partial image and the second partial image at a set of points in a neighborhood of the at least one chord.

19. The media of claim 14, further comprising identifying a chord between the first section and the second section for each pair of the subset and reconstructing the first partial image and the second partial image of each pair in a neighborhood of the identified chord to generate a two-dimensional first partial image and a two-dimensional second partial image and performing image registration using the two-dimensional first partial image and the two-dimensional second partial image.

20. The media of claim 14, further comprising identifying a chord surface between the first section and the second section for each pair of the subset and reconstructing the first partial image and the second partial image of each pair on the identified chord surface to generate a three-dimensional first partial image a three-dimensional second partial image and performing image registration using the three-dimensional first partial image and the three-dimensional second partial image.