Search

Objective. Plastic scintillation dosimeters (PSDs) are highly suitable for real-time dosimetry on the MR-linac. For optimal performance, the primary signal (scintillation) needs to be separated from secondary optical effects (Cerenkov, fluorescence and optical fiber attenuation). This requires a spectral separation approach and careful calibration. Currently, the 'classic' calibration is a multi-step procedure using both kV and MV x-ray sources, requiring an uninterrupted optical connection between the dosimeter and read-out system, complicating efficient use of PSDs. Therefore, we present a more time-efficient and more practical novel calibration technique for PSDs suitable for MR-linac dosimetry. Approach. The novel calibration relies on prior spectral information combined with two 10 × 10 cm 2 field irradiations on the 1.5 T MR-linac. Performance of the novel calibration technique was evaluated focusing on its reproducibility, performance characteristics (repeatability, linearity, dose rate dependency, output factors, angular response and detector angle dependency) and IMRT deliveries. To investigate the calibration stability over time, prior spectral information up to 315 days old was used. To quantify the time efficiency, each step of the novel and classic calibration was timed. Main results. The novel calibration showed a high reproducibility with a maximum relative standard deviation of 0.2%. The novel method showed maximum differences of 1.2% compared to the gold-standard calibration, while reusing old classic calibrations after reconnecting fibers showed differences up to 3.0%. The novel calibration improved time efficiency from 105 to 30 min compared to the classic method. Significance. The novel calibration method showed a gain in time efficiency and practicality while preserving the dosimetric accuracy. Therefore, this method can replace the traditional method for PSDs suitable for MR-linac dosimetry, using prior spectral information of up to a year. This novel calibration facilitates reconnecting the detector to the read-out system which would lead to unacceptable dosimetric results with the classic calibration method., (Creative Commons Attribution license.)

Medical image registration is an integral part of various clinical applications including image guidance, motion tracking, therapy assessment and diagnosis. We present a robust approach for mono-modal and multi-modal medical image registration. To this end, we propose the novel shape operator based local image distance (SOLID) which estimates the similarity of images by comparing their second-order curvature information. Our similarity metric is rigorously tailored to be suitable for comparing images from different medical imaging modalities or image contrasts. A critical element of our method is the extraction of local features using higher-order shape information, enabling the accurate identification and registration of smaller structures. In order to assess the efficacy of the proposed similarity metric, we have implemented a variational image registration algorithm that relies on the principle of matching the curvature information of the given images. The performance of the proposed algorithm has been evaluated against various alternative state-of-the-art variational registration algorithms. Our experiments involve mono-modal as well as multi-modal and cross-contrast co-registration tasks in a broad variety of anatomical regions. Compared to the evaluated alternative registration methods, the results indicate a very favorable accuracy, precision and robustness of the proposed SOLID method in various highly challenging registration tasks., (Creative Commons Attribution license.)

Objective . Intrafraction motion is a major concern for the safety and effectiveness of high dose stereotactic body radiotherapy (SBRT) in the upper abdomen. In this study, the impact of the intrafraction motion on the delivered dose was assessed in a patient group that underwent MR-guided radiotherapy for upper abdominal malignancies with an abdominal corset. Approach . Fast online 2D cine MRI was used to extract tumor motion during beam-on time. These tumor motion profiles were combined with linac log files to reconstruct the delivered dose in 89 fractions of MR-guided SBRT in twenty patients. Aside the measured tumor motion, motion profiles were also simulated for a wide range of respiratory amplitudes and drifts, and their subsequent dosimetric impact was calculated in every fraction. Main results . The average (SD) D 99% of the gross tumor volume (GTV), relative to the planned D 99% , was 0.98 (0.03). The average (SD) relative D 0.5 cc of the duodenum, small bowel and stomach was 0.99 (0.03), 1.00 (0.03), and 0.97 (0.05), respectively. No correlation of respiratory amplitude with dosimetric impact was observed. Fractions with larger baseline drifts generally led to a larger uncertainty of dosimetric impact on the GTV and organs at risk (OAR). The simulations yielded that the delivered dose is highly dependent on the direction of on baseline drift. Especially in anatomies where the OARs are closely abutting the GTV, even modest LR or AP drifts can lead to substantial deviations from the planned dose. Significance . The vast majority of the fractions was only modestly impacted by intrafraction motion, increasing our confidence that MR-guided SBRT with abdominal compression can be safely executed for patients with abdominal tumors, without the use of gating or tracking strategies., (Creative Commons Attribution license.)

The treatment of cancer with proton radiation therapy was first suggested in 1946 followed by the first treatments in the 1950s. As of 2020, almost 200 000 patients have been treated with proton beams worldwide and the number of operating proton therapy (PT) facilities will soon reach one hundred. PT has long moved from research institutions into hospital-based facilities that are increasingly being utilized with workflows similar to conventional radiation therapy. While PT has become mainstream and has established itself as a treatment option for many cancers, it is still an area of active research for various reasons: the advanced dose shaping capabilities of PT cause susceptibility to uncertainties, the high degrees of freedom in dose delivery offer room for further improvements, the limited experience and understanding of optimizing pencil beam scanning, and the biological effect difference compared to photon radiation. In addition to these challenges and opportunities currently being investigated, there is an economic aspect because PT treatments are, on average, still more expensive compared to conventional photon based treatment options. This roadmap highlights the current state and future direction in PT categorized into four different themes, 'improving efficiency', 'improving planning and delivery', 'improving imaging', and 'improving patient selection'., (© 2021 Institute of Physics and Engineering in Medicine.)

Image-guided radiation therapy (IGRT) allows radiation dose deposition with a high degree of geometric accuracy. Previous studies have demonstrated that such therapies may benefit from the employment of deformable image registration (DIR) algorithms, which allow both the automatic tracking of anatomical changes and accumulation of the delivered radiation dose over time. In order to ensure patient care and safety, however, the estimated deformations must be subjected to stringent quality assurance (QA) measures. In the present study we propose to extend the state-of-the-art methodology for QA of DIR algorithms by a set of novel biomechanical criteria. The proposed biomechanical criteria imply the calculation of the normal and shear mechanical stress, which would occur within the observed tissues as a result of the estimated deformations. The calculated stress is then compared to plausible physiological limits, providing thus the anatomical plausibility of the estimated deformations. The criteria were employed for the QA of three DIR algorithms in the context of abdominal conebeam computed tomography and magnetic resonance radiotherapy guidance. An initial evaluation of organ boundary alignment capabilities indicated that all three algorithms perform similarly. However, an analysis of the deformations within the organ boundaries with respect to the proposed biomechanical QA criteria revealed different degrees of anatomical plausibility. Additionally, it was demonstrated that violations of these criteria are also indicative of errors within the dose accumulation process. The proposed QA criteria, therefore, provide a tissue-dependent assessment of the anatomical plausibility of the deformations estimated by DIR algorithms, showcasing potential in ensuring patient safety for future adaptive IGRT treatments.

This work aimed to quantify any principal magnetic field (B 0 ) inhomogeneity and changes in MR image geometric distortion with continuous linac gantry rotation on an Elekta Unity MR-linac. This situation occurs for around a second between treatment beams during current image guided radiotherapy treatment and would occur frequently in foreseeable real-time adaptive radiotherapy treatment. Pixel by pixel maps of B 0 inhomogeneity were obtained via repeated high temporal resolution pulse sequences with the linac gantry static at 36 gantry angles spaced ten degrees apart, and in continuous rotation at both 1 and 2 rpm. Individual B 0 maps were subtracted from average maps across all data and the residual peak to peak inhomogeneity was calculated for each. The bulk geometric shift and change in physical extent of a 10 cm diameter spherical flood phantom during continuous linac gantry rotation at 1 and 2 rpm was compared to the static gantry case for two pulse sequences: the real-time clinical monitoring bFFE sequence and a non-clinical EPI sequence, chosen for its susceptibility to geometric distortion. The peak to peak inhomogeneity in the deviation-from-average ppm maps, plotted against gantry angle with the gantry in continuous rotation at 1 and 2 rpm were negligibly different from equivalent data obtained with the gantry static. The real-time clinical monitoring pulse sequence was shown to give negligible geometric distortion during continuous gantry motion, whilst a non-clinical EPI sequence showed bulk shifts of the order of one pixel and gantry angle dependent changes in extent, demonstrating the sensitivity of the chosen method. MR imaging on the Elekta Unity MR-Linac with the gantry in continuous motion is negligibly different from the static gantry case with current clinical pulse sequences. Real-time tracking and treatment plan adaptation using MR images obtained with the linac gantry in motion is possible.

The out-of-field surface dose contribution due to backscattered or ejected electrons, focused by the magnetic field, is evaluated in this work. This electron streaming effect (ESE) can contribute to out-of-field skin doses in orthogonal magnetic resonance guided radiation therapy machines. Using the EGSnrc Monte Carlo package, a phantom is set-up along the central axis of an incident 10 [Formula: see text] 10 cm 2 7 MV FFF photon beam. The phantom exit or entry surface is inclined with respect to the magnetic field, and an out-of-field water panel is positioned 10 cm away from, and centered on, the isocenter. The doses from streaming backscattered or ejected electrons, for either a 0.35 T or 1.5 T magnetic field, are evaluated in the out-of-field water panel for surface inclines of 10, 30, and 45°. The magnetic field focuses electrons emitted from the inclined phantom. Dose distributions at the surface of the out-of-field water panel are sharper in the 1.5 T magnetic field as compared to 0.35 T. The maximum doses for the 0.35 T simulations are 23.2%, 37.8%, and 39.0% for the respective 10, 30, and 45° simulations. For 1.5 T, for the same angles, the maximum values are 17.1%, 29.8%, and 35.8%. Dose values drop to below 2% within the first 1 cm of the out-of-field water phantom. The phantom thickness is an important variable in the magnitude of the ESE dose. The ESE can produce large out-of-field skin doses and must be a consideration in treatment planning in the MRgRT work-flow. Treatments often include multiple beams which will serve to spread out the effect, and many beams, such as anterior-posterior, will reduce the skin dose due to the ESE. A 1 cm thick shielding of either a bolus placed on the patient or mounted on the present RF coils would greatly reduce the ESE dose contributions. Further exploration of the capabilities of treatment planning systems to screen for this effect is required.

The purpose of this study was to develop a method enabling synthetic computed tomography (sCT) generation of the whole abdomen using magnetic resonance imaging (MRI) scans of pediatric patients with abdominal tumors. The proposed method relies on an automatic atlas-based segmentation of bone and lungs followed by an MRI intensity to synthetic Hounsfield unit conversion. Separate conversion algorithms were used for bone, lungs and soft-tissue. Rigidly registered CT and T2-weighted MR images of 30 patients in treatment position and with the same field of view were used for the evaluation of the atlas and the conversion algorithms. The dose calculation accuracy of the generated sCTs was verified for volumetric modulated arc therapy (VMAT) and pencil beam scanning (PBS). VMAT and PBS plans were robust optimized on an internal target volume (ITV) against a patient set-up uncertainty of 5 mm. Average differences between CT and sCT dose calculations for the ITV V 95% were 0.5% (min 0.0%; max 5.0%) and 0.0% (min  -0.1%; max 0.1%) for VMAT and PBS dose distributions, respectively. Average differences for the mean dose to the organs at risk were  <0.2% (min  -0.6%; max 1.2%) and  <0.2% (min  -2.0%; max 2.6%) for VMAT and PBS dose distributions, respectively. Results show that MRI-only photon and proton dose calculations are feasible for children with abdominal tumors.

The purpose of this study is to investigate the attenuation characteristics of a novel radiofrequency (RF) coil, which is the first coil that is solely dedicated to MR guided radiotherapy with a 1.5 T MR-linac. Additionally, we investigated the impact of the treatment beam on the MRI performance of this RF coil. First, the attenuation characteristics of the RF coil were characterized. Second, we investigated the impact of the treatment beam on the MRI performance of the RF coil. We additionally demonstrated the ability of the anterior coil to attenuate returning electrons and thereby reducing the dose to the skin at the distal side of the treatment beam. Intensity modulated radiation therapy simulation of a clinically viable treatment plan for spinal bone metastasis shows a decrease of the dose to the planned tumor volume of 1.8% as a result of the MR coil around the patient. Ionization chamber and film measurements show that the anterior and posterior coil attenuate the beam homogeneously by 0.4% and 2.2%, respectively. The impact of the radiation resulted in a slight drop of the time-course signal-to-noise ratio and was dependent on imaging parameters. However, we could not observe any image artifacts resulting from this irradiation in any situation. In conclusion, the investigated MR-coil can be utilized for treatments with the 1.5 T-linac system. However, there is still room for improvement when considering both the dosimetric and imaging performance of the coil.

A magnetic resonance (MR)-only radiotherapy workflow can reduce cost, radiation exposure and uncertainties introduced by CT-MRI registration. A crucial prerequisite is generating the so called pseudo-CT (pCT) images for accurate dose calculation and planning. Many pCT generation methods have been proposed in the scope of photon radiotherapy. This work aims at verifying for the first time whether a commercially available photon-oriented pCT generation method can be employed for accurate intensity-modulated proton therapy (IMPT) dose calculation. A retrospective study was conducted on ten prostate cancer patients. For pCT generation from MR images, a commercial solution for creating bulk-assigned pCTs, called MR for Attenuation Correction (MRCAT), was employed. The assigned pseudo-Hounsfield Unit (HU) values were adapted to yield an increased agreement to the reference CT in terms of proton range. Internal air cavities were copied from the CT to minimise inter-scan differences. CT- and MRCAT-based dose calculations for opposing beam IMPT plans were compared by gamma analysis and evaluation of clinically relevant target and organ at risk dose volume histogram (DVH) parameters. The proton range in beam's eye view (BEV) was compared using single field uniform dose (SFUD) plans. On average, a [Formula: see text] mm) gamma pass rate of 98.4% was obtained using a [Formula: see text] dose threshold after adaptation of the pseudo-HU values. Mean differences between CT- and MRCAT-based dose in the DVH parameters were below 1 Gy ([Formula: see text]). The median proton range difference was [Formula: see text] mm, with on average 96% of all BEV dose profiles showing a range agreement better than 3 mm. Results suggest that accurate MR-based proton dose calculation using an automatic commercial bulk-assignment pCT generation method, originally designed for photon radiotherapy, is feasible following adaptation of the assigned pseudo-HU values.

Combining magnetic-resonance imaging (MRI) and proton therapy (PT) using pencil-beam scanning (PBS) may improve image-guided radiotherapy. We aimed at assessing the impact of a magnetic field on PBS-PT plan quality and robustness. Specifically, the robustness against anatomical changes and positioning errors in an MRI-guided scenario with a 30 cm radius 1.5 T magnetic field was studied for prostate PT. Five prostate cancer patients with three consecutive CT images (CT1-3) were considered. Single-field uniform dose PBS-PT plans were generated on the segmented CT1 with Monte-Carlo-based treatment planning software for inverse optimization. Plans were optimized at 90° gantry angle without B-field (no B), with  ±1.5 T B-field (B and minus B), as well as at 81° gantry angle and  +1.5 T (B G81). Plans were re-calculated on aligned CT2 and CT3 to study the impact of anatomical changes. Dose distributions were compared in terms of changes in DVH parameters, proton range and gamma-index pass-rates. To assess the impact of positioning errors, DVH parameters were compared for  ±5 mm CT1 patient shifts in anterior-posterior (AP) and left-right (LR) direction. Proton beam deflection considerably reduced robustness against inter-fractional changes for the B scenario. Range agreement, gamma-index pass-rates and PTV V95% were significantly lower compared to no B. Improved robustness was obtained for minus B and B G81, the latter showing only minor differences to no B. The magnetic field introduced slight dosimetric changes under LR shifts. The impact of AP shifts was considerably larger, and equivalent for scenarios with and without B-field. Results suggest that robustness equivalent to PT without magnetic field can be achieved by adaptation of the treatment parameters, such as B-field orientation (minus B) with respect to the patient and/or gantry angle (B G81). MRI-guided PT for prostate cancer might thus be implemented without compromising robustness compared to state-of-the-art CT-guided PT.

Reduction of motion uncertainty by applying adaptive radiotherapy strategies depends largely on the temporal behavior of this motion. To fully optimize adaptive strategies, insight into target motion is needed. The purpose of this study was to analyze stability and evolution in time of motion uncertainty of both the gross tumor volume (GTV) and clinical target volume (CTV) for patients with rectal cancer. We scanned 16 patients daily during one week, on a 1.5 T MRI scanner in treatment position, prior to each radiotherapy fraction. Single slice sagittal cine MRIs were made at the beginning, middle, and end of each scan session, for one minute at 2 Hz temporal resolution. GTV and CTV motion were determined by registering a delineated reference frame to time-points later in time. The 95th percentile of observed motion (dist95%) was taken as a measure of motion. The stability of motion in time was evaluated within each cine-MRI separately. The evolution of motion was investigated between the reference frame and the cine-MRIs of a single scan session and between the reference frame and the cine-MRIs of several days later in the course of treatment. This observed motion was then converted into a PTV-margin estimate. Within a one minute cine-MRI scan, motion was found to be stable and small. Independent of the time-point within the scan session, the average dist95% remains below 3.6 mm and 2.3 mm for CTV and GTV, respectively 90% of the time. We found similar motion over time intervals from 18 min to 4 days. When reducing the time interval from 18 min to 1 min, a large reduction in motion uncertainty is observed. A reduction in motion uncertainty, and thus the PTV-margin estimate, of 71% and 75% for CTV and tumor was observed, respectively. Time intervals of 15 and 30 s yield no further reduction in motion uncertainty compared to a 1 min time interval.

The application of thermal modelling for hyperthermia and thermal ablation is severely hampered by lack of information about perfusion and vasculature. However, recently, with the advent of sophisticated angiography and dynamic contrast enhanced (DCE) imaging techniques, it has become possible to image small vessels and blood perfusion bringing the ultimate goal of patient specific thermal modelling closer within reach. In this study dynamic contrast enhanced multi-slice CT imaging techniques are employed to investigate the feasibility of this concept for regional hyperthermia treatment of the prostate. The results are retrospectively compared with clinical thermometry data of a patient group from an earlier trial. Furthermore, the role of the prostate vasculature in the establishment of the prostate temperature distribution is studied. Quantitative 3D perfusion maps of the prostate were constructed for five patients using a distributed-parameter tracer kinetics model to analyse dynamic CT data. CT angiography was applied to construct a discrete vessel model of the pelvis. Additionally, a discrete vessel model of the prostate vasculature was constructed of a prostate taken from a human corpse. Three thermal modelling schemes with increasing inclusion of the patient specific physiological information were used to simulate the temperature distribution of the prostate during regional hyperthermia. Prostate perfusion was found to be heterogeneous and T3 prostate carcinomas are often characterized by a strongly elevated tumour perfusion (up to 70-80 ml 100 g(-1) min(-1)). This elevated tumour perfusion leads to 1-2 degrees C lower tumour temperatures than thermal simulations based on a homogeneous prostate perfusion. Furthermore, the comparison has shown that the simulations with the measured perfusion maps result in consistently lower prostate temperatures than clinically achieved. The simulations with the discrete vessel model indicate that significant pre-heating takes place in the prostate capsule vasculature which forms a possible explanation for the discrepancy. Pre-heating in the larger pelvic vessels is very moderate, approximately 0.1-0.3 degrees C. In conclusion, perfusion imaging provides important input for thermal modelling and can be used to obtain a lower limit on the prostate and tumour temperature in regional hyperthermia. However, it is not sufficient to calculate in detail the prostate temperature distribution in individual patients. The prostate vasculature plays such a crucial role that a patient specific discrete vessel model of the prostate vasculature is required.