IMRT credentialing for prospective trials using institutional virtual phantoms: results of a joint European Organization for the Research and Treatment of Cancer and Radiological Physics Center project.

BACKGROUND AND PURPOSE Intensity-modulated radiotherapy (IMRT) credentialing for a EORTC study was performed using an anthropomorphic head phantom from the Radiological Physics Center (RPC; RPC(PH)). Institutions were retrospectively requested to irradiate their institutional phantom (INST(PH)) using the same treatment plan in the framework of a Virtual Phantom Project (VPP) for IMRT credentialing. MATERIALS AND METHODS CT data set of the institutional phantom and measured 2D dose matrices were requested from centers and sent to a dedicated secure EORTC uploader. Data from the RPC(PH) and INST(PH) were thereafter centrally analyzed and inter-compared by the QA team using commercially available software (RIT; ver.5.2; Colorado Springs, USA). RESULTS Eighteen institutions participated to the VPP. The measurements of 6 (33%) institutions could not be analyzed centrally. All other centers passed both the VPP and the RPC ±7%/4 mm credentialing criteria. At the 5%/5 mm gamma criteria (90% of pixels passing), 11(92%) as compared to 12 (100%) centers pass the credentialing process with RPC(PH) and INST(PH) (p = 0.29), respectively. The corresponding pass rate for the 3%/3 mm gamma criteria (90% of pixels passing) was 2 (17%) and 9 (75%; p = 0.01), respectively. CONCLUSIONS IMRT dosimetry gamma evaluations in a single plane for a H&N prospective trial using the INST(PH) measurements showed agreement at the gamma index criteria of ±5%/5 mm (90% of pixels passing) for a small number of VPP measurements. Using more stringent, criteria, the RPC(PH) and INST(PH) comparison showed disagreement. More data is warranted and urgently required within the framework of prospective studies.

Monte Carlo Dose Calculation on Deforming Anatomy

This article presents the implementation and validation of a dose calculation approach for deforming anatomical objects. Deformation is represented by deformation vector fields leading to deformed voxel grids representing the different deformation scenarios. Particle transport in the resulting deformed voxels is handled through the approximation of voxel surfaces by triangles in the geometry implementation of the Swiss Monte Carlo Plan framework. The focus lies on the validation methodology which uses computational phantoms representing the same physical object through regular and irregular voxel grids. These phantoms are chosen such that the new implementation for a deformed voxel grid can be compared directly with an established dose calculation algorithm for regular grids. Furthermore, separate validation of the aspects voxel geometry and the density changes resulting from deformation is achieved through suitable design of the validation phantom. We show that equivalent results are obtained with the proposed method and that no statistically significant errors are introduced through the implementation for irregular voxel geometries. This enables the use of the presented and validated implementation for further investigations of dose calculation on deforming anatomy.

A 20-channel receive-only mouse array coil for a 3 T clinical MRI system

A 20-channel phased-array coil for MRI of mice has been designed, constructed, and validated with bench measurements and high-resolution accelerated imaging. The technical challenges of designing a small, high density array have been overcome using individual small-diameter coil elements arranged on a cylinder in a hexagonal overlapping design with adjacent low impedance preamplifiers to further decouple the array elements. Signal-to-noise ratio (SNR) and noise amplification in accelerated imaging were simulated and quantitatively evaluated in phantoms and in vivo mouse images. Comparison between the 20-channel mouse array and a length-matched quadrature driven small animal birdcage coil showed an SNR increase at the periphery and in the center of the phantom of 3- and 1.3-fold, respectively. Comparison with a shorter but SNR-optimized birdcage coil (aspect ratio 1:1 and only half mouse coverage) showed an SNR gain of twofold at the edge of the phantom and similar SNR in the center. G-factor measurements indicate that the coil is well suited to acquire highly accelerated images.

High-accuracy patient-to-image registration for the facilitation of image-guided robotic microsurgery on the head

Image-guided microsurgery requires accuracies an order of magnitude higher than today's navigation systems provide. A critical step toward the achievement of such low-error requirements is a highly accurate and verified patient-to-image registration. With the aim of reducing target registration error to a level that would facilitate the use of image-guided robotic microsurgery on the rigid anatomy of the head, we have developed a semiautomatic fiducial detection technique. Automatic force-controlled localization of fiducials on the patient is achieved through the implementation of a robotic-controlled tactile search within the head of a standard surgical screw. Precise detection of the corresponding fiducials in the image data is realized using an automated model-based matching algorithm on high-resolution, isometric cone beam CT images. Verification of the registration technique on phantoms demonstrated that through the elimination of user variability, clinically relevant target registration errors of approximately 0.1 mm could be achieved.

Generalized eMC implementation for Monte Carlo dose calculation of electron beams from different machine types

The electron Monte Carlo (eMC) dose calculation algorithm available in the Eclipse treatment planning system (Varian Medical Systems) is based on the macro MC method and uses a beam model applicable to Varian linear accelerators. This leads to limitations in accuracy if eMC is applied to non-Varian machines. In this work eMC is generalized to also allow accurate dose calculations for electron beams from Elekta and Siemens accelerators. First, changes made in the previous study to use eMC for low electron beam energies of Varian accelerators are applied. Then, a generalized beam model is developed using a main electron source and a main photon source representing electrons and photons from the scattering foil, respectively, an edge source of electrons, a transmission source of photons and a line source of electrons and photons representing the particles from the scrapers or inserts and head scatter radiation. Regarding the macro MC dose calculation algorithm, the transport code of the secondary particles is improved. The macro MC dose calculations are validated with corresponding dose calculations using EGSnrc in homogeneous and inhomogeneous phantoms. The validation of the generalized eMC is carried out by comparing calculated and measured dose distributions in water for Varian, Elekta and Siemens machines for a variety of beam energies, applicator sizes and SSDs. The comparisons are performed in units of cGy per MU. Overall, a general agreement between calculated and measured dose distributions for all machine types and all combinations of parameters investigated is found to be within 2% or 2 mm. The results of the dose comparisons suggest that the generalized eMC is now suitable to calculate dose distributions for Varian, Elekta and Siemens linear accelerators with sufficient accuracy in the range of the investigated combinations of beam energies, applicator sizes and SSDs.

Macro Monte Carlo for dose calculation of proton beams

Although the Monte Carlo (MC) method allows accurate dose calculation for proton radiotherapy, its usage is limited due to long computing time. In order to gain efficiency, a new macro MC (MMC) technique for proton dose calculations has been developed. The basic principle of the MMC transport is a local to global MC approach. The local simulations using GEANT4 consist of mono-energetic proton pencil beams impinging perpendicularly on slabs of different thicknesses and different materials (water, air, lung, adipose, muscle, spongiosa, cortical bone). During the local simulation multiple scattering, ionization as well as elastic and inelastic interactions have been taken into account and the physical characteristics such as lateral displacement, direction distributions and energy loss have been scored for primary and secondary particles. The scored data from appropriate slabs is then used for the stepwise transport of the protons in the MMC simulation while calculating the energy loss along the path between entrance and exit position. Additionally, based on local simulations the radiation transport of neutrons and the generated ions are included into the MMC simulations for the dose calculations. In order to validate the MMC transport, calculated dose distributions using the MMC transport and GEANT4 have been compared for different mono-energetic proton pencil beams impinging on different phantoms including homogeneous and inhomogeneous situations as well as on a patient CT scan. The agreement of calculated integral depth dose curves is better than 1% or 1 mm for all pencil beams and phantoms considered. For the dose profiles the agreement is within 1% or 1 mm in all phantoms for all energies and depths. The comparison of the dose distribution calculated using either GEANT4 or MMC in the patient also shows an agreement of within 1% or 1 mm. The efficiency of MMC is up to 200 times higher than for GEANT4. The very good level of agreement in the dose comparisons demonstrate that the newly developed MMC transport results in very accurate and efficient dose calculations for proton beams.

Determining the optical properties of a gelatin-TiO2 phantom at 780 nm

Tissue phantoms play a central role in validating biomedical imaging techniques. Here we employ a series of methods that aim to fully determine the optical properties, i.e., the refractive index n, absorption coefficient μa, transport mean free path ℓ∗, and scattering coefficient μs of a TiO2 in gelatin phantom intended for use in optoacoustic imaging. For the determination of the key parameters μa and ℓ∗, we employ a variant of time of flight measurements, where fiber optodes are immersed into the phantom to minimize the influence of boundaries. The robustness of the method was verified with Monte Carlo simulations, where the experimentally obtained values served as input parameters for the simulations. The excellent agreement between simulations and experiments confirmed the reliability of the results. The parameters determined at 780 nm are n=1.359(±0.002), μ′s=1/ℓ∗=0.22(±0.02) mm-1, μa= 0.0053(+0.0006-0.0003) mm-1, and μs=2.86(±0.04) mm-1. The asymmetry parameter g obtained from the parameters ℓ∗ and μ′s is 0.93, which indicates that the scattering entities are not bare TiO2 particles but large sparse clusters. The interaction between the scattering particles and the gelatin matrix should be taken into account when developing such phantoms.

Effect of multislice acquisition on T1 and T2 measurements of articular cartilage at 3T

PURPOSE: The aim of this study was to investigate the effect of magnetization transfer on multislice T(1) and T(2) measurements of articular cartilage. MATERIALS AND METHODS: A set of phantoms with different concentrations of collagen and contrast agent (Gd-DTPA(2-)) were used for the in vitro study. A total of 20 healthy knees were used for the in vivo study. T(1) and T(2) measurements were performed using fast-spin-echo inversion-recovery (FSE-IR) sequence and multi-spin-echo (MSE) sequence, respectively, in both in vitro and in vivo studies. We investigated the difference in T(1) and T(2) values between that measured by single-slice acquisition and that measured by multislice acquisition. RESULTS: Regarding T(1) measurement, a large drop of T(1) in all slices and also a large interslice variation in T(1) were observed when multislice acquisition was used. Regarding T(2) measurement, a substantial drop of T(2) in all slices was observed; however, there was no apparent interslice variation when multislice acquisition was used. CONCLUSION: This study demonstrated that the adaptation of multislice acquisition technique for T(1) measurement using FSE-IR methodology is difficult and its use for clinical evaluation is problematic. In contrast, multislice acquisition for T(2) measurement using MSE was clinically applicable if inaccuracies caused by multislice acquisition were taken into account. J. Magn. Reson. Imaging 2007;26:109-117. (c) 2007 Wiley-Liss, Inc.

Conversion of CT numbers into tissue parameters for Monte Carlo dose calculations: a multi-centre study

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

Comparative dose measurements by spiral tomography for preimplant diagnosis: the Scanora machine versus the Cranex Tome radiography unit

Objective. The purpose of this study was to determine the dose profile of the Cranex Tome radiography unit and compare it with that of the Scanora machine.Study design. The radiation dose delivered by the Cranex Tome radiography unit during the cross-sectional mode was determined. Single tooth gaps in regions 3 (16) and 30 (46) were simulated. Dosimetry was carried out with 2 phantoms, a head and neck phantom and a full-body phantom loaded with 142 thermoluminescent dosimeters (TLD) and 280 TLD, respectively; all locations corresponded to radiosensitive organs or tissues. The recorded local mean organ doses were compared with those measured in another study evaluating the Scanora machine.Results. Generally, dose values from the Cranex Tome radiography unit reached only 50% to 60% of the values measured for the Scanora machine. The effective dose was calculated as 0.061 mSv and 0.04 mSv for tooth regions 3 (16) and 30 (46), respectively. Corresponding values for the Scanora machine were 0.117 mSv and 0.084 mSv.Conclusion. Cross-sectional imaging in the molar region of the upper and the lower jaw can be performed with the Cranex Tome unit, which delivers only approximately half of the dose that the Scanora machine delivers.

Hypothetical mortality risk associated with spiral computed tomography of the maxilla and mandible

In the present study, dose measurements have been conducted following examination of the maxilla and mandible with spiral computed tomography (CT). The measurements were carried out with 2 phantoms, a head and neck phantom and a full body phantom. The analysis of applied thermoluminescent dosimeters yielded radiation doses for organs and tissues in the head and neck region between 0.6 and 16.7 mGy when 40 axial slices and 120 kV/165 mAs were used as exposure parameters. The effective dose was calculated as 0.58 and 0.48 mSv in the maxilla and mandible, respectively. Tested methods for dose reduction showed a significant decrease of radiation dose from 40 to 65%. Based on these results, the mortality risk was estimated according to calculation models recommended by the Committee on the Biological Effects of Ionizing Radiations and by the International Commission on Radiological Protection. Both models resulted in similar values. The mortality risk ranges from 46.2 x 10.6 for 20-year-old men to 11.2 x 10(-6) for 65-year-old women. Using 2 methods of dose reduction, the mortality risk decreased by approximately 50 to 60% to 19.1 x 10(-6) for 20-year-old men and 5.5 x 10(-6) for 65-year-old women. It can be concluded that a CT scan of the maxillofacial complex causes a considerable radiation dose when compared with conventional radiographic examinations. Therefore, a careful indication for this imaging technique and dose reduction methods should be considered in daily practice.

Effect of X-ray tube parameters, iodine concentration, and patient size on image quality in pulmonary computed tomography angiography: a chest-phantom-study

OBJECTIVES: The aim of this phantom study was to evaluate the contrast-to-noise ratio (CNR) in pulmonary computed tomography (CT)-angiography for 300 and 400 mg iodine/mL contrast media using variable x-ray tube parameters and patient sizes. We also analyzed the possible strategies of dose reduction in patients with different sizes. MATERIALS AND METHODS: The segmental pulmonary arteries were simulated by plastic tubes filled with 1:30 diluted solutions of 300 and 400 mg iodine/mL contrast media in a chest phantom mimicking thick, intermediate, and thin patients. Volume scanning was done with a CT scanner at 80, 100, 120, and 140 kVp. Tube current-time products (mAs) varied between 50 and 120% of the optimal value given by the built-in automatic dose optimization protocol. Attenuation values and CNR for both contrast media were evaluated and compared with the volume CT dose index (CTDI(vol)). Figure of merit, calculated as CNR/CTDIvol, was used to quantify image quality improvement per exposure risk to the patient. RESULTS: Attenuation of iodinated contrast media increased both with decreasing tube voltage and patient size. A CTDIvol reduction by 44% was achieved in the thin phantom with the use of 80 instead of 140 kVp without deterioration of CNR. Figure of merit correlated with kVp in the thin phantom (r = -0.897 to -0.999; P < 0.05) but not in the intermediate and thick phantoms (P = 0.09-0.71), reflecting a decreasing benefit of tube voltage reduction on image quality as the thickness of the phantom increased. Compared with the 300 mg iodine/mL concentration, the same CNR for 400 mg iodine/mL contrast medium was achieved at a lower CTDIvol by 18 to 40%, depending on phantom size and applied tube voltage. CONCLUSIONS: Low kVp protocols for pulmonary embolism are potentially advantageous especially in thin and, to a lesser extent, in intermediate patients. Thin patients profit from low voltage protocols preserving a good CNR at a lower exposure. The use of 80 kVp in obese patients may be problematic because of the limitation of the tube current available, reduced CNR, and high skin dose. The high CNR of the 400 mg iodine/mL contrast medium together with lower tube energy and/or current can be used for exposure reduction.

Rapid estimation of cartilage T2 based on double echo at steady state (DESS) with 3 Tesla

The double-echo-steady-state (DESS) sequence generates two signal echoes that are characterized by a different contrast behavior. Based on these two contrasts, the underlying T2 can be calculated. For a flip-angle of 90 degrees , the calculated T2 becomes independent of T1, but with very low signal-to-noise ratio. In the present study, the estimation of cartilage T2, based on DESS with a reduced flip-angle, was investigated, with the goal of optimizing SNR, and simultaneously minimizing the error in T2. This approach was validated in phantoms and on volunteers. T2 estimations based on DESS at different flip-angles were compared with standard multiecho, spin-echo T2. Furthermore, DESS-T2 estimations were used in a volunteer and in an initial study on patients after cartilage repair of the knee. A flip-angle of 33 degrees was the best compromise for the combination of DESS-T2 mapping and morphological imaging. For this flip angle, the Pearson correlation was 0.993 in the phantom study (approximately 20% relative difference between SE-T2 and DESS-T2); and varied between 0.429 and 0.514 in the volunteer study. Measurements in patients showed comparable results for both techniques with regard to zonal assessment. This DESS-T2 approach represents an opportunity to combine morphological and quantitative cartilage MRI in a rapid one-step examination.

Quality controls for two heel bone ultrasounds used in the Swiss Evaluation of the Methods of Measurement of Osteoporotic Fracture Risk Study

Because of the important morbidity and mortality associated with osteoporosis, it is essential to detect subjects at risk by screening methods, such as bone quantitative ultrasounds (QUSs). Several studies showed that QUS could predict fractures. None, however, compared prospectively different QUS devices, and few data of quality controls (QCs) have been published. The Swiss Evaluation of the Methods of Measurement of Osteoporotic Fracture Risk is a prospective multicenter study that compared three QUSs for the assessment of hip fracture risk in a population of 7609 women age >/=70 yr. Because the inclusion phase lasted 20 mo, and because 10 centers participated in this study, QC became a major issue. We therefore developed a QC procedure to assess the stability and precision of the devices, and for their cross-calibration. Our study focuses on the two heel QUSs. The water bath system (Achilles+) had a higher precision than the dry system (Sahara). The QC results were highly dependent on temperature. QUS stability was acceptable, but Sahara must be calibrated regularly. A sufficient homogeneity among all the Sahara devices could be demonstrated, whereas significant differences were found among the Achilles+ devices. For speed of sound, 52% of the differences among the Achilles+ was explained by the water s temperature. However, for broadband ultrasound attenuation, a maximal difference of 23% persisted after adjustment for temperature. Because such differences could influence measurements in vivo, it is crucial to develop standardized phantoms to be used in prospective multicenter studies.

Lung cancer screening with CT: evaluation of radiologists and different computer assisted detection software (CAD) as first and second readers for lung nodule detection at different dose levels

OBJECTIVES To find the best pairing of first and second reader at highest sensitivity for detecting lung nodules with CT at various dose levels. MATERIALS AND METHODS An anthropomorphic lung phantom and artificial lung nodules were used to simulate screening CT-examination at standard dose (100 mAs, 120 kVp) and 8 different low dose levels, using 120, 100 and 80 kVp combined with 100, 50 and 25 mAs. At each dose level 40 phantoms were randomly filled with 75 solid and 25 ground glass nodules (5-12 mm). Two radiologists and 3 different computer aided detection softwares (CAD) were paired to find the highest sensitivity. RESULTS Sensitivities at standard dose were 92%, 90%, 84%, 79% and 73% for reader 1, 2, CAD1, CAD2, CAD3, respectively. Combined sensitivity for human readers 1 and 2 improved to 97%, (p1=0.063, p2=0.016). Highest sensitivities--between 97% and 99.0%--were achieved by combining any radiologist with any CAD at any dose level. Combining any two CADs, sensitivities between 85% and 88% were significantly lower than for radiologists combined with CAD (p<0.03). CONCLUSIONS Combination of a human observer with any of the tested CAD systems provide optimal sensitivity for lung nodule detection even at reduced dose at 25 mAs/80 kVp.