PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Comput Assist Tomogr. Author manuscript; available in PMC 2010 April 27.
Published in final edited form as:
PMCID: PMC2860287
NIHMSID: NIHMS192596

Demonstration of Dose and Scatter Reductions for Interior Computed Tomography

Deepak Bharkhada, MS,*§ Hengyong Yu, PhD,|| Robert Dixon, PhD,|| Yuchuan Wei, PhD,*§ J. Jeffrey Carr, MD,||# J. Daniel Bourland, PhD,§ Ryan Best, Richard Hogan,|| and Ge Wang, PhD*§||

Abstract

With continuing developments in computed tomography (CT) technology and its increasing use of CT imaging, the ionizing radiation dose from CT is becoming a major public concern particularly for high-dose applications such as cardiac imaging. We recently proposed a novel interior tomography approach for x-ray dose reduction that is very different from all the previously proposed methods. Our method only uses the projection data for the rays passing through the desired region of interest. This method not only reduces x-ray dose but scatter as well. In this paper, we quantify the reduction in the amount of x-ray dose and scattered radiation that could be achieved using this method. Results indicate that interior tomography may reduce the x-ray dose by 18% to 58% and scatter to the detectors by 19% to 59% as the FOV is reduced from 50 to 8.6 cm.

Keywords: computed tomography (CT), cardiac CT, interior tomography, Farmer chamber, pencil chamber, x-ray dose reduction, x-ray scatter reduction, CTDI

Tremendous progress in computed tomography (CT) technology such as the introduction of multi-slice CT and subsecond rotation times have increased the clinical uses of CT. The number of CT scans per year in the United States has increased from 3 million in 1980 to 62 million in 2006.1,2 Computed tomography is among the largest source of radiation exposure in radiology. It has been estimated that although CT studies constitute only 4% of all radiological procedures, they account for 40% of the radiation dose delivered.3 Cancer risk associated with CT is now becoming a growing concern. It has been estimated that 1.5% to 2% of all cancer cases may be associated with CT.1 Cancer risk is especially a concern for high-dose applications such as cardiac CT.4 Computed tomography technology is rapidly evolving, and multi-detector CT scanners with up to 320 slices have now become available. Cardiac CT imaging is also progressing considerably, with technologies such as controlled cardiac CT5,6 and dynamic volumetric cardiac CT7,8 and scanning trajectories such as saddle curve9 and composite circling.10 Hence, dramatic progress in CT technology, especially with multi-detector CT, is foreseeable in the near future.11 As a result, future CT imaging procedures are likely to improve the diagnosis of cardiovascular disease, including coronary artery disease, risk factor assessment, and other cardiovascular diseases such as congenital heart disease in infants and children. Therefore, it is necessary to find ways for effective dose reduction, especially for high-dose applications like cardiac CT.

Most of the proposed dose-reduction strategies have focused on modulation of scan parameters, modification of scanner geometry, and use of prospective gating.12,13 Although these methods can reduce dose, significant additional dose reduction can be achieved by developing acquisition and corresponding reconstruction methods that require less raw projection data. The first landmark work along this direction is the well-known fan-beam half-scan formula.14 Since then, several imaging conditions have been proposed that have enabled further reduction in the amount of required projection data.1517 However, these conditions do not apply to the case where the x-ray beam is entirely limited to the region of interest (ROI) to be reconstructed. Such a condition could reduce the dose significantly for interior imaging applications such as cardiac CT. Last year, our group18 proposed an improved condition in this direction. As illustrated in Figure 1, we proved that the image on an interior line segment can be exactly reconstructed as long as it is within the scan field of view and x-ray attenuation coefficients are known for some portion of the segment. This line segment can be completely within the object, and the projection data necessary for reconstruction are fully truncated. The interior tomography method was also independently proposed by Kudo et al.19,20 Figure 2 compares the difference between the data acquisition configurations of conventional CT and interior tomography. The filtered back-projection (FBP) algorithm used in current CT scanners requires the object to be completely within x-ray beam in all views, but interior tomography only needs the ROI to be within the x-ray beam in all views. Owing to reduced x-ray scan field of view (SFOV), the interior tomography method not only decreases the x-ray dose but also the x-ray scatter.

FIGURE 1
Exact interior reconstruction settings. A, One known region adjacent to the reconstructible region. B, Two known regions adjacent to the reconstructible region, where f(x) is the image and g(x) the Hilbert transform.
FIGURE 2
Difference between conventional CT and interior tomography. A, Conventional CT with the x-ray beam containing the object support. B, Interior tomography with the x-ray beam only covering the ROI, where 1 is for an object support, 2 for an ROI, and 3 for ...

Because our interior tomography approach is for interior imaging and uses localized projection data, the external markers on the human body cannot serve as known regions for reconstruction. We recently tested the feasibility of using blood intensity in aorta and/or bone intensity in vertebra as the known information for cardiac imaging applications and obtained very promising results.21 Figure 3 compares an FBP reconstruction image using nontruncated projection data with an interior tomography reconstruction image using only projection data for the rays passing through the ROI. Details of the reconstruction method can be found in Bharkhada et al.21 The known information for blood can be obtained by scanning a blood sample. Alternatively, the known information can be obtained from existing previous CT scans or from low-dose CT scans obtained immediately before the interior tomography scanning session. We believe that our interior tomography method can significantly reduce x-ray dose and scatter for interior CT imaging applications such as cardiac and dental.

FIGURE 3
Comparison of reconstructed images for conventional CT and interior tomography. A, Original image showing the ROI. B, Image reconstructed using FBP reconstruction from nontruncated projection data set. C, Region of interest image reconstructed with interior ...

Radiation dose is described as the amount of energy deposited per unit mass, say, in a patient's body, as a result of ionizing radiation exposure. The SI unit of radiation dose is the gray, and 1 Gy = 1 J/kg (energy per mass). The commonly used radiation dose parameter in CT is the CT dose index (CTDI). In practice, CTDI100 (in Gy or mGy) is the measured parameter, using a pencil ionization chamber to integrate the radiation exposure of a single axial scan for a length of 100 mm.22 Recently, a novel and more practical method has been proposed and validated,23,24 which uses a 0.6-cm3 Farmer chamber (small volume ionization chamber) to measure CTDIL, the accumulated dose for any scan length including CTDI100. The Farmer ionization chamber is placed in a standard polymethyl methacrylate (PMMA) phantom of either 16 or 32 cm in diameter. CTDI100 is measured on the central and peripheral axes in the phantom. CTDIw is the weighted average of the CTDI100 measurements at the center and the peripheral locations of the phantom given as

CTDIw=23CTDI100(periphery)+13CTDI100(center).

CTDIw reflects the mean absorbed dose over the central plane (z = 0) of the phantom at the center of a 10-cm scan length for a pitch P = 1. Computed tomographic dose index volume (CTDIvol) islikewisethe mean doseon the central x-y plane at z = 0 for a pitch P and is derived from CTDIw as CTDIvol = p–1CTDIw.22 The expected dose reduction, using interior tomography, is determined from dose measurements using the Farmer chamber.

In the diagnostic x-ray energy range, the predominant source of scattered photons and electrons is through the Compton effect. In this interaction, a photon interacts with a loosely bound or free electron within the attenuating medium. This results in release of a recoil Compton electron and deflection of the incident photon from its original path with reduced energy.25 The scattered photons add to the signal along different x-ray paths and present themselves as noise in the reconstructed image. Because interior tomography requires reduced SFOV, the scattered photons that would have come from the regions outside the reduced SFOV but within the nominal SFOV would be absent, thus reducing the contributions of scattered radiation. It is well known that scatter is reduced with reduction in radiation field size. The expected reduction in scattered radiation is determined from scattered radiation measurements using the diagnostic pancake-style ionization chamber.

The main goal of this article was to quantify the x-ray dose and scatter radiation reduction achieved as a result of using our interior tomography method.

METHODS

The CTDI100 Acquisition Method Using a Pencil Chamber

Using a pencil chamber of active length [ell], a single axial rotation is made about its center with the phantom and table held stationary, measuring the integral of the single rotation axial dose profile over (2,2). For the determination of CTDI100, the active length is =100mm. The charge q collected is given by

q=1Nk22f(z)dz

where Nk is the chamber calibration factor in mGy/nC (as given by an accredited dosimetry calibration laboratory) and

CTDIL=1N×TL2L2f(z)dz

where N × T is the actual detector length. The chamber is vented, and pressure and temperature corrections are also applied.

CTDI100 Using Farmer Chamber

With the small chamber method, the accumulated dose at z = 0 is measured by integrating the current from the ionization chamber located at a fixed point in the phantom at the midpoint (z = 0) of the scanned length L = vt0, whereas the phantom and chamber are translated by the couch through the beam plane at velocity v during the helical acquisition of total time t0. Using the Farmer chamber guarantees that the scan length is always identical to the integration length L of the single rotation dose profile f(z). For the determination of CTDI100, this scan length must be 100 mm. Then, CTDI100 = pDL (0), where DL (0) is the dose measured by the Farmer chamber and p is an arbitrary acquisition pitch.24

Experimental Validation of Farmer Chamber

We measured the CTDI100 at the center and the periphery of the PMMA body phantom 32 cm in diameter. The phantom length is 25 cm, which is greater than the standard phantom of 14 cm length. The additional length provides better accuracy in the measurement of CTDI100. A 10-cm-long Victoreen 500-200 pencil ionization chamber (Victoreen Instruments, Mödling, Austria) and a NE 2505/3A (0.6 cm3) Farmer chamber (Nuclear Enterprises Ltd, Reading, England) are used to make the measurement. In a GE Lightspeed 8-slice CT scanner, a single axial rotation is used for the measurements with the pencil chamber, and a helical series is used for the Farmer chamber. The measurements are made at the same kVp and mA for both the chambers. Details for the scan protocols are presented in the results.

Interior Tomography Dose Measurements

The same PMMA cylindrical body phantom and the Farmer chamber NE 2505/3A (0.6 cm3) are used for interior tomography dose measurements. Helical trajectory is used for a GE Lightspeed 16-slice scanner with a large (50 cm) SFOV. To create small SFOVs and thus simulate the interior tomography method, collimator lead plates 0.0625 inch thick are placed in front of the standard xy CT collimator to produce various openings (field sizes). This results in small SFOVs while the large SFOV is still set on the CT scanner. Thus, the same bowtie filter is used for small or large SFOV x-ray beams. The measurements for each SFOV are made at 3 different voltages: 80, 100, and 120 kVp. To measure the SFOV size, we exposed Kodak X-omat V x-ray film (Kodak) at the isocenter using a stationary gantry. The SFOV was determined from the measured width on the blackened film.

Interior Tomography Scatter Measurements

Scatter measurement is also carried out using the PMMA body phantom. The measurements are made for the same kVp, mA, and SFOVs as used for dose estimation on the GE Lightspeed 16-slice scanner. However, these measurements are made using a large-volume (15 cm3) diagnostic pancake-style ionization chamber (Fluke Biomedical 96035B). This chamber is used because it has greater sensitive volume and produces more current/charge that is suitable for low-dose measurements. In addition, the pencil chamber is too long for scatter measurements. To measure the scattered radiation reaching the detectors, the pancake chamber was placed just outside the primary beam approximately 180 degrees from the x-ray source, slightly off the detector z-position less than 3 cm. A stationary x-ray beam was used. The charge q collected is multiplied by the chamber calibration factor (milligray per nanocoulomb) to obtain the scattered radiation dose in milligray.

RESULTS

To demonstrate that the CTDI100 measurements made with the Farmer chamber are the same as those obtained with the commonly used pencil chamber, CTDI100 was measured at both the peripheral and the central axes of the phantom with both ionization chambers, using CT irradiation parameters of 120 kVp, 300 mA, rotation τ = 1 seconds, and N × T = 8 × 1.25 = 10 mm. We also calculated CTDIvol from these measurements. An axial scan was used for measurements with the pencil chamber. A helical scan with pitch p = 0.875, table velocity

v=p×N×Tτ=8.75mmss

for a total exposure time of t0 = 11.5 seconds, and an acquisition scan length L = vt0 = 100.625 mm were used for measurements with the Farmer chamber. The accredited dosimetry calibration laboratory calibration factor for the pencil chamber is 2.735 mGy/nC, and the calibration factor for the Farmer chamber is 41.69 mGy/nC. Table 1 compares the measured CTDI100 and CTDIvol values using the pencil chamber and the Farmer chamber. It is easy to see that values for CTDI100 and CTDIvol measured with both chambers agree within 1%.

TABLE 1
Comparison Between Pencil Chamber and Farmer Chamber

To quantify the dose reduction achieved using the interior tomography method, we next measured CTDI100 values at the central and peripheral axes of the phantom with the Farmer chamber using a helical trajectory. CTDIvol values were calculated from these values at various kVp and for 5 different SFOVs, 1 large (50 cm) SFOV, and 4 small SFOVs (8.6, 13.4, 16.8, and 20.6 cm) using the extra lead plate collimation to define the 4 different small SFOV apertures. These SFOVs were derived from the fan beam widths measured as increased film density from 4 sets of x-ray films at the isocenter. The measured CTDI100 and CTDIvol values along with the percentage dose reduction for the small SFOVs as compared to the 50 cm SFOV are presented in Tables 24. From these tables, it can be seen that as the SFOV decreases, the percentage dose reduction increases. In addition, the percentage dose reduction is greater at the periphery than at the center as one would expect. Although the actual dose increases with increasing kVp, no significant differences are observed in percentage dose reduction for the same small SFOV. In other words, percentage dose reduction is a function of size of the SFOV but seems to be independent of the kVp.

TABLE 2
Percentage Reduction of CTDI100 Values and CTDIvol for Small SFOVs as Compared to the 50-cm SFOV at 80 kVp
TABLE 4
Percentage Reduction of CTDI100 Values and CTDIvol for Small SFOVs as Compared to the 50-cm SFOV at 120 kVp

To quantify the scatter reduction at the detectors achieved using the interior tomography method, we measured scattered radiation as previously described in the methods at the same SFOVs, kVp values, and mA. The calibration factor of the diagnostic pancake ionization chamber used is 0.1682 mGy/nC. The measured scattered radiation along with the percentage reduction for the small SFOVs as compared to the 50 cm SFOV are presented in Table 5. From this table, we can see that as the SFOV decreases, the percentage of scatter reduction increases. Although the actual scattered radiation increases with increasing kVp, no significant differences are observed in percentage reduction in scattered radiation for the same small SFOV. In other words, percentage reduction in scattered radiation is a function of size of the SFOV, but it seems to be independent of the kVp.

TABLE 5
Percentage Reduction in Scatter Radiation to the Detectors for Small SFOVs as Compared to the 50-cm SFOV at Approximately 180 Degrees from the X-Ray Source

DISCUSSION

Progress in CT technology has led to an enormous increase in the use of CT for noninvasive imaging. Popular high-dose CT applications such as cardiac imaging, while providing clinically meaningful data, have become a growing concern owing to the ionizing nature of x-rays and their potential risks to cause cancer. Several directions have been proposed to reduce CT x-ray dose. We believe that our recently proposed interior tomography method, which uses only the localized projection data to reconstruct the ROI, has a significant potential to reduce CT x-ray dose. As an added advantage, it can also reduce scattered radiation, which could provide improved reconstructed image quality. The reduction in x-ray dose and scattered radiation that could be achieved using our interior tomography method was quantified by measuring x-ray dose and scattered radiation for several SFOVs and by comparing them to the measurements for a large reference SFOV.

Our results indicate that as the x-ray beam SFOV is decreased, the dose at the center and the periphery also decrease, thus leading to an overall dose reduction. As expected, the dose reduction is greater at the periphery than at the center. Similarly, as the SFOV was decreased, the amount of scattered radiation also decreased. Although the actual dose and scattered radiation increase with the kVp, the percentage reduction in dose and scattered radiation seems to be a function of SFOV only and independent of kVp.

The measured dose and scattered radiation for a few different SFOVs give us an estimate of the x-ray dose and scatter radiation reduction potential of our interior tomography method. In addition, our measurements are limited to a circular SFOV that is adequate for cardiac imaging owing to the shape and location of the heart. For other interior applications, the SFOV could potentially be different in different views. A model for predicting x-ray dose and scattered radiation reduction using interior tomography is considered beyond the scope of this paper but will provide additional insights and capability to predict dose and scattered radiation reduction for any desired SFOV.

In conclusion, our interior tomography method has a great potential for x-ray dose reduction. The results indicate that it can reduce x-ray dose in the range of 18% to 58% and scatter in the range of 19% to 59%. Thus, we have demonstrated that it gives a significant dose reduction for interior applications such as cardiac CT and thus helps to make the use of CT safer. At the same time, reduced scatter radiation could provide improved reconstructed image quality.

TABLE 3
Percentage Reduction of CTDI100 Values and CTDIvol for Small SFOVs as Compared to the 50-cm SFOV at 100 kVp

ACKNOWLEDGMENTS

The authors thank Thomas Payne in the Radiology Department of Wake Forest University School of Medicine for his assistance with measurement of scattered radiation.

This work was partially supported by NIH/NIBIB grants (EB002667, EB004287, and EB007288A) and a grant from Toshiba Medical Research Institute USA, Inc.

REFERENCES

1. Brenner DJ, Hall EJ. Computed tomography—an increasing source of radiation exposure. N Engl J Med. 2007;357:2277–2284. [PubMed]
2. IMV . 2006 Market Summary Report. IMV Medical Information Division; Des Plains, IL: 2006.
3. Shrimpton PC, Edyvean S. CT scanner dosimetry. Br J Radiol. 1998;71:1. [PubMed]
4. Cody DD, Mahesh M. AAPM/RSNA physics tutorials AAPM/RSNA physics tutorial for residents: technologic advances in multidetector CT with a focus on cardiac imaging. Radiographics. 2007;27:1829–1837. [PubMed]
5. Bai EW, Wang C, Liu Y, et al. Controlled Cardiac Computed Tomography. Internation Journal of Biomedical Imaging, 2006. 2006:11. Article ID 12819.
6. Bharkhada D, et al. A comparative study on interpolation methods for controlled cardiac CT. Int J Imaging Syst Technol. 2007;17:91–98.
7. Bharkhada D, Yu HY, Wang G. Knowledge-based dynamic volumetric cardiac CT with saddle curve trajectory. J Comput Assist Tomogr. 2008;32:942–950. [PMC free article] [PubMed]
8. Wang G, Zhao S, Heuscher D. A knowledge-based cone-beam x-ray CT algorithm for dynamic volumetric cardiac imaging. Med Phys. 2002;29:1807–1822. [PubMed]
9. Pack JD, Noo F, Kudo H. Investigation of saddle trajectories for cardiac CT imaging in cone-beam geometry. Phys Med Biol. 2004;49:2317–2336. [PubMed]
10. Yu HY, Wang G. Cone-beam composite-circling scan and exact image reconstruction for a quasi-short object. Int J Biomed Imaging. 2007;2007:10. article ID 87319. [PMC free article] [PubMed]
11. Budoff M, Achenbach S, Blumenthal RS, et al. Assessment of coronary artery disease by cardiac computed tomography: a scientific statement from the American Heart Association Committee on Cardiovascular Imaging and Intervention, Council on Cardiovascular Radiology and Intervention, and Committee on Cardiac Imaging, Council on Clinical Cardiology. Circulation. 2006;114:1761–1791. [PubMed]
12. Kalra MK, Maher MM, Toth TL, et al. Strategies for CT radiation dose optimization. Radiology. 2004;230:619–628. [PubMed]
13. Paul JF, Abada HT. Strategies for reduction of radiation dose in cardiac multislice CT. Eur Radiol. 2007;17:2028–2037. [PubMed]
14. Parker DL. Optimal short scan convolution reconstruction for fan beam CT. Med Phys. 1982;9:254–257. [PubMed]
15. Defrise M, et al. Truncated Hilbert transform and image reconstruction from limited tomographic data. Inverse Problems. 2006;22:1037–1053.
16. Noo F, Clackdoyle R, Pack JD. A two-step Hilbert transform method for 2D image reconstruction. Phys Med Biol. 2004;49:3903–3923. [PubMed]
17. Noo F, Defrise M, Clackdoyle R, et al. Image reconstruction from fan-beam projections on less than a short scan. Phys Med Biol. 2002;47:2525–2546. [PubMed]
18. Ye YB, Yu HY, Wei YC, et al. A general local reconstruction approach based on truncated Hilbert transform. Int J Biomed Imaging. 2007;2007:8. Article ID 63634. [PMC free article] [PubMed]
19. Kudo H, Courdurier M, Noo F, et al. Tiny a priori knowledge solves the interior problem in computed tomography. Phys Med Biol. 2008;53:2207. [PubMed]
20. Courdurier M, Noo F, Defrise M, et al. Solving the interior problem of computed tomography using a priori knowledge. Inverse Problems. 2008;24:065001. [PMC free article] [PubMed]
21. Bharkhada D, Yu HY, Ge SP, et al. Cardiac CT radiation dose reduction using interior reconstruction algorithm with the aorta and vertebra as known information. J Comput Assist Tomogr. 2007 To appear. [PMC free article] [PubMed]
22. Morin RL, Gerber TC, McCollough CH. Radiation dose in computed tomography of the heart. Circulation. 2003;107:917–922. [PubMed]
23. Dixon RL. A new look at CT dose measurement: beyond CTDI. Med Phys. 2003;30:1272. [PubMed]
24. Dixon RL, Ballard AC. Experimental validation of a versatile system of CT dosimetry using a conventional ion chamber: beyond CTDI100. Med Phys. 2007;34:3399. [PubMed]
25. Hendee WR, Ritenour ER. Medical Imaging Physics. 4th ed Wiley-Liss; New York: 2002.