Localized motion artifact reduction in projection imaging

• US 5,287,276 A
• Filed: 10/22/1992
• Issued: 02/15/1994
• Est. Priority Date: 11/19/1990
• Status: Expired due to Term

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Crawford, Carl R., Ritchie, Cameron J.
• Primary Examiner(s): Envall, Jr., Roy N.
• Assistant Examiner(s): Huntley, David
• Application Number:
  US07/964,335
• Time in Patent Office: 481 Days
• Field of Search: 364/413.13, 364/413.14, 364/413.15, 364/413.16, 364/413.17, 364/413.19, 128/653.2
• US Class Current: 378/4
• CPC Class Codes: A61B 5/113 occurring during breathing
• CPC Class Codes: A61B 5/721 using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: A61B 6/5264 due to motion
• CPC Class Codes: A61B 6/589 Setting distance between source unit and patient
• 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 2211/412 Dynamic
• CPC Class Codes: G06T 2211/421 Filtered back projection [FBP]

Abstract

An x-ray CT scanner acquires projection data from a series of projections during a scan of a patient'"'"'s chest Movement of the patient'"'"'s chest due to respiration is also sensed during the scan and this acquired motion data is employed along with an warping function which models chest motion to calculate factors which correct the acquired projection data and reduce motion artifacts in an image produced by back projecting the acquired projection data.

Claims

1. A method for producing an image of a moving patient composed of a plurality of voxels, each voxel having reference coordinate defined at a reference time while the patient is in motion, and a displaced coordinate which at the reference time may equal the reference coordinate, but at other times may not equal the reference coordinate, the steps comprising:

performing a scan in which a plurality of projections of a projection data set are acquired at different times and throughout a range of different projection angles (θ

), each projection measuring physical characteristic of the voxels at their displaced coordinates;

acquiring a motion parameter with each projection of the projection set, each acquired motion parameter being indicative of the movement of the patient as the projection set is acquired; and

back projecting the set of projections to produce the image of the patient using a back projection formula which is modified to account for the displacement of the voxels from their reference coordinates at the time of the projection is acquired, the modification being according to a predetermined warping function that relates the displaced coordinates of the voxels to their reference coordinates as a function of the motion parameter;

whereby artifacts produced in the image by movement of the patient from projection-to-projection during the scan are reduced.

2. The method as recited in claim 1 in which the motion parameter acquired with each set of projection data indicates the position of the anterior chest wall of the patient.

3. The method as recited in claim 1 wherein the step of backprojecting includes for each voxel at a reference coordinate and for each projection:

a) employing the warping function to determine from the motion parameter the displaced coordinates of that voxel at the time of the acquisition of the projection; and

b) backprojecting the voxel at the reference coordinate but employing the displaced coordinate to identify data from the projection used for the backprojection.

4. The method as recited in claim 3 in which the motion parameter acquired with each projection indicates the position of the anterior chest wall of the patient and the warping function provides a displacement coordinate for each voxel based on the position of the anterior chest wall, the reference coordinate of the voxel, and a motion reference value which indicates the position of the anterior chest wall of the patient at a predetermined point in the patient'"'"'s respiratory cycle.

5. The method as recited in claim 1 in which the projection data is x-ray attenuation values.

6. The method as recited in claim 5 in which the motion of the patient is due to respiration and the warping function is calculated based on a geometric model of the motion of the patient'"'"'s chest cavity during respiration and the acquired motion parameter which indicates a position of the patient as the projection set is acquired.

7. The method as recited in claim 1 in which the set of projection data is acquired with parallel x-ray beams.

8. The method as recited in claim 1 in which the set of projection data is acquired with a curved array of x-ray detectors that receive a fan-beam of x-rays.

9. The method as recited in claim 1 in which the set of projection data is acquired with a flat array of x-ray detectors that receive a fan-beam of x-rays.