Cone beam CT apparatus using truncated projections and a previously acquired 3D CT image

• US 7,440,535 B2
• Filed: 03/31/2005
• Issued: 10/21/2008
• Est. Priority Date: 04/21/2004
• Status: Expired due to Fees

• Current Assignee: Koninklijke Philips Electronics N.V. (Koninklijke Philips N.V.)
• Original Assignee: Koninklijke Philips Electronics N.V. (Koninklijke Philips N.V.)
• Inventors: Netsch, Thomas, Rose, Georg, Schomberg, Hermann
• Primary Examiner(s): Glick; Edward J
• Assistant Examiner(s): Taningco; Alexander H
• Application Number:
  US11/568,115
• Publication Number:
  US 20070195923A1
• Time in Patent Office: 1,300 Days
• Field of Search: 378/4, 378/901
• US Class Current: 378/4
• CPC Class Codes: A61B 6/027 characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
• 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/432 Truncation
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

A 3D image of a region of an object is computed from truncated cone beam projection data acquired with an x-ray device and a prior CT image representing a larger region of the object. The truncated projection data are extrapolated to derive pseudoprojection data associated with projection directions outside the detector, and an intermediate CT image is reconstructed based on the truncated projection data completed with the pseudoprojection data. The prior CT image is then registered with the intermediate CT image. Forward projection data associated with projection directions outside the detector are computed from the truncated projection data and the registered prior CT image. The 3D image is finally reconstructed based on the truncated projection data completed with the forward projection data.

Claims

1. An apparatus for reconstructing a 3D image of at least a first region of an object from truncated cone beam projection data acquired by means of an x-ray device having a source-detector assembly moved along an acquisition trajectory such that a cone beam between the source and the detector encompasses said first region, the apparatus comprising:

extrapolation means for processing the truncated projection data to derive pseudo-projection data associated with projection directions outside the detector;

first 3D reconstruction means for reconstructing an intermediate CT image covering at least the first region of the object based on the truncated projection data completed with the pseudo-projection data;

registration means for computing registration data representing a spatial transformation selected to register a prior CT image, representing a second region of the object including the first region and larger than the first region, with the intermediate CT image;

forward projection means for processing the truncated projection data and the registered prior CT image to derive forward projection data associated with projection directions outside the detector; and

second 3D reconstruction means for synthesizing a reconstructed CT image covering at least the first region of the object based on the truncated projection data completed with the forward projection data.

2. The apparatus as claimed in claim 1, wherein the registration means comprise means for determining a template of voxels of the intermediate CT image which are representative of contours present in the intermediate CT image, and means for determining a spatial transformation of the prior CT image which maximizes a local correlation measurement, at the voxels of the template, between the intermediate CT image and the transformed prior CT image.

3. The apparatus as claimed in claim 2, wherein the means for determining the template comprise means for determining a set of voxels of the intermediate CT image having respective x-ray attenuation values above a threshold, and means for selecting the voids of the template whereby each void of the template has a respective local neighborhood including at least a predefined number of voxels of said set.

4. The apparatus as claimed in claim 3, wherein the means for selecting the voxels of the template are arranged to select voids of the intermediate CT image having maximum edge strength values among the voids having local neighborhoods including at least the predefined number of voids of said set.

5. The apparatus as claimed in claims 2, wherein the maximized local correlation measurement is of the form

6. The apparatus as claimed in claim 1, wherein the forward projection means comprise first computing means for obtaining first forward projection data with a coarser resolution than the truncated cone beam projection data by means of the registered prior CT image, and second computing means for interpolating the first forward projection data to provide the forward projection data with the same spatial resolution as the truncated projection data.

7. The apparatus as claimed in any one of the preceding claims, wherein the forward projection means comprise linear transformation means to fit the forward projection data with the truncated cone beam projection data in a region located along the border of the detector.

8. A computer readable storage medium containing instructions which, when executed by a computer for reconstructing a 3D image of at least a first region of an object from truncated cone beam projection data acquired by means of an x-ray device having a source-detector assembly moved along an acquisition trajectory such that a cone beam between the source and the detector encompasses said first region, cause the computer to perform the step of:

extrapolating the truncated projection data to derive pseudo-projection data associated with projection directions outside the detector;

reconstructing an intermediate 3D CT image covering at least the first region of the object based on the truncated projection data completed with the pseudo-projection data;

registering a prior CT image, representing a second region of the object including the first region and larger than the first region, with the intermediate CT image;

processing the truncated projection data and the registered prior CT image to derive forward projection data associated with projection directions outside the detector; and

synthesizing a reconstructed CT image covering at least the first region of the object based on the truncated projection data completed with the forward projection data.

9. The computer readable storage medium as claimed in claim 8, wherein the instructions for registering the prior CT image with the intermediate CT image comprise instructions for determining a template of voxels of the intermediate CT image which are representative of contours present in the intermediate CT image, and instructions for determining a spatial transformation of the prior CT image which maximizes a local correlation measurement, at the voxels of the template, between the intermediate CT image and the transformed prior CT image.

10. The computer readable storage medium as claimed in claim 9, wherein the instructions for determining the template comprise instructions for determining a set of voxels of the intermediate CT image having respective x-ray attenuation values above a threshold, and instructions for selecting the voxels of the template whereby each voxel of the template has a respective local neighborhood including at least a predefined number of voxels of said set.

11. The computer readable storage medium as claimed in claim 10, wherein the instructions for selecting the voxels of the template select voxels of the intermediate CT image having maximum edge strength values among the voxels having local neighborhoods including at least the predefined number of voxels of said set.

12. The computer readable storage medium as claimed in claim 8, wherein the instructions to derive forward projection data comprise instructions for obtaining first forward projection data with a coarser resolution than the truncated cone beam projection data by means of the registered prior CT image, and instructions for interpolating the first forward projection data to provide the forward projection data with the same spatial resolution as the truncated projection data.

13. The computer readable storage medium as claimed in claim 8, wherein the instructions to derive forward projection data comprise instructions to fit the forward projection data with the truncated cone beam projection data in a region located along the border of the detector.