In silico investigation of factors affecting the MV imaging performance of a novel water-equivalent EPID

PurposeA Geant4 model of a novel, water-equivalent electronic portal imaging device (EPID) prototype for radiotherapy imaging and dosimetry utilising an array of plastic scintillating fibres (PSFs) has been developed. Monte Carlo (MC) simulations were performed to quantify the PSF-EPID imaging performance and to investigate design aspects affecting performance for optimisation.MethodsUsing the Geant4 model, the PSF-EPID’s imaging performance for 6 MV photon beams was quantified in terms of its modulation transfer function (MTF), noise power spectrum (NPS) and detective quantum efficiency (DQE). Model parameters, including fibre dimensions, optical cladding reflectivity and scintillation yield, were varied to investigate impact on imaging performance.ResultsThe MC-calculated DQE(0) for the reference PSF-EPID geometry employing 30 mm fibres was approximately nine times greater than values reported for commercial EPIDs. When using 10 mm long fibres, the PSF-EPID DQE(0) was still approximately three times greater than that of a commercial EPID. Increased fibre length, cladding reflectivity and scintillation yield produced the greatest decreases in NPS and increases in DQE.ConclusionsThe potential to develop an optimised next-generation water-equivalent EPID with MV imaging performance at least comparable to commercial EPIDs has been demonstrated. Factors most important for optimising prototype design include fibre length, cladding reflectivity and scintillation yield.     

A new cone-beam computed tomography dosimetry method providing optimal measurement positions: A Monte Carlo study

PurposeThe conventional weighted computed tomography dose index (CTDI_w) may not be suitable for cone-beam computed tomography (CBCT) dosimetry because a cross-sectional dose distribution is angularly inhomogeneous owing to partial angle irradiations. This study was conducted to develop a new dose metric (f(0)^CB_w) for CBCT dosimetry to determine a more accurate average dose in the central cross-sectional plane of a cylindrical phantom using Monte Carlo simulations.MethodsFirst, cross-sectional dose distributions of cylindrical polymethyl methacrylate phantoms over a wide range of phantom diameters (8–40 cm) were calculated for various CBCT scan protocols. Then, by obtaining linear least-squares fits of the full datasets of the cross-sectional dose distributions, the optimal radial positions, which represented measurement positions for the average phantom dose, were determined. Finally, the f(0)^CB_w method was developed by averaging point doses at the optimal radial positions of the phantoms. To demonstrate its validity, the relative differences between the average doses and each dose index value were estimated for the devised f(0)^CB_w, conventional CTDI_w, and Haba’s CTDI_w methods, respectively.ResultsThe relative differences between the average doses and each dose index value were within 4.1%, 16.7%, and 11.9% for the devised, conventional CTDI_w, and Haba’s CTDI_w methods, respectively.ConclusionsThe devised f(0)^CB_w value was calculated by averaging four “point doses” at 90° intervals and the optimal radial positions of the cylindrical phantom. The devised method can estimate the average dose more accurately than the previously developed CTDI_w methods for CBCT dosimetry.     

Generation of clinical 177Lu SPECT/CT images based on Monte Carlo simulation with GATE

PurposePatient-specific dosimetry in MRT relies on quantitative imaging, pharmacokinetic assessment and absorbed dose calculation. The DosiTest project was initiated to evaluate the uncertainties associated with each step of the clinical dosimetry workflow through a virtual multicentric clinical trial. This work presents the generation of simulated clinical SPECT datasets based on GATE Monte Carlo modelling with its corresponding experimental CT image, which can subsequently be processed by commercial image workstations.MethodsThis study considers a therapy cycle of 6.85 GBq ¹⁷⁷Lu-labelled DOTATATE derived from an IAEA-Coordinated Research Project (E23005) on “Dosimetry in Radiopharmaceutical therapy for personalised patient treatment”. Patient images were acquired on a GE Infinia-Hawkeye 4 gamma camera using a medium energy (ME) collimator. Simulated SPECT projections were generated based on experimental time points and validated against experimental SPECT projections using flattened profiles and gamma index. The simulated projections were then incorporated into the patient SPECT/CT DICOM envelopes for processing and their reconstruction within a commercial image workstation.ResultsGamma index passing rate (2% − 1 pixel criteria) between 95 and 98% and average gamma between 0.28 and 0.35 among different time points revealed high similarity between simulated and experimental images. Image reconstruction of the simulated projections was successful on HERMES and Xeleris workstations, a major step forward for the initiation of a multicentric virtual clinical dosimetry trial based on simulated SPECT/CT images.ConclusionsRealistic ¹⁷⁷Lu patient SPECT projections were generated in GATE. These modelled datasets will be circulated to different clinical departments to perform dosimetry in order to assess the uncertainties in the entire dosimetric chain.     

Monte Carlo kilovoltage X-ray tube simulation: A statistical analysis and compact simulation method

IntroductionMonte Carlo (MC) simulations are a powerful tool for improving image quality in X-ray based imaging modalities. An accurate X-ray source model is essential to MC modeling for CBCT but can be difficult to implement on a GPU while maintaining efficiency and memory limitations. A statistical analysis of the photon distribution from a MC X-ray tube simulation is conducted in hopes of building a compact source model.Materials & methodsMC simulations of an X-ray tube were carried out using BEAMnrc. The resulting photons were sorted into four categories: primary, scatter, off-focal radiation (OFR), and both (scatter and OFR). A statistical analysis of the photon components (energy, position, direction) was completed. A novel method for a compact (memory efficient) representation of the PHSP data was implemented and tested using different statistical based linear transformations (PCA, ZCA, ICA), as well as a geometrical transformation.ResultsThe statistical analysis showed all photon groupings had strong correlations between position and direction, with the largest correlation in the primary data. The novel method was successful in compactly representing the primary (error < 2%) and scatter (error < 6%) photon groupings by reducing the component correlations.Discussion & conclusionStatistical linear transforms provide a method of reducing the memory required to accurately simulate an X-ray source in a GPU MC system. If all photon types are required, the proposed method reduces the memory requirements by 3.8 times. When only primary and scatter data is needed, the memory requirement is reduced from gigabytes to kilobytes.     

Quantifying the performance of two different types of commercial software programs for 3D patient dose reconstruction for prostate cancer patients: Machine log files vs. machine log files with EPID images

We clarified the reconstructed 3D dose difference between two different commercial software programs (Mobius3D v2.0 and PerFRACTION v1.6.4).Five prostate cancer patients treated with IMRT (74 Gy/37 Fr) were studied. Log files and cine EPID images were acquired for each fraction. 3D patient dose was reconstructed using log files (Mobius3D) or log files with EPID imaging (PerFRACTION). The treatment planning dose was re-calculated on homogeneous and heterogeneous phantoms, and log files and cine EPID images were acquired. Measured doses were compared with the reconstructed point doses in the phantom. Next, we compared dosimetric metrics (mean dose for PTV, rectum, and bladder) calculated by Mobius3D and PerFRACTION for all fractions from five patients.Dose difference at isocenter between measurement and reconstructed dose for two software programs was within 3.0% in both homogeneous and heterogeneous phantoms. Moreover, the dose difference was larger using skip arc plan than that using full arc plan, especially for PerFRACTION (e.g., dose difference at isocenter for PerFRACTION: 0.34% for full arc plan vs. −4.50% for skip arc plan in patient 1).For patients, differences in dosimetric parameters were within 1% for almost all fractions. PerFRACTION had wider range of dose difference between first fraction and the other fractions than Mobius3D (e.g., maximum difference: 0.50% for Mobius3D vs. 1.85% for PerFRACTION), possibly because EPID may detect some types of MLC positioning errors such as miscalibration errors or mechanical backlash which cannot be detected by log files, or that EPID data might include image acquisition failure and image noise.     

A Monte Carlo software bench for simulation of spectral k-edge CT imaging: Initial results

PurposeSpectral Computed Tomography (SCT) systems equipped with photon counting detectors (PCD) are clinically desired, since such systems provide not only additional diagnostic information but also radiation dose reductions by a factor of two or more. The current unavailability of clinical PCDs makes a simulation of such systems necessary.MethodsIn this paper, we present a Monte Carlo-based simulation of a SCT equipped with a PCD. The aim of this development is to facilitate research on potential clinical applications. Our MC simulator takes into account scattering interactions within the scanned object and has the ability to simulate scans with and without scatter and a wide variety of imaging parameters. To demonstrate the usefulness of such a MC simulator for development of SCT applications, a phantom with contrast targets covering a wide range of clinically significant iodine concentrations is simulated. With those simulations the impact of scatter and exposure on image quality and material decomposition results is investigated.ResultsOur results illustrate that scatter radiation plays a significant role in visual as well as quantitative results. Scatter radiation can reduce the accuracy of contrast agent concentration by up to 15%.ConclusionsWe present a reliable and robust software bench for simulation of SCTs equipped with PCDs.     

Monte Carlo dose in a prosthesis phantom based on exact geometry vs streak artefact contaminated CT data as benchmarked against Gafchromic film measurements

PurposeIn radiotherapy, accurate calculation of patient radiation dose is very important for good clinical outcome. In the presence of metallic implants, the dose calculation accuracy could be compromised by metal artefacts generated in computed tomography (CT) images of patients. This study investigates the influence of metal-induced CT artefacts on MC dose calculations in a pelvic prosthesis phantom.MethodsA pelvic phantom containing unilateral Ti prosthesis was CT-scanned and accurate Hounsfield unit (HU) values were assigned to known materials of the phantom as opposed to HU values produced through the artefact CT images of the phantom. Using the DOSXYZnrc MC code, dose calculations were computed in the phantom model constructed from the original CT images containing the artefacts and artefact-free images made from the exact geometry of the phantom with known materials. The dose calculations were benchmarked against Gafchromic EBT3 film measurements using 15 MeV electron and 10 MV photon beams.ResultsThe average deviations between film and MC dose data decreased from 3 ± 2% to 1 ± 1% and from about 6 ± 2% to 3 ± 1% for the artefact and artefact-free phantom models against film data for the electron and photon fields, respectively.ConclusionsFor the Ti prosthesis phantom, the presence of metal-induced CT artefacts could cause dose inaccuracies of about 3%. Construction of an artefact-free phantom model made from the exact geometry of the phantom with known materials to overcome the effect of artefacts is advantageous compared to using CT data directly of which the exact tissue composition is not well-known.     

Towards component-based validation of GATE: Aspects of the coincidence processor

GATE is public domain software widely used for Monte Carlo simulation in emission tomography. Validations of GATE have primarily been performed on a whole-system basis, leaving the possibility that errors in one sub-system may be offset by errors in others. We assess the accuracy of the GATE PET coincidence generation sub-system in isolation, focusing on the options most closely modeling the majority of commercially available scanners.Independent coincidence generators were coded by teams at Toshiba Medical Research Unit (TMRU) and UC Davis. A model similar to the Siemens mCT scanner was created in GATE. Annihilation photons interacting with the detectors were recorded. Coincidences were generated using GATE, TMRU and UC Davis code and results compared to “ground truth” obtained from the history of the photon interactions. GATE was tested twice, once with every qualified single event opening a time window and initiating a coincidence check (the “multiple window method”), and once where a time window is opened and a coincidence check initiated only by the first single event to occur after the end of the prior time window (the “single window method”). True, scattered and random coincidences were compared. Noise equivalent count rates were also computed and compared.The TMRU and UC Davis coincidence generators agree well with ground truth. With GATE, reasonable accuracy can be obtained if the single window method option is chosen and random coincidences are estimated without use of the delayed coincidence option. However in this GATE version, other parameter combinations can result in significant errors.     

A Monte Carlo method to account for sampling error in multi-species indicators

The usefulness of biodiversity indicators strongly increases if accompanied by measures of uncertainty. In the case of indicators that combine population indices of species, however, the inclusion of the uncertainty of the species indices has shown to be hard to realize, usually due to imperfections in monitoring programmes. Missing values and time series of different lengths preclude the use of analytical approaches, whereas bootstrapping across sites requires the raw abundance data on the site level, which may not always be available. Sometimes bootstrapping across species rather than sites is opted for, but this approach ignores the uncertainty attached to species indices. We developed a method to account for sampling error of species indices in the calculation of multi-species indicators based on Monte Carlo simulation of annual species indices. The construction of confidence intervals enables various trend assessments, like testing for linear or smooth trends, testing for changes between two time points, testing the significance of a suspected change-point and testing for differences between two multi-species indicators. Here, we compare our method with conventional methods and illustrate the benefits of our approach using Dutch breeding bird indicators.     

Image quality characteristics for virtual monoenergetic images using dual-layer spectral detector CT: Comparison with conventional tube-voltage images

PurposeTo investigate the image quality characteristics for virtual monoenergetic images compared with conventional tube-voltage image with dual-layer spectral CT (DLCT).MethodsHelical scans were performed using a first-generation DLCT scanner, two different sizes of acrylic cylindrical phantoms, and a Catphan phantom. Three different iodine concentrations were inserted into the phantom center. The single-tube voltage for obtaining virtual monoenergetic images was set to 120 or 140 kVp. Conventional 120- and 140-kVp images and virtual monoenergetic images (40–200-keV images) were reconstructed from slice thicknesses of 1.0 mm. The CT number and image noise were measured for each iodine concentration and water on the 120-kVp images and virtual monoenergetic images. The noise power spectrum (NPS) was also calculated.ResultsThe iodine CT numbers for the iodinated enhancing materials were similar regardless of phantom size and acquisition method. Compared with the iodine CT numbers of the conventional 120-kVp images, those for the monoenergetic 40-, 50-, and 60-keV images increased by approximately 3.0-, 1.9-, and 1.3-fold, respectively. The image noise values for each virtual monoenergetic image were similar (for example, 24.6 HU at 40 keV and 23.3 HU at 200 keV obtained at 120 kVp and 30-cm phantom size). The NPS curves of the 70-keV and 120-kVp images for a 1.0-mm slice thickness over the entire frequency range were similar.ConclusionVirtual monoenergetic images represent stable image noise over the entire energy spectrum and improved the contrast-to-noise ratio than conventional tube voltage using the dual-layer spectral detector CT.     

Measurement-based model of a wide-bore CT scanner for Monte Carlo dosimetric calculations with GMCTdospp software

The aim of this work was to create a model of a wide-bore Siemens Somatom Sensation Open CT scanner for use with GMCTdospp, which is an EGSnrc-based software tool dedicated for Monte Carlo calculations of dose in CT examinations.The method was based on matching spectrum and filtration to half value layer and dose profile, and thus was similar to the method of Turner et al. (Med. Phys. 36, pp. 2154–2164). Input data on unfiltered beam spectra were taken from two sources: the TASMIP model and IPEM Report 78. Two sources of HVL data were also used, namely measurements and documentation. Dose profile along the fan-beam was measured with Gafchromic RTQA-1010 (QA+) film. Two-component model of filtration was assumed: bow-tie filter made of aluminum with 0.5 mm thickness on central axis, and flat filter made of one of four materials: aluminum, graphite, lead, or titanium.Good agreement between calculations and measurements was obtained for models based on the measured values of HVL. Doses calculated with GMCTdospp differed from the doses measured with pencil ion chamber placed in PMMA phantom by less than 5%, and root mean square difference for four tube potentials and three positions in the phantom did not exceed 2.5%. The differences for models based on HVL values from documentation exceeded 10%. Models based on TASMIP spectra and IPEM78 spectra performed equally well.     

Monte Carlo simulation of damascene electroplating: effects of additives

The influence of additives on the filling process of via holes in damascene electroplating is investigated with use of a kinetic Monte Carlo method. The basic system is the solid-by-solid model for crystal growth which includes the vacancy formation during the growth of thin film. Three kinds of additives are included in the model to control the local surface growth rate. Inhibitors and levelers have the effect of preventing the deposition, while accelerators increase the local growth rate. Levelers are modeled to stick to the tips of the surface. We performed a series of simulations by changing the parameters which characterize the additives to see their influence on the filling mechanism. It is shown that void-free filling is possible by the combination of the effects of the additives.     

Evaluation of three scatter correction methods based on estimation of photopeak scatter spectrum in SPECT imaging: A simulation study

Three practical methods for scatter correction of Tc-99m SPECT images are evaluated. Among these, two methods, three-energy window (TEW) methods using the trapezoidal and triangular approximations, have been described previously by investigators, and a new approximation is offered in this work. The SIMIND (SIMulation of Imaging Nuclear Detectors) Monte Carlo program is used to simulate a line source placed at on-axis and 5 cm off-axis locations, a cold-sphere/hot-background phantom, a hot-sphere/cold-background phantom, and a more clinically realistic NCAT (Nonuniform Rational B-spline-based CArdiac-Torso) phantom. For evaluation of these methods, the scatter line-spread functions and scatter fractions for the on- and off-axis line source, image contrast, signal-to-noise ratio and relative noise for the cold spheres, and recovery coefficient for the hot spheres of different diameters are compared. For the NCAT phantom, a line profile through a slice of the reconstructed image is considered before and after scatter correction, and also image contrast defined by this profile is used to compare the correction methods. The results of this study indicate that for the line source simulation the scatter fractions obtained from the proposed method are a better estimation of true scatter fractions. Also, for both the sphere simulation and NCAT simulation, the proposed method improves the image contrast as compared to the two other methods.     

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.