Whole-body FDG PET-MR oncologic imaging: pitfalls in clinical interpretation related to inaccurate MR-based attenuation correction

Simultaneous data collection for positron emission tomography and magnetic resonance imaging (PET/MR) is now a reality. While the full benefits of concurrently acquiring PET and MR data and the potential added clinical value are still being evaluated, initial studies have identified several important potential pitfalls in the interpretation of fluorodeoxyglucose (FDG) PET/MRI in oncologic whole-body imaging, the majority of which being related to the errors in the attenuation maps created from the MR data. The purpose of this article was to present such pitfalls and artifacts using case examples, describe their etiology, and discuss strategies to overcome them. Using a case-based approach, we will illustrate artifacts related to (1) Inaccurate bone tissue segmentation; (2) Inaccurate air cavities segmentation; (3) Motion-induced misregistration; (4) RF coils in the PET field of view; (5) B0 field inhomogeneity; (6) B1 field inhomogeneity; (7) Metallic implants; (8) MR contrast agents.

Image Guided Radiation Therapy Using Synthetic Computed Tomography Images in Brain Cancer

PURPOSE

The development of synthetic computed tomography (CT) (synCT) derived from magnetic resonance (MR) images supports MR-only treatment planning. We evaluated the accuracy of synCT and synCT-generated digitally reconstructed radiographs (DRRs) relative to CT and determined their performance for image guided radiation therapy (IGRT).

METHODS AND MATERIALS

Magnetic resonance simulation (MR-SIM) and CT simulation (CT-SIM) images were acquired of an anthropomorphic skull phantom and 12 patient brain cancer cases. SynCTs were generated using fluid attenuation inversion recovery, ultrashort echo time, and Dixon data sets through a voxel-based weighted summation of 5 tissue classifications. The DRRs were generated from the phantom synCT, and geometric fidelity was assessed relative to CT-generated DRRs through bounding box and landmark analysis. An offline retrospective analysis was conducted to register cone beam CTs (n=34) to synCTs and CTs using automated rigid registration in the treatment planning system. Planar MV and KV images (n=37) were rigidly registered to synCT and CT DRRs using an in-house script. Planar and volumetric registration reproducibility was assessed and margin differences were characterized by the van Herk formalism.

RESULTS

Bounding box and landmark analysis of phantom synCT DRRs were within 1 mm of CT DRRs. Absolute planar registration shift differences ranged from 0.0 to 0.7 mm for phantom DRRs on all treatment platforms and from 0.0 to 0.4 mm for volumetric registrations. For patient planar registrations, the mean shift differences were 0.4 ± 0.5 mm (range, -0.6 to 1.6 mm), 0.0 ± 0.5 mm (range, -0.9 to 1.2 mm), and 0.1 ± 0.3 mm (range, -0.7 to 0.6 mm) for the superior-inferior (S-I), left-right (L-R), and anterior-posterior (A-P) axes, respectively. The mean shift differences in volumetric registrations were 0.6 ± 0.4 mm (range, -0.2 to 1.6 mm), 0.2 ± 0.4 mm (range, -0.3 to 1.2 mm), and 0.2 ± 0.3 mm (range, -0.2 to 1.2 mm) for the S-I, L-R, and A-P axes, respectively. The CT-SIM and synCT derived margins were <0.3 mm different.

CONCLUSION

DRRs generated by synCT were in close agreement with CT-SIM. Planar and volumetric image registrations to synCT-derived targets were comparable with CT for phantom and patients. This validation is the next step toward MR-only planning for the brain.

Spatial Precision in Magnetic Resonance Imaging-Guided Radiation Therapy: The Role of Geometric Distortion

Because magnetic resonance imaging-guided radiation therapy (MRIgRT) offers exquisite soft tissue contrast and the ability to image tissues in arbitrary planes, the interest in this technology has increased dramatically in recent years. However, intrinsic geometric distortion stemming from both the system hardware and the magnetic properties of the patient affects MR images and compromises the spatial integrity of MRI-based radiation treatment planning, given that for real-time MRIgRT, precision within 2 mm is desired. In this article, we discuss the causes of geometric distortion, describe some well-known distortion correction algorithms, and review geometric distortion measurements from 12 studies, while taking into account relevant imaging parameters. Eleven of the studies reported phantom measurements quantifying system-dependent geometric distortion, while 2 studies reported simulation data quantifying magnetic susceptibility-induced geometric distortion. Of the 11 studies investigating system-dependent geometric distortion, 5 reported maximum measurements less than 2 mm. The simulation studies demonstrated that magnetic susceptibility-induced distortion is typically smaller than system-dependent distortion but still nonnegligible, with maximum distortion ranging from 2.1 to 2.6 mm at a field strength of 1.5 T. As expected, anatomic landmarks containing interfaces between air and soft tissue had the largest distortions. The evidence indicates that geometric distortion reduces the spatial integrity of MRI-based radiation treatment planning and likely diminishes the efficacy of MRIgRT. Better phantom measurement techniques and more effective distortion correction algorithms are needed to achieve the desired spatial precision.

Characterization of 3D geometric distortion of magnetic resonance imaging scanners commissioned for radiation therapy planning

OBJECTIVE

To develop a method for the assessment and characterization of 3D geometric distortion as part of routine quality assurance for MRI scanners commissioned for Radiation Therapy planning.

MATERIALS AND METHODS

In this study, the in-plane and through-plane geometric distortions on a 1.5T GE MRI-SIM unit are characterized and the 2D and 3D correction algorithms provided by the vendor are evaluated. We used a phantom developed by GE Healthcare that covers a large field of view of 500mm, and consists of layers of foam embedded with a matrix of ellipsoidal markers. An in-house Java-based software module was developed to automatically assess the geometric distortion by calculating the center of each marker using the center of mass method, correcting of gross rotation errors and comparing the corrected positions with a CT gold standard data set. Spatial accuracy of typical pulse sequences used in RT planning was assessed (2D T1/T2 FSE, 3D CUBE, T1 SPGR) using the software. The accuracy of vendor specific geometric distortion correction (GDC) algorithms was quantified by measuring distortions before and after the application of the 2D and 3D correction algorithms.

RESULTS

Our algorithm was able to accurately calculate geometric distortion with sub-pixel precision. For all typical MR sequences used in Radiotherapy, the vendor's GDC was able to substantially reduce the distortions. Our results showed also that the impact of the acquisition produced a maximum variation of 0.2mm over a radial distance of 200mm. It has been shown that while the 2D correction algorithm remarkably reduces the in-plane geometric distortion, 3D geometric distortion further reduced the geometric distortion by correcting both in-plane and through-plane distortions in all acquisitions.

CONCLUSION

The presented methods represent a valuable tool for routine quality assurance of MR applications that require stringent spatial accuracy assessment such as radiotherapy. The phantom used in this study provides three dimensional arrays of control points. These tools and the detailed results can be also used for developing new geometric distortion correction algorithms or improving the existing ones.

Continuous table acquisition MRI for radiotherapy treatment planning: distortion assessment with a new extended 3D volumetric phantom

PURPOSE

Accurate geometry is required for radiotherapy treatment planning (RTP). When considering the use of magnetic resonance imaging (MRI) for RTP, geometric distortions observed in the acquired images should be considered. While scanner technology and vendor supplied correction algorithms provide some correction, large distortions are still present in images, even when considering considerably smaller scan lengths than those typically acquired with CT in conventional RTP. This study investigates MRI acquisition with a moving table compared with static scans for potential geometric benefits for RTP.

METHODS

A full field of view (FOV) phantom (diameter 500 mm; length 513 mm) was developed for measuring geometric distortions in MR images over volumes pertinent to RTP. The phantom consisted of layers of refined plastic within which vitamin E capsules were inserted. The phantom was scanned on CT to provide the geometric gold standard and on MRI, with differences in capsule location determining the distortion. MRI images were acquired with two techniques. For the first method, standard static table acquisitions were considered. Both 2D and 3D acquisition techniques were investigated. With the second technique, images were acquired with a moving table. The same sequence was acquired with a static table and then with table speeds of 1.1 mm/s and 2 mm/s. All of the MR images acquired were registered to the CT dataset using a deformable B-spline registration with the resulting deformation fields providing the distortion information for each acquisition.

RESULTS

MR images acquired with the moving table enabled imaging of the whole phantom length while images acquired with a static table were only able to image 50%-70% of the phantom length of 513 mm. Maximum distortion values were reduced across a larger volume when imaging with a moving table. Increased table speed resulted in a larger contribution of distortion from gradient nonlinearities in the through-plane direction and an increased blurring of capsule images, resulting in an apparent capsule volume increase by up to 170% in extreme axial FOV regions. Blurring increased with table speed and in the central regions of the phantom, geometric distortion was less for static table acquisitions compared to a table speed of 2 mm/s over the same volume. Overall, the best geometric accuracy was achieved with a table speed of 1.1 mm/s.

CONCLUSIONS

The phantom designed enables full FOV imaging for distortion assessment for the purposes of RTP. MRI acquisition with a moving table extends the imaging volume in the z direction with reduced distortions which could be useful particularly if considering MR-only planning. If utilizing MR images to provide additional soft tissue information to the planning CT, standard acquisition sequences over a smaller volume would avoid introducing additional blurring or distortions from the through-plane table movement.

Initial clinical experience with a radiation oncology dedicated open 1.0T MR-simulation

The purpose of this study was to describe our experience with 1.0T MR-SIM including characterization, quality assurance (QA) program, and features necessary for treatment planning. Staffing, safety, and patient screening procedures were developed. Utilization of an external laser positioning system (ELPS) and MR-compatible couchtop were illustrated. Spatial and volumetric analyses were conducted between CT-SIM and MR-SIM using a stereotactic QA phantom with known landmarks and volumes. Magnetic field inhomogeneity was determined using phase difference analysis. System-related, in-plane distortion was evaluated and temporal changes were assessed. 3D distortion was characterized for regions of interest (ROIs) 5-20 cm away from isocenter. American College of Radiology (ACR) recommended tests and impact of ELPS on image quality were analyzed. Combined ultrashort echotime Dixon (UTE/Dixon) sequence was evaluated. Amplitude-triggered 4D MRI was implemented using a motion phantom (2-10 phases, ~ 2 cm excursion, 3-5 s periods) and a liver cancer patient. Duty cycle, acquisition time, and excursion were evaluated between maximum intensity projection (MIP) datasets. Less than 2% difference from expected was obtained between CT-SIM and MR-SIM volumes, with a mean distance of &lt; 0.2 mm between landmarks. Magnetic field inhomogeneity was &lt; 2 ppm. 2D distortion was &lt; 2 mm over 28.6-33.6 mm of isocenter. Within 5 cm radius of isocenter, mean 3D geometric distortion was 0.59 ± 0.32 mm (maximum = 1.65 mm) and increased 10-15 cm from isocenter (mean = 1.57 ± 1.06 mm, maximum = 6.26 mm). ELPS interference was within the operating frequency of the scanner and was characterized by line patterns and a reduction in signal-to-noise ratio (4.6-12.6% for TE = 50-150 ms). Image quality checks were within ACR recommendations. UTE/Dixon sequences yielded detectability between bone and air. For 4D MRI, faster breathing periods had higher duty cycles than slow (50.4% (3 s) and 39.4% (5 s), p &lt; 0.001) and ~fourfold acquisition time increase was measured for ten-phase versus two-phase. Superior-inferior object extent was underestimated 8% (6 mm) for two-phase as compared to ten-phase MIPs, although &lt; 2% difference was obtained for ≥ 4 phases. 4D MRI for a patient demonstrated acceptable image quality in ~ 7 min. MR-SIM was integrated into our workflow and QA procedures were developed. Clinical applicability was demonstrated for 4D MRI and UTE imaging to support MR-SIM for single modality treatment planning.

Evaluating organ delineation, dose calculation and daily localization in an open-MRI simulation workflow for prostate cancer patients

BACKGROUND

This study describes initial testing and evaluation of a vertical-field open Magnetic Resonance Imaging (MRI) scanner for the purpose of simulation in radiation therapy for prostate cancer. We have evaluated the clinical workflow of using open MRI as a sole modality for simulation and planning. Relevant results related to MRI alignment (vs. CT) reference dataset with Cone-Beam CT (CBCT) for daily localization are presented.

METHODS

Ten patients participated in an IRB approved study utilizing MRI along with CT simulation with the intent of evaluating the MRI-simulation process. Differences in prostate gland volume, seminal vesicles, and penile bulb were assessed with MRI and compared to CT. To evaluate dose calculation accuracy, bulk-density-assignments were mapped onto respective MRI datasets and treated IMRT plans were re-calculated. For image localization purposes, 400 CBCTs were re-evaluated with MRI as the reference dataset and daily shifts compared against CBCT-to-CT registration. Planning margins based on MRI/CBCT shifts were computed using the van Herk formalism.

RESULTS

Significant organ contour differences were noted between MRI and CT. Prostate volumes were on average 39.7% (p = 0.002) larger on CT than MRI. No significant difference was found in seminal vesicle volumes (p = 0.454). Penile bulb volumes were 61.1% higher on CT, without statistical significance (p = 0.074). MRI-based dose calculations with assigned bulk densities produced agreement within 1% with heterogeneity corrected CT calculations. The differences in shift positions for the cohort between CBCT-to-CT registration and CBCT-to-MRI registration are -0.15 ± 0.25 cm (anterior-posterior), 0.05 ± 0.19 cm (superior-inferior), and -0.01 ± 0.14 cm (left-right).

CONCLUSIONS

This study confirms the potential of using an open-field MRI scanner as primary imaging modality for prostate cancer treatment planning simulation, dose calculations and daily image localization.

The rotating biplanar linac-magnetic resonance imaging system

We have successfully built linac-magnetic resonance imaging (MR) systems based on a linac waveguide placed between open MR planes (perpendicular) or through the central opening of one of the planes (parallel) to improve dosimetric properties. It rotates on a gantry to irradiate at any angle. Irradiation during MR imaging and automatic 2-dimensional MR image-based target tracking and automatic beam steering to the moving target have been demonstrated with our systems. The functioning whole-body system (0.6-T MR and 6-MV linac) has been installed in an existing clinical vault without removing the walls or the ceiling and without the need of a helium exhaust vent.

MRI-based treatment planning with pseudo CT generated through atlas registration

PURPOSE

To evaluate the feasibility and accuracy of magnetic resonance imaging (MRI)-based treatment planning using pseudo CTs generated through atlas registration.

METHODS

A pseudo CT, providing electron density information for dose calculation, was generated by deforming atlas CT images previously acquired on other patients. The authors tested 4 schemes of synthesizing a pseudo CT from single or multiple deformed atlas images: use of a single arbitrarily selected atlas, arithmetic mean process using 6 atlases, and pattern recognition with Gaussian process (PRGP) using 6 or 12 atlases. The required deformation for atlas CT images was derived from a nonlinear registration of conjugated atlas MR images to that of the patient of interest. The contrasts of atlas MR images were adjusted by histogram matching to reduce the effect of different sets of acquisition parameters. For comparison, the authors also tested a simple scheme assigning the Hounsfield unit of water to the entire patient volume. All pseudo CT generating schemes were applied to 14 patients with common pediatric brain tumors. The image similarity of real patient-specific CT and pseudo CTs constructed by different schemes was compared. Differences in computation times were also calculated. The real CT in the treatment planning system was replaced with the pseudo CT, and the dose distribution was recalculated to determine the difference.

RESULTS

The atlas approach generally performed better than assigning a bulk CT number to the entire patient volume. Comparing atlas-based schemes, those using multiple atlases outperformed the single atlas scheme. For multiple atlas schemes, the pseudo CTs were similar to the real CTs (correlation coefficient, 0.787-0.819). The calculated dose distribution was in close agreement with the original dose. Nearly the entire patient volume (98.3%-98.7%) satisfied the criteria of chi-evaluation (<2% maximum dose and 2 mm range). The dose to 95% of the volume and the percentage of volume receiving at least 95% of the prescription dose in the planning target volume differed from the original values by less than 2% of the prescription dose (root-mean-square, RMS < 1%). The PRGP scheme did not perform better than the arithmetic mean process with the same number of atlases. Increasing the number of atlases from 6 to 12 often resulted in improvements, but statistical significance was not always found.

CONCLUSIONS

MRI-based treatment planning with pseudo CTs generated through atlas registration is feasible for pediatric brain tumor patients. The doses calculated from pseudo CTs agreed well with those from real CTs, showing dosimetric accuracy within 2% for the PTV when multiple atlases were used. The arithmetic mean process may be a reasonable choice over PRGP for the synthesis scheme considering performance and computational costs.

A simple and efficient methodology to improve geometric accuracy in gamma knife radiation surgery: implementation in multiple brain metastases

PURPOSE

To propose, verify, and implement a simple and efficient methodology for the improvement of total geometric accuracy in multiple brain metastases gamma knife (GK) radiation surgery.

METHODS AND MATERIALS

The proposed methodology exploits the directional dependence of magnetic resonance imaging (MRI)-related spatial distortions stemming from background field inhomogeneities, also known as sequence-dependent distortions, with respect to the read-gradient polarity during MRI acquisition. First, an extra MRI pulse sequence is acquired with the same imaging parameters as those used for routine patient imaging, aside from a reversal in the read-gradient polarity. Then, "average" image data are compounded from data acquired from the 2 MRI sequences and are used for treatment planning purposes. The method was applied and verified in a polymer gel phantom irradiated with multiple shots in an extended region of the GK stereotactic space. Its clinical impact in dose delivery accuracy was assessed in 15 patients with a total of 96 relatively small (<2 cm) metastases treated with GK radiation surgery.

RESULTS

Phantom study results showed that use of average MR images eliminates the effect of sequence-dependent distortions, leading to a total spatial uncertainty of less than 0.3 mm, attributed mainly to gradient nonlinearities. In brain metastases patients, non-eliminated sequence-dependent distortions lead to target localization uncertainties of up to 1.3 mm (mean: 0.51 ± 0.37 mm) with respect to the corresponding target locations in the "average" MRI series. Due to these uncertainties, a considerable underdosage (5%-32% of the prescription dose) was found in 33% of the studied targets.

CONCLUSIONS

The proposed methodology is simple and straightforward in its implementation. Regarding multiple brain metastases applications, the suggested approach may substantially improve total GK dose delivery accuracy in smaller, outlying targets.

Integrated image reconstruction and gradient nonlinearity correction

PURPOSE

To describe a model-based reconstruction strategy for routine magnetic resonance imaging that accounts for gradient nonlinearity (GNL) during rather than after transformation to the image domain, and demonstrate that this approach reduces the spatial resolution loss that occurs during strictly image-domain GNL-correction.

METHODS

After reviewing conventional GNL-correction methods, we propose a generic signal model for GNL-affected magnetic resonance imaging acquisitions, discuss how it incorporates into contemporary image reconstruction platforms, and describe efficient nonuniform fast Fourier transform-based computational routines for these. The impact of GNL-correction on spatial resolution by the conventional and proposed approaches is investigated on phantom data acquired at varying offsets from gradient isocenter, as well as on fully sampled and (retrospectively) undersampled in vivo acquisitions.

RESULTS

Phantom results demonstrate that resolution loss that occurs during GNL-correction is significantly less for the proposed strategy than for the standard approach at distances >10 cm from isocenter with a 35 cm field-of-view gradient coil. The in vivo results suggest that the proposed strategy better preserves fine anatomical detail than retrospective GNL-correction while offering comparable geometric correction.

CONCLUSION

Accounting for GNL during image reconstruction allows geometric distortion to be corrected with less spatial resolution loss than is typically observed with the conventional image domain correction strategy.

A dual model HU conversion from MRI intensity values within and outside of bone segment for MRI-based radiotherapy treatment planning of prostate cancer

PURPOSE

The lack of electron density information in magnetic resonance images (MRI) poses a major challenge for MRI-based radiotherapy treatment planning (RTP). In this study the authors convert MRI intensity values into Hounsfield units (HUs) in the male pelvis and thus enable accurate MRI-based RTP for prostate cancer patients with varying tissue anatomy and body fat contents.

METHODS

T1/T2*-weighted MRI intensity values and standard computed tomography (CT) image HUs in the male pelvis were analyzed using image data of 10 prostate cancer patients. The collected data were utilized to generate a dual model HU conversion technique from MRI intensity values of the single image set separately within and outside of contoured pelvic bones. Within the bone segment local MRI intensity values were converted to HUs by applying a second-order polynomial model. This model was tuned for each patient by two patient-specific adjustments: MR signal normalization to correct shifts in absolute intensity level and application of a cutoff value to accurately represent low density bony tissue HUs. For soft tissues, such as fat and muscle, located outside of the bone contours, a threshold-based segmentation method without requirements for any patient-specific adjustments was introduced to convert MRI intensity values into HUs. The dual model HU conversion technique was implemented by constructing pseudo-CT images for 10 other prostate cancer patients. The feasibility of these images for RTP was evaluated by comparing HUs in the generated pseudo-CT images with those in standard CT images, and by determining deviations in MRI-based dose distributions compared to those in CT images with 7-field intensity modulated radiation therapy (IMRT) with the anisotropic analytical algorithm and 360° volumetric-modulated arc therapy (VMAT) with the Voxel Monte Carlo algorithm.

RESULTS

The average HU differences between the constructed pseudo-CT images and standard CT images of each test patient ranged from -2 to 5 HUs and from 22 to 78 HUs in soft and bony tissues, respectively. The average local absolute value differences were 11 HUs in soft tissues and 99 HUs in bones. The planning target volume doses (volumes 95%, 50%, 5%) in the pseudo-CT images were within 0.8% compared to those in CT images in all of the 20 treatment plans. The average deviation was 0.3%. With all the test patients over 94% (IMRT) and 92% (VMAT) of dose points within body (lower than 10% of maximum dose suppressed) passed the 1 mm and 1% 2D gamma index criterion. The statistical tests (t- and F-tests) showed significantly improved (p ≤ 0.05) HU and dose calculation accuracies with the soft tissue conversion method instead of homogeneous representation of these tissues in MRI-based RTP images.

CONCLUSIONS

This study indicates that it is possible to construct high quality pseudo-CT images by converting the intensity values of a single MRI series into HUs in the male pelvis, and to use these images for accurate MRI-based prostate RTP dose calculations.

Estimation of geometrically undistorted B(0) inhomogeneity maps

Geometric accuracy of MRI is one of the main concerns for its use as a sole image modality in precision radiation therapy (RT) planning. In a state-of-the-art scanner, system level geometric distortions are within acceptable levels for precision RT. However, subject-induced B0 inhomogeneity may vary substantially, especially in air-tissue interfaces. Recent studies have shown distortion levels of more than 2 mm near the sinus and ear canal are possible due to subject-induced field inhomogeneity. These distortions can be corrected with the use of accurate B0 inhomogeneity field maps. Most existing methods estimate these field maps from dual gradient-echo (GRE) images acquired at two different echo-times under the assumption that the GRE images are practically undistorted. However distortion that may exist in the GRE images can result in estimated field maps that are distorted in both geometry and intensity, leading to inaccurate correction of clinical images. This work proposes a method for estimating undistorted field maps from GRE acquisitions using an iterative joint estimation technique. The proposed method yields geometrically corrected GRE images and undistorted field maps that can also be used for the correction of images acquired by other sequences. The proposed method is validated through simulation, phantom experiments and applied to patient data. Our simulation results show that our method reduces the root-mean-squared error of the estimated field map from the ground truth by ten-fold compared to the distorted field map. Both the geometric distortion and the intensity corruption (artifact) in the images caused by the B0 field inhomogeneity are corrected almost completely. Our phantom experiment showed improvement in the geometric correction of approximately 1 mm at an air-water interface using the undistorted field map compared to using a distorted field map. The proposed method for undistorted field map estimation can lead to improved geometric distortion correction at air-tissue interfaces, especially in low readout-bandwidth acquisitions, thus making them suitable for clinical use in precision RT without increasing the treatment planning margin.

Toward magnetic resonance-only simulation: segmentation of bone in MR for radiation therapy verification of the head

PURPOSE

To develop a practical method to localize bones in magnetic resonance (MR) images, to create "computed tomography-like" MR images (ctMRI) that could be used for radiation therapy verification, and to generate MR-based digitally reconstructed radiographs (DRR).

METHODS AND MATERIALS

Using T1-weighted MR images, an air mask was derived from the manual contouring of all airways within the head and neck region using axial images at 6 anatomic levels. Compact bone, spongy bone, and soft tissue masks were then automatically generated using the statistical data derived from MR intensities and the air mask. ctMRI were then generated by mapping the MR intensities of the voxels within these masks into the CT number ranges of corresponding tissues. MR-based DRRs created from ctMRI were quantitatively evaluated using the co-registered MR and CT head images of 20 stereotactic radiosurgery patients. Ten anatomical points, positioned on the skull segmented using a threshold of 300 HU, in CT and ctMRI, were used to determine the differences in distance between MR-based DRRs and CT-based DRRs, and to evaluate the geometric accuracy of ctMRI and MR-based DRRs.

RESULTS

The bony structures were identified on ctMRI and were visible in the MR-based DRRs. From the 20 patient cases, the mean geometric difference and standard deviation between the 10 anatomical points on MR-based and CT-based DRRs was -0.05 ± 0.85 mm, respectively. This included uncertainty in image fusion. The maximum distance difference was 1.88 mm.

CONCLUSIONS

A practical method was developed to segment bone from MR images. The ctMRI created can be used for radiation treatment verification when MR-only simulation is performed. MR-based DRRs can be used in place of CT-based DRRs.

Assessment and characterization of the total geometric uncertainty in Gamma Knife radiosurgery using polymer gels

PURPOSE

This work proposes and implements an experimental methodology, based on polymer gels, for assessing the total geometric uncertainty and characterizing its contributors in Gamma Knife (GK) radiosurgery.

METHODS

A treatment plan consisting of 26, 4-mm GK single shot dose distributions, covering an extended region of the Leksell stereotactic space, was prepared and delivered to a polymer gel filled polymethyl methacrylate (PMMA) head phantom (16 cm diameter) used to accurately reproduce every link in the GK treatment chain. The center of each shot served as a "control point" in the assessment of the GK total geometric uncertainty, which depends on (a) the spatial dose delivery uncertainty of the PERFEXION GK unit used in this work, (b) the spatial distortions inherent in MR images commonly used for target delineation, and (c) the geometric uncertainty contributor associated with the image registration procedure performed by the Leksell GammaPlan (LGP) treatment planning system (TPS), in the case that registration is directly based on the apparent fiducial locations depicted in each MR image by the N-shaped rods on the Leksell localization box. The irradiated phantom was MR imaged at 1.5 T employing a T2-weighted pulse sequence. Four image series were acquired by alternating the frequency encoding axis and reversing the read gradient polarity, thus allowing the characterization of the MR-related spatial distortions.

RESULTS

MR spatial distortions stemming from main field (B0) inhomogeneity as well as from susceptibility and chemical shift phenomena (also known as sequence dependent distortions) were found to be of the order of 0.5 mm, while those owing to gradient nonlinearities (also known as sequence independent distortions) were found to increase with distance from the MR scanner isocenter extending up to 0.47 mm at an Euclidean distance of 69.6 mm. Regarding the LGP image registration procedure, the corresponding average contribution to the total geometric uncertainty ranged from 0.34 to 0.80 mm. The average total geometric uncertainty, which also includes the GK spatial dose delivery uncertainty, was found equal to (0.88 ± 0.16), (0.88 ± 0.26), (1.02 ± 0.09), and (1.15 ± 0.24) mm for the MR image series acquired with the read gradient polarity (direction) set toward right, left, posterior, and anterior, respectively.

CONCLUSIONS

The implemented methodology seems capable of assessing the total geometric uncertainty, as well as of characterizing its contributors, ascribed to the entire GK treatment delivery (i.e., from MR imaging to GK dose delivery) for an extended region of the Leksell stereotactic space. Results obtained indicate that the selection of both the frequency encoding axis and the read gradient polarity during MRI acquisition may affect the magnitude as well as the spatial components of the total geometric uncertainty.

Patient-induced susceptibility effect on geometric distortion of clinical brain MRI for radiation treatment planning on a 3T scanner

Concerns about the geometric accuracy of MRI in radiation therapy (RT) have been present since its invention. Although modern scanners typically have system levels of geometric accuracy that meet requirements of RT, subject-specific distortion is variable, and methods to in vivo assess and control patient-induced geometric distortion are not yet resolved. This study investigated the nature and magnitude of the subject-induced susceptibility effect on geometric distortions in clinical brain MRI, and tested the feasibility of in vivo quality control using field inhomogeneity mapping. For 19 consecutive patients scanned on a dedicated 3T MR scanner, B0 field inhomogeneity maps were acquired and analyzed to determine subject-induced distortions. For 3D T1 weighted images frequency-encoded with a bandwidth of 180 Hz/pixel, 86.9% of the estimated displacements were <0.5 mm, 97.4% <1 mm, and only 0.1% of displacements > 2 mm. The maximum displacement was <4 mm. The greatest distortions were observed at the interfaces with air at the sinuses. Displacements decayed to less than 1 mm over a distance of 8 mm. Metal surgical wires generated smaller distortions, with an averaged maximum displacement of 0.76 mm. Repeat acquisition of the field maps in 17 patients revealed a within-subject standard deviation of 0.25 ppm, equivalent to 0.22 mm displacement in the frequency-encoding direction in the 3D T1 weighted images. Susceptibility-induced voxel displacements in the brain are generally small, but should be monitored for precision RT. These effects are manageable at 3T and lower fields, and the methods applied can be used to monitor for potential local errors in individual patients, as well as to correct for local distortions as needed.

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.

Implications of tissue magnetic susceptibility-related distortion on the rotating magnet in an MR-linac design

PURPOSE

One of the recently published concepts that combine the soft-tissue imaging capabilities of MRI with external beam radiotherapy involves the rigid coupling of a linac with a rotating biplanar low-field MR imaging system. While such a system would prevent possible image distortion resulting from relative motion between the magnet and the linac, the rotation of the magnet around the patient can itself introduce possibilities for image distortion that need to be addressed. While there are straightforward techniques in the literature for correcting distortions from gradient nonlinearities and nonuniform magnetic fields during image reconstruction, the correction of distortions related to tissue magnetic susceptibility is more complex. This work investigates the extent of this latter distortion type under the regime of a rotating magnetic field.

METHODS

CT images covering patient anatomy in the head, lung, and male pelvic regions were obtained and segmented into components of air, bone, and soft tissue. Each of these three components was assigned bulk magnetic susceptibility values in accordance with those found in the literature. A finite-difference algorithm was then implemented to solve for magnetic field distortion maps should the anatomies be placed in the uniform polarizing field of an MR system. The algorithm was repeated multiple times as the polarizing field was rotated axially about the virtual patient in 15 degrees increments. In this way, a map of maximum distortion, and the range of distortion as the magnetic field is rotated about each anatomical region could be determined. The consequence of these susceptibility distortions in terms of geometric signal shift was calculated for 0.2 T, as well as another low-field system (0.5 T), and a higher field 1.5 T system for comparison, using the assumption of a frequency encoding gradient strength of 5 mT/m.

RESULTS

At 0.2 T, the susceptibility-related distortion was limited to less than 0.5 mm given an encoding gradient strength of 5 mT/m or higher. To maintain this same level of geometric accuracy, the 0.5 T system would require a moderately higher minimum gradient strength of 11 mT/m, and at a typical MR field strength of 1.5 T this minimum gradient strength would increase to 33 mT/m. The influence of magnetic susceptibility on mean frequency shift as the field orientation was rotated was also investigated and found to account for less than half a millimeter at 1.5 T, and negligible for low-field systems.

CONCLUSIONS

A study of three sites (head, lung, and prostate) that are vulnerable to magnetic susceptibility-related distortions were studied, and showed that in the context of a rotating polarizing magnet, low-field systems can maintain geometric accuracy of 0.5 mm with at most moderate limitations on sequence parameters. This conclusion will likely apply only to endogenous tissues, as implanted materials such as titanium can create field distortions much in excess of what may normally be induced in the body. Items containing such materials (hip prostheses, for example) will require individual scrutiny.