- Open Access
Evaluating organ delineation, dose calculation and daily localization in an open-MRI simulation workflow for prostate cancer patients
Radiation Oncologyvolume 10, Article number: 37 (2015)
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.
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.
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).
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.
Magnetic resonance imaging (MRI) is a supplement to computed tomography (CT) in radiation oncology where superior soft tissue contrast is required for anatomical delineation. MRI simulation in radiation therapy (RT) could potentially be a routine modality in the future primarily because MRI offers the ability to image tumors and surrounding healthy tissues with significantly better soft tissue contrast than CT (lung being a notable exception). MRI is an established standard for target and organ-at-risk delineation for brain and spinal cord cancers and is used in conjunction with CT for treatment planning [1-3]. In prostate treatment in particular, MRI as a visualization tool for structures has been investigated [4-7].
MRI-only workflows have been introduced with the primary endpoint of eliminating the need for CT simulation (CT-SIM). Devic  provides an excellent review of the current use of MRI in radiation oncology, including MRI simulation. Recently, MRI-only simulators have been introduced, incorporating features commonly found in CT-SIM: flat table tops to accommodate immobilization devices and external laser systems for patient marking and alignment. These two major additions facilitate better agreement in patient positioning between MRI and radiation therapy delivery. As MRI technology continues to evolve, new radiation therapy systems are utilizing MRI as an on-board localization technology, such as the three cobalt source unit from ViewRay (Oakwood Village, OH) , the MRI-LINAC from Elekta (Stockholm, SE) , or the MRI-on-rails solution being installed at Princess Margaret Hospital in Toronto, CA .
Nevertheless, CT-SIM remains the gold standard for treatment planning because CT voxels demonstrate a tissue’s electron density, a factor critical to determining the absorbed dose, as well excellent geometric accuracy and less stringent safety considerations. Electron density information is not easily accessible from MRI; and in an MRI-only simulation workflow treatment planners need to assign bulk densities to multiple contours that are not currently generated in a CT workflow. MRI treatment planning is currently a topic of considerable interest in radiation oncology [12-15]. Despite its promise, issues related to geometric distortions (stemming from magnetic field inhomogeneities [16,17] and patient-induced artifacts [18,19]), dose calculation [20-22] and the use of MRI datasets for image-guidance need to be optimized to facilitate accurate and efficient clinical workflow.
This study’s goal was to evaluate the potential of vertical-field open design MRI scanner as a primary modality for simulation, planning, and IGRT-based localization for treatment of patients with prostate cancers. Dose calculation accuracy was evaluated, including the use of bulk density assignments to overcome limitations in electron density conversion of the MRI signal. A unique aspect of our study centers on using MRI as a reference dataset for CBCT-MRI, 3D-3D image registration for image-guided RT. Such data is important for routine clinical use of MRI, and has not been previously published. Additionally, for an open platform MRI simulator, we investigated issues related to image acquisition sequences, and the possible pitfalls of improper imaging and its effect on anatomical contouring.
Ten prostate cancer patients enrolled on an Institutional Review Board approved study (to be in compliance of the Helsinki declaration) had a MRI acquired along with standard CT simulation. MRIs were obtained from a one Tesla (1.0 T) open MRI scanner (Panorama, Philips Healthcare, Best, Netherlands) that has been developed as a simulation device for RT procedures. Figure 1 shows our workflow for CT and MRI simulation. Whole pelvis acquisitions were performed using T1-weighted Fast Field Echo, T2-weighted Turbo Spin Echo and 3D balanced Turbo Field Echo (bTFE) sequences. Table 1 shows acquisition parameters for these scans. T1 and T2 scans were full field of view (FOV) in order to use them for treatment planning. The bTFE scans were ‘coned-down’ to smaller regions of interest to help keep acquisition time down and allow for improved organ delineation. Of note is that these acquisition parameters are optimized for RT planning by the manufacturer, and are not applicable to diagnostic imaging protocols. A flat table top insert (Civco, Orange City, IA) was used to reproduce typical patient positioning.
The MRI simulation process shares many aspects with CT simulation. The anatomy of interest (pelvis, head, etc.) is externally targeted and then moved to the magnet isocenter. A three plane survey (analogous to CT scouts) is performed with the largest FOV to determine the acquisition volume. Acquisition time for three whole pelvis sequences totals twenty minutes.
All anatomical delineation and planning were done in Eclipse® (Varian Medical Systems, Palo Alto, CA). First, all MRIs were rigidly registered to CT images focusing on a match of the prostate. Secondly, organ delineation was performed/reviewed by a single observer for the T1-MRI and CT images. Retrospectively, another observer contoured the T2 and bTFE volumes. All contours were reviewed by an attending physician.
All patients were planned using intensity modulated radiation therapy (IMRT) with seven or nine fields. The clinical target volume (CTV) was defined as prostate or as prostate plus proximal seminal vesicles based on disease staging. For patient treatment, the prostate volume was based on the CT volume instead of the MRI drawn prostate volume. The planned target volume (PTV) was an expansion of the CTV by 1 cm in all directions except 0.6 cm posteriorly. IMRT dose constraints followed QUANTEC  guidelines. Plans were prescribed so the dose to ninety-five percent of the volume (D95) was 75.6 Gy in 42 fractions.
Retrospective dose calculations were performed on T1-MRI images using bulk-density-assignment for each of the structures (body contour, PTV, seminal vesicles, bladder, rectum, penile bulb and bony structure), by assigning bony tissue to 480 Hounsfield units (HU) and soft tissue to 0 HU (water equivalent). The bone value (480 HU) was determined by contouring bony tissue on CT and measuring the average CT number. Soft tissue is comprised of fat and muscle which have HU values slightly lower and higher than water, respectively, so we chose to represent all soft tissue with water density. No air pockets were present in the treated area. Had they been present, a value of −1000 HU would have been used. Bulk-density-assignment for dose calculation using MRI datasets in the context of prostate cancer have been used by other investigators [13,21]. MRI datasets cannot be used to generate dose plans in Eclipse v11. Instead, unit-density CT datasets were registered to MRI datasets and MRI contours copied to CT to create workable datasets for dose calculation. The same beam angles, leaf segment positions and monitor units were calculated to isocenter. We evaluated the impact of bulk-density-assignment on IMRT optimization using a few representative datasets. Inverse planning was performed using bulk-density-assignment and plans were compared to those in which actual CT densities were used. For both datasets, the identical beam arrangements, plan and optimization parameters were used for inverse planning. The calculated DVHs for the prostate gland, bladder and rectum were compared.
CBCT-based treatment localization
IMRT delivery was performed on Trilogy® and Truebeam® linear accelerators (Varian Medical Systems, Palo Alto, CA). All patients underwent daily CBCT localization and shifts were recorded. Retrospectively, another observer used MRI datasets as the primary reference sets, allowing for virtual realignment of all CBCT images (~400 datasets). MRI-CBCT registration was performed manually prioritizing alignment of the prostate-rectum interface correctly, which is standard practice. By analyzing these shifts and using the van Herk formalism , we calculated the CTV-to-PTV margins needed for the CT and MRI reference sets. This formalism calculates the margins needed to deliver 95% of the prescribed dose to 90% of the patient population. It is calculated using the cohort’s inter-patient uncertainty (Σ) and the inter-fraction uncertainty (σ) in the formula .
Margin = 2.5 Σ + 0.7 σ. We also calculated the group’s mean error (M). The inter-patient uncertainty (Σ) is calculated by averaging all daily shifts and then taking the standard deviation of each patient’s average value. The inter-fraction uncertainty (σ) is calculated by taking the standard deviation of all daily shifts and then taking the root mean square (RMS) of each patient’s standard deviation.
Volumetric differences in contours from MRI and CT datasets are summarized in Table 2. Significant differences were noted between MRI and CT. Prostate and penile bulb volumes were on average 39.7% (p = 0.002) and 61.1% (p = 0.074) larger on CT than on MRI respectively. The prostate volume difference agrees well with values from Roach et al.  Seminal vesicle volumes were underestimated on the CT scan by 1.7% (p = 0.454) compared to the volumes from MRI. CTV (prostate plus seminal vesicles) were on average 27.5% (p = 0.002) smaller on the MRI than they were on CT.
Inverse-optimized plans using either bulk-density-assignment or CT-pixel densities have dose parameters that differ on average by less than 0.5% which was characterized for head and neck IMRT cases by Karotki et al.  We analyzed the dose to 99% and 95% of the MRI-simulation drawn prostate (D99) and (D95) along with the global maximum dose for the MRI- and CT-based plans. Patient results are presented in Table 3. Patient eight could not have bulk-density-assignment because a hip prosthetic did not allow for accurate contouring of femoral heads and surrounding bony tissue. The global maximum dose differed by 1.01% between MRI and CT plans for the patient cohort, and both the D99 and D95 were within 0.2% of the CT-based dose calculation.
CBCT-based treatment localization
Patient 8 was not evaluated for this section due to metal artifacts impacting CBCT image quality. The averages and standard deviations of daily shifts used for treatment localization for each patient are summarized in Table 4. The differences in shift positions for the entire cohort between CBCT-to-CT registration and CBCT-to-MRI registration were −0.15 ± 0.25 cm in the AP direction, 0.07 ± 0.19 cm in the SI direction and −0.01 ± 0.14 cm in the LR direction.
The margins for the CBCT-CT registration were 0.87 cm (AP), 0.64 cm (SI) and 0.69 cm (LR) compared to the CBCT-MRI registration values of 1.06 cm (AP), 0.66 cm (SI) and 0.68 cm (LR). These margin calculations show that CBCT-MRI registration is within two millimeters of those between CBCT-CT. The largest margin difference is in the AP direction. This is due to bowel preparation issues at time of MRI simulation. This issue is expanded on in the discussion section. Imaging modality did not impact daily localization. MRI performs as well as CT as a reference image for registration. Figure 2 shows the capability of T1-MRI for providing an image reference set for daily CBCT by comparing the respective CT/MRI registration to the same CBCT slice. Some CBCT-MRI registrations were performed multiple times in order to measure the reproducibility, with the registration uncertainty (standard deviation) being 0.08, 0.09 and 0.07 cm in the AP, SI and LR directions, respectively.
We found clear visual improvement in anatomical delineation from MRI as seen by the agreement of our volume differences compared to the literature . The penile bulb is often difficult to discern on CT. The difference, both spatial and volumetric, between the MRI volume and the CT volume are stark. Consequently, if using only CT, we would not be able to offer reduction of the collateral radiation damage to these structures during radiation planning and delivery, which may alter the quality of life of patients treated with radiation.
A novel aspect of our study was using MRI as the reference dataset for 3D-3D CBCT-MRI matching for image-guided RT. Such data has not been previously reported for the prostate. Since patients underwent CT and MRI simulation, we could compare the shifts and calculated planning margins obtained with the both datasets. Looking closely at each patient, differences in shifts between MRI and CT registration (Table 4) are consistent overall for each modality. The only registration difference greater than 3 mm between CT and MRI is for patient 7 in the AP direction. This is due to non-compliance of bowel preparation at the MRI and CT simulations. Due to non-compliance, the prostate position relative to bony anatomy changed significantly. Consequently this dramatically impacts registering the CT and MRI images, particularly in the AP direction, which resulted in systematic differences (>6 mm) between MRI and CT-generated shifts. If we exclude this outlier data from the analysis, the margins are nearly identical between MR and CT, as observed in Table 4. A previous 3D-3D CBCT-MRI matching study performed by Buhl et al.  reported on CBCT shifts compared to reference MRI images of the brain and showed mean and standard deviation values of 0.8 ± 0.6 mm, 1.5 ± 1.2 mm and 1.2 ± 1.2 mm differences in the AP, SI and LR directions, respectively. These compare very well to our patient cohort mean and standard deviation values of 1.5 ± 2.5 mm in the AP direction, 0.7 ± 1.9 mm in the SI direction and −0.1 ± 1.4 mm in the LR direction.
One disadvantage of the study was that we could not simulate both MRI and CT during the same appointment session. Since patients had to return for another image acquisition there were instances where rectal/bladder filling was inconsistent between sessions. This created a source of uncertainty during the image registration process since changes in bladder and rectum filling could deform other organs such as the seminal vesicles or prostate. In the clinical implementation of a MRI-only workflow, this source of error would be mitigated.
During our study, we spent a great deal of effort outlining the femoral heads and pelvic bones. MRI simulation and treatment planning relies on correctly delineating all structures since bulk electron densities are assigned. We found normal T1 and T2 image acquisitions provided substandard image quality to properly delineate these bones. bTFE acquisition provided adequate image quality at the expense of additional scan time. Inverting the T1 image grayscale gave a clearer outline of the edge of the femoral heads and other bony structures compared to T1 and T2 images. Image comparisons can be seen in Figure 3. Future utilization of ultra-fast (UTE)  and water/fat separating multi-echo (mDIXON)  sequences are likely necessary for proper identification of bony anatomy. Image registration to the daily CBCT was found to be straightforward and uncomplicated. Aligning the capsule edge, seminal vesicles, rectum interface and the femoral heads was no more difficult than with a CT reference image.
MRI simulation is not without limitations. The first consideration with MRI simulation is whether to use open or closed bore systems. Closed bore systems allow for greater magnet strength, useful for MRI spectroscopy and reducing overall scan time. Open-bore systems consist of 2 horizontal magnets. This open design allows for greater patient comfort and allows simulation therapists increased access to patients, just as in CT simulation. However, as open-bore systems tend to use lower strength magnets which can adversely impact the signal-to-noise ratio. Our open-MRI system uses rigid radiofrequency (RF) coils with integrated solenoid technology and receiver elements perpendicular to the body’s long axis. These are mounted to the patient table to improve the signal-to-noise ratio. Using these coils allows our 1.0 T vertical system to maximize image quality. Enders et al.  found out that in a blinded randomized study evaluating vertical and horizontal field image quality, closed-bore systems had better image quality, though open-bore systems still had an average image quality grade of ‘moderate’, which is acceptable in non-diagnostic environments. The last consideration in bore style is the need to accommodate immobilization devices that are increasingly utilized in radiation oncology departments in order to fixate patients and minimize patient motion during simulation and subsequent treatment. Open-bore systems are more accommodating to a wider arrange of immobilization devices.
Magnetic field distortions may cause geometric aberrations that can result in inaccurate identification of the patient’s skin contour, a crucial component to radiation treatment planning [16,17]. It is important that geometric distortions of each MRI scanner used for simulation be characterized so that associated uncertainties across the FOV are understood [29,30]. Geometric distortions are classified as system distortions or object-related distortions. Geometric distortion is evaluated at our institution on a daily basis for all three axes using a vendor provided 2-dimensional geometric distortion phantom with a field of view (FOV) encompassing imaging sequence (~35 cm x 40 cm). The difference in center of mass position between the acquired image and the ideal position map is measured. Deviations are displayed as isocurves that delineate 2 mm distortion within the primary FOV up to 6 mm of distortion at the FOV edge. The plot is visually inspected to insure the 2 mm isocurve does not advance past the peripheral imaging region. Object-related distortions are another concern and cannot be evaluated with the use of a phantom. Multiple proposals exist for assessing and minimizing object-related distortion [18,19,31-33], but the issue remains uncorrected in clinical settings. Distortions due to tissue susceptibility scale almost linearly with the B0 field strength, and can be considered to be minimal (<0.5 mm) at a field strength of 1.0 T . Relative to CT, uncorrected MRI geometric uncertainties can degrade contouring accuracy while the enhanced soft-tissue contrast improves contouring precision.
Electron density information will always be necessary for accurate dose calculation. For this study only two densities were assigned, but in other areas of the body where air pockets or other tissue densities exist, this could be a bottleneck for treatment planning efficiency without intelligent autosegmentation tools. For MRI-only simulation, we will need to quickly and accurately determine a tissue’s electron density. One possibility is to create so-called “synthetic CTs”, which are currently being researched to infer electron density information directly from MRI [21,34]. For metallic implants such as hip prostheses, CT can use scatter correction algorithms to correct for this presence , whereas MRI is compromised by the metal artifact and could represent a contraindication for MRI-only simulation. Additional safety concerns need to be thoroughly investigated and due diligence performed; any medical product introduced into a high magnetic field environment must be properly vetted. Examples include implanted medical devices and infusion catheters. Our clinic found it necessary to complete MRI screening forms with the patients three times prior to simulation in order to guarantee that incompatible materials do not enter the magnetic field. Patients that are non-compliant in this regard cannot currently be simulated in a MRI-only simulation workflow.
This study confirms that MRI has the potential to be a primary imaging modality for treatment planning simulation, dose calculations, and daily image localization. Readily available radiation oncology directed MR simulation can increase the population of patients that take advantage of more precise tumor delineation given the improved soft tissue contrast with MRI. Consequently it is possible to develop more accurate treatment plans, which have the potential to positively influence patient outcomes. Ultimately this hypothesis needs to be tested, and in this regard, we are currently performing a prospective clinical trial to correlate patient outcomes (toxicities and tumor control) with dose distributions generated using MRI-defined targets.
Written informed consent was obtained from the patients for the publication of this report and any accompanying images.
Datta NR, David R, Gupta RK, Lal P. Implications of contrast-enhanced CT-based and MRI-based target volume delineations in radiotherapy treatment planning for brain tumors. J Cancer Res Ther. 2008;4(1):9–13.
Lemort M, Canizares AC, Kampouridis S. Advances in imaging head and neck tumours. Curr Opin Oncol. 2006;18(3):234–9.
Prabhakar R, Haresh KP, Ganesh T, Joshi RC, Julka PK, Rath GK. Comparison of computed tomography and magnetic resonance based target volume in brain tumors. J Cancer Res Ther. 2007;3(2):121–3.
McLaughlin PW, Troyer S, Berri S, Narayana V, Meirowitz A, Roberson PL, et al. Functional anatomy of the prostate: implications for treatment planning. Int J Radiat Oncol Biol Phys. 2005;63(2):479–91.
Debois M, Oyen R, Maes F, Verswijvel G, Gatti G, Bosmans H, et al. The contribution of magnetic resonance imaging to the three-dimensional treatment planning of localized prostate cancer. Int J Radiat Oncol Biol Phys. 1999;45(4):857–65.
Hricak H. MR imaging and MR spectroscopic imaging in the pre-treatment evaluation of prostate cancer. Br J Radiol. 2005;78 Suppl 2:S103–11.
McLaughlin PW, Narayana V, Meirovitz A, Troyer S, Roberson PL, Gonda Jr R, et al. Vessel-sparing prostate radiotherapy: dose limitation to critical erectile vascular structures (internal pudendal artery and corpus cavernosum) defined by MRI. Int J Radiat Oncol Biol Phys. 2005;61(1):20–31.
Devic S. MRI simulation for radiotherapy treatment planning. Med Phys. 2012;39(11):6701–11.
Mutic S, Dempsey JF. The ViewRay system: magnetic resonance-guided and controlled radiotherapy. Semin Radiat Oncol. 2014;24(3):196–9.
Lagendijk JJ, Raaymakers BW, van Vulpen M. The magnetic resonance imaging-linac system. Semin Radiat Oncol. 2014;24(3):207–9.
Stanescu T, Wachowicz K, Jaffray DA. Characterization of tissue magnetic susceptibility-induced distrotions for MRIgRT. Med Phys. 2012;39(12):7185–93.
Klein E, Chen Z, Chetty IJ, Dogan N. Oncology Scan – Physics. Int J Radiat Oncol Biol Phys. 2012;84(4):871–3.
Karlsson M, Karlsson MG, Nyholm T, Amies C, Zackrisson B. Dedicated magnetic resonance imaging in the radiotherapy clinic. Int J Radiat Oncol Biol Phys. 2009;74(2):644–51.
Lee YK, Bollet M, Charles-Edwards G, Flower MA, Leach MO, McNair H, et al. Radiotherapy treatment planning of prostate cancer using magnetic resonance imaging alone. Radiother Oncol. 2003;66(2):203–16.
Prabhakar R, Julka PK, Ganesh T, Munshi A, Joshi RC, Rath GK. Feasibility of using MRI alone for 3D radiation treatment planning in brain tumors. Jpn J Clin Oncol. 2007;37(6):405–11.
Janke A, Zhao H, Cowin GJ, Galloway GJ, Doddrell DM. Use of spherical harmonic deconvolution methods to compensate for nonlinear gradient effects on MRI images. Magn Reson Med. 2004;52(1):115–22.
Karger CP, Höss A, Bendl R, Canda V, Schad L. Accuracy of device-specific 2D and 3D image distortion correction algorithms for magnetic resonance imaging of the head provided by a manufacturer. Phys Med Biol. 2006;51(12):N253–61.
Pasquier D, Lartiqau E. MRI prostate radiation therapy planning: when the patient distorts his own image (Regarding Lambert et al., Radiother Oncol 2011; 98:330–334). Radiother Oncol. 2012;102(1):163.
Wang H, Balter J, Cao Y. Patient-induced susceptibility effect on geometric distortion of clinical brain MRI for radiation treatment planning on a 3 T scanner. Phys Med Biol. 2013;58(3):465–77.
Hsu SH, Cao Y, Huang K, Feng M, Balter JM. Investigation of a method for generating synthetic CT models from MRI scans of the head and neck for radiotherapy. Phys Med Biol. 2013;58(23):8419–35.
Karotki A, Mah K, Meijer G, Meltsner M. Comparison of bulk electron density and voxel-based electron density treatment planning. J Appl Clin Med Phys. 2011;12(4):97–104.
Johansson A, Karlsson M, Nyholm T. CT subsititute derived from MRI sequences with ultrashort echo time. Med Phys. 2011;38(5):2708–14.
Marks LB, Yorke ED, Jackson A, Ten Haken RK, Constine LS, Eisbruch A, et al. Use of normal tissue complication probability models in the clinic. Int J Radiat Oncol Biol Phys. 2010;76(3 Suppl):S10–9.
van Herk M, Remeijer P, Rasch C, Lebesque JV. The probability of correct target dosage: dose-population histograms for deriving treatment margins in radiotherapy. Int J Radiat Oncol Biol Phys. 2000;47(4):1121–35.
Roach 3rd M, Faillace-Akazawa P, Malfatti C, Holland J, Hricak H. Prostate volumes defined by magnetic resonance imaging and computerized tomographic scans for three-dimensional conformal radiotherapy. Int J Radiat Oncol Biol Phys. 1996;35(5):1011–8.
Buhl SK, Duun-Christensen AK, Kristensen BH, Behrens CF. Clinical evaluation of 3D/3D MRI-CBCT automatching on brain tumors for online patient setup verification – A step towards MRI-based treatment planning. Acta Oncol. 2010;49(7):1085–91.
Helle M, Stehning C, Traughber MS, Schadewaldt N, Schulz H, Renisch S, et al. Comparison of sequences for MR-based cortical bone imaging and tissue classification in the pelvis at 3.0 T with subsequent generation of electron density maps and digitally reconstructed radiographs. ISMRM proceedings 2013. www.ismrm.org/13/session76.htm.
Enders J, Rief M, Zimmermann E, Asbach P, Diederichs G, Wetz C, et al. High-field open versus short-bore magnetic resonance imaging of the spine: a randomized controlled comparison of image quality. PLoS One. 2013;8(12):e83427.
Wang D, Doddrell DM, Cowin G. A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging. Magn Reson Imaging. 2004;22(4):529–42.
Moerland MA, Beersma R, Bhagwandien R, Wijrdeman HK, Bakker CJ. Analysis and correction of geometric distortions in 1.5 T magnetic resonance images for use in radiotherapy treatment planning. Phys Med Biol. 1995;40:1651–4.
Baldwin LN, Wachowicz K, Fallone BG. A two-step scheme for distortion rectification of magnetic resonance images. Med Phys. 2009;36(9):3917–26.
Cao Y, Wang H, Balter J. Characterization of Patient-Induced Geometric Distortions in Clincal Brain MRI on a 3 T MR Simulator. Pract Radiat Oncol. 2013;3(2 Suppl 1):S7.
Wachowicz K, Stanescu T, Thomas SD, Fallone BG. Implications of tissue magnetic susceptibility-related distortion on the rotating magnet in an MR-linac design. Med Phys. 2010;37(4):1714–21.
Kim J, Glide-Hurst C, Doemer A, Wen N, Movsas B, Chetty IJ. Implementation of a novel algorithm for generating synthetic CT images from magnetic resonance imaging datasets for prostate cancer radiation therapy. Int J Radiat Oncol Biol Phys 2014 Nov 7. [Epub ahead of print].
Li H, Noel C, Chen H, Harold Li H, Low D, Moore K, et al. Clinical evaluation of commercial orthopedic metal artifact reduction tool for CT simulations in radiation therapy. Med Phys. 2012;39(12):7507–17.
The authors declare that they have no competing interests.
All authors were all involved in the conception and design of the project. All authors were involved in the data acquisition, analysis or interpretation. All authors were involved in drafting or revising the manuscript. All authors read and approved the final manuscript.