Embodiments of the present invention relate to possibilities of creating a CT reconstruction of an object that has a high resolution in an object region that is of particular interest, said resolution being lower in other object regions which adjoin or comprise said object region of interest.
In the non-destructive testing of objects, or in the non-destructive examination of patients by means of X-ray computed tomography, the achievable resolution is limited essentially by two factors. They include, firstly, the finite expansion of the radiation source, i.e. of that area from which the radiation used for tomography is emitted (for example the focal spot of an X-ray tube), and, secondly, the finite expansion of the detector elements. With a finite expansion of one of these two components, the idealized way of looking at an idealized “X-ray beam” of an infinitesimal expansion, which perspective underlies many reconstruction algorithms, is no longer fulfilled by the object to be examined.
The pixels used for sampling the individual projection images in an X-ray-sensitive detector, such as, for example, those of an electronic flat-panel image converter (for example of a CCD comprising a radiation-converting coating, or of a directly converting semiconductor detector, or the like) naturally have an intermediate distance that is finite in each case, whereby the resolution is also limited.
Since this typically corresponds to the practical circumstances, it shall be assumed below, without prejudice to the generality, that the expansion of the radiation source is smaller than the limitation of the resolution that is caused by the finite distance of the detector pixels.
In order to be able to reconstruct in an artefact-free manner, by means of standard algorithms, a series of X-ray projection images, i.e. a plurality of recordings, or X-rays, obtained from different perspectives by means of an extensive one- or two-dimensional detector, the object may be fully contained within the horizontal extensions in each of the recorded projections (if the object projects beyond the detector at the top and/or at the bottom, this will not lead to any artefacts in the reconstruction). In other words, for three-dimensional reconstruction (CT reconstruction), the object may be fully imaged in the horizontal extension on each two-dimensional shadow image (projection). Therefore, the horizontal extension describes, in this context, that orientation of the object relative to which the perspective is changed by the rotation. In the vertical direction which is perpendicular thereto, the object may be imaged in an incomplete manner. With a point-shaped radiation source and detector having finite dimensions, the object consequently cannot be positioned at any distance from the detector, since otherwise the geometrical projection of the object would protrude beyond the detector. The more severe the violation of the condition of fully imaging the object in each projection, the more intense the interferences in the CT sectional images or in the three-dimensional reconstructions of the object examined will be, said interferences being caused by the image reconstruction algorithm. The share of artefacts in the reconstructed layers or models thus increasingly impairs the diagnostic conclusion that can be drawn from the images, until, in extreme cases, said images contain no more useful information.
On the one hand, the resolution is thus limited by the apparatus used and, in particular, by the resolution of the sensor used. On the other hand, the condition of completely imaging the object in the individual projections, which condition is set up by the image reconstruction algorithm, limits the spatial resolution since, in this manner, the optical magnification of the object to be examined—caused by geometric variations of the distances between the detector and the object—on the detector is capped.
This is relevant particularly to geometrically expanded objects, such as in the non-destructive testing of the material of motors or similarly large components, wherein fine details can no longer be resolved with the available sensor resolution when, as was described above, the entire object is to be imaged, in each projection view, onto the sensor having a finite size.
The possibility of obtaining a complete set of projection data in a magnified view by not imaging the complete object per projection is limited since, in this case, the image reconstruction algorithm produces considerable artefacts. This is due to the fact that the image reconstruction algorithm depends on containing complete information about the object from all perspectives. However, this is not the case if the edge of object is not imaged in each the projections. However, this edge contributes, with a view rotated by 90°, to the absorption coefficient and, thus, to the entire X-ray absorption, so that the CT reconstruction algorithm produces considerable artefacts in reconstruction, whose sizes increase as the proportion of the region that is not imaged at the edge of the object increases.
Therefore, it is useful to be able to examine partial regions of an object by means of CT methods in a manner that is complete and free from artefacts, even with objects that have large geometrical expansions.
Some embodiments of the present invention enable highly detailed CT imaging of an object region of interest within an object in that, initially, a first projection data set of the object is produced at a first, low resolution, and in that subsequently, a second projection data set, which comprises only the object region of interest, is detected at a second, higher resolution. By combining the first and second projection data sets thus obtained in accordance with a combination rule, a CT reconstruction of the first region of the object comprising the first resolution may be obtained, it being possible to reconstruct that object region of interest that is arranged within or at the edge of the first region at the second, higher resolution. In this context it is possible, in particular, to reconstruct, by detecting the object twice, a complete data set that remains free from image artefacts within the object region of interest.
In accordance with some embodiments, the first and second regions are recorded simultaneously in that, for example, two detectors of different geometric expansions and spatial resolutions are used. The high-resolution detector may be arranged upstream from the low-resolution detector, for example, it being possible to substitute those regions of the low-resolution detector that are shadowed by the high-resolution detector by using the image data of the high-resolution detector to supplement the missing data in the image of the low-resolution detector.
In some embodiments, the object is moved, or the distance between the detector and the object is varied, so that the object may be recorded by means of a detector with different magnification factors, so as to obtain, by subsequently combining the two projection data sets, the representation of the object region of interest (ROI=region of interest), said representation having a high resolution and being free from artefacts.
In addition, in accordance with some embodiments, the concept of duplicate recording or of duplicate production of projection data sets of different resolutions may also be implemented in connection with different X-ray sources or with X-radiation of different wavelengths. Since the absorption cross-section of X-radiation is dependent on the material and energy, a maximum image contrast may be achieved, depending on the object investigated, for different energies of the X-radiation used. Thus, if the object to be examined contains different materials, different X-ray energies may result in that the achievable resolution is identically high for both materials despite the different absorption properties.
In the context of the present application or invention, the term resolution is therefore not only to be understood as spatial resolution, but as a term describing how much information may be obtained from the data or the individual projection recordings. An illustrative example of this may be the fact that, even though a detector with randomly high spatial resolution is used, the information obtained is only approximately zero if X-radiation of such low energy is used that it is almost entirely absorbed by the object. In the extreme case of very high X-ray energy, where hardly any absorption takes place, this is also true, of course.
If two projection recordings with different resolutions are produced, different possibilities or combination rules will result in accordance with which both projection recordings, or the different projection data sets consisting of several projection recordings, may be combined such that an image which is free from artefacts and has the high resolution of the second projection recording will result in the object region of interest.
A second possibility of combination and an associated second combination rule result when the data in the second projection recording are processed such that only information about the high-frequency image portions within the object region of interest is taken from said data. The low-frequency image information (low spatial frequencies) is already contained in the first projection recordings. If this pre-processing of the second projection data set is performed such that all of the image information or almost all of the image information above a cutoff frequency, which is defined by the maximally imageable position information in the lower-resolution images of the first projection data set, is contained therein, both projection data sets may be combined, by simply adding the corresponding recordings, such that the resolution of the second projection recordings is available within the object region of interest without artefacts or similar image interferences resulting from the separate pre-processing.
In some embodiments, pre-transformation is used for the second projection recordings so as to extract those higher-frequency image portions associated with the object region of interest which correspond, in a retro-projected image, to a wavelet decomposition of the recording. Due to the positional- and frequency-locality property of the wavelet base functions, the above criterion can be ensured. This means that the representation or wavelet decomposition can be selected such that the wavelets contribute to the reconstruction only within the object region of interest, but describe, within said object region of interest, the relevant information, namely the high-frequency image portions. In this manner, it is possible to obtain combination recordings which are provided, within the object region of interest, by the desired high resolution in an artefact-free manner. In particular the high-frequency and low-frequency image portions may be processed and reconstructed separately, it being possible to combine the representations after the reconstruction by means of addition.
It is readily possible for a person skilled in the art to draw, from the wavelet representation of the reconstructed recordings that is to be achieved, the conclusion as to which filters and/or operations of a pre-transformation are to be applied to the individual projection recordings.
According to an embodiment, a method of producing a CT reconstruction of an object including a high-resolution object region of interest may have the steps of: producing a first projection data set of a first region of the object that encloses the object region of interest and includes at least one projection recording of a first resolution; producing a second projection data set for the object region of interest that includes at least one second projection recording of a second, higher resolution; pre-transforming the second projection data set to acquire a pre-transformed second projection data set including, within the object region of interest, only image frequencies above a predetermined cutoff frequency, which is equal to or smaller than an upper imaging frequency, defined by the first resolution, within the first projection data set; and combining the first and the pre-transformed second projection data sets in accordance with a combination rule so as to acquire a CT reconstruction of the first region of the object including the first resolution, and of the object region of interest including the second, higher resolution.
According to another embodiment, a device for producing a CT reconstruction of an object including a high-resolution object region of interest may have: a detector configured to produce a first projection data set of a first region of the object that encloses the object region of interest and includes at least one projection recording of a first resolution; produce a second projection data set of an object region of interest including at least a second projection recording of a second, higher resolution; and a combiner configured to pre-transform the second projection data set to acquire a pre-transformed second projection data set including, within the object region of interest, only image frequencies above a predetermined cutoff frequency, which is equal to or smaller than an upper imaging frequency, defined by the first resolution, within the first projection data set, and to combine the first and the pre-transformed second projection data sets in accordance with a combination rule so as to acquire a CT reconstruction of the first region of the object including the first resolution, and of the object region of interest including the second, higher resolution.
According to another embodiment, a computer program has a program code for performing the method of producing a CT reconstruction of an object including a high-resolution object region of interest, wherein the method may have the steps of: producing a first projection data set of a first region of the object that encloses the object region of interest and includes at least one projection recording of a first resolution; producing a second projection data set for the object region of interest that includes at least one second projection recording of a second, higher resolution; pre-transforming the second projection data set to acquire a pre-transformed second projection data set including, within the object region of interest, only image frequencies above a predetermined cutoff frequency, which is equal to or smaller than an upper imaging frequency, defined by the first resolution, within the first projection data set; and combining the first and the pre-transformed second projection data sets in accordance with a combination rule so as to acquire a CT reconstruction of the first region of the object including the first resolution, and of the object region of interest including the second, higher resolution, when the program runs on a computer.
Other elements, features, steps, characteristics and advantages of the present invention will become more apparent from the following detailed description of the preferred embodiments with reference to the attached drawings.
Embodiments of the present invention will be detailed subsequently referring to the appended drawings, in which:
As was already described above, for the application of standard CT image reconstruction methods, the object (in this case replaced by the spatial expansion W of the measuring field 4) should be fully contained, in the horizontal expansion, within the individual projection recordings of the detector 6. A projection data set consisting of several projection recordings is typically obtained in that either the object and/or the measuring field 4 is rotated relative to the arrangement of the X-ray source 2 and the detector 6, it being possible, in principle, to freely select the angle increments Δα1 of the rotation, as will be explained below. Alternatively, the X-ray source 2 and the detector 6 may be rotated about the center of rotation (center of the measuring field 4), as, for example, in a computer tomograph for application in human medicine. If the above requirement of complete imaging of the object in the projections is met, this will mean that an effective spatial resolution Δeffective will depend on the size of the object and, thus, on the maximally applicable magnification, since effective sampling (within the object) by a detector (pixel distance Δ) decreases as the magnification factor increases:
Δeffective=Δ/M (1)
The magnification factor M in this context is defined by the distances of the axis of rotation (which typically has the object arranged centrally thereon) from the detector and/or from the source
wherein b is the distance of the source from the detector, and a is the distance of the source from the center of rotation.
Let W be the diameter of the measuring field and, thus, the maximum object size that guarantees that each projection is fully imaged to the detector. The upper limit for the magnification results from this as
wherein D is the width of the input face of the detector, which is typically specified by the design of the flat-panel image converter installed within the CT system and which cannot be adapted to different problems.
For a flat-panel image converter with ND detector elements per column or row (the validity of this consideration not being limited to square detectors), the following applies:
This means that in the event of complete imaging of the object in the above-described sense, the resolution will be effectively limited by W/ND even with a sufficiently small focal spot. For example, a cylindrical casting having a diameter of 200 mm can be imaged, while using a 1,024-column detector, no better than at a resolution of approx. 195 μm in the sectional planes.
On the basis of the preceding considerations, it is impossible to select a magnification factor that is sufficiently large that, for each projection, only an inner sub-region of an object is recorded within the measuring field 4, but is recorded with a larger magnification to make up for it, since this would lead to intense image artefacts which in extreme cases would render useful interpretation of the recordings impossible.
Generally, within the context of the present description, elements that are identical or similar in function are provided with the same reference numerals, it being possible to interchange their descriptions with regard to the individual figures.
a shows the situation of how a projection data set of the entire object (measuring field) 4 is initially created by means of one and the same detector 6, said projection data set having a first resolution determined by the above-described geometric conditions.
b illustrates how a second projection data set may be obtained, using the same detector 6, for an object region of interest, which here is arranged at the center of the measuring field 4 without prejudice to the generality. Alternative embodiments or alternative applications of the concept may naturally also image, at a high resolution, such object regions that are not located at the center of rotation.
As is depicted in
The detector is connected to a combination means 8 that enables combining the projection data sets, as will be set forth later on.
In other words,
a and 2b therefore show two possible operating modes of a detection means 7, which enables creating a first projection data set having a first resolution, and a second projection data set having a second, higher resolution.
The combination means 8 enables combining the first and second projection data sets in accordance with a combination rule so as to obtain a CT reconstruction of the first region 4 of the object with the first resolution, and of the object region 10 of interest with the second, higher resolution.
Both data sets may be combined by suitable means such that within the object region of interest (ROI), in this case, for example, the central region, the desired resolution is achieved without the image quality in the volume data set, which represents the entire object, comprising, after the CT reconstruction, artefacts that result from the incompleteness of the data.
Embodiments of suitable combinations or suitable combination rules (means for combining) will be explained below in more detail with reference to
In the approach suggested in
This means that with a flat-panel image converter of a specific size, said converter being predefined by the CT system available, the user-selectable parameter regarding the number of angular positions on the full circle from which projection images are acquired may be selected.
With an assumed angle increment of Δα, radial sampling at the distance RM from the center of rotation 12 results directly from the arc length:
Δr=RM·Δα (6)
The angular increment follows directly from the number of angular positions on the full circle:
If radial sampling at the edge of the measuring field (the maximum distance of two rays in successive projections) is equated with the sampling distance Δeffective within the projection from (1), the following results for the number of projection angles:
In the simplest case, the number of the acquired projections will be identical in both measurements described above. However, in order to reduce the overall time that may be taken to perform both measurements and, thus, to increase the test efficiency of a computer tomography system designed in accordance with the present invention, the number of projections may be reduced, in the first measurement (which fully images the object), to such an extent that the minimum requirements for angle sampling are met only within the sub-region to be reconstructed at a high resolution.
This means that the fact that there is an interest only in the image information within the object region of interest can additionally be accounted for in that, even in the lower-resolution recording, the information having the resolution that is maximally possible in view of the reduced resolution is recorded only for the object region of interest.
The first detector 20, which has a lower spatial resolution, may be arranged, for example, upstream or downstream from the second detector 22 having the higher spatial resolution. If useful for the image reconstruction or for combining the images, the image information which is missing for the second detector 22 by shadowing on the projection recordings of the first detector 20, may be readily substituted by the image information of the second detector 22. An example of such a substitution would be, for example, to sum the intensities of the high-resolution pixels in a weighted manner or to combine them in any other manner so as to achieve an intensity value corresponding to a detector pixel of the first detector 20.
The arrangement shown here by way of example is not limited to square geometries; rather, any detector shapes may be randomly combined with one another. Also, the second detector 22 may be arranged such that it does not fully cover an even number of pixels of the first detector 20, as is indicated here for simplicity's sake.
In addition, in accordance with some embodiments, the second detector 22 may be arranged to be movable relative to the first detector 20, as is indicated here, for example, by a directional arrow 24, so that the second detector 22 may be moved in two dimensions (x and y directions) relative to the first detector 20. With setups that are not fully rotationally symmetric, i.e. wherein the second detector 22—as is already shown here—is not arranged at the center of the first detector 20, this may be used for having the second detector follow the region of interest or the object region of interest upon rotation of the object. In addition, the second detector 22 may also be movable in three dimensions so as to be able to vary, additionally, the effective resolution of the second detector 22 independently of the first detector 20.
Additionally, a second detector arranged upstream from a first detector acts as a prefilter for X-radiation, and ensures, at a constant tube voltage or X-ray energy, that the beam qualities and, thus, the contrast with which different materials are imaged differ between the two projections of different spatial resolutions. In this manner, a two-spectra data set may additionally be recorded which can be evaluated using the known methods of material analysis.
Alternatively or additionally, as is indicated by the supplementation step 28, a second intermediate image 30 may be produced from the second projection recording in that image information outside the object region of interest recorded directly by means of the second detector 22 is supplemented with a resolution corresponding to the resolution of the second detector 22. The image information, useful for supplementation, for generating the second intermediate image may be obtained from the image information of the first projection recording 20. In both cases it is useful, in order to produce the intermediate image (third intermediate image), which may be handed over to a standard CT algorithm, to suitably combine the first projection recording 20 and the second projection recording 22. If the first intermediate image 26 is directly combined with the second projection recording 22, the above may be effected by a summation 32, for example. Alternatively, the image information that is already contained within the second projection recording 22 may be fully removed from the interpolated intermediate image representation 26 of the first projection recording prior to summation, so that there will be no “overexposure”, in terms of intensity, of the image region in question.
In other words,
In this manner, a back projection of the projection data set 50 having the first, low resolution may initially be performed separately from a back projection of the second projection data set 52 having the second, higher resolution. Provided that in the first projection data set 50, the object comprising the object region of interest is fully imaged, a CT reconstruction of the first projection data set may be performed using conventional CT image reconstruction algorithms. They include, for example, applying globally effective filters to the individual projection recordings so as to obtain filtered projection recordings 54. By means of back-projecting the filtered projection recordings, a low-resolution CT reconstruction 56 of the data may occur. Irrespective thereof, the projection recordings associated with the second, higher-resolution projection data set 52 may be filtered using a location- and frequency-local filter (corresponding to wavelet decomposition, for example) on the basis of said second projection data set 52, so as to obtain a high-frequency-filtered representation of the projection recordings of the higher resolution 58.
By means of the back projection 60, an intermediate CT reconstruction of the object is then produced within the object region of interest, which intermediate CT reconstruction contains only information about high-frequency image portions (intermediate back projection or intermediate reconstruction 62). If the intermediate CT reconstruction contains only image portions of a frequency that is above a cutoff frequency defined by the finite resolution of the first projection recordings, the intermediate CT reconstructions 56 and 62, which have been produced separately, may be combined by means of a combiner 64 so as to obtain, after the combination step has been effected, a CT reconstruction 66 of the object wherein the object region of interest comprising the second, higher resolution is contained without image artefacts having been produced by the image reconstruction. During combining 56 and 62, 56 is possibly treated with an operator associated with the filters, applied to the projection, for reconstructing the low-frequency portions, and 62 is treated with an operator associated with the filters, applied to the projection, for reconstructing the high-frequency portions, before the two intermediate CT reconstructions which have been pre-processed in this manner are added.
In other words,
With the CT reconstruction, any CT reconstruction algorithms may be used as have been known so far or will be developed in the future.
In a pre-transformation step 70, the second projection data set, or the projection recording of the second projection data set, is pre-transformed so as to obtain a pre-transformed second projection data set 58 which merely comprises image frequencies above a predetermined frequency. This frequency may be, e.g., the cutoff frequency of the first projection data set, or any other frequency matched or tuned to the cutoff frequency. The pre-transformation may be high-pass filtering, for example.
However, any other pre-transformations which result in that only image frequencies above a predetermined cutoff frequency are contained within the pre-transformed second projection data set are also applicable. Contained merely within the second projection data set evidently also means that the lower frequencies need not be fully removed, but that it is sufficient to suppress them with regard to the high-frequency image portions.
Parallel therewith, the first projection data set may optionally be subject to low-pass filtering adapted to the cutoff frequency of the high-pass filtering of the second projection data set. Alternatively, said low-pass filtering may also be dispensed with, however, since, due to the finite detector resolution with which the first projection data set was obtained, the data of the first projection data set per se have a cutoff frequency above which no image frequencies can be contained within the projection recordings.
In accordance with some embodiments of the present invention, while combining the first and the pre-transformed second projection data sets in accordance with a combination rule, a CT reconstruction 66 is eventually produced which in the first region of the object comprises the first resolution, and in the object region of interest comprises the second, higher resolution. In accordance with the embodiment described in
The CT reconstruction 66 is then obtained in that a CT reconstruction algorithm is applied to the third projection data set.
In other words, in
The method illustrated in
The difference consists in that on the basis of the first projection data set and the pre-transformed second projection data set, separate intermediate CT reconstructions are initially produced which are subsequently copied to become the final CT reconstruction 66. In other words, a first intermediate CT reconstruction 90 is produced from the first projection data set 50 or its pre-transformed representation 82. By analogy therewith, a second intermediate CT reconstruction 92 is produced from the pre-transformed second projection data set 58, said second intermediate CT reconstruction 92 obviously comprising, even in the 3-dimensional representation, only image frequencies above the predetermined cutoff frequency. The CT reconstruction is then obtained by combining the first intermediate CT reconstruction 90 and the second intermediate CT reconstruction 92. In accordance with an embodiment of the invention, these two intermediate CT reconstructions 90 and 92 are combined by means of a wavelet synthesis.
I.e., the method suggested in
The method described in
In other words, it is therefore possible to enable—from two or more recordings of the same object with different geometric magnifications—reconstruction of large measuring regions with such a high spatial resolution as cannot be achieved with a single measurement (multi-scan CT, multi-scan ROI CT).
Although the inventive embodiments have been discussed mainly in the context of material tests, it is readily possible to apply them to diagnosing methods in humans. When used for non-destructive testing of materials, the method may be used, in particular, for testing large components that have been produced by means of any production methods and from any materials for which, in their inner sub-regions, particularly high demands are placed on the spatial resolution that is effective in the testing. As an example, one might mention the testing of light metal casting pistons or turbine screws for power plants or generators. The testing precision is greatly improved by the above-described concept, both with regard to the sensitivity to cracks, to voids and pores, and with regard to the specificity in separating correct and erroneous positive defect findings.
In addition, applying the inventive concept enables simultaneously, i.e. with only one single CT construction, testing the further regions of the object with a low spatial resolution, it being possible to supplement the testing with the low spatial resolution by further testing by means of a method with a lower penetration depth.
Depending on the circumstances, the inventive method for CT reconstruction of an object may be implemented in hardware or in software. The implementation may be on a digital storage medium, in particular a disc or a CD having electronically readable control signals which can cooperate with a programmable computer system such that the inventive method for CT reconstruction of an object is performed. Generally, the inventive thus also consists in a computer program product having a program code, stored on a machine-readable carrier, for performing the inventive method, when the computer program product runs on a computer. In other words, the invention may thus be realized as a computer program having a program code for performing the method, when the computer program runs on a computer.
While this invention has been described in terms of several embodiments, there are alterations, permutations, and equivalents which fall within the scope of this invention. It should also be noted that there are many alternative ways of implementing the methods and compositions of the present invention. It is therefore intended that the following appended claims be interpreted as including all such alterations, permutations and equivalents as fall within the true spirit and scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
102008021639.9 | Apr 2008 | DE | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/EP09/03149 | 4/30/2009 | WO | 00 | 6/3/2010 |