Christopher M. Rank

MRI-based treatment plan simulation and adaptation for ion radiotherapy using a classification-based approach

BackgroundIn order to benefit from the highly conformal irradiation of tumors in ion radiotherapy, sophisticated treatment planning and simulation are required. The purpose of this study was to investigate the potential of MRI for ion radiotherapy treatment plan simulation and adaptation using a classification-based approach.MethodsFirstly, a voxelwise tissue classification was applied to derive pseudo CT numbers from MR images using up to 8 contrasts. Appropriate MR sequences and parameters were evaluated in cross-validation studies of three phantoms. Secondly, ion radiotherapy treatment plans were optimized using both MRI-based pseudo CT and reference CT and recalculated on reference CT. Finally, a target shift was simulated and a treatment plan adapted to the shift was optimized on a pseudo CT and compared to reference CT optimizations without plan adaptation.ResultsThe derivation of pseudo CT values led to mean absolute errors in the range of 81 - 95 HU. Most significant deviations appeared at borders between air and different tissue classes and originated from partial volume effects. Simulations of ion radiotherapy treatment plans using pseudo CT for optimization revealed only small underdosages in distal regions of a target volume with deviations of the mean dose of PTV between 1.4 - 3.1% compared to reference CT optimizations. A plan adapted to the target volume shift and optimized on the pseudo CT exhibited a comparable target dose coverage as a non-adapted plan optimized on a reference CT.ConclusionsWe were able to show that a MRI-based derivation of pseudo CT values using a purely statistical classification approach is feasible although no physical relationship exists. Large errors appeared at compact bone classes and came from an imperfect distinction of bones and other tissue types in MRI. In simulations of treatment plans, it was demonstrated that these deviations are comparable to uncertainties of a target volume shift of 2 mm in two directions indicating that especially applications for adaptive ion radiotherapy are interesting.

MRI-based simulation of treatment plans for ion radiotherapy in the brain region

PURPOSE To test the potential of MRI-based treatment plan simulation for ion radiotherapy in the brain region. MATERIALS AND METHODS A classification-based tissue segmentation method based on discriminant analysis was employed to derive so-called pseudo CT numbers from MR images of three patients with lesions in the head region undergoing ion radiotherapy. Treatment plans for ions, and for comparison purposes also for photons, were subsequently optimized and simulated using both MRI-based pseudo CT and a standard X-ray-based reference CT. RESULTS Pseudo CTs revealed mean absolute errors in CT number in the range of 141-165 HU. While soft tissue was in good agreement with reference CT values, large deviations appeared at air cavities and bones as well as at interfaces of different tissue types. In simulations of ion treatment plans, pseudo CT optimizations showed small underdosages of target volumes with deviations in the PTV mean dose of 0.4-2.0% in comparison to reference CT optimizations. In contrast, the PTV mean dose in photon treatment plans differed by no more than 0.2%. CONCLUSIONS The main challenge in deriving pseudo CT numbers from MRI was the correct assignment of air and compact bone. In this study, the impact of deviations on simulations of ion and photon treatment plans in the brain region was small, however for more complicated morphologies a further improvement of the classification method including MR imaging of compact bone is required.

4D respiratory motion-compensated image reconstruction of free-breathing radial MR data with very high undersampling.

To develop four‐dimensional (4D) respiratory time‐resolved MRI based on free‐breathing acquisition of radial MR data with very high undersampling.

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.

MR–Consistent Simultaneous Reconstruction of Attenuation and Activity for Non–TOF PET/MR

Attenuation correction (AC) is required for accurate quantification of the reconstructed activity distribution in positron emission tomography (PET). For simultaneous PET/magnetic resonance (MR), however, AC is challenging, since the MR images do not provide direct information on the attenuating properties of the underlying tissue. Standard MR-based AC does not account for the presence of bone and thus leads to an underestimation of the activity distribution. To improve quantification for non-time-of-flight PET/MR, we propose an algorithm which simultaneously reconstructs activity and attenuation distribution from the PET emission data using available MR images as anatomical prior information. The MR information is used to derive voxel-dependent expectations on the attenuation coefficients. The expectations are modeled using Gaussian-like probability functions. An iterative reconstruction scheme incorporating the prior information on the attenuation coefficients is used to update attenuation and activity distribution in an alternating manner. We tested and evaluated the proposed algorithm for simulated 3D PET data of the head and the pelvis region. Activity deviations were below 5% in soft tissue and lesions compared to the ground truth whereas standard MR-based AC resulted in activity underestimation values of up to 12%.

MLAA-based attenuation correction of flexible hardware components in hybrid PET/MR imaging

BackgroundAccurate PET quantification demands attenuation correction (AC) for both patient and hardware attenuation of the 511 keV annihilation photons. In hybrid PET/MR imaging, AC for stationary hardware components such as patient table and MR head coil is straightforward, employing CT-derived attenuation templates. AC for flexible hardware components such as MR-safe headphones and MR radiofrequency (RF) surface coils is more challenging. Registration-based approaches, aligning CT-based attenuation templates with the current patient position, have been proposed but are not used in clinical routine. Ignoring headphone or RF coil attenuation has been shown to result in regional activity underestimation values of up to 18%.We propose to employ the maximum-likelihood reconstruction of attenuation and activity (MLAA) algorithm to estimate the attenuation of flexible hardware components. Starting with an initial attenuation map not including flexible hardware components, the attenuation update of MLAA is applied outside the body outline only, allowing to estimate hardware attenuation without modifying the patient attenuation map. Appropriate prior expectations on the attenuation coefficients are incorporated into MLAA. The proposed method is investigated for non-TOF PET phantom and 18F-FDG patient data acquired with a clinical PET/MR device, using headphones or RF surface coils as flexible hardware components.ResultsAlthough MLAA cannot recover the exact physical shape of the hardware attenuation maps, the overall attenuation of the hardware components is accurately estimated. Therefore, the proposed algorithm significantly improves PET quantification. Using the phantom data, local activity underestimation when neglecting hardware attenuation was reduced from up to 25% to less than 3% under- or overestimation as compared to reference scans without hardware present or to CT-derived AC. For the patient data, we found an average activity underestimation of 7.9% evaluated in the full brain and of 6.1% for the abdominal region comparing the uncorrected case with MLAA.ConclusionsMLAA is able to provide accurate estimations of the attenuation of flexible hardware components and can therefore be used to significantly improve PET quantification. The proposed approach can be readily incorporated into clinical workflow.

Investigation of the halo-artifact in 68Ga-PSMA-11-PET/MRI

Objectives Combined positron emission tomography (PET) and magnetic resonance imaging (MRI) targeting the prostate-specific membrane antigen (PSMA) with a 68Ga-labelled PSMA-analog (68Ga-PSMA-11) is discussed as a promising diagnostic method for patients with suspicion or history of prostate cancer. One potential drawback of this method are severe photopenic (halo-) artifacts surrounding the bladder and the kidneys in the scatter-corrected PET images, which have been reported to occur frequently in clinical practice. The goal of this work was to investigate the occurrence and impact of these artifacts and, secondly, to evaluate variants of the standard scatter correction method with regard to halo-artifact suppression. Methods Experiments using a dedicated pelvis phantom were conducted to investigate whether the halo-artifact is modality-, tracer-, and/or concentration-dependent. Furthermore, 31 patients with history of prostate cancer were selected from an ongoing 68Ga-PSMA-11-PET/MRI study. For each patient, PET raw data were reconstructed employing six different variants of PET scatter correction: absolute scatter scaling, relative scatter scaling, and relative scatter scaling combined with prompt gamma correction, each of which was combined with a maximum scatter fraction (MaxSF) of MaxSF = 75% or MaxSF = 40%. Evaluation of the reconstructed images with regard to halo-artifact suppression was performed both quantitatively using statistical analysis and qualitatively by two independent readers. Results The phantom experiments did not reveal any modality-dependency (PET/MRI vs. PET/CT) or tracer-dependency (68Ga vs. 18F-FDG). Patient- and phantom-based data indicated that halo-artifacts derive from high organ-to-background activity ratios (OBR) between bladder/kidneys and surrounding soft tissue, with a positive correlation between OBR and halo size. Comparing different variants of scatter correction, reducing the maximum scatter fraction from the default value MaxSF = 75% to MaxSF = 40% was found to efficiently suppress halo-artifacts in both phantom and patient data. In 1 of 31 patients, reducing the maximum scatter fraction provided new PET-based information changing the patient’s diagnosis. Conclusion Halo-artifacts are particularly observed for 68Ga-PSMA-11-PET/MRI due to 1) the biodistribution of the PSMA-11-tracer resulting in large OBRs for bladder and kidneys and 2) inaccurate scatter correction methods currently used in clinical routine, which tend to overestimate the scatter contribution. If not compensated for, 68Ga-PSMA-11 uptake pathologies may be masked by halo-artifacts leading to false-negative diagnoses. Reducing the maximum scatter fraction was found to efficiently suppress halo-artifacts.

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.

Improved clinical workflow for simultaneous whole-body PET/MRI using high-resolution CAIPIRINHA-accelerated MR-based attenuation correction

PURPOSE To explore the value and reproducibility of a novel magnetic resonance based attenuation correction (MRAC) using a CAIPIRINHA-accelerated T1-weighted Dixon 3D-VIBE sequence for whole-body PET/MRI compared to the clinical standard. METHODS The PET raw data of 19 patients from clinical routine were reconstructed with standard MRAC (MRACstd) and the novel MRAC (MRACcaipi), a prototype CAIPIRINHA accelerated Dixon 3D-VIBE sequence, both acquired in 19 s/bed position. Volume of interests (VOIs) for liver, lung and all voxels of the total image stack were created to calculate standardized uptake values (SUVmean) followed by inter-method agreement (Passing-Bablok regression, Bland-Altman analysis). A voxel-wise SUV comparison per patient was performed for intra-individual correlation between MRACstd and MRACcaipi. Difference images (MRACstd-MRACcaipi) of attenuation maps and SUV images were calculated. The image quality of in/opposed-phase water and fat images obtained from MRACcaipi was assessed by two readers on a 5-point Likert-scale including intra-class coefficients for inter-reader agreement. RESULTS SUVmean correlations of VOIs demonstrated high linearity (0.95<Spearmans rho<1, p<0.0001, respectively), substantiated by voxel-wise SUV scatter-plots (1.79×108 pixels). Outliers could be explained by different physiological conditions between the scans such as different segmentation of air-containing tissue, lungs, kidneys, metal implants, diaphragm edge or small air bubbles in the gastrointestinal tracts that moved between MRAC acquisitions. Nasal sinuses and the trachea were better segmented in MRACcaipi. High-resolution T1w Dixon 3D VIBE images were acquired in all cases and could be used for PET/MRI fusion. MRACcaipi images were of high diagnostic quality (4.2±0.8) with 0.92-0.96 intra-class correlation. CONCLUSIONS The novel prototype MRACcaipi extends the value for attenuation correction by providing a high spatial resolution DIXON-based dataset suited for diagnostic assessment towards time-efficient whole-body PET/MRI.