US6275561B1 - Computer tomagraphy method with helicoidal scanning of an examination area - Google Patents

Computer tomagraphy method with helicoidal scanning of an examination area Download PDF

Info

Publication number
US6275561B1
US6275561B1 US09/380,972 US38097299A US6275561B1 US 6275561 B1 US6275561 B1 US 6275561B1 US 38097299 A US38097299 A US 38097299A US 6275561 B1 US6275561 B1 US 6275561B1
Authority
US
United States
Prior art keywords
rotation
axis
examination zone
data
radiation source
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
US09/380,972
Inventor
Per-Erik Danielsson
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
US Philips Corp
Original Assignee
US Philips Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from PCT/SE1998/000029 external-priority patent/WO1998030980A1/en
Application filed by US Philips Corp filed Critical US Philips Corp
Priority to US09/380,972 priority Critical patent/US6275561B1/en
Priority claimed from PCT/IB1999/000027 external-priority patent/WO1999036885A1/en
Assigned to U.S. PHILIPS CORPORATION reassignment U.S. PHILIPS CORPORATION ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: DANIELSSON, PER-ERIK
Application granted granted Critical
Publication of US6275561B1 publication Critical patent/US6275561B1/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/027Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral

Definitions

  • the invention relates to a computer tomography method which involves helical scanning of an examination zone by means of a scanning unit which includes a radiation source and a detector unit, an object present in the examination zone and the scanning unit simultaneously rotating about an axis of rotation relative to one another and performing a motion parallel to the direction of the axis of rotation, resulting in a relative motion in the form of a helix, and also involves a reconstruction of the spatial distribution of the absorption within the examination zone from the measuring data acquired by the detector unit.
  • the invention also relates to a computer tomography apparatus for performing a method of this kind.
  • a method and a computer tomography apparatus of the kind set forth are known from DE-OS 195 45 778 (Tam).
  • This known method enables the scanning of an extended examination zone in the direction of the axis of rotation by means of a cone beam and the reconstruction of the absorption distribution in the examination zone also in case the object present therein, for example a patient, is longer than the part of the examination zone for which data has been acquired.
  • each group comprising a plurality of planes which extend parallel to the axis of rotation and each of which contains a respective fan beam
  • the invention utilizes exclusively measuring data acquired while the radiation source irradiates the points within the examination zone from an angular range of exactly 180° as seen from the relevant point itself; the rays associated with the measuring data then pass exactly through the area between two adjacent turns of the helix. On the one hand this angular range suffices to enable exact reconstruction while on the other hand it avoids the use of redundant measuring data.
  • the type of rebinning performed in relation to these measuring data is essential (rebinning is to be understood to mean the resorting of the measuring data from the sequence yielded by the acquisition as well as the re-interpolation of the measuring data on a different grid).
  • Rebinning is performed with data having been acquired along rays which emanate from the same source point and form a fan beam extending in a plane parallel to the rotation axis.
  • the subsequent processing steps i.e. the preferably one-dimensional filtering and the reconstruction, are thus significantly facilitated.
  • each group contains only mutually parallel planes so that the further processing is significantly facilitated. Notably an outstanding image quality is thus obtained.
  • the rebinning is performed on a virtual detector which extends perpendicularly to the planes associated with the relevant group and has a rectangular surface. The re-interpolation on an equidistant measuring point grid as required for the subsequent processing steps is thus significantly facilitated.
  • the absorption distribution could in principle be reconstructed from the filtered data of different groups by means of so-called generalized projections as described in the document published by Schaller et al. in SPIE, Vol. 3032, 32, pp. 213 to 224.
  • a preferred type of reconstruction is realized by back projection of the filtered data in conformity with claim 3 .
  • the filtering operation could, also be performed, for example, by subjecting the data produced by the rebinning operation to a convolution with an appropriate filter kernel.
  • the filtering operation defined in claim 4 requires less calculation time.
  • Claim 6 discloses a computer tomography apparatus for performing the method according to the invention and claim 7 defines an advantageous embodiment thereof.
  • the shape of the collimator arrangement and/or the shape of the detector unit in this embodiment ensure that each point within the examination zone “sees” the radiation source, upon emerging from the radiation beam generated thereby, at an angle which has been shifted through exactly 180°( ⁇ ) with respect to the angle upon entering of the radiation beam.
  • the advantage of this step resides in the fact that all measuring data required for an exact reconstruction is measured (and no other data). Thus, removal or co-weighting of redundant measuring data is not necessary.
  • FIG. 1 is a diagrammatic representation of a computer tomography apparatus according to the invention
  • FIG. 2 shows the helical scanning path described relative to one another by the scanning unit and an object present in the examination zone
  • FIG. 3 is a perspective view of the scanning unit and the examination zone
  • FIG. 4 is a plan view of FIG. 3,
  • FIG. 5 is a developed view of the detector unit
  • FIG. 6 illustrates various possibilities for the arrangement of the detector unit
  • FIG. 7 is a developed view of a detector unit situated on a cylinder around the radiation source
  • FIG. 8 is a perspective view of fans associated with a group of rebinned data
  • FIG. 9 shows a flow chart illustrating the processing of the measuring data
  • FIG. 10 is a side elevation of the detector units with some rays
  • FIG. 11 is a plan view of the arrangement shown in FIG. 10, and
  • FIG. 12 shows the same plan view as FIG. 11, be it after partial rebinning of measuring data.
  • the computer tomography apparatus shown in FIG. 1 includes a gantry 1 which is capable of rotation about an axis of rotation 14 which extends parallel to the z direction. To this end, the gantry is driven at a preferably constant angular speed by a motor 2 .
  • a radiation source S for example an X-ray tube, which is provided with a collimator arrangement 3 which forms a conical radiation beam 4 from the radiation produced by the radiation source S, i.e. a radiation beam having a finite dimension in the direction of the z axis as well as in the direction perpendicular thereto (i.e. in the x-y plane of the Cartesian co-ordinate system shown in FIG. 1 ).
  • the radiation beam irradiates an examination zone 13 or (not shown) an object, for example a patient arranged on a patient table.
  • the X-ray beam 4 is incident on a two-dimensional detector unit 16 which is attached to the gantry 1 and comprises a plurality of detector rows, each of which comprises a plurality of detector elements.
  • Each detector element detects a ray from the radiation beam 4 in each radiation source position.
  • the detector unit 16 may be arranged on an arc of circle which coincides with the circular path of the radiation source S during the rotation.
  • the angle of aperture K max of the radiation beam 4 determines the diameter of the examination zone 13 which is concentric with the axis of rotation 14 and in which the object to be examined must be present during the acquisition of the measuring values by means of the detector unit.
  • the examination zone 13 or a patient who is arranged, for example on a patient table present therein, can be shifted parallel to the direction of the axis of rotation 14 or the z axis by means of a motor 5 .
  • the measuring data then acquired by the detector unit 16 is applied to an image processing computer 10 which derives therefrom the distribution of the absorption of the emitted radiation in the part of the examination zone 13 covered by the cone beam 4 and reproduces it, for example on a monitor 11 .
  • the two motors 2 and 5 , the image processing computer 10 , the radiation source S and the transfer of the measuring data from the detector unit 16 to the image processing computer 10 are controlled by means of an appropriate control unit 7 .
  • the control unit 7 controls the motors 2 and 5 in such a manner that the ratio of the speed v of the examination zone 13 to the angular speed of the gantry 1 is constant.
  • the radiation source S and the examination zone then move along a helical path relative to one another. It is then irrelevant in principle whether the scanning unit or the examination zone performs the rotation or translation; only the relative motion is of importance.
  • the radiation source S (and the detector unit 16 connected thereto via the gantry) move along the helical path 17 in FIG. 2, whereas the examination zone 13 (or the object present therein) which is not shown in FIG. 2 is stationary.
  • the conical radiation beam 4 emitted by the source S is incident on the detector unit 16 arranged to the other side of the examination zone.
  • the conical beam 4 is depicted as if it contained a respective fan beam in a plurality of planes parallel to the axis of rotation 14 (or to the z direction). All of said flat fan beams emanate from the relevant position of the radiation source S or intersect one another in this position.
  • FIG. 3 is a perspective view of the radiation source S, symbolized by the point wherefrom the radiation emanates, the detector unit 16 and the cylindrical examination zone 13 having the radius r. Also shown is a cylinder 12 which concentrically encloses the examination zone 13 and has a radius R; the helical scanning path ( 17 in FIG. 2) is situated on this cylinder which, therefore, will also be referred to as the helix cylinder hereinafter.
  • the detector unit 16 consisting of a mosaic of detector elements as indicated in FIG.
  • FIG. 3 also shows two rays 18 a and 18 b which are incident on the lower edge and the upper edge, respectively, of the detector unit and pass through the boundary of the examination zone 13 at the points Q 1 and Q 2 .
  • the point Q on the ray between these points is the point on the ray which is situated at the shortest distance from the axis of rotation 14 .
  • the rays 18 a and 18 b strike the points Q 1 -Q-Q 2 upon entering the radiation beam and upon leaving the radiation beam. Thus, they are detected by the detector unit at different instants.
  • FIG. 4 is a plan view of the arrangement of FIG. 3, i.e. a view taken parallel to the z axis or the axis of rotation 14 .
  • the radius r may also be larger than R/2 (but always smaller than R) or smaller as is shown in FIG. 1 .
  • the projection of the two rays 18 a and 18 b on the x-y plane encloses an angle K with respect to the x axis, i.e. a ray emanating from S and passing through the axis of rotation 14 .
  • the lower ray 18 a becomes the upper ray 18 b after the examination zone 13 has been rotated through the angle ⁇ +2K and shifted proportionally in the z direction.
  • the point Q 1 is situated between the radiation source S and Q 2 upon entering of the radiation cone, exactly the reverse situation occurs upon leaving.
  • each point in the examination zone belongs to one II-line and one II-line only. Therefore, each point is irradiated from an angular range 180° viewed from this-point itself. This suffices (and is necessary) to enable reconstruction of each point within the examination zone 13 which has entered and emerged from the cone beam.
  • the detector unit 16 thus delivers the measuring values required for exact reconstruction, but no redundant measuring values so that the reconstruction is considerably simplified.
  • FIG. 5 is a developed view of a detector unit 16 from the helix cylinder 12 in the plane of drawing.
  • FIG. 5 shows the curves for the points Q 1 , Q and Q 2 .
  • FIG. 5 also shows the curves for the points on another II-line with the points P 1 , P, P 2 whose projection in the z direction coincides with the upper ray 18 b in FIG. 4 upon entering the radiation beam (i.e. when their connecting line to the radiation source intersects the lower edge of the detector unit, see FIG. 3 ).
  • the development of the detector unit need not necessarily have the shape of a parallelogram as shown in FIG. 5 .
  • the measuring data from the detector elements which are situated outside the parallelogram on the detector of FIG. 5 can be ignored.
  • the detector unit can also describe a helical arc of circle 16 b or 16 a about the radiation source S which is tangent to the examination zone 13 or the helix cylinder.
  • the detector unit may also have a flat surface 16 c or be shaped arbitrarily.
  • FIG. 7 shows the developed views of the detector units 16 and 16 a in the plane of drawing which are tangent to the helix cylinder 12 on an arc of circle around the radiation source. It appears that the height of this development, i.e. the dimension in the direction of the z axis, varies and analogously do the dimensions of the detector rows, i.e. in conformity with the function h/cosK, where K is the angle enclosed by the projection of a ray in the x-y plane relative to the z axis (for example, see FIG. 4 ).
  • each measuring value from each detector element is first divided by a reference value and the resultant quotient is logarithmized.
  • the measuring data thus formed represents the line integral of the absorption of the radiation along a ray connecting the radiation source to the relevant detector element.
  • the subsequent processing steps serve to determine the spatial distribution of the absorption from these line integrals of the absorption.
  • a rebinning operation is performed. After the rebinning, or before that, the measuring data is weighted in that this data is mutliplied by a factor which corresponds to the cosine of the angle enclosed by the ray (for example, 18) associated with the measuring data relative to a plane intersecting the axis of rotation at right angles.
  • This weighting step can be omitted in those cases where the distance between two turns of the helix is small in comparison with their radius. Therefore, this step is not shown separately in FIG. 9 .
  • FIG. 10 is a side elevation of the arrangement shown in FIG. 3 .
  • FIG. 10 shows six rays of the conical beam 4 ; three rays 401 . . . 403 thereof strike the upper edge while three rays 411 . . . 413 strike the lower edge of the detector unit.
  • the rays 402 and 412 pass through the axis of rotation 14 whereas the rays 401 , 403 , and 411 , 413 pass the axis of rotation to the left and to the right, respectively.
  • Each time two of these rays constitute the edge rays of a fan beam whose rays are situated in a plane parallel to the z axis or axis of rotation 14 , for example the rays 401 and 411 , the rays 402 and 412 and the rays 403 and 413 .
  • FIG. 11 is a plan view of the arrangement shown in FIG. 10 . Because the rays 401 . . . 403 and 411 . . . 413 are situated in the same plane, extending perpendicularly to the plane of drawing of FIG. 11, they appear as a single ray in FIG. 11 .
  • the planes defined by these fan beams intersect one another in the radiation source position S I .
  • the plane containing the fan beams 401 , 411 encloses an angle +K 1 relative to the central plane containing the axis of rotation, whereas an angle ⁇ K 1 exists between the central plane and the planes of the fan beam 403 , 413 .
  • the fan beam 402 , 412 is identical to the central plane (being the x-z plane in FIG. 11 ).
  • FIG. 11 shows two further radiation source positions S I ⁇ K1 and S I+K1 to both sides of the central radiation source position S I , and a respective fan beam which emanates from this source position and passes through the axis of rotation 14 by way of the fan beams 420 and 430 , respectively.
  • the fan beams 420 and 403 , 413 extend parallel to one another, like the fan beams 401 , 411 and 430 .
  • the fan beams from different radiation source positions being situated in planes parallel to one another (and to the axis of rotation 14 ), or the measuring data associated with the rays constituting these fan beams, are combined so as to form a respective group.
  • fan beams and associated measuring values for which the sum of the angle characterizing the radiation source position (I or I ⁇ K 1 or I+K 1 ) and the fan angle K (being the angle enclosed by the plane of the fan beam relative to a plane containing the axis of rotation 14 ; these are, for example the angles ⁇ K 1 and +K 1 , respectively, in FIG. 11) is constant.
  • FIG. 12 shows such a set of fan beams, generated in different radiation source positions and situated in parallel planes extending through the examination zone 13 .
  • a virtual detector 72 is arranged perpendicularly to the planes in which said fan beams are situated and in the axis of rotation.
  • the dimensions of this virtual detector in the x-y plane correspond to the diameter of the examination zone 13 (2r).
  • the dimensions of the virtual detector in the z direction amount to h/2. This is because it can be demonstrated that the upper and lower edge rays of all fan beams coincide exactly with the upper edge and the lower edge (in the z direction) of this flat, virtual and exactly rectangular detector.
  • FIG. 8 is a perspective representation of the configuration, it being assumed that the radiation source moves along the helical path 17 .
  • the upper and the left-hand side of the virtual detector 72 are denoted by solid lines, whereas the lower side is represented by a dotted line.
  • the upper edge and the lower edge of the virtual detector 72 are extended by way of dashed lines 72 so as to form a rectangle 720 ; the helical path 17 extends from the lower right-hand corner point to the upper corner point thereof.
  • the fan beams emanating from different radiation source positions on the helical path and extending in mutually parallel planes, are represented as triangles whose upper and lower sides denote the upper and the lower edge ray of the fan beam and whose vertical side is intended to denote the position of the column of the detector unit which is struck by the relevant fan beam.
  • the upper and the lower edge rays enter the upper and the lower edge of the plane virtual detector 72 .
  • the fan beams situated to the right and to the left of the center are intercepted by detector columns which are situated higher and lower, respectively, in the z direction than the column which detects the fan beam which emanates from the central radiation source position and passes through the axis of rotation 14 .
  • each fan beam covers a column of the virtual detector.
  • all fan beams are assigned to a respective one of the groups so that, possibly after interpolation, each group contains fan beams situated in parallel planes which perpendicularly intersect the virtual detector 72 associated with the relevant group.
  • the second part of the rebinning operation is performed.
  • This is a further interpolation operation which is required because of the following.
  • the fan beams cover a vertical strip or a column in the virtual detector and the rays associated with a fan beam can be incident on the virtual detector at equidistant points, but the columns or vertical strips are situated at different distances from one another (for reasons of geometry they are situated nearer to one another at the outer side than at the inner side). Therefore, in the step 103 the data produced by the step 102 is interpolated in such a manner that the associated rays and the associated line integrals of the absorption are obtained for a regular Cartesian grid on the virtual detector. This completes rebinning with a parallel beam geometry on a rectangular detector surface with regularly distributed grid points, so that the subsequent processing is significantly facilitated.
  • the rebinning operation carried out in the steps 102 and 103 thus produces groups of measuring data and the associated rays, which would occur if a flat, rectangular detector in a plane containing the axis of rotation 14 (i.e. the virtual detector) were to acquire the measuring data of a radiation source extending along a segment of the line 17 and emitting fan beams extending perpendicularly to the detector plane and parallel to the axis of rotation 14 .
  • a one-dimensional filtering operation is performed in the step 104 .
  • a simple, one-dimensional, location-independent, preferably ramp-like filter is required in the line direction (the line direction extends vertically in the representation of FIG. 12, or in the longitudinal direction of the rectangle 720 , 72 in the perspective view of FIG. 8, and hence perpendicularly to the axis of ration 14 ).
  • the filtering operation can in principle be carried out by subjecting the data resulting from the rebinning operation to a convolution with a suitable one-dimensional filter kernel.
  • step 104 A simpler possibility, however, consists in subjecting the data produced by the rebinning operation to a Fourier transformation first in the step 104 .
  • step 105 the data thus transformed in the spatial frequency domain is subjected to a ramp-like filtering operation (in the line direction) during which the damping decreases linearly as the value of the frequency increases.
  • step 106 the data thus filtered in the spatial frequency domain is subjected to an inverse Fourier transformation, yielding filtered projection data.
  • the virtual detector on which the rebinning is based is flat and contains the axis of rotation.
  • the puncture points of the individual rays and the (virtual) detector then no longer describe lines perpendicular to the axis of rotation, but possibly curved lines.
  • the measuring data which is associated, within the fan beams associated with the same group, with the corresponding rays (for example, the each time upper rays, the upper rays but one the lower rays).
  • the filtered data of each group is backprojected into the spatial domain. The filtered data is thus assigned to the voxels in the examination zone (i.e.
  • Each voxel thus receives contributions from different groups of rays (or the associated data) which, in respect of this voxel, enclose an angle of 180° relative to one another.

Abstract

The invention relates to a CT method which involves helical scanning of an examination zone by means of a cone beam. Measuring data for the three-dimensional tomographic imaging of objects of arbitrary length are now acquired completely and without redundancy. The reconstruction requires only a one-dimensional filtering operation. Particularly simple processing steps are enabled by a rebinning operation during which fan beams which are situated in planes extending parallel to one another and to the axis of rotation are combined.

Description

This application is a continuation-in-part of application Ser. No. 09/228,219, filed Jan. 11, 1999.
BACKGROUND OF THE INVENTION
1. Field of the Invention
The invention relates to a computer tomography method which involves helical scanning of an examination zone by means of a scanning unit which includes a radiation source and a detector unit, an object present in the examination zone and the scanning unit simultaneously rotating about an axis of rotation relative to one another and performing a motion parallel to the direction of the axis of rotation, resulting in a relative motion in the form of a helix, and also involves a reconstruction of the spatial distribution of the absorption within the examination zone from the measuring data acquired by the detector unit. The invention also relates to a computer tomography apparatus for performing a method of this kind.
2. Description of the Related Art
A method and a computer tomography apparatus of the kind set forth are known from DE-OS 195 45 778 (Tam). This known method enables the scanning of an extended examination zone in the direction of the axis of rotation by means of a cone beam and the reconstruction of the absorption distribution in the examination zone also in case the object present therein, for example a patient, is longer than the part of the examination zone for which data has been acquired.
However, for this method it is necessary to define in advance the so-called region of interest (ROI) and at the beginning and at the end of the ROI (in relation to the direction of rotation) an additional scan of the examination zone must be performed along a circular path which extends perpendicularly to the axis of rotation. The reconstruction of the absorption distribution may commence only after the scanning of the examination zone has been completed. The change over from a circular scan to a helical scan and back to a circular scan of the examination zone necessitates an abrupt acceleration or deceleration of either the scanning unit or the object in the examination zone; this could cause unsharpness. It is a further drawback that the scanning region must be defined in advance.
SUMMARY OF THE INVENTION
Therefore, it is an object of the present invention to provide a method of the kind set forth which does not require additional circular scanning motions and enables the reconstruction already during the acquisition of the measuring data. This object is achieved by a method of the kind set forth in that it includes the following steps:
a) using exclusively the measuring data for the reconstruction which are obtained from rays which pass exactly through the area between two adjacent turns of the helix, the points within the examination zone thereby being irradiated over an angular range of exactly 180 degrees—as seen from the point itself,
b) rebinning the measuring data and the associated rays so as to form a number of groups, each group comprising a plurality of planes which extend parallel to the axis of rotation and each of which contains a respective fan beam,
c) filtering the data of each group formed by the rebinning
d) reconstructing of the spatial distribution of the absorption from the filtered data of different groups.
The invention utilizes exclusively measuring data acquired while the radiation source irradiates the points within the examination zone from an angular range of exactly 180° as seen from the relevant point itself; the rays associated with the measuring data then pass exactly through the area between two adjacent turns of the helix. On the one hand this angular range suffices to enable exact reconstruction while on the other hand it avoids the use of redundant measuring data. The type of rebinning performed in relation to these measuring data is essential (rebinning is to be understood to mean the resorting of the measuring data from the sequence yielded by the acquisition as well as the re-interpolation of the measuring data on a different grid). Rebinning is performed with data having been acquired along rays which emanate from the same source point and form a fan beam extending in a plane parallel to the rotation axis. The subsequent processing steps, i.e. the preferably one-dimensional filtering and the reconstruction, are thus significantly facilitated.
In principle there are various possibilities for grouping the data emanating from the above-mentioned fan beams. For most of these groups, however, it will then be necessary to weight the measuring data with suitably chosen weighting factors (which may be dependent on the type of detector and on the type of rebinning). This necessity is eliminated, however, in the preferred version which is disclosed in claim 2 and in which each group contains only mutually parallel planes so that the further processing is significantly facilitated. Notably an outstanding image quality is thus obtained. In the further version according to claim 5, for each group the rebinning is performed on a virtual detector which extends perpendicularly to the planes associated with the relevant group and has a rectangular surface. The re-interpolation on an equidistant measuring point grid as required for the subsequent processing steps is thus significantly facilitated.
The absorption distribution could in principle be reconstructed from the filtered data of different groups by means of so-called generalized projections as described in the document published by Schaller et al. in SPIE, Vol. 3032, 32, pp. 213 to 224. A preferred type of reconstruction, however, is realized by back projection of the filtered data in conformity with claim 3.
The filtering operation could, also be performed, for example, by subjecting the data produced by the rebinning operation to a convolution with an appropriate filter kernel. The filtering operation defined in claim 4, however, requires less calculation time.
Claim 6 discloses a computer tomography apparatus for performing the method according to the invention and claim 7 defines an advantageous embodiment thereof. The shape of the collimator arrangement and/or the shape of the detector unit in this embodiment ensure that each point within the examination zone “sees” the radiation source, upon emerging from the radiation beam generated thereby, at an angle which has been shifted through exactly 180°(Ξ) with respect to the angle upon entering of the radiation beam. The advantage of this step resides in the fact that all measuring data required for an exact reconstruction is measured (and no other data). Thus, removal or co-weighting of redundant measuring data is not necessary.
BRIEF DESCRIPTION OF THE DRAWING
The invention will be described in detail hereinafter with reference to the drawings. Therein:
FIG. 1 is a diagrammatic representation of a computer tomography apparatus according to the invention,
FIG. 2 shows the helical scanning path described relative to one another by the scanning unit and an object present in the examination zone,
FIG. 3 is a perspective view of the scanning unit and the examination zone,
FIG. 4 is a plan view of FIG. 3,
FIG. 5 is a developed view of the detector unit,
FIG. 6 illustrates various possibilities for the arrangement of the detector unit,
FIG. 7 is a developed view of a detector unit situated on a cylinder around the radiation source,
FIG. 8 is a perspective view of fans associated with a group of rebinned data,
FIG. 9 shows a flow chart illustrating the processing of the measuring data,
FIG. 10 is a side elevation of the detector units with some rays,
FIG. 11 is a plan view of the arrangement shown in FIG. 10, and
FIG. 12 shows the same plan view as FIG. 11, be it after partial rebinning of measuring data.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
The computer tomography apparatus shown in FIG. 1 includes a gantry 1 which is capable of rotation about an axis of rotation 14 which extends parallel to the z direction. To this end, the gantry is driven at a preferably constant angular speed by a motor 2. On the gantry there is mounted a radiation source S, for example an X-ray tube, which is provided with a collimator arrangement 3 which forms a conical radiation beam 4 from the radiation produced by the radiation source S, i.e. a radiation beam having a finite dimension in the direction of the z axis as well as in the direction perpendicular thereto (i.e. in the x-y plane of the Cartesian co-ordinate system shown in FIG. 1).
The radiation beam irradiates an examination zone 13 or (not shown) an object, for example a patient arranged on a patient table. After having traversed the examination zone 13, the X-ray beam 4 is incident on a two-dimensional detector unit 16 which is attached to the gantry 1 and comprises a plurality of detector rows, each of which comprises a plurality of detector elements. Each detector element detects a ray from the radiation beam 4 in each radiation source position. The detector unit 16 may be arranged on an arc of circle which coincides with the circular path of the radiation source S during the rotation.
The angle of aperture Kmax of the radiation beam 4 (the angle of aperture is to be understood to mean the angle enclosed by a ray of the beam 4 which is situated at the edge in the x-y plane relative to the ray which intersects the axis of rotation 14 at right angles) then determines the diameter of the examination zone 13 which is concentric with the axis of rotation 14 and in which the object to be examined must be present during the acquisition of the measuring values by means of the detector unit. The examination zone 13, or a patient who is arranged, for example on a patient table present therein, can be shifted parallel to the direction of the axis of rotation 14 or the z axis by means of a motor 5. The measuring data then acquired by the detector unit 16 is applied to an image processing computer 10 which derives therefrom the distribution of the absorption of the emitted radiation in the part of the examination zone 13 covered by the cone beam 4 and reproduces it, for example on a monitor 11. The two motors 2 and 5, the image processing computer 10, the radiation source S and the transfer of the measuring data from the detector unit 16 to the image processing computer 10 are controlled by means of an appropriate control unit 7.
The control unit 7 controls the motors 2 and 5 in such a manner that the ratio of the speed v of the examination zone 13 to the angular speed of the gantry 1 is constant. The radiation source S and the examination zone then move along a helical path relative to one another. It is then irrelevant in principle whether the scanning unit or the examination zone performs the rotation or translation; only the relative motion is of importance.
Therefore, it has been assumed in FIG. 2 that the radiation source S (and the detector unit 16 connected thereto via the gantry) move along the helical path 17 in FIG. 2, whereas the examination zone 13 (or the object present therein) which is not shown in FIG. 2 is stationary. The conical radiation beam 4 emitted by the source S is incident on the detector unit 16 arranged to the other side of the examination zone. The conical beam 4 is depicted as if it contained a respective fan beam in a plurality of planes parallel to the axis of rotation 14 (or to the z direction). All of said flat fan beams emanate from the relevant position of the radiation source S or intersect one another in this position.
FIG. 3 is a perspective view of the radiation source S, symbolized by the point wherefrom the radiation emanates, the detector unit 16 and the cylindrical examination zone 13 having the radius r. Also shown is a cylinder 12 which concentrically encloses the examination zone 13 and has a radius R; the helical scanning path (17 in FIG. 2) is situated on this cylinder which, therefore, will also be referred to as the helix cylinder hereinafter. The detector unit 16, consisting of a mosaic of detector elements as indicated in FIG. 1 which may be arranged in columns (parallel to the axis of rotation 14) and rows, is situated on the circumference of the helix cylinder 12 between two successive turns of the helix 12; its dimensions in the z direction thus correspond to the pitch h of the turns of the helix. In this case it is assumed that the radiation source S and the detector unit 16 are stationary whereas the examination zone 13 with the object present therein is shifted in the direction of the axis of rotation 14 and at the same time rotates counterclockwise about this axis of rotation 14; the cylinder 13, therefore, would be shifted out of the helix cylinder 12 in the upwards direction.
FIG. 3 also shows two rays 18 a and 18 b which are incident on the lower edge and the upper edge, respectively, of the detector unit and pass through the boundary of the examination zone 13 at the points Q1 and Q2. The point Q on the ray between these points is the point on the ray which is situated at the shortest distance from the axis of rotation 14. The rays 18 a and 18 b strike the points Q1-Q-Q2 upon entering the radiation beam and upon leaving the radiation beam. Thus, they are detected by the detector unit at different instants.
FIG. 4 is a plan view of the arrangement of FIG. 3, i.e. a view taken parallel to the z axis or the axis of rotation 14. The angle of aperture Kmax of the cone beam in the direction of the x-y plane then amounts to 45°, meaning that the radius r=R/2. However, the radius r may also be larger than R/2 (but always smaller than R) or smaller as is shown in FIG. 1.
The projection of the two rays 18 a and 18 b on the x-y plane encloses an angle K with respect to the x axis, i.e. a ray emanating from S and passing through the axis of rotation 14. The lower ray 18 a becomes the upper ray 18 b after the examination zone 13 has been rotated through the angle Ξ+2K and shifted proportionally in the z direction. However, whereas the point Q1 is situated between the radiation source S and Q2 upon entering of the radiation cone, exactly the reverse situation occurs upon leaving.
This means that the points Q1 and Q2, and all other points on the line through Q1 and Q2, have been irradiated or projected onto the detector unit 16 in an angular range of exactly 180°, viewed from the relevant point. A line with points simultaneously entering and simultaneously leaving the examination zone 13 in this manner is called a II-line. The line 18 in FIG. 3 and FIG. 4 is such a II-line and it will be evident that a II-line is any line intercounting two points on the same turn of the helical scanning path 17. It can be demonstrated that each point in the examination zone belongs to one II-line and one II-line only. Therefore, each point is irradiated from an angular range 180° viewed from this-point itself. This suffices (and is necessary) to enable reconstruction of each point within the examination zone 13 which has entered and emerged from the cone beam. The detector unit 16 thus delivers the measuring values required for exact reconstruction, but no redundant measuring values so that the reconstruction is considerably simplified.
FIG. 5 is a developed view of a detector unit 16 from the helix cylinder 12 in the plane of drawing. The development is formed as a parallelogram having sides extending parallel to the z direction; the upper and the lower side enclose an angle M relative to the axis of rotation in conformity with the slope of the helix, which angle M can be calculated from the relation tanM=h/2ΞR. In this respect it is assumed that the speeds of translation and rotation (or the angular speed) are constant and that a complete rotation about the axis of rotation 14 takes place in the same period of time in which a displacement over the distance h takes place in the z direction.
During the passage of the cone beam by a point, its projection on the detector unit 16 continuously changes position. Starting at the lower edge (or at the lower detector row) of the detector unit, it describes a curve on the detector unit which terminates at the upper edge. FIG. 5 shows the curves for the points Q1, Q and Q2. FIG. 5 also shows the curves for the points on another II-line with the points P1, P, P2 whose projection in the z direction coincides with the upper ray 18 b in FIG. 4 upon entering the radiation beam (i.e. when their connecting line to the radiation source intersects the lower edge of the detector unit, see FIG. 3). These points are situated nearer to the radiation source S than the points Q1-Q-Q2 during their passage through the cone beam, and they pass the cone beam during a rotation through the angle Ξ-2K about the axis of rotation 14. During their passage through the cone beam 4 these points nevertheless also “see” the radiation source S from an angle of exactly 180°. The greater the distance between two points on the same II-line, the greater the difference will be between the two curves described by these two points on the detector.
The development of the detector unit need not necessarily have the shape of a parallelogram as shown in FIG. 5. Use could also be made of a larger, for example rectangular detector if the collimator 3 (FIG. 1) limits the conical radiation beam of the radiation source S in such a manner that the development of the area of the detector unit struck by the radiation beam has exactly the shape shown in FIG. 5. Instead of this step, or in combination therewith, the measuring data from the detector elements which are situated outside the parallelogram on the detector of FIG. 5 can be ignored.
It is not necessary either for the detector elements to be situated on the circumference of the helix cylinder 12 (FIG. 3). As is indicated in FIG. 6, diagrammatically representing a parallel projection of the arrangement of FIG. 3 in the direction of the z axis, the detector unit can also describe a helical arc of circle 16 b or 16 a about the radiation source S which is tangent to the examination zone 13 or the helix cylinder. The detector unit may also have a flat surface 16 c or be shaped arbitrarily. For all these versions it is essential only that the edges of the detector unit (or the respective area of the detector unit on which the radiation beam 4 is incident) coincide with the central projection of two segments of a turn of the helical scanning path 17, or that each point “sees” the radiation source in an angular range of exactly 180° during its passage.
FIG. 7 shows the developed views of the detector units 16 and 16 a in the plane of drawing which are tangent to the helix cylinder 12 on an arc of circle around the radiation source. It appears that the height of this development, i.e. the dimension in the direction of the z axis, varies and analogously do the dimensions of the detector rows, i.e. in conformity with the function h/cosK, where K is the angle enclosed by the projection of a ray in the x-y plane relative to the z axis (for example, see FIG. 4).
The further processing by the image processing computer 10 of the data acquired by the multi-row detector unit 16 will be described in detail hereinafter with reference to the flow chart shown in FIG. 9. After the initialization (block 100), each measuring value from each detector element is first divided by a reference value and the resultant quotient is logarithmized. The measuring data thus formed represents the line integral of the absorption of the radiation along a ray connecting the radiation source to the relevant detector element. The subsequent processing steps serve to determine the spatial distribution of the absorption from these line integrals of the absorption.
To this end, first a rebinning operation is performed. After the rebinning, or before that, the measuring data is weighted in that this data is mutliplied by a factor which corresponds to the cosine of the angle enclosed by the ray (for example, 18) associated with the measuring data relative to a plane intersecting the axis of rotation at right angles. This weighting step, however, can be omitted in those cases where the distance between two turns of the helix is small in comparison with their radius. Therefore, this step is not shown separately in FIG. 9.
During the rebinning operation, in a first step 102, there are formed groups of fan beams which are situated in planes parallel to one another and parallel to the axis of rotation 14, or groups of measuring data associated with the rays constituting these fan beams. This will first be described with reference to FIG. 10 which is a side elevation of the arrangement shown in FIG. 3. FIG. 10 shows six rays of the conical beam 4; three rays 401 . . . 403 thereof strike the upper edge while three rays 411 . . . 413 strike the lower edge of the detector unit. The rays 402 and 412 pass through the axis of rotation 14 whereas the rays 401, 403, and 411, 413 pass the axis of rotation to the left and to the right, respectively. Each time two of these rays constitute the edge rays of a fan beam whose rays are situated in a plane parallel to the z axis or axis of rotation 14, for example the rays 401 and 411, the rays 402 and 412 and the rays 403 and 413.
FIG. 11 is a plan view of the arrangement shown in FIG. 10. Because the rays 401 . . . 403 and 411 . . . 413 are situated in the same plane, extending perpendicularly to the plane of drawing of FIG. 11, they appear as a single ray in FIG. 11. The planes defined by these fan beams intersect one another in the radiation source position SI. The plane containing the fan beams 401, 411 encloses an angle +K1 relative to the central plane containing the axis of rotation, whereas an angle −K1 exists between the central plane and the planes of the fan beam 403, 413. The fan beam 402, 412 is identical to the central plane (being the x-z plane in FIG. 11).
FIG. 11 shows two further radiation source positions SI−K1 and SI+K1 to both sides of the central radiation source position SI, and a respective fan beam which emanates from this source position and passes through the axis of rotation 14 by way of the fan beams 420 and 430, respectively. It appears that the fan beams 420 and 403, 413 extend parallel to one another, like the fan beams 401, 411 and 430. According to the invention in the step 102 the fan beams from different radiation source positions being situated in planes parallel to one another (and to the axis of rotation 14), or the measuring data associated with the rays constituting these fan beams, are combined so as to form a respective group. Thus, combined in one group are those fan beams (and associated measuring values) for which the sum of the angle characterizing the radiation source position (I or I−K1 or I+K1) and the fan angle K (being the angle enclosed by the plane of the fan beam relative to a plane containing the axis of rotation 14; these are, for example the angles −K1 and +K1, respectively, in FIG. 11) is constant.
In practice only discrete values occur for the angle I of the radiation source positions or, due to the finite dimensions of the detector elements, of the fan angle K. These discrete values deviate from one another by increments dI and dK; it may be that dI≠dK. Because of this inequality, the sum of these two angles cannot always have exactly the same value for different radiation source positions; fan beams may occur for which the sum is slightly larger or slightly smaller than the angle I associated with the central radiation source position SI, i.e. the associated fan beams are not situated in parallel planes. In that case interpolation can be applied so as to derive measuring data corresponding to a fan beam which is situated in an exactly parallel plane or for which the sum is exactly α, using measuring data corresponding to these original fan beams which deviate slightly from the angle α.
FIG. 12 shows such a set of fan beams, generated in different radiation source positions and situated in parallel planes extending through the examination zone 13. A virtual detector 72 is arranged perpendicularly to the planes in which said fan beams are situated and in the axis of rotation. The dimensions of this virtual detector in the x-y plane correspond to the diameter of the examination zone 13 (2r). The dimensions of the virtual detector in the z direction amount to h/2. This is because it can be demonstrated that the upper and lower edge rays of all fan beams coincide exactly with the upper edge and the lower edge (in the z direction) of this flat, virtual and exactly rectangular detector.
FIG. 8 is a perspective representation of the configuration, it being assumed that the radiation source moves along the helical path 17. The upper and the left-hand side of the virtual detector 72 are denoted by solid lines, whereas the lower side is represented by a dotted line. The upper edge and the lower edge of the virtual detector 72 are extended by way of dashed lines 72 so as to form a rectangle 720; the helical path 17 extends from the lower right-hand corner point to the upper corner point thereof. The fan beams, emanating from different radiation source positions on the helical path and extending in mutually parallel planes, are represented as triangles whose upper and lower sides denote the upper and the lower edge ray of the fan beam and whose vertical side is intended to denote the position of the column of the detector unit which is struck by the relevant fan beam.
Even though the radiation source positions to the left and to the right of the center are situated higher and lower (measured in the direction of the axis of rotation), respectively, than the central radiation source position, the upper and the lower edge rays enter the upper and the lower edge of the plane virtual detector 72. This is due to the fact that the fan beams situated to the right and to the left of the center are intercepted by detector columns which are situated higher and lower, respectively, in the z direction than the column which detects the fan beam which emanates from the central radiation source position and passes through the axis of rotation 14.
It will be evident from the foregoing that each fan beam covers a column of the virtual detector. During the step 102 all fan beams (situated in planes parallel to the axis of rotation) are assigned to a respective one of the groups so that, possibly after interpolation, each group contains fan beams situated in parallel planes which perpendicularly intersect the virtual detector 72 associated with the relevant group.
After all measuring data and the associated fan beams have thus been detected for at least one group, in the step 103 the second part of the rebinning operation is performed. This is a further interpolation operation which is required because of the following. The fan beams cover a vertical strip or a column in the virtual detector and the rays associated with a fan beam can be incident on the virtual detector at equidistant points, but the columns or vertical strips are situated at different distances from one another (for reasons of geometry they are situated nearer to one another at the outer side than at the inner side). Therefore, in the step 103 the data produced by the step 102 is interpolated in such a manner that the associated rays and the associated line integrals of the absorption are obtained for a regular Cartesian grid on the virtual detector. This completes rebinning with a parallel beam geometry on a rectangular detector surface with regularly distributed grid points, so that the subsequent processing is significantly facilitated.
The rebinning operation carried out in the steps 102 and 103 thus produces groups of measuring data and the associated rays, which would occur if a flat, rectangular detector in a plane containing the axis of rotation 14 (i.e. the virtual detector) were to acquire the measuring data of a radiation source extending along a segment of the line 17 and emitting fan beams extending perpendicularly to the detector plane and parallel to the axis of rotation 14.
Subsequently, a one-dimensional filtering operation is performed in the step 104. For the given rebinning operation, merely a simple, one-dimensional, location-independent, preferably ramp-like filter is required in the line direction (the line direction extends vertically in the representation of FIG. 12, or in the longitudinal direction of the rectangle 720, 72 in the perspective view of FIG. 8, and hence perpendicularly to the axis of ration 14). The filtering operation can in principle be carried out by subjecting the data resulting from the rebinning operation to a convolution with a suitable one-dimensional filter kernel.
A simpler possibility, however, consists in subjecting the data produced by the rebinning operation to a Fourier transformation first in the step 104. In the step 105 the data thus transformed in the spatial frequency domain is subjected to a ramp-like filtering operation (in the line direction) during which the damping decreases linearly as the value of the frequency increases. In the step 106 the data thus filtered in the spatial frequency domain is subjected to an inverse Fourier transformation, yielding filtered projection data.
It is not essential that the virtual detector on which the rebinning is based is flat and contains the axis of rotation. The puncture points of the individual rays and the (virtual) detector then no longer describe lines perpendicular to the axis of rotation, but possibly curved lines. In that case it is necessary to use for the filtering operation the measuring data which is associated, within the fan beams associated with the same group, with the corresponding rays (for example, the each time upper rays, the upper rays but one the lower rays). During the next step 107 the filtered data of each group is backprojected into the spatial domain. The filtered data is thus assigned to the voxels in the examination zone (i.e. it is “spread” among these voxels) which were struck by the associated beam during the acquisition. Each voxel thus receives contributions from different groups of rays (or the associated data) which, in respect of this voxel, enclose an angle of 180° relative to one another.
This already yields a complete reconstruction of a part of the examination zone as soon as the measuring data produced by one pass of the scanning unit through an angle of 180°+2Kmax, relative to the examination zone, has been processed in conformity with the steps 101-107. The region thus reconstructed can be displayed on the monitor immediately thereafter; during this display measuring data can still be acquired and processed in conformity with the steps 101 to 107. The acquisition of the measuring data and the reconstruction of the absorption distribution in the examination zone can then be terminated at an arbitrary instant (block 108).
All references cited herein are incorporated herein by reference in their entirety and for all purposes to the same extent as if each individual publication or patent or patent application was specifically and individually indicated to be incorporated by reference in its entirety for all purposes.

Claims (7)

What is claimed is:
1. A computer tomography method which involves helical scanning of an examination zone by means of a scanning unit which includes a radiation source (S) and a detector unit (16), an object present in the examination zone (13) and the scanning unit simultaneously rotating about an axis of rotation relative to one another and performing a motion parallel to the direction of the axis of rotation, resulting in a relative motion in the form of a helix (17), and also involves a reconstruction of the spatial distribution of the absorption within the examination zone from the measuring data acquired by the detector unit, characterized in that it includes the following steps
a) using exclusively the measuring data for the reconstruction which are obtained from rays which pass exactly through the area between two adjacent turns of the helix, the points within the examination zone thereby being irradiated over an angular range of exactly 180 degrees—as seen from the point itself,
b) rebinning the measuring data and the associated rays so as to form a number of groups, each group comprising a plurality of planes which extend parallel to the axis of rotation and each of which contains a respective fan beam,
c) filtering the data of each group, formed by the rebinning,
d) reconstructing the spatial distribution of the absorption from the filtered data of different groups.
2. A computer tomography method as claimed in claim 1, characterized in that the measuring data is rebinned in such a manner that the planes containing a respective fan beam in each group extend parallel to one another and to the axis of rotation.
3. A computer tomography method as claimed in claim 1, characterized in that the reconstruction step includes the backprojection of the filtered data of a plurality of groups.
4. A computer tomography method as claimed in claim 2, characterized in that the filtering operation includes the following steps:
a) one-dimensional Fourier transformation of the data of each group in the direction perpendicular to the axis of rotation,
b) application of a ramp filter to the values yielded by the Fourier transformation,
c) inverse Fourier transformation of the filtered data.
5. A computer tomography method as claimed in claim 1, characterized in that the rebinning is performed on a respective flat virtual detector which extends perpendicularly to the planes of each group and contains the axis of rotation.
6. A computer tomography apparatus for carrying out the method claimed in claim 1, provided with a scanning unit which includes a radiation source and a detector unit connected thereto, a drive device for simultaneously rotating an object present in the examination zone and the scanning unit relative to one another about an axis of rotation and for causing them to perform a motion in the direction of the axis of rotation, and also with a reconstruction unit for reconstructing the spatial distribution of the absorption within the examination zone from the measuring data acquired by the detector unit, characterized in that it includes
a) means for using exclusively measuring data for the reconstruction which has been acquired during irradiation of the points within the examination zone from an angular range of 180°, the radiation source,
b) means for rebinning the measuring data so as to form a plurality of groups, each group containing a plurality of planes which extend parallel to the axis of rotation and in which a respective fan beam is situated,
c) means for the one-dimensional filtering of the data produced by the rebinning operation of each group in the direction perpendicular to the axis of rotation,
d) means for the reconstruction of the spatial distribution of the absorption from the filtered data of different groups.
7. A computer tomography apparatus as claimed in claim 6, characterized in that a collimator arrangement and/or detector unit is constructed in such a manner that all connecting lines between the radiation source and the angula range, mutually offset in the direction of rotation, of the area of the detector unit struck by the radiation beam, or the detector unit, intersect two segments of a turn of the helix which neighbor one another in the direction of rotation and along which the radiation source and the object are movable relative to one another.
US09/380,972 1998-01-13 1999-01-12 Computer tomagraphy method with helicoidal scanning of an examination area Expired - Fee Related US6275561B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US09/380,972 US6275561B1 (en) 1998-01-13 1999-01-12 Computer tomagraphy method with helicoidal scanning of an examination area

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
SE9800029 1998-01-13
PCT/SE1998/000029 WO1998030980A1 (en) 1997-01-14 1998-01-13 Technique and arrangement for tomographic imaging
US09/228,219 US6240157B1 (en) 1997-01-14 1999-01-11 Technique and arrangement for tomographic imaging
PCT/IB1999/000027 WO1999036885A1 (en) 1998-01-13 1999-01-12 Computer tomagraphy method with helicoidal scanning of an examination area
US09/380,972 US6275561B1 (en) 1998-01-13 1999-01-12 Computer tomagraphy method with helicoidal scanning of an examination area

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US09/228,219 Continuation-In-Part US6240157B1 (en) 1997-01-14 1999-01-11 Technique and arrangement for tomographic imaging

Publications (1)

Publication Number Publication Date
US6275561B1 true US6275561B1 (en) 2001-08-14

Family

ID=85557410

Family Applications (1)

Application Number Title Priority Date Filing Date
US09/380,972 Expired - Fee Related US6275561B1 (en) 1998-01-13 1999-01-12 Computer tomagraphy method with helicoidal scanning of an examination area

Country Status (1)

Country Link
US (1) US6275561B1 (en)

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030072406A1 (en) * 2001-09-26 2003-04-17 Xiaochun Yang Versatile cone-beam imaging apparatus and method
US20030202637A1 (en) * 2001-09-26 2003-10-30 Xiaochun Yang True 3D cone-beam imaging method and apparatus
WO2003094736A1 (en) * 2002-05-10 2003-11-20 Research Foundation Of The University Of Central Florida, Incorporated Ct image reconstruction method
US20040028173A1 (en) * 2002-08-06 2004-02-12 Koninklijke Philips Electronics N.V. Reconstruction method for tilted-gantry computed tomography
US20040131142A1 (en) * 2002-10-18 2004-07-08 Herbert Bruder Imaging method for a multi-slice spiral CT scan with 3D reconstruction, and a computed tomography unit for carrying out this method
US20040149923A1 (en) * 2001-05-11 2004-08-05 Beekman Frederik Johannes Method for obtaining a tomographic image, including apparatus
WO2004095066A1 (en) * 2003-04-22 2004-11-04 Philips Intellectual Property & Standards Gmbh Attenuation map generation from pet scans
US20040232348A1 (en) * 2001-05-11 2004-11-25 Beekman Frederik Johannes Method of obtaining a tomographic image
US20040240719A1 (en) * 2003-05-09 2004-12-02 Achim Gruebnau Method for producing images in spiral computed tomography, and a spiral CT unit
US20050123092A1 (en) * 2001-12-14 2005-06-09 Mistretta Charles A. Virtual spherical anode computed tomography
US20050243962A1 (en) * 2002-04-11 2005-11-03 Michael Grass Computed tomography method
FR2872659A1 (en) * 2004-06-30 2006-01-06 Gen Electric METHOD AND APPARATUS FOR DIRECT RECONSTRUCTION IN TOMOSYNTHESIS IMAGING
US20060018439A1 (en) * 2004-07-21 2006-01-26 Xiangyang Tang Method and apparatus for 3D reconstruction of images
US20060083344A1 (en) * 2004-10-15 2006-04-20 General Electric Company Methods and apparatus for reconstruction in helical cone beam volumetric CT
US20060140335A1 (en) * 2003-02-14 2006-06-29 Heuscher Dominic J System and method for helical cone-beam computed tomography with exact reconstruction
US20060233294A1 (en) * 2003-02-19 2006-10-19 Koninklijke Philips Electronics N.V. Computerized tomography method with helical relative movement and conical beam
US20060262895A1 (en) * 2005-05-20 2006-11-23 Ge Medical Systems Global Technology Company, Llc X-ray CT image reconstruction method and X-ray CT system
CN101882319A (en) * 2010-06-24 2010-11-10 西北工业大学 Cone beam CT (Computed Tomography) fast reconstructing method based on minimum three-dimensional convex hull
US10795457B2 (en) 2006-12-28 2020-10-06 D3D Technologies, Inc. Interactive 3D cursor
US11228753B1 (en) 2006-12-28 2022-01-18 Robert Edwin Douglas Method and apparatus for performing stereoscopic zooming on a head display unit
US11275242B1 (en) 2006-12-28 2022-03-15 Tipping Point Medical Images, Llc Method and apparatus for performing stereoscopic rotation of a volume on a head display unit
US11315307B1 (en) 2006-12-28 2022-04-26 Tipping Point Medical Images, Llc Method and apparatus for performing rotating viewpoints using a head display unit

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5323007A (en) * 1992-02-07 1994-06-21 Univ. Of Chicago Development Corp. Argonne National Laboratories Method of recovering tomographic signal elements in a projection profile or image by solving linear equations
US5744802A (en) * 1995-10-25 1998-04-28 Adac Laboratories Image generation from limited projections in positron emission tomography using multi-slice rebinning
US5828718A (en) * 1996-11-27 1998-10-27 Analogic Corporation Method and apparatus for helical computed tomography scanning with asymmetric detector system
US5848117A (en) * 1996-11-27 1998-12-08 Analogic Corporation Apparatus and method for computed tomography scanning using halfscan reconstruction with asymmetric detector system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5323007A (en) * 1992-02-07 1994-06-21 Univ. Of Chicago Development Corp. Argonne National Laboratories Method of recovering tomographic signal elements in a projection profile or image by solving linear equations
US5744802A (en) * 1995-10-25 1998-04-28 Adac Laboratories Image generation from limited projections in positron emission tomography using multi-slice rebinning
US5828718A (en) * 1996-11-27 1998-10-27 Analogic Corporation Method and apparatus for helical computed tomography scanning with asymmetric detector system
US5848117A (en) * 1996-11-27 1998-12-08 Analogic Corporation Apparatus and method for computed tomography scanning using halfscan reconstruction with asymmetric detector system

Cited By (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040149923A1 (en) * 2001-05-11 2004-08-05 Beekman Frederik Johannes Method for obtaining a tomographic image, including apparatus
US7145153B2 (en) * 2001-05-11 2006-12-05 Universitair Medisch Centrum Utrecht Method of obtaining a tomographic image
US20040232348A1 (en) * 2001-05-11 2004-11-25 Beekman Frederik Johannes Method of obtaining a tomographic image
US6944259B2 (en) 2001-09-26 2005-09-13 Massachusetts Institute Of Technology Versatile cone-beam imaging apparatus and method
US6865246B2 (en) 2001-09-26 2005-03-08 Massachusetts Institute Of Technology True 3D cone-beam imaging method and apparatus
US20030072406A1 (en) * 2001-09-26 2003-04-17 Xiaochun Yang Versatile cone-beam imaging apparatus and method
US20030202637A1 (en) * 2001-09-26 2003-10-30 Xiaochun Yang True 3D cone-beam imaging method and apparatus
US7333588B2 (en) * 2001-12-14 2008-02-19 Wisconsin Alumni Research Foundation Virtual spherical anode computed tomography
US20050123092A1 (en) * 2001-12-14 2005-06-09 Mistretta Charles A. Virtual spherical anode computed tomography
US20050243962A1 (en) * 2002-04-11 2005-11-03 Michael Grass Computed tomography method
US7142628B2 (en) * 2002-04-11 2006-11-28 Koninklijke Philips Electronics N.V. Computed tomography method
WO2003094736A1 (en) * 2002-05-10 2003-11-20 Research Foundation Of The University Of Central Florida, Incorporated Ct image reconstruction method
US6751283B2 (en) 2002-08-06 2004-06-15 Koninklijke Philips Electronics, N.V. Reconstruction method for tilted-gantry computed tomography
US20040028173A1 (en) * 2002-08-06 2004-02-12 Koninklijke Philips Electronics N.V. Reconstruction method for tilted-gantry computed tomography
US20040131142A1 (en) * 2002-10-18 2004-07-08 Herbert Bruder Imaging method for a multi-slice spiral CT scan with 3D reconstruction, and a computed tomography unit for carrying out this method
US7058156B2 (en) * 2002-10-18 2006-06-06 Siemens Aktiengesellschaft Imaging method for a multi-slice spiral CT scan with 3D reconstruction, and a computed tomography unit for carrying out this method
US20060140335A1 (en) * 2003-02-14 2006-06-29 Heuscher Dominic J System and method for helical cone-beam computed tomography with exact reconstruction
US7187747B2 (en) * 2003-02-19 2007-03-06 Koninklijke Philips Electronics N.V. Computerized tomography method with helical relative movement and conical beam
US20060233294A1 (en) * 2003-02-19 2006-10-19 Koninklijke Philips Electronics N.V. Computerized tomography method with helical relative movement and conical beam
US20060243914A1 (en) * 2003-04-22 2006-11-02 Koehler Thomas Attenuation map generation from pet scans
WO2004095066A1 (en) * 2003-04-22 2004-11-04 Philips Intellectual Property & Standards Gmbh Attenuation map generation from pet scans
US20040240719A1 (en) * 2003-05-09 2004-12-02 Achim Gruebnau Method for producing images in spiral computed tomography, and a spiral CT unit
FR2872659A1 (en) * 2004-06-30 2006-01-06 Gen Electric METHOD AND APPARATUS FOR DIRECT RECONSTRUCTION IN TOMOSYNTHESIS IMAGING
US20060018439A1 (en) * 2004-07-21 2006-01-26 Xiangyang Tang Method and apparatus for 3D reconstruction of images
US7583777B2 (en) * 2004-07-21 2009-09-01 General Electric Company Method and apparatus for 3D reconstruction of images
US20060083344A1 (en) * 2004-10-15 2006-04-20 General Electric Company Methods and apparatus for reconstruction in helical cone beam volumetric CT
US7539281B2 (en) * 2004-10-15 2009-05-26 General Electric Company Methods and apparatus for reconstruction in helical cone beam volumetric CT
US20060262895A1 (en) * 2005-05-20 2006-11-23 Ge Medical Systems Global Technology Company, Llc X-ray CT image reconstruction method and X-ray CT system
US7623615B2 (en) * 2005-05-20 2009-11-24 Ge Medical Systems Global Technology Company, Llc X-ray CT image reconstruction method and X-ray CT system
US11036311B2 (en) 2006-12-28 2021-06-15 D3D Technologies, Inc. Method and apparatus for 3D viewing of images on a head display unit
US10795457B2 (en) 2006-12-28 2020-10-06 D3D Technologies, Inc. Interactive 3D cursor
US10936090B2 (en) 2006-12-28 2021-03-02 D3D Technologies, Inc. Interactive 3D cursor for use in medical imaging
US10942586B1 (en) 2006-12-28 2021-03-09 D3D Technologies, Inc. Interactive 3D cursor for use in medical imaging
US11016579B2 (en) 2006-12-28 2021-05-25 D3D Technologies, Inc. Method and apparatus for 3D viewing of images on a head display unit
US11228753B1 (en) 2006-12-28 2022-01-18 Robert Edwin Douglas Method and apparatus for performing stereoscopic zooming on a head display unit
US11275242B1 (en) 2006-12-28 2022-03-15 Tipping Point Medical Images, Llc Method and apparatus for performing stereoscopic rotation of a volume on a head display unit
US11315307B1 (en) 2006-12-28 2022-04-26 Tipping Point Medical Images, Llc Method and apparatus for performing rotating viewpoints using a head display unit
US11520415B2 (en) 2006-12-28 2022-12-06 D3D Technologies, Inc. Interactive 3D cursor for use in medical imaging
CN101882319B (en) * 2010-06-24 2011-11-30 西北工业大学 Cone beam CT fast reconstructing method based on minimum three-dimensional convex hull
CN101882319A (en) * 2010-06-24 2010-11-10 西北工业大学 Cone beam CT (Computed Tomography) fast reconstructing method based on minimum three-dimensional convex hull

Similar Documents

Publication Publication Date Title
US6275561B1 (en) Computer tomagraphy method with helicoidal scanning of an examination area
US6269141B1 (en) Computer tomography apparatus with a conical radiation beam and a helical scanning trajectory
JP3547455B2 (en) Reconstruction of CT images of slices with turning motion
US7106825B2 (en) Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system
EP1643446B1 (en) Image reconstruction method and x-ray computed tomography apparatus
US6078638A (en) Pixel grouping for filtering cone beam detector data during 3D image reconstruction
US6285733B1 (en) Computed tomography method utilizing a conical radiation beam
US6018561A (en) Mask boundary correction in a cone beam imaging system using simplified filtered backprojection image reconstruction
US6130930A (en) Exact region of interest cone beam imaging without circle scans
US20090297009A1 (en) Method of reconstructing an image function from radon data
EP0892966A1 (en) Technique and arrangement for tomographic imaging
JP2000081318A (en) Scanning and data collecting method for three- dimensional computer tomography imaging and imaging system
JP4553894B2 (en) System and method for helical cone beam computed tomography with accurate reconstruction
JP2000081318A5 (en)
WO2001006931A1 (en) Volumetric ct image reconstruction
US6904117B2 (en) Tilted gantry helical cone-beam Feldkamp reconstruction for multislice CT
US6426989B2 (en) Computed tomography method
JP4342164B2 (en) Computed tomography equipment
Bruder et al. Single-slice rebinning reconstruction in spiral cone-beam computed tomography
JP4557321B2 (en) Image reconstruction device
JP2001057976A (en) Method and apparatus for stereoscopic image reconstruction, and ct scanner
JP2001054517A (en) Computer tomographic image forming device and operating method
JP4347061B2 (en) Sequential computed tomography method
US6333960B1 (en) Exact region of interest cone beam imaging without circle scans
JP4146907B2 (en) Computer tomography method including spiral scanning process of diagnostic area

Legal Events

Date Code Title Description
AS Assignment

Owner name: U.S. PHILIPS CORPORATION, NEW YORK

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:DANIELSSON, PER-ERIK;REEL/FRAME:010335/0905

Effective date: 19990125

FPAY Fee payment

Year of fee payment: 4

FPAY Fee payment

Year of fee payment: 8

REMI Maintenance fee reminder mailed
LAPS Lapse for failure to pay maintenance fees
STCH Information on status: patent discontinuation

Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362

FP Lapsed due to failure to pay maintenance fee

Effective date: 20130814