High correlation between radiation dose estimates for 256-slice CT obtained by highly parallelized hybrid Monte Carlo computation and solid-state metal-oxide semiconductor field-effect transistor measurements in physical anthropomorphic phantoms.

PURPOSE: Accurate, patient-specific radiation dosimetry for CT scanning is critical to optimize radiation doses and balance dose against image quality. While Monte Carlo (MC) simulation is often used to estimate doses from CT, comparison of estimates to experimentally measured values is lacking for advanced CT scanners incorporating novel design features. We aimed to compare radiation dose estimates from MC simulation to doses measured in physical anthropomorphic phantoms using metal-oxide semiconductor field-effect transistors (MOSFETs) in a 256-slice CT scanner. METHODS: Fifty MOSFETs were placed in organs within tissue-equivalent anthropomorphic adult and pediatric radiographic phantoms, which were scanned using a variety of chest, cardiac, abdomen, brain, and whole-body protocols on a 256-slice system. MC computations were performed on voxelized CT reconstructions of the phantoms using a highly parallel MC tool developed specifically for diagnostic X-ray energies and rapid computation. Doses were compared between MC estimates and physical measurements. RESULTS: The average ratio of MOSFET to MC dose in the in-field region was close to 1 (range, 0.96-1.12; mean ± SD, 1.01 ± 0.04), indicating outstanding agreement between measured and simulated doses. The difference between measured and simulated doses tended to increase with distance from the in-field region. The error in the MC simulations due to the limited number of simulated photons was less than 1%. The errors in the MOSFET dose determinations in the in-field region for a single scan were mainly due to the calibration method and were typically about 6% (8% if the error in the reading of the ionization chamber that was used for the MOSFET calibration was included). CONCLUSIONS: Radiation dose estimation using a highly parallelized MC method is strongly correlated with experimental measurements in physical adult and infant anthropomorphic phantoms for a wide range of scans performed on a 256-slice CT scanner. Incorporation into CT scanners of radiation-dose distribution estimation, employing the scanner's reconstructed images of the patient, may offer the potential for accurate patient-specific CT dosimetry.

Calibration and error analysis of metal-oxide-semiconductor field-effect transistor dosimeters for computed tomography radiation dosimetry.

PURPOSE: Metal-oxide-semiconductor field-effect transistors (MOSFETs) serve as a helpful tool for organ radiation dosimetry and their use has grown in computed tomography (CT). While different approaches have been used for MOSFET calibration, those using the commonly available 100 mm pencil ionization chamber have not incorporated measurements performed throughout its length, and moreover, no previous work has rigorously evaluated the multiple sources of error involved in MOSFET calibration. In this paper, we propose a new MOSFET calibration approach to translate MOSFET voltage measurements into absorbed dose from CT, based on serial measurements performed throughout the length of a 100-mm ionization chamber, and perform an analysis of the errors of MOSFET voltage measurements and four sources of error in calibration. METHODS: MOSFET calibration was performed at two sites, to determine single calibration factors for tube potentials of 80, 100, and 120 kVp, using a 100-mm-long pencil ion chamber and a cylindrical computed tomography dose index (CTDI) phantom of 32 cm diameter. The dose profile along the 100-mm ion chamber axis was sampled in 5 mm intervals by nine MOSFETs in the nine holes of the CTDI phantom. Variance of the absorbed dose was modeled as a sum of the MOSFET voltage measurement variance and the calibration factor variance, the latter being comprised of three main subcomponents: ionization chamber reading variance, MOSFET-to-MOSFET variation and a contribution related to the fact that the average calibration factor of a few MOSFETs was used as an estimate for the average value of all MOSFETs. MOSFET voltage measurement error was estimated based on sets of repeated measurements. The calibration factor overall voltage measurement error was calculated from the above analysis. RESULTS: Calibration factors determined were close to those reported in the literature and by the manufacturer (~3 mV/mGy), ranging from 2.87 to 3.13 mV/mGy. The error σV of a MOSFET voltage measurement was shown to be proportional to the square root of the voltage V: σV=cV where c = 0.11 mV. A main contributor to the error in the calibration factor was the ionization chamber reading error with 5% error. The usage of a single calibration factor for all MOSFETs introduced an additional error of about 5-7%, depending on the number of MOSFETs that were used to determine the single calibration factor. The expected overall error in a high-dose region (~30 mGy) was estimated to be about 8%, compared to 6% when an individual MOSFET calibration was performed. For a low-dose region (~3 mGy), these values were 13% and 12%. CONCLUSIONS: A MOSFET calibration method was developed using a 100-mm pencil ion chamber and a CTDI phantom, accompanied by an absorbed dose error analysis reflecting multiple sources of measurement error. When using a single calibration factor, per tube potential, for different MOSFETs, only a small error was introduced into absorbed dose determinations, thus supporting the use of a single calibration factor for experiments involving many MOSFETs, such as those required to accurately estimate radiation effective dose.

A pilot study of patient-specific cardiovascular MDCT dose maps and their utility in estimating patient-specific organ and effective doses in obese patients.

BACKGROUND: Estimates of effective dose (E) for cardiovascular CT are obtained from a scanner-provided dose metric, the dose-length product (DLP), and a conversion factor. These estimates may not adequately represent the risk of a specific scan to obese adults. OBJECTIVE: Our objective was to create dose maps sensitive to patient size and anatomy in the irradiated region from a patient's own CT images and compare measured E (EDoseMap) to doses determined from standard DLP conversion (EDLP) in obese adults. METHODS: 21 obese patients (mean body mass index, 39 kg/m(2)) underwent CT of the pulmonary veins, thoracic aorta, or coronary arteries. DLP values were converted to E. A Monte Carlo tool was used to simulate X-ray photon interaction with virtual phantoms created from each patient's image set. Organ doses were determined from dose maps. EDoseMap was computed as a weighted sum of organ doses multiplied by tissue-weighting factors. RESULTS: EDLP (mean ± SD, 5.7 ± 3.3 mSv) was larger than EDoseMap (3.4 ± 2.4 mSv) (difference = 2.3; P < .001). CONCLUSION: Dose maps derived from patient CT images yielded lower effective doses than DLP conversion methods. Considering over all patient size, organ size, and tissue composition could lead to better dose metrics for obese patients.

Projection extension for region of interest imaging in cone-beam CT.

RATIONALE AND OBJECTIVES: For 3D X-ray imaging during interventions, changes of the imaged object are often restricted to a small part of the field of view, suggesting region of interest (ROI) imaging by irradiating this area only. In this article, we present a novel method for extension of truncated projections in order to avoid truncation artifacts in C-arm based 3D ROI imaging. MATERIALS AND METHODS: The method makes use of prior knowledge by combining forward projections of a previously acquired, nontruncated 3D reference image with the truncated ROI projections. Rigid registration between the two datasets is achieved by using a technique based on local cross-correlation. To account for a gray value mismatch between the two data sets due to, e.g., differing beam quality and different contributions of scattered radiation, a linear gray level transformation is applied to the forward-projected reference data. RESULTS: The performance of different gray value transformation schemes is systematically assessed by means of numerical simulations. For various simulated scenarios, the best performing transformation has been identified, providing practical guidelines for selecting a scheme depending on the origin of the gray-level mismatch. Experiments prove the high performance of the developed method. CONCLUSION: The presented technique enables almost artifact-free 3D ROI imaging during interventions. This actually allows for repeated scans at low dose and enables intraprocedural imaging of large objects even with a small detector. However, applicability of the method is limited to scenarios where direct access to a reference image, e.g., a prior CT scan, is available.

Performance simulation of an x-ray detector for spectral CT with combined Si and Cd[Zn]Te detection layers.

The most obvious problem in obtaining spectral information with energy-resolving photon counting detectors in clinical computed tomography (CT) is the huge x-ray flux present in conventional CT systems. At high tube voltages (e.g. 140 kVp), despite the beam shaper, this flux can be close to 109 Mcps mm⁻² in the direct beam or in regions behind the object, which are close to the direct beam. Without accepting the drawbacks of truncated reconstruction, i.e. estimating missing direct-beam projection data, a photon-counting energy-resolving detector has to be able to deal with such high count rates. Sub-structuring pixels into sub-pixels is not enough to reduce the count rate per pixel to values that today's direct converting Cd[Zn]Te material can cope with (≤ 10 Mcps in an optimistic view). Below 300 µm pixel pitch, x-ray cross-talk (Compton scatter and K-escape) and the effect of charge diffusion between pixels are problematic. By organising the detector in several different layers, the count rate can be further reduced. However this alone does not limit the count rates to the required level, since the high stopping power of the material becomes a disadvantage in the layered approach: a simple absorption calculation for 300 µm pixel pitch shows that the required layer thickness of below 10 Mcps/pixel for the top layers in the direct beam is significantly below 100 µm. In a horizontal multi-layer detector, such thin layers are very difficult to manufacture due to the brittleness of Cd[Zn]Te. In a vertical configuration (also called edge-on illumination (Ludqvist et al 2001 IEEE Trans. Nucl. Sci. 48 1530-6, Roessl et al 2008 IEEE NSS-MIC-RTSD 2008, Conf. Rec. Talk NM2-3)), bonding of the readout electronics (with pixel pitches below 100 µm) is not straightforward although it has already been done successfully (Pellegrini et al 2004 IEEE NSS MIC 2004 pp 2104-9). Obviously, for the top detector layers, materials with lower stopping power would be advantageous. The possible choices are, however, quite limited, since only 'mature' materials, which operate at room temperature and can be manufactured reliably should reasonably be considered. Since GaAs is still known to cause reliability problems, the simplest choice is Si, however with the drawback of strong Compton scatter which can cause considerable inter-pixel cross-talk. To investigate the potential and the problems of Si in a multi-layer detector, in this paper the combination of top detector layers made of Si with lower layers made of Cd[Zn]Te is studied by using Monte Carlo simulated detector responses. It is found that the inter-pixel cross-talk due to Compton scatter is indeed very high; however, with an appropriate cross-talk correction scheme, which is also described, the negative effects of cross-talk are shown to be removed to a very large extent.

Directional view interpolation for compensation of sparse angular sampling in cone-beam CT.

In flat detector cone-beam computed tomography and related applications, sparse angular sampling frequently leads to characteristic streak artifacts. To overcome this problem, it has been suggested to generate additional views by means of interpolation. The practicality of this approach is investigated in combination with a dedicated method for angular interpolation of 3-D sinogram data. For this purpose, a novel dedicated shape-driven directional interpolation algorithm based on a structure tensor approach is developed. Quantitative evaluation shows that this method clearly outperforms conventional scene-based interpolation schemes. Furthermore, the image quality trade-offs associated with the use of interpolated intermediate views are systematically evaluated for simulated and clinical cone-beam computed tomography data sets of the human head. It is found that utilization of directionally interpolated views significantly reduces streak artifacts and noise, at the expense of small introduced image blur.