Pelvic organ motion and dosimetric implications during horizontal patient rotation for prostate radiation therapy.

PURPOSE: Gantry-free radiation therapy systems utilizing patient rotation would be simpler and more cost effective than the conventional gantry-based systems. Such a system could enable the expansion of radiation therapy to meet global demand and reduce capital costs. Recent advances in adaptive radiation therapy could potentially be applied to correct for gravitational deformation during horizontal patient rotation. This study aims to quantify the pelvic organ motion and the dosimetric implications of horizontal rotation for prostate intensity-modulated radiation therapy (IMRT) treatments. METHODS: Eight human participants who previously received prostate radiation therapy were imaged in a clinical magnetic resonance imaging (MRI) scanner using a bespoke patient rotation system (PRS). The patients were imaged every 45 degrees during a full roll rotation (0-360 degrees). Whole pelvic bone, prostate, rectum, and bladder motion were compared to the supine position using dice similarity coefficient (DSC) and mean absolute surface distance (MASD). Prostate centroid motion was compared in the left-right (LR), superior-inferior (SI), and anterior-posterior (AP) direction prior to and following pelvic bone-guided rigid registration. Seven-field prostate IMRT treatment plans were generated for each patient rotation angles under three adaption scenarios: No plan adaption, rigid planning target volume (PTV)-guided alignment to the prostate, and plan re-optimization. Prostate, rectum, and bladder doses were compared for each adaption scenario. RESULTS: Pelvic bone motion within the PRS of up to 53 mm relative to the supine position was observed for some participants. Internal organ motion was greatest at the 180-degree PRS couch angle (prone), with prostate centroid motion range < 2 mm LR, 0 mm to 14 mm SI, and -11 mm to 4 mm AP. Rotation with no adaption of the treatment plan resulted in an underdose to the PTV -- in some instances up to 75% (D95%: 78 ± 0.3 Gy at supine to 20 ± 15.0 Gy at the 225-degree PRS couch angle). Bladder dose was reduced during the rotation by up to 98% (V60 Gy: 15.0 ± 9.4% supine to 0.3 ± 0.5% at the 225-degree PRS couch angle). In some instances, the rectum dose increased during rotation (V60Gy: 20.0 ± 4.5% supine to 25.0 ± 15.0% at the 135-degree PRS couch angle). Rigid PTV-guided alignment resulted in PTV coverage which, though statistically lower (P < 0.05 for all D95% values), was within 1 Gy of the supine plans. Plan re-optimization resulted in a statistically equivalent PTV coverage compared to the supine plans (P > 0.05 for all D95% metrics and all within ±0.4 Gy). For both rigid PTV-guided alignment and plan re-optimization, rectum dose volume metrics were reduced compared to the supine position between the 90- and 225-degree PRS couch angles (P < 0.05). Bladder dose volume metrics were not impacted by rotation. CONCLUSION: Pelvic bone and internal organ motion are present during patient rotation. Rigid PTV-guided alignment to the prostate will be a requirement if prostate IMRT is to be safely delivered using patient rotation. Plan re-optimization for each PRS couch angle to account for anatomical deformations further improves the PTV coverage.

Towards MR-guided electron therapy: Measurement and simulation of clinical electron beams in magnetic fields.

PURPOSE: In the current era of MRI-linac radiotherapy, dose optimization with arbitrary dose distributions is a reality. For the first time, we present new and targeted experiments and modeling to aid in evaluating the potential dose improvements offered with an electron beam mode during MRI-linac radiotherapy. METHODS: Small collimated (1 cm diameter and 1.5 × 1.5 cm2 square) electron beams (6, 12 and 20 MeV) from a clinical linear accelerator (Varian Clinac 2100C) are incident perpendicular and parallel to the strong and localized magnetic fields (0-0.7 T) generated by a permanent magnet device. Gafchromic EBT3 film is placed inside a slab phantom to measure two-dimensional dose distributions. A benchmarked and comprehensive Monte Carlo model (Geant4) is established to directly compare with experiments. RESULTS: With perpendicular fields a 5% narrowing of the beam FWHM and a 10 mm reduction in the 15% isodose penetration is seen for the 20 MeV beam. In the inline setup the penumbral width is reduced by up to 20%, and a local central dose enhancement of 100% is observed. Monte Carlo simulations are in agreement with the measured dose distributions (2% or 2 mm). CONCLUSION: A new range of experiments have been performed to offer insight into how an electron beam mode could offer additional choices in MRI-linac radiotherapy. The work extends on historic studies to bring a successful unified experimental and Monte Carlo modeling approach for studying small field electron beam dosimetry inside magnetic fields. The results suggest further work, particularly on the inline magnetic field scenario.

Anatomical deformation due to horizontal rotation: towards gantry-free radiation therapy.

Gantry-free radiation therapy systems may be simpler and more cost effective, particularly for MRI-guided photon or hadron therapy. This study aims to understand and quantify anatomical deformations caused by horizontal rotation with scan sequences sufficiently short to facilitate integration into an MRI-guided workflow. Rigid and non-rigid pelvic deformations due to horizontal rotation were quantified for a cohort of 8 healthy volunteers using a bespoke patient rotation system and a clinical MRI scanner. For each volunteer a reference scan was acquired at 0° followed by sequential faster scans in 45° increments through to 360°. All fast scans were registered to the 0° image via a three-step process: first, images were aligned using MR visible couch markers. Second, the scans were pre-processed then rigidly registered to the 0° image. Third, the rigidly registered scans were non-rigidly registered to the 0° image to assess soft tissue deformation. The residual differences after rigid and non-rigid registration were determined from the transformation matrix and the deformation vector field, respectively. The rigid registration yielded mean rotations of ⩽2.5° in all cases. The average 3D translational magnitudes range was 5.8 ± 2.9 mm-30.0 ± 11.0 mm. Translations were most significant in the left-right (LR) direction. Smaller translations were observed in the anterior-posterior (AP) and superior-inferior (SI) directions. The maximum deformation magnitudes range was: 10.0 ± 0.9 mm-28.0 ± 2.8 mm and average deformation magnitudes range: 2.3 ± 0.6 mm-7.5 ± 1.0 mm. Average non-rigid deformation magnitude was correlated with BMI (correlation coefficient 0.84, p  = 0.01). Rigid pelvic deformations were most significant in the LR direction but could be accounted for with on-line adjustments. Non-rigid deformations can be significant and will need to be accounted for in order to facilitate the delivery of gantry-free therapy with an automated patient rotation system.

Technical Note: Real-time image-guided adaptive radiotherapy of a rigid target for a prototype fixed beam radiotherapy system.

PURPOSE: Fixed beam radiotherapy systems utilize couch movement and rotation instead of gantry rotation in order to simplify linear accelerator design. We investigate the ability to deliver fixed beam treatments with the same level of clinical accuracy as conventional (rotating beam) treatments using real-time image guidance to maintain this accuracy in the presence of rigid target motion. METHODS: A prototype fixed beam radiotherapy system was built using a standard linac with the beam fixed in the vertical position and a computer controlled rotation stage that rotated a rigid phantom about the superior-inferior axis. Kilovoltage Intrafraction Monitoring (KIM) and real-time beam adaptation with MLC tracking was applied to a five-field IMRT treatment plan with motion introduced to the phantom. The same IMRT treatment was also delivered with real-time adaptation using the conventional rotating beam geometry. Film dosimetry was used to measure the dose delivered with a fixed beam compared to a rotating beam, as well as to compare treatments delivered with and without real-time adaptation. RESULTS: The dose distributions were found to be equivalent between the fixed beam and rotating beam geometry for real-time adaptive radiotherapy using KIM and MLC tracking beam adaptation. Gamma analysis on the films showed agreement >98% using a 2%/2 mm criteria with adaptation for static shifts and periodic motion. CONCLUSIONS: Fixed beam treatments with real-time beam adaptation are dosimetrically equivalent to conventional treatments with a rotating beam, even in the presence of rigid target motion. This suggests that, for a rigid target, the high clinical accuracy of real-time adaptive radiotherapy can be achieved with simpler beam geometry.

An interdimensional correlation framework for real-time estimation of six degree of freedom target motion using a single x-ray imager during radiotherapy.

Increasing evidence suggests that intrafraction tumour motion monitoring needs to include both 3D translations and 3D rotations. Presently, methods to estimate the rotation motion require the 3D translation of the target to be known first. However, ideally, translation and rotation should be estimated concurrently. We present the first method to directly estimate six-degree-of-freedom (6DoF) motion from the target's projection on a single rotating x-ray imager in real-time. This novel method is based on the linear correlations between the superior-inferior translations and the motion in the other five degrees-of-freedom. The accuracy of the method was evaluated in silico with 81 liver tumour motion traces from 19 patients with three implanted markers. The ground-truth motion was estimated using the current gold standard method where each marker's 3D position was first estimated using a Gaussian probability method, and the 6DoF motion was then estimated from the 3D positions using an iterative method. The 3D position of each marker was projected onto a gantry-mounted imager with an imaging rate of 11 Hz. After an initial 110° gantry rotation (200 images), a correlation model between the superior-inferior translations and the five other DoFs was built using a least square method. The correlation model was then updated after each subsequent frame to estimate 6DoF motion in real-time. The proposed algorithm had an accuracy (±precision) of -0.03 ± 0.32 mm, -0.01 ± 0.13 mm and 0.03 ± 0.52 mm for translations in the left-right (LR), superior-inferior (SI) and anterior-posterior (AP) directions respectively; and, 0.07 ± 1.18°, 0.07 ± 1.00° and 0.06 ± 1.32° for rotations around the LR, SI and AP axes respectively on the dataset. The first method to directly estimate real-time 6DoF target motion from segmented marker positions on a 2D imager was devised. The algorithm was evaluated using 81 motion traces from 19 liver patients and was found to have sub-mm and sub-degree accuracy.

Towards real-time MRI-guided 3D localization of deforming targets for non-invasive cardiac radiosurgery.

Radiosurgery to the pulmonary vein antrum in the left atrium (LA) has recently been proposed for non-invasive treatment of atrial fibrillation (AF). Precise real-time target localization during treatment is necessary due to complex respiratory and cardiac motion and high radiation doses. To determine the 3D position of the LA for motion compensation during radiosurgery, a tracking method based on orthogonal real-time MRI planes was developed for AF treatments with an MRI-guided radiotherapy system. Four healthy volunteers underwent cardiac MRI of the LA. Contractile motion was quantified on 3D LA models derived from 4D scans with 10 phases acquired in end-exhalation. Three localization strategies were developed and tested retrospectively on 2D real-time scans (sagittal, temporal resolution 100 ms, free breathing). The best-performing method was then used to measure 3D target positions in 2D-2D orthogonal planes (sagittal-coronal, temporal resolution 200-252 ms, free breathing) in 20 configurations of a digital phantom and in the volunteer data. The 3D target localization accuracy was quantified in the phantom and qualitatively assessed in the real data. Mean cardiac contraction was ⩽ 3.9 mm between maximum dilation and contraction but anisotropic. A template matching approach with two distinct template phases and ECG-based selection yielded the highest 2D accuracy of 1.2 mm. 3D target localization showed a mean error of 3.2 mm in the customized digital phantoms. Our algorithms were successfully applied to the 2D-2D volunteer data in which we measured a mean 3D LA motion extent of 16.5 mm (SI), 5.8 mm (AP) and 3.1 mm (LR). Real-time target localization on orthogonal MRI planes was successfully implemented for highly deformable targets treated in cardiac radiosurgery. The developed method measures target shifts caused by respiration and cardiac contraction. If the detected motion can be compensated accordingly, an MRI-guided radiotherapy system could potentially enable completely non-invasive treatment of AF.

Motion prediction in MRI-guided radiotherapy based on interleaved orthogonal cine-MRI.

In-room cine-MRI guidance can provide non-invasive target localization during radiotherapy treatment. However, in order to cope with finite imaging frequency and system latencies between target localization and dose delivery, tumour motion prediction is required. This work proposes a framework for motion prediction dedicated to cine-MRI guidance, aiming at quantifying the geometric uncertainties introduced by this process for both tumour tracking and beam gating. The tumour position, identified through scale invariant features detected in cine-MRI slices, is estimated at high-frequency (25 Hz) using three independent predictors, one for each anatomical coordinate. Linear extrapolation, auto-regressive and support vector machine algorithms are compared against systems that use no prediction or surrogate-based motion estimation. Geometric uncertainties are reported as a function of image acquisition period and system latency. Average results show that the tracking error RMS can be decreased down to a [0.2; 1.2] mm range, for acquisition periods between 250 and 750 ms and system latencies between 50 and 300 ms. Except for the linear extrapolator, tracking and gating prediction errors were, on average, lower than those measured for surrogate-based motion estimation. This finding suggests that cine-MRI guidance, combined with appropriate prediction algorithms, could relevantly decrease geometric uncertainties in motion compensated treatments.

Dose enhancement in radiotherapy of small lung tumors using inline magnetic fields: A Monte Carlo based planning study.

PURPOSE: To report on significant dose enhancement effects caused by magnetic fields aligned parallel to 6 MV photon beam radiotherapy of small lung tumors. Findings are applicable to future inline MRI-guided radiotherapy systems. METHODS: A total of eight clinical lung tumor cases were recalculated using Monte Carlo methods, and external magnetic fields of 0.5, 1.0, and 3 T were included to observe the impact on dose to the planning target volume (PTV) and gross tumor volume (GTV). Three plans were 6 MV 3D-CRT plans while 6 were 6 MV IMRT. The GTV's ranged from 0.8 to 16 cm(3), while the PTV's ranged from 1 to 59 cm(3). In addition, the dose changes in a 30 cm diameter cylindrical water phantom were investigated for small beams. The central 20 cm of this phantom contained either water or lung density insert. RESULTS: For single beams, an inline magnetic field of 1 T has a small impact in lung dose distributions by reducing the lateral scatter of secondary electrons, resulting in a small dose increase along the beam. Superposition of multiple small beams leads to significant dose enhancements. Clinically, this process occurs in the lung tissue typically surrounding the GTV, resulting in increases to the D98% (PTV). Two isolated tumors with very small PTVs (3 and 6 cm(3)) showed increases in D98% of 23% and 22%. Larger PTVs of 13, 26, and 59 cm(3) had increases of 9%, 6%, and 4%, describing a natural fall-off in enhancement with increasing PTV size. However, three PTVs bounded to the lung wall showed no significant increase, due to lack of dose enhancement in the denser PTV volume. In general, at 0.5 T, the GTV mean dose enhancement is around 60% lower than that at 1 T, while at 3 T, it is 5%-60% higher than 1 T. CONCLUSIONS: Monte Carlo methods have described significant and predictable dose enhancement effects in small lung tumor plans for 6 MV radiotherapy when an external inline magnetic field is included. Results of this study indicate that future clinical inline MRI-guided radiotherapy systems will be able to deliver a dosimetrically superior treatment to small (PTV < 15 cm(3)), isolated lung tumors over non-MRI-Linac systems. This increased efficacy coincides with the reimbursement in the United States of lung CT screening and the likely rapid growth in the number of patients with small lung tumors to be treated with radiotherapy.

MagicPlate-512: A 2D silicon detector array for quality assurance of stereotactic motion adaptive radiotherapy.

PURPOSE: Spatial and temporal resolutions are two of the most important features for quality assurance instrumentation of motion adaptive radiotherapy modalities. The goal of this work is to characterize the performance of the 2D high spatial resolution monolithic silicon diode array named "MagicPlate-512" for quality assurance of stereotactic body radiation therapy (SBRT) and stereotactic radiosurgery (SRS) combined with a dynamic multileaf collimator (MLC) tracking technique for motion compensation. METHODS: MagicPlate-512 is used in combination with the movable platform HexaMotion and a research version of radiofrequency tracking system Calypso driving MLC tracking software. The authors reconstruct 2D dose distributions of small field square beams in three modalities: in static conditions, mimicking the temporal movement pattern of a lung tumor and tracking the moving target while the MLC compensates almost instantaneously for the tumor displacement. Use of Calypso in combination with MagicPlate-512 requires a proper radiofrequency interference shielding. Impact of the shielding on dosimetry has been simulated by (GEANT)4 and verified experimentally. Temporal and spatial resolutions of the dosimetry system allow also for accurate verification of segments of complex stereotactic radiotherapy plans with identification of the instant and location where a certain dose is delivered. This feature allows for retrospective temporal reconstruction of the delivery process and easy identification of error in the tracking or the multileaf collimator driving systems. A sliding MLC wedge combined with the lung motion pattern has been measured. The ability of the MagicPlate-512 (MP512) in 2D dose mapping in all three modes of operation was benchmarked by EBT3 film. RESULTS: Full width at half maximum and penumbra of the moving and stationary dose profiles measured by EBT3 film and MagicPlate-512 confirm that motion has a significant impact on the dose distribution. Motion, no motion, and motion with MLC tracking profiles agreed within 1 and 0.4 mm, respectively, for all field sizes tested. Use of electromagnetic tracking system generates a fluctuation of the detector baseline up to 10% of the full scale signal requiring a proper shielding strategy. MagicPlate-512 is also able to reconstruct the dose variation pulse-by-pulse in each pixel of the detector. An analysis of the dose transients with motion and motion with tracking shows that the tracking feedback algorithm used for this experiment can compensate effectively only the effect of the slower transient components. The fast changing components of the organ motion can contribute only to discrepancy of the order of 15% in penumbral region while the slower components can change the dose profile up to 75% of the expected dose. CONCLUSIONS: MagicPlate-512 is shown to be, potentially, a valid alternative to film or 2D ionizing chambers for quality assurance dosimetry in SRS or SBRT. Its high spatial and temporal resolutions allow for accurate reconstruction of the profile in any conditions with motion and with tracking of the motion. It shows excellent performance to reconstruct the dose deposition in real time or retrospectively as a function of time for detailed analysis of the effect of motion in a specific pixel or area of interest.

Determining appropriate imaging parameters for kilovoltage intrafraction monitoring: an experimental phantom study.

Kilovoltage intrafraction monitoring (KIM) utilises the kV imager during treatment for real-time tracking of prostate fiducial markers. However, its effectiveness relies on sufficient image quality for the fiducial tracking task. To guide the performance characterisation of KIM under different clinically relevant conditions, the effect of different kV parameters and patient size on image quality, and quantification of MV scatter from the patient to the kV detector panel were investigated in this study. Image quality was determined for a range of kV acquisition frame rates, kV exposure, MV dose rates and patient sizes. Two methods were used to determine image quality; the ratio of kV signal through the patient to the MV scatter from the patient incident on the kilovoltage detector, and the signal-to-noise ratio (SNR). The effect of patient size and frame rate on MV scatter was evaluated in a homogeneous CIRS pelvis phantom and marker segmentation was determined utilising the Rando phantom with embedded markers. MV scatter incident on the detector was shown to be dependent on patient thickness and frame rate. The segmentation code was shown to be successful for all frame rates above 3 Hz for the Rando phantom corresponding to a kV to MV ratio of 0.16 and an SNR of 1.67. For a maximum patient dimension less than 36.4 cm the conservative kV parameters of 5 Hz at 1 mAs can be used to reduce dose while retaining image quality, where the current baseline kV parameters of 10 Hz at 1 mAs is shown to be adequate for marker segmentation up to a patient dimension of 40 cm. In conclusion, the MV scatter component of image quality noise for KIM has been quantified. For most prostate patients, use of KIM with 10 Hz imaging at 1 mAs is adequate however image quality can be maintained and imaging dose reduced by altering existing acquisition parameters.

Proton beam deflection in MRI fields: Implications for MRI-guided proton therapy.

PURPOSE: This paper investigates, via magnetic modeling and Monte Carlo simulation, the ability to deliver proton beams to the treatment zone inside a split-bore MRI-guided proton therapy system. METHODS: Field maps from a split-bore 1 T MRI-Linac system are used as input to geant4 Monte Carlo simulations which model the trajectory of proton beams during their paths to the isocenter of the treatment area. Both inline (along the MRI bore) and perpendicular (through the split-bore gap) orientations are simulated. Monoenergetic parallel and diverging beams of energy 90, 195, and 300 MeV starting from 1.5 and 5 m above isocenter are modeled. A phase space file detailing a 2D calibration pattern is used to set the particle starting positions, and their spatial location as they cross isocenter is recorded. No beam scattering, collimation, or modulation of the proton beams is modeled. RESULTS: In the inline orientation, the radial symmetry of the solenoidal style fringe field acts to rotate the protons around the beam's central axis. For protons starting at 1.5 m from isocenter, this rotation is 19° (90 MeV) and 9.8° (300 MeV). A minor focusing toward the beam's central axis is also seen, but only significant, i.e., 2 mm shift at 150 mm off-axis, for 90 MeV protons. For the perpendicular orientation, the main MRI field and near fringe field act as the strongest to deflect the protons in a consistent direction. When starting from 1.5 m above isocenter shifts of 135 mm (90 MeV) and 65 mm (300 MeV) were observed. Further to this, off-axis protons are slightly deflected toward or away from the central axis in the direction perpendicular to the main deflection direction. This leads to a distortion of the phase space pattern, not just a shift. This distortion increases from zero at the central axis to 10 mm (90 MeV) and 5 mm (300 MeV) for a proton 150 mm off-axis. In both orientations, there is a small but subtle difference in the deflection and distortion pattern between protons fired parallel to the beam axis and those fired from a point source. This is indicative of the 3D spatially variant nature of the MRI fringe field. CONCLUSIONS: For the first time, accurate magnetic and Monte Carlo modeling have been used to assess the transport of generic proton beams toward a 1 T split-bore MRI. Significant rotation is observed in the inline orientation, while more complex deflection and distortion are seen in the perpendicular orientation. The results of this study suggest that due to the complexity and energy-dependent nature of the magnetic deflection and distortion, the pencil beam scanning method will be the only choice for delivering a therapeutic proton beam inside a potential MRI-guided proton therapy system in either the inline or perpendicular orientation. Further to this, significant correction strategies will be required to account for the MRI fringe fields.

Radiotherapy beyond cancer: target localization in real-time MRI and treatment planning for cardiac radiosurgery.

PURPOSE: Atrial fibrillation (AFib) is the most common cardiac arrhythmia that affects millions of patients world-wide. AFib is usually treated with minimally invasive, time consuming catheter ablation techniques. While recently noninvasive radiosurgery to the pulmonary vein antrum (PVA) in the left atrium has been proposed for AFib treatment, precise target location during treatment is challenging due to complex respiratory and cardiac motion. A MRI linear accelerator (MRI-Linac) could solve the problems of motion tracking and compensation using real-time image guidance. In this study, the authors quantified target motion ranges on cardiac magnetic resonance imaging (MRI) and analyzed the dosimetric benefits of margin reduction assuming real-time motion compensation was applied. METHODS: For the imaging study, six human subjects underwent real-time cardiac MRI under free breathing. The target motion was analyzed retrospectively using a template matching algorithm. The planning study was conducted on a CT of an AFib patient with a centrally located esophagus undergoing catheter ablation, representing an ideal case for cardiac radiosurgery. The target definition was similar to the ablation lesions at the PVA created during catheter treatment. Safety margins of 0 mm (perfect tracking) to 8 mm (untracked respiratory motion) were added to the target, defining the planning target volume (PTV). For each margin, a 30 Gy single fraction IMRT plan was generated. Additionally, the influence of 1 and 3 T magnetic fields on the treatment beam delivery was simulated using Monte Carlo calculations to determine the dosimetric impact of MRI guidance for two different Linac positions. RESULTS: Real-time cardiac MRI showed mean respiratory target motion of 10.2 mm (superior-inferior), 2.4 mm (anterior-posterior), and 2 mm (left-right). The planning study showed that increasing safety margins to encompass untracked respiratory motion leads to overlapping structures even in the ideal scenario, compromising either normal tissue dose constraints or PTV coverage. The magnetic field caused a slight increase in the PTV dose with the in-line MRI-Linac configuration. CONCLUSIONS: The authors' results indicate that real-time tracking and motion compensation are mandatory for cardiac radiosurgery and MRI-guidance is feasible, opening the possibility of treating cardiac arrhythmia patients completely noninvasively.

Quality assurance for the clinical implementation of kilovoltage intrafraction monitoring for prostate cancer VMAT.

PURPOSE: Kilovoltage intrafraction monitoring (KIM) is a real-time 3D tumor monitoring system for cancer radiotherapy. KIM uses the commonly available gantry-mounted x-ray imager as input, making this method potentially more widely available than dedicated real-time 3D tumor monitoring systems. KIM is being piloted in a clinical trial for prostate cancer patients treated with VMAT (NCT01742403). The purpose of this work was to develop clinical process and quality assurance (QA) practices for the clinical implementation of KIM. METHODS: Informed by and adapting existing guideline documents from other real-time monitoring systems, KIM-specific QA practices were developed. The following five KIM-specific QA tests were included: (1) static localization accuracy, (2) dynamic localization accuracy, (3) treatment interruption accuracy, (4) latency measurement, and (5) clinical conditions accuracy. Tests (1)-(4) were performed using KIM to measure static and representative patient-derived prostate motion trajectories using a 3D programmable motion stage supporting an anthropomorphic phantom with implanted gold markers to represent the clinical treatment scenario. The threshold for system tolerable latency is <1 s. The tolerances for all other tests are that both the mean and standard deviation of the difference between the programmed trajectory and the measured data are <1 mm. The (5) clinical conditions accuracy test compared the KIM measured positions with those measured by kV/megavoltage (MV) triangulation from five treatment fractions acquired in a previous pilot study. RESULTS: For the (1) static localization, (2) dynamic localization, and (3) treatment interruption accuracy tests, the mean and standard deviation of the difference are <1.0 mm. (4) The measured latency is 350 ms. (5) For the tests with previously acquired patient data, the mean and standard deviation of the difference between KIM and kV/MV triangulation are <1.0 mm. CONCLUSIONS: Clinical process and QA practices for the safe clinical implementation of KIM, a novel real-time monitoring system using commonly available equipment, have been developed and implemented for prostate cancer VMAT.

DMLC tracking and gating can improve dose coverage for prostate VMAT.

PURPOSE: To assess and compare the dosimetric impact of dynamic multileaf collimator (DMLC) tracking and gating as motion correction strategies to account for intrafraction motion during conventionally fractionated prostate radiotherapy. METHODS: A dose reconstruction method was used to retrospectively assess the dose distributions delivered without motion correction during volumetric modulated arc therapy fractions for 20 fractions of five prostate cancer patients who received conventionally fractionated radiotherapy. These delivered dose distributions were compared with the dose distributions which would have been delivered had DMLC tracking or gating motion correction strategies been implemented. The delivered dose distributions were constructed by incorporating the observed prostate motion with the patient's original treatment plan to simulate the treatment delivery. The DMLC tracking dose distributions were constructed using the same dose reconstruction method with the addition of MLC positions from Linac log files obtained during DMLC tracking simulations with the observed prostate motions input to the DMLC tracking software. The gating dose distributions were constructed by altering the prostate motion to simulate the application of a gating threshold of 3 mm for 5 s. RESULTS: The delivered dose distributions showed that dosimetric effects of intrafraction prostate motion could be substantial for some fractions, with an estimated dose decrease of more than 19% and 34% from the planned CTVD99% and PTV D95% values, respectively, for one fraction. Evaluation of dose distributions for DMLC tracking and gating deliveries showed that both interventions were effective in improving the CTV D99% for all of the selected fractions to within 4% of planned value for all fractions. For the delivered dose distributions the difference in rectum V65% for the individual fractions from planned ranged from -44% to 101% and for the bladder V65% the range was -61% to 26% from planned. The application of tracking decreased the maximum rectum and bladder V65% difference to 6% and 4%, respectively. CONCLUSIONS: For the first time, the dosimetric impact of DMLC tracking and gating to account for intrafraction motion during prostate radiotherapy has been assessed and compared with no motion correction. Without motion correction intrafraction prostate motion can result in a significant decrease in target dose coverage for a small number of individual fractions. This is unlikely to effect the overall treatment for most patients undergoing conventionally fractionated treatments. Both DMLC tracking and gating demonstrate dose distributions for all assessed fractions that are robust to intrafraction motion.

Electron contamination modeling and reduction in a 1 T open bore inline MRI-linac system.

PURPOSE: A potential side effect of inline MRI-linac systems is electron contamination focusing causing a high skin dose. In this work, the authors reexamine this prediction for an open bore 1 T MRI system being constructed for the Australian MRI-Linac Program. The efficiency of an electron contamination deflector (ECD) in purging electron contamination from the linac head is modeled, as well as the impact of a helium gas region between the deflector and phantom surface for lowering the amount of air-generated contamination. METHODS: Magnetic modeling of the 1 T MRI was used to generate 3D magnetic field maps both with and without the presence of an ECD located immediately below the MLC's. Forty-seven different ECD designs were modeled and for each the magnetic field map was imported into Geant4 Monte Carlo simulations including the linac head, ECD, and a 30 × 30 × 30 cm(3) water phantom located at isocenter. For the first generation system, the x-ray source to isocenter distance (SID) will be 160 cm, resulting in an 81.2 cm long air gap from the base of the ECD to the phantom surface. The first 71.2 cm was modeled as air or helium gas, with the latter encased between two windows of 50 µm thick high density polyethlyene. 2D skin doses (at 70 µm depth) were calculated across the phantom surface at 1 × 1 mm(2) resolution for 6 MV beams of field size of 5 × 5, 10 × 10, and 20 × 20 cm(2). RESULTS: The skin dose was predicted to be of similar magnitude as the generic systems modeled in previous work, 230% to 1400% of D(max) for 5 × 5 to 20 × 20 cm(2), respectively. Inclusion of the ECD introduced a nonuniformity to the MRI imaging field that ranged from ∼20 to ∼140 ppm while the net force acting on the ECD ranged from ∼151 N to ∼1773 N. Various ECD designs were 100% efficient at purging the electron contamination into the ECD magnet banks; however, a small percentage were scattered back into the beam and continued to the phantom surface. Replacing a large portion of the extended air-column between the ECD and phantom surface with helium gas is a key element as it significantly minimized the air-generated contamination. When using an optimal ECD and helium gas region, the 70 µm skin dose is predicted to increase moderately inside a small hot spot over that of the case with no magnetic field present for the jaw defined square beams examined here. These increases include from 12% to 40% of [Formula: see text] for 5 × 5 cm(2), 18% to 55% of D(max) for 10 × 10 cm(2), and from 23% to 65% of D(max) for 20 × 20 cm(2). CONCLUSIONS: Coupling an efficient ECD and helium gas region below the MLCs in the 160 cm isocenter MRI-linac system is predicted to ameliorate the impact electron contamination focusing has on skin dose increases. An ECD is practical as its impact on the MRI imaging distortion is correctable, and the mechanical forces acting on it manageable from an engineering point of view.

Estimation of effective imaging dose for kilovoltage intratreatment monitoring of the prostate position during cancer radiotherapy.

Kilovoltage intratreatment monitoring (KIM) is a novel real-time localization modality where the tumor position is continuously measured during intensity modulated radiation therapy (IMRT) or intensity modulated arc therapy (IMAT) by a kilovoltage (kV) x-ray imager. Adding kV imaging during therapy adds radiation dose. The additional effective dose is quantified for prostate radiotherapy and compared to dose from other localization modalities. The software PCXMC 2.0 was used to calculate the effective dose delivered to a phantom as a function of imager angle and field size for a Varian On-Board Imager. The average angular effective dose was calculated for a field size of 6 cm × 6 cm. The average angular effective dose was used in calculations for different treatment scenarios. Treatment scenarios considered were treatment type and fractionation. For all treatment scenarios, (i.e. conventionally fractionated and stereotactic body radiotherapy (SBRT), IMRT and IMAT), the total KIM dose at 1 Hz ranged from 2-10 mSv. This imaging dose is less than the Navotek radioactive implant dose (64 mSv) and a standard SBRT cone beam computed tomography pretreatment scan dose (22 mSv) over an entire treatment regime. KIM delivers an acceptably low effective dose for daily use as a real-time image-guidance method for prostate radiotherapy.

Electron contamination modeling and skin dose in 6 MV longitudinal field MRIgRT: Impact of the MRI and MRI fringe field.

PURPOSE: In recent times, longitudinal field MRI-linac systems have been proposed for 6 MV MRI-guided radiotherapy (MRIgRT). The magnetic field is parallel with the beam axis and so will alter the transport properties of any electron contamination particles. The purpose of this work is to provide a first investigation into the potential effects of the MR and fringe magnetic fields on the electron contamination as it is transported toward a phantom, in turn, providing an estimate of the expected patient skin dose changes in such a modality. METHODS: Geant4 Monte Carlo simulations of a water phantom exposed to a 6 MV x-ray beam were performed. Longitudinal magnetic fields of strengths between 0 and 3 T were applied to a 30 × 30 × 20 cm(3) phantom. Surrounding the phantom there is a region where the magnetic field is at full MRI strength, consistent with clinical MRI systems. Beyond this the fringe magnetic field entering the collimation system is also modeled. The MRI-coil thickness, fringe field properties, and isocentric distance are varied and investigated. Beam field sizes of 5 × 5, 10 × 10, 15 × 15 and 20 × 20 cm(2) were simulated. Central axis dose, 2D virtual entry skin dose films, and 70 µm skin depth doses were calculated using high resolution scoring voxels. RESULTS: In the presence of a longitudinal magnetic field, electron contamination from the linear accelerator is encouraged to travel almost directly toward the patient surface with minimal lateral spread. This results in a concentration of electron contamination within the x-ray beam outline. This concentration is particularly encouraged if the fringe field encompasses the collimation system. Skin dose increases of up to 1000% were observed for certain configurations and increases above Dmax were common. In nonmagnetically shielded cases, electron contamination generated from the jaw faces and air column is trapped and propagated almost directly to the phantom entry region, giving rise to intense dose hot spots inside the x-ray treatment field. These range up to 1000% or more of Dmax at the CAX, depending on field size, isocenter, and coil thickness. In the case of a fully magnetically shielded collimation system and the lowest MRI field of 0.25 T, the entry skin dose is expected to increase to at least 40%, 50%, 65%, and 80% of Dmax for 5 × 5, 10 × 10, 15 × 15, and 20 × 20 cm(2), respectively. CONCLUSIONS: Electron contamination from the linac head and air column may cause considerable skin dose increases or hot spots at the beam central axis on the entry side of a phantom or patient in longitudinal field 6 MV MRIgRT. This depends heavily on the properties of the magnetic fringe field entering the linac beam collimation system. The skin dose increase is also related to the MRI-coil thickness, the fringe field, and the isocenter distance of the linac. The results of this work indicate that the properties of the MRI fringe field, electron contamination production, and transport must be considered carefully during the design stage of a longitudinal MRI-linac system.

Real-time profiling of respiratory motion: baseline drift, frequency variation and fundamental pattern change.

To precisely ablate tumor in radiation therapy, it is important to locate the tumor position in real time during treatment. However, respiration-induced tumor motions are difficult to track. They are semi-periodic and exhibit variations in baseline, frequency and fundamental pattern (oscillatory amplitude and shape). In this study, we try to decompose the above-mentioned components from discrete observations in real time. Baseline drift, frequency (equivalently phase) variation and fundamental pattern change characterize different aspects of respiratory motion and have distinctive clinical indications. Furthermore, smoothness is a valid assumption for each one of these components in their own spaces, and facilitates effective extrapolation for the purpose of estimation and prediction. We call this process 'profiling' to reflect the integration of information extraction, decomposition, processing and recovery. The proposed method has three major ingredients: (1) real-time baseline and phase estimation based on elliptical shape tracking in augmented state space and Poincaré sectioning principle; (2) estimation of the fundamental pattern by unwarping the observation with phase estimate from the previous step; (3) filtering of individual components and assembly in the original temporal-displacement signal space. We tested the proposed method with both simulated and clinical data. For the purpose of prediction, the results are comparable to what one would expect from a human operator. The proposed approach is fully unsupervised and data driven, making it ideal for applications requiring economy, efficiency and flexibility.

Accuracy in the localization of thoracic and abdominal tumors using respiratory displacement, velocity, and phase.

UNLABELLED: Current four dimensional (4D) computed tomography (CT) reconstruction techniques are retrospectively created based on either the phase or displacement of the respiratory signal. Both techniques have known limitations which cause clinically significant motion artifacts in 4D CT images. These artifacts, which appear as undefined or irregular boundaries in the 4D CT images, cause systematic errors in patient contouring and dose calculations. The aim of this work was to evaluate the reproducibility of tumor position as a function of displacement, phase, and velocity of the respiratory signal, respectively, in order to determine the optimum parameter or combination of parameters to use in order to minimize artifacts in 4D CT images or to accurately deliver radiation to relevant structures during treatment. METHOD AND MATERIALS: Estimated tumor centroid position and respiratory signal data were acquired with the Cyberknife Synchrony system for 26 thoracic radiotherapy patients (52 fractions). A reference respiratory cycle was calculated for each patient. Displacement, phase, and velocity of ten data points were calculated from this reference respiratory cycle, where each point represents an image bin. The corresponding tumor position was then sorted into these image bins if the phase, displacement, simultaneous displacement and phase, or simultaneous displacement and velocity of the respiratory signal were within tolerances of 0.5 mm for displacement and 0.5 mm/s for velocity, respectively, from the corresponding data of the reference cycle for each image bin. RESULTS: The mean of the standard deviations of tumor positions over all bins and all fractions for the superior-inferior direction were 2.13 +/- 1.01 mm for phase sorting, 1.20 +/- 0.76 mm for displacement sorting, 1.20 +/- 0.71 mm for simultaneous displacement and phase sorting, and 1.10 +/- 0.71 mm for simultaneous displacement and velocity sorting, with maximum deviations of 43.0, 16.1, 15.5, and 14.1 mm for each scenario, respectively. The same trend was observed for the anterior-posterior and left-right directions. A linear dependence was observed between the mean of the standard deviations of tumor positions over all fractions as a function of the velocity of the respiratory signal at each bin for all the sorting scenarios. A substantially larger gradient for the phase sorting scenario, compared to the other scenarios, suggests that tumor localization will become increasingly less accurate as the velocity of the tumor increases during a breathing cycle, e.g., if the amplitude of motion increases while the period of the respiratory cycle stays constant or during mid inhale or exhale phases of the respiratory cycle. CONCLUSION: This study illustrates that position of a tumor can be determined more accurately if displacement and velocity are used simultaneously as sorting parameters for 4D CT images or during treatment. A real-time displacement and velocity based 4D CT image sorting method may therefore produce fewer and smaller artifacts in 4D CT images than current retrospective sorting methods.

Prospective displacement and velocity-based cine 4D CT.

Four dimensional (4D) computed tomography (CT) image sorting is currently a retrospective procedure. Mismatches in displacement and/or phase of a patient's respiratory signal, corresponding with two dimensional images taken at subsequent couch positions, become visible as artifacts in reconstructed 4D CT images. These artifacts appear as undefined or irregular boundaries in the 4D CT images and cause systematic errors in patient contouring and dose calculations. In addition, the substantially higher dose required for 4D CT, compared with 3D CT, is of concern. To minimize these problems, we developed a prospective respiratory displacement and velocity based cine 4D CT (PDV CT) method to trigger image acquisition if the displacement and velocity of the respiratory signal occurred within predetermined tolerances simultaneously. The use of velocity avoids real-time phase estimation. Real-time image acquisition ensures data sufficiency, while avoiding the need for redundant data. This may potentially result in a lower dose to the patient. PDV CT was compared with retrospective 4D CT acquisition methods, using respiratory signals of 24 lung cancer patients (103 sessions) under free breathing conditions. Image acquisition was simulated for each of these sessions from the respiratory signal. The root mean square (RMS) of differences between displacements and velocities of the respiratory signal corresponding to subsequent images was calculated in order to evaluate the image-sorting accuracy of each method. Patient dose reductions of 22 to 50% were achieved during image acquisition depending on the model parameters chosen. The mean RMS differences over all sessions and image bins show that PDV CT produces similar results to retrospective displacement sorting overall, although improvements of the RMS difference up to 20% were achieved depending on the model parameters chosen. Velocity RMS differences improved between 30 and 45% when compared with retrospective phase sorting. The efficiency in acquisition compared with retrospective phase sorting varied from approximately 10% for displacement and velocity tolerances of 1 mm and 4 mm/s, respectively, to 80 to 93% for 4 mm and 4 mm/s. The lower variation in the displacement and velocity of the respiratory signal in each image bin indicates that PDV CT could be a valuable tool for reducing artifacts in 4D CT images and lowering patient dose, although the cost may be increased acquisition time.