SlicerRT: radiation therapy research toolkit for 3D Slicer

Unsupervised model for structure segmentation applied to brain computed tomography

This article presents an unsupervised method for segmenting brain computed tomography scans. The proposed methodology involves image feature extraction and application of similarity and continuity constraints to generate segmentation maps of the anatomical head structures. Specifically designed for real-world datasets, this approach applies a spatial continuity scoring function tailored to the desired number of structures. The primary objective is to assist medical experts in diagnosis by identifying regions with specific abnormalities. Results indicate a simplified and accessible solution, reducing computational effort, training time, and financial costs. Moreover, the method presents potential for expediting the interpretation of abnormal scans, thereby impacting clinical practice. This proposed approach might serve as a practical tool for segmenting brain computed tomography scans, and make a significant contribution to the analysis of medical images in both research and clinical settings.

Enabling extracorporeal ultrasound imaging with the da Vinci robot for transoral robotic surgery: a feasibility study

PURPOSE

Transoral robotic surgery (TORS) is a challenging procedure due to its small workspace and complex anatomy. Ultrasound (US) image guidance has the potential to improve surgical outcomes, but an appropriate method for US probe manipulation has not been defined. This study evaluates using an additional robotic (4th) arm on the da Vinci Surgical System to perform extracorporeal US scanning for image guidance in TORS.

METHODS

A stereoscopic imaging system and da Vinci-compatible US probe attachment were developed to enable control of the extracorporeal US probe from the surgeon console. The prototype was compared to freehand US by nine operators in three tasks on a healthy volunteer: (1) identification of the common carotid artery, (2) carotid artery scanning, and (3) identification of the submandibular gland. Operator workload and user experience were evaluated using a questionnaire.

RESULTS

The robotic US tasks took longer than freehand US tasks (2.09x longer; p = 0.001 ) and had higher operator workload (2.12x higher; p = 0.004 ). However, operator-rated performance was closer (avg robotic/avg freehand = 0.66; p = 0.017 ), and scanning performance measured by MRI-US average Hausdorff distance provided no statistically significant difference.

CONCLUSION

Extracorporeal US scanning for intraoperative US image guidance is a convenient approach for providing the surgeon direct control over the US image plane during TORS, with little modification to the existing operating room workflow. Although more time-consuming and higher operator workload, several methods have been identified to address these limitations.

To degas or not to degas? The effectiveness of pharmaceuticals in reducing the impact of bowel volume variations during liver SBRT treatment

The indications for liver stereotactic body radiation therapy (SBRT) continue to expand in the management of liver cancer due to the improved rates of local control with acceptable normal tissue toxicity. Changes in internal anatomy, such as the bowel, may negatively impact the precision of treatment delivery of SBRT liver treatment by influencing daily image matching. Institutions have developed various approaches to promoting bowel volume consistency. One such strategy is the administration of pharmaceuticals. The administration of pharmaceuticals, such as Simethicone, has been adopted by the Princess Alexandra Hospital Radiation Oncology Department (ROPAIR) as a method to promote consistency in the amount of bowel gas observed in liver cancer patients. This case series examines a group of patients treated at ROPAIR with liver SBRT to determine whether current practices effectively reduce the impact of bowel volume variations for liver cancer patients. Initial observations from this hypothesis generating research suggest potential improved consistency of the small bowel's anatomical position for liver SBRT patients who were administered Simethicone (Bowel bag dice similarity coefficient - Simethicone group = 0.79-0.92, Standard group = 0.24-0.93). However, it appeared that this strategy alone may not be entirely effective achieving consistency in the amount of bowel gas present throughout the duration of treatment. Further investigation into the refinement of liver SBRT pre-treatment preparation is therefore recommended.

Improvement in male pelvis magnetic resonance image contouring following radiologist-delivered training

INTRODUCTION

The magnetic resonance linear accelerator (MRL) combines both magnetic resonance imaging and a linear accelerator, allowing for daily treatment adaptation. This study aimed to assess the impact of radiologist-delivered training in magnetic resonance (MR) contouring of relevant structures within the male pelvis.

METHODS

Two radiation oncologists, two radiation oncology registrars and seven radiation therapists completed contouring on 10 male pelvis MR datasets both pre- and post-training. A 2-hour MR anatomy training session was delivered by a radiologist, who also provided the 'gold standard' contours. The pre- and post-training contours were compared against the gold standard with Dice similarity coefficient (DSC) and Hausdorff distances calculated; and the pre- and post-confidence scores and timing were compared.

RESULTS

The improvement in DSC were significant in prostate, rectum and seminal vesicles, with a post-training median DSC of 0.87 ± 0.06, 0.92 ± 0.04 and 0.80 ± 0.14, respectively. The median Hausdorff improved with a median of 1.46 ± 0.78 mm, 0.52 ± 0.32 mm and 1.11 ± 0.86 mm for prostate, rectum and seminal vesicles, respectively. Bladder concordance was high both pre- and post-training. Urethra contours improved post-training, however, remained difficult to contour with a median post-DSC of 0.51 ± 0.24. Overall, confidence scoring improved (P < 0.001) and timing decreased by an average of 4.4 ± 16.4 min post-training.

CONCLUSION

Radiologist-delivered training improved concordance of male pelvis contouring on MR datasets. Further work is required in the identification of urethra on MRs. These findings are of importance in the MRL adaptive workflow.

An investigation into the impact of volumetric rescanning and fractionation treatment on dose homogeneity in liver cancer proton therapy

The Pencil Beam Scanning (PBS) technique in modern particle therapy offers a highly conformal dose distribution but poses challenges due to the interplay effect, an interaction between respiration-induced organ movement and PBS. This study evaluates the effectiveness of different volumetric rescanning strategies in mitigating this effect in liver cancer proton therapy. We used a Geant4-based Monte Carlo simulation toolkit, 'TOPAS,' and an image registration toolbox, 'Elastix,' to calculate 4D dose distributions from 5 patients' four-dimensional computed tomography (4DCT). We analyzed the homogeneity index (HI) value of the Clinical Tumor Volume (CTV) at different rescan numbers and treatment times. Our results indicate that dose homogeneity stabilizes at a low point after a week of treatment, implying that both rescanning and fractionation treatments help mitigate the interplay effect. Notably, an increase in the number of rescans doesn't significantly reduce the mean dose to normal tissue but effectively prevents high localized doses to tissue adjacent to the CTV. Rescanning techniques, based on statistical averaging, require no extra equipment or patient cooperation, making them widely accessible. However, the number of rescans, tumor location, diaphragm movement, and treatment fractionation significantly influence their effectiveness. Therefore, deciding the number of rescans should involve considering the number of beams, treatment fraction size, and total delivery time to avoid unnecessary treatment extension without significant clinical benefits. The results showed that 2-3 rescans are more clinically suitable for liver cancer patients undergoing proton therapy.

Histopathologically confirmed radiation-induced damage of the brain - an in-depth analysis of radiation parameters and spatio-temporal occurrence

BACKGROUND

Radiation-induced damage (RID) after radiotherapy (RT) of primary brain tumors and metastases can be challenging to clinico-radiographically distinguish from tumor progression. RID includes pseudoprogression and radiation necrosis; the latter being irreversible and often associated with severe symptoms. While histopathology constitutes the diagnostic gold standard, biopsy-controlled clinical studies investigating RID remain limited. Whether certain brain areas are potentially more vulnerable to RID remains an area of active investigation. Here, we analyze histopathologically confirmed cases of RID in relation to the temporal and spatial dose distribution.

METHODS

Histopathologically confirmed cases of RID after photon-based RT for primary or secondary central nervous system malignancies were included. Demographic, clinical, and dosimetric data were collected from patient records and treatment planning systems. We calculated the equivalent dose in 2 Gy fractions (EQD22) and the biologically effective dose (BED2) for normal brain tissue (α/β ratio of 2 Gy) and analyzed the spatial and temporal distribution using frequency maps.

RESULTS

Thirty-three patients were identified. High-grade glioma patients (n = 18) mostly received one normofractionated RT series (median cumulative EQD22 60 Gy) to a large planning target volume (PTV) (median 203.9 ccm) before diagnosis of RID. Despite the low EQD22 and BED2, three patients with an accelerated hyperfractionated RT developed RID. In contrast, brain metastases patients (n = 15; 16 RID lesions) were often treated with two or more RT courses and with radiosurgery or fractionated stereotactic RT, resulting in a higher cumulative EQD22 (median 162.4 Gy), to a small PTV (median 6.7 ccm). All (n = 34) RID lesions occurred within the PTV of at least one of the preceding RT courses. RID in the high-grade glioma group showed a frontotemporal distribution pattern, whereas, in metastatic patients, RID was observed throughout the brain with highest density in the parietal lobe. The cumulative EQD22 was significantly lower in RID lesions that involved the subventricular zone (SVZ) than in lesions without SVZ involvement (median 60 Gy vs. 141 Gy, p = 0.01).

CONCLUSIONS

Accelerated hyperfractionated RT can lead to RID despite computationally low EQD22 and BED2 in high-grade glioma patients. The anatomical location of RID corresponded to the general tumor distribution of gliomas and metastases. The SVZ might be a particularly vulnerable area.

A systematic approach for calibrating a Monte Carlo code to a treatment planning system for obtaining dose, LET, variable proton RBE and out-of-field dose

Objective. While integration of variable relative biological effectiveness (RBE) has not reached full clinical implementation, the importance of having the ability to recalculate proton treatment plans in a flexible, dedicated Monte Carlo (MC) code cannot be understated . Here we provide a step-wise method for calibrating dose from a MC code to a treatment planning system (TPS), to obtain required parameters for calculating linear energy transfer (LET), variable RBE and in general enabling clinical realistic research studies beyond the capabilities of a TPS.Approach. Initially, Pristine Bragg peaks (PBP) were calculated in both the Eclipse TPS and the FLUKA MC code. A rearranged Bortfeld energy-range relation was applied to the initial energy of the beam to fine-tune the range of the MC code at 80% dose level distal to the PBP. The energy spread was adapted by dividing the TPS range by the MC range for dose level 80%-20% distal to the PBP. Density and relative proton stopping power were adjusted by comparing the TPS and MC for different Hounsfield units. To find the relationship of dose per primary particle from the MC to dose per monitor unit in the TPS, integration was applied to the area of the Bragg curve. The calibration was validated for spread-out Bragg peaks (SOBP) in water and patient treatment plans. Following the validation, variable RBE were calculated using established models.Main results.The PBPs ranges were within ±0.3mm threshold, and a maximum of 5.5% difference for the SOBPs was observed. The patient validation showed excellent dose agreement between the TPS and MC, with the greatest differences for the lung tumor patient.Significance. Aprocedure for calibrating a MC code to a TPS was developed and validated. The procedure enables MC-based calculation of dose, LET, variable RBE, advanced (secondary) particle tracking and more from treatment plans.

A Tool to Integrate Electrophysiological Mapping for Cardiac Radioablation of Ventricular Tachycardia

Purpose

Cardiac radioablation is an emerging therapy for recurrent ventricular tachycardia. Electrophysiology (EP) data, including electroanatomic maps (EAM) and electrocardiographic imaging (ECGI), provide crucial information for defining the arrhythmogenic target volume. The absence of standardized workflows and software tools to integrate the EP maps into a radiation planning system limits their use. This study developed a comprehensive software tool to enable efficient utilization of the mapping for cardiac radioablation treatment planning.

Methods and Materials

The tool, HeaRTmap, is a Python-scripted plug-in module on the open-source 3D Slicer software platform. HeaRTmap is able to import EAM and ECGI data and visualize the maps in 3D Slicer. The EAM is translated into a 3D space by registration with cardiac magnetic resonance images (MRI) or computed tomography (CT). After the scar area is outlined on the mapping surface, the tool extracts and extends the annotated patch into a closed surface and converts it into a structure set associated with the anatomic images. The tool then exports the structure set and the images as The Digital Imaging and Communications in Medicine Standard in Radiotherapy for a radiation treatment planning system to import. Overlapping the scar structure on simulation CT, a transmural target volume is delineated for treatment planning.

Results

The tool has been used to transfer Ensite NavX EAM data into the Varian Eclipse treatment planning system in radioablation on 2 patients with ventricular tachycardia. The ECGI data from CardioInsight was retrospectively evaluated using the tool to derive the target volume for a patient with left ventricular assist device, showing volumetric matching with the clinically used target with a Dice coefficient of 0.71.

Conclusions

HeaRTmap smoothly fuses EP information from different mapping systems with simulation CT for accurate definition of radiation target volume. The efficient integration of EP data into treatment planning potentially facilitates the study and adoption of the technique.

Prior information guided auto-segmentation of clinical target volume of tumor bed in postoperative breast cancer radiotherapy

BACKGROUND

Accurate delineation of clinical target volume of tumor bed (CTV-TB) is important but it is also challenging due to surgical effects and soft tissue contrast. Recently a few auto-segmentation methods were developed to improve the process. However, those methods had comparatively low segmentation accuracy. In this study the prior information was introduced to aid auto-segmentation of CTV-TB based on a deep-learning model.

METHODS

To aid the delineation of CTV-TB, the tumor contour on preoperative CT was transformed onto postoperative CT via deformable image registration. Both original and transformed tumor contours were used for prior information in training an auto-segmentation model. Then, the CTV-TB contour on postoperative CT was predicted by the model. 110 pairs of preoperative and postoperative CT images were used with a 5-fold cross-validation strategy. The predicted contour was compared with the clinically approved contour for accuracy evaluation using dice similarity coefficient (DSC) and Hausdorff distance.

RESULTS

The average DSC of the deep-learning model with prior information was improved than the one without prior information (0.808 vs. 0.734, P < 0.05). The average DSC of the deep-learning model with prior information was higher than that of the traditional method (0.808 vs. 0.622, P < 0.05).

CONCLUSIONS

The introduction of prior information in deep-learning model can improve segmentation accuracy of CTV-TB. The proposed method provided an effective way to automatically delineate CTV-TB in postoperative breast cancer radiotherapy.

Influence of air mapping errors on the dosimetric accuracy of prostate CBCT-guided online adaptive radiation therapy

PURPOSE

CBCT-guided online adaptive radiotherapy (oART) plans presently utilize daily synthetic CTs (sCT) that are automatically generated using deformable registration algorithms. These algorithms may have poor performance at reproducing variable volumes of gas present during treatment. Therefore, we have analyzed the air mapping error between the daily CBCTs and the corresponding sCT and explored its dosimetric effect on oART plan calculation.

METHODS

Abdominopelvic air volume was contoured on both the daily CBCT images and the corresponding synthetic images for 207 online adaptive pelvic treatments. Air mapping errors were tracked over all fractions. For two case studies representing worst case scenarios, dosimetric effects of air mapping errors were corrected in the sCT images using the daily CBCT air contours, then recalculating dose. Dose volume histogram statistics and 3D gamma passing rates were used to compare the original and air-corrected sCT-based dose calculations.

RESULTS

All analyzed patients showed observable air pocket contour differences between the sCT and the CBCT images. The largest air volume difference observed in daily CBCT images for a given patient was 276.3 cc, a difference of more than 386% compared to the sCT. For the two case studies, the largest observed change in DVH metrics was a 2.6% reduction in minimum PTV dose, with all other metrics varying by less than 1.5%. 3D gamma passing rates using 1%/1 mm criteria were above 90% when comparing the uncorrected and corrected dose distributions.

CONCLUSION

Current CBCT-based oART workflow can lead to inaccuracies in the mapping of abdominopelvic air pockets from daily CBCT to the sCT images used for the optimization and calculation of the adaptive plan. Despite the large observed mapping errors, the dosimetric effects of such differences on the accuracy of the adapted plan dose calculation are unlikely to cause differences greater than 3% for prostate treatments.

Novel tongue-positioning device to reduce tongue motions during radiation therapy for head and neck cancer: Geometric and dosimetric evaluation

This study aimed to assess the performance of a tongue-positioning device in interfractional tongue position reproducibility by cone-beam computed tomography (CBCT). Fifty-two patients treated with radiation therapy (RT) while using a tongue positioning device were included in the study. All patients were treated with 28 or 30 fractions using the volumetric modulated arc therapy technique. CBCT images were acquired at the 1st, 7th, 11th, 15th, 19th, 23th, and 27th fractions. Tongues on planning computed tomography (pCT) and CBCT images were contoured in the treatment planning system. Geometric differences in the tongue between pCT and CBCT were assessed by the Dice similarity coefficient (DSC) and averaged Hausdorff distance (AHD). Two-dimensional in vivo measurements using radiochromic films were performed in 13 patients once a week during sessions. The planned dose distributions were compared with the measured dose distributions using gamma analysis with criteria of 3%/3 mm. In all patients, the mean DSC at the 1st fraction (pCT versus 1st CBCT) was 0.80 while the mean DSC at the 27th fraction (pCT versus 27th CBCT) was 0.77 with statistical significance (p-value = 0.015). There was no statistically significant difference in DSC between the 1st fraction and any other fraction, except for the 27th fraction. There was statistically significant difference in AHD between the 1st fraction and the 19th, 23th, and 27th fractions (p-value < 0.05). In vivo measurements showed an average gamma passing rate of 90.54%. There was no significant difference between measurements at the 1st week and those at other weeks. The tongue geometry during RT was compared between pCT and CBCT. In conclusion, the novel tongue-positioning device was found to minimize interfractional variations in position and shape of the tongue.

Combining ultrasound radiomics, complete blood count, and serum biochemical biomarkers for diagnosing intestinal disorders in cats using machine learning

This retrospective analytical observational cohort study aimed to model and predict the classification of feline intestinal diseases from segmentations of a transverse section from small intestine ultrasound (US) image, complete blood count (CBC), and serum biochemical profile data using a variety of machine-learning approaches. In 149 cats from three institutions, images were obtained from cats with biopsy-confirmed small cell epitheliotropic lymphoma (lymphoma), inflammatory bowel disease (IBD), no pathology ("healthy"), and other conditions (warrant a biopsy for further diagnosis). CBC, blood serum chemistry, small intestinal ultrasound, and small intestinal biopsy were obtained within a 2-week interval. CBC and serum biomarkers and radiomic features were combined for modeling. Four classification schemes were investigated: (1) normal versus abnormal; (2) warranting or not warranting a biopsy; (3) lymphoma, IBD, healthy, or other conditions; and (4) lymphoma, IBD, or other conditions. Two feature selection methods were used to identify the top 3, 5, 10, and 20 features, and six machine learning models were trained. The average (95% CI) performance of models for all combinations of features, numbers of features, and types of classifiers was 0.886 (0.871-0.912) for Model 1 (normal vs. abnormal), 0.751 (0.735-0.818) for Model 2 (biopsy vs. no biopsy), 0.504 (0.450-0.556) for Model 3 (lymphoma, IBD, healthy, or other), and 0.531 (0.426-0.589), for Model 4 (lymphoma, IBD, or other). Our findings suggest model accuracies above 0.85 can be achieved in Model 1 and 2, and that including CBC and biochemistry data with US radiomics data did not significantly improve accuracy in our models.

Extension of matRad with a modified microdosimetric kinetic model for carbon ion treatment planning: Comparison with Monte Carlo calculation

BACKGROUND

Treatment planning is essential for in silico particle therapy studies. matRad is an open-source research treatment planning system (TPS) based on the local effect model, which is a type of relative biological effectiveness (RBE) model.

PURPOSE

This study aims to implement a microdosimetric kinetic model (MKM) in matRad and develop an automation algorithm for Monte Carlo (MC) dose recalculation using the TOPAS code. In addition, we provide the developed MKM extension as open-source tool for users.

METHODS

Carbon beam data were generated using TOPAS MC pencil beam irradiation. We parameterized the TOPAS MC beam data with a double-Gaussian fit and modeled the integral depth doses and lateral spot profiles in the range of 100-430 MeV/u. To implement the MKM, the specific energy data table for Z = 1-6 and integrated depth-specific energy data were acquired based on the Kiefer-Chatterjee track structure and TOPAS MC simulation, respectively. Generic data were integrated into matRad, and treatment planning was performed based on these data. The optimized plan parameters were automatically converted into MC simulation input. Finally, the matRad TPS and TOPAS MC simulations were compared using the RBE-weighted dose calculation results. A comparison was made for three geometries: homogeneous water phantom, inhomogeneous phantom, and patient.

RESULTS

The RBE-weighted dose (DRBE ) distribution agreed with TOPAS MC within 1.8% for all target sizes for the homogeneous phantom. For the inhomogeneous phantom, the relative difference in the range of 80% of the prescription dose in the distal fall-off region (R80) between the matRad TPS and TOPAS MC was 0.6% (1.1 mm). DRBE between the TPS and the MC was within 4.0%. In the patient case, the difference in the dose-volume histogram parameters for the target volume between the TPS and the MC was less than 2.7%. The relative difference in R80 was 0.7% (1.2 mm).

CONCLUSIONS

The MKM was successfully implemented in matRad TPS, and the RBE-weighted dose was comparable to that of TOPAS MC. The MKM-implemented matRad was released as an open-source tool. Further investigations with MC simulations can be conducted using this tool, providing a good option for carbon ion research.

AAPM WGDCAB Report 372: A joint AAPM, ESTRO, ABG, and ABS report on commissioning of model-based dose calculation algorithms in brachytherapy

The introduction of model-based dose calculation algorithms (MBDCAs) in brachytherapy provides an opportunity for a more accurate dose calculation and opens the possibility for novel, innovative treatment modalities. The joint AAPM, ESTRO, and ABG Task Group 186 (TG-186) report provided guidance to early adopters. However, the commissioning aspect of these algorithms was described only in general terms with no quantitative goals. This report, from the Working Group on Model-Based Dose Calculation Algorithms in Brachytherapy, introduced a field-tested approach to MBDCA commissioning. It is based on a set of well-characterized test cases for which reference Monte Carlo (MC) and vendor-specific MBDCA dose distributions are available in a Digital Imaging and Communications in Medicine-Radiotherapy (DICOM-RT) format to the clinical users. The key elements of the TG-186 commissioning workflow are now described in detail, and quantitative goals are provided. This approach leverages the well-known Brachytherapy Source Registry jointly managed by the AAPM and the Imaging and Radiation Oncology Core (IROC) Houston Quality Assurance Center (with associated links at ESTRO) to provide open access to test cases as well as step-by-step user guides. While the current report is limited to the two most widely commercially available MBDCAs and only for 192 Ir-based afterloading brachytherapy at this time, this report establishes a general framework that can easily be extended to other brachytherapy MBDCAs and brachytherapy sources. The AAPM, ESTRO, ABG, and ABS recommend that clinical medical physicists implement the workflow presented in this report to validate both the basic and the advanced dose calculation features of their commercial MBDCAs. Recommendations are also given to vendors to integrate advanced analysis tools into their brachytherapy treatment planning system to facilitate extensive dose comparisons. The use of the test cases for research and educational purposes is further encouraged.

Commissioning dose computation model for proton source in pencil beam scanning therapy by convolution neural networks

Objective. Proton source model commissioning (PSMC) is critical for ensuring accurate dose calculation in pencil beam scanning (PBS) proton therapy using Monte Carlo (MC) simulations. PSMC aims to match the calculated dose to the delivered dose. However, commissioning the 'nominal energy' and 'energy spread' parameters in PSMC can be challenging, as these parameters cannot be directly obtained from solving equations. To efficiently and accurately commission the nominal energy and energy spread in a proton source model, we developed a convolution neural network (CNN) named 'PSMC-Net.'Methods. The PSMC-Net was trained separately for 33 energies (E, 70-225 MeV with a step of 5 MeV plus 226.09 MeV). For eachE, a dataset was generated consisting of 150 source model parameters (15 nominal energies ∈ [E,E+ 1.5 MeV], ten spreads ∈ [0, 1]) and the corresponding 150 MC integrated depth doses (IDDs). Of these 150 data pairs, 130 were used for training the network, 10 for validation, and 10 for testing.Results. The source model, built by 33 measured IDDs and 33 PSMC-Nets (cost 0.01 s), was used to compute the MC IDDs. The gamma passing rate (GPRs, 1 mm/1%) between MC and measured IDDs was 99.91 ± 0.12%. However, when no commissioning was made, the corresponding GPR was reduced to 54.11 ± 22.36%, highlighting the tremendous significance of our CNN commissioning method. Furthermore, the MC doses of a spread-out Bragg peak and 20 patient PBS plans were also calculated, and average 3D GPRs (2 mm/2% with a 10% threshold) were 99.89% and 99.96 ± 0.06%, respectively.Significance. We proposed a nova commissioning method of the proton source model using CNNs, which made the PSMC process easy, efficient, and accurate.

Creation of waterproof, TLD probes for dose measurements to validate image-based radiopharmaceutical therapy dosimetry workflow

Voxel-level dosimetry based on nuclear medicine images offers patient-specific personalization of radiopharmaceutical therapy (RPT) treatments. Clinical evidence is emerging demonstrating improvements in treatment precision in patients when voxel-level dosimetry is used compared to MIRD. Voxel-level dosimetry requires absolute quantification of activity concentrations in the patient, but images from SPECT/CT scanners are not quantitative and require calibration using nuclear medicine phantoms. While phantom studies can validate a scanner's ability to recover activity concentrations, these studies provide only a surrogate for the true metric of interest: absorbed doses. Measurements using thermoluminescent dosimeters (TLDs) are a versatile and accurate method of measuring absorbed dose. In this work, a TLD probe was manufactured that can fit into currently available nuclear medicine phantoms for the measurement of absorbed dose of RPT agents. Next, 748 MBq of I-131 was administered to a 16 ml hollow source sphere placed in a 6.4 L Jaszczak phantom in addition to six TLD probes, each holding 4 TLD-100 1 × 1 × 1 mm TLD-100 (LiF:Mg,Ti) microcubes. The phantom then underwent a SPECT/CT scan in accordance with a standard SPECT/CT imaging protocol for I-131. The SPECT/CT images were then input into a Monte Carlo based RPT dosimetry platform named RAPID and a three dimensional dose distribution in the phantom was estimated. Additionally, a GEANT4 benchmarking scenario (denoted 'idealized') was created using a stylized representation of the phantom. There was good agreement for all six probes, the differences between measurement and RAPID ranged between -5.5% and 0.9%. The difference between the measured and the idealized GEANT4 scenario was calculated and ranged from -4.3% and -20.5%. This work demonstrates good agreement between TLD measurements and RAPID. In addition, it introduces a novel TLD probe that can be easily introduced into clinical nuclear medicine workflows to provide QA of image-based dosimetry for RPT treatments.

IA PULMONARY VEIN ATLAS FOR RADIOTHERAPY PLANNING

BACKGROUND AND PURPOSE

Cardiac arrhythmia is a recognised potential complication of thoracic radiotherapy, but the responsible cardiac substructures for arrhythmogenesis have not been identified. Arrhythmogenic tissue is commonly located in the pulmonary veins (PVs) of cardiology patients with arrhythmia, however these structures are not currently considered organs-at-risk during radiotherapy planning. A standardised approach to their delineation was developed and evaluated.

MATERIALS AND METHODS

The gross and radiological anatomy relevant to atrial fibrillation was derived from cardiology and radiology literature by a multidisciplinary team. A region of interest and contouring instructions for radiotherapy computed tomography scans were iteratively developed and subsequently evaluated. Radiation oncologists (n=5) and radiation technologists (n=2) contoured the PVs on the four-dimensional planning datasets of five patients with locally advanced lung cancer treated with 1.8-2.75 Gy fractions. Contours were compared to reference contours agreed by the researchers using geometric and dosimetric parameters.

RESULTS

The mean dose to the PVs was 35% prescription dose. Geometric and dosimetric similarity of the observer contours with reference contours was fair, with an overall mean Dice of 0.80 ± 0.02. The right superior PV (mean DSC 0.83 ± 0.02) had better overlap than the left (mean DSC 0.80 ± 0.03), but the inferior PVs were equivalent (mean DSC of 0.78). The mean difference in mean dose was 0.79 Gy ± 0.71 (1.46% ± 1.25).

CONCLUSION

A PV atlas with multidisciplinary approval led to reproducible delineation for radiotherapy planning, supporting the utility of the atlas in future clinical radiotherapy cardiotoxicity research encompassing arrhythmia endpoints.

CBCT-based synthetic CT generated using CycleGAN with HU correction for adaptive radiotherapy of nasopharyngeal carcinoma

This study aims to utilize a hybrid approach of phantom correction and deep learning for synthesized CT (sCT) images generation based on cone-beam CT (CBCT) images for nasopharyngeal carcinoma (NPC). 52 CBCT/CT paired images of NPC patients were used for model training (41), validation (11). Hounsfield Units (HU) of the CBCT images was calibrated by a commercially available CIRS phantom. Then the original CBCT and the corrected CBCT (CBCT_cor) were trained separately with the same cycle generative adversarial network (CycleGAN) to generate SCT1 and SCT2. The mean error and mean absolute error (MAE) were used to quantify the image quality. For validations, the contours and treatment plans in CT images were transferred to original CBCT, CBCT_cor, SCT1 and SCT2 for dosimetric comparison. Dose distribution, dosimetric parameters and 3D gamma passing rate were analyzed. Compared with rigidly registered CT (RCT), the MAE of CBCT, CBCT_cor, SCT1 and SCT2 were 346.11 ± 13.58 HU, 145.95 ± 17.64 HU, 105.62 ± 16.08 HU and 83.51 ± 7.71 HU, respectively. Moreover, the average dosimetric parameter differences for the CBCT_cor, SCT1 and SCT2 were 2.7% ± 1.4%, 1.2% ± 1.0% and 0.6% ± 0.6%, respectively. Using the dose distribution of RCT images as reference, the 3D gamma passing rate of the hybrid method was significantly better than the other methods. The effectiveness of CBCT-based sCT generated using CycleGAN with HU correction for adaptive radiotherapy of nasopharyngeal carcinoma was confirmed. The image quality and dose accuracy of SCT2 were outperform the simple CycleGAN method. This finding has great significance for the clinical application of adaptive radiotherapy for NPC.

PyRaDiSe: A Python package for DICOM-RT-based auto-segmentation pipeline construction and DICOM-RT data conversion

BACKGROUND AND OBJECTIVE

Despite fast evolution cycles in deep learning methodologies for medical imaging in radiotherapy, auto-segmentation solutions rarely run in clinics due to the lack of open-source frameworks feasible for processing DICOM RT Structure Sets. Besides this shortage, available open-source DICOM RT Structure Set converters rely exclusively on 2D reconstruction approaches leading to pixelated contours with potentially low acceptance by healthcare professionals. PyRaDiSe, an open-source, deep learning framework independent Python package, addresses these issues by providing a framework for building auto-segmentation solutions feasible to operate directly on DICOM data. In addition, PyRaDiSe provides profound DICOM RT Structure Set conversion and processing capabilities; thus, it applies also to auto-segmentation-related tasks, such as dataset construction for deep learning model training.

METHODS

The PyRaDiSe package follows a holistic approach and provides DICOM data handling, deep learning model inference, pre-processing, and post-processing functionalities. The DICOM data handling allows for highly automated and flexible handling of DICOM image series, DICOM RT Structure Sets, and DICOM registrations, including 2D-based and 3D-based conversion from and to DICOM RT Structure Sets. For deep learning model inference, extending given skeleton classes is straightforwardly achieved, allowing for employing any deep learning framework. Furthermore, a profound set of pre-processing and post-processing routines is included that incorporate partial invertibility for restoring spatial properties, such as image origin or orientation.

RESULTS

The PyRaDiSe package, characterized by its flexibility and automated routines, allows for fast deployment and prototyping, reducing efforts for auto-segmentation pipeline implementation. Furthermore, while deep learning model inference is independent of the deep learning framework, it can easily be integrated into famous deep learning frameworks such as PyTorch or Tensorflow. The developed package has successfully demonstrated its capabilities in a research project at our institution for organs-at-risk segmentation in brain tumor patients. Furthermore, PyRaDiSe has shown its conversion performance for dataset construction.

CONCLUSIONS

The PyRaDiSe package closes the gap between data science and clinical radiotherapy by enabling deep learning segmentation models to be easily transferred into clinical research practice. PyRaDiSe is available on https://github.com/ubern-mia/pyradise and can be installed directly from the Python Package Index using pip install pyradise.

Large scale crowdsourced radiotherapy segmentations across a variety of cancer anatomic sites

Clinician generated segmentation of tumor and healthy tissue regions of interest (ROIs) on medical images is crucial for radiotherapy. However, interobserver segmentation variability has long been considered a significant detriment to the implementation of high-quality and consistent radiotherapy dose delivery. This has prompted the increasing development of automated segmentation approaches. However, extant segmentation datasets typically only provide segmentations generated by a limited number of annotators with varying, and often unspecified, levels of expertise. In this data descriptor, numerous clinician annotators manually generated segmentations for ROIs on computed tomography images across a variety of cancer sites (breast, sarcoma, head and neck, gynecologic, gastrointestinal; one patient per cancer site) for the Contouring Collaborative for Consensus in Radiation Oncology challenge. In total, over 200 annotators (experts and non-experts) contributed using a standardized annotation platform (ProKnow). Subsequently, we converted Digital Imaging and Communications in Medicine data into Neuroimaging Informatics Technology Initiative format with standardized nomenclature for ease of use. In addition, we generated consensus segmentations for experts and non-experts using the Simultaneous Truth and Performance Level Estimation method. These standardized, structured, and easily accessible data are a valuable resource for systematically studying variability in segmentation applications.

Dosiomics and radiomics to predict pneumonitis after thoracic stereotactic body radiotherapy and immune checkpoint inhibition

Introduction

Pneumonitis is a relevant side effect after radiotherapy (RT) and immunotherapy with checkpoint inhibitors (ICIs). Since the effect is radiation dose dependent, the risk increases for high fractional doses as applied for stereotactic body radiation therapy (SBRT) and might even be enhanced for the combination of SBRT with ICI therapy. Hence, patient individual pre-treatment prediction of post-treatment pneumonitis (PTP) might be able to support clinical decision making. Dosimetric factors, however, use limited information and, thus, cannot exploit the full potential of pneumonitis prediction.

Methods

We investigated dosiomics and radiomics model based approaches for PTP prediction after thoracic SBRT with and without ICI therapy. To overcome potential influences of different fractionation schemes, we converted physical doses to 2 Gy equivalent doses (EQD2) and compared both results. In total, four single feature models (dosiomics, radiomics, dosimetric, clinical factors) were tested and five combinations of those (dosimetric+clinical factors, dosiomics+radiomics, dosiomics+dosimetric+clinical factors, radiomics+dosimetric+clinical factors, radiomics+dosiomics+dosimetric+clinical factors). After feature extraction, a feature reduction was performed using pearson intercorrelation coefficient and the Boruta algorithm within 1000-fold bootstrapping runs. Four different machine learning models and the combination of those were trained and tested within 100 iterations of 5-fold nested cross validation.

Results

Results were analysed using the area under the receiver operating characteristic curve (AUC). We found the combination of dosiomics and radiomics features to outperform all other models with AUCradiomics+dosiomics, D = 0.79 (95% confidence interval 0.78-0.80) and AUCradiomics+dosiomics, EQD2 = 0.77 (0.76-0.78) for physical dose and EQD2, respectively. ICI therapy did not impact the prediction result (AUC ≤ 0.5). Clinical and dosimetric features for the total lung did not improve the prediction outcome.

Conclusion

Our results suggest that combined dosiomics and radiomics analysis can improve PTP prediction in patients treated with lung SBRT. We conclude that pre-treatment prediction could support clinical decision making on an individual patient basis with or without ICI therapy.

Evaluation of MRI-only based online adaptive radiotherapy of abdominal region on MR-linac

PURPOSE

A hybrid magnetic resonance linear accelerator (MRL) can perform magnetic resonance imaging (MRI) with high soft-tissue contrast to be used for online adaptive radiotherapy (oART). To obtain electron densities needed for the oART dose calculation, a computed tomography (CT) is often deformably registered to MRI. Our aim was to evaluate an MRI-only based synthetic CT (sCT) generation as an alternative to the deformed CT (dCT)-based oART in the abdominal region.

METHODS

The study data consisted of 57 patients who were treated on a 0.35 T MRL system mainly for abdominal tumors. Simulation MRI-CT pairs of 43 patients were used for training and validation of a prototype convolutional neural network sCT-generation algorithm, based on HighRes3DNet, for the abdominal region. For remaining test patients, sCT images were produced from simulation MRIs and daily MRIs. The dCT-based plans were re-calculated on sCT with identical calculation parameters. The sCT and dCT were compared in terms of geometric agreement and calculated dose.

RESULTS

The mean and one standard deviation of the geometric agreement metrics over dCT-sCT-pairs were: mean error of 8 ± 10 HU, mean absolute error of 49 ± 10 HU, and Dice similarity coefficient of 55 ± 12%, 60 ± 5%, and 82 ± 15% for bone, fat, and lung tissues, respectively. The dose differences between the sCT and dCT-based dose for planning target volumes were 0.5 ± 0.9%, 0.6 ± 0.8%, and 0.5 ± 0.8% at D_2% , D_50% , and D_98% in physical dose and 0.8 ± 1.4%, 0.8 ± 1.2%, and 0.6 ± 1.1% in biologically effective dose (BED). For organs-at-risk, the dose differences of all evaluated dose-volume histogram points were within [-4.5%, 7.8%] and [-1.1 Gy, 3.5 Gy] in both physical dose and BED.

CONCLUSIONS

The geometric agreement metrics were within typically reported values and most average relative dose differences were within 1%. Thus, an MRI-only sCT-based approach is a promising alternative to the current clinical practice of the abdominal oART on MRL.

Applying Multi-Metric Deformable Image Registration for Dose Accumulation in Combined Cervical Cancer Radiotherapy

(1) Purpose: Challenges remain in dose accumulation for cervical cancer radiotherapy combined with external beam radiotherapy (EBRT) and brachytherapy (BT) as there are many large and complex organ deformations between different treatments. This study aims to improve deformable image registration (DIR) accuracy with the introduction of multi-metric objectives for dose accumulation of EBRT and BT. (2) Materials and methods: Twenty cervical cancer patients treated with EBRT (45-50 Gy/25 fractions) and high-dose-rate BT (≥20 Gy in 4 fractions) were included for DIR. The multi-metric DIR algorithm included an intensity-based metric, three contour-based metrics, and a penalty term. Nonrigid B-spine transformation was used to transform the planning CT images from EBRT to the first BT, with a six-level resolution registration strategy. To evaluate its performance, the multi-metric DIR was compared with a hybrid DIR provided by commercial software. The DIR accuracy was measured by the Dice similarity coefficient (DSC) and Hausdorff distance (HD) between deformed and reference organ contours. The accumulated maximum dose of 2 cc (D_2cc) of the bladder and rectum was calculated and compared to simply addition of D_2cc from EBRT and BT (ΔD_2cc). (3) Results: The mean DSC of all organ contours for the multi-metric DIR were significantly higher than those for the hybrid DIR (p ≤ 0.011). In total, 70% of patients had DSC > 0.8 using the multi-metric DIR, while 15% of patients had DSC > 0.8 using the commercial hybrid DIR. The mean ΔD_2cc of the bladder and rectum for the multi-metric DIR were 3.25 ± 2.29 and 3.54 ± 2.02 Gy_EQD2, respectively, whereas those for the hybrid DIR were 2.68 ± 2.56 and 2.32 ± 3.25 Gy_EQD2, respectively. The multi-metric DIR resulted in a much lower proportion of unrealistic D_2cc than the hybrid DIR (2.5% vs. 17.5%). (4) Conclusions: Compared with the commercial hybrid DIR, the introduced multi-metric DIR significantly improved the registration accuracy and resulted in a more reasonable accumulated dose distribution.

Advantages of a Training Course for Surgical Planning in Virtual Reality for Oral and Maxillofacial Surgery: Crossover Study

BACKGROUND

As an integral part of computer-assisted surgery, virtual surgical planning (VSP) leads to significantly better surgery results, such as for oral and maxillofacial reconstruction with microvascular grafts of the fibula or iliac crest. It is performed on a 2D computer desktop screen (DS) based on preoperative medical imaging. However, in this environment, VSP is associated with shortcomings, such as a time-consuming planning process and the requirement of a learning process. Therefore, a virtual reality (VR)-based VSP application has great potential to reduce or even overcome these shortcomings due to the benefits of visuospatial vision, bimanual interaction, and full immersion. However, the efficacy of such a VR environment has not yet been investigated.

OBJECTIVE

This study aimed to demonstrate the possible advantages of a VR environment through a substep of VSP, specifically the segmentation of the fibula (calf bone) and os coxae (hip bone), by conducting a training course in both DS and VR environments and comparing the results.

METHODS

During the training course, 6 novices were taught how to use a software application in a DS environment (3D Slicer) and in a VR environment (Elucis) for the segmentation of the fibula and os coxae, and they were asked to carry out the maneuvers as accurately and quickly as possible. Overall, 13 fibula and 13 os coxae were segmented for each participant in both methods (VR and DS), resulting in 156 different models (78 fibula and 78 os coxae) per method (VR and DS) and 312 models in total. The individual learning processes in both environments were compared using objective criteria (time and segmentation performance) and self-reported questionnaires. The models resulting from the segmentation were compared mathematically (Hausdorff distance and Dice coefficient) and evaluated by 2 experienced radiologists in a blinded manner.

RESULTS

A much faster learning curve was observed for the VR environment than the DS environment (β=.86 vs β=.25). This nearly doubled the segmentation speed (cm³/min) by the end of training, leading to a shorter time (P<.001) to reach a qualitative result. However, there was no qualitative difference between the models for VR and DS (P=.99). The VR environment was perceived by participants as more intuitive and less exhausting, and was favored over the DS environment.

CONCLUSIONS

The more rapid learning process and the ability to work faster in the VR environment could save time and reduce the VSP workload, providing certain advantages over the DS environment.

Intra- and interobserver assessments of intestinal wall thickness and segmentations from transverse sections of feline abdominal ultrasound images

Measurements of intestinal wall thicknesses from ultrasound imaging (US) are routinely used to support diagnoses of intestinal disorders in cats, however published studies describing observer agreement are currently lacking. The aim of this retrospective, observer agreement study was to quantify inter- and intraobserver repeatability and agreement in the measurement of intestinal wall layer thicknesses and the segmentation of transverse sections of small intestines in US images of 20 cats. Intestinal wall layer thickness measurements of the mucosa, submucosa, muscularis, serosa layer, and total thickness of these layers were performed on five cats with small cell epitheliotropic lymphoma, five with inflammatory bowel disease, and 10 with other conditions. Thickness measurements and the segmentation encompassing the serosa layer were obtained from five observers four times non-sequentially. The average standard deviation in thickness measurements (95% confidence interval) in the mucosa, submucosa, muscularis, serosa, and total thickness were 0.35 (0.07-0.95), 0.24 (0.07-0.52), 0.22 (0.06-0.49), 0.20 (0.05-0.49), and 0.57 (0.11-1.60) mm, respectively. The average intraclass correlation coefficients, which estimates the degree of consistency in thickness measurements and segmentation areas for each observer, ranged from 0.355 to 0.870 and 0.850 to 0.993, respectively. The interclass correlation coefficient, which estimates the degree of consistency when measuring a thickness or segmentation area over all observers ranged from 0.115 to 0.753, and 0.811 to 0.902, respectively. The overall average Dice Coefficient, which estimates the extent of overlap of the segmentations for all observers was 0.957 (0.933 to 0.972). Our results suggest segmentations of small intestines have a higher interobserver agreement than measurements of intestinal wall thicknesses.

CrossMoDA 2021 challenge: Benchmark of cross-modality domain adaptation techniques for vestibular schwannoma and cochlea segmentation

Domain Adaptation (DA) has recently been of strong interest in the medical imaging community. While a large variety of DA techniques have been proposed for image segmentation, most of these techniques have been validated either on private datasets or on small publicly available datasets. Moreover, these datasets mostly addressed single-class problems. To tackle these limitations, the Cross-Modality Domain Adaptation (crossMoDA) challenge was organised in conjunction with the 24th International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI 2021). CrossMoDA is the first large and multi-class benchmark for unsupervised cross-modality Domain Adaptation. The goal of the challenge is to segment two key brain structures involved in the follow-up and treatment planning of vestibular schwannoma (VS): the VS and the cochleas. Currently, the diagnosis and surveillance in patients with VS are commonly performed using contrast-enhanced T1 (ceT₁) MR imaging. However, there is growing interest in using non-contrast imaging sequences such as high-resolution T2 (hrT₂) imaging. For this reason, we established an unsupervised cross-modality segmentation benchmark. The training dataset provides annotated ceT₁ scans (N=105) and unpaired non-annotated hrT₂ scans (N=105). The aim was to automatically perform unilateral VS and bilateral cochlea segmentation on hrT₂ scans as provided in the testing set (N=137). This problem is particularly challenging given the large intensity distribution gap across the modalities and the small volume of the structures. A total of 55 teams from 16 countries submitted predictions to the validation leaderboard. Among them, 16 teams from 9 different countries submitted their algorithm for the evaluation phase. The level of performance reached by the top-performing teams is strikingly high (best median Dice score - VS: 88.4%; Cochleas: 85.7%) and close to full supervision (median Dice score - VS: 92.5%; Cochleas: 87.7%). All top-performing methods made use of an image-to-image translation approach to transform the source-domain images into pseudo-target-domain images. A segmentation network was then trained using these generated images and the manual annotations provided for the source image.

3D slicer-based calculation of hematoma irregularity index for predicting hematoma expansion in intracerebral hemorrhage

BACKGROUND

Irregular hematoma is considered as a risk sign of hematoma expansion. The aim of this study was to quantify hematoma irregularity with computed tomography based on 3D Slicer.

METHODS

Patients with spontaneous intracerebral hemorrhage who underwent an initial and subsequent non-contrast computed tomography (CT) at a single medical center between January 2019 to January 2020 were retrospectively identified. The Digital Imaging and Communication in Medicine (DICOM) standard images were loaded into the 3D Slicer, and the surface area (S) and volume (V) of hematoma were calculated. The hematoma irregularity index (HII) was defined as [Formula: see text]. Logistic regression analyses and receiver operating characteristic (ROC) curve analysis were performed to assess predictive performance of HII.

RESULTS

The enrolled patients were divided into those with hematoma enlargement (n = 36) and those without the enlargement (n = 57). HII in hematoma expansion group was 130.4 (125.1-140.0), and the index in non-enlarged hematoma group was 118.6 (113.5-122.3). There was significant difference in HII between the two groups (P < 0.01). Multivariate logistic regression analysis revealed that the HII was significantly associated with hematoma expansion before (odds ratio = 1.203, 95% confidence interval [CI], 1.115-1.298; P < 0.001) and after adjustment for age, hematoma volume, Glasgow Coma Scale score (odds ratio = 1.196, 95% CI, 1.102-1.298, P < 0.001). The area under the ROC curve was 0.86 (CI, 0.78-0.93, P < 0.01), and the best cutoff of HII for predicting hematoma growth was 123.8.

CONCLUSION

As a quantitative indicator of irregular hematoma, HII can be calculated using the 3D Slicer. And the HII was independently correlated with hematoma expansion.

Utilizing 3D Slicer to incorporate tomographic images into GATE Monte Carlo simulation for personalized dosimetry in yttrium-90 radioembolization

PURPOSE

Monte Carlo (MC) simulation is an important technique that can help design advanced and challenging experimental setups. GATE (Geant4 application for tomographic emission) is a useful simulation toolkit for applications in nuclear medicine. Transarterial radioembolization is a treatment for liver cancer, where microspheres embedded with yttrium-90 (⁹⁰ Y) are administered intra-arterially to the tumor. Personalized dosimetry for this treatment may provide higher dosimetry accuracy compared to the conventional partition model (PM) calculation. However, incorporation of three-dimensional tomographic input data into MC simulation is an intricate process. In this article, 3D Slicer, free and open-source software, was utilized for the incorporation of patient tomographic images into GATE to demonstrate the feasibility of personalized dosimetry in hepatic radioembolization with ⁹⁰ Y.

METHODS

In this article, the steps involved in importing, segmenting, and registering tomographic images using 3D Slicer were thoroughly described, before importing them into GATE for MC simulation. The absorbed doses estimated using GATE were then compared with that of PM. SlicerRT, a 3D Slicer extension, was then used to visualize the isodose from the MC simulation.

RESULTS

A workflow diagram consisting of all the steps taken in the utilization of 3D Slicer for personalized dosimetry in ⁹⁰ Y radioembolization has been presented in this article. In comparison to the MC simulation, the absorbed doses to the tumor and normal liver were overestimated by PM by 105.55% and 20.23%, respectively, whereas for lungs, the absorbed dose estimated by PM was underestimated by 25.32%. These values were supported by the isodose distribution obtained via SlicerRT, suggesting the presence of beta particles outside the volumes of interest. These findings demonstrate the importance of personalized dosimetry for a more accurate absorbed dose estimation compared to PM.

CONCLUSION

The methodology provided in this study can assist users (especially students or researchers who are new to MC simulation) in navigating intricate steps required in the importation of tomographic data for MC simulation. These steps can also be utilized for other radiation therapy related applications, not necessarily limited to internal dosimetry.

The application of multiple metrics in deformable image registration for target volume delineation of breast tumor bed

BACKGROUND AND PURPOSE

For postoperative breast cancer patients, deformable image registration (DIR) is challenged due to the large deformations and non-correspondence caused by tumor resection and clip insertion. To deal with it, three metrics (fiducial-, region-, and intensity-based) were jointly used in DIR algorithm for improved accuracy.

MATERIALS AND METHODS

Three types of metrics were combined to form a single-objective function in DIR algorithm. Fiducial-based metric was used to minimize the distance between the corresponding point sets of two images. Region-based metric was used to improve the overlap between the corresponding areas of two images. Intensity-based metric was used to maximize the correlation between the corresponding voxel intensities of two images. The two CT images, one before surgery and the other after surgery, were acquired from the same patient in the same radiotherapy treatment position. Twenty patients who underwent breast-conserving surgery and postoperative radiotherapy were enrolled in this study.

RESULTS

For target registration error, the difference between the proposed and the conventional registration methods was statistically significant for soft tissue (2.06 vs. 7.82, p = 0.00024 < 0.05) and body boundary (3.70 vs. 6.93, p = 0.021 < 0.05). For visual assessment, the proposed method achieved better matching result for soft tissue and body boundary.

CONCLUSIONS

Comparing to the conventional method, the registration accuracy of the proposed method was significantly improved. This method provided a feasible way for target volume delineation of tumor bed in postoperative radiotherapy of breast cancer patients.

SlicerHeart: An open-source computing platform for cardiac image analysis and modeling

Cardiovascular disease is a significant cause of morbidity and mortality in the developed world. 3D imaging of the heart's structure is critical to the understanding and treatment of cardiovascular disease. However, open-source tools for image analysis of cardiac images, particularly 3D echocardiographic (3DE) data, are limited. We describe the rationale, development, implementation, and application of SlicerHeart, a cardiac-focused toolkit for image analysis built upon 3D Slicer, an open-source image computing platform. We designed and implemented multiple Python scripted modules within 3D Slicer to import, register, and view 3DE data, including new code to volume render and crop 3DE. In addition, we developed dedicated workflows for the modeling and quantitative analysis of multi-modality image-derived heart models, including heart valves. Finally, we created and integrated new functionality to facilitate the planning of cardiac interventions and surgery. We demonstrate application of SlicerHeart to a diverse range of cardiovascular modeling and simulation including volume rendering of 3DE images, mitral valve modeling, transcatheter device modeling, and planning of complex surgical intervention such as cardiac baffle creation. SlicerHeart is an evolving open-source image processing platform based on 3D Slicer initiated to support the investigation and treatment of congenital heart disease. The technology in SlicerHeart provides a robust foundation for 3D image-based investigation in cardiovascular medicine.

Interobserver variability of Gross Tumour Volume delineation for colorectal liver metastases using CT and MRI

Purpose

The purpose of this study was to evaluate the interobserver variability in the contouring of the gross tumor volume (GTV) on magnetic resonance (MR) imaging and computed tomography (CT) for colorectal liver metastases in the setting of SABR.

Methods and Materials

Three expert radiation oncologists contoured 10 GTV volumes on 3 MR imaging sequences and on the CT image data set. Three metrics were chosen to evaluate the interobserver variability: the conformity index, the DICE coefficient, and the maximum Hausdorff distance (HDmax). Statistical analysis of the results was performed using a 1-sided permutation test.

Results

For all 3 metrics, the MR liver acquisition volume acquisition (MR LAVA) showed the lowest interobserver variability. Analysis showed a significant difference (P < .01) in the mean DICE, an overlap metric, for MR LAVA (0.82) and CT (0.74). The HDmax that highlights boundary errors also showed a significant difference (P = .04) with MR LAVA having a lower mean HDmax (7.2 mm) compared with CT (5.7 mm). The mean HDmax for both MR single shot fast spin echo (SSFSE) (19.3 mm) and diffusion weighted image (9.5 mm) showed large interobserver variability with MR SSFSE having a mean HDmax of 19.3 mm. A volume comparison between MR LAVA and CT showed a significantly higher volume for small GTVs (<5 cm³) when using MR LAVA for contouring in comparison to CT.

Conclusions

This study reported the lowest interobserver variability for the MR LAVA, thus indicating the benefit of using MR to complement CT when contouring GTV for colorectal liver metastases.

Validation of an established deep learning auto-segmentation tool for cardiac substructures in 4D radiotherapy planning scans

•

Cardiotoxicity is a common complication of lung cancer radiotherapy.

•

Segmentation of cardiac substructures is time-consuming and challenging.

•

Deep learning segmentation tools can perform this task in 3D and 4D scans.

•

Performance is high when assessed geometrically, dosimetrically and clinically.

•

Auto-segmentation tools may accelerate clinical workflows and enable research.

Background

Emerging data suggest that dose-sparing several key cardiac regions is prognostically beneficial in lung cancer radiotherapy. The cardiac substructures are challenging to contour due to their complex geometry, poor soft tissue definition on computed tomography (CT) and cardiorespiratory motion artefact. A neural network was previously trained to generate the cardiac substructures using three-dimensional radiotherapy planning CT scans (3D-CT). In this study, the performance of that tool on the average intensity projection from four-dimensional (4D) CT scans (4D-AVE), now commonly used in lung radiotherapy, was evaluated.

Materials and Methods

The 4D-AVE of n=20 patients completing radiotherapy for lung cancer 2015–2020 underwent manual and automated cardiac substructure segmentation. Manual and automated substructures were compared geometrically and dosimetrically. Two senior clinicians also qualitatively assessed the auto-segmentation tool’s output.

Results

Geometric comparison of the automated and manual segmentations exhibited high levels of similarity across parameters, including volume difference (11.8% overall) and Dice similarity coefficient (0.85 overall), and were consistent with 3D-CT performance. Differences in mean (median 0.2 Gy, range −1.6–0.3 Gy) and maximum (median 0.4 Gy, range −2.2–0.9 Gy) doses to substructures were generally small. Nearly all structures (99.5 %) were deemed to be appropriate for clinical use without further editing.

Conclusions

Cardiac substructure auto-segmentation using a deep learning-based tool trained on a 3D-CT dataset was feasible on the 4D-AVE scan, meaning this tool is suitable for use on 4D-CT radiotherapy planning scans. Application of this tool would increase the practicality of routine clinical cardiac substructure delineation, and enable further cardiac radiation effects research.

Pneumonitis after Stereotactic Thoracic Radioimmunotherapy with Checkpoint Inhibitors: Exploration of the Dose-Volume-Effect Correlation

Simple Summary

Stereotactic body radiation therapy (SBRT) is widely applied for treatment of early stage lung cancer and pulmonary metastases. Modern immune checkpoint blockade (ICB) is progressively used in cancer treatment. Pneumonitis is a relevant side effect of both thoracic SBRT and ICB. Currently, it remains unclear whether we can presume the same radiation dose–volume–effect correlations and dose constraints for safe application of SBRT + ICB. We present a dose–volume–effect correlation analysis method using pneumonitis contours and dose–volume histograms (DVH). We showed dosimetric differences for pneumonitis volumes between SBRT + ICB and SBRT alone. We found a large extent of pneumonitis, even bilateral and apart from the radiation field for combined SBRT + ICB. We noticed a shift in pneumonitis DVHs towards lower doses and a trend towards decreased areas under the curve (AUC) for SBRT + ICB. This provides a direction for re-evaluation and potential adaptation of lung dose constraints for combined SBRT and ICB.

Abstract

Thoracic stereotactic body radiation therapy (SBRT) is extensively used in combination with immune checkpoint blockade (ICB). While current evidence suggests that the occurrence of pneumonitis as a side effect of both treatments is not enhanced for the combination, the dose–volume correlation remains unclear. We investigate dose–volume–effect correlations for pneumonitis after combined SBRT + ICB. We analyzed patient clinical characteristics and dosimetric data for 42 data sets for thoracic SBRT with ICB treatment (13) and without (29). Dose volumes were converted into 2 Gy equivalent doses (EQD2), allowing for dosimetric comparison of different fractionation regimes. Pneumonitis volumes were delineated and corresponding DVHs were analyzed. We noticed a shift towards lower doses for combined SBRT + ICB treatment, supported by a trend of smaller areas under the curve (AUC) for SBRT+ ICB (median AUC 1337.37 vs. 5799.10, p = 0.317). We present a DVH-based dose–volume–effect correlation method and observed large pneumonitis volumes, even with bilateral extent in the SBRT + ICB group. We conclude that further studies using this method with enhanced statistical power are needed to clarify whether adjustments of the radiation dose constraints are required to better estimate risks of pneumonitis after the combination of SBRT and ICB.

The dose accumulation and the impact of deformable image registration on dose reporting parameters in a moving patient undergoing proton radiotherapy

INTRODUCTION

Potential changes in patient anatomy during proton radiotherapy may lead to a deviation of the delivered dose. A dose estimate can be computed through a deformable image registration (DIR) driven dose accumulation. The present study evaluates the accumulated dose uncertainties in a patient subject to an inadvertent breathing associated motion.

MATERIALS AND METHODS

A virtual lung tumour was inserted into a pair of single participant landmark annotated computed tomography images depicting opposite breathing phases, with the deep inspiration breath-hold the planning reference and the exhale the off-reference geometry. A novel Monte Carlo N-Particle, Version 6 (MCNP6) dose engine was developed, validated and used in treatment plan optimization. Three DIR methods were compared and used to transfer the exhale simulated dose to the reference geometry. Dose conformity and homogeneity measures from International Committee on Radioactivity Units and Measurements (ICRU) reports 78 and 83 were evaluated on simulated dose distributions registered with different DIR algorithms.

RESULTS

The MCNP6 dose engine handled patient-like geometries in reasonable dose calculation times. All registration methods were able to align image associated landmarks to distances, comparable to voxel sizes. A moderate deterioration of ICRU measures was encountered in comparing doses in on and off-reference anatomy. There were statistically significant DIR driven differences in ICRU measures, particularly a 10% difference in the relative D_98% for planning tumour volume and in the 3 mm/3% gamma passing rate.

CONCLUSIONS

T he dose accumulation over two anatomies resulted in a DIR driven uncertainty, important in reporting the associated ICRU measures for quality assurance.

Unsupervised pseudo CT generation using heterogenous multicentric CT/MR images and CycleGAN: Dosimetric assessment for 3D conformal radiotherapy

PURPOSE

Absorbed dose calculation in magnetic resonance-guided radiation therapy (MRgRT) is commonly based on pseudo CT (pCT) images. This study investigated the feasibility of unsupervised pCT generation from MRI using a cycle generative adversarial network (CycleGAN) and a heterogenous multicentric dataset. A dosimetric analysis in three-dimensional conformal radiotherapy (3DCRT) planning was also performed.

MATERIAL AND METHODS

Overall, 87 T1-weighted and 102 T2-weighted MR images alongside with their corresponding computed tomography (CT) images of brain cancer patients from multiple centers were used. Initially, images underwent a number of preprocessing steps, including rigid registration, novel CT Masker, N4 bias field correction, resampling, resizing, and rescaling. To overcome the gradient vanishing problem, residual blocks and mean squared error (MSE) loss function were utilized in the generator and in both networks (generator and discriminator), respectively. The CycleGAN was trained and validated using 70 T1 and 80 T2 randomly selected patients in an unsupervised manner. The remaining patients were used as a holdout test set to report final evaluation metrics. The generated pCTs were validated in the context of 3DCRT.

RESULTS

The CycleGAN model using masked T2 images achieved better performance with a mean absolute error (MAE) of 61.87 ± 22.58 HU, peak signal to noise ratio (PSNR) of 27.05 ± 2.25 (dB), and structural similarity index metric (SSIM) of 0.84 ± 0.05 on the test dataset. T1-weighted MR images used for dosimetric assessment revealed a gamma index of 3%, 3 mm, 2%, 2 mm and 1%, 1 mm with acceptance criteria of 98.96% ± 1.1%, 95% ± 3.68%, 90.1% ± 6.05%, respectively. The DVH differences between CTs and pCTs were within 2%.

CONCLUSIONS

A promising pCT generation model capable of handling heterogenous multicenteric datasets was proposed. All MR sequences performed competitively with no significant difference in pCT generation. The proposed CT Masker proved promising in improving the model accuracy and robustness. There was no significant difference between using T1-weighted and T2-weighted MR images for pCT generation.

Experimental validation of 4D log file-based proton dose reconstruction for interplay assessment considering amplitude-sorted 4DCTs

PURPOSE

The unpredictable interplay between dynamic proton therapy delivery and target motion in the thorax can lead to severe dose distortions. A fraction-wise four-dimensional (4D) dose reconstruction workflow allows for the assessment of the applied dose after patient treatment while considering the actual beam delivery sequence extracted from machine log files, the recorded breathing pattern and the geometric information from a 4D computed tomography scan (4DCT). Such an algorithm capable of accounting for amplitude-sorted 4DCTs was implemented and its accuracy as well as its sensitivity to input parameter variations was experimentally evaluated.

METHODS

An anthropomorphic thorax phantom with a movable insert containing a target surrogate and a radiochromic film was irradiated with a monoenergetic field for various 1D target motion forms (sin, sin4) and peak-to-peak amplitudes (5/10/15/20/30 mm). The measured characteristic film dose distributions were compared to the respective sections in the 4D reconstructed doses using a 2D γ-analysis (3mm, 3%); γ-pass rates were derived for different dose grid resolutions (1mm/3mm) and deformable image registrations (DIR, automatic/manual) applied during the 4D dose reconstruction process. In an additional analysis, the sensitivity of reconstructed dose distributions against potential asynchronous timing of the motion and machine log files was investigated for both a monoenergetic field and more realistic 4D robustly optimized fields by artificially introduced offsets of ± 1/5/25/50/250 ms. The resulting dose distributions with asynchronized log files were compared to those with synchronized log files by means of a 3D γ-analysis (1mm, 1%) and the evaluation of absolute dose differences.

RESULTS

The induced characteristic interplay patterns on the films were well reproduced by the 4D dose reconstruction with 2D γ-pass rates ≥95% for almost all cases with motion magnitudes ≤15 mm. In general, the 2D γ-pass rates showed a significant decrease for larger motion amplitudes and increase when using a finer dose grid resolution but were not affected by the choice of motion form (sin, sin4). There was also a trend, though not statistically significant, towards the manually defined DIR for better quality of the reconstructed dose distributions in the area imaged by the film. The 4D dose reconstruction results for the monoenergetic as well as the 4D robustly optimized fields were robust against small asynchronies between motion and machine log files of up to 5 ms, which is in the order of potential network latencies.

CONCLUSIONS

We have implemented a 4D log file-based proton dose reconstruction that accounts for amplitude-sorted 4DCTs. Its accuracy was proven to be clinically acceptable for target motion magnitudes of up to 15 mm. Particular attention should be paid to the synchronization of the log file generating systems as the reconstructed dose distribution may vary with log file asynchronies larger than those caused by realistic network delays. This article is protected by copyright. All rights reserved.

Dosimetric assessment of patient dose calculation on a deep learning-based synthesized computed tomography image for adaptive radiotherapy

Purpose

Dose computation using cone beam computed tomography (CBCT) images is inaccurate for the purpose of adaptive treatment planning. The main goal of this study is to assess the dosimetric accuracy of synthetic computed tomography (CT)‐based calculation for adaptive planning in the upper abdominal region. We hypothesized that deep learning‐based synthetically generated CT images will produce comparable results to a deformed CT (CTdef) in terms of dose calculation, while displaying a more accurate representation of the daily anatomy and therefore superior dosimetric accuracy.

Methods

We have implemented a cycle‐consistent generative adversarial networks (CycleGANs) architecture to synthesize CT images from the daily acquired CBCT image with minimal error. CBCT and CT images from 17 liver stereotactic body radiation therapy (SBRT) patients were used to train, test, and validate the algorithm.

Results

The synthetically generated images showed increased signal‐to‐noise ratio, contrast resolution, and reduced root mean square error, mean absolute error, noise, and artifact severity. Superior edge matching, sharpness, and preservation of anatomical structures from the CBCT images were observed for the synthetic images when compared to the CTdef registration method. Three verification plans (CBCT, CTdef, and synthetic) were created from the original treatment plan and dose volume histogram (DVH) statistics were calculated. The synthetic‐based calculation shows comparatively similar results to the CTdef‐based calculation with a maximum mean deviation of 1.5%.

Conclusions

Our findings show that CycleGANs can produce reliable synthetic images for the adaptive delivery framework. Dose calculations can be performed on synthetic images with minimal error. Additionally, enhanced image quality should translate into better daily alignment, increasing treatment delivery accuracy.

Proof-of-concept of DosiTest: A virtual multicentric clinical trial for assessing uncertainties in molecular radiotherapy dosimetry

Clinical dosimetry in molecular radiotherapy (MRT) is a multi-step procedure, prone to uncertainties at every stage of the dosimetric workflow. These are difficult to assess, especially as some are complex or even impossible to measure experimentally. The DosiTest project was initiated to assess the variability associated with clinical dosimetry, by setting up a 'virtual' multicentric clinical dosimetry trial based on Monte Carlo (MC) modelling. A reference patient model with a realistic geometry and activity input for a specific tracer is considered. Reference absorbed dose rate distribution maps are generated at various time-points from MC modelling, combining precise information on density and activity distributions (voxel wise). Then, centre-specific calibration and patient SPECT/CT datasets are modelled, on which the clinical centres can perform clinical (i.e. image-based) dosimetry. The results of this dosimetric analysis can be benchmarked against the reference dosimetry to assess the variability induced by implementing different clinical dosimetry approaches. The feasibility of DosiTest is presented here for a clinical situation of therapeutic administration of ¹⁷⁷Lu-DOTATATE (Lutathera®) peptide receptor radionuclide therapy (PRRT). From a real patient dataset composed of 5 SPECT/CT images and associated calibrations, we generated the reference absorbed dose rate images with GATE. Then, simulated SPECT/CT image generation based on GATE was performed, both for a calibration phantom and virtual patient images. Based on this simulated dataset, image-based dosimetry could be performed, and compared with reference dosimetry. The good agreement, between real and simulated images, and between reference and image-based dosimetry established the proof of concept of DosiTest.

Transient Enlargement in Meningiomas Treated with Stereotactic Radiotherapy

Simple Summary

Accurate assessment of treatment efficacy is a prerequisite for the improvement in therapeutic outcomes in clinical trials. However, it is very challenging to accurately track the size of meningiomas after radiotherapy, because of their complex shapes and often slow growth. Measuring the whole tumor volume as opposed to simple diameter measurements to assess treatment efficacy, therefore, is very promising but little is known on expected volumetric changes of meningiomas following radiotherapy. Therefore, in this study, we meticulously investigated volumetric changes in meningiomas following radiotherapy incorporating volumetric measurements from 468 MRI studies and evaluated newly proposed RANO volumetric criteria in the context of radiotherapy. We found that temporary tumor enlargement after radiotherapy overall was rare but occurred significantly more frequently after high than after low single doses of radiation, represented an important differential diagnosis to tumor progression and would have skewed results in a clinical trial if not accounted for.

Abstract

To investigate the occurrence of pseudoprogression/transient enlargement in meningiomas after stereotactic radiotherapy (RT) and to evaluate recently proposed volumetric RANO meningioma criteria for response assessment in the context of RT. Sixty-nine meningiomas (benign: 90%, atypical: 10%) received stereotactic RT from January 2005–May 2018. A total of 468 MRI studies were segmented longitudinally during a median follow-up of 42.3 months. Best response and local control were evaluated according to recently proposed volumetric RANO criteria. Transient enlargement was defined as volumetric increase ≥20% followed by a subsequent regression ≥20%. The mean best volumetric response was −23% change from baseline (range, −86% to +19%). According to RANO, the best volumetric response was SD in 81% (56/69), MR in 13% (9/69) and PR in 6% (4/69). Transient enlargement occurred in only 6% (4/69) post RT but would have represented 60% (3/5) of cases with progressive disease if not accounted for. Transient enlargement was characterized by a mean maximum volumetric increase of +181% (range, +24% to +389 %) with all cases occurring in the first year post-RT (range, 4.1–10.3 months). Transient enlargement was significantly more frequent with SRS or hypofractionation than with conventional fractionation (25% vs. 2%, p = 0.015). Five-year volumetric control was 97.8% if transient enlargement was recognized but 92.9% if not accounted for. Transient enlargement/pseudoprogression in the first year following SRS and hypofractionated RT represents an important differential diagnosis, especially because of the high volumetric control achieved with stereotactic RT. Meningioma enlargement during subsequent post-RT follow-up and after conventional fractionation should raise suspicion for tumor progression.

A Subset of Intracanalicular Vestibular Schwannomas Demonstrates Minimal Growth Over a 10-Year Period

OBJECTIVE

Vestibular schwannomas (VS) commonly undergo magnetic resonance imaging (MRI) surveillance, but long-term data to support the ideal frequency is limited. Herein, we aim to investigate intracanalicular VS growth predictors and long-term growth rates (GR).

STUDY DESIGN

Retrospective chart review.

SETTING

Two tertiary care centers.

PATIENTS

Sporadic intracanalicular VS with initial conservative management and at least two sequential MRIs.

INTERVENTION

Serial MRI.

MAIN OUTCOME MEASURES

VS were categorized by baseline internal auditory canal tertile sublocalization (fundus, midpoint, porus) and size (≤100, 100-200, >200 mm3). Throughout follow-up, volumetric GR (mm3/yr) were determined (baseline-3 yrs, 3-5 yrs, 5-10 yrs) and treatment rates were assessed.

RESULTS

Ninety-nine intracanalicular VS were identified (mean follow-up of 6.1 ± 4.5 yrs). Mean GR before 5-year follow-up were comparable for baseline tertile involvement and size. After 5-year follow-up, mean GR of VS involving the fundus at baseline were lower than those involving the midpoint and fundus (6.17 ± 21.16 and 119.74 ± 117.57 mm3/yr, respectively; p = 0.034). Mean GR of VS with less than or equal to 100 mm3 at baseline (-7.29 ± 25.44 mm3/yr) were lower than those with 100 to 200 mm3 (86.55 ± 103.99 mm3/yr; p = 0.011) and more than 200 mm3 (45.70 ± 35.71 mm3/yr; p = 0.031). Vestibular schwannomas involving the midpoint and fundus had greater treatment rates compared with VS involving only the fundus (p < 0.001).

CONCLUSIONS

Baseline tertile involvement and size may predict long-term intracanalicular VS growth where fundal tumors or those less than or equal to 100 mm3 exhibit little long-term growth. Extending surveillance after 5-year follow-up may be reasonable for fundal VS.

Treatment Planning Study for Microbeam Radiotherapy Using Clinical Patient Data

Microbeam radiotherapy (MRT) is a novel, still preclinical dose delivery technique. MRT has shown reduced normal tissue effects at equal tumor control rates compared to conventional radiotherapy. Treatment planning studies are required to permit clinical application. The aim of this study was to establish a dose comparison between MRT and conventional radiotherapy and to identify suitable clinical scenarios for future applications of MRT. We simulated MRT treatment scenarios for clinical patient data using an inhouse developed planning algorithm based on a hybrid Monte Carlo dose calculation and implemented the concept of equivalent uniform dose (EUD) for MRT dose evaluation. The investigated clinical scenarios comprised fractionated radiotherapy of a glioblastoma resection cavity, a lung stereotactic body radiotherapy (SBRT), palliative bone metastasis irradiation, brain metastasis radiosurgery and hypofractionated breast cancer radiotherapy. Clinically acceptable treatment plans were achieved for most analyzed parameters. Lung SBRT seemed the most challenging treatment scenario. Major limitations comprised treatment plan optimization and dose calculation considering the tissue microstructure. This study presents an important step of the development towards clinical MRT. For clinical treatment scenarios using a sophisticated dose comparison concept based on EUD and EQD2, we demonstrated the capability of MRT to achieve clinically acceptable dose distributions.

Contouring cardiac substructures on average intensity projection 4D-CT for lung cancer radiotherapy: a proposal of a heart valve contouring atlas

BACKGROUND AND PURPOSE

A detailed contouring atlas of the heart valves is lacking. Existing heart contouring atlases have not been evaluated on average intensity projection four-dimensional non-contrast computed tomography (AVE-4D-CT) scans, routinely used for organ-at-risk delineation in lung cancer radiotherapy. We aimed to develop the heart valve contouring atlas and to assess inter-observer variation in delineation of the heart, its substructures, and coronary arteries on AVE-4D-CT scans, along with its impact on radiotherapy doses.

MATERIALS AND METHODS

A heart valve contouring atlas was developed. Five radiation oncologists and four cardiologists delineated the valves according to this atlas, and the remaining heart substructures according to the existing atlases, on AVE-4D-CT scans of ten patients who underwent radio(chemo)therapy for NSCLC. The observer contours were then compared to the collectively defined "reference" contours. Spatial variation was assessed using the Sørensen-Dice similarity coefficient (DSC), directed average Hausdorff distance (DAH), directed Hausdorff distance (HD), and the mean distance to agreement (MDA). The effect of spatial variation on radiotherapy doses was assessed using the patients' treatment plans.

RESULTS

Inter-observer contour overlap (mean DSC) was 0.68, 0.49, 0.45 and 0.45, and inter-observer contour separation (mean DAH) was 2.1, 3.4, 2.6 and 2.9 mm for the pulmonic (PV), aortic (AV), mitral (MV) and tricuspid valve (TV), respectively. Mean HD was higher for TV and MV (13.3 and 11.7mm) than for AV and PV (7.8 and 7mm). The highest mean MDA of 3.1mm was found for AV, and the lowest (1.9mm) for PV. Inter-observer agreement was the lowest for the coronary arteries, but statistically significant dose variation was found mainly in the left ventricular septal and anterior segments.

CONCLUSION

Our atlas enables reproducible delineation of the heart valves. Delineation of the heart and its substructures on AVE-4D-CT scans is feasible, with inter-observer variability similar to that reported on conventional non-contrast CT scans.

Oesophageal IGRT considerations for SBRT of LA-NSCLC: barium-enhanced CBCT and interfraction motion

Background

To determine the optimal volume of barium for oesophageal localisation on cone-beam CT (CBCT) for locally-advanced non-small cell lung cancers (NSCLC) and quantify the interfraction oesophageal movement relative to tumour.

Methods

Twenty NSCLC patients with mediastinal and/or hilar disease receiving radical radiotherapy were recruited. The first five patients received 25 ml of barium prior to their planning CT and alternate CBCTs during treatment. Subsequent five patient cohorts, received 15 ml, 10 ml and 5 ml. Six observers contoured the oesophagus on each of the 107 datasets and consensus contours were created. Overall 642 observer contours were generated and interobserver contouring reproducibility was assessed. The kappa statistic, dice coefficient and Hausdorff Distance (HD) were used to compare barium-enhanced CBCTs and non-enhanced CBCTs. Oesophageal displacement was assessed using the HD between consensus contours of barium-enhanced CBCTs and planning CTs.

Results

Interobserver contouring reproducibility was significantly improved in barium-enhanced CBCTs compared to non-contrast CBCTs with minimal difference between barium dose levels. Only 10 mL produced a significantly higher kappa (0.814, p = 0.008) and dice (0.895, p = 0.001). The poorer the reproducibility without barium, the greater the improvement barium provided. The median interfraction HD between consensus contours was 4 mm, with 95% of the oesophageal displacement within 15 mm.

Conclusions

10 mL of barium significantly improves oesophageal localisation on CBCT with minimal image artifact. The oesophagus moves substantially and unpredictably over a course of treatment, requiring close daily monitoring in the context of hypofractionation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-021-01946-8.

Profeel-An open source dosimetry data visualization and analysis software

BACKGROUND AND OBJECTIVES

This article presents Profeel, a Matlab (MathWorks Inc., MA) based open source dosimetry data visualization and analysis software. Profeel aims to support quality assurance, dosimetry and research in the field of radiotherapy by providing an environment to visualize, process and analyse measured and simulated dosimetry data from several data sources used in radiotherapy practice and research.

METHODS

The processing and analysis tools are based on routinely used dosimetry analysis methods, such as gamma analysis, different data normalizations and data filtering. Additionally the Profeel performs an automatic 1 dimensional profile and percentage depth dose analysis in accordance with International Electrotechnical Commission definitions. All data can be operated by user created custom functions and lower dimensionality data can be extracted from volume doses and dose planes.

RESULTS

Profeel supports data import in all 3 dimensions and offers an intuitive user interface to perform data visualization, processing and analysis between simulated and measured data. Profeel and its source code are distributed free of charge under the General Public Licence (GPL).

CONCLUSIONS

Profeel has shown to be an agile tool for fulfilling various needs of several researchers and since Profeel is under constant development and is an open source project, community needs, issues and bug reports are taken into account in the development.

Segmentation of vestibular schwannoma from MRI, an open annotated dataset and baseline algorithm

Automatic segmentation of vestibular schwannomas (VS) from magnetic resonance imaging (MRI) could significantly improve clinical workflow and assist patient management. We have previously developed a novel artificial intelligence framework based on a 2.5D convolutional neural network achieving excellent results equivalent to those achieved by an independent human annotator. Here, we provide the first publicly-available annotated imaging dataset of VS by releasing the data and annotations used in our prior work. This collection contains a labelled dataset of 484 MR images collected on 242 consecutive patients with a VS undergoing Gamma Knife Stereotactic Radiosurgery at a single institution. Data includes all segmentations and contours used in treatment planning and details of the administered dose. Implementation of our automated segmentation algorithm uses MONAI, a freely-available open-source framework for deep learning in healthcare imaging. These data will facilitate the development and validation of automated segmentation frameworks for VS and may also be used to develop other multi-modal algorithmic models.

Development and Validation of an Efficient MRI Radiomics Signature for Improving the Predictive Performance of 1p/19q Co-Deletion in Lower-Grade Gliomas

Simple Summary

Low-grade gliomas (LGG) with the 1p/19q co-deletion mutation have been proven to have a better survival prognosis and response to treatment than individuals without the mutation. Identifying this mutation has a vital role in managing LGG patients; however, the current diagnostic gold standard, including the brain-tissue biopsy or the surgical resection of the tumor, remains highly invasive and time-consuming. We proposed a model based on the eXtreme Gradient Boosting (XGBoost) classifier to detect 1p/19q co-deletion mutation using non-invasive medical images. The performance of our model achieved 87% and 82.8% accuracy on the training and external test set, respectively. Significantly, the prediction was based on only seven optimal wavelet radiomics features extracted from brain Magnetic Resonance (MR) images. We believe that this model can address clinicians in the rapid diagnosis of clinical 1p/19q co-deletion mutation, thereby improving the treatment prognosis of LGG patients.

Abstract

The prognosis and treatment plans for patients diagnosed with low-grade gliomas (LGGs) may significantly be improved if there is evidence of chromosome 1p/19q co-deletion mutation. Many studies proved that the codeletion status of 1p/19q enhances the sensitivity of the tumor to different types of therapeutics. However, the current clinical gold standard of detecting this chromosomal mutation remains invasive and poses implicit risks to patients. Radiomics features derived from medical images have been used as a new approach for non-invasive diagnosis and clinical decisions. This study proposed an eXtreme Gradient Boosting (XGBoost)-based model to predict the 1p/19q codeletion status in a binary classification task. We trained our model on the public database extracted from The Cancer Imaging Archive (TCIA), including 159 LGG patients with 1p/19q co-deletion mutation status. The XGBoost was the baseline algorithm, and we combined the SHapley Additive exPlanations (SHAP) analysis to select the seven most optimal radiomics features to build the final predictive model. Our final model achieved an accuracy of 87% and 82.8% on the training set and external test set, respectively. With seven wavelet radiomics features, our XGBoost-based model can identify the 1p/19q codeletion status in LGG-diagnosed patients for better management and address the drawbacks of invasive gold-standard tests in clinical practice.

Clinically Oriented Target Contour Evaluation Using Geometric and Dosimetric Indices Based on Simple Geometric Transformations

Purpose:

In radiotherapy, geometric indices are often used to evaluate the accuracy of contouring. However, the ability of geometric indices to identify the error of contouring results is limited primarily because they do not consider the clinical background. The purpose of this study is to investigate the relationship between geometric and clinical dosimetric indices.

Methods:

Four different types of targets were selected (C-shaped target, oropharyngeal cancer, metastatic spine cancer, and prostate cancer), and the translation, scaling, rotation, and sine function transformation were performed with the software Python to introduce systematic and random errors. The transformed contours were regarded as reference contours. Dosimetric indices were obtained from the original dose distribution of the radiotherapy plan. The correlations between geometric and dosimetric indices were quantified by linear regression.

Results:

The correlations between the geometric and dosimetric indices were inconsistent. For systematic errors, and with the exception of the sine function transformation (R²: 0.023-0.04, P > 0.05), the geometric transformations of the C-shaped target were correlated with the D98% and D_mean (R²: 0.689-0.988), 80% of which were P < 0.001. For the random errors, the correlations obtained by the all targets were R² > 0.384, P < 0.05. The Wilcoxon signed-rank test was used to compare the spatial direction resolution capability of geometric indices in different directions of the C-shaped target (with systematic errors), and the results showed only the volumetric geometric indices with P < 0.05.

Conclusions:

Clinically, an assessment of the contour accuracy of the region-of-interest is not feasible based on geometric indices alone. Dosimetric indices should be added to the evaluations of the accuracy of the delineation results, which can be helpful for explaining the clinical dose response relationship of delineation more comprehensively and accurately.

Deformable image registration to assist clinical decision for radiotherapy treatment adaptation for head and neck cancer patients

Head and neck (H&N) cancer patients often present anatomical and geometrical changes in tumors and organs at risk (OARs) during radiotherapy treatment. These changes may result in the need to adapt the existing treatment planning, using an expert's subjective opinion, for offline adaptive radiotherapy and a new treatment planning before each treatment, for online adaptive radiotherapy. In the present study, a fast methodology is proposed to assist in planning adaptation clinical decision using tumor and parotid glands percentage volume changes during treatment. The proposed approach was applied to 40 Η&Ν cases, with one planning Computed Tomography (pCT) image and CBCT scans for 6 weeks of treatment per case. Deformable registration was used for each patient's pCT image alignment to its weekly CBCT. The calculated transformations were used to align each patient's anatomical structures to the weekly anatomy. Clinical target volume (CTV) and parotid gland volume percentage changes were calculated in each case. The accuracy of the achieved image alignment was validated qualitatively and quantitatively. Furthermore, statistical analysis was performed to test if there is a statistically significant correlation between CTV and parotid glands volume percentage changes. Average MDA for CTV and parotid glands between corresponding structures defined by an expert in CBCTs and automatically calculated through registration was 1.4 ± 0.1 mm and 1.5 ± 0.1 mm, respectively. The mean registration time of the first CBCT image registration for 40 cases was lower than 3.4 min. Five patients show more than 20% tumor volume change. Six patients show more than 30% parotid glands volume change. Ten out of 40 patients proposed for planning adaptation. All the statistical tests performed showed no correlation between CTV/parotid glands percentage volume changes. The aim to assist in clinical decision making on a fast and automatic way was achieved using the proposed methodology, thereby reducing workload in clinical practice.

Estimation of delivered dose to lung tumours considering setup uncertainties and breathing motion in a cohort of patients treated with stereotactic body radiation therapy

INTRODUCTION

Dose-response relationships for local control of lung tumours treated with stereotactic body radiotherapy (SBRT) have proved ambiguous, however, these have been based on the prescribed or planned dose. Delivered dose to the target may be a better predictor for local control. In this study, the probability of the delivered minimum dose to the clinical target volume (CTV) in relation to the prescribed dose was estimated for a cohort of patients, considering geometrical uncertainties.

MATERIALS AND METHODS

Delivered doses were retrospectively simulated for 50 patients treated with SBRT for lung tumours, comparing two image-guidance techniques: pre-treatment verification computed tomography (IG1) and online cone-beam computed tomography (IG2). The prescribed dose was typically to the 67% isodose line of the treatment plan. Simulations used in-house software that shifted the static planned dose according to a breathing motion and sampled setup/matching errors. Each treatment was repeatedly simulated, generating a multiplicity of dose-volume histograms (DVH). From these, tumour-specific and population-averaged statistics were derived.

RESULTS

For IG1, the probability that the minimum CTV dose (D_98%) exceeded 100% of the prescribed dose was 90%. With IG2, this probability increased to 99%.

CONCLUSIONS

Doses below the prescribed dose were delivered to a considerably larger part of the population prior to the introduction of online soft-tissue image-guidance. However, there is no clear evidence that this impacts local control, when compared to previous published data.