R.H.N. Tijssen

MRI-based tumor motion characterization and gating schemes for radiation therapy of pancreatic cancer

BACKGROUND AND PURPOSE To characterize pancreatic tumor motion and to develop a gating scheme for radiotherapy in pancreatic cancer. MATERIALS AND METHODS Two cine MRIs of 60s each were performed in fifteen pancreatic cancer patients, one in sagittal direction and one in coronal direction. A Minimum Output Sum of Squared Error (MOSSE) adaptive correlation filter was used to quantify tumor motion in craniocaudal, lateral and anteroposterior directions. To develop a gating scheme, stability of the breathing phases was examined and a gating window assessment was created, incorporating tumor motion, treatment time and motion margins. RESULTS The largest tumor motion was found in craniocaudal direction, with an average peak-to-peak amplitude of 15mm (range 6-34mm). Amplitude of the tumor in the anteroposterior direction was on average 5mm (range 1-13mm). The least motion was seen in lateral direction (average 3mm, range 2-5mm). The end exhale position was the most stable position in the breathing cycle and tumors spent more time closer to the end exhale position than to the end inhale position. On average, a margin of 25% of the maximum craniocaudal breathing amplitude was needed to achieve full target coverage with a duty cycle of 50%. When reducing the duty cycle to 50%, a margin of 5mm was sufficient to cover the target in 11 out of 15 patients. CONCLUSION Gated delivery for radiotherapy of pancreatic cancer is best performed around the end exhale position as this is the most stable position in the breathing cycle. Considerable margin reduction can be established at moderate duty cycles, yielding acceptable treatment efficiency. However, motion patterns and amplitude do substantially differ between individual patients. Therefore, individual treatment strategies should be considered for radiotherapy in pancreatic cancer.

Optimizing 4-Dimensional Magnetic Resonance Imaging Data Sampling for Respiratory Motion Analysis of Pancreatic Tumors

PURPOSE To determine the optimum sampling strategy for retrospective reconstruction of 4-dimensional (4D) MR data for nonrigid motion characterization of tumor and organs at risk for radiation therapy purposes. METHODS AND MATERIALS For optimization, we compared 2 surrogate signals (external respiratory bellows and internal MRI navigators) and 2 MR sampling strategies (Cartesian and radial) in terms of image quality and robustness. Using the optimized protocol, 6 pancreatic cancer patients were scanned to calculate the 4D motion. Region of interest analysis was performed to characterize the respiratory-induced motion of the tumor and organs at risk simultaneously. RESULTS The MRI navigator was found to be a more reliable surrogate for pancreatic motion than the respiratory bellows signal. Radial sampling is most benign for undersampling artifacts and intraview motion. Motion characterization revealed interorgan and interpatient variation, as well as heterogeneity within the tumor. CONCLUSIONS A robust 4D-MRI method, based on clinically available protocols, is presented and successfully applied to characterize the abdominal motion in a small number of pancreatic cancer patients.

Image-driven, model-based 3D abdominal motion estimation for MR-guided radiotherapy.

Respiratory motion introduces substantial uncertainties in abdominal radiotherapy for which traditionally large margins are used. The MR-Linac will open up the opportunity to acquire high resolution MR images just prior to radiation and during treatment. However, volumetric MRI time series are not able to characterize 3D tumor and organ-at-risk motion with sufficient temporal resolution. In this study we propose a method to estimate 3D deformation vector fields (DVFs) with high spatial and temporal resolution based on fast 2D imaging and a subject-specific motion model based on respiratory correlated MRI. In a pre-beam phase, a retrospectively sorted 4D-MRI is acquired, from which the motion is parameterized using a principal component analysis. This motion model is used in combination with fast 2D cine-MR images, which are acquired during radiation, to generate full field-of-view 3D DVFs with a temporal resolution of 476 ms. The geometrical accuracies of the input data (4D-MRI and 2D multi-slice acquisitions) and the fitting procedure were determined using an MR-compatible motion phantom and found to be 1.0-1.5 mm on average. The framework was tested on seven healthy volunteers for both the pancreas and the kidney. The calculated motion was independently validated using one of the 2D slices, with an average error of 1.45 mm. The calculated 3D DVFs can be used retrospectively for treatment simulations, plan evaluations, or to determine the accumulated dose for both the tumor and organs-at-risk on a subject-specific basis in MR-guided radiotherapy.

4D MR imaging using robust internal respiratory signal.

The purpose of this study is to investigate the feasibility of using internal respiratory (IR) surrogates to sort four-dimensional (4D) magnetic resonance (MR) images. The 4D MR images were constructed by acquiring fast 2D cine MR images sequentially, with each slice scanned for more than one breathing cycle. The 4D volume was then sorted retrospectively using the IR signal. In this study, we propose to use multiple low-frequency components in the Fourier space as well as the anterior body boundary as potential IR surrogates. From these potential IR surrogates, we used a clustering algorithm to identify those that best represented the respiratory pattern to derive the IR signal. A study with healthy volunteers was performed to assess the feasibility of the proposed IR signal. We compared this proposed IR signal with the respiratory signal obtained using respiratory bellows. Overall, 99% of the IR signals matched the bellows signals. The average difference between the end inspiration times in the IR signal and bellows signal was 0.18 s in this cohort of matching signals. For the acquired images corresponding to the other 1% of non-matching signal pairs, the respiratory motion shown in the images was coherent with the respiratory phases determined by the IR signal, but not the bellows signal. This suggested that the IR signal determined by the proposed method could potentially correct the faulty bellows signal. The sorted 4D images showed minimal mismatched artefacts and potential clinical applicability. The proposed IR signal therefore provides a feasible alternative to effectively sort MR images in 4D.

An optimization framework to maximize signal-to-noise ratio in simultaneous multi-slice body imaging

Parallel imaging is essential for the acceleration of abdominal and pelvic 2D multi‐slice imaging, in order to reduce scan time and mitigate motion artifacts. Controlled Aliasing In Parallel Imaging Results IN Higher Acceleration (CAIPIRINHA) accelerated imaging has been shown to increase the signal‐to‐noise ratio (SNR) significantly compared with in‐plane parallel imaging with similar acceleration. We hypothesize that for CAIPIRINHA‐accelerated abdominal imaging the consistency of image quality and SNR is more difficult to achieve due to the subject‐specific coil sensitivity profiles, caused by (1) flexible coil placement; (2) variations in anatomy; and (3) variations in scan coverage along the superior–inferior direction. To test this, a mathematical framework is introduced that calculates the (retained) SNR for in‐plane and simultaneous multi‐slice (SMS)‐accelerated acquisitions. Moreover, this framework was used to optimize the sampling pattern by maximizing the local SNR within a region of interest (ROI) through non‐linear, RF‐induced CAIPIRINHA slice shifts. The framework was evaluated on 14 healthy subjects and the optimized sampling pattern was compared with in‐plane acceleration and CAIPIRINHA acceleration with linear slice shifts, which are primarily used in brain imaging. We demonstrate that the field of view (FOV) in the superior–inferior direction, the coil positioning and the individual anatomy have a large impact on the image SNR (changes up to 50% for varying coil positions and 40% differences between subjects) and image artifacts for simultaneous multi‐slice acceleration. Consequently, sampling patterns have to be optimized for acquisitions employing different FOVs and ideally on an individual basis. Optimization of the sampling pattern, which exploits non‐linear shifts between slices, showed a considerable SNR increase (10–30%) for higher acceleration factors. The framework outlined in this article can be used to optimize sampling patterns for a broad range of accelerated body acquisitions on an individual basis. Copyright

Effect of intra-fraction motion on the accumulated dose for free-breathing MR-guided stereotactic body radiation therapy of renal-cell carcinoma

Stereotactic body radiation therapy (SBRT) has shown great promise in increasing local control rates for renal-cell carcinoma (RCC). Characterized by steep dose gradients and high fraction doses, these hypo-fractionated treatments are, however, prone to dosimetric errors as a result of variations in intra-fraction respiratory-induced motion, such as drifts and amplitude alterations. This may lead to significant variations in the deposited dose. This study aims to develop a method for calculating the accumulated dose for MRI-guided SBRT of RCC in the presence of intra-fraction respiratory variations and determine the effect of such variations on the deposited dose. For this, RCC SBRT treatments were simulated while the underlying anatomy was moving, based on motion information from three motion models with increasing complexity: (1) STATIC, in which static anatomy was assumed, (2) AVG-RESP, in which 4D-MRI phase-volumes were time-weighted, and (3) PCA, a method that generates 3D volumes with sufficient spatio-temporal resolution to capture respiration and intra-fraction variations. Five RCC patients and two volunteers were included and treatments delivery was simulated, using motion derived from subject-specific MR imaging. Motion was most accurately estimated using the PCA method with root-mean-squared errors of 2.7, 2.4, 1.0 mm for STATIC, AVG-RESP and PCA, respectively. The heterogeneous patient group demonstrated relatively large dosimetric differences between the STATIC and AVG-RESP, and the PCA reconstructed dose maps, with hotspots up to [Formula: see text] of the D99 and an underdosed GTV in three out of the five patients. This shows the potential importance of including intra-fraction motion variations in dose calculations.

Intrafraction motion quantification and planning target volume margin determination of head-and-neck tumors using cine magnetic resonance imaging

PURPOSE To quantify intrafractional motion to determine population-based radiotherapy treatment margins for head-and-neck tumors. METHODS Cine MR imaging was performed in 100 patients with head-and-neck cancer on a 3T scanner in a radiotherapy treatment setup. MR images were analyzed using deformable image registration (optical flow algorithm) and changes in tumor contour position were used to calculate the tumor motion. The tumor motion was used together with patient setup errors (450 patients) to calculate population-based PTV margins. RESULTS Tumor motion was quantified in 84 patients (12/43/29 nasopharynx/oropharynx/larynx, 16 excluded). The mean maximum (95th percentile) tumor motion (swallowing excluded) was: 2.3 mm in superior, 2.4 mm in inferior, 1.8 mm in anterior and 1.7 mm in posterior direction. PTV margins were: 2.8 mm isotropic for nasopharyngeal tumors, 3.2 mm isotropic for oropharyngeal tumors and 4.3 mm in inferior-superior and 3.2 mm in anterior-posterior for laryngeal tumors, for our institution. CONCLUSIONS Intrafractional head-and-neck tumor motion was quantified and population-based PTV margins were calculated. Although the average tumor motion was small (95th percentile motion <3.0 mm), tumor motion varied considerably between patients (0.1-12.0 mm). The intrafraction motion expanded the CTV-to-PTV with 1.7 mm for laryngeal tumors, 0.6 mm for oropharyngeal tumors and 0.2 mm for nasopharyngeal tumors.

Characterization of the first RF coil dedicated to 1.5 T MR guided radiotherapy

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.

Nuts and bolts of 4D-MRI for radiotherapy

Magnetic resonance imaging (MRI) is increasingly being used in the radiotherapy workflow because of its superior soft tissue contrast and high flexibility in contrast. In addition to anatomical and functional imaging, MRI can also be used to characterize the physiologically induced motion of both the tumor and organs-at-risk. Respiratory-correlated 4D-MRI has gained large interest as an alternative to 4D-CT for the characterization of respiratory motion throughout the thorax and abdomen. These 4D-MRI data sets consist of three spatial dimensions and the respiratory phase or amplitude over the fourth dimension (opposed to time-resolved 4D-MRI that represents time in the fourth dimension). Over the last 15 years numerous methods have been presented in literature. This review article provides a comprehensive overview of the various 4D-MRI techniques, and describes the differences in MRI data acquisition and 4D data set generation from a methodological point of view. The current status and future perspective of these techniques are highlighted, and the requirements for safe introduction into the clinic (e.g. method validation) are discussed.

Characterization of imaging latency for real-time MRI-guided radiotherapy

Hybrid MR-linac systems can use fast dynamic MR sequences for tumor tracking and adapt the radiation treatment in real-time. For this the imaging latency must be as short as possible. This work describes how different acquisition parameters influence this latency. First, the latency was measured for Cartesian readouts with phase encode orderings linear, reverse-linear, and high-low. Second, the latency was measured for radial readouts with linear and golden angle profile orderings. To reduce the latency, a spatio-temporal (k-t) filter that suppresses the k-space center of earlier acquired spokes was implemented for the golden angle sequence. For Cartesian readouts a high-low ordering achieved a three times lower latency compared to a linear ordering with our sampling parameters. For radial readouts the filter was able to reduce the acquisition latency from half the acquisition time to a quarter of the acquisition time. The filter did not compromise the signal-to-noise ratio and the artifact power.