A single centre intercomparison between commercial treatment planning systems for 90Y radioembolization using virtual and experimental phantoms.

INTRODUCTION: Dedicated Treatment Planning Systems (TPSs) were developed to personalize 90Y-transarterial radioembolization. This study evaluated the agreement among four commercial TPSs assessing volumes of interest (VOIs) volumes and dose metrics. METHODS: A homogeneous (EH) and an anthropomorphic phantom with hot and cold inserts (EA) filled with 99mTc-pertechnetate were acquired with a SPECT/CT scanner. Their virtual versions (VH and VA, respectively) and a phantom with activity inside a single voxel (VK) were generated by an in-house MATLAB script. Images and delineated VOIs were imported into the TPSs to compute voxel-based absorbed dose distributions with various dose deposition approaches: local deposition method (LDM) and dose kernel convolution (DKC) with/without local density correction (LDC). VOI volumes and mean absorbed doses were assessed against their median value across TPSs. Dose-volume histograms (DVHs) and VK-derived dose profiles were evaluated. RESULTS: Small (<2.1 %) and large (up to 42.4 %) relative volume differences were observed on large (>500 ml) and small VOIs, respectively. Mean absorbed doses relative differences were < 3 % except for small VOIs with steep dose gradients (up to 89.1 % in the VA Cold Sphere VOI). Within the same TPS, LDC negligibly affected the mean absorbed dose, while DKC and LDM showed differences up to 63 %. DHVs were mostly overlapped in experimental phantoms, with some differences in the virtual versions. Dose profiles agreed within 1 %. CONCLUSION: TPSs showed an overall good agreement except for small VOI volumes and mean absorbed doses of VOIs with steep dose gradients. These discrepancies should be considered in the dosimetry uncertainty assessment, thus requiring an appropriate harmonization.

A novel tool for predicting the dose distribution of non-sealed 188 Re (Rhenium) resin in non-melanoma skin cancers (NMSC) patients.

BACKGROUND: High-dose rate brachytherapy using a non-sealed 188 Rhenium resin (188 Re) is a recently approved treatment option for non-melanoma skin cancer (NMSC). The treatment goal is to deliver a personalized absorbed dose to the deepest point of neoplastic infiltration corresponding to the minimal target dose. The treatment consists of the application of a 188 Re-based resin over a plastic foil placed on the target skin surface. However, there is no treatment planning tool to assess the 188 Re activity needed for a personalized treatment. PURPOSE: The paper aims to present a novel Monte Carlo (MC)-based tool for 188 Re-based resin activity and dose calculation, experimentally validated using Gafchromic EBT3 films. METHODS: MC simulations were carried out using FLUKA modeling density and composition of 188 Re resin. The MC-based look up table (LUT) was incorporated in an ad hoc developed tool. The proposed tool allows the personalized calculation of treatment parameters (i.e., activity to be dispensed, the treatment duration, and dose volume histograms), according to the target dimension. The proposed tool was compared using Bland-Altman analysis to the previous calculation approaches conducted using VARSKIN in a retrospective cohort of 76 patients. The tool was validated in ad hoc experimental set ups using a stack of calibrated Gafchromic EBT3 films covered by a plastic film and exposed using a homogenous activity distribution of 188 Re eluate and a heterogeneous activity distribution of 188 Re resin mimic the patient treatment. RESULTS: The agreement between the proposed tool and VARSKIN was evaluated on the investigated cohort with median range of target area, target depth, and treatment time equal to 4.8 [1.0-60.1] cm2 , 1.1 [0.2-3.0] mm, and 70 [21-285] min, with a median range of target dose (Gy) of 23.5 [10-54.9]. The calculated minimal target doses, ranged from 1% to 10% for intermediate target depths (1.2 ± 0.7 mm), while showing significant differences in the estimation of superficial (maximal) target doses. The agreement between MC calculation and measurements at different plans in a stack of Gafchromic EBT3 films was within 10% for both the homogenous and heterogeneous activity distribution of 188 Re. Worst agreements were observed for absorbed doses lower than 0.3 Gy. CONCLUSIONS: Our results support the implementation of our MC-based tool in the practical routine for calculating the 188 Re resin activity and treatment parameters necessary for obtaining the prescribed minimal target dose.

A novel figure of merit to investigate 68Ga PET/CT image quality based on patient weight and lesion size using Q.Clear reconstruction algorithm: A phantom study.

INTRODUCTION: Q.Clear is a Bayesian penalised-likelihood algorithm that uses a ß-value for positron emission tomography(PET)/computed tomography(CT) image reconstruction(IR). Our study proposes a novel figure of merit, named CRBV, to compare the Q.Clear performances using 68Ga PET/CT image with the ordered-subset-expectation-maximization(OSEM) algorithm and to identify the optimal ß-values for these images using two phantoms mimicking normal and overweight patients. METHODS: NEMA IQ phantom with or without a ring of water-filled plastic bags (NEMAstd and NEMAow, respectively) was acquired and reconstructed with OSEM and Q.Clear at various ß-values and minutes/bed position(min/bp). Contrast recovery(CR), background variability(BV) and CRBV were calculated. Highest CRBV values were used to identify optimal ß-value ranges. RESULTS: Q.Clear with 250 ≤ ß ≤ 800 improved CRBV compared to OSEM for all the investigated spheres and acquisition setups. Outside of this range, Q.Clear still outperformed OSEM with few exceptions depending on spheres diameters and phantoms(e.g.,ß-value = 1600 for diameters ≤ 17 mm using the NEMAow phantom). Regarding the CRBV performance for IR optimization, for the 4 min/bp NEMAstd IR, ß-values = 300 ÷ 350 allowed to simultaneously optimize all diameters(except for the 10 mm); for the NEMAow IR, ß-values = 350 ÷ 500 were needed for diameters > 20 mm, while ß-values = 200 ÷ 250 were selected for the remaining diameters. For the 2 min/bp, ß-value = 500 was suitable for diameters > 17 mm in both NEMAstd and NEMAow IR, while for smaller diameters ß-value = 200 and ß-values = 250 ÷ 350 were obtained for NEMAstd and NEMAow, respectively. CONCLUSION: Almost all tested ß-values of Q.Clear improved the CRBV compared to OSEM. In both phantoms, simulating normal and over-weight patients, optimal ß-values were found according to lesion sizes and investigated acquisition times.

A novel tool for motion-related dose inaccuracies reduction in 99mTc-MAA SPECT/CT images for SIRT planning.

INTRODUCTION: In Selective Internal Radiation Therapy (SIRT), 90Y is administered to primary/secondary hepatic lesions. An accurate pre-treatment planning using 99mTc-MAA SPECT/CT allows the assessment of its feasibility and of the activity to be injected. Unfortunately, SPECT/CT suffers from patient-specific respiratory motion which causes artifacts and absorbed dose inaccuracies. In this study, a data-driven solution was developed to correct the respiratory motion. METHODS: The tool realigns the barycenter of SPECT projection images and shifts them to obtain a fine registration with the attenuation map. The tool was validated using a modified dynamic phantom with several breathing patterns. We compared the absorbed dose distributions derived from uncorrected(Dm)/corrected(Dc) images with static ones(Ds) in terms of γ-passing rates, 210 Gy isodose volumes, dose-volume histograms and percentage differences of mean doses (i.e., ΔD¯m and ΔD¯c, respectively). The tool was applied to twelve SIRT patients and the Bland-Altman analysis was performed on mean doses. RESULTS: In the phantom study, the agreement between Dc and Ds was higher (γ-passing rates generally > 90%) than Dm and Ds. The isodose volumes in Dc were closer than Dm to Ds, with differences up to 10% and 30% respectively. A reduction from a median ΔD¯m = -19.3% to ΔD¯c = -0.9%, from ΔD¯m = -42.8% to ΔD¯c = -7.0% and from ΔD¯m = 1586% to ΔD¯c = 47.2% was observed in liver-, tumor- and lungs-like structures. The Bland-Altman analysis on patients showed variations (±50 Gy) and (±4 Gy) between D¯c and D¯m of tumor and lungs, respectively. CONCLUSION: The proposed tool allowed the correction of 99mTc-MAA SPECT/CT images, improving the accuracy of the absorbed dose distribution.

Overview of commercial treatment planning systems for targeted radionuclide therapy.

INTRODUCTION: Targeted Radionuclide Therapy (TRT) is a branch of cancer medicine dealing with the therapeutic use of radioisotopes associated with biological vectors accumulating in the tumors/targets, indicated as Molecular Radiotherapy (MRT), or directly injected into the arteries that supply blood to liver tumour vasculature, indicated as Selective RT (SRT). The aim of this work is to offer a panoramic view on the increasing number of commercially-available TRT treatment planning systems (TPSs). MATERIALS AND METHODS: A questionnaire was sent to manufacturers' representatives. Academic software were not considered. Questions were grouped as follows: general information, clinical workflow, calibration procedure, image processing/reconstruction, image registration and segmentation tools, time-activity curve (TAC) fitting and absorbed dose calculation. RESULTS: All software reported have CE-marking. TPSs were divided between SRT-dedicated software [4] and MRT [5] dosimetry software. In SRT, since no kinetic process is involved, absorbed dose calculation does not require TAC fitting, and image registration is not fully developed in all TPS. All software requires a radionuclide-specific calibration. In SRT, a relative image calibration can be obtained by scaling the counts to a known activity. Automated VOI contouring and rigid/deformable propagation between different acquisitions time-points is implemented in most TPSs, although DICOM export is rare. Different TAC fits are available depending on the number of time-points. Voxel S-value and Local deposition methods are the most frequent dosimetric approaches; dose-voxel kernel convolution and semi-Monte Carlo method are also available. CONCLUSIONS: Available TPSs allows performing personalized dosimetry in clinical practice. Individual variations in methodology/algorithms must be considered in the standardisation/harmonization processes.

A novel tool for assessing the correlation of internal/external markers during SGRT guided stereotactic ablative radiotherapy treatments.

INTRODUCTION: An in-house developed tool was implemented and validated to investigate the skin surface, hepatic dome, and target displacement for stereotactic ablative radiotherapy (SABR) of thoracic/abdominal lesions using a Surface Guided Radiation Therapy (SGRT) system combined with 4D- images. MATERIALS AND METHODS: Fourteen consecutive patients with tumors near the hepatic dome undergoing SABR treatments were analyzed. For each patient, a planning 4D-CT and five 4D-CBCT images were acquired. The C-RAD technology was also used to register/monitor the position of the skin reference point (SRP) as an external marker representative of patient breathing. The 4D images were imported in the developed tool, and the absolute maximum height (Pmax,dome) of the hepatic dome on the ten respiratory phases was semi-automatically detected. Similarly, the contour of the skin surface was extracted in correspondence with the SRP position. The tool has been validated using an ad hoc modified moving phantom with pre-selected amplitudes and numbers of cycles. The Pearson correlation coefficients and Bland-Altman plots were calculated. RESULTS: There was a strong correlation between the skin motion amplitude based on 4D-CBCT and the C-RAD in all the patients (0.90 ± 0.08). Similarly, the mean ± SD of Pearson correlation coefficients of skin and Pmax,dome movements registered by 4D-CT and 4D-CBCT were 0.90 ± 0.05 and 0.94 ± 0.05, respectively. The mean ± SD of Pearson correlation coefficients comparing the skin and Pmax,dome displacements within each imaging modality were 0.88 ± 0.05 and 0.90 ± 0.05 for 4D-CT and 4D-CBCT, respectively. The SRP displacement during the set-up imaging and the treatment delivery were similar in all the investigated patients. Similar results were obtained for the ad hoc modified phantom in the preliminary validation phase. CONCLUSION: The strong correlation between the tumor/ hepatic dome and skin displacements confirms that the SGRT approach can be considered appropriate for intra- and inter-fraction motion management in SABR therapy.

Comprehensive Analysis of Radiomic Datasets by RadAR.

Quantitative analysis of biomedical images, referred to as radiomics, is emerging as a promising approach to facilitate clinical decisions and improve patient stratification. The typical radiomic workflow includes image acquisition, segmentation, feature extraction, and analysis of high-dimensional datasets. While procedures for primary radiomic analyses have been established in recent years, processing the resulting radiomic datasets remains a challenge due to the lack of specific tools for doing so. Here we present RadAR (Radiomics Analysis with R), a new software to perform comprehensive analysis of radiomic features. RadAR allows users to process radiomic datasets in their entirety, from data import to feature processing and visualization, and implements multiple statistical methods for analysis of these data. We used RadAR to analyze the radiomic profiles of more than 850 patients with cancer from publicly available datasets and showed that it was able to recapitulate expected results. These results demonstrate RadAR as a reliable and valuable tool for the radiomics community. SIGNIFICANCE: A new computational tool performs comprehensive analysis of high-dimensional radiomic datasets, recapitulating expected results in the analysis of radiomic profiles of >850 patients with cancer from independent datasets.

A comparative study of physical image quality in digital and synthetic mammography from commercially available mammography systems.

We present a comparison between full field digital mammography and synthetic mammography, performed on several mammography systems from four different manufacturers. The analysis is carried out on both the digital and synthetic images of two commercially available mammography phantoms, and focuses on a set of objective metrics that encode the geometrical appearance of imaging features of diagnostic interest. In particular, we measured sizes and contrasts of several clusters of microcalcification specks, shapes and contrasts of circular masses, and the power spectrum of background regions mimicking the heterogeneous texture of the breast parenchyma. Despite the potential issues of tomosynthesis in terms of image blurring, the synthetic images do not highlight any globally significant differences in the rendering of the details of interest, when compared to the original digital mammograms: relative contrasts are generally preserved, as well as the geometry of broad structures. We conclude that, as far as the considered objective metrics are concerned, the image quality of synthetic mammography does not exhibit significant differences with respect to the one of full field digital mammography, for all the considered systems.

Fully automated VMAT treatment planning for advanced-stage NSCLC patients.

PURPOSE: To develop a fully automated procedure for multicriterial volumetric modulated arc therapy (VMAT) treatment planning (autoVMAT) for stage III/IV non-small cell lung cancer (NSCLC) patients treated with curative intent. MATERIALS AND METHODS: After configuring the developed autoVMAT system for NSCLC, autoVMAT plans were compared with manually generated clinically delivered intensity-modulated radiotherapy (IMRT) plans for 41 patients. AutoVMAT plans were also compared to manually generated VMAT plans in the absence of time pressure. For 16 patients with reduced planning target volume (PTV) dose prescription in the clinical IMRT plan (to avoid violation of organs at risk tolerances), the potential for dose escalation with autoVMAT was explored. RESULTS: Two physicians evaluated 35/41 autoVMAT plans (85%) as clinically acceptable. Compared to the manually generated IMRT plans, autoVMAT plans showed statistically significant improved PTV coverage (V95% increased by 1.1% ± 1.1%), higher dose conformity (R50 reduced by 12.2% ± 12.7%), and reduced mean lung, heart, and esophagus doses (reductions of 0.9 Gy ± 1.0 Gy, 1.5 Gy ± 1.8 Gy, 3.6 Gy ± 2.8 Gy, respectively, all p < 0.001). To render the six remaining autoVMAT plans clinically acceptable, a dosimetrist needed less than 10 min hands-on time for fine-tuning. AutoVMAT plans were also considered equivalent or better than manually optimized VMAT plans. For 6/16 patients, autoVMAT allowed tumor dose escalation of 5-10 Gy. CONCLUSION: Clinically deliverable, high-quality autoVMAT plans can be generated fully automatically for the vast majority of advanced-stage NSCLC patients. For a subset of patients, autoVMAT allowed for tumor dose escalation.