WO2005023114A2 - Computer tomography method using a cone-shaped bundle of rays - Google Patents
Computer tomography method using a cone-shaped bundle of rays Download PDFInfo
- Publication number
- WO2005023114A2 WO2005023114A2 PCT/IB2004/051523 IB2004051523W WO2005023114A2 WO 2005023114 A2 WO2005023114 A2 WO 2005023114A2 IB 2004051523 W IB2004051523 W IB 2004051523W WO 2005023114 A2 WO2005023114 A2 WO 2005023114A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- radiation source
- trajectory
- measured values
- partial
- rays
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000002591 computed tomography Methods 0.000 title claims abstract description 38
- 230000005855 radiation Effects 0.000 claims abstract description 129
- 230000003797 telogen phase Effects 0.000 claims description 26
- 238000001914 filtration Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 4
- 230000003205 diastolic effect Effects 0.000 claims description 2
- 230000000875 corresponding effect Effects 0.000 description 9
- 239000013598 vector Substances 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 230000029058 respiratory gaseous exchange Effects 0.000 description 5
- 206010073306 Exposure to radiation Diseases 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
- A61B6/527—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion using data from a motion artifact sensor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4021—Arrangements for generating radiation specially adapted for radiation diagnosis involving movement of the focal spot
- A61B6/4028—Arrangements for generating radiation specially adapted for radiation diagnosis involving movement of the focal spot resulting in acquisition of views from substantially different positions, e.g. EBCT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Definitions
- the invention relates to a computer tomography method in which an examination area, in which a patient may be located, is irradiated by a cone-shaped bundle of rays and a radiation source which emits the bundle of rays moves relative to the examination area. Measured values which depend on the intensity of the radiation transmitted through the examination area are acquired by a suitable detector unit, and a CT image of the examination area is reconstructed from these measured values. Moreover, the invention relates to a computer tomography scanner for carrying out the method and to a computer program for controlling the computer tomography scanner.
- the radiation source location of the radiation source moves relative to the examination area on a circular or helix-like trajectory.
- the circular trajectory has the disadvantage that the regions of the examination area which do not lie in the plane of the trajectory cannot be reconstructed using exact methods. Although these regions can be reconstructed using approximative methods, this leads to artefacts in the CT images. In this connection, a data record acquired along a circular trajectory is not complete.
- a data record is complete when, for each point in the examination area that is to be reconstructed, each plane containing this point intersects the trajectory at least once, and when this point is passed through by rays which come from at least one of these points of intersection and the absorption of these rays is measured. If, on the other hand, the radiation source location moves on a helix-like trajectory during the acquisition, the acquired measured values are complete and the examination area can be reconstructed using exact reconstruction methods, for example using the reconstruction method published in the article "Analysis of an exact inversion algorithm for spiral cone-beam CT", A. Katsevich, Phys. Med. Biol. 47, 2583-2597, 2002 (El).
- this object is achieved by a computer tomography method comprising the following steps: • producing a cone-shaped bundle of rays by means of a radiation source, said bundle of rays passing through an examination area and coming from a radiation source location, • producing a relative movement between the radiation source location on the one hand and the examination area on the other hand along an overall trajectory, wherein the relative movement comprises a rotation about an axis of rotation along a first, closed partial trajectory, a rotation about the axis of rotation along at least one second, closed partial trajectory and a movement along at least one third partial trajectory which connects the first and the at least one second partial trajectories to one another, • acquiring measured values which depend on the intensity in the bundle of rays on the other side of the examination area, by means of a detector unit during the relative movement, • reconstructing a CT image of the examination area from the measured values.
- each plane containing this point intersects the overall trajectory at least once. It is thus possible with this overall trajectory, unlike in the case of a circular overall trajectory, to acquire a complete data record if the respective point is passed through by rays which come from at least one of the respective points of intersection.
- the examination area which lies between the first and the at least one second partial trajectory can thus be reconstructed using exact reconstruction methods, and this leads to an improved image quality.
- the radiation source location can, unlike in the case of the helixlike trajectory, rotate along the first and the at least one second trajectory about the axis of rotation until a sufficient amount of measured values to reconstruct the trajectory planes has been acquired.
- the radiation source location can rotate relative to the examination area along the first and the at least one second trajectory until a desired amount of measured values which correlate with a rest phase of the object has been acquired for the reconstruction.
- the third partial trajectory is also closed so that the radiation source location can also rotate along this partial trajectory until a desired amount of measured values has been acquired for the reconstruction.
- the movement of the radiation source location relative to the examination area on closed partial trajectories allows simpler technical implementation of the movement since no reversal of the rotation movement of the radiation source is required and in particular the radiation source can rotate with a constant angular velocity.
- the first and the at least one second partial trajectory in each case form a circle, which makes it possible to further simplify implementation of the method by means of a computer tomography scanner and leads to a reduced computational complexity during the reconstruction.
- Claim 4 describes a reconstruction using a filtered back-projection, which leads to a further improvement in the image quality while having a low degree of computational complexity. A further improvement in image quality can be achieved with the embodiments as claimed in claims 5 and 6.
- measured values are essentially acquired only when the object is in a rest phase.
- the radiation source location moves along the first and the at least one second partial trajectory until the radiation source locations on the first and the at least one second trajectory that belong to the acquired measured values have scanned a predefined angular range or predefined angular ranges.
- These angular ranges may be selected such that a periodically moving object can be reconstructed with good image quality. They may be selected, for example, such that from each point in the examination area that is to be reconstructed the radiation source location can be seen during the rest phases over a continuous or total angular range of 180°.
- Claim 8 provides that a dynamic collimator arrangement collimates the cone- shaped bundle of rays such that for each radiation source location the cone-shaped bundle of rays is delimited by two adjacent partial trajectories selected from the first and the at least one second partial trajectory. This leads, for example, in respect of a human patient, to a reduction of the exposure to radiation.
- a computer tomography scanner for carrying out the method according to the invention is described in claim 9.
- the embodiment as claimed in claim 10 allows a movement of the radiation source location of the radiation source, in particular a focal spot of an X-ray source, and thus a movement of the radiation source location relative to the examination area along the at least one third partial trajectory, without an object located in the examination area having to be moved in the direction of the axis of rotation, for example by moving a patient table.
- the acquisition process is therefore more pleasant for a patient, and an undesired movement of a patient relative to a patient table becomes more unlikely, so that the image quality is further improved. Since the focal spot of an X-ray source can be moved more quickly than a physical object, this preferred embodiment furthermore leads to more rapid acquisition of the measured values.
- the preferred embodiment of a computer tomography scanner as claimed in claim 11 allows the detection and recording of the movement of an object. It is thus possible to acquire measured values essentially only when the object is in a rest phase or to use for reconstruction only measured values which have been acquired during a rest phase. This leads to a further improvement in image quality.
- the computer tomography scanner has a dynamic collimator arrangement which is adapted such that for each radiation source location the cone-shaped bundle of rays is delimited by two adjacent partial trajectories selected from the first and the at least one second partial trajectory. As mentioned above, this leads to a reduction in the exposure to radiation.
- Claim 13 defines a computer program for controlling a computer tomography scanner as claimed in claim 9.
- Fig. 1 shows a computer tomography scanner according to the invention which can be used to carry out the method according to the invention.
- Fig. 2 shows an overall trajectory according to the invention on which the radiation source location can move relative to the examination area.
- Fig. 3 shows a further overall trajectory according to the invention on which the radiation source location can move relative to the examination area.
- Fig. 4 shows a further overall trajectory on which the radiation source location can move relative to the examination area.
- Fig. 5 shows a flowchart of the method according to the invention.
- Fig. 6 shows an overall trajectory according to the invention with a radiation source location and a bundle of rays.
- FIG. 7 shows an overall trajectory with a radiation source location and a virtual, planar detector surface.
- Fig. 8 shows the virtual, planar detector surface onto which the second partial trajectory is projected, viewed from a point on a first partial trajectory.
- Figs. 9-11 show the virtual, planar detector surface onto which the second partial trajectory is projected, with filter lines and directions, viewed from a point on the first partial trajectory.
- Fig. 12 shows the virtual, planar detector surface onto which the first and second partial trajectories are projected, with filter lines and directions, viewed from a point on a third partial trajectory.
- Fig. 13 shows the virtual, planar detector surface onto which the first and second partial trajectories are projected, viewed from a point on the third partial trajectory.
- the computer tomography scanner shown in Fig. 1 comprises a gantry 1 which can rotate about an axis of rotation 14 that runs parallel to the z direction of the coordinate system shown in Fig. 1.
- the gantry 1 is driven by a motor 2 at a preferably constant but adjustable angular velocity.
- Attached to the gantry 1 is a radiation source S, for example an X-ray source.
- the latter is provided with a collimator arrangement 3 which screens out, from the radiation generated by the radiation source S at the radiation source location, a cone-shaped bundle of rays 4, that is to say a bundle of rays which has a finite extent other than zero both in the z direction and in a direction pe ⁇ endicular thereto (that is to say in a plane pe ⁇ endicular to the axis of rotation 14).
- the bundle of rays 4 penetrates a cylindrical examination area 13 in which a periodically moving object (not shown) is located in this embodiment according to the invention.
- This object is in this case a beating heart which carries out intrinsic movements and under some circumstances is also moved back and forth on account of breathing movements of the patient.
- non-moving objects or other periodically moving body organs such as the liver or brain, periodically moving parts of body organs or periodically moving technical objects may also be irradiated.
- the computer tomography method according to the invention can also be carried out when an object which is not moving is located in the examination area 13.
- objects within which other objects are moving may be irradiated.
- the propagation of a contrast agent in the human body can be observed using the computer tomography scanner (perfusion examination). For this pu ⁇ ose, the body is generally irradiated over a relatively long period of time with low-dose radiation.
- the bundle of rays 4 impinges on a detector unit 16 which is attached to the gantry 1 and has a detector surface that comprises a large number of detector elements, said detector elements being arranged in rows and columns in the manner of a matrix in this embodiment.
- the detector columns preferably run parallel to the axis of rotation 14.
- the detector rows are located in planes pe ⁇ endicular to the axis of rotation and in this embodiment are arranged in an arc of a circle about the radiation source S (focus-centered detector). In other embodiments, however, they may also be shaped differently, e.g. describe an arc of a circle about the axis of rotation 14 or be rectilinear.
- Each detector element that is struck by the bundle of rays 4 supplies in each position of the radiation source location a measured value for a ray of the bundle of rays 4.
- the opening angle of the bundle of rays designated ma ⁇ , determines the diameter of the object cylinder within which the object to be examined is located during the acquisition of the measured values.
- the opening angle is defined as the angle enclosed by a ray lying in a plane pe ⁇ endicular to the axis of rotation 14 at the edge of the bundle of rays and a plane running through the radiation source S and the axis of rotation 14.
- the examination area 13 or the object or the patient table may be displaced by means of a motor 5 parallel to the axis of rotation 14 or the z axis.
- the gantry may also be displaced in this direction.
- the X-ray source has a focal spot that can be moved parallel to the axis of rotation, from which focal spot the X-rays come.
- the radiation source location can thus be displaced parallel to the axis of rotation.
- the radiation source S and the focal spot of the X-ray source and the detector unit 16 can describe an overall trajectory 20 relative to the examination area 13, which overall trajectory runs on an imaginary cylindrical surface. Examples of such overall trajectories are shown in Figs. 2 to 4.
- first closed partial trajectory 21, at least one second closed partial trajectory 25, 26, 29 and at least one third partial trajectory 23, 27, 31, 32 which connects the first partial trajectory 21 and the at least one second partial trajectory 25, 26, 29 to one another.
- the first partial trajectory 21 and the at least one second partial trajectory 25, 26, 29 are arranged around the axis of rotation 14.
- two circular first and second partial trajectories 21, 25 with the same diameter and the same axis of rotation 14 are connected by a third elliptical partial trajectory 23.
- two circular partial trajectories 21, 25 with the same axis of rotation 14 and the same diameter are connected to one another by an ellipse section 27.
- a further overall trajectory according to the invention is shown in Fig.
- the heart movement is recorded by an electrocardiograph 8 in a known manner.
- the chest area of a patient is connected to the electrocardiograph 8 via electrodes (not shown).
- the movement of the object may be monitored in other ways.
- the movement information could for example be obtained from the values measured by the detector unit 16 themselves, so that there is no need to record the movement by means of an additional device such as an electrocardiograph.
- a kymogram is first created from the measured values, and the movement can be derived from said kymogram in a known manner.
- the measured values acquired by the detector unit 16 are fed to a reconstruction computer 10 which is connected to the detector unit 16 for example via a contactless data transmission system (not shown). Moreover, an electrocardiogram is transmitted from the electrocardiograph 8 to the reconstruction computer 10.
- the reconstruction computer 10 reconstructs a CT image of the examination area 13, in particular an abso ⁇ tion distribution, and displays it for example on a monitor 11.
- the two motors 2 and 5, the reconstruction computer 10, the radiation source S, the electrocardiograph 8 and the transfer of the measured values from the detector unit 16 to the reconstruction computer 10 are controlled by the control unit 7. Moreover, the control unit 7 controls transmission of the electrocardiogram from the electrocardiograph 8 to the reconstruction computer 10. Fig.
- the gantry 1 rotates at an angular velocity which in this example of embodiment is constant. However, it may also vary, for example as a function of time or of the radiation source position or of the position of the radiation source location.
- the radiation of the radiation source S is switched on and the radiation source location is moved relative to the examination area 13 along the overall trajectory 20 by moving the patient table or the focal spot.
- the electrocardiograph 8 is activated so that an electrocardiogram is recorded at the same time as the measured values are being acquired.
- the distance between the first partial trajectory 21 and the second partial trajectory 25 is selected such that the object to be reconstructed can be arranged completely in the part of the examination area which is delimited by these two partial trajectories 21, 25, since, as already mentioned above, the acquired data record is complete only in respect of this part.
- this distance should be selected such that the projections of the first partial trajectory 21 and of the second partial trajectory 25 along the rays emitted by the radiation source S fall on the detector unit.
- the distance between these two partial trajectories 21 and 25 and/or the extent of the detector unit in the direction of the axis of rotation 14 should be selected such that the part 33 of the cone-shaped bundle of rays 4 shown by way of example i Fig.
- the radiation source S is in this way controlled by the control unit 7, to which the electrocardiogram is continuously transmitted by the electrocardiograph 8. If a temporally precise switching on and off of the radiation source is technically problematic, so that acquisition of the measured values when the object is in a rest phase is not possible, the radiation source is controlled such that at least the majority, in particular more than 80% or 90%, of the measured values are acquired while the object is in a rest phase, so that measured values are essentially acquired only in the rest phases. As an alternative, the measured values may be acquired independently of the object movement, although in this case only measured values acquired while the object was in a rest phase are taken into account in the reconstruction described below.
- the radiation source S or more precisely the radiation source location B, moves relative to the examination area along the first partial trajectory 21.
- the radiation source location B rotates along this partial trajectory 21 until measured values have been acquired during one or more rest phases of the object over as large an angular range as possible.
- the radiation source S moves on the partial trajectory 21 until measured values which correlate with the one or more rest phases of the object have been acquired over the entire angular range of 360°.
- the patient table or focal spot of the X-ray source is moved parallel to the axis of rotation 14, so that the radiation source location advances relative to the examination area 13 along an ellipse 23.
- the gantry 1 may furthermore rotate at a constant angular velocity without the rotation movement having to be accelerated or slowed down.
- the overall trajectory 20 according to the invention can therefore be realized relatively easily with known computer tomography scanners. While the radiation source location is moving relative to the examination area 13 along the third partial trajectory 23 also, measured values are acquired only when the object is in a rest phase. The radiation source location moves on this ellipse 23 until measured values have been acquired over as large an angular range as possible, in particular over the entire angular range of 360°, along the partial trajectory 23.
- the patient table and focal spot move so that the radiation source location is guided along the third partial trajectory 23 toward the second partial trajectory 25.
- the patient table and focal spot end their movement parallel to the axis of rotation, so that the radiation source location is now moving along the second partial trajectory 25.
- the measured values are acquired only when the object is in a rest phase, and measurement is continued until measured values have been acquired over as large an angular range as possible along the second partial trajectory 25, in particular over an angular range of 360°.
- the sequence in which the partial trajectories are passed through is not important.
- the measured values could therefore also be acquired, for example, first along the second partial trajectory 25 or along the third partial trajectory 23. If measured values are acquired a number of times from one and the same position of the radiation source location, these measurement values which are redundant per se may be used to increase the signal-to-noise ratio. For this pu ⁇ ose, measured values whose rays come from the same position of the radiation source location and strike the same detector element may be averaged, for example arithmetically, to form one measured value. Once the measured values have been acquired, a CT image is reconstructed from these measured values. This could take place for example using a known iterative reconstruction method.
- an iteration step may comprise the following steps: a) simulating the acquisition along the overall trajectory to produce a simulated set of measured values; b) difference formation between the simulated and the corresponding measured values that were actually measured; c) back-projection of the differences along the rays belonging to the respective differences and measured values into the examination area.
- ART Algebraic Reconstruction Technique
- the reconstruction is carried out using a filtered back-projection.
- the measured values are firstly partially derived. This partial derivation is explained in more detail below.
- the first partial trajectory 21, y, , and the second partial trajectory 25, y 2 can be described by the following equations:
- the angular positions s, and s 2 in each case assume values of between 0 and 2 ⁇ .
- the variable R denotes the radius of the first 21 and second 23 circular partial trajectories, and z 0 is the z position at which the trajectory plane of the second partial trajectory 25 and the axis of rotation 14 intersect.
- the elliptical third partial trajectory 23 can be described by the following equation:
- the angular position s ⁇ of the ellipse 23 may assume values of between 0 and 2 ⁇ . If instead of a complete ellipse only an ellipse section 27 is used as the third partial trajectory, then the angular position s 3 assumes a correspondingly smaller range of values.
- the angular position s assumes values of between 0 and ⁇ .
- the respectively associated measured values D f (y l (s l ), ⁇ ) can be described by
- D f (y l (s l ),Q) ldlf(y, (s l ) + l ⁇ ) .
- the index assumes the values 1, 2 and 3
- ⁇ is a unit vector which points in the different ray directions of the rays coming from the respective radiation source location.
- the function f(x) denotes the abso ⁇ tion value at the point x in the examination area 13.
- the partial derivation in step 103 can then be described by the following equation: Parallel rays which come from different positions of the radiation source location are therefore derived on the basis of the angular position s, of the radiation source location B.
- measured values whose rays run parallel and come from the same partial trajectory but from different radiation source locations can be combined to form in each case one group.
- the measured values of a group are in each case derived on the basis of the angular position s, of the radiation source location on the respective partial trajectory. This derivation may be carried out, for example, using a Fourier transform or the finite element method.
- the measured values are projected along the respective rays onto a virtual, planar detector surface 35.
- the virtual detector surface 35 contains the axis of rotation 14 and rotates with the radiation source S about the axis of rotation 14.
- a ray of the radiation source S is always pe ⁇ endicular on the axis of rotation 14 and the virtual detector surface 35.
- the virtual detector surface 35 does not move in the direction of the axis of rotation 14 and is selected to be so large that all the rays of the bundle of rays 33 shown in Fig. 6 and described above penetrate the virtual detector surface 35.
- the size of the virtual detector surface 35 must therefore be selected such that all the rays which come from one point on the overall trajectory 20 and pass through the examination area and leave through the region delimited by the partial trajectories 21 and 25 are detected. In this example of embodiment this detector surface 35 is rectangular.
- the derived measured values projected onto the virtual detector surface 35 are then filtered in step 105.
- the filter lines and directions along which the measured values are filtered on the detector surface 35 depend on which partial trajectory 21, 23, 25 the radiation source location was located during acquisition of the respective measured value, in which direction it is moving on this partial trajectory and on the position of the measured value on the virtual detector surface 35.
- the detector surface 35 is divided into a number of regions with the aid of projections of the first partial trajectory 21 and second partial trajectory 25 along the rays coming from the respective radiation source location.
- the detector surface 35 is divided into three regions 45, 47, 49 by the projection 37 of the second partial trajectory 25 onto the detector surface 35 along the rays coming from the first partial trajectory 21 and by a straight line 43 which contacts this projection 37 at a tangent at one point and runs pe ⁇ endicular to the axis of rotation 14. These regions are shown in Fig. 8. In this figure and also in Figs. 9 to 13, in each case the virtual detector surface 35 can be seen, wherein the direction from left to right in these figures corresponds essentially to the direction in which the radiation source location advances.
- the direction from left to right corresponds to the component of the advance movement of the radiation source location which is oriented pe ⁇ endicular to the axis of rotation.
- a measured value lies or is arranged in a certain region on the virtual detector surface 35, this means that the measured value originally acquired by the detector unit 16 has been projected into this region in step 104.
- a filter line 51 (Fig. 9) is assigned to a measured value D whose ray comes from the first partial trajectory 21 and which lies in the region 47 delimited by the straight line 43, said filter line running through this measured value and parallel to the straight line 43 and pe ⁇ endicular to the axis of rotation 14.
- Two filter lines 55, 57 or 59, 61 are assigned to a measured value D whose ray comes from the first partial trajectory and which lies in one of the regions 45, 49 that are delimited by the straight line 43 and the projection 37 of the second partial trajectory 25, said filter lines in each case running at a tangent to the projection 37 (cf. Fig. 10 and Fig. 11).
- the two tangents of the projection 37 are therefore assigned as filter lines 55, 57 or 59, 61 to such a measured value D, said filter lines running through this measured value D.
- the filter directions 56, 58 and 60, 62 belonging to these filter lines 55, 57 and 59, 61 are likewise shown in Fig. 10 and Fig. 11.
- all measured values whose rays come from the first partial trajectory 21, regardless of their arrangement on the planar, virtual detector surface 35 could be filtered along filter lines which run pe ⁇ endicular to the axis of rotation 14. That is to say, in Fig. 9, the filter lines 51 and filter directions 52 would then be distributed over the entire detector surface 35.
- filter lines which is described in more detail below, is of course possible only when measured values are arranged along the filter lines.
- the opening angle of the cone-shaped bundle of rays 4 must therefore be selected to be so large that measured values are acquired and arranged also in the region of the virtual detector surface that is delimited by the projection 37 and the upper edge of the detector surface 35.
- the filter lines and filter directions of those measured values whose rays come from the second partial trajectory 25 correspond to the filter lines and filter directions for measured values whose rays come from a radiation source position on the first partial trajectory 21. That is to say, in order to obtain the corresponding filter lines and filter directions for these measured values, in Figs.
- Filter lines of the first set run, for a measured value D which lies in the regions 67 or 69, at a tangent to the projection 37 of the second partial trajectory 25 onto the detector surface 35, with the point of contact of this filter line 64 lying to the right of the respective measured value 80 in Fig. 12. If, on the other hand, the measured value lies in one of the regions 63 or 65, then the filter line of the first set runs at a tangent to the projection 39 of the first partial trajectory 21 onto the detector surface 35, with the point of contact of this filter line with the projection 39 lying to the left of the measured value in Fig. 12.
- the filter directions of the filter lines of the first set run essentially from right to left.
- the filter line of the second set runs at a tangent to the projection 37, with the point of contact of this filter line with the projection 37 in Fig. 12 lying to the left of the measured value.
- the filter line of the second set runs at a tangent to the projection 39, with the point of contact of this filter line 66 with the projection 39 being to the right of the measured value.
- the filter directions of the filter lines of the second set run from right to left for measured values whose rays come from the section of the third partial trajectory 23 on which the radiation source location moves in the direction of the second partial trajectory 25. In Fig. 12 they run from left to right for measured values whose rays come from the section of the third partial trajectory 23 on which the radiation source location moves in the direction of the first partial trajectory 21.
- measured values may be arranged in regions 82, 84, which lie between the projections 37, 39, but not in one of the regions 63, 65, 67, 69.
- the measured values which lie in these regions 82, 84 are in this example of embodiment not taken into account in the filtering and back-projection. If the overall trajectory has a number of second partial trajectories 26, 29, as shown in Fig.
- each part of the examination area 13 that is delimited by two adjacent partial trajectories selected from the first and second partial trajectories can then be considered individually, with the method according to the invention being carried out for each of these parts. That is to say, in step 105, in the case of one first partial trajectory 21 and two second partial trajectories 26 and 29, two virtual detector surfaces 35 may be defined which are each passed through by the bundle of rays delimited by the adjacent partial trajectories 21, 29 and 29, 26. The respectively adjacent partial trajectories 21, 29 or 29, 26 are then projected onto each of these virtual detector surfaces 35 and then, as described above, filter lines and filter directions can be defined on the respective virtual detector surface. In another embodiment, if there are a number of second partial trajectories 26,
- (y 1 (5,), ⁇ (y j (5 , ,),x)) is the respectively filtered measured value whose ray comes from y, (s t ) and runs through the point x in the examination area 13.
- the redundancy factor ⁇ may take account of redundancies arising during filtering or reconstruction. Examples of possible redundancy factors are mentioned below.
- e ⁇ (y , (s t ), x) is a vector which describes the course of the respective filter line and is defined by the following equation:
- the K angle is the angle between the ray belonging to the measured value that is to be filtered and the ray belonging to the respective measured value that is to be multiplied by the first weighting factor and lies on the filter line.
- the filtered measured values may in each case be determined at points on the virtual detector surface which lie between the measured values lying on the filter line.
- each measured value lying on the respective filter line is multiplied by a second weighting factor which is equal to the K angle distance d ⁇ or ⁇ between two adjacent measured values lying on the filter line.
- the measured values are arranged on a filter line in an equidistant manner so that the K angle distance between two adjacent measured values is constant for a filter line. If the K angle distances between a measured value and the two adjacent measured values are different, then the K angle distance in the filter direction can be taken as the weighting factor for example, or else the weighting factor can be determined by arithmetic averaging of the two K angle distances. If in other embodiments the K angle distance is the same for all filter lines, the multiplication by the second weighting factor could be omitted. The summation of the weighted measured values lying on the respective filter line is repeated for all filter lines of a measured value, so that for one measured value a number of sums equal to the number of filter lines is determined.
- the result is just one sum of weighted measured values lying on the filter line, since only one filter line is assigned to this measured value.
- two filter lines are assigned to a measured value which lies for example in the region 49 on the detector surface 35 in Fig. 8 or Fig. 10, so that after filtering two corresponding sums are produced for this measured value. If a number of such sums are produced for a measured value, they are then added in accordance with equation (6). The resulting sum is then multiplied by a redundancy factor and finally forms the filtered measured value.
- the redundancy factor depends on which radiation source location belongs to the respective measured value and where the measured value is arranged on the virtual detector surface. If the ray belonging to a measured value comes from the first partial trajectory 21, then the redundancy factor is 1/2 when the measured value is arranged in the region 47 in Fig. 8, and is 1/4 when the measured value is arranged in one of the regions 45 or 49 in Fig. 8. The same applies in respect of measured values whose rays come from the second partial trajectory 25, wherein, as already described above, the virtual detector surface 35 in Fig. 8 has to be reflected along a horizontal axis that is arranged centrally on the detector surface, in order to obtain the corresponding regions 45, 47 and 49. If the ray belonging to a measured value comes from the third partial trajectory
- the redundancy factor for this measured value is 1/4. If the third trajectory forms only a half-ellipse, as shown in Fig. 3, then the redundancy factor for rays which come from the third trajectory is 1/2.
- the multiplication by the weighting factors and the summation or filtering of the measured values in accordance with equation (6) may be carried out in a known manner using a Fourier transform. Once all the measured values on the virtual, planar detector surface 35 have been filtered along the filter lines assigned to them, the filtered measured values are used to reconstruct the abso ⁇ tion distribution in the examination area by means of a back-projection.
- f(x) denotes the value to be reconstructed, in particular the abso ⁇ tion value, at the point x in the examination area.
- f(x) denotes the value to be reconstructed, in particular the abso ⁇ tion value, at the point x in the examination area.
- a voxel V(x) is determined within a predefinable field of view FOV in the examination area 13 between the first partial trajectory 21 and the second partial trajectory 25.
- a radiation source location which has not yet been taken into account in this back-projection is selected on the overall trajectory (or a projection).
- a projection is in this connection the set of all the rays which come from a radiation source location.
- a specific radiation source location is selected only if the part 33 (Fig. 6) of the cone-shaped bundle of rays 4 coming from this radiation source location that is delimited by the two adjacent first and second partial trajectories 21, 25 passes through the voxel V(x).
- step 107 only projections which have been acquired during a rest phase of the object should be selected. This may be done with the aid of the electrocardiogram recorded by the electrocardiograph 8. Even if in step 102 projections have been acquired only during rest phases of the object, this can be checked in step 107 by means of the electrocardiogram, that is to say a check can be made as to whether the projections have actually been acquired during the rest phases of the object. Projections which have not been acquired during a rest phase may continue to be disregarded in the following back-projection.
- step 108 a check is made as to whether there is, in the projection selected in step 107, a measured value whose ray runs through the center of the voxel V(x). If such a measured value or ray cannot be found, the point at which a central ray would strike the detector surface 35 is determined. The associated measured value is then calculated by inte ⁇ olating the measured values of adjacent rays.
- the measured value which can be assigned to the ray passing through the voxel, or the measured value obtained by inte ⁇ olation, is multiplied in step 109 by a weighting factor which becomes smaller as the distance of the radiation source location from the voxel position x in the FOV increases, and in particular is selected to be proportional to the reciprocal value of this distance or equal to said reciprocal value.
- the measured value multiplied by this weighting factor is added to the voxel ⁇ (x). Initially, all voxels V(x) are set to the same initial value, for example to the value 0.
- step 111 a check is made as to whether all the projections which can be taken into account for the back-projection, as mentioned above in relation to step 107, have been considered. If this is not the case, the flowchart branches to step 107. Otherwise, it continues with step 112.
- step 112 a check is made as to whether all voxels V(x) in the FOV have been passed through. If this is not the case, the method continues with step 106. If, on the other hand, all voxels V(x) in the FOV have been passed though, then the abso ⁇ tion in the entire FOV has been determined and the reconstruction method is terminated (step 113).
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Physiology (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Mathematical Optimization (AREA)
- Pulmonology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
Claims
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006525236A JP4606414B2 (en) | 2003-09-05 | 2004-08-23 | Computed tomography method using conical ray bundle |
US10/570,088 US7403587B2 (en) | 2003-09-05 | 2004-08-23 | Computer tomography method using a cone-shaped bundle of rays |
EP04744809A EP1663004A2 (en) | 2003-09-05 | 2004-08-23 | Computer tomography method using a cone-shaped bundle of rays |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP03102708.9 | 2003-09-05 | ||
EP03102708 | 2003-09-05 |
Publications (3)
Publication Number | Publication Date |
---|---|
WO2005023114A2 true WO2005023114A2 (en) | 2005-03-17 |
WO2005023114A3 WO2005023114A3 (en) | 2005-07-07 |
WO2005023114A8 WO2005023114A8 (en) | 2006-06-29 |
Family
ID=34259205
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IB2004/051523 WO2005023114A2 (en) | 2003-09-05 | 2004-08-23 | Computer tomography method using a cone-shaped bundle of rays |
Country Status (4)
Country | Link |
---|---|
US (1) | US7403587B2 (en) |
EP (1) | EP1663004A2 (en) |
JP (1) | JP4606414B2 (en) |
WO (1) | WO2005023114A2 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007151849A (en) * | 2005-12-06 | 2007-06-21 | Ge Medical Systems Global Technology Co Llc | X-ray ct imaging method and x-ray ct apparatus |
WO2007119124A1 (en) * | 2006-04-19 | 2007-10-25 | Philips Intellectual Property & Standards Gmbh | Cone beam computed tomography with multiple partial scanning trajectories |
WO2008075244A2 (en) * | 2006-12-15 | 2008-06-26 | Koninklijke Philips Electronics N.V. | Imaging system for imaging an object in an examination zone |
WO2009083850A1 (en) * | 2007-12-21 | 2009-07-09 | Koninklijke Philips Electronics, N.V. | Dynamic collimation in cone beam computed tomography to reduce patient exposure |
US8401144B2 (en) | 2008-08-07 | 2013-03-19 | Koninklijke Philips Electronics N.V. | Method and apparatus for correcting artifacts in circular CT scans |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007505674A (en) * | 2003-09-18 | 2007-03-15 | コニンクリユケ フィリップス エレクトロニクス エヌ.ブイ. | Imaging processing method and apparatus |
EP1819276A1 (en) * | 2004-10-15 | 2007-08-22 | Philips Intellectual Property & Standards GmbH | Computer tomography method |
EP1987772A4 (en) * | 2006-02-20 | 2012-08-01 | Shimadzu Corp | Tomographic device and processing program |
US7848479B1 (en) * | 2007-06-29 | 2010-12-07 | University Of Central Florida Research Foundation, Inc. | Image reconstruction for a general circle-plus trajectory |
JP2011502683A (en) * | 2007-11-19 | 2011-01-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Imaging apparatus, imaging method, and computer program for determining image of region of interest |
EP2234541B1 (en) * | 2007-12-20 | 2016-05-04 | Koninklijke Philips N.V. | Stereo tube attenuation filter |
US8139709B2 (en) * | 2008-09-15 | 2012-03-20 | University Of Utah Research Foundation | Staggered circular scans for CT imaging |
US7831011B2 (en) * | 2008-11-21 | 2010-11-09 | General Electric Co. | Computed tomography method and system |
US20100202583A1 (en) * | 2009-02-03 | 2010-08-12 | Ge Wang | Systems and Methods for Exact or Approximate Cardiac Computed Tomography |
US7881426B2 (en) * | 2009-02-26 | 2011-02-01 | Morpho Detection, Inc. | Method and system for performing a scan of an object |
DE102009020400B4 (en) * | 2009-05-08 | 2016-04-21 | Siemens Aktiengesellschaft | Method and apparatus for image determination from x-ray projections taken when traversing a trajectory |
JP5507682B2 (en) | 2009-06-30 | 2014-05-28 | アナロジック コーポレーション | Efficient quasi-rigid 3D image reconstruction algorithm for CT scanner |
DE102009042922B4 (en) * | 2009-09-24 | 2019-01-24 | Siemens Healthcare Gmbh | Method and apparatus for image determination from x-ray projections taken when traversing a trajectory |
EP2315178A1 (en) * | 2009-10-12 | 2011-04-27 | Siemens Aktiengesellschaft | Reconstruction of 3D image datasets from X-ray cone-beam data |
US20120099768A1 (en) | 2010-10-20 | 2012-04-26 | Medtronic Navigation, Inc. | Method and Apparatus for Reconstructing Image Projections |
JP5674546B2 (en) * | 2011-04-27 | 2015-02-25 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Measuring apparatus, program, tomography system, tomography apparatus and network system |
BR112014008090A2 (en) * | 2011-10-03 | 2017-04-18 | Fei Co | computed tomography imaging process and system, and computer readable storage media |
CN105093342B (en) * | 2014-05-14 | 2017-11-17 | 同方威视技术股份有限公司 | Spiral ct system and method for reconstructing |
EP3133554B1 (en) * | 2015-08-17 | 2019-07-10 | FEI Company | Novel acquisition and processing of data in a tomographic imaging apparatus |
US10383590B2 (en) * | 2015-09-28 | 2019-08-20 | General Electric Company | Methods and systems for adaptive scan control |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5187659A (en) * | 1990-09-04 | 1993-02-16 | General Electric Company | Cone beam scanning trajectories for three-dimensional computerized tomography data acquisition where object is larger than the field of view |
US5375156A (en) * | 1992-03-31 | 1994-12-20 | Siemens Medical Systems, Inc. | Method and apparatus for 3-D computer tomography |
US5278884A (en) * | 1992-12-18 | 1994-01-11 | General Electric Company | Complete 3D CT data acquisition using practical scanning paths on the surface of a sphere |
US5611026A (en) * | 1992-12-21 | 1997-03-11 | General Electric Company | Combining a priori data with partial scan data to project three dimensional imaging of arbitrary objects with computerized tomography |
US6321245B1 (en) * | 1997-04-02 | 2001-11-20 | International Business Machines Corporation | Method and system for performing fast division using non linear interpolation |
US5802134A (en) | 1997-04-09 | 1998-09-01 | Analogic Corporation | Nutating slice CT image reconstruction apparatus and method |
US6014419A (en) | 1997-11-07 | 2000-01-11 | Hu; Hui | CT cone beam scanner with fast and complete data acquistion and accurate and efficient regional reconstruction |
US6009142A (en) * | 1998-03-31 | 1999-12-28 | Siemens Corporate Research, Inc. | Practical cone beam image reconstruction using local regions-of-interest |
US6084937A (en) | 1998-07-27 | 2000-07-04 | Siemens Corporate Research, Inc. | Adaptive mask boundary correction in a cone beam imaging system |
US6018561A (en) | 1998-07-27 | 2000-01-25 | Siemens Corporate Research, Inc. | Mask boundary correction in a cone beam imaging system using simplified filtered backprojection image reconstruction |
US6275560B1 (en) * | 1998-12-22 | 2001-08-14 | General Electric Company | Cardiac gated computed tomography system |
JP2000197627A (en) * | 1999-01-05 | 2000-07-18 | Hitachi Medical Corp | X-ray ct device |
JP4378812B2 (en) * | 1999-11-19 | 2009-12-09 | 株式会社島津製作所 | Cone beam type radiation CT system |
JP4407005B2 (en) * | 1999-11-25 | 2010-02-03 | 株式会社島津製作所 | Cone beam type radiation CT system |
DE10021219A1 (en) * | 2000-04-29 | 2001-10-31 | Philips Corp Intellectual Pty | Computer tomography procedure |
US6504892B1 (en) | 2000-10-13 | 2003-01-07 | University Of Rochester | System and method for cone beam volume computed tomography using circle-plus-multiple-arc orbit |
US6324245B1 (en) | 2000-10-17 | 2001-11-27 | Siemens Corporation Research, Inc. | Method and apparatus for identifying and correcting image inaccuracies caused by simplified processing of masked cone beam projection data |
US6330298B1 (en) | 2000-10-17 | 2001-12-11 | Siemens Corporate Research, Inc. | Method and apparatus for simplifying the correction of image inaccuracies caused by processing of masked cone beam projection data |
US6771733B2 (en) * | 2001-08-16 | 2004-08-03 | University Of Central Florida | Method of reconstructing images for spiral and non-spiral computer tomography |
US6459756B1 (en) | 2001-10-30 | 2002-10-01 | Siemens Corporate Research, Inc. | System and method for providing normalization correction for image reconstruction in a reduced pitch spiral scan cone beam computed tomography imaging system |
-
2004
- 2004-08-23 JP JP2006525236A patent/JP4606414B2/en not_active Expired - Fee Related
- 2004-08-23 WO PCT/IB2004/051523 patent/WO2005023114A2/en active Application Filing
- 2004-08-23 EP EP04744809A patent/EP1663004A2/en not_active Withdrawn
- 2004-08-23 US US10/570,088 patent/US7403587B2/en not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
None |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007151849A (en) * | 2005-12-06 | 2007-06-21 | Ge Medical Systems Global Technology Co Llc | X-ray ct imaging method and x-ray ct apparatus |
WO2007119124A1 (en) * | 2006-04-19 | 2007-10-25 | Philips Intellectual Property & Standards Gmbh | Cone beam computed tomography with multiple partial scanning trajectories |
JP2009534079A (en) * | 2006-04-19 | 2009-09-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Cone beam computed tomography with multiple partial scan trajectories |
US7822171B2 (en) | 2006-04-19 | 2010-10-26 | Koninklijke Philips Electronics N.V. | Cone beam computed tomography with multiple partial scanning trajectories |
WO2008075244A2 (en) * | 2006-12-15 | 2008-06-26 | Koninklijke Philips Electronics N.V. | Imaging system for imaging an object in an examination zone |
WO2008075244A3 (en) * | 2006-12-15 | 2008-08-14 | Koninkl Philips Electronics Nv | Imaging system for imaging an object in an examination zone |
WO2009083850A1 (en) * | 2007-12-21 | 2009-07-09 | Koninklijke Philips Electronics, N.V. | Dynamic collimation in cone beam computed tomography to reduce patient exposure |
US8213568B2 (en) | 2007-12-21 | 2012-07-03 | Koninklijke Philips Electronics N.V. | Dynamic collimation in cone beam computed tomography to reduce patient exposure |
US8401144B2 (en) | 2008-08-07 | 2013-03-19 | Koninklijke Philips Electronics N.V. | Method and apparatus for correcting artifacts in circular CT scans |
Also Published As
Publication number | Publication date |
---|---|
US20070019776A1 (en) | 2007-01-25 |
JP2007503905A (en) | 2007-03-01 |
WO2005023114A3 (en) | 2005-07-07 |
US7403587B2 (en) | 2008-07-22 |
EP1663004A2 (en) | 2006-06-07 |
JP4606414B2 (en) | 2011-01-05 |
WO2005023114A8 (en) | 2006-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7403587B2 (en) | Computer tomography method using a cone-shaped bundle of rays | |
US6466640B1 (en) | Computed tomography system and method | |
JP5019193B2 (en) | Reconstruction method and X-ray computed tomography apparatus for determining image data values at reconstruction points in a CT image related to a scan target | |
US5991356A (en) | Radiation tomography method and apparatus | |
EP1800264B1 (en) | Image reconstruction with voxel dependent interpolation | |
CN103390284B (en) | The CT image reconstruction measured in field in extension | |
US20080267455A1 (en) | Method for Movement Compensation of Image Data | |
US7672423B2 (en) | Short scan cardiac CT on a quasi axial trajectory | |
JP4440588B2 (en) | CT image forming apparatus and CT apparatus for subject moving periodically | |
JP2009534079A (en) | Cone beam computed tomography with multiple partial scan trajectories | |
JP2008546464A (en) | Suppression of band artifact in cardiac CT scan | |
JP2008519636A (en) | Computed tomography method for inspection of periodic moving objects | |
JP5539719B2 (en) | Image forming system | |
Koken et al. | Aperture weighted cardiac reconstruction for cone-beam CT | |
US20100232663A1 (en) | Computed tomography reconstruction for two tilted circles | |
JP2008539931A (en) | Reconstruction method of helical cone beam CT | |
US7426257B2 (en) | Computer tomography method for a periodically moving object | |
JP2002034970A (en) | Method and device for spiral reconstitution in multi-slice ct scan | |
US20060120586A1 (en) | Method and system for extracting information about the cardiac cycle from CT projection data | |
US6999550B2 (en) | Method and apparatus for obtaining data for reconstructing images of an object | |
JP2011502683A (en) | Imaging apparatus, imaging method, and computer program for determining image of region of interest | |
US7292717B2 (en) | Computer tomography method using redundant measured values | |
US7436925B2 (en) | Computer tomography method for objects moving periodically | |
Hoppe et al. | Accurate image reconstruction using real C-arm data from a Circle-plus-arc trajectory | |
Manzke et al. | Image reconstruction with multi-slice CT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated states |
Kind code of ref document: A2 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A2 Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2004744809 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2007019776 Country of ref document: US Ref document number: 10570088 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2006525236 Country of ref document: JP |
|
WWP | Wipo information: published in national office |
Ref document number: 2004744809 Country of ref document: EP |
|
CFP | Corrected version of a pamphlet front page | ||
CR1 | Correction of entry in section i |
Free format text: IN PCT GAZETTE 11/2005 UNDER (71) REPLACE "FOR AE, AG, AL...ZM, ZW ONLY" BY "FOR ALL DESIGNATED STATES EXCEPT DE, US" |
|
WWP | Wipo information: published in national office |
Ref document number: 10570088 Country of ref document: US |