Exact region of interest cone beam imaging without circle scans

• US 6,333,960 B1
• Filed: 06/30/1999
• Issued: 12/25/2001
• Est. Priority Date: 03/22/1999
• Status: Expired due to Term

• Current Assignee: Siemens Corp. (Siemens AG)
• Original Assignee: Siemens Corporate Research Incorporated (Siemens AG)
• Inventors: Tam, Kwok
• Primary Examiner(s): BRUCE, DAVID VERNON
• Application Number:
  US09/343,770
• Time in Patent Office: 909 Days
• Field of Search: 378/4, 378/15, 378/901
• US Class Current: 378/15
• CPC Class Codes: A61B 6/027 characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
• CPC Class Codes: A61B 6/032 Transmission computed tomography [CT]
• CPC Class Codes: A61B 6/4085 Cone-beams
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

A scanning and data acquisition method and apparatus for three dimensional (3D) computerized tomography (CT) imaging of an ROI in an object, wherein data is calculated on a plurality of line segments L formed in cone beam projection data acquired on a detector at a plurality of source positions. The endpoints of the line segments L formed in the data acquired at each of the source positions is determined by a data combination mask having its bounds formed by cone beam projection onto the plane of the detector of portions of the source scan path that are above and below the source position that acquired the cone beam data in which the data for the line segments L are being calculated. When processing cone beam data near the top or bottom edges of the ROI, a first group of the line segments L has both of its endpoints determined by the boundaries of the data combination mask, and a second group of the line segments L has one of its endpoints determined by a boundary of the data combination mask, and another of its endpoints determined by the x-axis of the mask, so as to avoid data corruption from x-rays that passed outside the ROI. The line integral data calculated for both groups of the line segments L formed in the mask are then processed using an image reconstruction technique for providing a contribution to a 3D image reconstruction of the ROI in the object.

Claims

1. A method for performing three dimensional computerized tomographic imaging of a region-of-interest (ROI) in an object using a cone beam source of radiation energy, comprising the steps of:

defining a source scanning trajectory along a longitudinal axis as a scan path that encircles the ROI in the object and is traversed by the cone beam source;

using the cone beam source, fixed relative to an area detector with both source and detector movably positioned relative to the object, for applying radiation energy towards the object from a plurality of source positions along the scan path as said source traverses the scan path, said applying causing said area detector to acquire cone beam projection data corresponding to respective portions of the object at each of said source positions;

masking the cone beam projection data acquired at each of said source positions with a mask having upper and lower boundaries formed by cone beam projections onto the plane of the detector of portions of the source scan path that are above and below, respectively, the source position that acquired the data being masked;

calculating data for each of a plurality of line segments L formed in the masked cone beam projection data acquired at each of said source positions, said mask being used during said calculating step to divide said plurality of line segments L into separate groups, each group having a different technique for limiting the length of the line segments L therein; and

reconstruction processing the data calculated for said line segments L having their length determined by the mask as defined above, for reconstructing a 3D image of the ROI in the object;

wherein the calculating step includes a step for determining a point (C₀) on an axial line crossing the mask which intersects a projection of the source position that acquired the data, the point C_o is determined by;

a) drawing a first line (606) which is tangent to the scan path (604) and perpendicular to its longitudinal axis;

b) drawing a second line which is parallel to the first line, and which passes through the source position (S_i) which acquired the data being masked; and

c) defining as point C_o, that point where the second line intersects the detector (612).

2. The method of claim 1, wherein during said calculating step, those of said line segments L formed in the masked cone beam data which cross said axial line to one side of C₀ have their length limited using one technique, and those of said line segments L which cross said axial line to the other side of C₀ have their length limited using another technique.

3. The method of claim 2, wherein said one technique comprises determining both endpoints of the line segments L using the upper and lower boundaries of the mask, and said other technique determines one of the endpoints of the line segments L using one of the upper or lower boundaries of the mask, and its other end point is determined using the axial line of the mask.

4. The method of claim 1, wherein:

said calculating step calculates line integral derivative data for each of said line segments L, and

said reconstruction processing step comprises accumulating the derivative data in Radon space, and inversion processing of the accumulated Radon data for reconstructing an image of the ROI in the object.

5. The method of claim 1, wherein:

said calculating step calculates line integral derivative data for each of said line segments L, and

said reconstruction processing step comprises 2D backprojecting the line integral derivatives calculated for a plurality of said line segments L onto a 2D space corresponding to the plane of the detector, and 3D backprojecting the 2D backprojection results into a 3D space, thereby reconstructing a 3D image of the ROI in the object.

6. The method of claim 5, wherein after completion of the 2D backprojecting step, a 1D d/dt filtering step is performed on the results in said 2D space along lines in the direction of a projection of the scan path.

7. The method of claim 6, wherein the 3D backprojecting step comprises performing a weighted 3D backprojection of the 1D d/dt filtering results at discrete points along said lines onto a plurality of sample points in a 3D object volume.