Moti Paudel

Evaluation of a commercial MRI Linac based Monte Carlo dose calculation algorithm with geant 4

PURPOSE This paper provides a comparison between a fast, commercial, in-patient Monte Carlo dose calculation algorithm (GPUMCD) and geant4. It also evaluates the dosimetric impact of the application of an external 1.5 T magnetic field. METHODS A stand-alone version of the Elekta™ GPUMCD algorithm, to be used within the Monaco treatment planning system to model dose for the Elekta™ magnetic resonance imaging (MRI) Linac, was compared against GEANT4 (v10.1). This was done in the presence or absence of a 1.5 T static magnetic field directed orthogonally to the radiation beam axis. Phantoms with material compositions of water, ICRU lung, ICRU compact-bone, and titanium were used for this purpose. Beams with 2 MeV monoenergetic photons as well as a 7 MV histogrammed spectrum representing the MRI Linac spectrum were emitted from a point source using a nominal source-to-surface distance of 142.5 cm. Field sizes ranged from 1.5 × 1.5 to 10 × 10 cm(2). Dose scoring was performed using a 3D grid comprising 1 mm(3) voxels. The production thresholds were equivalent for both codes. Results were analyzed based upon a voxel by voxel dose difference between the two codes and also using a volumetric gamma analysis. RESULTS Comparisons were drawn from central axis depth doses, cross beam profiles, and isodose contours. Both in the presence and absence of a 1.5 T static magnetic field the relative differences in doses scored along the beam central axis were less than 1% for the homogeneous water phantom and all results matched within a maximum of ±2% for heterogeneous phantoms. Volumetric gamma analysis indicated that more than 99% of the examined volume passed gamma criteria of 2%-2 mm (dose difference and distance to agreement, respectively). These criteria were chosen because the minimum primary statistical uncertainty in dose scoring voxels was 0.5%. The presence of the magnetic field affects the dose at the interface depending upon the density of the material on either sides of the interface. This effect varies with the field size. For example, at the water-lung interface a 33.94% increase in dose was observed (relative to the Dmax), by both GPUMCD and GEANT4 for the field size of 2 × 2 cm(2) (compared to no B-field case), which increased to 47.83% for the field size of 5 × 5 cm(2) in the presence of the magnetic field. Similarly, at the lung-water interface, the dose decreased by 19.21% (relative to Dmax) for a field size of 2 × 2 cm(2) and by 30.01% for 5 × 5 cm(2) field size. For more complex combinations of materials the dose deposition also becomes more complex. CONCLUSIONS The GPUMCD algorithm showed good agreement against GEANT4 both in the presence and absence of a 1.5 T external magnetic field. The application of 1.5 T magnetic field significantly alters the dose at the interfaces by either increasing or decreasing the dose depending upon the density of the material on either side of the interfaces.

Experimental evaluation of a GPU-based Monte Carlo dose calculation algorithm in the Monaco treatment planning system

A new GPU-based Monte Carlo dose calculation algorithm (GPUMCD), developed by the vendor Elekta for the Monaco treatment planning system (TPS), is capable of modeling dose for both a standard linear accelerator and an Elekta MRI linear accelerator. We have experimentally evaluated this algorithm for a standard Elekta Agility linear accelerator. A beam model was developed in the Monaco TPS (research version 5.09.06) using the commissioned beam data for a 6 MV Agility linac. A heterogeneous phantom representing several scenarios - tumor-in-lung, lung, and bone-in-tissue - was designed and built. Dose calculations in Monaco were done using both the current clinical Monte Carlo algorithm, XVMC, and the new GPUMCD algorithm. Dose calculations in a Pinnacle TPS were also produced using the collapsed cone convolution (CCC) algorithm with heterogeneity correction. Calculations were compared with the measured doses using an ionization chamber (A1SL) and Gafchromic EBT3 films for 2×2 cm2,5×5 cm2, and 10×2 cm2 field sizes. The percentage depth doses (PDDs) calculated by XVMC and GPUMCD in a homogeneous solid water phantom were within 2%/2 mm of film measurements and within 1% of ion chamber measurements. For the tumor-in-lung phantom, the calculated doses were within 2.5%/2.5 mm of film measurements for GPUMCD. For the lung phantom, doses calculated by all of the algorithms were within 3%/3 mm of film measurements, except for the 2×2 cm2 field size where the CCC algorithm underestimated the depth dose by ∼5% in a larger extent of the lung region. For the bone phantom, all of the algorithms were equivalent and calculated dose to within 2%/2 mm of film measurements, except at the interfaces. Both GPUMCD and XVMC showed interface effects, which were more pronounced for GPUMCD and were comparable to film measurements, whereas the CCC algorithm showed these effects poorly. PACS number(s): 87.53.Bn, 87.55.dh, 87.55.km.A new GPU‐based Monte Carlo dose calculation algorithm (GPUMCD), developed by the vendor Elekta for the Monaco treatment planning system (TPS), is capable of modeling dose for both a standard linear accelerator and an Elekta MRI linear accelerator. We have experimentally evaluated this algorithm for a standard Elekta Agility linear accelerator. A beam model was developed in the Monaco TPS (research version 5.09.06) using the commissioned beam data for a 6 MV Agility linac. A heterogeneous phantom representing several scenarios — tumor‐in‐lung, lung, and bone‐in‐tissue — was designed and built. Dose calculations in Monaco were done using both the current clinical Monte Carlo algorithm, XVMC, and the new GPUMCD algorithm. Dose calculations in a Pinnacle TPS were also produced using the collapsed cone convolution (CCC) algorithm with heterogeneity correction. Calculations were compared with the measured doses using an ionization chamber (A1SL) and Gafchromic EBT3 films for 2×2 cm2,5×5 cm2, and 10×2 cm2 field sizes. The percentage depth doses (PDDs) calculated by XVMC and GPUMCD in a homogeneous solid water phantom were within 2%/2 mm of film measurements and within 1% of ion chamber measurements. For the tumor‐in‐lung phantom, the calculated doses were within 2.5%/2.5 mm of film measurements for GPUMCD. For the lung phantom, doses calculated by all of the algorithms were within 3%/3 mm of film measurements, except for the 2×2 cm2 field size where the CCC algorithm underestimated the depth dose by ∼5% in a larger extent of the lung region. For the bone phantom, all of the algorithms were equivalent and calculated dose to within 2%/2 mm of film measurements, except at the interfaces. Both GPUMCD and XVMC showed interface effects, which were more pronounced for GPUMCD and were comparable to film measurements, whereas the CCC algorithm showed these effects poorly. PACS number(s): 87.53.Bn, 87.55.dh, 87.55.km

SU-E-T-416: Experimental Evaluation of a Commercial GPU-Based Monte Carlo Dose Calculation Algorithm

Purpose: A new commercial GPU-based Monte Carlo dose calculation algorithm (GPUMCD) developed by the vendor Elekta™ to be used in the Monaco Treatment Planning System (TPS) is capable of modeling dose for both a standard linear accelerator and for an Elekta MRI-Linear accelerator (modeling magnetic field effects). We are evaluating this algorithm in two parts: commissioning the algorithm for an Elekta Agility linear accelerator (the focus of this work) and evaluating the algorithm’s ability to model magnetic field effects for an MRI-linear accelerator. Methods: A beam model was developed in the Monaco TPS (v.5.09.06) using the commissioned beam data for a 6MV Agility linac. A heterogeneous phantom representing tumor-in-lung, lung, bone-in-tissue, and prosthetic was designed/built. Dose calculations in Monaco were done using the current clinical algorithm (XVMC) and the new GPUMCD algorithm (1 mm3 voxel size, 0.5% statistical uncertainty) and in the Pinnacle TPS using the collapsed cone convolution (CCC) algorithm. These were compared with the measured doses using an ionization chamber (A1SL) and Gafchromic EBT3 films for 2×2 cm2, 5×5 cm2, and 10×10 cm2 field sizes. Results: The calculated central axis percentage depth doses (PDDs) in homogeneous solid water were within 2% compared to measurements for XVMC and GPUMCD. For tumor-in-lung and lung phantoms, doses calculated by all of the algorithms were within the experimental uncertainty of the measurements (±2% in the homogeneous phantom and ±3% for the tumor-in-lung or lung phantoms), except for 2×2 cm2 field size where only the CCC algorithm differs from film by 5% in the lung region. The analysis for bone-in-tissue and the prosthetic phantoms are ongoing. Conclusion: The new GPUMCD algorithm calculated dose comparable to both the XVMC algorithm and to measurements in both a homogeneous solid water medium and the heterogeneous phantom representing lung or tumor-in-lung for 2×2 cm2-10×10 cm2 field sizes. Funding support was obtained from Elekta.

SU‐E‐T‐203: Comparison of a Commercial MRI‐Linear Accelerator Based Monte Carlo Dose Calculation Algorithm and Geant4

Purpose: An MRI-linear accelerator is currently being developed by the vendor Elekta™. The treatment planning system that will be used to model dose for this unit uses a Monte Carlo dose calculation algorithm, GPUMCD, that allows for the application of a magnetic field. We tested this radiation transport code against an independent Monte-Carlo toolkit Geant4 (v.4.10.01) both with and without the magnetic field applied. Methods: The setup comprised a 6 MeV mono-energetic photon beam emerging from a point source impinging on a homogeneous water phantom at 100 cm SSD. The comparisons were drawn from the percentage depth doses (PDD) for three different field sizes (1.5 x 1.5 cm2, 5 x 5 cm2, 10 x 10 cm2) and dose profiles at various depths. A 1.5 T magnetic field was applied perpendicular to the direction of the beam. The transport thresholds were kept the same for both codes. Results: All of the normalized PDDs and profiles agreed within ± 1 %. In the presence of the magnetic field, PDDs rise more quickly reducing the depth of maximum dose. Near the beam exit point in the phantom a hot spot is created due to the electron return effect. This effect is more pronounced for the larger field sizes. Profiles selected parallel to the external field show no effect, however, the ones selected perpendicular to the direction of the applied magnetic field are shifted towards the direction of the Lorentz force applied by the magnetic field on the secondary electrons. It is observed that these profiles are not symmetric which indicates a lateral build up of the dose. Conclusion: There is a good general agreement between the PDDs/profiles calculated by both algorithms thus far. We are proceeding towards clinically relevant comparisons in a heterogeneous phantom for polyenergetic beams. Funding for this work has been provided by Elekta.

The dosimetric impact of gadolinium-based contrast media in GBM brain patient plans for a MRI-Linac

Dosimetric effects of gadolinium based contrast media (Gadovist) were evaluated for the Elekta MRI linear accelerator using the research version of the Monaco treatment planning system (TPS). In order to represent a gadolinium uptake, the contrast was manually assigned to a phantom as well as to the gross tumour volume (GTV) of 6 glioblastoma multiforme (GBM) patients. A preliminary estimate of the dose enhancement, due to gadolinium, was performed using the phantom irradiated with a single beam. A more complicated assessment was performed for the GBM patients using a 7 field IMRT technique. The material table in Monaco was modified in order to identify the presence of a non-biological material. The dose distribution was modelled using GPUMCD (MC algorithm in Monaco) for an unmodified (or default) material table (DMT) as well as for a modified (or custom) material table (CMT) for both the phantom and patients. Various concentrations ranging between 8 and 157 mg ml-1 were used to represent the gadolinium uptake in the patients GTV. It was assumed that the gadolinium concentration remained the same for the entire course of radiation treatment. Results showed that at the tissue-Gadovist interface, inside the phantom, dose scored using the DMT was 7% lower compared to that using the CMT for 157 mg ml-1 concentration of gadolinium. Dosimetric differences in the case of the patient study were measured using the DVH parameters. D 50% was higher by 6% when the DMT was used compared to the CMT for dose modelling for a gadolinium concentration of 157 mg ml-1. This difference decreased gradually with decreasing concentration of gadolinium. It was concluded that dosimetric differences can be quantified in Monaco if the tumour-gadolinium concentration is more than 23 mg ml-1. If the gadolinium concentration is lower than 23 mg ml-1, then a correction for the presence of gadolinium may not be necessary in the TPS.

WE‐G‐18A‐07: Clinical Evaluation of Normalized Metal Artifact Reduction in KVCT Using MVCT Prior Images (MVCT‐NMAR) Technique in Radiotherapy

PURPOSE To evaluate the metal artifacts in diagnostic kVCT images of patients that are corrected using a normalized metal artifact reduction method with MVCT prior images, MVCT-NMAR. METHODS An MVCTNMAR algorithm was developed and applied to five patients: three with bilateral hip prostheses, one with unilateral hip prosthesis and one with dental fillings. The corrected images were evaluated for visualization of tissue structures and their interfaces, and for radiotherapy dose calculations. They were also compared against the corresponding images corrected by a commercial metal artifact reduction technique, O-MAR, on a Phillips™ CT scanner. RESULTS The use of MVCT images for correcting kVCT images in the MVCT-NMAR technique greatly reduces metal artifacts, avoids secondary artifacts, and makes patient images more useful for correct dose calculation in radiotherapy. These improvements are significant over the commercial correction method, provided the MVCT and kVCT images are correctly registered. The remaining and the secondary artifacts (soft tissue blurring, eroded bones, false bones or air pockets, CT number cupping within the metal) present in O-MAR corrected images are removed in the MVCT-NMAR corrected images. Large dose reduction is possible outside the planning target volume (e.g., 59.2 Gy in comparison to 52.5 Gy in pubic bone) when these MVCT-NMAR corrected images are used in TomoTherapy™ treatment plans, as the corrected images no longer require directional blocks for prostate plans in order to avoid the image artifact regions. CONCLUSION The use of MVCT-NMAR corrected images in radiotherapy treatment planning could improve the treatment plan quality for cancer patients with metallic implants. Moti Raj Paudel is supported by the Vanier Canada Graduate Scholarship, the Endowed Graduate Scholarship in Oncology and the Dissertation Fellowship at the University of Alberta. The authors acknowledge the CIHR operating grant number MOP 53254.

SU-E-I-92: Is Photon Starvation Preventing Metal Artifact Reduction Algorithm From Working in KVCT?

PURPOSE High density/high atomic number metallic objects create shading and streaking metal artifacts in the CT image that can cause inaccurate delineation of anatomical structures or inaccurate radiation dose calculation. A modified iterative maximum-likelihood polychromatic algorithm for CT (mIMPACT) that models the energy response of detectors, photon interaction processes and beam polychromaticity has successfully reduced metal artifacts in MVCT. Our extension of mIMPACT in kVCT did not significantly reduce metal artifacts for high density metal like steel. We hypothesize that photon starvation may result in the measured data in a commercial kVCT imaging beam. METHODS We measured attenuation of a range of steel plate thicknesses, sandwiched between two 12cm thick solid water blocks, using a Phillips Big Bore CTTM scanner in scout acquisition mode with 120kVp and 200mAs. The transmitted signal (y) was normalized to the air scan signal (y0 ) to get attenuation [i.e., ln(y/y0 )] data for a detector. RESULTS Below steel plate thickness of 13.4mm, the variations in measured attenuation as a function of view number are characterized by a quantum noise and show increased attenuation with metal thickness. On or above this thickness the attenuation shows discrete levels in addition to the quantum noise. Some views have saturated attenuation value. The histograms of the measured attenuation for up to 36.7mm of steel show this trend. The detector signal is so small that the quantization levels in the analog to digital (A-to-D) converter are visible, a clear indication of photon starvation. CONCLUSION Photons reaching the kVCT detector after passing through a thick metal plate are either so low in number that the signal measured has large quantum noise, or are completely absorbed inside the plate creating photon starvation. This is un-interpretable by the mIMPACT algorithm and cannot reduce metal artifacts in kVCT for certain realistic thicknesses of steel hip implants. Moti Raj Paudel is supported by the Vanier Canada Graduate Scholarship, the Endowed Graduate Scholarship in Oncology and the Dissertation Fellowship at the University of Alberta. The authors acknowledge the CIHR operating grant number MOP 53254.

SU-DD-A3-04: Evaluation of Metal Artifact Reduction Using MVCT and Model Based Image Reconstruction

Purpose: To evaluate the performance of a model based image reconstruction in reducing metal artifacts in MVCT system in the image of a phantom representing bilateral hip prostheses, and to compare with traditionally reconstructed image.Method and Materials: A cylindrical plexiglass (19 cm diameter) phantom containing two steel inserts (2 cm diameter) was scanned using bremsstrahlung radiation from 6MeV electron beam passing through 4 cm thick solid water in Varian Clinac 2300C with bench‐top MVCT system. Iterative maximum‐likelihood polychromatic algorithm for CT (IMPACT) modified to include pair and triplet production for the MV spectrum was used with air, plexiglass, bone, and iron as base substances. Detector calibrated signal without beam‐hardening correction was used to get filtered back‐projected (FBP) image and used as an initial image in IMPACT. The final image at 1.25 MeV was obtained after 130 iterations. The second image was reconstructed using FBP with traditional signal and beam‐hardening corrections. The quantitative analysis included calculation of average attenuation coefficient in various regions of interest and their comparisons with theoretical value. Results: Visual inspection of the traditional image shows that the region between two steel inserts is dark with additional white streaks. Attenuation coefficients in ROIs fluctuate largely and deviate from theoretical values. In the iteratively reconstructed image, metal artifacts are remarkably reduced leaving behind only faint shadings in between the inserts. The average attenuation coefficients in ROIs were close to theoretical values. Image profile through the steel inserts shows restored uniform background between them. Conclusion: This experiment emphasizes the importance of model based image reconstruction and MVCT system for the metal artifact reduction. The beam‐hardening correction applied in conventional image reconstruction creates severe darkening and white streaks in the image. The complete removal of the metal artifacts requires further modeling and strategies which will be employed in future studies.