Dual energy CT (DECT) acquires an object at two different detected spectra w1(E) and w2(E). Rawdata-based techniques can be applied whenever the integration lines (rays) of both spectra are identical. Then, the rawdata q1 and q2 undergo a decomposition process of type p = p(q1, q2). For spiral cone-beam dual source CT (DSCT) and for micro-DSCT the integration lines are disjunct and therefore inconsistent. Hence, one typically uses image-based subtraction techniques that are inferior to the rawdata-based methods. Our technique provides consistent rawdata in two steps. Starting from standard reconstructions of both scans, polychromatic forward projections of the volumes 1 and 2 yield the desired consistent rawdata q1 and q2. The second step computes the decomposition p(q1, q2) and reconstructs these sinograms using filtered backprojection. All computations are done in parallel beam geometry. To avoid loosing spatial resolution the linear terms of p are handled in image domain. Only the higher order terms undergo the two-step process. The algorithm was run on a Cell Broadband Engine (CBE) at 3.2 GHz in a Playstation 3 (PS3) computer (Sony Computer Entertainment, Japan). To evaluate our approach we decompose patient rawdata acquired with a SOMATOM Definition clinical DSCT scanner (Siemens Medical Solutions, Forchheim, Germany). The PS3 processes 35 slices per second (5122 pixels, 512 views per half rotation, two forward projections, decomposition, filtering, one backprojection); a typical volume of thousand slices is processed in half a minute.
[1]
A. Coleman,et al.
A beam-hardening correction using dual-energy computed tomography.
,
1985,
Physics in medicine and biology.
[2]
A. Macovski,et al.
Energy-selective reconstructions in X-ray computerised tomography
,
1976,
Physics in medicine and biology.
[3]
Michael Knaup,et al.
Hyperfast parallel-beam and cone-beam backprojection using the cell general purpose hardware.
,
2007,
Medical physics.
[4]
W. Kalender,et al.
An algorithm for noise suppression in dual energy CT material density images.
,
1988,
IEEE transactions on medical imaging.
[5]
W. Kalender,et al.
Evaluation of a prototype dual-energy computed tomographic apparatus. I. Phantom studies.
,
1986,
Medical physics.
[6]
Lellery Storm,et al.
Photon cross sections from 1 keV to 100 MeV for elements Z=1 to Z=100
,
1970
.
[7]
G. Barnes,et al.
Semiempirical model for generating tungsten target x-ray spectra.
,
1991,
Medical physics.
[8]
W. Kalender,et al.
Empirical Dual Energy Calibration (EDEC) for Cone-Beam Computed Tomography
,
2006,
2006 IEEE Nuclear Science Symposium Conference Record.