Patent Application
20070025499
The present application hereby claims priority under 35 U.S.C. §119 on German patent application number DE 10 2005 034 876.9 filed Jul. 26, 2005, the entire contents of which is hereby incorporated herein by reference.
The invention generally relates to a method for compiling computer tomographic representations. For example it may relate to one using a CT system with at least two angularly offset ray sources, wherein at least a first ray cone with a larger fan angle and a second ray cone with a smaller fan angle scan an object circularly or spirally and the first ray cone generates a first dataset A because of its absorption in the object and the second ray cone generates a dataset B, the B data of the smaller ray cone being supplemented with other data at the edge to give an expanded dataset B+ for reconstruction of the CT representation, the expanded B+data of the first ray cone and the A data of the larger ray cone is subjected to a convolution operation to give B+′ and A′ data, and a back projection being carried out from the convoluted B+′ and A′ data.
To compile computer tomographic representations using a CT with a plurality of angularly offset ray sources, it is generally known to combine data from the individual ray sources and thereby carry out the reconstruction, i.e. a convolution of the data with subsequent back projection, for 2D or 3D reconstruction.
For reconstruction from a dataset of a ray cone which does not fully scan the object, it is furthermore known to supplement the peripherally existing sharp boundaries with corresponding extrapolations of data in order to reduce artifacts which have been created by a sharp delimitation of the data. In this regard, reference is made for example to document DE 198 54 917 A1 which describes such an extrapolation of peripheral data, albeit in a CT with a single-tube/single-detector system.
It has been found that both a simple data supplement in a CT system with a plurality of ray cones, respectively from the other ray cone, can lead to artifacts in the reconstruction of computer tomographic representations.
A method is provided which, in at least one embodiment, allows a data supplement between ray cones in a CT with a plurality of offset ray cones but reduces or even avoids the artifacts occurring in the prior art.
In at least one embodiment, the inventors have discovered that it can be advantageous for the reconstruction to carry out the employed convolution with datasets which are respectively supplemented with data of another ray cone, but that the supplementary data should at least partially be removed again after the convolution and before the actual back projection.
In particular, such data removal is advantageous whenever the supplementary data come from complementary data of the other respective ray cone since, with these data, the filter direction which is normally carried out in the scan direction is incorrect for these complementary rays.
According to this basic concept, in at least one embodiment the inventors propose a method for compiling computer tomographic representations using a CT system with at least two angularly offset ray sources, wherein at least a first ray cone with a larger fan angle and a second ray cone with a smaller fan angle scan an object circularly or spirally and the first ray cone generates a first dataset A because of its absorption in the object and the second ray cone generates a dataset B, the B data of the smaller ray cone being supplemented with other data at the edge to give an expanded dataset B+ for reconstruction of the CT representation, the expanded B+ data of the first ray cone and the A data of the larger ray cone is subjected to a convolution operation to give B+′ and A′ data, and a back projection being carried out from the convoluted B+′ and A′ data, wherein the dataset B is supplemented with data of the dataset A and supplementary data are removed from the dataset B+′ after the convolution but before the back projection.
This type of data supplement and subsequent data reduction reduces artifacts otherwise occurring in the CT representation.
In principle, either only data of the dataset A with the same ray direction may be used for supplementing the dataset B or data with the same ray direction and with a complementary ray direction may also be used. In respect of the data to be removed again after the convolution and before the back projection, on the one hand all the supplementary data may be removed from the dataset B+′, or alternatively only the supplementary which come from complementary rays of the dataset A may be removed from the dataset B+′.
The reconstruction carried out may be a voxel-wise 3D reconstruction, or a planar 2D reconstruction may be carried out, in which case the data removed from the dataset B+′ must be replaced by data from the convoluted dataset A′ for the 2D reconstruction.
A different inclination of the rays relative to the z axis may furthermore be neglected when supplementing the dataset B with data of the dataset A. It is likewise possible to use interpolated data from the dataset A for supplementing the dataset B.
It is furthermore possible for weighting to achieve a smooth transition to be carried out in the supplemented dataset B+ at the transition between the data of the dataset B and the supplementary data from the dataset A.
According to at least one embodiment of the invention, different procedures may then be carried out with the result of the convolution of the original datasets A and B, either by recombining these datasets to give a common dataset before the back projection or by carrying out one back projection per dataset and generating the common image with this result.
In the latter case, a back projection may respectively be carried out separately with the dataset A′ and with the dataset B+′ (−), optionally with a further data supplement, and a common image may be generated by weighting. Transitional weighting can in this case advantageously be carried out in the transition range of the measurement field of the smaller ray cone B.
Corresponding to the former case, a common 2D back projection is thus carried out with the A′ data and the − optionally postprocessed—B+′ (−) data, to which end the two datasets may be combined by weighting to give a dataset AUB before the back projections. Here, there is also the possibility of transitional weighting in the transition range of the data of the measurement field of the smaller ray cone B.
For the case of 3D reconstruction, a common 3D back projection may be carried out voxel-wise with the A′ data and the—optionally postprocessed—B+′ (−) data, and the two datasets may be combined by weighting to give a new dataset before the voxel-wise back projections.
According to at least one embodiment of the method described above, the Inventors also propose to improve a computer tomography system known per se for the tomographic representation of an object with at least two angularly offset ray sources, which form a first ray cone with a larger fan angle and a second ray cone with a smaller fan angle, scan an object circularly or spirally and generate detector output data, the first ray cone generating a first dataset A because of its absorption in the object and the second ray cone generating a dataset B, which are processed with the aid of stored computer programs or program modules in a computation and control unit to give tomographic representations of the object, so that the stored computer programs or program modules also comprise program code which implements the method steps described above.
The invention will be explained in more detail below with reference to the example embodiments with the aid of the figures, only the features necessary for understanding the example embodiments of the invention being represented, the following reference numerals being used in the figures: 1: CT system with two X-ray tubes; 2: first tube; 3: first detector; 4: second tube; 5: second detector; 6: CT housing; 7: opening in the CT system; 8: displaceable patient support; 9: system axis, 10: control and computation unit; 11: large ray cone of the A system; 12: large measurement field of the A system; 13: small ray cone of the B system; 14: small measurement field of the B system; 15: image; A: data of the A system; B: data of the B system; DA detector of the A system; DB detector of the B system; FA focus of the A system; FB focus of the B system; P: patient; Prg1-Prgn: computer programs; α: projection angle; βA: fan angle of the A system; βB: fan angle of the B system; I-VI: method steps.
In detail:
A patient P, who can be moved along the system axis 9 through an opening 7 at the centre of rotation of the tube/detector combination, lies on a displaceable patient support 8, the actual scanning taking place either sequentially and circularly or, with a continuous forward movement of the patient, in a spiral-shaped profile relative to the patient.
The computer tomographic system 1 is controlled by a control and computation unit 10, in which the computer programs Prg1 to Prgn employed to carry out the method are stored and can be used in operation.
A section through the gantry of such a system of
The first ray cone 11 has a fan angle βA which is designed to be substantially greater than the fan angle βB of the second ray cone 13. The measurement field 12 associated with the focus FA is correspondingly much larger than the measurement field 14 of the focus FB.
If objects which exceed the smaller measurement field 14 of the B system are studied with such a CT, then strong edge effects occur when convoluting the data of the B system for the reconstruction because the attenuation at the edge has not decreased to zero, but instead forms a sharp edge. As already mentioned above, a data supplement must in any case be carried out before the convolution. In the case according to the invention, this data supplement is carried out with data which come from the other ray cone 11, i.e. the A system.
Here, for each missing ray of the B detector, it is possible to look for a ray of the A detector so that, on the one hand, the direction of the ray coincides or so that a complementary ray then with the correct direction is used and, on the other hand, the z position of the rays must coincide as well as possible with the z positions of the rays from the B detector. Since the z position for an inclined ray is not the same over the entire ray length, the z position of the rays at a point where their distance from the z axis is least may for example be used for comparison here.
The different inclination of the rays relative to the z axis may be neglected in this method, although it should in any case be pointed out that the data supplement in the method presented here is in principle carried out approximately. Thus, the data used may be correspondingly interpolated if there are no exact supplementary rays.
In order to avoid hard transitions between the two data ranges, i.e. the original data range and the supplementary data range, it may be expedient to use weighting with the aid of a smooth transition function, for example a cosine function.
In principle, it is possible to apply the method according to at least one embodiment of the invention both to 2D reconstructions and to 3D reconstructions, in the case of a 2D reconstruction it being necessary for the data removed after the convolution to be replaced again, while in 3D reconstruction it is also possible to carry out the reformatting without supplementary data since compensation is provided here via the normalization which takes place anyway.
Owing to the smaller measurement field in column <B>, originally there are in fact only restricted data B there between the dotted lines which show the limits of the B system. In the first step, the data B not present in the B system are replaced by data A of the A system—represented by the two arrows between columns—and the dataset B+ is formed. A convolution is subsequently carried out for both sides, the result of which is represented by the datasets A′ and B+′ in row III. If the data lying outside the actual measurement field of the B system are now used for the back projection, then the resulting CT image would comprise increased artifacts—as already explained above. It is therefore more favorable for the convoluted data lying outside the B range now to be cut out again, so as to give the data B+′ (−) represented in row IV. According to at least one embodiment of the invention, the already convoluted data A′ of the A system can now be used in the range outside the B system in order thereby to supplement the data B+′ (−), so as to give the dataset B+′ (−) +, of row V which, for example, may be used for back projection in the 2D reconstruction.
Since both the convoluted dataset A′ from the detector data of the further detector DA, i.e. of the more widely fanned ray cone 11, and of the convoluted and supplemented dataset B+′ (−) are now available, it is possible to proceed with the back projection in different ways.
In one variant shown in
Another variant of a 2D reconstruction according to at least one embodiment of the invention is shown in
A procedure similar to that described here for 2D reconstruction can also be applied to voxel-wise 3D reconstruction. This method according to at least one embodiment of the invention may be used in both reconstruction methods, which are based on a scan with at least two ray cones having different fan widths and scan fields of different sizes.
It is to be understood that the features of the embodiments of the invention as mentioned above may be used not only in the respectively indicated combination, but also in other combinations or individually, without departing from the scope of the invention.
Further, elements and/or features of different example embodiments may be combined with each other and/or substituted for each other within the scope of this disclosure and appended claims.
Still further, any one of the above-described and other example features of the present invention may be embodied in the form of an apparatus, method, system, computer program and computer program product. For example, of the aforementioned methods may be embodied in the form of a system or device, including, but not limited to, any of the structure for performing the methodology illustrated in the drawings.
Even further, any of the aforementioned methods may be embodied in the form of a program. The program may be stored on a computer readable media and is adapted to perform any one of the aforementioned methods when run on a computer device (a device including a processor). Thus, the storage medium or computer readable medium, is adapted to store information and is adapted to interact with a data processing facility or computer device to perform the method of any of the above mentioned embodiments.
The storage medium may be a built-in medium installed inside a computer device main body or a removable medium arranged so that it can be separated from the computer device main body. Examples of the built-in medium include, but are not limited to, rewriteable non-volatile memories, such as ROMs and flash memories, and hard disks. Examples of the removable medium include, but are not limited to, optical storage media such as CD-ROMs and DVDs; magneto-optical storage media, such as MOs; magnetism storage media, including but not limited to floppy disks (trademark), cassette tapes, and removable hard disks; media with a built-in rewriteable non-volatile memory, including but not limited to memory cards; and media with a built-in ROM, including but not limited to ROM cassettes; etc. Furthermore, various information regarding stored images, for example, property information, may be stored in any other form, or it may be provided in other ways.
Example embodiments being thus described, it will be obvious that the same may be varied in many ways. Such variations are not to be regarded as a departure from the spirit and scope of the present invention, and all such modifications as would be obvious to one skilled in the art are intended to be included within the scope of the following claims.
| Number | Date | Country | Kind |
|---|---|---|---|
| 10 2005 034 876.9 | Jul 2005 | DE | national |
