Isodose surface differences: A novel tool for the comparison of dose distributions.

BACKGROUND: Comparing dose distributions is a routine task in radiotherapy, mainly in patient-specific quality assurance (PSQA). Currently, the evaluation of the dose distributions is being performed mainly with statistical methods, which could underestimate the clinical importance of the spotted differences, as per the literature. PURPOSE: This paper aims to provide proof-of-concept for a novel dose distribution comparison method based on the difference of the isodose surfaces. The new method connects acceptance tolerance to QA limitations (equipment capabilities) and integrates a clinical approach into the analysis procedure. METHODS: The distance of dose points from the isocenter can be used as a function to define the shape of an isodose surface expressed as a histogram. Isodose surface differences (ISD) are defined as the normalized differences of reference and evaluated surface histograms plotted against their corresponding isodose. Acceptance tolerances originate from actual QA tolerances and are presented clinically intuitively. The ISD method was compared to the gamma index using intentionally erroneous VMAT and IMRT plans. RESULTS: Results revealed that the ISD method is sensitive to all errors induced in the plans. Discrepancies are presented per isodose, enabling the evaluation of the plan in two regions representing PTV and Normal Tissue. ISD manages to flag errors that would remain undetected under the gamma analysis. CONCLUSION: The ISD method is a meaningful, QA-related, registration-free, and clinically oriented technique of dose distribution evaluation. This method can be used either as a standalone or an auxiliary tool to the well-established evaluation procedures, overcoming significant limitations reported in the literature.

Dosimetric accuracy of the Convolution algorithm for Leksell Gamma Plan radiosurgery treatment planning: Evaluation in the presence of clinically relevant inhomogeneities.

PURPOSE: The Leksell Gamma Plan Convolution algorithm (LGP-Convolution) has not been widely adopted. This mainly stems from the higher calculated beam-on times relative to the standard ray tracing-based LGP-TMR10 dose calculation algorithm. This study aims to evaluate the accuracy of the LGP-Convolution in scenarios where the treated lesions are in the vicinity of or encompassed by bone and/or air inhomogeneities. METHODS: The solid water dosimetry phantom provided by the vendor was modified to include bone and air inhomogeneities. Two treatment planning scenarios were investigated involving a single shot and multiple shots, respectively. Treatment planning and dose prescription were performed using the LGP-Convolution algorithm. Triple channel film dosimetry was performed using GafChromic EBT3 films calibrated in terms of absorbed dose to water in a 60 Co beam. Monte Carlo (MC) simulation dosimetry was also performed in the inhomogeneous experimental geometry using the EGSnrc MC platform and a previously validated sector-based phase-space source model. MC simulations were also employed to determine correction factors required for converting EBT3 measurements at points within the bone and air inhomogeneities from dose-to-water values to the corresponding dose to medium values. RESULTS AND CONCLUSIONS: EBT3 dose to medium correction factors ranged with field size (4, 8, or 16 mm) within 0.941-0.946 for bone and 0.745-0.749 for air inhomogeneities. An excellent agreement was found between the LGP-Convolution calculations with corresponding EBT3 and MC dose to medium results at all measurement points, except those located inside the air inhomogeneity. The latter is of no clinical importance and excluding them yielded gamma index passing rates of nearly 100% for 3% local dose difference and 1 mm distance-to-agreement criteria. The excellent agreement observed between LGP-Convolution calculations and film as well as MC results of dose to medium indicates that the latter is the quantity reported by the LGP-Convolution.

Target localization accuracy in frame-based stereotactic radiosurgery: Comparison between MR-only and MR/CT co-registration approaches.

PURPOSE: In frame-based Gamma Knife (GK) stereotactic radiosurgery two treatment planning workflows are commonly employed; one based solely on magnetic resonance (MR) images and the other based on magnetic resonance/computed tomography (MR/CT) co-registered images. In both workflows, target localization accuracy (TLA) can be deteriorated due to MR-related geometric distortions and/or MR/CT co-registration uncertainties. In this study, the overall TLA following both clinical workflows is evaluated for cases of multiple brain metastases. METHODS: A polymer gel-filled head phantom, having the Leksell stereotactic headframe attached, was CT-imaged and irradiated by a GK Perfexion unit. A total of 26 4-mm shots were delivered at 26 locations directly defined in the Leksell stereotactic space (LSS), inducing adequate contrast in corresponding T2-weighted (T2w) MR images. Prescribed shot coordinates served as reference locations. An additional MR scan was acquired to implement the "mean image" distortion correction technique. The TLA for each workflow was assessed by comparing the radiation-induced target locations, identified in MR images, with corresponding reference locations. Using T1w MR and CT images of 15 patients (totaling 81 lesions), TLA in clinical cases was similarly assessed, considering MR-corrected data as reference. For the MR/CT workflow, both global and region of interest (ROI)-based MR/CT registration approaches were studied. RESULTS: In phantom measurements, the MR-corrected workflow demonstrated unsurpassed TLA (median offset of 0.2 mm) which deteriorated for MR-only and MR/CT workflows (median offsets of 0.8 and 0.6 mm, respectively). In real-patient cases, the MR-only workflow resulted in offsets that exhibit a significant positive correlation with the distance from the MR isocenter, reaching 1.1 mm (median 0.6 mm). Comparable results were obtained for the MR/CT-global workflow, although a maximum offset of 1.4 mm was detected. TLA was improved with the MR/CT-ROI workflow resulting in median/maximum offsets of 0.4 mm/1.1 mm. CONCLUSIONS: Subpixel TLA is achievable in all workflows. For the MR/CT workflow, a ROI-based MR/CT co-registration approach could considerably increase TLA and should be preferred instead of a global registration.

Impact of systematic MLC positional uncertainties on the quality of single-isocenter multi-target VMAT-SRS treatment plans.

PURPOSE: To study the impact of systematic MLC leaf positional uncertainties (stemming from mechanical inaccuracies or sub-optimal MLC modeling) on the quality of intracranial single-isocenter multi-target VMAT-SRS treatment plans. An estimation of appropriate tolerance levels is attempted. METHODS: Five patients, with three to four metastases and at least one target lying in close proximity to organs-at-risk (OARs) were included in this study. A single-isocenter multi-arc VMAT plan per patient was prepared, which served as the reference for dosimetric impact evaluation. A range of leaf offsets was introduced (±0.03 mm up to ±0.30 mm defined at the MLC plane) to both leaf banks, by varying the leaf offset MLC modeling parameter in Monaco for all the prepared plans, in order to simulate projected leaf offsets of ±0.09 mm up to ±0.94 mm at the isocenter plane, respectively. For all offsets simulated and cases studied, dose distributions were re-calculated and compared with the corresponding reference ones. An experimental dosimetric procedure using the SRS mapCHECK diode array was also performed to support the simulation study results and investigate its suitability to detect small systematic leaf positional errors. RESULTS: Projected leaf offsets of ±0.09 mm were well-tolerated with respect to both target dosimetry and OAR-sparing. A linear relationship was found between D95% percentage change and projected leaf offset (slope: 12%/mm). Impact of projected offset on target dosimetry was strongly associated with target volume. In two cases, plans that could be considered potentially clinically unacceptable (i.e., clinical dose constraint violation) were obtained even for projected offsets as small as 0.19 mm. The performed experimental dosimetry check can detect potential small systematic leaf errors. CONCLUSIONS: Plan quality indices and dose-volume metrics are very sensitive to systematic sub-millimeter leaf positional inaccuracies, projected at the isocenter plane. Acceptable and tolerance levels in systematic MLC uncertainties need to be tailored to VMAT-SRS spatial and dosimetric accuracy requirements.

Dosimetric impact of rotational errors on the quality of VMAT-SRS for multiple brain metastases: Comparison between single- and two-isocenter treatment planning techniques.

PURPOSE: In the absence of a 6D couch and/or assuming considerable intrafractional patient motion, rotational errors could affect target coverage and OAR-sparing especially in multiple metastases VMAT-SRS cranial cases, which often involve the concurrent irradiation of off-axis targets. This work aims to study the dosimetric impact of rotational errors in such applications, under a comparative perspective between the single- and two-isocenter treatment techniques. METHODS: Ten patients (36 metastases) were included in this study. Challenging cases were only considered, with several targets lying in close proximity to OARs. Two multiarc VMAT plans per patient were prepared, involving one and two isocenters, serving as the reference plans. Different degrees of angular offsets at various orientations were introduced, simulating rotational errors. Resulting dose distributions were evaluated and compared using commonly employed dose-volume and plan quality indices. RESULTS: For single-isocenter plans and 1° rotations, plan quality indices, such as coverage, conformity index and D95% , deteriorated significantly (>5%) for distant targets from the isocenter (at> 4-6 cm). Contrarily, for two-isocenter plans, target distances to nearest isocenter were always shorter (≤4 cm), and, consequently, 1° errors were well-tolerated. In the most extreme case considered (2° around all axes) conformity index deteriorated by on-average 7.2%/cm of distance to isocenter, if one isocenter is used, and 2.6%/cm, for plans involving two isocenters. The effect is, however, strongly associated with target volume. Regarding OARs, for single-isocenter plans, significant increase (up to 63%) in Dmax and D0.02cc values was observed for any angle of rotation. Plans that could be considered clinically unacceptable were obtained even for the smallest angle considered, although rarer for the two-isocenter planning approach. CONCLUSION: Limiting the lesion-to-isocenter distance to ≤4 cm by introducing additional isocenter(s) appears to partly mitigate severe target underdosage, especially for smaller target sizes. If OAR-sparing is also a concern, more stringent rotational error tolerances apply.

Quality assurance protocol for superficial and deep hyperthermia systems established by the Hellenic Association of Medical Physicists (HAMP) in cooperation with the Hellenic Society of Oncologic Hyperthemia (HSOH): A study based on European Society for Hyperthermic Oncology (ESHO) quality assurance guidelines.

PURPOSE: During the last years hyperthermia is a developing therapeutic modality in Greece. Quality assurance (QA) procedures are essential for ensuring the correct operation of the hyperthermia system and therefore the selective heating of the tumor with minimum toxicity to the surrounding healthy tissues. The European Society for Hyperthermic Oncology (ESHO) has proposed QA guidelines for superficial as well as deep hyperthermia systems. The purpose of this study was to describe the adapted QA protocol for superficial and deep hyperthermia systems established in Greece. METHODS: A working group was created by the Hellenic Association of Medical Physicists (HAMP) for the proposal of QA guidelines for superficial and deep hyperthermia systems. A review of the protocol proposed by ESHO, together with the existing protocols in other European and International centers, as well as protocols suggested by European or International organizations, was performed. Then, a protocol was suggested, describing procedures for QA according to the current technology and the existing equipment used in Greece. RESULTS: A protocol describing the procedures for QA of superficial and deep hyperthermia systems was proposed. These procedures aim to evaluate the correct operation of the device, the thermometric system, the generator, the incorporated power meter and the applicators. It will also ensure the electrical safety of the devices. CONCLUSIONS: The proposed protocol, applied by medical physicists in Greece, will ensure an efficient treatment with safety and minimum adverse effects. This protocol has been approved by the Hellenic Society of Oncologic Hyperthermia and the Hellenic Association of Medical Physicists.

A Web simulation of medical image reconstruction and processing as an educational tool.

Web educational resources integrating interactive simulation tools provide students with an in-depth understanding of the medical imaging process. The aim of this work was the development of a purely Web-based, open access, interactive application, as an ancillary learning tool in graduate and postgraduate medical imaging education, including a systematic evaluation of learning effectiveness. The pedagogic content of the educational Web portal was designed to cover the basic concepts of medical imaging reconstruction and processing, through the use of active learning and motivation, including learning simulations that closely resemble actual tomographic imaging systems. The user can implement image reconstruction and processing algorithms under a single user interface and manipulate various factors to understand the impact on image appearance. A questionnaire for pre- and post-training self-assessment was developed and integrated in the online application. The developed Web-based educational application introduces the trainee in the basic concepts of imaging through textual and graphical information and proceeds with a learning-by-doing approach. Trainees are encouraged to participate in a pre- and post-training questionnaire to assess their knowledge gain. An initial feedback from a group of graduate medical students showed that the developed course was considered as effective and well structured. An e-learning application on medical imaging integrating interactive simulation tools was developed and assessed in our institution.

Dosimetry in 1.5 T MR-Linacs: Monte Carlo determination of magnetic field correction factors and investigation of the air gap effect.

BACKGROUND: In Magnetic Resonance-Linac (MR-Linac) dosimetry formalisms, a new correction factor, kB,Q , has been introduced to account for corresponding changes to detector readings under the beam quality, Q, and the presence of magnetic field, B. PURPOSE: This study aims to develop and implement a Monte Carlo (MC)-based framework for the determination of kB,Q correction factors for a series of ionization chambers utilized for dosimetry protocols and dosimetric quality assurance checks in clinical 1.5 T MR-Linacs. Their dependencies on irradiation setup conditions are also investigated. Moreover, to evaluate the suitability of solid phantoms for dosimetry checks and end-to-end tests, changes to the detector readings due to the presence of small asymmetrical air gaps around the detector's tip are quantified. METHODS: Phase space files for three irradiation fields of the ELEKTA Unity 1.5 T/7 MV flattening-filter-free MR-Linac were provided by the manufacturer and used as source models throughout this study. Twelve ionization chambers (three farmer-type and nine small-cavity detectors, from three manufacturers) were modeled (including their dead volume) using the EGSnrc MC code package. kB,Q values were calculated for the 10 × 10 cm2 irradiation field and for four cardinal orientations of the detectors' axes with respect to the 1.5 T magnetic field. Potential dependencies of kB,Q values with respect to field size, depth, and phantom material were investigated by performing additional simulations. Changes to the detectors' readings due to the presence of small asymmetrical air gaps (0.1 up to 1 mm) around the chambers' sensitive volume in an RW3 solid phantom were quantified for three small-cavity chambers and two orientations. RESULTS: For both parallel (to the magnetic field) orientations, kB,Q values were found close to unity. The maximum correction needed was 1.1%. For each detector studied, the kB,Q values calculated for the two parallel orientations agreed within uncertainties. Larger corrections (up to 5%) were calculated when the detectors were oriented perpendicularly to the magnetic field. Results were compared with corresponding ones found in the literature, wherever available. No considerable dependence of kB,Q with respect to field size (down to 3 × 3 cm2 ), depth, or phantom material was noticed, for the detectors investigated. As compared to the perpendicular one, in the parallel to the magnetic field orientation, the air gap effect is minimized but is still considerable even for the smallest air gap considered (0.1 mm). CONCLUSION: For the 10 × 10 cm2 field, magnetic field correction factors for 12 ionization chambers and four orientations were determined. For each detector, the kB,Q value may be also applied for dosimetry procedures under different irradiation parameters provided that the orientation is taken into account. Moreover, if solid phantoms are used, even the smallest asymmetrical air gap may still bias small-cavity chamber response. This work substantially expands the availability and applicability of kB,Q correction factors that are detector- and orientation-specific, enabling more options in MR-Linac dosimetry checks, end-to-end tests, and quality assurance protocols.

Dose-response dependencies of OSL dosimeters in conventional linacs and 1.5T MR-linacs: an experimental and Monte Carlo study.

Objective. This work focuses on the optically stimulated luminescence dosimetry (OSLD) dose-response characterization, with emphasis on 1.5T MR-Linacs.Approach. Throughout this study, the nanoDots OSLDs (Landauer, USA) were considered. In groups of three, the mean OSLD response was measured in a conventional linac and an MR-Linac under various irradiation conditions to investigate (i) dose-response linearity with and without the 1.5T magnetic field, (ii) signal fading rate and its dependencies, (iii) beam quality, detector orientation and dose rate dependencies in a conventional linac, (iii) potential MR imaging related effects on OSLD response and (iv) detector orientation dependence in an MR-Linac. Monte Carlo calculations were performed to further quantify angular dependence after rotating the detector around its central axis parallel to the magnetic field, and determine the magnetic field correction factors,kB,Q,for all cardinal detector orientations.Main results. OSLD dose-response supralinearity in an MR-Linac setting was found to agree within uncertainties with the corresponding one in a conventional linac, for the axial detector orientation investigated. Signal fading rate does not depend on irradiation conditions for the range of 3-30 d considered. OSLD angular (orientation) dependence is more pronounced under the presence of a magnetic field. OSLDs irradiated with and without real-time T2w MR imaging enabled during irradiation yielded the same response within uncertainties.kB,Qvalues were determined for all three cardinal orientations. Corrections needed reached up to 6.4%. However, if OSLDs are calibrated in the axial orientation and then irradiated in an MR-Linac placed again in the axial orientation (perpendicular to the magnetic field), then simulations suggest thatkB,Qcan be considered unity within uncertainties, irrespective of the incident beam angle.Significance. This work contributes towards OSLD dose-response characterization and relevant correction factors availability. OSLDs are suitable for QA checks in MR-based beam gating applications andin vivodosimetry in MR-Linacs.

Scattered Radiation Distribution Utilizing Three Different Cone-Beam Computed Tomography Devices for Maxillofacial Diagnostics: A Research Study.

This study aimed to estimate scattered radiation and its spatial distribution around three cone-beam computed tomography (CBCT) devices, in order to determine potential positions for an operator to stand if they needed to be inside the CBCT room. The following devices were tested: Morita Accuitomo (CBCT1), Newtom Giano HR (CBCT2), Newtom VGi (CBCT3). Scattered radiation measurements were performed using different kVp, mA, and Field of View (FOV) options. An anthropomorphic phantom (NATHANIA) was placed inside the X-ray gantry to simulate clinical conditions. Scattered measurements were taken with the Inovision model 451P Victoreen ionization chamber once placed at fixed distances from each irradiation isocenter, away from the primary beam. A statistically significant (p < 0.001) difference was found in the mean value of the scattered radiation estimations between the CBCT devices. Scattered radiation was reduced with a different rate for each CBCT device as distance was increased. For CBCT1 the reduction was 0.047 µGy, for CBCT2 it was 0.036 µGy, and for CBCT3 it was 0.079 µGy, for every one meter from the X-ray gantry. Therefore, at certain distances from the central X-ray, the scattered radiation was below the critical level of 1 mGy, which is defined by the radiation protection guidelines as the exposure radiation limit of the general population. Consequently, an operator could stay inside the room accompanying the patient being scanned, if necessary.

Surgical Gesture Recognition in Laparoscopic Tasks Based on the Transformer Network and Self-Supervised Learning.

In this study, we propose a deep learning framework and a self-supervision scheme for video-based surgical gesture recognition. The proposed framework is modular. First, a 3D convolutional network extracts feature vectors from video clips for encoding spatial and short-term temporal features. Second, the feature vectors are fed into a transformer network for capturing long-term temporal dependencies. Two main models are proposed, based on the backbone framework: C3DTrans (supervised) and SSC3DTrans (self-supervised). The dataset consisted of 80 videos from two basic laparoscopic tasks: peg transfer (PT) and knot tying (KT). To examine the potential of self-supervision, the models were trained on 60% and 100% of the annotated dataset. In addition, the best-performing model was evaluated on the JIGSAWS robotic surgery dataset. The best model (C3DTrans) achieves an accuracy of 88.0%, a 95.2% clip level, and 97.5% and 97.9% (gesture level), for PT and KT, respectively. The SSC3DTrans performed similar to C3DTrans when training on 60% of the annotated dataset (about 84% and 93% clip-level accuracies for PT and KT, respectively). The performance of C3DTrans on JIGSAWS was close to 76% accuracy, which was similar to or higher than prior techniques based on a single video stream, no additional video training, and online processing.

On the potential of 2D ion chamber arrays for high-dose rate remote afterloading brachytherapy quality assurance.

Objective. To investigate the potential of 2D ion chamber arrays to serve as a standalone tool for the verification of source strength, positioning and dwell time, within the framework of192Ir high-dose rate brachytherapy device quality assurance (QA).Approach.A commercially available ion chamber array was used. Fitting of a 2D Lorentzian peak function to experimental data from a multiple source dwell position irradiation on a frame-by-frame basis, facilitated tracking of the source center orthogonal projection on the array plane. For source air kerma strength verification, Monte Carlo simulation was employed to obtain a chamber array- and source-specific correction factor of calibration with a 6 MV photon beam. This factor converted the signal measured by each ion chamber element to air kerma in free space. A source positioning correction was also applied to lift potential geometry mismatch between experiment and Monte Carlo simulation.Main results.Spatial and temporal accuracy of source movement was verified within 0.5 mm and 0.02 s, respectively, in compliance with the test endpoints recommended by international professional societies. The source air kerma strength was verified experimentally within method uncertainties estimated as 1.44% (k = 1). The source positioning correction method employed did not introduce bias to experimental results of irradiations where source positioning was accurate. Development of a custom jig attachable to the chamber array for accurate and reproducible experimental set up would improve testing accuracy and obviate the need for source positioning correction in air kerma strength verification.Significance.Delivery of a single irradiation plan, optimized based on results of this work, to a 2D ion chamber array can be used for concurrent testing of source position, dwell time and air kerma strength, and the procedure can be expedited through automation. Chamber arrays merit further study in treatment planning QA and real time,in vivodose verification.

Monte Carlo-Based Radiobiological Investigation of the Most Optimal Ion Beam Forming SOBP for Particle Therapy.

Proton (p) and carbon (C) ion beams are in clinical use for cancer treatment, although other particles such as He, Be, and B ions have more recently gained attention. Identification of the most optimal ion beam for radiotherapy is a challenging task involving, among others, radiobiological characterization of a beam, which is depth-, energy-, and cell type- dependent. This study uses the FLUKA and MCDS Monte Carlo codes in order to estimate the relative biological effectiveness (RBE) for several ions of potential clinical interest such as p, 4He, 7Li, 10Be, 10B, and 12C forming a spread-out Bragg peak (SOBP). More specifically, an energy spectrum of the projectiles corresponding to a 5-cm SOBP at a depth of 8 cm was used. All secondary particles produced by the projectiles were considered and RBE was determined based on radiation-induced Double Strand Breaks (DSBs), as calculated by MCDS. In an attempt to identify the most optimal ion beam, using the latter data, biological optimization was performed and the obtained depth-dose distributions were inter-compared. The results showed that 12C ions are more effective inside the SOBP region, which comes at the expense of higher dose values at the tail (i.e., after the SOBP). In contrast, p beams exhibit a higher DSOPB/DEntrance ratio, if physical doses are considered. By performing a biological optimization in order to obtain a homogeneous biological dose (i.e., dose × RBE) in the SOBP, the corresponding advantages of p and 12C ions are moderated. 7Li ions conveniently combine a considerably lower dose tail and a DSOPB/DEntrance ratio similar to 12C. This work contributes towards identification of the most optimal ion beam for cancer therapy. The overall results of this work suggest that 7Li ions are of potential interest, although more studies are needed to demonstrate the relevant advantages. Future work will focus on studying more complex beam configurations.

Quantification of Nanoscale Dose Enhancement in Gold Nanoparticle-Aided External Photon Beam Radiotherapy.

The recent progress in Nanotechnology has introduced Gold Nanoparticles (AuNPs) as promising radiosensitizing agents in radiation oncology. This work aims to estimate dose enhancement due to the presence of AuNPs inside an irradiated water region through Monte Carlo calculations. The GATE platform was used to simulate 6 MV photon histories generated from a TrueBeam® linear accelerator with and without a Flattening Filter (FF) and model AuNPs clusters. The AuNPs size, concentration and distribution pattern were examined. To investigate different clinical irradiation conditions, the effect of field size, presence of FF and placement of AuNPs in water were evaluated. The range of Dose Enhancement Factors (DEF = DoseAu/DoseWater) calculated in this study is 0.99 ± 0.01-1.26 ± 0.02 depending on photon beam quality, distance from AuNPs surface, AuNPs size and concentration and pattern of distribution. The highest DEF is reported for irradiation using un-flattened photon beams and at close distances from AuNPs. The obtained findings suggest that dose deposition could be increased in regions that represent whole cells or subcellular targets (mitochondria, cell nucleus, etc.). Nevertheless, further and consistent research is needed in order to make a step toward AuNP-aided radiotherapy in clinical practice.

Geometric distortion assessment in 3T MR images used for treatment planning in cranial Stereotactic Radiosurgery and Radiotherapy.

Magnetic Resonance images (MRIs) are employed in brain Stereotactic Radiosurgery and Radiotherapy (SRS/SRT) for target and/or critical organ localization and delineation. However, MRIs are inherently distorted, which also impacts the accuracy of the Magnetic Resonance Imaging/Computed Tomography (MRI/CT) co-registration process. In this phantom-based study, geometric distortion is assessed in 3T T2-weighted images (T2WIs), while the efficacy of an MRI distortion correction technique is also evaluated. A homogeneous polymer gel-filled phantom was CT-imaged before being irradiated with 26 4-mm Gamma Knife shots at predefined locations (reference control points). The irradiated phantom was MRI-scanned at 3T, implementing a T2-weighted protocol suitable for SRS/SRT treatment planning. The centers of mass of all shots were identified in the 3D image space by implementing an iterative localization algorithm and served as the evaluated control points for MRI distortion detection. MRIs and CT images were spatially co-registered using a mutual information algorithm. The inverse transformation matrix was applied to the reference control points and compared with the corresponding MRI-identified ones to evaluate the overall spatial accuracy of the MRI/CT dataset. The mean image distortion correction technique was implemented, and resulting MRI-corrected control points were compared against the corresponding reference ones. For the scanning parameters used, increased MRI distortion (>1mm) was detected at areas distant from the MRI isocenter (>5cm), while median radial distortion was 0.76mm. Detected offsets were slightly higher for the MRI/CT dataset (0.92mm median distortion). The mean image distortion correction improves geometric accuracy, but residual distortion cannot be considered negligible (0.51mm median distortion). For all three datasets studied, a statistically significant positive correlation between detected spatial offsets and their distance from the MRI isocenter was revealed. This work contributes towards the wider adoption of 3T imaging in SRS/SRT treatment planning. The presented methodology can be employed in commissioning and quality assurance programmes of corresponding treatment workflows.

Early Prediction of Planning Adaptation Requirement Indication Due to Volumetric Alterations in Head and Neck Cancer Radiotherapy: A Machine Learning Approach.

BACKGROUND: During RT cycles, the tumor response pattern could affect tumor coverage and may lead to organs at risk of overdose. As such, early prediction of significant volumetric changes could therefore reduce potential radiation-related adverse effects. Nevertheless, effective machine learning approaches based on the radiomic features of the clinically used CBCT images to determine the tumor volume variations due to RT not having been implemented so far. METHODS: CBCT images from 40 HN cancer patients were collected weekly during RT treatment. From the obtained images, the Clinical Target Volume (CTV) and Parotid Glands (PG) regions of interest were utilized to calculate 104 delta-radiomics features. These features were fed on a feature selection and classification procedure for the early prediction of significant volumetric alterations. RESULTS: The proposed framework was able to achieve 0.90 classification performance accuracy while detecting a small subset of discriminative characteristics from the 1st week of RT. The selected features were further analyzed regarding their effects on temporal changes in anatomy and tumor response modeling. CONCLUSION: The use of machine learning algorithms offers promising perspectives for fast and reliable early prediction of large volumetric deviations as a result of RT treatment, exploiting hidden patterns in the overall anatomical characteristics.

Assessment of Radiation-Induced Bladder and Bowel Cancer Risks after Conventionally and Hypo-Fractionated Radiotherapy for the Preoperative Management of Rectal Carcinoma.

Preoperative management of rectal carcinoma can be performed by employing either conventionally or hypo-fractionated Radiotherapy (CFRT or HFRT, respectively), delivered by Intensity Modulated Radiotherapy (IMRT) or Volumetric Modulated Arc Therapy (VMAT) plans, employing 6 MV or 10 MV photon beams. This study aims to dosimetrically and radiobiologically compare all available approaches, with emphasis on the risk of radiation-induced second cancer to the bladder and bowel. Computed Tomography (CT) scans and relevant radiotherapy contours from 16 patients were anonymized and analyzed retrospectively. For each case, CFRT of 25 × 2 Gy and HFRT of 5 × 5 Gy were both considered. IMRT and VMAT plans using 6 MV and 10 MV photons were prepared. Plan optimization was performed, considering all clinically used plan quality indices and dose-volume constraints for the critical organs. Resulting dose distributions were analyzed and compared. Moreover, the Lifetime Attributable Risk (LAR) for developing radiation-induced bladder and bowel malignancies were assessed using a non-linear mechanistic model, assuming patient ages at treatment of 45, 50, 55 and 60 years. All 128 plans created were clinically acceptable. Risk of second bladder cancer reached 0.26% for HFRT (5 × 5 Gy) and 0.19% for CFRT (25 × 2 Gy) at the age of 45. Systematically higher risks were calculated for HFRT (5 × 5 Gy) as compared to CFRT (25 × 2 Gy), with 6 MV photons resulting in slightly increased LAR, as well. Similar or equal bowel cancer risks were calculated for all techniques and patient ages investigated (range 0.05-0.14%). This work contributes towards radiotherapy treatment protocol selection criteria for the preoperative irradiation of rectal carcinoma. However, more studies are needed to establish the associated radiation-induced risk of each RT protocol.

Dose perturbation in the radiotherapy of breast cancer patients implanted with the Magna-Site: a Monte Carlo study.

External beam radiation therapy (RT) is often offered to breast cancer patients after surgical mastectomy followed by breast reconstruction with silicone implants. In some cases, the RT is administered while the patient is still implanted with a temporary tissue expander including a high-density metallic port, which is expected to affect the planned dose distribution. This work uses Monte Carlo (MC) simulation in order to evaluate the aforementioned effect when the McGhan Style 133 Tissue Expander with the Magna-Site injection port is used. Simulations have been performed on a patient model built using the actual CT images of the patient for two irradiation schemes, involving two tangential photon beams of 6 MV and 18 MV respectively. MC results show that the presence of the Magna-Site within the two irradiation fields leads to an overall reduction of absorbed dose for points lying in the shadow of the metallic port (relative to each of the opposing beams). The relative reduction compared to dose results without the expander in place ranges from 7% to 13% for the 6 MV beam and is around 6% for the 18 MV photon beam. However, in the close vicinity of the metallic port, increased absorbed doses are observed, due to the increase of secondary electrons emerging from the metallic part of the insert.

Lung cancer histology classification from CT images based on radiomics and deep learning models.

Adenocarcinoma (AC) and squamous cell carcinoma (SCC) are frequent reported cases of non-small cell lung cancer (NSCLC), responsible for a large fraction of cancer deaths worldwide. In this study, we aim to investigate the potential of NSCLC histology classification into AC and SCC by applying different feature extraction and classification techniques on pre-treatment CT images. The employed image dataset (102 patients) was taken from the publicly available cancer imaging archive collection (TCIA). We investigated four different families of techniques: (a) radiomics with two classifiers (kNN and SVM), (b) four state-of-the-art convolutional neural networks (CNNs) with transfer learning and fine tuning (Alexnet, ResNet101, Inceptionv3 and InceptionResnetv2), (c) a CNN combined with a long short-term memory (LSTM) network to fuse information about the spatial coherency of tumor's CT slices, and (d) combinatorial models (LSTM + CNN + radiomics). In addition, the CT images were independently evaluated by two expert radiologists. Our results showed that the best CNN was Inception (accuracy = 0.67, auc = 0.74). LSTM + Inception yielded superior performance than all other methods (accuracy = 0.74, auc = 0.78). Moreover, LSTM + Inception outperformed experts by 7-25% (p < 0.05). The proposed methodology does not require detailed segmentation of the tumor region and it may be used in conjunction with radiological findings to improve clinical decision-making. Lung cancer histology classification from CT images based on CNN + LSTM.

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.