Patent Application
20100286928
This application claims the benefit of DE 10 2009 020 400.8 filed May 8, 2009, which is hereby incorporated by reference.
The present embodiments relate to a method and a device for determining attenuation coefficients for an object using a movable X-ray source and a detector.
X-ray methods are standard techniques in medical engineering. With simple X-ray photographs, X-ray radiation is transmitted through an object that is to be examined and subsequently recorded by a detector. The recording or projection represents information about the attenuation of the transmitted X-ray beams on a path through the object. The attenuation of the X-ray radiation is dependent on the density of the object parts penetrated by the radiation. The density in turn yields information about the properties of the object, which are typically presented in visual form for diagnostic purposes. In the case of an X-ray photograph, the intensity registered by the detector is dependent on the overall composition of the object along the path traveled by the X-ray beam (i.e., information integrated over a distance is available). Consequently, attenuation coefficients of the object as a function of all three space coordinates are not obtained from a single X-ray photograph. For a three-dimensional image, therefore, a plurality of X-ray photographs are taken from different recording positions, and a three-dimensional image is reconstructed from the plurality of X-ray photographs. In medical engineering, this approach is adopted in computed tomography (CT). Within the framework of computed tomography, the X-ray source and X-ray detector travel along a path or trajectory, and recordings are taken along the trajectory. The recordings are used to reconstruct a three-dimensional image of attenuation coefficients, which relate to the density.
Image reconstruction is a complex, mathematical problem. Two groups of methods have become established for three-dimensional image construction: approximate and exact methods. In this context, theoretically exact methods are methods which, mathematically, include no approximations; the numerical conversion may introduce errors. The approximate methods (e.g., the Feldkamp algorithm) initially had the advantage of the significantly less complex numerical conversion. Thanks to skillful formulations of the mathematical problem, which have been proposed in the last several years, there is now available a theoretically exact formulation, which can be converted numerically with a realistic amount of effort. This is described, for example, in U.S. Pat. No. 6,771,733 B2. U.S. Pat. No. 6,771,733 B2 discloses a reconstruction formula (formula 10), which is well suited to the numerical conversion of an exact method. For implementation, this reconstruction formula generally makes use of a further transformation according to the path or trajectory used. The corresponding formula for a spiral path is expressed, for example, in the cited publication as formula 29.
In respect of the numerical conversion, however, difficulties continue to exist. One challenge that remains is finding a reasonable compromise between the number of projections recorded and the image quality, the image quality generally being higher, the more images that are recorded. However, it is also desirable to limit the number of recordings made in order to limit the exposure to radiation of patients being examined.
The present embodiments may obviate one or more of the drawbacks or limitations in the related art. For example, in one embodiment, the number of projections in approximate and exact reconstruction methods may be reduced without a deterioration in image quality.
In the present embodiments, attenuation coefficients for an object are determined using a movable X-ray source and a detector. In this context, the term detector is to be understood in a broad sense (e.g., the detector may be a detector system having a suitable arrangement of detectors or detector surfaces). Either both the X-ray source and the detector may be movable (e.g., CT applications) or only the X-ray source may be movable (e.g., tomosynthesis in mammography). In one embodiment, the attenuation coefficients are determined using projections (e.g., recordings of X-ray radiation transmitted through the object), a path or trajectory being specified for the movable X-ray source. Positions are specified on the trajectory (e.g., equidistantly) for the numeric calculation of a derivative of projections recorded by the detector. A plurality of scanning positions (e.g., 2 or 3) is specified for each of the specified positions. The movable X-ray source travels along the trajectory, a projection being recorded for each of the specified scanning positions. In one embodiment, the trajectory may be followed at a constant speed. At a constant speed, the system used for the recordings is subject to comparatively low mechanical stresses. A constant speed may be used for the duration of the treatment, because the trajectory can be completed more quickly due to the lower mechanical loads. Accordingly, different time intervals occur between recordings of projections (e.g., shorter time intervals in the case of adjacent recordings that are assigned to the same plurality of specified scanning positions than if the adjacent recordings do not belong to the same plurality of scanning positions). The projections obtained are used for the numeric calculation of the projection derivative in relation to the trajectory, the projections recorded for the associated plurality of scanning positions being used for each of the positions. The calculated projection derivatives are used for the reconstruction using a theoretically exact or approximate rule or formula for determining the attenuation coefficients for the object. In one embodiment, further acts may be provided (e.g., for preprocessing the projections prior to the reconstruction using exact or approximate formulas). Examples of such further acts would be, for example, resorting (e.g., binning) and filtering. The acts in the numeric derivative calculation and reconstruction may also be linked to one another in such a way that a complete calculation of all derivatives is not performed before the start of the back-projection. These are measures known to the person of skill in the art, for which the person of skill in the art will make provision as a matter of routine as an embodiment of the teaching according to the present embodiments. In one embodiment, a usable theoretically exact formula or mathematical rule can be obtained, for example, by a transformation of formula 10 from U.S. Pat. No. 6,771,733 B2 for a specific trajectory (e.g., spiral, circle, or line). However, the present embodiments may also be used with approximate methods. For example, tomosynthesis is usually not exact by construction (due to the incomplete scanning path), but reconstruction methods that make use of a derivative along the trajectory are used.
Both in the case of some approximate methods and in the case of the known theoretically exact formulas that are suitable for a numeric conversion for the reconstruction for determining attenuation coefficients in three dimensions, a derivative of projections appears according to a parameter, which parameterizes the trajectory followed by the recording source. This derivative plays an important role for the numeric conversion of the formulas for the back-projection. After positions for determining derivatives of the projections recorded by the detector are specified, and scanning positions are assigned to the specified positions, the quality of the numeric derivative may be optimized by suitable specification of the scanning positions. As a result, the same quality may be achieved with fewer recordings overall. As a consequence, images may be recorded that have the same quality yet expose the patient to a lower dose of radiation.
The present embodiments also include a device through which the methods described above are implemented. The device may use software, hardware, firmware or a combination thereof in order to realize the present embodiments.
The individual acts of the methods described above may also be realized using corresponding functional software modules.
During the scanning of the patient P, dose-dependent signals acquired by the detector 5 are transmitted to a computing unit 7 via a data/control line 6. From the measured raw data, the spatial structure of a scanned region of the patient P is calculated in terms of absorption values in a known manner with the aid of known methods that are stored in program modules P1 to Pn. In one embodiment, a theoretically exact reconstruction method is applied.
In other respects, the CT scanner is likewise operated and controlled using the computing unit 7 and a keyboard 9. The computing unit 7 may be a processor, computer, server or other device, configured with data stored in non-transitory computer readable media, for performing the acts described herein. Computed data may be output via a monitor 8 or a printer (not shown).
The majority of the efficiently usable approximate and exact reconstruction methods (e.g. those described in U.S. Pat. No. 6,771,733 B2) require the calculation of the derivative along the scanning path.
In the case of exact three-dimensional reconstruction methods, there are essentially three methods that are used for calculating the derivative along the scanning path: the “direct scheme”; the “chain rule scheme;” and the “new scheme” (cf. F. Noo, S. Hoppe, F. Dennerlein, G. Lauritsch, and J. Hornegger. “A New Scheme for View-Dependent Data Differentiation in Fan-Beam and Cone-Beam Computed Tomography.” Physics in Medicine and Biology 52.17 (2007): 5393-5414.). The direct scheme is easy to implement and yields good results provided the scanning rate along the scanning path is high. It has been demonstrated, however, that if the scanning rate is too low the direct scheme fails to deliver usable results. The chain rule scheme generally copes better with a lower scanning rate. This does not apply unreservedly, however, since the chain rule scheme behaves like the direct scheme with certain types of scanning paths, such as, for example, a linear scanning path. Consequently, the chain rule scheme cannot be used without restriction when the scanning rate is too low. The new scheme was developed in order to solve this problem, so that usable results can still be achieved even with lower scanning rates. However, it is likely that a higher scanning rate will lead to more precise and consequently better results, even if a higher scanning rate is not essential for the use of the new scheme.
This is the starting point for the present embodiments. A departure is made from constant scanning in order to enable the number of scanning points to be reduced. Points on a trajectory are specified (e.g., equidistantly). A plurality of scanning points are specified for each of the specified points on the trajectory. The point of the trajectory for which the derivative is calculated does not necessarily coincide with one of the plurality of scanning points. This may be the case for an uneven number of scanning points per plurality of scanning points.
In contrast to the conventional scheme, in which only the interval between scanning points can be varied, it is possible in one embodiment to optimize the ratio of the number of scanning points to image quality by setting at least three variables such as, for example, a spacing of the points for which the derivative is calculated, a number of scanning points per derivative calculation and a spacing of the scanning points.
Curve A of
The invention is not limited to the cases described in the present embodiments. In particular, other recording protocols than are shown in
While the present invention has been described above by reference to various embodiments, it should be understood that many changes and modifications can be made to the described embodiments. It is therefore intended that the foregoing description be regarded as illustrative rather than limiting, and that it be understood that all equivalents and/or combinations of embodiments are intended to be included in this description.
| Number | Date | Country | Kind |
|---|---|---|---|
| 10 2009 020 400.8 | May 2009 | DE | national |
