Radiation dose estimates based on Monte Carlo simulation for spiral breast computed tomography imaging in a large cohort of patients.

BACKGROUND: Spiral breast computed tomography (BCT) equipped with a photon-counting detector (PCD) is a new radiological modality allowing for the compression-free acquisition of high-resolution 3-D datasets of the breast. Optimized dose exposu04170/re setups according to breast size were previously proposed but could not effectively be applied in a clinical environment due to ambiguity in measuring breast size. PURPOSE: This study aims to report the standard radiation dose values in a large cohort of patients examined with BCT, and to provide a mathematical model to estimate radiation dose based on morphological features of the breast. METHODS: This retrospective study was conducted on 1657 BCT examinations acquired between 2018 and 2021 from 829 participants (57 ± 10 years, all female). Applying a dedicated breast tissue segmentation algorithm and Monte Carlo (MC) simulation, mean absorbed dose (MAD), mean glandular dose (MGD), mean skin dose (MSD), maximum glandular dose (maxGD), and maximum skin dose (maxSD) were calculated and related to morphological features such as breast volume, effective diameter, breast length, skin volume, and glandularity. Effective dose (ED) was calculated by applying the corresponding beam and tissue weighting factors, 1 Sv/Gy and 0.12 per breast. Relevant morphological features predicting dose values were identified based on the Spearman's rank correlation coefficient. Exponential or bi-exponential models predicting the dose values as a function of morphological features were fitted by using a non-linear least squares (LS) method. The models were validated by assessing R2 and residual standard error (RSE). RESULTS: The most relevant morphological features for radiation dose estimation were the breast volume (correlation coefficient: -0.8), diameter (-0.7), and length (-0.6). The glandularity presented a weak-positive correlation (0.4) with MGD and maxGD due to the inhomogeneous distribution of the glandularity and absorbed dose in the 3-D breast volume. The standard MGDs were calculated to be 7.3 ± 0.7, 6.5 ± 0.3, and 5.9 ± 0.3 mGy, MADs to 7.6 ± 0.8, 6.8 ± 0.3, and 6.2 ± 0.3 mGy, maxSDs to 19.9 ± 1.6, 19.5 ± 0.5, and 18.9 ± 0.5 mGy, and EDs to 0.88 ± 0.08, 0.78 ± 0.04, and 0.72 ± 0.04 mSv for small, medium, and large breasts with average breast lengths of 5.9 ± 1.6, 8.7 ± 1.3, and 12.2 ± 2.0 cm, respectively. The estimated glandularity - 23.1 ± 16.9, 12.5 ± 11.4, and 6.9 ± 7.3% from small to large breasts. The mathematical models were able to estimate the MAD, MGD, MSD, and maxSD as a function of each morphological feature with only upto 0.5 mGy RSE. CONCLUSION: We presented the typical morphological features and standard dose values according to the breast size acquired from a large patient cohort. We established radiation dose estimation models allowing accurate estimation of dose values including MGD with an acceptable RSE based on each of the easily measured morphological features of the breast. Clinicians could use the breast length to operate as a dosimetric alert of the scanner prior to a BCT scan. Radiation exposure for BCT was lower than diagnostic mammography (MG) and cone-beam breast CT (BCT).

Investigation of spectral performance for single-scan contrast-enhanced breast CT using photon-counting technology: A phantom study.

PURPOSE: Contrast-enhanced imaging of the breast is frequently used in breast MRI and has recently become more common in mammography. The purpose of this study was to make single-scan contrast-enhanced imaging feasible for photon-counting breast CT (pcBCT) and to assess the spectral performance of a pcBCT scanner by evaluating iodine maps and virtual non-contrast (VNC) images. METHODS: We optimized the settings of a pcBCT to maximize the signal-to-noise ratio between iodinated contrast agent and breast tissue. Therefore, an electronic energy threshold dividing the x-ray spectrum used into two energy bins was swept from 23.17 keV to 50.65 keV. Validation measurements were performed by placing syringes with contrast agent (2.5 mg/ml to 40 mg/ml) in phantoms with 7.5 cm and 12 cm in diameter. Images were acquired at different tube currents and reconstructed with 300 µm isotropic voxel size. Iodine maps and VNC images were generated using image-based material decomposition. Iodine concentrations and CT values were measured for each syringe and compared to the known concentrations and reference CT values. RESULTS: Maximal signal-to-noise ratios were found at a threshold position of 32.59 keV. Accurate iodine quantification (average root mean square error of 0.56 mg/ml) was possible down to a concentration of 2.5 mg/ml for all tube currents investigated. The enhancement has been sufficiently removed in the VNC images, so they can be interpreted as unenhanced CT images. Only minor changes of CT values compared to a conventional CT scan were observed. Noise was increased by the decomposition by a factor of 2.62 and 4.87 (7.5 cm and 12 cm phantoms) but did not compromise the accuracy of the iodine quantification. CONCLUSIONS: Accurate iodine quantification and generation of VNC images can be achieved using contrast-enhanced pcBCT from a single CT scan in the absence of temporal or spatial misalignment. Using iodine maps and VNC images, pcBCT has the potential to reduce dose, shorten examination and reading time, and to increase cancer detection rates.

Metal Artifact Reduction in X-ray Computed Tomography Using Computer-Aided Design Data of Implants as Prior Information.

OBJECTIVES: The performance of metal artifact reduction (MAR) methods in x-ray computed tomography (CT) suffers from incorrect identification of metallic implants in the artifact-affected volumetric images. The aim of this study was to investigate potential improvements of state-of-the-art MAR methods by using prior information on geometry and material of the implant. MATERIALS AND METHODS: The influence of a novel prior knowledge-based segmentation (PS) compared with threshold-based segmentation (TS) on 2 MAR methods (linear interpolation [LI] and normalized-MAR [NORMAR]) was investigated. The segmentation is the initial step of both MAR methods. Prior knowledge-based segmentation uses 3-dimensional registered computer-aided design (CAD) data as prior knowledge to estimate the correct position and orientation of the metallic objects. Threshold-based segmentation uses an adaptive threshold to identify metal. Subsequently, for LI and NORMAR, the selected voxels are projected into the raw data domain to mark metal areas. Attenuation values in these areas are replaced by different interpolation schemes followed by a second reconstruction. Finally, the previously selected metal voxels are replaced by the metal voxels determined by PS or TS in the initial reconstruction. First, we investigated in an elaborate phantom study if the knowledge of the exact implant shape extracted from the CAD data provided by the manufacturer of the implant can improve the MAR result. Second, the leg of a human cadaver was scanned using a clinical CT system before and after the implantation of an artificial knee joint. The results were compared regarding segmentation accuracy, CT number accuracy, and the restoration of distorted structures. RESULTS: The use of PS improved the efficacy of LI and NORMAR compared with TS. Artifacts caused by insufficient segmentation were reduced, and additional information was made available within the projection data. The estimation of the implant shape was more exact and not dependent on a threshold value. Consequently, the visibility of structures was improved when comparing the new approach to the standard method. This was further confirmed by improved CT value accuracy and reduced image noise. CONCLUSIONS: The PS approach based on prior implant information provides image quality which is superior to TS-based MAR, especially when the shape of the metallic implant is complex. The new approach can be useful for improving MAR methods and dose calculations within radiation therapy based on the MAR corrected CT images.

Dose Optimization for Computed Tomography Localizer Radiographs for Low-Dose Lung Computed Tomography Examinations.

INTRODUCTION: Recent studies have shown a substantial reduction of radiation dose from computed tomography (CT) scans down to 0.1 mSv for lung cancer screening and cardiac examinations, when applying optimization techniques. Hence, CT localizer radiographs (LRs) might now be considered a significant contributor to the total dose of the CT examination. We investigated in our study the potential for reducing dose of the LRs by adapting the patient-specific acquisition parameters of the LR. MATERIALS AND METHODS: Localizer radiographs covering the lungs were acquired on 2 clinical scanners (64 slices, conventional detector [CD]; 96 slices, fully integrated detector [ID]) for 3 semianthropomorphic phantoms, representing a slim, a normal, and an obese adult. Starting at 120-kV tube voltage and 250-mA current were reduced until the image quality of the LR, and thereby the accuracy of the automatic exposure control was compromised; this was defined as a deviation of measured attenuation values in the center of the LR of more than 5% from the reference values measured at the highest tube voltage and current. Subsequent Monte Carlo calculations on anthropomorphic phantoms were performed to calculate organ and effective dose values for the respective optimal settings. In addition, effective dose values normalized to CTDIvol for tube voltages ranging from 60 to 160 kV were determined for the different combinations of phantom sizes, sexes, and LR views to evaluate dose efficiency. RESULTS: For the CD scanner, the optimal LR settings depended strongly on phantom size. Higher tube voltage and current were necessary for the larger phantoms. The ID scanner showed uncompromised LR quality for all phantoms using the lowest possible tube voltage-tube current combination of 80 kV and 20 mA. Depending on patient size and LR direction, effective dose values for the optimal settings ranged from 6 to 53 µSv and 3 to 11 µSv for the CD and ID scanner, respectively. For the example of an anterior-posterior LR on a normal patient, using the optimal settings instead of the standard settings on the ID scanner reduced LR dose from 53 µSv (120 kV, 30 mA) to 10 µSv (80 kV, 20 mA). The simulations for the different tube voltages show that effective dose and CTDIvol behave similarly for different views and patient sizes. However, the tube voltage level itself impacts the relationship between CTDIvol and effective dose, by up to a factor of 2. DISCUSSION: Dose from LRs may contribute significantly to the total effective dose of low-dose CT examinations such as lung cancer screening. Optimal LR settings can reduce LR dose substantially, but adaptations have to consider scanner characteristics, detector technology, and patient size. Thus, for low-dose CT examinations, such as cardiac examinations and lung cancer screening, LR optimization may result in a significant dose reduction and thereby in a substantial reduction of total dose.

Technical feasibility proof for high-resolution low-dose photon-counting CT of the breast.

X-ray computed tomography (CT) has been proposed and evaluated multiple times as a potentially alternative method for breast imaging. All efforts shown so far have been criticized and partly disapproved because of their limited spatial resolution and higher patient dose when compared to mammography. Our concept for a dedicated breast CT (BCT) scanner therefore aimed at novel apparatus and detector design to provide high spatial resolution of about 100 µm and average glandular dose (AGD) levels of 5 mGy or below. Photon-counting technology was considered as a solution to reach these goals. The complete concept was previously evaluated and confirmed by simulations and basic experiments on laboratory setups. We here present measurements of dose, technical image quality parameters and surgical specimen results on such a scanner. For comparison purposes, the specimens were also imaged with digital mammography (DM) and breast tomosynthesis (BT) apparatus. Results show that photon-counting BCT (pcBCT) at 5 mGy AGD offers sufficiently high 3D spatial resolution for reliable detectability of calcifications and soft tissue delineation. KEY POINTS: • Photon-counting detector technology allows for spatial resolution better than 100 µm. • pcBCT allows for dose levels in the screening mammography range. • pcBCT provides the highest quality imaging of microcalcifications.

Performance of Photon-Counting Breast Computed Tomography, Digital Mammography, and Digital Breast Tomosynthesis in Evaluating Breast Specimens.

RATIONALE AND OBJECTIVES: This study compared a novel photon-counting breast computed tomography (pcBCT) system with digital mammography (DM) and digital breast tomosynthesis (DBT) systems. For this reason, surgical specimens were examined with all three techniques and rated by three observers. MATERIALS AND METHODS: A total of 30 surgical specimens were investigated with DM, DBT, and pcBCT; the associated images were shown to three experienced radiologists. Findings (22 microcalcifications and 23 mass lesions) were recorded and compared to the results of the pathological examination. Sensitivity and specificity for detection of microcalcifications and lesions were calculated and displayed using receiver operating characteristic curves. RESULTS: Sensitivity for microcalcifications was 82% for DM, 70% for DBT, and 85% for pcBCT. Specificity for microcalcifications was 71% for DM, 75% for DBT, and 83% for pcBCT. Sensitivity for lesions was 45% for DM, 62% for DBT, and 65% for pcBCT. Specificity for lesions was 76% for DM, 62% for DBT, and 76% for pcBCT. CONCLUSIONS: pcBCT showed a comparable or superior performance compared to the clinically approved DM and DBT systems. Mass lesion detectability can be increased further by the use of contrast media.

High Atomic Number Contrast Media Offer Potential for Radiation Dose Reduction in Contrast-Enhanced Computed Tomography.

OBJECTIVES: Spectral optimization of x-ray computed tomography (CT) has led to substantial radiation dose reduction in contrast-enhanced CT studies using standard iodinated contrast media. The purpose of this study was to analyze the potential for further dose reduction using high-atomic-number elements such as hafnium and tungsten. As in previous studies, spectra were determined for which the patient dose necessary to provide a given contrast-to-noise ratio (CNR) is minimized. MATERIALS AND METHODS: We used 2 different quasi-anthropomorphic phantoms representing the liver cross-section of a normal adult and an obese adult patient with the lateral widths of 360 and 460 mm and anterior-posterior heights of 200 and 300 mm, respectively. We simulated and measured on 2 different scanners with x-ray spectra from 80 to 140 kV and from 70 to 150 kV, respectively. We determined the contrast for iodine-, hafnium-, and tungsten-based contrast media, the noise, and 3-dimensional dose distributions at all available tube voltages by measurements and by simulations. The dose-weighted CNR was determined as optimization parameter. RESULTS: Simulations and measurements were in good agreement regarding their dependence on energy for all parameters investigated. Hafnium provided the best performance for normal and for obese patient phantoms, indicating a dose reduction potential of 30% for normal and 50% for obese patients at 120 kV compared with iodine; this advantage increased further with higher kV values. Dose-weighted CNR values for tungsten were always slightly below the hafnium results. Iodine proved to be the superior choice at voltage values of 80 kV and below. DISCUSSION: Hafnium and tungsten both seem to be candidates for contrast-medium-enhanced CT of normal and obese adult patients with strongly reduced radiation dose at unimpaired image quality. Computed tomography examinations of obese patients will decrease in dose for higher kV values.

Comparative investigation of the detective quantum efficiency of direct and indirect conversion detector technologies in dedicated breast CT.

PURPOSE: To investigate the dose saving potential of direct-converting CdTe photon-counting detector technology for dedicated breast CT. MATERIALS AND METHODS: We analyzed the modulation transfer function (MTF), the noise power spectrum (NPS) and the detective quantum efficiency (DQE) of two detector technologies, suitable for breast CT (BCT): a flat-panel energy-integrating detector with a 70 µm and a 208 µm thick gadolinium oxysulfide (GOS) and a 150 µm thick cesium iodide (CsI) scintillator and a photon-counting detector with a 1000 µm thick CdTe sensor. RESULTS: The measurements for GOS scintillator thicknesses of 70 µm and 208 µm delivered 10% pre-sampled MTF values of 6.6 mm(-1) and 3.2 mm(-1), and DQE(0) values of 23% and 61%. The 10% pre-sampled MTF value for the 150 µm thick CsI scintillator 6.9 mm(-1), and the DQE(0) value was 49%. The CdTe sensor reached a 10% pre-sampled MTF value of 8.5 mm(-1) and a DQE(0) value of 85%. CONCLUSION: The photon-counting CdTe detector technology allows for significant dose reduction compared to the energy-integrating scintillation detector technology used in BCT today. Our comparative evaluation indicates that a high potential dose saving may be possible for BCT by using CdTe detectors, without loss of spatial resolution.

Generating and using patient-specific whole-body models for organ dose estimates in CT with increased accuracy: feasibility and validation.

The estimation of patient dose using Monte Carlo (MC) simulations based on the available patient CT images is limited to the length of the scan. Software tools for dose estimation based on standard computational phantoms overcome this problem; however, they are limited with respect to taking individual patient anatomy into account. The purpose of this study was to generate whole-body patient models in order to take scattered radiation and over-scanning effects into account. Thorax examinations were performed on three physical anthropomorphic phantoms at tube voltages of 80 kV and 120 kV; absorbed dose was measured using thermoluminescence dosimeters (TLD). Whole-body voxel models were built as a combination of the acquired CT images appended by data taken from widely used anthropomorphic voxel phantoms. MC simulations were performed both for the CT image volumes alone and for the whole-body models. Measured and calculated dose distributions were compared for each TLD chip position; additionally, organ doses were determined. MC simulations based only on CT data underestimated dose by 8%-15% on average depending on patient size with highest underestimation values of 37% for the adult phantom at the caudal border of the image volume. The use of whole-body models substantially reduced these errors; measured and simulated results consistently agreed to better than 10%. This study demonstrates that combined whole-body models can provide three-dimensional dose distributions with improved accuracy. Using the presented concept should be of high interest for research studies which demand high accuracy, e.g. for dose optimization efforts.

Experimental validation of a single shaped filter approach for CT using variable source-to-filter distance for examination of arbitrary object diameters.

The purpose of this study was to validate the use of a single shaped filter (SF) for computed tomography (CT) using variable source-to-filter distance (SFD) for the examination of different object diameters.A SF was designed by performing simulations with the purpose of achieving noise homogeneity in the reconstructed volume and dose reduction for arbitrary phantom diameters. This was accomplished by using a filter design method thats target is to achieve a homogeneous detector noise, but also uses a correction factor for the filtered back projection process. According to simulation results, a single SF designed for one of the largest phantom diameters meets the requirements for all diameters when SFD can be adjusted. To validate these results, a SF made of aluminium alloy was manufactured. Measurements were performed on a CT scanner with polymethyl methacrylate (PMMA) phantoms of diameters from 40-100 mm. The filter was positioned at SFDs ranging from 97-168 mm depending on the phantom diameter. Image quality was evaluated for the reconstructed volume by assessing CT value accuracy, noise homogeneity, contrast-to-noise ratio weighted by dose (CNRD) and spatial resolution. Furthermore, scatter distribution was determined with the use of a beam-stop phantom. Dose was measured for a PMMA phantom with a diameter of 100 mm using a calibrated ionization chamber.The application of a single SF at variable SFD led to improved noise uniformity and dose reduction: noise homogeneity was improved from 15% down to about 0%, and dose was reduced by about 37%. Furthermore, scatter dropped by about 32%, which led to reduced cupping artifacts and improved CT value accuracy. Spatial resolution and CNRD was not affected by the SF.By means of a single SF with variable SFD designed for CT, significant dose reduction can be achieved and image quality can be improved by reducing noise inhomogeneity as well as scatter-induced artifacts.

A quality assurance framework for the fully automated and objective evaluation of image quality in cone-beam computed tomography.

PURPOSE: Thousands of cone-beam computed tomography (CBCT) scanners for vascular, maxillofacial, neurological, and body imaging are in clinical use today, but there is no consensus on uniform acceptance and constancy testing for image quality (IQ) and dose yet. The authors developed a quality assurance (QA) framework for fully automated and time-efficient performance evaluation of these systems. In addition, the dependence of objective Fourier-based IQ metrics on direction and position in 3D volumes was investigated for CBCT. METHODS: The authors designed a dedicated QA phantom 10 cm in length consisting of five compartments, each with a diameter of 10 cm, and an optional extension ring 16 cm in diameter. A homogeneous section of water-equivalent material allows measuring CT value accuracy, image noise and uniformity, and multidimensional global and local noise power spectra (NPS). For the quantitative determination of 3D high-contrast spatial resolution, the modulation transfer function (MTF) of centrally and peripherally positioned aluminum spheres was computed from edge profiles. Additional in-plane and axial resolution patterns were used to assess resolution qualitatively. The characterization of low-contrast detectability as well as CT value linearity and artifact behavior was tested by utilizing sections with soft-tissue-equivalent and metallic inserts. For an automated QA procedure, a phantom detection algorithm was implemented. All tests used in the dedicated QA program were initially verified in simulation studies and experimentally confirmed on a clinical dental CBCT system. RESULTS: The automated IQ evaluation of volume data sets of the dental CBCT system was achieved with the proposed phantom requiring only one scan for the determination of all desired parameters. Typically, less than 5 min were needed for phantom set-up, scanning, and data analysis. Quantitative evaluation of system performance over time by comparison to previous examinations was also verified. The maximum percentage interscan variation of repeated measurements was less than 4% and 1.7% on average for all investigated quality criteria. The NPS-based image noise differed by less than 5% from the conventional standard deviation approach and spatially selective 10% MTF values were well comparable to subjective results obtained with 3D resolution pattern. Determining only transverse spatial resolution and global noise behavior in the central field of measurement turned out to be insufficient. CONCLUSIONS: The proposed framework transfers QA routines employed in conventional CT in an advanced version to CBCT for fully automated and time-efficient evaluation of technical equipment. With the modular phantom design, a routine as well as an expert version for assessing IQ is provided. The QA program can be used for arbitrary CT units to evaluate 3D imaging characteristics automatically.

Assessment of patient dose from CT localizer radiographs.

PURPOSE: With recently introduced technical innovations for CT systems, the dose of CT scan acquisitions has been substantially reduced; even effective dose values below 1 mSv have been reported. Due to this development, dose of the localizer radiograph may contribute substantially to dose of the whole CT examination. Since there are only limited data in the literature regarding patient dose for the different types of localizer radiographs, patient dose values were estimated in our study by measurements and Monte Carlo simulations and compared to dose values of typical CT examinations. METHODS: First, dose distributions were measured in anthropomorphic phantoms for three different body regions (head, thorax, abdomen-pelvic) and three positions of the x-ray tube (AP, PA, and lateral views); measured values were compared to simulated data using Monte Carlo techniques for validation purposes. Second, organ and effective dose values for the various investigated localizer radiograph scenarios were calculated and compared with published dose values for standard CT and low-dose CT examinations. RESULTS: For the anthropomorphic phantom, deviations of the dose values between measured and calculated results were in the range of 15%. Organ and effective dose values showed a strong dependence on the tube position. The largest differences were observed for chest localizer radiographs in the female phantom for the dose to the breast (AP: 1.01 mGy vs PA: 0.24 mGy). Overall effective dose values were in the range of 0.04-0.42 mSv per localizer radiograph acquisition. CONCLUSIONS: In view of the technical dose-reducing innovations in CT, localizer radiographs may substantially contribute to the total dose of the whole CT examination, particularly in the case of dedicated low-dose scans used, e.g., for young patients or screening purposes. Optimization of dose in localizer radiographs should be pursued further in the same way as it was done in CT.

Effect of shaped filter design on dose and image quality in breast CT.

The purpose of this study was to investigate the effect of shaped filters specifically designed for dedicated breast computed tomography (CT) scanners on dose and image quality. Optimization of filter shape and material in fan direction was performed using two different design methods, one aiming at homogeneous noise distributions in the CT images and the other aiming at a uniform dose distribution in the breast. The optimal filter thickness as a function of fan angle was determined iteratively to fulfil the above mentioned criteria for each breast diameter. Different filter materials (aluminium, copper, carbon, polytetrafluoroethylene) and breast phantoms with diameters between 80-180 mm were investigated. Noise uniformity in the reconstructed images, obtained from CT simulations based on ray-tracing methods, and dose in the breast, calculated with a Monte Carlo software tool, were used as figure of merit. Furthermore, CT-value homogeneity, the distribution of noise in cone direction, spatial resolution from centre to periphery and the contrast-to-noise ratio weighted by dose (CNRD) were evaluated. In addition, the decrease of scatter due to shaped filters was investigated. Since only few or one filter are practical in clinical CT systems, the effects of one shaped filter for different breast diameters were also investigated. In this case the filter, designed for the largest breast diameter, was simulated at variable source-to-filter distances depending on breast diameter. With the filter design method aiming at uniform noise distribution best results were obtained for aluminium as the filter material. Noise uniformity improved from 20} down to 5} and dose was reduced by about 30-40} for all breast diameters. No decrease of noise uniformity in cone direction, CT-value homogeneity, spatial resolution and the CNRD was detected with the shaped filter. However, a small improvement of CNRD was observed. Furthermore, a scatter reduction of about 20-30} and a more homogeneous scatter distribution were reached which led to reduced cupping artefacts. The simulations with one shaped filter at variable source-to-filter distance resulted in nearly homogeneous noise distributions and comparable dose reduction for all breast diameters. In conclusion, by means of shaped filters designed for breast CT, significant dose reduction can be achieved at unimpaired image quality. One shaped filter designed for the largest breast diameter used with variable source-to-filter distance appears to be the best solution for breast CT.

Fast on-site Monte Carlo tool for dose calculations in CT applications.

PURPOSE: Monte Carlo (MC) simulation is an established technique for dose calculation in diagnostic radiology. The major drawback is its high computational demand, which limits the possibility of usage in real-time applications. The aim of this study was to develop fast on-site computed tomography (CT) specific MC dose calculations by using a graphics processing unit (GPU) cluster. METHODS: GPUs are powerful systems which are especially suited to problems that can be expressed as data-parallel computations. In MC simulations, each photon track is independent of the others; each launched photon can be mapped to one thread on the GPU, thousands of threads are executed in parallel in order to achieve high performance. For further acceleration, the authors considered multiple GPUs. The total computation was divided into different parts which can be calculated in parallel on multiple devices. The GPU cluster is an MC calculation server which is connected to the CT scanner and computes 3D dose distributions on-site immediately after image reconstruction. To estimate the performance gain, the authors benchmarked dose calculation times on a 2.6 GHz Intel Xeon 5430 Quad core workstation equipped with two NVIDIA GeForce GTX 285 cards. The on-site calculation concept was demonstrated for clinical and preclinical datasets on CT scanners (multislice CT, flat-detector CT, and micro-CT) with varying geometry, spectra, and filtration. To validate the GPU-based MC algorithm, the authors measured dose values on a 64-slice CT system using calibrated ionization chambers and thermoluminesence dosimeters (TLDs) which were placed inside standard cylindrical polymethyl methacrylate (PMMA) phantoms. RESULTS: The dose values and profiles obtained by GPU-based MC simulations were in the expected good agreement with computed tomography dose index (CTDI) measurements and reference TLD profiles with differences being less than 5%. For 10(9) photon histories simulated in a 256 × 256 × 12 voxel thorax dataset with voxel size of 1.36 × 1.36 × 3.00 mm(3), calculation times of about 70 and 24 min were necessary with single-core and multiple-core central processing unit (CPU) solutions, respectively. Using GPUs, the same MC calculations were performed in 1.27 min (single card) and 0.65 min (two cards) without a loss in quality. Simulations were thus speeded up by factors up to 55 and 36 compared to single-core and multiple-core CPU, respectively. The performance scaled nearly linearly with the number of GPUs. Tests confirmed that the proposed GPU-based MC tool can be easily adapted to different types of CT scanners and used as service providers for fast on-site dose calculations. CONCLUSIONS: The Monte Carlo software package provides fast on-site calculation of 3D dose distributions in the CT suite which makes it a practical tool for any type of CT-specific application.

Development and evaluation of a phantom for dynamic contrast-enhanced imaging.

OBJECTIVES: Dynamic contrast-enhanced imaging allows assessing functional information in addition to morphology using various modalities. Several applications have been established in clinical practice; however, there is no standard with respect to injection protocols or postprocessing algorithms. The purpose of this study was to develop a phantom for generating reproducible contrast-enhancement curves and providing a standard for comparison of different protocols and modalities in dynamic imaging. MATERIALS AND METHODS: Our experimental setup consists of a peristaltic pump to generate a water flow through the phantom and a contrast injection pump. The phantom holds a sequence of layers allowing for assessment of perfusion, signal-to-noise ratio, and spatiotemporal resolution; the latter is the spatial resolution of structures with temporally changing contrast. Reproducibility was evaluated by the functional parameters time to peak, mean transit time, and peak enhancement by 24 scans over 4 weeks on a clinical computed tomography scanner. In addition, the area under the curve was evaluated for different injection durations at constant injection volume. Spatiotemporal resolution was assessed by spatial profiles on perfused bore patterns and compared for standard reconstructions, smooth reconstructions, and highly constrained backprojection for local reconstruction (HYPR LR). RESULTS: The phantom showed good reproducibility in repeated measurements, with maximal deviations of 4% for time to peak, 9% for mean transit time, and 8% for peak enhancement. Area under the curve was constant within 3.5% for different injection protocols. For the static case, HYPR LR maintained spatial resolution. For dynamic objects, however, HYPR LR reduced spatial resolution dependent on temporal dynamics by up to 19% for highest dynamics, which was still superior to smooth reconstructions (27%). CONCLUSIONS: The proposed phantom showed good reproducibility and therefore allows for comparing injection protocols or modalities in dynamic imaging. Assessment of spatiotemporal resolution under measurement conditions provides means for assessing postprocessing methods and reconstruction techniques in dynamic imaging.

Dosimetry concepts for scanner quality assurance and tissue dose assessment in micro-CT.

PURPOSE: At present, no established methods exist for dosimetry in micro computed tomography (micro-CT). The purpose of this study was therefore to investigate practical concepts for both dosimetric scanner quality assurance and tissue dose assessment for micro-CT. METHODS: The computed tomography dose index (CTDI) was adapted to micro-CT and measurements of the CTDI both free in air and in the center of cylindrical polymethyl methacrylate (PMMA) phantoms of 20 and 32 mm diameter were performed in a 6 month interval with a 100 mm pencil ionization chamber calibrated for low tube voltages. For tissue dose assessment, z-profile measurements using thermoluminescence dosimeters (TLDs) were performed and both profile and CTDI measurements were compared to Monte Carlo (MC) dose calculations to validate an existing MC tool for use in micro-CT. The consistency of MC calculations and TLD measurements was further investigated in two mice cadavers. RESULTS: CTDI was found to be a reproducible quantity for constancy tests on the micro-CT system under study, showing a linear dependence on tube voltage and being by definition proportional to mAs setting and z-collimation. The CTDI measured free in air showed larger systematic deviations after the 6 month interval compared to the CTDI measured in PMMA phantoms. MC calculations were found to match CTDI measurements within 3% when using x-ray spectra measured at our micro-CT installation and better than 10% when using x-ray spectra calculated from semi-empirical models. Visual inspection revealed good agreement for all z-profiles. The consistency of MC calculations and TLD measurements in mice was found to be better than 10% with a mean deviation of 4.5%. CONCLUSIONS: Our results show the CTDI implemented for micro-CT to be a promising candidate for dosimetric quality assurance measurements as it linearly reflects changes in tube voltage, mAs setting, and collimation used during the scan, encouraging further studies on a variety of systems. For tissue dose assessment, MC calculations offer an accurate and fast alternative to TLD measurements allowing for dose calculations specific to any geometry and scan protocol.

Iterative reconstruction methods in X-ray CT.

Iterative reconstruction (IR) methods have recently re-emerged in transmission x-ray computed tomography (CT). They were successfully used in the early years of CT, but given up when the amount of measured data increased because of the higher computational demands of IR compared to analytical methods. The availability of large computational capacities in normal workstations and the ongoing efforts towards lower doses in CT have changed the situation; IR has become a hot topic for all major vendors of clinical CT systems in the past 5 years. This review strives to provide information on IR methods and aims at interested physicists and physicians already active in the field of CT. We give an overview on the terminology used and an introduction to the most important algorithmic concepts including references for further reading. As a practical example, details on a model-based iterative reconstruction algorithm implemented on a modern graphics adapter (GPU) are presented, followed by application examples for several dedicated CT scanners in order to demonstrate the performance and potential of iterative reconstruction methods. Finally, some general thoughts regarding the advantages and disadvantages of IR methods as well as open points for research in this field are discussed.

A formulation of tissue- and water-equivalent materials using the stoichiometric analysis method for CT-number calibration in radiotherapy treatment planning.

Tissue- and water-equivalent materials (TEMs) are widely used in quality assurance and calibration procedures, both in radiodiagnostics and radiotherapy. In radiotherapy, particularly, the TEMs are often used for computed tomography (CT) number calibration in treatment planning systems. However, currently available TEMs may not be very accurate in the determination of the calibration curves due to their limitation in mimicking radiation characteristics of the corresponding real tissues in both low- and high-energy ranges. Therefore, we are proposing a new formulation of TEMs using a stoichiometric analysis method to obtain TEMs for the calibration purposes. We combined the stoichiometric calibration and the basic data method to compose base materials to develop TEMs matching standard real tissues from ICRU Report 44 and 46. First, the CT numbers of six materials with known elemental compositions were measured to get constants for the stoichiometric calibration. The results of the stoichiometric calibration were used together with the basic data method to formulate new TEMs. These new TEMs were scanned to validate their CT numbers. The electron density and the stopping power calibration curves were also generated. The absolute differences of the measured CT numbers of the new TEMs were less than 4 HU for the soft tissues and less than 22 HU for the bone compared to the ICRU real tissues. Furthermore, the calculated relative electron density and electron and proton stopping powers of the new TEMs differed by less than 2% from the corresponding ICRU real tissues. The new TEMs which were formulated using the proposed technique increase the simplicity of the calibration process and preserve the accuracy of the stoichiometric calibration simultaneously.

High-resolution spiral CT of the breast at very low dose: concept and feasibility considerations.

OBJECTIVE: Mammography, today's standard imaging approach, has deficits with respect to the superimposition of anatomical structures. Dedicated CT of the breast so far indicated that it can provide superior soft-tissue imaging, but that it still has significant limitations with respect to spatial resolution and dose. We have assessed novel dedicated breast CT technology. METHODS: Based on simulations and measurements we developed novel technology which uses direct-conversion CdTe material and photon-counting electronics with 100 µm detector element size for close to 100% dose efficiency. We assessed the potential for the imaging of microcalcifications of 100 to 200 µm diameter and soft-tissue lesions of 1 to 5 mm diameter by simulations at dose levels between 1 and 6 mGy. RESULTS: Microcalcifications of 150 µm and soft-tissue lesions of 2 mm diameter were found to be clearly detectable at an average glandular dose of 3 mGy. Separate displays are required for high-resolution microcalcification and for low-resolution soft-tissue analysis. Total CT data acquisition time will be below 10 s. CONCLUSION: Dedicated breast CT may eventually provide comprehensive diagnostic assessment of microcalcifications and soft-tissue structures at dose levels equivalent to or below those of two-view screening mammography.

Semiempirical analysis of materials' elemental composition to formulate tissue-equivalent materials: a preliminary study.

Tissue-equivalent materials are used for simplifying quality control and quality assurance procedures, both in diagnostic and therapeutic radiology. Important information to formulate a tissue-equivalent material is elemental composition of its base materials. However, this information is not easily obtained. Therefore we propose a stoichiometric analysis method to investigate the elemental composition of the base materials that can potentially be used for manufacturing tissue-equivalent materials. In this technique, we combined the stoichiometric calibration and the basic data method to obtain the elemental composition of materials from measured computer tomography (CT) numbers. The elemental composition, with the maximum number of the elements of the material in question up to the available number of different tube voltages at the CT scanner, was analysed using the proposed approach. We tested eight different cylinders in this study. The estimated elemental compositions of unspecified materials in the cylinders were evaluated by comparing the calculated and the simulated CT numbers to the measured ones; the results showed good correlation with maximum absolute differences of 1.9 and 3.7 HU, respectively. The accuracy of the stoichiometric analysis method to estimate the elemental composition was influenced by the accuracy of the measured CT numbers. The method proposed allows for determining the elemental composition of the base materials which can then be applied further to formulate tissue-equivalent materials.