Method for computed tomography motion estimation and compensation

• US 8,224,056 B2
• Filed: 12/15/2009
• Issued: 07/17/2012
• Est. Priority Date: 12/15/2009
• Status: Active Grant

• Current Assignee: General Electric Company
• Original Assignee: General Electronics Company
• Inventors: Pack, Jed Douglas, Edic, Peter Michael, Claus, Bernhard Erich Hermann, Iatrou, Maria, Miller, James Vradenburg
• Primary Examiner(s): Wu, Jingge
• Assistant Examiner(s): Drennan, Barry
• Application Number:
  US12/638,693
• Publication Number:
  US 20110142313A1
• Time in Patent Office: 945 Days
• Field of Search: None
• US Class Current: 382/131
• 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/10081 Computed x-ray tomography [CT]
• CPC Class Codes: G06T 2207/30048 Heart; Cardiac
• CPC Class Codes: G06T 2207/30101 Blood vessel; Artery; Vein; Vascular
• CPC Class Codes: G06T 2207/30172 Centreline of tubular or elongated structure
• CPC Class Codes: G06T 2211/412 Dynamic
• CPC Class Codes: G06T 7/11 Region-based segmentation
• CPC Class Codes: G06T 7/174 involving the use of two or more images
• CPC Class Codes: G06T 7/246 using feature-based methods, e.g. the tracking of corners or segments

Abstract

A method and system for motion estimation and compensation are disclosed. Initially, a set of one or more initial images is reconstructed using acquired imaging data. Further, one or more regions of interest are identified in this set of reconstructed initial images. At least a set of filters is applied to the identified regions of interest to generate a sequence of filtered images. Particularly, each of the filtered images in the generated sequence of filtered images includes data acquired near a different reference point. Subsequently, a motion path corresponding to each region of interest is determined based on one or more correspondences in the sequence of filtered images.

Claims

1. A method for estimating motion of an object of interest, comprising:

reconstructing a set of one or more initial images using acquired data;

identifying one or more regions of interest in the set of reconstructed initial images;

applying at least a set of directional Fourier filters to the identified regions of interest to generate a sequence of filtered images constrained in Fourier space, wherein each of the filtered images in the sequence comprises data acquired near a different reference point; and

estimating a motion path for each region of interest based on one or more correspondences in the sequence of filtered images.

2. The method of claim 1, wherein the region of interest comprises at least a portion of the object of interest.

3. The method of claim 2, wherein the object of interest comprises a coronary artery, a coronary valve, a heart chamber, the whole heart, or a lung nodule.

4. The method of claim 1, wherein reconstructing the initial image is performed using a short scan technique, a Feldkamp-Davis-Kress reconstruction technique or a “

Butterfly”

reconstruction.

5. The method of claim 1, wherein identifying the one or more regions of interest comprises employing a segmentation technique, a matched-filter based approach, an Eigen analysis based approach, or combinations thereof.

6. The method of claim 1, wherein the reference point comprises one of a particular instant of time, a view angle corresponding to a source location, or an angular space in a Fourier domain.

7. The method of claim 1, wherein identifying the one or more regions of interest comprises generating a representation of a centerline path for at least a portion of the object of interest.

8. The method of claim 7, wherein applying the set of filters comprises designing a set of filters that suppress the influence of data that was not acquired near a particular reference point.

9. The method of claim 1, wherein estimating a motion path, comprises:

identifying the displacements of one or more portions of an object of interest contained in a region of interest in the sequence of filtered images; and

determining a motion path based on the identified displacements.

10. The method of claim 1, wherein reconstructing the set of one or more initial images comprises reconstructing two or more initial images.

11. The method of claim 10, wherein estimating the motion path of a portion of an object of interest comprises:

selecting one or more pairs of frames from the sequence of filtered images wherein the two frames of each selected pair are generated by applying substantially the same filter to two different initial reconstructions;

identifying the relative displacement of a portion of the object of interest between the two frames in each selected pair of frames; and

determining the motion path based on the identified relative displacements.

12. The method of claim 1 wherein estimating a motion path comprises one of using a model for the motion that encourages smoothness, incorporates confidence information for each component of each displacement based on directional properties of the filter, or combinations thereof.

13. The method of claim 12, wherein estimating a motion path further comprises using a string model, a normalized correlation, or a combination thereof.

14. The method of claim 1, wherein the set of filters operate in a Fourier domain, a spatial domain, a wavelet domain, a time-frequency domain, a three-dimensional domain, a two-dimensional domain, a projection domain, or combinations thereof.

15. The method of claim 1, wherein applying the set of filters comprises designing the set of filters, wherein an impulse response of each of the set of filters comprises a high spatial compactness, a high angular compactness, a correspondence with a particular reference point, or combinations thereof.

16. The method of claim 1, further comprising compensating for motion of the object of interest.

17. The method of claim 16, wherein compensating for motion of the object of interest comprises warping a reconstruction grid using the estimated motion path before backprojecting the acquired data.

18. An imaging system, comprising:

a data acquisition system (DAS) for acquiring an initial image; and

an image processor for;

reconstructing a set of one or more initial images using acquired data;

identifying one or more regions of interest in the set of reconstructed initial images;

applying at least a set of filters to the identified regions of interest to generate a sequence of filtered images, wherein each of the filtered images in the sequence comprises data acquired near a different reference point; and

estimating a motion path for each region of interest based on one or more correspondences in the sequence of filtered images;

wherein each filter in the set of filters comprises a filter constrained to lie in a conjugate pair of opposing sections in Fourier space and in a band of Fourier space centered at the origin.

19. The imaging system of claim 18, wherein the reference point comprises one of a particular instant of time, a view angle corresponding to a source location, or an angular space in a Fourier domain.

20. The imaging system of claim 18, wherein the image processor further compensates for the estimated motion.

21. The imaging system of claim 18, wherein the data acquisition system comprises a computed tomography (CT) system, a multi-source imaging system, multi-detector imaging system, an X-Ray system, positron emission tomography (PET) scanner, a single photon emission computed tomography (SPECT) scanner, a microscope, a digital camera, a charge coupled device, or combinations thereof.

22. A computed tomography (CT) system, comprising:

at least one radiation source, for generating a conical beam for imaging an object, wherein the at least one radiation source and the object move relative to each other around a determined axis;

a detector assembly coupled to the at least one radiation source; and

a computing device coupled to detector assembly for estimating motion of the object,

wherein the computing device includes one or more image processors for;

reconstructing a set of one or more initial images using acquired data;

identifying one or more regions of interest in the set of reconstructed initial images;

applying at least a set of directional Fourier filters to the identified regions of interest to generate a sequence of filtered images constrained in Fourier space, wherein each of the filtered images in the sequence comprises data acquired near a different reference point; and

estimating a motion path for each region of interest based on one or more correspondences in the sequence of filtered images.

23. The CT system of claim 22, wherein the image processor further compensates for the estimated motion.

24. The CT system of claim 22, wherein the reference point comprises one of a particular instant of time, a view angle corresponding to a source location, or an angular space in a Fourier domain.

25. A method for estimating motion of an object of interest, comprising:

reconstructing a set of one or more initial images using acquired data;

identifying one or more regions of interest in the set of reconstructed initial images;

applying at least a set of filters to the identified regions of interest to generate a sequence of filtered images, wherein each of the filtered images in the sequence comprises data acquired near a different reference point; and

estimating a motion path for each region of interest based on one or more correspondences in the sequence of filtered images;

wherein reconstructing the set of one or more initial images comprises reconstructing two or more initial images; and

wherein estimating the motion path of a portion of the object of interest comprises;

selecting one or more pairs of frames from the sequence of filtered images, wherein the two frames of each selected pair are generated by applying substantially the same filter to two different initial reconstructions;

identifying the relative displacement of a portion of the object of interest between the two frames in each selected pair of frames; and

determining the motion path based on the identified relative displacements.