Three-Dimensional Margin as a Predictor of Local Tumor Progression after Microwave Ablation: Intraprocedural versus 4-8-Week Postablation Assessment.

PURPOSE: To evaluate the prognostic accuracy of intraprocedural and 4-8-week (current standard) post-microwave ablation zone (AZ) and margin assessments for prediction of local tumor progression (LTP) using 3-dimensional (3D) software. MATERIALS AND METHODS: Data regarding 100 colorectal liver metastases (CLMs) in 75 patients were collected from 2 prospective fluorodeoxyglucose positron emission tomography (PET)/computed tomography (CT)-guided microwave ablation (MWA) trials. The target CLMs and theoretical 5- and 10-mm margins were segmented and registered intraprocedurally and at 4-8 weeks after MWA contrast-enhanced CT (or magnetic resonance [MR] imaging) using the same methodology and 3D software. Tumor and 5- and 10-mm minimal margin (MM) volumes not covered by the AZ were defined as volumes of insufficient coverage (VICs). The intraprocedural and 4-8-week post-MWA VICs were compared as predictors of LTP using receiver operating characteristic curve analysis. RESULTS: The median follow-up time was 19.6 months (interquartile range, 7.97-36.5 months). VICs for 5- and 10-mm MMs were predictive of LTP at both time assessments. The highest accuracy for the prediction of LTP was documented with the intra-ablation 5-mm VIC (area under the curve [AUC], 0.78; 95% confidence interval, 0.66-0.89). LTP for a VIC of 6-10-mm margin category was 11.4% compared with 4.3% for >10-mm margin category (P < .001). CONCLUSIONS: A 3D 5-mm MM is a critical endpoint of thermal ablation, whereas optimal local tumor control is noted with a 10-mm MM. Higher AUCs for prediction of LTP were achieved for intraprocedural evaluation than for the 4-8-week postablation 3D evaluation of the AZ.

The OpenGATE ecosystem for Monte Carlo simulation in medical physics.

This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.

3D margin assessment predicts local tumor progression after ablation of colorectal cancer liver metastases.

OBJECTIVE: To determine the feasibility and prognostic value of 3D measuring of the ablation margins using a dedicated image registration software. METHODS: This retrospective study included 104 colorectal liver metastases in 68 consecutive patients that underwent microwave ablation between 08/2012 and 08/2019. The minimal ablation margin (MM) was measured in 2D using anatomic landmarks on contrast enhanced CT(CECT) 4-8 weeks post-ablation, and in 3D using an image registration software and immediate post-ablation CECT. Local tumor progression (LTP) was assessed by imaging up to 24 months after ablation. A blinded interventional radiologist provided feedback on the possibility of additional ablation after examining the 3D-margin measurements. RESULTS: The 3D-margin assessment was completed in 79/104 (76%) tumors without the need for target manipulation. In 25/104 (24%) tumors, manipulation was required due to image misregistration. LTP was observed in 40/104 (38.5%) tumors: 92.5% vs 7.5% for those with margin <5mm vs ≥5mm, respectively (p = 0.0001). The 2D and 3D-assessments identified margin <5mm in 17/104 (16%), and in 74/104 (71%) ablated tumors, respectively (p < 0.01). The sensitivity and specificity of the 3D software for predicting LTP was 93% (37/40) and 42% (27/64), respectively. Additional ablation to achieve a MM of 5 mm would have been offered in 26/37 cases if the 3D-margin assessment was available intraoperatively. CONCLUSION: Image registration software can measure ablation margins and detect MM under 5 mm intraoperatively, with significantly higher sensitivity than the 2D technique using landmarks on the post-ablation CECT. The identification of a margin under 5 mm is strongly associated with LTP.

Advanced Monte Carlo simulations of emission tomography imaging systems with GATE.

Built on top of the Geant4 toolkit, GATE is collaboratively developed for more than 15 years to design Monte Carlo simulations of nuclear-based imaging systems. It is, in particular, used by researchers and industrials to design, optimize, understand and create innovative emission tomography systems. In this paper, we reviewed the recent developments that have been proposed to simulate modern detectors and provide a comprehensive report on imaging systems that have been simulated and evaluated in GATE. Additionally, some methodological developments that are not specific for imaging but that can improve detector modeling and provide computation time gains, such as Variance Reduction Techniques and Artificial Intelligence integration, are described and discussed.

PETSTEP: Generation of synthetic PET lesions for fast evaluation of segmentation methods.

PURPOSE: This work describes PETSTEP (PET Simulator of Tracers via Emission Projection): a faster and more accessible alternative to Monte Carlo (MC) simulation generating realistic PET images, for studies assessing image features and segmentation techniques. METHODS: PETSTEP was implemented within Matlab as open source software. It allows generating three-dimensional PET images from PET/CT data or synthetic CT and PET maps, with user-drawn lesions and user-set acquisition and reconstruction parameters. PETSTEP was used to reproduce images of the NEMA body phantom acquired on a GE Discovery 690 PET/CT scanner, and simulated with MC for the GE Discovery LS scanner, and to generate realistic Head and Neck scans. Finally the sensitivity (S) and Positive Predictive Value (PPV) of three automatic segmentation methods were compared when applied to the scanner-acquired and PETSTEP-simulated NEMA images. RESULTS: PETSTEP produced 3D phantom and clinical images within 4 and 6 min respectively on a single core 2.7 GHz computer. PETSTEP images of the NEMA phantom had mean intensities within 2% of the scanner-acquired image for both background and largest insert, and 16% larger background Full Width at Half Maximum. Similar results were obtained when comparing PETSTEP images to MC simulated data. The S and PPV obtained with simulated phantom images were statistically significantly lower than for the original images, but led to the same conclusions with respect to the evaluated segmentation methods. CONCLUSIONS: PETSTEP allows fast simulation of synthetic images reproducing scanner-acquired PET data and shows great promise for the evaluation of PET segmentation methods.

32P brachytherapy conformal source model RIC-100 for high-dose-rate treatment of superficial disease: Monte Carlo calculations, diode measurements, and clinical implementation.

PURPOSE: A novel (32)P brachytherapy source has been in use at our institution intraoperatively for temporary radiation therapy of the spinal dura and other localized tumors. We describe the dosimetry and clinical implementation of the source. METHODS AND MATERIALS: Dosimetric evaluation for the source was done with a complete set of MCNP5 Monte Carlo calculations preceding clinical implementation. In addition, the depth dose curve and dose rate were measured by use of an electron field diode to verify the Monte Carlo calculations. Calibration procedures using the diode in a custom-designed phantom to provide an absolute dose calibration and to check dose uniformity across the source area for each source before treatment were established. RESULTS: Good agreement was established between the Monte Carlo calculations and diode measurements. Quality assurance measurements results are provided for about 100 sources used to date. Clinical source calibrations were usually within 10% of manufacturer specifications. Procedures for safe handling of the source are described. DISCUSSION: Clinical considerations for using the source are discussed.

Influence of photon energy cuts on PET Monte Carlo simulation results.

PURPOSE: The purpose of this work is to study the influence of photon energy cuts on the results of positron emission tomography (PET) Monte Carlo (MC) simulations. METHODS: MC simulations of PET scans of a box phantom and the NEMA image quality phantom are performed for 32 photon energy cut values in the interval 0.3-350 keV using a well-validated numerical model of a PET scanner. The simulations are performed with two MC codes, egs_pet and GEANT4 Application for Tomographic Emission (GATE). The effect of photon energy cuts on the recorded number of singles, primary, scattered, random, and total coincidences as well as on the simulation time and noise-equivalent count rate is evaluated by comparing the results for higher cuts to those for 1 keV cut. To evaluate the effect of cuts on the quality of reconstructed images, MC generated sinograms of PET scans of the NEMA image quality phantom are reconstructed with iterative statistical reconstruction. The effects of photon cuts on the contrast recovery coefficients and on the comparison of images by means of commonly used similarity measures are studied. RESULTS: For the scanner investigated in this study, which uses bismuth germanate crystals, the transport of Bi X(K) rays must be simulated in order to obtain unbiased estimates for the number of singles, true, scattered, and random coincidences as well as for an unbiased estimate of the noise-equivalent count rate. Photon energy cuts higher than 170 keV lead to absorption of Compton scattered photons and strongly increase the number of recorded coincidences of all types and the noise-equivalent count rate. The effect of photon cuts on the reconstructed images and the similarity measures used for their comparison is statistically significant for very high cuts (e.g., 350 keV). The simulation time decreases slowly with the increase of the photon cut. CONCLUSIONS: The simulation of the transport of characteristic x rays plays an important role, if an accurate modeling of a PET scanner system is to be achieved. The simulation time decreases slowly with the increase of the cut which, combined with the accuracy loss at high cuts, means that the usage of high photon energy cuts is not recommended for the acceleration of MC simulations.

Pencil beam approach for correcting the energy dependence artifact in film dosimetry for IMRT verification.

The higher sensitivity to low-energy scattered photons of radiographic film compared to water can lead to significant dosimetric error when the beam quality varies significantly within a field. Correcting for this artifact will provide greater accuracy for intensity modulated radiation therapy (IMRT) verification dosimetry. A procedure is developed for correction of the film energy-dependent response by creating a pencil beam kernel within our treatment planning system to model the film response specifically. Film kernels are obtained from EGSnrc Monte Carlo simulations of the dose distribution from a 1 mm diameter narrow beam in a model of the film placed at six depths from 1.5 to 40 cm in polystyrene and solid water phantoms. Kernels for different area phantoms (50 x 50 cm2 and 25 x 25 cm2 polystyrene and 30 x 30 cm2 solid water) are produced. The Monte Carlo calculated kernel is experimentally verified with film, ion chamber and thermoluminescent dosimetry (TLD) measurements in polystyrene irradiated by a narrow beam. The kernel is then used in convolution calculations to, predict the film response in open and IMRT fields. A 6 MV photon beam and Kodak XV2 film in a polystyrene phantom are selected to test the method as they are often used in practice and can result in large energy-dependent artifacts. The difference in dose distributions calculated with the film kernel and the water kernel is subtracted from film measurements to obtain a practically film artifact free IMRT dose distribution for the Kodak XV2 film. For the points with dose exceeding 5 cGy (11% of the peak dose) in a large modulated field and a film measurement inside a large polystyrene phantom at depth of 10 cm, the correction reduces the fraction of pixels for which the film dose deviates from dose to water by more than 5% of the mean film dose from 44% to 6%.

Validation of GATE Monte Carlo simulations of the GE Advance/Discovery LS PET scanners.

The recently developed GATE (GEANT4 application for tomographic emission) Monte Carlo package, designed to simulate positron emission tomography (PET) and single photon emission computed tomography (SPECT) scanners, provides the ability to model and account for the effects of photon noncollinearity, off-axis detector penetration, detector size and response, positron range, photon scatter, and patient motion on the resolution and quality of PET images. The objective of this study is to validate a model within GATE of the General Electric (GE) Advance/Discovery Light Speed (LS) PET scanner. Our three-dimensional PET simulation model of the scanner consists of 12 096 detectors grouped into blocks, which are grouped into modules as per the vendor's specifications. The GATE results are compared to experimental data obtained in accordance with the National Electrical Manufactures Association/Society of Nuclear Medicine (NEMA/SNM), NEMA NU 2-1994, and NEMA NU 2-2001 protocols. The respective phantoms are also accurately modeled thus allowing us to simulate the sensitivity, scatter fraction, count rate performance, and spatial resolution. In-house software was developed to produce and analyze sinograms from the simulated data. With our model of the GE Advance/Discovery LS PET scanner, the ratio of the sensitivities with sources radially offset 0 and 10 cm from the scanner's main axis are reproduced to within 1% of measurements. Similarly, the simulated scatter fraction for the NEMA NU 2-2001 phantom agrees to within less than 3% of measured values (the measured scatter fractions are 44.8% and 40.9 +/- 1.4% and the simulated scatter fraction is 43.5 +/- 0.3%). The simulated count rate curves were made to match the experimental curves by using deadtimes as fit parameters. This resulted in deadtime values of 625 and 332 ns at the Block and Coincidence levels, respectively. The experimental peak true count rate of 139.0 kcps and the peak activity concentration of 21.5 kBq/cc were matched by the simulated results to within 0.5% and 0.1% respectively. The simulated count rate curves also resulted in a peak NECR of 35.2 kcps at 10.8 kBq/cc compared to 37.6 kcps at 10.0 kBq/cc from averaged experimental values. The spatial resolution of the simulated scanner matched the experimental results to within 0.2 mm.

Predicting energy response of radiographic film in a 6 MV x-ray beam using Monte Carlo calculated fluence spectra and absorbed dose.

The advantage of radiographic film is that it allows two-dimensional, high-resolution dose measurement. While there is concern over its photon energy dependence, these problems are considered acceptable within small fields, where the scatter component is small. The application of film dosimetry to intensity modulated radiotherapy (IMRT) raises additional concern since the primary fluence may vary significantly within the field. The varying primary fluence in combination with a large scatter fraction, present for large fields and large depths, causes the spectrum at various points within the IMRT field to differ from the spectrum in the uniform fields typically used for calibrating the film. As a result, significant artifacts are introduced in the measured dose distribution. The purpose of this work is to quantify and develop a method to correct for these artifacts. Two approaches based on Monte Carlo (MC) simulations are examined. In the first method, the film artifact, as quantified by film and ion chamber output measurements in uniform square fields, is derived from the MC calculated ratio of absorbed doses to film and to water. In the second method, the measured film artifact is correlated with MC calculated photon spectra, revealing a strong correlation between the measured artifact and the "scatter"-to-"primary" ratio, defined by the ratio of the number of photons below to the number of photons above 0.1 MeV, independent of field size and depth. These methods are evaluated in high- and low-dose regions of a large intensity-modulated field created with a central block. The spectral approach is also tested with a clinical IMRT field. The absorbed dose method accurately corrects the measured film dose in the open part of the field and in points under the block and outside the field. The dose error is reduced from as much as 16% of the open field dose to less than 1%, as verified with an ion chamber. The spectral method accurately corrects the measured film dose in the open region of the centrally blocked field, but does not fully correct for the film artifact for points under the block and outside the field, where the spectrum is substantially different. Applied to the clinical field, the corrected film measurement shows good agreement with data obtained with a two-dimensional diode array.