Scaled radiography reconstruction

Abstract

The invention relates to off-center detector 3D X-ray or proton radiography reconstruction. Redundancy weighting with a steep weighting function around the iso-axis typically leads to artifacts in the reconstruction, for example, if inconsistencies between two nominal redundant projections occur, e.g. due to slightly incorrect detector calibration or scatter correction, etc. With the present invention, an approach is presented for overcoming or mitigating these problems.

Claims

1. An apparatus for off-center detector 3D X-ray reconstruction or for proton radiography reconstruction, comprising a processing unit arranged for obtaining acquisition data measured by an off-center X-ray detector or by a proton detector, wherein the acquisition data includes a first portion and a second portion, wherein the first and second portion both cover at least a detector overlap region, determining a first reference value indicative of measured line integral values for the first portion in the detector overlap region and determining a second reference value of measured line integral values for the second portion in the detector overlap region, scaling at least one of the first portion and the second portion of the acquisition data based on a ratio of the first and second reference values, and performing reconstruction using the result of the scaling.

2. The apparatus for proton radiography reconstruction according to claim 1, further comprising the proton detector.

3. The apparatus according to claim 2, wherein the proton detector is arranged for performing multiple measurements with different offsets.

4. The apparatus for off-center detector 3D X-ray reconstruction according to claim 1, further comprising the off-center 3D X-ray detector.

5. A method for off-center detector 3D X-ray reconstruction or for proton radiography reconstruction, comprising the steps of obtaining acquisition data measured by an off-center X-ray detector or by a proton detector, wherein the acquisition data includes a first portion and a second portion, wherein the first and second portion both cover at least a detector overlap region, determining a first reference value indicative of measured line integral values for the first portion in the detector overlap region and determining a second reference value of measured line integral values for the second portion in the detector overlap region, scaling at least one of the first portion and the second portion of the acquisition data based on the first and second reference values, and performing reconstruction using the result of the scaling.

6. The method according to claim 5, wherein only the second portion of the acquisition data is scaled with a ratio of the first reference value and the second reference value.

7. The method according to claim 5, wherein the first portion is scaled with a ratio of the first reference value and a base value and the second portion is scaled with a ratio of the base value and the second reference value.

8. The method according to claim 5, wherein the first portion is scaled with a ratio of the first reference value and a first base value obtained by forward projection of a planning CT of a first arc and the second portion is scaled with a ratio of a second base value corresponding to the first base value recomputed for angles of a second arc and the second reference value.

9. The method for according to claim 5, wherein the acquisition data is obtained in a course of a rotation of 360° with a fixed detector off-set or by two scans of more than 180° and shifting the detector from one side to another side.

10. The method according to claim 5, wherein the performing of the reconstruction includes a redundancy weighting.

11. The method for according to claim 5, wherein the acquisition data is obtained by realizing an arbitrary trajectory around an examination object.

12. A software product stored on a non-transitory computer-readable medium for off-center detector 3D X-ray reconstruction or proton radiography reconstruction, the software product comprising program code means for causing a computer to carry out the steps of the method as claimed in claim 5 when the software product is run on the computer.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) In the following drawings:

(2) FIG. 1 shows an illustration of overlap and redundant data regions for short scan and 360° acquisitions,

(3) FIG. 2 shows a reconstruction result of conventional processing and of a processing according to the invention in comparison to a reference for a simulation with a drift of the magnitude of the line integral values of +/−5%,

(4) FIG. 3 shows schematically a system including an apparatus according to an embodiment of the invention, and

(5) FIG. 4 shows a schematic flow diagram illustrating a method according to an embodiment of the invention.

DETAILED DESCRIPTION OF EMBODIMENTS

(6) In the context of the invention, it is proposed to use the (quasi-) redundant rays measured with different detector offset to force equality in a certain projection region to mitigate the artifacts produced by inconsistencies, as discussed above and explained with respect to embodiments as indicated below.

(7) In a certain area in the region where (quasi-) redundant data are measured, the average (or an indicative value thereof) of the measured line integral values for each projection is logged.

(8) In an embodiment, the pre-processing and back-projection of the first half of the complete acquisition is unchanged, while for the second half, the measured line integral on the region of interest is compared to that with the projection containing the (quasi-) redundant data. The measured line integral values of the second half acquisition are scaled with the ratio of the two line integral values.

(9) The selection of areas of redundant data is illustrated in FIG. 1 and reconstruction results are shown in FIG. 2.

(10) FIG. 1 shows an illustration of overlap and redundant data regions for short scan and 360° acquisitions.

(11) The first portion 1 and the second portion 2 of the acquisition data overlap in a detector overlap region 3, while this depends on the size of the physical detector and the offset.

(12) In case of two short scans, the detector overlap region 3 is used for calculating the line integral values.

(13) In case of a 360° acquisition, there is a small(er) region 4 of quasi-redundant data, which may be used for calculating the line integral values.

(14) FIG. 2 shows a reconstruction result of conventional processing and of a processing according to the invention in comparison to a reference for a simulation with a drift of the magnitude of the line integral values of +/−5%.

(15) In the upper portion of FIG. 2, axial views of the reference (on the left), a conventional reconstruction (in the middle) and a reconstruction according to the invention (on the right) are shown, while the lower portions shows corresponding coronal views. The reconstructions results are for a +/−5% drift of magnitude in the line integral values for a 360° acquisition with a 40 mm detector overlap (L/W 30/150 HU, single slice).

(16) FIG. 3 shows schematically a system including an apparatus 10 according to an embodiment of the invention.

(17) The apparatus 10 for reconstruction is coupled to an off-center detector X-ray apparatus 20 and a memory 30, while the X-ray apparatus 20 and the memory 30 are also coupled.

(18) The X-ray apparatus 20 is arranged to obtain the raw data of the measurement by carrying out the X-ray operation and may transfer the raw data (directly or after some preliminary processing) either to the reconstruction apparatus 10 or to the memory 30.

(19) The apparatus 10 for reconstruction thus may obtain the data on which the reconstruction is based from the X-ray apparatus 20 or from the memory 30.

(20) The apparatus 10 is arranged to perform the reconstruction according to the present invention and is furthermore arranged to forward the reconstruction data to a display 40 for display and/or to the memory 30 for storage (and later display).

(21) FIG. 4 shows a schematic flow diagram illustrating a method according to an embodiment of the invention

(22) In an obtainment step 100 acquisition data is obtained, e.g. directly from an X-ray apparatus including an off-center X-ray detector (a concept with which the skilled person is well familiar, so that no further explanation is needed) or from a memory or the like in which such data is stored.

(23) The acquisition data includes a first portion and a second portion, wherein the first and second portion both cover at least a detector overlap region, as it is illustrated in FIG. 1.

(24) In a following determination step 101, a first reference value indicative of measured line integral values for the first portion in the detector overlap region is determined, while, in parallel or separately, also a second reference value of measured line integral values for the second portion in the detector overlap region is determined.

(25) Based on the first and second reference values, in a scaling step 102 at least one of the first portion and the second portion of the acquisition are scaled.

(26) Using the result of the scaling, in a reconstruction step 103, the reconstruction is performed, while this reconstruction process as such corresponds to the conventional approaches.

(27) While the invention has been illustrated and described in detail in the drawings and foregoing description, such illustration and description are to be considered illustrative or exemplary and not restrictive; the invention is not limited to the disclosed embodiments.

(28) The invention can, in particular, be used for or implemented in the context of C-arm systems, RT on-board imager, Spect/CT systems, etc.

(29) Furthermore, while it is realized that the above discussion of exemplary embodiments refers to the case of an off-center 3D X-ray reconstruction, it will be appreciated by the skilled person that the above discussion similarly applies to, for example, PT radiography (image stitching with different detector offsets at same angle) or Proton CT.

(30) Other variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed invention, from a study of the drawings, the disclosure, and the appended claims.

(31) In the claims, the word “comprising” does not exclude other elements or steps, and the indefinite article “a” or “an” does not exclude a plurality.

(32) A single processor, device or other unit may fulfill the functions of several items recited in the claims. The mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measures cannot be used to advantage.

(33) Operations like obtaining or acquiring data, determining or calculating values, scaling and performing reconstruction can be implemented as program code means of a computer program and/or as dedicated hardware.

(34) A computer program may be stored and/or distributed on a suitable medium, such as an optical storage medium or a solid-state medium, supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems.

(35) Any reference signs in the claims should not be construed as limiting the scope.