Andreas Wetscherek

MRI-guided lung SBRT: Present and future developments

Stereotactic body radiotherapy (SBRT) is rapidly becoming an alternative to surgery for the treatment of early-stage non-small cell lung cancer patients. Lung SBRT is administered in a hypo-fractionated, conformal manner, delivering high doses to the target. To avoid normal-tissue toxicity, it is crucial to limit the exposure of nearby healthy organs-at-risk (OAR). Current image-guided radiotherapy strategies for lung SBRT are mostly based on X-ray imaging modalities. Although still in its infancy, magnetic resonance imaging (MRI) guidance for lung SBRT is not exposure-limited and MRI promises to improve crucial soft-tissue contrast. Looking beyond anatomical imaging, functional MRI is expected to inform treatment decisions and adaptations in the future. This review summarises and discusses how MRI could be advantageous to the different links of the radiotherapy treatment chain for lung SBRT: diagnosis and staging, tumour and OAR delineation, treatment planning, and inter- or intrafractional motion management. Special emphasis is placed on a new generation of hybrid MRI treatment devices and their potential for real-time adaptive radiotherapy.

Respiratory motion compensation for simultaneous PET/MR based on highly undersampled MR data

PURPOSE Positron emission tomography (PET) of the thorax region is impaired by respiratory patient motion. To account for motion, the authors propose a new method for PET/magnetic resonance (MR) respiratory motion compensation (MoCo), which uses highly undersampled MR data with acquisition times as short as 1 min/bed. METHODS The proposed PET/MR MoCo method (4D jMoCo PET) uses radial MR data to estimate the respiratory patient motion employing MR joint motion estimation and image reconstruction with temporal median filtering. Resulting motion vector fields are incorporated into the system matrix of the PET reconstruction. The proposed approach is evaluated for the thorax region utilizing a PET/MR simulation with 1 min MR acquisition time and simultaneous PET/MR measurements of six patients with MR acquisition times of 1 and 5 min and radial undersampling factors of 11.2 and 2.2, respectively. Reconstruction results are compared to 3D PET, 4D gated PET and a standard MoCo method (4D sMoCo PET), which performs iterative image reconstruction and motion estimation sequentially. Quantitative analysis comprises the parameters SUVmean, SUVmax, full width at half-maximum/lesion volume, contrast and signal-to-noise ratio. RESULTS For simulated PET data, our quantitative analysis shows that the proposed 4D jMoCo PET approach with temporal filtering achieves the best quantification accuracy of all tested reconstruction methods with a mean absolute deviation of 2.3% when compared to the ground truth. For measured PET patient data, the mean absolute deviation of 4D jMoCo PET using a 1 min MR acquisition for motion estimation is 2.1% relative to the 5 min MR acquisition. This demonstrates a robust behavior even in case of strong undersampling at MR acquisition times as short as 1 min. In contrast, 4D sMoCo PET shows considerable reduction of quantification accuracy for the 1 min MR acquisition time. Relative to 3D PET, the proposed 4D jMoCo PET approach with temporal filtering yields an average increase of SUVmean, SUVmax, and contrast of 29.9% and 13.8% for simulated and measured PET data, respectively. CONCLUSIONS Employing artifact-robust motion estimation enables PET/MR respiratory MoCo with MR acquisition times as short as 1 min/bed improving PET image quality and quantification accuracy.

T2-Weighted 4D Magnetic Resonance Imaging for Application in Magnetic Resonance-Guided Radiotherapy Treatment Planning

Objectives The aim of this study was to develop and verify a method to obtain good temporal resolution T2-weighted 4-dimensional (4D-T2w) magnetic resonance imaging (MRI) by using motion information from T1-weighted 4D (4D-T1w) MRI, to support treatment planning in MR-guided radiotherapy. Materials and Methods Ten patients with primary non–small cell lung cancer were scanned at 1.5 T axially with a volumetric T2-weighted turbo spin echo sequence gated to exhalation and a volumetric T1-weighted stack-of-stars spoiled gradient echo sequence with golden angle spacing acquired in free breathing. From the latter, 20 respiratory phases were reconstructed using the recently developed 4D joint MoCo-HDTV algorithm based on the self-gating signal obtained from the k-space center. Motion vector fields describing the respiratory cycle were obtained by deformable image registration between the respiratory phases and projected onto the T2-weighted image volume. The resulting 4D-T2w volumes were verified against the 4D-T1w volumes: an edge-detection method was used to measure the diaphragm positions; the locations of anatomical landmarks delineated by a radiation oncologist were compared and normalized mutual information was calculated to evaluate volumetric image similarity. Results High-resolution 4D-T2w MRI was obtained. Respiratory motion was preserved on calculated 4D-T2w MRI, with median diaphragm positions being consistent with less than 6.6 mm (2 voxels) for all patients and less than 3.3 mm (1 voxel) for 9 of 10 patients. Geometrical positions were coherent between 4D-T1w and 4D-T2w MRI as Euclidean distances between all corresponding anatomical landmarks agreed to within 7.6 mm (Euclidean distance of 2 voxels) and were below 3.8 mm (Euclidean distance of 1 voxel) for 355 of 470 pairs of anatomical landmarks. Volumetric image similarity was commensurate between 4D-T1w and 4D-T2w MRI, as mean percentage differences in normalized mutual information (calculated over all respiratory phases and patients), between corresponding respiratory phases of 4D-T1w and 4D-T2w MRI and the tie-phase of 4D-T1w and 3-dimensional T2w MRI, were consistent to 0.41% ± 0.37%. Four-dimensional T2w MRI displayed tumor extent, structure, and position more clearly than corresponding 4D-T1w MRI, especially when mobile tumor sites were adjacent to organs at risk. Conclusions A methodology to obtain 4D-T2w MRI that retrospectively applies the motion information from 4D-T1w MRI to 3-dimensional T2w MRI was developed and verified. Four-dimensional T2w MRI can assist clinicians in delineating mobile lesions that are difficult to define on 4D-T1w MRI, because of poor tumor-tissue contrast.

Regional Lung Ventilation Analysis Using Temporally Resolved Magnetic Resonance Imaging.

Objective We propose a computer-aided method for regional ventilation analysis and observation of lung diseases in temporally resolved magnetic resonance imaging (4D MRI). Methods A shape model–based segmentation and registration workflow was used to create an atlas-derived reference system in which regional tissue motion can be quantified and multimodal image data can be compared regionally. Model-based temporal registration of the lung surfaces in 4D MRI data was compared with the registration of 4D computed tomography (CT) images. A ventilation analysis was performed on 4D MR images of patients with lung fibrosis; 4D MR ventilation maps were compared with corresponding diagnostic 3D CT images of the patients and 4D CT maps of subjects without impaired lung function (serving as reference). Results Comparison between the computed patient-specific 4D MR regional ventilation maps and diagnostic CT images shows good correlation in conspicuous regions. Comparison to 4D CT-derived ventilation maps supports the plausibility of the 4D MR maps. Dynamic MRI-based flow-volume loops and spirograms further visualize the free-breathing behavior. Conclusions The proposed methods allow for 4D MR-based regional analysis of tissue dynamics and ventilation in spontaneous breathing and comparison of patient data. The proposed atlas-based reference coordinate system provides an automated manner of annotating and comparing multimodal lung image data.

Characterization of the diffusion coefficient of blood

To characterize the diffusion coefficient of human blood for accurate results in intravoxel incoherent motion imaging.

Tumour auto-contouring on 2d cine MRI for locally advanced lung cancer: A comparative study

Background and purpose Radiotherapy guidance based on magnetic resonance imaging (MRI) is currently becoming a clinical reality. Fast 2d cine MRI sequences are expected to increase the precision of radiation delivery by facilitating tumour delineation during treatment. This study compares four auto-contouring algorithms for the task of delineating the primary tumour in six locally advanced (LA) lung cancer patients. Material and methods Twenty-two cine MRI sequences were acquired using either a balanced steady-state free precession or a spoiled gradient echo imaging technique. Contours derived by the auto-contouring algorithms were compared against manual reference contours. A selection of eight image data sets was also used to assess the inter-observer delineation uncertainty. Results Algorithmically derived contours agreed well with the manual reference contours (median Dice similarity index: ⩾0.91). Multi-template matching and deformable image registration performed significantly better than feature-driven registration and the pulse-coupled neural network (PCNN). Neither MRI sequence nor image orientation was a conclusive predictor for algorithmic performance. Motion significantly degraded the performance of the PCNN. The inter-observer variability was of the same order of magnitude as the algorithmic performance. Conclusion Auto-contouring of tumours on cine MRI is feasible in LA lung cancer patients. Despite large variations in implementation complexity, the different algorithms all have relatively similar performance.

Super-resolution T2-weighted 4D MRI for image guided radiotherapy

Background and purpose The superior soft-tissue contrast of 4D-T2w MRI motivates its use for delineation in radiotherapy treatment planning. We address current limitations of slice-selective implementations, including thick slices and artefacts originating from data incompleteness and variable breathing. Materials and methods A method was developed to calculate midposition and 4D-T2w images of the whole thorax from continuously acquired axial and sagittal 2D-T2w MRI (1.5 × 1.5 × 5.0 mm3). The method employed image-derived respiratory surrogates, deformable image registration and super-resolution reconstruction. Volunteer imaging and a respiratory motion phantom were used for validation. The minimum number of dynamic acquisitions needed to calculate a representative midposition image was investigated by retrospectively subsampling the data (10–30 dynamic acquisitions). Results Super-resolution 4D-T2w MRI (1.0 × 1.0 × 1.0 mm3, 8 respiratory phases) did not suffer from data incompleteness and exhibited reduced stitching artefacts compared to sorted multi-slice MRI. Experiments using a respiratory motion phantom and colour-intensity projection images demonstrated a minor underestimation of the motion range. Midposition diaphragm differences in retrospectively subsampled acquisitions were <1.1 mm compared to the full dataset. 10 dynamic acquisitions were found sufficient to generate midposition MRI. Conclusions A motion-modelling and super-resolution method was developed to calculate high quality 4D/midposition T2w MRI from orthogonal 2D-T2w MRI.

The impact of 2D cine MR imaging parameters on automated tumor and organ localization for MR-guided real-time adaptive radiotherapy

Abstract 2D cine MR imaging may be utilized to monitor rapidly moving tumors and organs-at-risk for real-time adaptive radiotherapy. This study systematically investigates the impact of geometric imaging parameters on the ability of 2D cine MR imaging to guide template-matching-driven autocontouring of lung tumors and abdominal organs. Abdominal 4D MR images were acquired of six healthy volunteers and thoracic 4D MR images were obtained of eight lung cancer patients. At each breathing phase of the images, the left kidney and gallbladder or lung tumor, respectively, were outlined as volumes of interest. These images and contours were used to create artificial 2D cine MR images, while simultaneously serving as 3D ground truth. We explored the impact of five different imaging parameters (pixel size, slice thickness, imaging plane orientation, number and relative alignment of images as well as strategies to create training images). For each possible combination of imaging parameters, we generated artificial 2D cine MR images as training and test images. A template-matching algorithm used the training images to determine the tumor or organ position in the test images. Subsequently, a 3D base contour was shifted to the determined position and compared to the ground truth via centroid distance and Dice similarity coefficient. The median centroid distance between adapted and ground truth contours was 1.56 mm for the kidney, 3.81 mm for the gallbladder and 1.03 mm for the lung tumor (median Dice similarity coefficient: 0.95, 0.72 and 0.93). We observed that a decrease in image resolution led to a modest decrease in localization accuracy, especially for the small gallbladder. However, for all volumes of interest localization accuracy varied substantially more between subjects than due to the different imaging parameters. Automated tumor and organ localization using 2D cine MR imaging and template-matching-based autocontouring is robust against variation of geometric imaging parameters. Future work and optimization efforts of 2D cine MR imaging for real-time adaptive radiotherapy is needed to characterize the influence of sequence- and anatomical site-specific imaging contrast.

Texture analysis using proton density and T2 relaxation in patients with histological usual interstitial pneumonia (UIP) or nonspecific interstitial pneumonia (NSIP)

Objectives The purpose of our study was to assess proton density (PD) and T2 relaxation time of usual interstitial pneumonia (UIP) and nonspecific interstitial pneumonia (NSIP) and to evaluate their utility in differentiating the two patterns. Furthermore, we aim to investigate whether these two parameters could help differentiate active-inflammatory and stable-fibrotic lesions in NSIP. Methods 32 patients (mean age: 69 years; M:F, 1:1) with pathologically proven disease (UIP:NSIP, 1:1), underwent thoracic thin-section multislice CT scan and 1.5T MRI. A total of 437 regions-of-interest (ROIs) were classified at CT as advanced, moderate or mild alterations. Based on multi-echo single-shot TSE sequence acquired at five echo times, with breath-holding at end-expiration and ECG-triggering, entire lung T2 and PD maps were generated from each subject. The T2 relaxation time and the respective signal intensity were quantified by performing a ROI measurement on the T2 and PD maps in the corresponding CT selected areas of the lung. Results UIP and NSIP regional patterns could not be differentiated by T2 relaxation times or PD values alone. Overall, a strong positive correlation was found between T2 relaxation and PD in NSIP, r = 0.64, p<0.001; however, this correlation was weak in UIP, r = 0.20, p = 0.01. T2 relaxation showed significant statistical difference between active-inflammatory and stable-fibrotic NSIP regions at all levels, p<0.05, while for the analysis of ventral lesions PD proved no statistical difference, p>0.05. Conclusions T2 relaxation times and PD values may provide helpful quantitative information for differentiating NSIP from UIP pattern. These parameters have the potential to differentiate active-inflammatory and stable-fibrotic lesions in NSIP.

OC-0072: Respiratory time-resolved 4D MR imaging for RT applications with acquisition times below one minute

S33 ______________________________________________________________________________________________________ Conclusion: DIR-based registration methods showed that the vast majority of failures originated in the high dose target volumes and received full prescribed doses suggesting biological rather than technology-related causes of failure. Validated DIR-based registration is recommended for accurate failure characterization and a novel typologyindicative taxonomy is recommended for failure reporting in the IMRT era.