Dual energy CT for attenuation correction with PET/CT Ting Xia, Adam M. Alessio, and Paul E. Kinahan Citation: Medical Physics 41, 012501 (2014); doi: 10.1118/1.4828838 View online: http://dx.doi.org/10.1118/1.4828838 View Table of Contents: http://scitation.aip.org/content/aapm/journal/medphys/41/1?ver=pdfcov Published by the American Association of Physicists in Medicine Articles you may be interested in Evaluation of respiratory and cardiac motion correction schemes in dual gated PET/CT cardiac imaging Med. Phys. 41, 072504 (2014); 10.1118/1.4881099 CT head-scan dosimetry in an anthropomorphic phantom and associated measurement of ACR accreditationphantom imaging metrics under clinically representative scan conditions Med. Phys. 40, 081917 (2013); 10.1118/1.4815964 Voxel clustering for quantifying PET-based treatment response assessment Med. Phys. 40, 012401 (2013); 10.1118/1.4764900 Characterization of adaptive statistical iterative reconstruction algorithm for dose reduction in CT: A pediatric oncology perspective Med. Phys. 39, 5520 (2012); 10.1118/1.4745563 Single and dual energy attenuation correction in PET/CT in the presence of iodine based contrast agents Med. Phys. 35, 1959 (2008); 10.1118/1.2903476

Dual energy CT for attenuation correction with PET/CT Ting Xia Department of Bioengineering, University of Washington, Seattle, Washington 98105

Adam M. Alessio Department of Radiology, University of Washington, Seattle, Washington 98105

Paul E. Kinahana) Departments of Radiology and Bioengineering, University of Washington, Seattle, Washington 98105

(Received 8 May 2012; revised 15 October 2013; accepted for publication 15 October 2013; published 4 December 2013) Purpose: The authors evaluate the energy dependent noise and bias properties of monoenergetic images synthesized from dual-energy CT (DECT) acquisitions. These monoenergetic images can be used to estimate attenuation coefficients at energies suitable for positron emission tomography (PET) and single-photon emission computed tomography (SPECT) imaging. This is becoming more relevant with the increased use of quantitative imaging by PET/CT and SPECT/CT scanners. There are, however, potential variations in the noise and bias of synthesized monoenergetic images as a function of energy. Methods: The authors used analytic approximations and simulations to estimate the noise and bias of synthesized monoenergetic images of water-filled cylinders with different shapes and the NURBSbased cardiac-torso (NCAT) phantom from 40 to 520 keV, the range of SPECT and PET energies. The dual-kVp spectra were based on the GE Lightspeed VCT scanner at 80 and 140 kVp with added filtration of 0.5 mm Cu. The authors evaluated strategies of noise suppression with sinogram smoothing and dose minimization with reduction of tube currents at the two kVp settings. The authors compared the impact of DECT-based attenuation correction with single-kVp CT-based attenuation correction on PET quantitation for the NCAT phantom for soft tissue and high-Z materials of bone and iodine contrast enhancement. Results: Both analytic calculations and simulations displayed the expected minimum noise value for a synthesized monoenergetic image at an energy between the mean energies of the two spectra. In addition the authors found that the normalized coefficient of variation in the synthesized attenuation map increased with energy but reached a plateau near 160 keV, and then remained constant with increasing energy up to 511 keV and beyond. The bias was minimal, as the linear attenuation coefficients of the synthesized monoenergetic images were within 2.4% of the known true values across the entire energy range. Compared with no sinogram smoothing, sinogram smoothing can dramatically reduce noise in the DECT-derived attenuation map. Through appropriate selection of tube currents for high and low kVp scans, DECT can deliver roughly the same amount of radiation dose as that of a single kVp CT scan, but could be used for PET attenuation correction with reduced bias in contrast agent regions by a factor of ∼2.6 and slightly reduced RMSE for the total image. Conclusions: When DECT is used for attenuation correction at higher energies, there is a noise amplification that is dependent on the energy of the synthesized monoenergetic image of linear attenuation coefficients. Sinogram smoothing reduces the noise amplification in DECT-derived attenuation maps without increasing bias. With an appropriate selection of CT techniques, a DECT scan with the same radiation dose as a single CT scan can result in a PET image with improved quantitative accuracy. © 2014 American Association of Physicists in Medicine. [http://dx.doi.org/10.1118/1.4828838] Key words: dual energy CT, PET/CT, attenuation correction, noise and bias 1. INTRODUCTION Quantitative imaging by positron emission tomography (PET) and single-photon emission computed tomography (SPECT) is receiving increased attention for clinical applications.1–5 In quantitative PET and SPECT, correction for the effect of photon attenuation is of paramount importance.6 The xray computed tomography (CT) component in a PET/CT or SPECT/CT system not only provides precise anatomical localization of regions identified on the tracer uptake images, 012501-1

Med. Phys. 41 (1), January 2014

but is also used for attenuation correction of the PET or SPECT emission data.7–10 At the energies of x-ray CT, attenuation is due to Compton scatter and photoelectric absorption, while at SPECT energies, and particularly PET energies, Compton scatter is the dominant process for biological materials. A list of energies of interest are given in Table I, and Fig. 1 plots the mass attenuation coefficients of common materials over the energy ranges relevant to PET/CT and SPECT/CT imaging.

0094-2405/2014/41(1)/012501/11/$30.00

© 2014 Am. Assoc. Phys. Med.

012501-1

012501-2

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

TABLE I. Photon energies of common isotopes used in PET/CT and SPECT/CT Imaging. Mode PET SPECT SPECT SPECT SPECT SPECT

Isotope

Energy (keV)

All (i.e. F-18, C-11, N-13, etc.) Xe-133 Tc-99m I-123 Ga-67 In-111

511 80.9 140 160 185, 93, 300 171, 245

CT-based attenuation correction offers several benefits over transmission scanning with an isotope at or near the energy of the nuclear tracer isotope.9 CT scans are substantially faster than lengthy nuclear transmission scans and can provide attenuation correction factors that are orders of magnitude lower in noise. Moreover, CT scans contain no bias from the emission contamination present in postinjection imaging.9 CT attenuation correction factors are, however, susceptible to bias due to the need to transform the polyenergetic CT image to attenuation correction factors at the nuclear isotope energy, typically either at 140 or 511 keV. Three potential solutions are employed for transforming attenuation maps from CT energies to SPECT and PET energies:9 segmentation, linear scaling, and multilinear scaling. Segmentation methods can be used to separate the CT image into regions corresponding to different tissue types, which then are replaced with appropriate attenuation coefficients at the appropriate energy. This method has a potential source of error due to tissue misclassification and the assumption of uniform attenuation throughout a region.11 A second method uses linear scaling of the entire CT image with the ratio of attenuation coefficients of water (representing soft tissues) at the photon energies of CT and the target energy, offering a simple solution for the transformation.7 For bone, however, linear scaling is a poor approximation, since the photoelectric absorption to

012501-2

Compton scatter ratio changes differently for bone compared to soft tissue across the range of energies of interest.7 Finally, the most common method for transforming to SPECT or PET energies is Multilinear/Hybrid Scaling.8–12 In these methods, different scaling factors (for water and air and for water and bone, respectively) are used to calculate the attenuation values for CT numbers H for which −1000 < H < 0, and for H > 0. In practice, the multilinear scaling method and other hybrid methods have been shown to give reasonable results for low-Z biological materials. However, high-Z materials, such as contrast agents, can introduce significant bias.9, 13–16 Dual energy (It is more accurate to use the term dualkVp CT, since each CT acquisition is polyenergetic, not monoenergetic. However, the term dual-energy CT is already in wide use.) CT (DECT) is an accepted method to estimate material properties using either photoelectric and Compton components17 or physical basis materials such as plastic and aluminum.18 DECT has been proposed to remove the bias from the CT attenuation correction for SPECT (Refs. 19 and 20) and PET.9, 20 A challenge with the use of DECT is the significant noise amplification due to the poorly conditioned inverse problem of estimating the component sinograms,21 leading to excessive noise amplification. It has already been shown in the literature that the noise amplification is energy dependent in the diagnostic CT energy range,17, 22, 23 but the noise and bias at SPECT or PET energies are largely unknown. Another potential drawback of dual energy is the additional patient radiation dose required to acquire two separate energies and/or to reduce noise. To address these issues with DECT-based attenuation correction we have previously investigated alternate methods using iterative CT reconstruction24 and image-based strategies.25 In this study, we aim to (1) evaluate the energy dependent noise properties of monoenergetic images synthesized from dual-energy CT in the energy range appropriate for nuclear imaging (140-511 keV) and (2) minimize noise and radiation dose in DECT-based attenuation correction.

2. MATERIALS AND METHODS

F IG . 1. Mass attenuation coefficient as a function of energy for different materials. Medical Physics, Vol. 41, No. 1, January 2014

We used three approaches to study the noise and bias properties of using DECT for PET attenuation correction as a function of energy: analytical approximations to directly calculate sinogram noise, simulations of cylindrical water phantoms to measure sinogram noise, and simulations of cylindrical water phantoms and the NCAT phantom to measure noise and bias in the reconstructed PET images. In Sec. 2.A, we describe the theoretical prediction that uses an approximate analytical method. In Sec. 2.B, we outline the simulation tools. In Sec. 2.C, we measure the energy dependent noise of the DECT derived monoenergetic sinograms and corresponding reconstructed images for two different size water phantoms. A simple noise suppression and dose minimization strategy is proposed in Sec. 2.D. In Sec. 2.E, we evaluate the impact of DECT-based attenuation correction on PET image quantitation with the NCAT phantom.

012501-3

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

2.A. Theoretical prediction

We first used analytic approximations to estimate the variance of the synthesized monoenergetic sinogram. This

012501-3

approach was based on the work of Doost-Hoseini and others,18, 22, 23 where the variance of an entry in the monoenergetic sinogram is approximated by

2 σ12 2 σ2 + (f (E )k − f (E )k ) a m 1b b m 1a I12 I22 , 2 (k1a k2b − k1b k2a )

(fa (Em )k2b − fb (Em )k2a )2 σm2 =

where I1 and I2 are entries in the measured signals for the two kVp scans, σ 2 1 and σ 2 2 are the variances for the two kVp scans, fa (E) and fb (E) are attenuation coefficients of the basis materials at energy E, Em is the target energy, kia = (1 / Ii )(∂Ii / ∂Aa ) is the variation of measured signal intensity with basis material, and Aa is the line integral through the basis material a. It was shown that σ m reached its minimum at an optimum target energy Eopt , which was a function of object thickness.23 We used Eq. (1) to estimate the energydependent coefficient of variation (COV) for an object at the monoenergetic target energy Em as σm , (2) COVEm = μm where σ m is the noise of the object and μm is the linear attenuation coefficient of the object at energy Em . We calculated the normalized COVEm for the target energy Em ranging from 10 to 525 keV with 1 keV increments, where the normalized COVEm was the COVEm divided by the minimum value of COV at the optimal energy. We note that the derivation of Eq. (1) used a first order Taylor approximation and assumed monoenergetic x-rays,23 thus not accounting for beam-hardening effects. In order to confirm this analytic relationship, we used polyenergetic simulations to estimate the noise of the synthesized monoenergetic data, σ m , over the same energy range of 10–525 keV. 2.B. Simulation tools

2.B.1. CT and PET simulators

We used the Computer Assisted Tomography SIMulator (CATSIM) for simulation of the x-ray CT imaging.26 It includes a Monte Carlo-based simulation package called CATDOSE for estimation of the CT radiation dose. Both CATSIM and CATDOSE have been validated relative to GE Lightspeed Volumetric Computed Tomography (Lightspeed VCT, General Electric Medical Systems, Waukesha, WI) measurements for a series of test objects.26 The CT simulation models the effects of the x-ray source focal spot, polychromatic tube spectra, beam conditioning, bowtie filter, and beam hardening. Although fully described in Ref. 26, we provide a brief summary here of the properties modeled by CATSIM. CT raw data yi at sinogram index, i, are formed from the combination of a quantum noise process pi and electronic noise process zi , i.e., yi = pi + zi , following methods Medical Physics, Vol. 41, No. 1, January 2014

(1)

described by De Man et al.,27 similar to the methods described by others.28, 29 The quantum noise is modeled as    Ek · xik fC , (3) pi = k

where E refers to photon energy, k is the energy bin index, fC is a factor to convert from keV to the number of electrons, and xik is a Poisson random process with mean λik , xik ∼ P(λik ), where   1   liso μok + yik (4) λik = η · Nik · exp − S s o and η is the detector quantum efficiency (fraction of photons absorbed), Nik is the number of photons arriving at the detector without attenuation for energy bin k, s is the “tube-ofresponse” subsampling index for the total of S subsamples (where the tube-of-response between the source and detector is subsampled to better model the finite beam width), liso is the intersection length between the line with index is and the object element with index o, μok is the linear attenuation co is the estimated efficient of object at o at energy k, and yik scatter at sinogram index i with energy index k. The electronic noise zi is modeled as a Gaussian random process, zi ∼ N(d, σe 2 ), where d is the dark current, and σ e is the standard deviation of the electronic noise. The PET images were simulated using a simplified version of ASIM,30 which included the physical effects of attenuation, background noise, and photon counting. The noisy prompt coincidences p˜ follow a Poisson pseudorandom process with mean m, (p˜ ∼ P (m)), defined as m = αt · t + αb · b,

(5)

where t represents the true coincidences, b is the background noise including scattered and random coincidences, and α t and α b are global scale factors for true coincidences and background noise. 2.B.2. Test objects

Three test phantoms were simulated. The first was a 10 cm diameter and 10 cm long water cylinder. The second was a uniform 20 × 30 cm elliptical water cylinder, with a 10 cm axial dimension. The third was derived from the NURBS-based cardiac-torso (NCAT) phantom.31, 32 The NCAT phantom had three regions with focal FDG uptake: A 1.6-cm-diameter bone lesion, a 3-cm-diameter contrast enhanced liver lesion, and a

012501-4

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

3-cm-diameter soft tissue lesion. In terms of mass percentage, the iodine contrast is 1% in blood, and equivalent to 10.6 mg/ml in our simulation. The PET activity distribution was set to tracer uptake ratios of (background: liver: lung: lesion) = (1:3:0.5:6). The parameters α t and α b in Eq. (5) were chosen to match total typical clinical counts (10 × 106 ) in the sinogram with a scatter fraction of 30%. The two water phantoms allowed evaluation of energy dependent noise and bias in the DECT derived attenuation map for uniform soft tissues. The rotational symmetry of the first 10 cm diameter water cylinder allowed evaluation of variance in the central projection bin of the synthesized monoenergetic sinograms, which assumed that noise is uncorrelated across view angles. The second elliptical water phantom invoked slightly more realistic beam hardening effects and allowed evaluation of energy dependent noise and bias in the DECT derived attenuation maps. The elliptical phantom was also used for dose estimation and optimization in a torsolike object. The NCAT phantom was used for a more realistic study of the quantitative accuracy of the estimated PET tracer uptake with selected combinations of DECT acquisition parameters for comparison to single-kVp CT-based attenuation correction with similar radiation dose levels. 2.B.3. Dual-energy basis material decomposition

We used the basis material decomposition (BMD) technique,18, 33 with basis materials of aluminum and plastic (polyethylene). For calibration we followed the procedure of Chuang et al.33 and defined polyethylene (Pl, ρ = 0.93 g/cm3 ) and aluminum (Al, ρ = 2.72 g/cm3 ) step wedge phantoms.34 The calibration and measurement procedures were as follows: Calibration: 1. Define step wedges made of aluminum and plastic for calibration (Al: 0–8 steps, 4 mm/step; Plastic: 0–17 steps, 12.5 mm/step). 2. Generate sinograms from high and low kVp scans for the step wedges. 3. Calculate In(I0 /I) for each combination of aluminum and plastic steps. 4. Generate calibration coefficients {di }, {ej } for each material: Aα = d0 + d1 TH + d2 TL + d3 TH TL + d4 TH2 + d5 TL2 + d6 (TH TL )2 + d7 TH3 + d8 TL3 ,

(6)

Aβ = e0 + e1 TH + e2 TL + e3 TH TL + e4 TH2 + e5 TL2 + e6 (TH TL )2 e7 TH3 + e8 TL3 ,

(7)

where Aα and Aβ correspond to the thicknesses of the two materials, with units of cm, and TH = In(I0 /I)H and TL = In(I0 /I)L are for the high (H) and low (L) kVp scans. Measurement: 5. Generate high and low kVp sinogram data for the test object at the given settings (mAs). Medical Physics, Vol. 41, No. 1, January 2014

012501-4

6. Using the generated coefficients {di }, {ej }, obtain basis material sinograms of Aα and Aβ for the test object according to Eqs. (6) and (7). 7. Obtain the basis material attenuation sinograms at desired energy E, using known attenuation values versus E for each material:  (8) μα (E) · ds = Aα · μα (E),  μβ (E) · ds = Aβ · μβ (E).

(9)

8. Sum the attenuation sinograms of each material to obtain the total attenuation sinogram of the test object at the desired energy E.    (10) μtest (E) · ds = μα (E) · ds + μβ (E) · ds. 9. Reconstruct attenuation maps at desired energy E.

2.C. Evaluation of energy-dependent noise and bias

In this subsection, we evaluated energy-dependent noise and bias on the dual energy CT derived monoenergetic sinograms and corresponding reconstructed attenuation maps using simple cylindrical and elliptical water phantoms. We simulated CT step-and-shoot mode with a gantry rotation time of 0.5 s. The dual-energy spectra were based on the GE Lightspeed VCT scanner at 80 and 140 kVp, with added 0.5 mm copper filtration following the bowtie filter. Multiple x-ray tube currents were simulated at [64, 91, 110, 128, 181, 256, 362, 512, 724, and 1024] mA for each kVp scan. Synthesized monoenergetic data and images were formed at [50, 60, 70, 80, 90, 100, 140, 171, 185, 245, and 511] keV covering the range of photon energies for SPECT and PET imaging listed in Table I. The CT sinogram was 888 channels × 984 views. Default ideal kernel filter and water-only beam hardening correction in CATSIM package were applied. The CT slice thickness was 3.125 mm, comparable with the PET slice thickness (for studies in Sec. 2.E). All DECT-derived attenuation maps were reconstructed to the size of 128 × 128 pixels using filtered back projection (FBP) over a 50 cm field of view (FOV) to match PET image dimensions. We estimated the COV according to Eq. (2) in the sinogram for the 10 cm diameter cylinder by estimating the mean attenuation and the standard deviation for the central sinogram bin across view angles. We also estimated the COV in the image domain using Eq. (2) as follows: For the 10 cm diameter water cylinder, a centered 2D region of interest (ROI) with diameter of 6 cm was used for analysis. For the 20 × 30 cm elliptical water cylinder, a total of 18 nonoverlapping 2D ROIs, each with radius of 2 cm, were drawn on the cross-section of the phantom as shown in Fig. 2. The averaged mean and COV for all ROIs were calculated. The bias was calculated as the difference between the average of the mean pixel value over all ROIs and the true attenuation value of water.

012501-5

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

012501-5

2.E. Comparison of PET attenuation correction with dual- and single kVp-derived attenuation map

To study the impact of DECT-based attenuation correction on PET quantitation, we used the NCAT phantom with iodinebased contrast agent to compare both DECT and single kVp CT attenuation corrected PET images. 2.E.1. Single kVp CT-based PET attenuation correction

F IG . 2. Elliptical cylindrical phantom with 30 × 20 cm cross-section, 18 2D ROIs each with radius of 2 cm were drawn for evaluation.

2.D. DECT noise suppression and dose minimization for attenuation correction

We then evaluated the effects of noise suppression and dose minimization on the DECT derived attenuation data, again using the elliptical water cylinder. To further reduce DECT noise, the high and low kVp CT sinograms of the elliptical water cylinder acquired with the tube current of 64, 91, 110, 128, 181, 256, 362, 512, 724, and 1024 mA for each kVp scan were filtered by 2D 5 × 5 boxcar smoothing, followed by an adaptive trimmed mean filter (ATM) 35–37 before logarithmic conversion. Then, with basis material decomposition method and image reconstruction similar to the procedure performed in Sec. 2.C, we generated the sinogram smoothed DECT-derived attenuation map for each synthesized energy of 50, 60, 70, 80, 90, 100, 140, 171, 185, 245, and 511 keV. For each synthesized energy, we evaluated the COV in the reconstructed attenuation map as a function of tube currents for high and low kVp scans. With the simulation package CATDOSE, we estimated the total radiation dose for each combination of tube currents at high and low kVp scans for the 20 × 30 cm elliptical water cylinder. Comparing the dose distribution and COV distribution as a function of tube currents for the two kVp CT scans, we were able to estimate the parameters for the DECT acquisition for attenuation correction with the minimum radiation dose. This was equivalent to estimating the parameters for the DECT acquisition for attenuation correction with the minimum COV for the given radiation dose. For comparison, a single kVp CT scan of the 20 × 30 cm elliptical water cylinder was obtained. The CT technique was selected to be comparable to clinical low dose CTAC protocol38 with a tube voltage of 120 kVp, tube current of 100 mA, and gantry rotation time of 0.5 s. The spectra were conditioned with a bowtie only. The slice thickness was the same as the DECT images. Radiation dose was estimated through CATDOSE. Medical Physics, Vol. 41, No. 1, January 2014

Acquisitions of NCAT phantom were simulated with CATSIM, which uses analytic projectors for a step-and-shoot acquisition with CT tube voltage of 120 kVp (with a bowtie filter but no other filtration), tube current of 100 mA, and gantry rotation time of 0.5 s. The CT slice thickness was 3.125 mm, comparable with the PET slice thickness. The CT sinogram was 888 channels × 984 views. Default ideal kernel filter and water-only beam hardening correction in CATSIM package were applied. The CT image was reconstructed to 128 × 128 image matrix, using filtered back projection (FBP) over a 50 cm field of view (FOV) to match PET image dimensions. Then a modified version of the bilinear transformation method that relates the CT number to the attenuation coefficient at 511 keV was used to generate the attenuation map for PET.9, 39 In the modified method, the slopes for the transformation line of the soft-tissue and bone were estimated with the best fit using the known true values. With the generated attenuation map, both noise-free and noisy attenuation corrected PET images were reconstructed using weighted OSEM (16 subsets, 14 iterations) on 128 × 128 pixel matrix over a 50 cm FOV. 2.E.2. DECT-based PET attenuation correction

The data flow of DECT-based PET attenuation correction process is shown in Fig. 3. The CT tube currents for the 140 and 80 kVp scans were selected based on the optimal estimation result of Sec. 2.D for the elliptical water cylinder, with added 0.5 mm copper filtration following the bowtie filter. This combination yielded the same radiation dose as that of 120 kVp 50 mAs single-kVp scan. The CT sinograms for the high and low kVp scans were then smoothed with a 2D 5 × 5 boxcar followed by the adaptive trimmed mean filter (ATM) before logarithmic conversion for noise suppression. With basis material decomposition method, the synthesized attenuation map at 511 keV was reconstructed with FBP to 128 × 128 pixels over the 50 cm FOV. The DECT-derived attenuation map was used for PET attenuation correction. Both noise-free and noisy attenuation corrected PET images were reconstructed using weighted OSEM (16 subsets, 14 iterations) on a 128 × 128 pixel matrix over a 50 cm FOV. 2.E.3. Image analysis and comparison

The difference images between the true attenuation map and the single-kVp or DECT derived attenuation map, and

012501-6

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

012501-6

F IG . 3. Data flow for evaluating use of DECT for PET attenuation correction.

the corresponding noise-free and noisy PET image were compared for the NCAT phantom. For each of the single-kVp and DECT derived attenuation map, six 2D ROIs were drawn on the attenuation map in locations of interest. Six corresponding ROIs were also drawn on the reconstructed noise-free and noisy PET images for analysis. The total root-mean-square errors (RMSE) within the spatial extent of the object of the attenuation and attenuation corrected PET images for the NCAT phantom were also calculated. A horizontal profile through the center of the noisefree PET images containing the lesions with contrast and in bone was drawn for comparative analysis.

estimated as

3. RESULTS

3.C. Evaluation of energy-dependent noise and bias

3.A. Theoretical prediction

The normalized COV for the central bin of the synthesized attenuation sinogram as a function of energy for the 10 cm

The theoretical prediction [Eq. (1)] for the central sinogram bin for the 10 cm diameter water cylinder is shown as one of the curves in Fig. 4, where the COV is normalized to a minimum value of 1. It can be observed from Fig. 4 that with increasing energy the normalized coefficient of variation first decreases to reach its expected minimum value at an energy point between the mean energies of the high and low transmitted spectra. Then the COV increases until it plateaus near 160 keV, at a level of approximately three times the minimum relative COV, and then remains almost constant with increasing energy.

Apl = 0.0272 + 24.4142TH − 17.2192TL + 3.3312TH TL − 0.3397TH2 − 2.5445TL2 + 0.0001TH2 TL2 − 0.2267TH3 + 0.1373TL3 ,

(11)

AAl = −0.0088 − 5.1178TH + 4.7601TL − 1.7392TH TL + 0.4399TH2 + 1.1547TL2 − 0.0003TH2 TL2 + 0.0719TH3 − 0.0422TL3 .

(12)

3.B. Dual-energy basis material decomposition

Before the added copper filter (Cu, Z = 64, 0.5 mm), the mean energy of the 80 and 140 kVp spectra were 42.7 and 59.7 keV, respectively. After the added filter, the mean energy of the 80 and 140 kVp spectra were 57.1 and 77.0 keV, respectively. The coefficients {di }, {ej } describing the relationships between aluminum and polyethylene thickness and the highand low-kVp log-transformed projections for the central bin without sinogram smoothing [Eqs. (6) and (7)] were Medical Physics, Vol. 41, No. 1, January 2014

F IG . 4. Normalized noise vs synthesized energy of monoenergetic attenuation data. The trend of decreasing linear attenuation coefficient of water with increased energy is also shown in the plot (the vertical axis at right side). Dual energy CT was acquired at 80 kVp with 256 mA and 140 kVp with 110 mA. Further details are given in the text.

012501-7

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

012501-7

diameter water cylinder is shown in Fig. 4 and agrees closely with the theoretical prediction. The normalized COV for the reconstructed attenuation map of the 10 cm diameter water cylinder as a function of energy is also shown for comparison. In this case the trend is similar to the theoretical prediction and sinogram domain simulations, but the exact amplification factor was slightly lower. Also shown in Fig. 4 is the normalized COV for the reconstructed attenuation map for the elliptical water cylinder. It can be seen that the image domain normalized COV curves for the elliptical and 10 cm circular water phantoms are similar if no sinogram smoothing is applied. For all the combinations of tube currents for high and low kVp scans, the percentage bias for both phantoms as a function of synthesized energy were very small, less than 2.4% for diagnostic CT energy range, and 0.4% for SPECT and PET energies. 3.D. DECT noise suppression and dose minimization for attenuation correction

The total radiation dose (mGy) to the 20 × 30 cm diameter water phantom as a function of tube currents of the high and low kVp scans are shown in Fig. 5. As expected, with increasing tube currents, the radiation dose increased. Figure 6 shows the averaged COV (unitless) of the 18 2D ROIs in the reconstructed attenuation map with sinogram smoothing as a function of tube currents for high and low kVp scans for the elliptical water phantom. The values acquired with DECT tube currents corresponding to give the same radiation dose as that of a single kVp scan of 120 kVp, 50 mAs (1.63 mGy) are also highlighted. It can be observed that, in general, COV vs tube currents in Fig. 6 has an inverse rela-

F IG . 6. Isocontour plot of coefficient of variation of DECT-derived attenuation map of elliptical water cylinder at 511 keV (PET photon energy) acquired with different CT tube currents at high and low kVps, respectively. The values obtained with DECT that have the same radiation dose as that of a single CT scan of 120 kVp, 50 mAs are shown with the solid curve copied from Fig. 5.

tionship with that of dose vs tube currents in Fig. 5. Based on Fig. 6, one can estimate the optimal combination of the tube currents that gave the same radiation dose as a single kVp CT scan, but generate an attenuation map at 511 keV (PET photon energy) with the minimum COV, which are 110 mA (55 mAs) for the 140 kVp scan and 256 mA for the 80 kVp scan. Note that the dual energy simulation included more spectral filtration (+0.5 mm Cu) than the single energy study. At the optimal tube current combination, the image domain normalized COV vs synthesized energy for the elliptical water cylinder with sinogram smoothing is also shown in Fig. 4 for comparison, which shows a reduction in noise at higher photon energies. The trend of normalized COV vs energy for sinogram smoothing case is the same as that without sinogram smoothing case, but the noise amplification factor is reduced from 2.4 to 1.9 at SPECT and PET photon energies.

3.E. Comparison of PET attenuation correction with DECT and single kVp-derived attenuation map

F IG . 5. Isocontour plot of the radiation dose (mGy) to a 20 × 30 cm diameter water phantom acquired as a function of mAs for high vs low energy (kVp) scans. The tube current levels that result in the same radiation dose as a low-dose single CT scan (120 kVp, 50 mAs, 1.63 mGy) are shown as a solid curve. Medical Physics, Vol. 41, No. 1, January 2014

Both DECT and single kVp derived attenuation maps were generated with roughly comparable radiation doses. At fixed optimal tube currents for high and low kVp CT scans, Fig. 7 shows the comparison of attenuation maps of the NCAT phantom derived from DECT or single-kVp scans. Table II lists the bias in each ROI and the total object RMSE for the two attenuation maps. Both PET attenuation maps have bias in the contrast agent region, 5.4% for the DECT-derived attenuation map compared to 11.0% in the single-kVp-based map. The total normalized RMSE for the object support region in the attenuation map is 2.6% with single-kVp-based estimation and to 1.9% with DECT method.

012501-8

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

F IG . 7. Comparison of attenuation maps at 511 keV.

The corresponding reconstructed PET images are shown in Fig. 8. The first row of Fig. 8 shows (from left to right) the reconstructed noise-free PET image with perfect attenuation correction, the noisy PET image attenuation corrected with the DECT method, and the one attenuation corrected with single-kVp method. The second row of Fig. 8 shows the difference image between the reconstructed PET images and the noise-free truth. The ROIs and profile for comparison are also illustrated in Fig. 8. The ROI quantitation results are compared in Table III. Figure 9 shows the comparison of profiles in the reconstructed noise-free PET images attenuation corrected with DECT and single-kVp derived methods. It can be clearly observed that the PET image attenuation corrected with DECT method slightly reduced RMSE for the total image and reduced the bias in contrast agent by a factor of ∼2.6. 4. DISCUSSION Standard methods of CT-based attenuation correction for SPECT and PET have been shown to work well for low density biological objects.9, 16 In cases where there is PET or SPECT tracer uptake in bone or other confounding highZ materials (e.g., contrast agent, implants) more accurate methods such as dual-energy CT-based attenuation correction may provide significantly more accurate images.19, 25 DECT is known to be problematic because of the noise amplification and the additional patient radiation dose required to perform two scans and to reduce noise. To analyze the properties of DECT-based attenuation correction we focused on the noise and bias of synthesized monoenergetic attenuation maps

012501-8

F IG . 8. Comparison of PET images using DECT and single-kVp CT-based AC methods.

as a function of the synthesized energy and source currents. Then we studied DECT noise suppression and dose minimization for attenuation correction. Finally, we used the NCAT phantom with iodine contrast to compare the noise and bias properties in attenuation corrected PET images. There are different strategies for dual energy CT imaging, such as separate dual kVp acquisitions, fast-switching interlaced dual kVp acquistions, simultaneous dual source acquisitions or multienergy detector-based. Our evaluations were based on dual kVp acquisitions. Considering the proposed attenuation correction is based on dual-energy synthesized attenuation maps (as shown in Fig. 3), the general methodology could be applied to other dual energy acquisition strategies. We used 0.5 mm Cu as additional flat filters for both high and low kVp scan, since our previous study indicated that this shapes the spectra to be more dose efficient.38 Although there are advantages for using different filtration for the two sources and a dual-tube scanner, we used the same filtration for both kVp scans as a single-tube CT scanner is more likely to be combined with a SPECT or PET scanner. Figure 4 shows the normalized coefficient of variation vs synthesized energy for water phantoms in the sinogram and image domains. The normalized coefficient of variation (COV) provides relevant information on the energydependent noise level, considering that the linear attenuation coefficient of water is energy dependent. The theoretical prediction was performed in the sinogram domain, and it closely matches the simulation results of the 10 cm diameter water cylinder. It should be noted that this close match is independent of whether COV is normalized or not, and only the central bin is used for simulation of the sinogram COV. The normalized COV vs energy trend could vary among different

TABLE II. Comparison of bias and RMSE in the DECT and SCT derived attenuation maps at 511 keV for the NCAT phantom. The RMSE was normalized to the true mean value within the spatial extent of the object.

DECT bias Single-kVp bias

ROI 1 Contrast (hot)

ROI 2 Soft tissue (hot)

ROI 3 Bone (hot)

ROI 4 Spine (normal)

ROI 5 Background (normal)

ROI 6 Lung (normal)

Total Object RMSE

5.4% 11.0%

− 0.85% − 2.5%

− 0.69% 0.66%

0.46% − 1.5%

− 1.3% − 2.0%

2.5% 11.5%

1.9% 2.6%

Medical Physics, Vol. 41, No. 1, January 2014

012501-9

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

012501-9

TABLE III. Comparison of noisy PET images with DECT and single kVp-based attenuation correction (AC) with different ROIs in the NCAT phantom (OSEM reconstruction, 16 subsets, 14 iterations). The RMSE was normalized to the true mean value within the spatial extent of the object.

PET image truth value Bias in PET with DECT AC Bias in PET with single-kVp AC

ROI 1 Contrast (hot)

ROI 2 Soft tissue (hot)

ROI 3 Bone (hot)

ROI 4 Spine (normal)

ROI 5 Background (normal)

ROI 6 Lung (normal)

Total Object RMSE

6.0 6.9% 17.8%

5.9 − 1.2% − 5.4%

5.8 − 3.3% − 1.1%

1.0 0.61% − 16.7%

1.0 − 1.2% − 4.3%

0.5 11.7% 42.9%

... 21.0% 23.3%

channels due to the transmitted spectra differences, but the general trend would be similar. The theoretical prediction and sinogram simulation do not agree completely, possibly due to the monoenergetic approximations in the theoretical prediction. Figure 4 also shows that the optimal synthesized energy where the COV reaches its minimum is object dependent. The energy for the minimum point in the curves for the two water phantoms were slightly different, likely due to the beam hardening effects for different thickness of the objects. In general, however, the minimal COV occurs when the energy is between the mean energy of the transmitted spectra. At higher energies, Fig. 4 indicates that the normalized COV increases with increasing energy until it plateaus around 160 keV (i.e., near SPECT energies), and then remains almost constant with increasing energy till up to the PET energy of 511 keV. At the PET energy of 511 keV the normalized COV of the attenuation maps for both water phantoms increased by a factor of approximately 2.4 compared to the minimal COV. Similar results comparing noise properties at 70 and 511 keV were also found by Stenner et al.40 We then evaluated CT sinogram smoothing as a method to further improve the noise/dose trade-offs. This was motivated by the requirement to smooth the CT image to match the PET resolution for CT-based attenuation correction to avoid artifacts. This smoothing could be performed prior to, during, or after CT image reconstruction. Due to the noise correlations that are introduced, we believe that smoothing should be performed either during or before reconstruction. We have previously evaluated the impact of smoothing the

F IG . 9. Profile comparison of the reconstructed noise-free PET images (OSEM reconstruction, 16 subsets, 14 iterations). Medical Physics, Vol. 41, No. 1, January 2014

raw sinogram with simple boxcar filters, both symmetric and non-symmetric, on the attenuation map and attenuation corrected PET images. For the given size of sinogram used in this study (888 channels, 984 views), we found that 5 × 5 boxcar smoothing followed by application of an adaptive trimmed mean filter35 before logarithmic conversion was a reasonable choice to provide substantial noise reduction while not introducing resolution mismatch artifacts.37 This is similar to the approach proposed by Colsher et al.36 Using a 2D 5 × 5 boxcar followed by adaptive trimmed mean filter for sinogram smoothing reduced the noise amplification factor from 2.4 to 1.9. It is possible to use statistically principled DECT reconstruction methods24 to allow for further improvements in the trade-off between radiation dose and image noise in DECT. It is important to note, however, that the main image quality metric is not CT image noise, but rather the quantitative accuracy of the PET image, which is determined by both noise and bias. For the same 10 cm diameter water cylinder, the maximum normalized COV in the image domain is lower than that in the central bin of the sinogram domain, probably due to noise correlations in image reconstruction, type of image reconstruction and kernel filters used, and bowtie filter effects across the total FOV. Figure 4 only shows one representative tube current combination for high and low kVp scans, and similar results could be obtained for other tube current combinations. As a check, the bias of the linear attenuation coefficients of the synthesized monoenergetic images were within 2.4% of the known true values across the entire energy range for various tube current combinations. Figures 5 and 6 showed the results of radiation dose as a function of tube currents and the tube current settings optimized for DECT-based PET attenuation correction. The dose estimation and parameter optimization were based on the elliptical water phantom. The study on the NCAT phantombased PET quantitation used the optimal parameters suggested from Fig. 6. Though it would be more precise if the dose calculation and parameter optimization were based on the NCAT phantom, the estimation was based on the elliptical water phantom as a more readily estimable surrogate. Figures 7–9 showed that with the same CT radiation dose as single energy CT, DECT acquisition parameters can be optimized to generate a more accurate attenuation map (less bias) leading to reduced PET bias in regions containing CT contrast agent and to similar or reduced RMSE for the total PET image. By comparing Tables II and III it can been seen that single-kVp-based attenuation correction leads to an overestimation of attenuation coefficients by 11.0% and

012501-10

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

the corresponding PET values by 17.8% in areas with contrast enhancement. With DECT-based attenuation correction, the overestimation of attenuation coefficients was reduced to 5.4% and the corresponding PET values to 6.9% for contrast agent. In practice, the single kVp-based attenuation correction method could be modified to reduce PET bias in iodine contrast. However, it would lead to increased bias for bone. In short, the bilinear method is problematic if both bone and iodine enhanced-regions are required to be quantitatively accurate in PET images. In the estimation of the coefficients for the DECT material basis decomposition formula (equations 6-7 and 112), non-zero constant terms have been proposed.33 However, it is reasonable to assume that Aα and Aβ should be exactly zero when TH and TL are both zero. We evaluated both cases and found that while the error of synthesized attenuation map at 511 keV increased in some regions and decreased in other regions, the overall root-mean-square error remained unchanged. The DECT results here are based on the basis material pair {polyethylene, aluminum}, however there are other potential choices of basis material pairs.18, 41 The bias for the DECT-based attenuation maps would ideally be zero in Table II. It has been shown, however, that such bias is a function of the choice of basis materials.41 For example, we repeated the DECT study with a {water, iodine} material pair (data not shown). In this case the bias in contrast agent region of the attenuation map was reduced from 5.4% to −1.6%. However, the bias in the lung region increased from 2.5% to 8.4%. These errors could be due to the strong K-edge of iodine at 33.2 keV in combination with the nonlinear estimation process. The optimal choice of basis material was outside the scope of this study. The findings in this study suggest that there is room for additional CT radiation dose reduction for both single- and dual-kVp-based attenuation correction in PET/CT imaging. We used additional flat Cu filtering to remove low-energy photons and shape the spectrum to be more dose-efficient for dual kVp CT scans, and this additional spectral shaping can be used for single-kVp CT scans as well.38 With the single and dual- kVp CT protocols in this study, Fig. 8 implies that the noise in the PET scans, not the CT techniques, continues to dominate the noise in the final PET images. The results here agree with our previous study38 that CT radiation dose could be further reduced by extending the standard techniques for PET/CT imaging with minimal effect on the PET noise and bias. The results presented here are based on simulations using a wide range of parameters to predict the general scanner behaviors for DECT-based PET attenuation correction. Before translation to clinical application, these results need to be confirmed in phantom-based physical experiments and patient studies. 5. CONCLUSION We addressed the energy dependent noise and bias properties for DECT derived attenuation correction for SPECT Medical Physics, Vol. 41, No. 1, January 2014

012501-10

and PET imaging. Our results showed the expected minimum noise value for a synthesized monoenergetic image at an energy between the mean energies of the two transmitted spectra. In addition it was found that for increasing energy, the coefficient of variation of water reached a plateau near 160 keV, at a level of approximately 2.4 times higher than the minimum coefficient of variation, and then remained constant with increasing energy. The trend of coefficient of variation vs synthesized energy provided confidence in the use of DECT for attenuation correction in PET/CT and SPECT/CT. Analysis of the noise, bias, and radiation dose indicated that there are optimal ranges for the mAs of the high and low kVp scans. We provided a solution for noise suppression and dose optimization, and showed that sinogram smoothing could further reduce the noise amplification in DECT-derived attenuation maps. With an appropriate selection of CT techniques, a DECT scan with the same radiation dose as a single CT scan can result in a PET image with improved quantitative accuracy. ACKNOWLEDGMENTS The authors acknowledge the support of General Electric for the use of CatSim, in particular Dr. Jed Pack, Dr. Bruno De Man and Dr. Paul Fitzgerald, and Mr. Steve Kohlmyer and Mr. Souma Sengupta. The authors also acknowledge many useful discussions with Dr. Jeffrey Fessler and Dr. Paul Segars. This work was supported by NIH Grants R01-CA160253 and R01-CA115870. a) Electronic

mail: [email protected] S. Rosenthal, J. Cullom, W. Hawkins, S. C. Moore, B. M. Tsui, and M. Yester, “Quantitative SPECT imaging: A review and recommendations by the Focus Committee of the Society of Nuclear Medicine Computer and Instrumentation Council,” J. Nucl. Med. 36, 1489–1513 (1995). 2 W. A. Weber, M. Schwaiger, and N. Avril, “Quantitative assessment of tumor metabolism using FDG-PET imaging,” Nucl. Med. Biol. 27, 683– 687 (2000). 3 W. A. Weber, “Assessing tumor response to therapy,” J. Nucl. Med. 50(Suppl 1), 1S–10S (2009). 4 W. Jentzen, L. Freudenberg, and A. Bockisch, “Quantitative imaging of (124)I with PET/CT in pretherapy lesion dosimetry. Effects impairing image quantification and their corrections,” Q. J. Nucl. Med. Mol. Imaging 55, 21–43 (2011). 5 E. P. Ficaro and J. R. Corbett, “Advances in quantitative perfusion SPECT imaging,” J. Nucl. Cardiol. 11, 62–70 (2004). 6 S. R. Cherry, J. A. Sorenson, and M. E. Phelps, Physics in Nuclear Medicine, 3rd ed. (Saunders, Philadelphia, PA, 2003). 7 K. J. Lacroix, B. M. W. Tsui, B. H. Hasegawa, and J. K. Brown, “Investigation of the use of x-ray Ct images for attenuation compensation in spect,” IEEE Trans. Nucl. Sci. 41, 2793–2799 (1994). 8 P. E. Kinahan, D. W. Townsend, T. Beyer, and D. Sashin, “Attenuation correction for a combined 3D PET/CT scanner,” Med. Phys. 25, 2046– 2053 (1998). 9 P. E. Kinahan, B. H. Hasegawa, and T. Beyer, “X-ray-based attenuation correction for positron emission tomography/computed tomography scanners,” Semin. Nucl. Med. 33, 166–179 (2003). 10 T. Beyer, D. W. Townsend, T. Brun, P. E. Kinahan, M. Charron, R. Roddy, J. Jerin, J. Young, L. Byars, and R. Nutt, “A combined PET/CT scanner for clinical oncology,” J. Nucl. Med. 41, 1369–1379 (2000). 11 P. J. Robinson and L. Kreel, “Pulmonary tissue attenuation with computedtomography - Comparison of inspiration and expiration scans,” J. Comput. Assist. Tomogr. 3, 740–748 (1979). 1 M.

012501-11

Xia, Alessio, and Kinahan: DECT attenuation correction for PET/CT

12 S.

C. Blankespoor, X. Wu, K. Kalki, J. K. Brown, H. R. Tang, C. E. Cann and B. H. Hasegawa, “Attenuation correction of SPECT using X-ray CT on an emission-transmission CT system: Myocardial perfusion assessment,” IEEE Trans. Nucl. Sci. 43, 2263–2274 (1996). 13 M. Abella, A. M. Alessio, D. A. Mankoff, L. R. Macdonald, J. J. Vaquero, M. Desco, and P. E. Kinahan, “Accuracy of CT-based attenuation correction in PET/CT bone imaging,” Phys. Med. Biol. 57, 2477–2490 (2012). 14 G. Antoch, L. S. Freudenberg, T. Egelhof, J. Stattaus, W. Jentzen, J. F. Debatin, and A. Bockisch, “Focal tracer uptake: A potential artifact in contrast-enhanced dual-modality PET/CT scans,” J. Nucl. Med. 43, 1339– 1342 (2002). 15 B. Teimourian, M. R. Ay, M. S. Zafarghandi, P. Ghafarian, H. Ghadiri, and H. Zaidi, “A novel energy mapping approach for CT-based attenuation correction in PET,” Med. Phys. 39, 2078–2089 (2012). 16 M. R. Ay, M. Shirmohammad, S. Sarkar, A. Rahmim, and H. Zaidi, “Comparative assessment of energy-mapping approaches in CT-based attenuation correction for PET,” Mol. Imaging Biol. 13, 187–198 (2011). 17 R. E. Alvarez and A. Macovski, “Energy-selective reconstructions in x-ray computerized tomography,” Phys. Med. Biol. 21, 733–744 (1976). 18 L. A. Lehmann, R. E. Alvarez, A. Macovski, W. R. Brody, N. J. Pelc, S. J. Riederer, and A. L. Hall, “Generalized image combinations in dual KVP digital radiography,” Med. Phys. 8, 659–667 (1981). 19 B. H. Hasegawa, T. F. Lang, J. K. Brown, E. L. Gingold, S. M. Reilly, S. C. Blankespoor, S. C. Liew, B. M. W. Tsui, and C. Ramanathan, “Objectspecific attenuation correction of spect with correlated dual-energy x-ray Ct,” IEEE Trans. Nucl. Sci. 40, 1242–1252 (1993). 20 M. J. Guy, I. A. Castellano-Smith, M. A. Flower, G. D. Flux, R. J. Ott, and D. Visvikis, “DETECT - Dual energy transmission estimation CT - for improved attenuation correction in SPECT and PET,” IEEE Trans. Nucl. Sci. 45, 1261–1267 (1998). 21 W. A. Kalender, E. Klotz, and L. Kostaridou, “An algorithm for noise suppression in dual energy CT material density images,” IEEE Trans. Med. Imaging 7, 218–224 (1988). 22 R. E. E. Alvarez, “A comparison of noise and dose in conventional and energy selective computed tomography,” IEEE Trans. Nucl. Sci. 26, 2853– 2856 (1979). 23 A. M. Doost-Hoseini, “Signal processing in dual-energy computerized tomography (beam-filter),” Ph.D. dissertation, Stanford University, 1984. 24 J. Noh, J. A. Fessler, and P. E. Kinahan, “Statistical sinogram restoration in dual-energy CT for PET attenuation correction,” IEEE Trans. Med. Imaging 28, 1688–1702 (2009). 25 P. E. Kinahan, A. M. Alessio, and J. A. Fessler, “Dual energy CT attenuation correction methods for quantitative assessment of response to cancer therapy with PET/CT imaging,” Technol. Cancer Res. Treat. 5, 319–327 (2006). 26 B. De Man, S. Basu, N. Chandra, B. Dunham, P. Edic, M. Iatrou, S. Mcolash, P. Sainath, C. Shaughnessy, B. Tower, and E. Williams, “CATSIM: A new computer assisted tomography simulation environment,” Proc. SPIE 6510, 6510G (2007).

Medical Physics, Vol. 41, No. 1, January 2014

27 B.

012501-11

De Man, J. Nuyts, P. Dupont, G. Marchal, and P. Suetens, “An iterative maximum-likelihood polychromatic algorithm for CT,” IEEE Trans. Med. Imaging 20, 999–1008 (2001). 28 P. J. La Riviere, “Penalized-likelihood sinogram smoothing for low-dose CT,” Med. Phys. 32, 1676–1683 (2005). 29 B. R. Whiting, P. Massoumzadeh, O. A. Earl, J. A. O’Sullivan, D. L. Snyder, and J. F. Williamson, “Properties of preprocessed sinogram data in x-ray computed tomography,” Med. Phys. 33, 3290–3303 (2006). 30 C. Comtat, P. E. Kinahan, M. Defrise, C. Michel, C. Lartizien, and D. W. Townsend, “Simulating whole-body PET scanning with rapid analytical methods,” in Proceedings of IEEE Nuclear Science Symposium and Medical Imaging Conference, Seattle, WA (IEEE, New York, NY, 1999), Vol. 3, pp. 1260–1264. 31 W. P. Segars, “Development and application of the new dynamic NURBSbased cardiac torso (NCAT) phantom,” Ph.D. dissertation, The University of North Carolina, 2001. 32 W. P. Segars, M. Mahesh, T. J. Beck, E. C. Frey, and B. M. W. Tsui, “Realistic CT simulation using the 4D XCAT phantom,” Med. Phys. 35, 3800– 3808 (2008). 33 K.-S Chuang and H. K. Huang, “Comparison of four dual energy image decomposition methods,” Phys. Med. Biol. 33, 445–466 (1988). 34 T. Xia, A. M. Alessio, and P. E. Kinahan, “Noise and bias properties of monoenergetic images from DECT used for attenuation correction with PET/CT and SPECT/CT,” Proc. Soc. Photo-Opt. Instrum. Eng. 7622, 762225–762228 (2010). 35 J. Hsieh, “Adaptive streak artifact reduction in computed tomography resulting from excessive x-ray photon noise,” Med. Phys. 25, 2139–2147 (1998). 36 J. G. Colsher, H. Jiang, J.-B. Thibault, A. Lonn, T. Pan, S. J. Lokitz, and T. G. Turkington, “Ultra low dose CT for attenuation correction in PET/CT,” in Proceedings of IEEE Nuclear Science Symposium and Medical Imaging Conference, Dresden, Germany (IEEE, New York, NY, 2008), pp. 5506–5511. 37 T. Xia, A. Alessio, and P. E. Kinahan, “Limits of ultra-low dose CT attenuation correction for PET/CT,” in Proceedings of IEEE Nuclear Science Symposium and Medical Imaging Conference, Orlando, Florida (IEEE, New York, NY, 2009), pp. 3074–3079. 38 T. Xia, A. M. Alessio, B. D. Man, R. Manjeshwar, E. Asma, and P. E. Kinahan, “Ultra-low dose CT attenuation correction for PET/CT,” Phys. Med. Biol. 57, 309–328 (2012). 39 C. Burger, G. Goerres, S. Schoenes, A. Buck, A. H. R. Lonn, and G. K. von Schulthess, “PET attenuation coefficients from CT images: Experimental evaluation of the transformation of CT into PET 511-keV attenuation coefficients,” Eur. J. Nucl. Med. Mol Imaging 29, 922–927 (2002). 40 P. Stenner and M. Kachelriess, “Dual energy exposure control (DEEC) for computed tomography: Algorithm and simulation study,” Med. Phys. 35, 5054–5060 (2008). 41 K. L. Goh, S. C. Liew, and B. H. Hasegawa, “Correction of energydependent systematic errors in dual-energy x-ray CT using a basis material coefficients transformation method,” IEEE Trans. Nucl. Sci. 44, 2419– 2424 (1997).

CT.

The authors evaluate the energy dependent noise and bias properties of monoenergetic images synthesized from dual-energy CT (DECT) acquisitions. These...
1MB Sizes 0 Downloads 0 Views