Thermal noise variance of a receive radiofrequency coil as a respiratory motion sensor.

PURPOSE: Development of a passive respiratory motion sensor based on the noise variance of the receive coil array. METHODS: Respiratory motion alters the body resistance. The noise variance of an RF coil depends on the body resistance and, thus, is also modulated by respiration. For the noise variance monitoring, the noise samples were acquired without and with MR signal excitation on clinical 1.5/3 T MR scanners. The performance of the noise sensor was compared with the respiratory bellow and with the diaphragm displacement visible on MR images. Several breathing patterns were tested. RESULTS: The noise variance demonstrated a periodic, temporal modulation that was synchronized with the respiratory bellow signal. The modulation depth of the noise variance resulting from the respiration varied between the channels of the array and depended on the channel's location with respect to the body. The noise sensor combined with MR acquisition was able to detect the respiratory motion for every k-space read-out line. CONCLUSION: Within clinical MR systems, the respiratory motion can be detected by the noise in receive array. The noise sensor does not require careful positioning unlike the bellow, any additional hardware, and/or MR acquisition. Magn Reson Med 77:221-228, 2017. © 2016 Wiley Periodicals, Inc.

Real-time auto-adaptive margin generation for MLC-tracked radiotherapy.

In radiotherapy, abdominal and thoracic sites are candidates for performing motion tracking. With real-time control it is possible to adjust the multileaf collimator (MLC) position to the target position. However, positions are not perfectly matched and position errors arise from system delays and complicated response of the electromechanic MLC system. Although, it is possible to compensate parts of these errors by using predictors, residual errors remain and need to be compensated to retain target coverage. This work presents a method to statistically describe tracking errors and to automatically derive a patient-specific, per-segment margin to compensate the arising underdosage on-line, i.e. during plan delivery. The statistics of the geometric error between intended and actual machine position are derived using kernel density estimators. Subsequently a margin is calculated on-line according to a selected coverage parameter, which determines the amount of accepted underdosage. The margin is then applied onto the actual segment to accommodate the positioning errors in the enlarged segment. The proof-of-concept was tested in an on-line tracking experiment and showed the ability to recover underdosages for two test cases, increasing [Formula: see text] in the underdosed area about [Formula: see text] and [Formula: see text], respectively. The used dose model was able to predict the loss of dose due to tracking errors and could be used to infer the necessary margins. The implementation had a running time of 23 ms which is compatible with real-time requirements of MLC tracking systems. The auto-adaptivity to machine and patient characteristics makes the technique a generic yet intuitive candidate to avoid underdosages due to MLC tracking errors.

Moderately degenerated lumbar motion segments: Are they truly unstable?

The two main load bearing tissues of the intervertebral disc are the nucleus pulposus and the annulus fibrosus. Both tissues are composed of the same basic components, but differ in their organization and relative amounts. With degeneration, the clear distinction between the two tissues disappears. The changes in biochemical content lead to changes in mechanical behaviour of the intervertebral disc. The aim of the current study was to investigate if well-documented moderate degeneration at the biochemical and fibre structure level leads to instability of the lumbar spine. By taking into account biochemical and ultrastructural changes to the extracellular matrix of degenerating discs, a set of constitutive material parameters were determined that described the individual tissue behaviour. These tissue biomechanical models were then used to simulate dynamic behaviour of the degenerated spinal motion segment, which showed instability in axial rotation, while a stabilizing effect in the other two principle bending directions. When a shear load was applied to the degenerated spinal motion segment, no sign of instability was found. This study found that reported changes to the nucleus pulposus and annulus fibrosus matrix during moderate degeneration lead to a more stable spinal motion segment and that such biomechanical considerations should be incorporated into the general pathophysiological understanding of disc degeneration and how its progress could affect low back pain and its treatments thereof.

Towards real-time thermometry using simultaneous multislice MRI.

MR-guided thermal therapies, such as high-intensity focused ultrasound (MRgHIFU) and laser-induced thermal therapy (MRgLITT) are increasingly being applied in oncology and neurology. MRI is used for guidance since it can measure temperature noninvasively based on the proton resonance frequency shift (PRFS). For therapy guidance using PRFS thermometry, high temporal resolution and large spatial coverage are desirable. We propose to use the parallel imaging technique simultaneous multislice (SMS) in combination with controlled aliasing (CAIPIRINHA) to accelerate the acquisition. We compare this with the sensitivity encoding (SENSE) acceleration technique. Two experiments were performed to validate that SMS can be used to increase the spatial coverage or the temporal resolution. The first was performed in agar gel using LITT heating and a gradient-echo sequence with echo-planar imaging (EPI), and the second was performed in bovine muscle using HIFU heating and a gradient-echo sequence without EPI. In both experiments temperature curves from an unaccelerated scan and from SMS, SENSE, and SENSE/SMS accelerated scans were compared. The precision was quantified by a standard deviation analysis of scans without heating. Both experiments showed a good agreement between the temperature curves obtained from the unaccelerated, and SMS accelerated scans, confirming that accuracy was maintained during SMS acceleration. The standard deviations of the temperature measurements obtained with SMS were significantly smaller than when SENSE was used, implying that SMS allows for higher acceleration. In the LITT and HIFU experiments SMS factors up to 4 and 3 were reached, respectively, with a loss of precision of less than a factor of 3. Based on these results we conclude that SMS acceleration of PRFS thermometry is a valuable addition to SENSE, because it allows for a higher temporal resolution or bigger spatial coverage, with a higher precision.

On-line MR imaging for dose validation of abdominal radiotherapy.

For quality assurance and adaptive radiotherapy, validation of the actual delivered dose is crucial.Intrafractional anatomy changes cannot be captured satisfactorily during treatment with hitherto available imaging modalitites. Consequently, dose calculations are based on the assumption of static anatomy throughout the treatment. However, intra- and interfraction anatomy is dynamic and changes can be significant.In this paper, we investigate the use of an MR-linac as a dose tracking modality for the validation of treatments in abdominal targets where both respiratory and long-term peristaltic and drift motion occur.The on-line MR imaging capability of the modality provides the means to perform respiratory gating of both delivery and acquisition yielding a model-free respiratory motion management under free breathing conditions.In parallel to the treatment, the volumetric patient anatomy was captured and used to calculate the applied dose. Subsequently, the individual doses were warped back to the planning grid to obtain the actual dose accumulated over the entire treatment duration. Ultimately, the planned dose was validated by comparison with the accumulated dose.Representative for a site subject to breathing modulation, two kidney cases (25 Gy target dose) demonstrated the working principle on volunteer data and simulated delivery. The proposed workflow successfully showed its ability to track local dosimetric changes. Integration of the on-line anatomy information could reveal local dose variations -2.3-1.5 Gy in the target volume of a volunteer dataset. In the adjacent organs at risk, high local dose errors ranging from -2.5 to 1.9 Gy could be traced back.

On-line 3D motion estimation using low resolution MRI.

Image processing such as deformable image registration finds its way into radiotherapy as a means to track non-rigid anatomy. With the advent of magnetic resonance imaging (MRI) guided radiotherapy, intrafraction anatomy snapshots become technically feasible. MRI provides the needed tissue signal for high-fidelity image registration. However, acquisitions, especially in 3D, take a considerable amount of time. Pushing towards real-time adaptive radiotherapy, MRI needs to be accelerated without degrading the quality of information. In this paper, we investigate the impact of image resolution on the quality of motion estimations. Potentially, spatially undersampled images yield comparable motion estimations. At the same time, their acquisition times would reduce greatly due to the sparser sampling. In order to substantiate this hypothesis, exemplary 4D datasets of the abdomen were downsampled gradually. Subsequently, spatiotemporal deformations are extracted consistently using the same motion estimation for each downsampled dataset. Errors between the original and the respectively downsampled version of the dataset are then evaluated. Compared to ground-truth, results show high similarity of deformations estimated from downsampled image data. Using a dataset with (2.5 mm)3 voxel size, deformation fields could be recovered well up to a downsampling factor of 2, i.e. (5 mm)3. In a therapy guidance scenario MRI, imaging speed could accordingly increase approximately fourfold, with acceptable loss of estimated motion quality.

On the suitability of Elekta's Agility 160 MLC for tracked radiation delivery: closed-loop machine performance.

For motion adaptive radiotherapy, dynamic multileaf collimator tracking can be employed to reduce treatment margins by steering the beam according to the organ motion. The Elekta Agility 160 MLC has hitherto not been evaluated for its tracking suitability. Both dosimetric performance and latency are key figures and need to be assessed generically, independent of the used motion sensor. In this paper, we propose the use of harmonic functions directly fed to the MLC to determine its latency during continuous motion. Furthermore, a control variable is extracted from a camera system and fed to the MLC. Using this setup, film dosimetry and subsequent γ statistics are performed, evaluating the response when tracking (MRI)-based physiologic motion in a closed-loop. The delay attributed to the MLC itself was shown to be a minor contributor to the overall feedback chain as compared to the impact of imaging components such as MRI sequences. Delay showed a linear phase behaviour of the MLC employed in continuously dynamic applications, which enables a general MLC-characterization. Using the exemplary feedback chain, dosimetry showed a vast increase in pass rate employing γ statistics. In this early stage, the tracking performance of the Agility using the test bench yielded promising results, making the technique eligible for translation to tracking using clinical imaging modalities.

3D geometry analysis of the medial meniscus--a statistical shape modeling approach.

The geometry-dependent functioning of the meniscus indicates that detailed knowledge on 3D meniscus geometry and its inter-subject variation is essential to design well functioning anatomically shaped meniscus replacements. Therefore, the aim of this study was to quantify 3D meniscus geometry and to determine whether variation in medial meniscus geometry is size- or shape-driven. Also we performed a cluster analysis to identify distinct morphological groups of medial menisci and assessed whether meniscal geometry is gender-dependent. A statistical shape model was created, containing the meniscus geometries of 35 subjects (20 females, 15 males) that were obtained from MR images. A principal component analysis was performed to determine the most important modes of geometry variation and the characteristic changes per principal component were evaluated. Each meniscus from the original dataset was then reconstructed as a linear combination of principal components. This allowed the comparison of male and female menisci, and a cluster analysis to determine distinct morphological meniscus groups. Of the variation in medial meniscus geometry, 53.8% was found to be due to primarily size-related differences and 29.6% due to shape differences. Shape changes were most prominent in the cross-sectional plane, rather than in the transverse plane. Significant differences between male and female menisci were only found for principal component 1, which predominantly reflected size differences. The cluster analysis resulted in four clusters, yet these clusters represented two statistically different meniscal shapes, as differences between cluster 1, 2 and 4 were only present for principal component 1. This study illustrates that differences in meniscal geometry cannot be explained by scaling only, but that different meniscal shapes can be distinguished. Functional analysis, e.g. through finite element modeling, is required to assess whether these distinct shapes actually influence the biomechanical performance of the meniscus.

Proof of concept of MRI-guided tracked radiation delivery: tracking one-dimensional motion.

In radiotherapy one aims to deliver a radiation dose to a tumour with high geometrical accuracy while sparing organs at risk (OARs). Although image guidance decreases geometrical uncertainties, treatment of cancer of abdominal organs is further complicated by respiratory motion, requiring intra-fraction motion compensation to fulfil the treatment intent. With an ideal delivery system, the optimal method of intra-fraction motion compensation is to adapt the beam collimation to the moving target using a dynamic multi-leaf collimator (MLC) aperture. The many guidance strategies for such tracked radiation delivery tested up to now mainly use markers and are therefore invasive and cannot deal with target deformations or adaptations for OAR positions. We propose to address these shortcomings using the online MRI guidance provided by an MRI accelerator and present a first step towards demonstration of the technical feasibility of this proposal. The position of a phantom subjected to one-dimensional (1D) periodic translation was tracked using a fast 1D MR sequence. Real-time communication with the MR scanner and control of the MLC aperture were established. Based on the time-resolved position of the phantom, tracked radiation delivery to the phantom was realized. Dose distributions for various delivery conditions were recorded on a gafchromic film. Without motion a sharply defined dose distribution is obtained, whereas considerable blur occurs for delivery to a moving phantom. With compensation for motion, the sharpness of the dose distribution is nearly restored. The total latency in our motion management architecture is approximately 200 ms. Combination of the recorded phantom and aperture positions with the planned dose distribution enabled the reconstruction of the delivered dose in all cases, which illustrates the promise of online dose accumulation and confirms that latency compensation could further enhance our results. For a simple 1D tracked delivery scenario, the technical feasibility of MRI-guided tracked radiation delivery is confirmed. More generic tracking scenarios require advanced MRI, leading to increased acquisition time and more challenging image processing problems. Latency compensation is therefore an important subject of future investigations.

Towards inherently distortion-free MR images for image-guided radiotherapy on an MRI accelerator.

In MR-guided interventions, it is mandatory to establish a solid relationship between the imaging coordinate system and world coordinates. This is particularly important in image-guided radiotherapy (IGRT) on an MRI accelerator, as the interaction of matter with γ-radiation cannot be visualized. In conventional acquisitions, off-resonance effects cause discrepancies between coordinate systems. We propose to mitigate this by using only phase encoding and to reduce the longer acquisitions by under-sampling and regularized reconstruction. To illustrate the performance of this acquisition in the presence of off-resonance phenomena, phantom and in vivo images are acquired using spin-echo (SE) and purely phase-encoded sequences. Data are retrospectively under-sampled and reconstructed iteratively. We observe accurate geometries in purely phase-encoded images for all cases, whereas SE images of the same phantoms display image distortions. Regularized reconstruction yields accurate phantom images under high acceleration factors. In vivo images were reconstructed faithfully while using acceleration factors up to 4. With the proposed technique, inherently undistorted images with one-to-one correspondence to world coordinates can be obtained. It is a valuable tool in geometry quality assurance, treatment planning and online image guidance. Under-sampled acquisition combined with regularized reconstruction can be used to accelerate the acquisition while retaining geometrical accuracy.

Integrated megavoltage portal imaging with a 1.5 T MRI linac.

In this note, the feasibility of complementing our hybrid 1.5 T MRI linac (MRL) with a megavoltage (MV) portal imager is investigated. A standard aSi MV detector panel is added to the system and both qualitative and quantitative performances are determined. Simultaneous MR imaging and transmission imaging can be performed without mutual interference. The MV image quality is compromised by beam transmission and longer isocentre distance; still, the field edges and bony anatomy can be detected at very low dose levels of 0.4 cGy. MV imaging integrated with the MRL provides an independent and well-established position verification tool, a field edge check and a calibration for alignment of the coordinate systems of the MRI and the accelerator. The portal imager can also be a valuable means for benchmarking MRI-guided position verification protocols on a patient-specific basis in the introductory phase.

Towards MRI-guided linear accelerator control: gating on an MRI accelerator.

To boost the possibilities of image guidance in radiotherapy by providing images with superior soft-tissue contrast during treatment, we pursue diagnostic quality MRI functionality integrated with a linear accelerator. Large respiration-induced semi-periodic target excursions hamper treatment of cancer of the abdominal organs. Methods to compensate in real time for such motion are gating and tracking. These strategies are most effective in cases where anatomic motion can be visualized directly, which supports the use of an integrated MRI accelerator. We establish here an infrastructure needed to realize gated radiation delivery based on MR feedback and demonstrate its potential as a first step towards more advanced image guidance techniques. The position of a phantom subjected to one-dimensional periodic translation is tracked with the MR scanner. Real-time communication with the MR scanner and control of the radiation beam are established. Based on the time-resolved position of the phantom, gated radiation delivery to the phantom is realized. Dose distributions for dynamic delivery conditions with varying gating windows are recorded on gafchromic film. The similarity between dynamically and statically obtained dose profiles gradually increases as the gating window is decreased. With gating windows of 5 mm, we obtain sharp dose profiles. We validate our gating implementation by comparing measured dose profiles to theoretical profiles calculated using the knowledge of the imposed motion pattern. Excellent correspondence is observed. At the same time, we show that real-time on-line reconstruction of the accumulated dose can be performed using time-resolved target position information. This facilitates plan adaptation not only on a fraction-to-fraction scale but also during one fraction, which is especially valuable in highly accelerated treatment strategies. With the currently established framework and upcoming improvements to our prototype-integrated MRI accelerator, we will realize more intricate MRI-guided linear accelerator control in the near future.

Real-time correction of magnetic field inhomogeneity-induced image distortions for MRI-guided conventional and proton radiotherapy.

Image-guided radiotherapy has the potential to increase the success of treatment by decreasing uncertainties concerning tumour position and shape. We pursue integrated diagnostic quality MRI functionality with radiotherapy systems to boost the possibilities of image guidance by providing images with superior soft-tissue contrast during treatment. However, the use of MR images in radiotherapy can be hindered by geometrical distortions due to magnetic field inhomogeneity problems. A method for fast correction of these distortions is presented and implemented. Using a 20 cm square phantom containing a regular grid, a measure of residual deformation after correction is established. At very low gradient strength (which leads to large deformations) a maximum displacement of 2.9 mm is shown to be reduced to 0.63 mm. Next, the method is applied in vivo to the case of pelvic body contour extraction for prostate radiotherapy treatment planning. Here, again with low gradient strengths, distortions of up to 6 mm can be reduced to 2 mm. All results are provided within a lag time of 8 ms. We discuss implications of image distortions for MRI-guided photon and proton radiotherapy separately, since the dose-depth curves in these treatments are very different. We argue that, although field inhomogeneities cannot be prevented from occurring, distortion correction is not always necessary in practice. This work opens new possibilities for investigating on-line MRI-based plan adaptations and ultimately MRI-based treatment planning.

Integrating a 1.5 T MRI scanner with a 6 MV accelerator: proof of concept.

At the UMC Utrecht, The Netherlands, we have constructed a prototype MRI accelerator. The prototype is a modified 6 MV Elekta (Crawley, UK) accelerator next to a modified 1.5 T Philips Achieva (Best, The Netherlands) MRI system. From the initial design onwards, modifications to both systems were aimed to yield simultaneous and unhampered operation of the MRI and the accelerator. Indeed, the simultaneous operation is shown by performing diagnostic quality 1.5 T MRI with the radiation beam on. No degradation of the performance of either system was found. The integrated 1.5 T MRI system and radiotherapy accelerator allow simultaneous irradiation and MR imaging. The full diagnostic imaging capacities of the MRI can be used; dedicated sequences for MRI-guided radiotherapy treatments will be developed. This proof of concept opens the door towards a clinical prototype to start testing MRI-guided radiation therapy (MRIgRT) in the clinic.