A criterion for the reliable use of MRI-only radiotherapy
© Korsholm et al.; licensee BioMed Central Ltd. 2014
Received: 23 July 2013
Accepted: 23 December 2013
Published: 9 January 2014
MRI-only radiotherapy will eliminate the systematic registration errors introduced when transferring MRI information to the CT. However, challenges concerning the missing information on electron density, necessary for dose calculation and patient setup on bony anatomy are introduced. This study presents a possible statistical approach to evaluate, if deviations based on MRI-only radiotherapy as compared to the CT based radiotherapy are acceptable.
18 head-and-neck, 21 prostate, 10 vesica and 8 pelvic patients were included in the study. Data from each patient contained a CT and a T2-weighted MRI scan, a structure set and a clinically approved CT based treatment plan, which was re-calculated with identical parameters on the density corrected MRI scans. A statistical analysis including a 95% confidence interval was performed in clinically relevant DVH points.
The mean differences in the investigated DVH points were in the order of 1.5% for the PTV and up to 4.2% for organs at risk. In addition, a proposed criterion of 2% dose difference in the PTV coverage for 95% of the patients was fulfilled for all diagnostic groups for a bulk segmented MRI in the DVH points, Dmedian and D2%, while only head-and-neck and prostate further fulfilled the criterion in D98%.
Here, we suggested a method for establishing a reliable use of MRI-only radiotherapy. A population-based study comparing CT based dose calculations with those obtained on a suggested segmentation of MRI should be initiated and acceptable deviations in clinically relevant DVH points should be established. Such a population-based approach could form a part of the clinical commissioning of MRI-only radiotherapy.
The gold standard for radiotherapy (RT) planning is computed tomography (CT). However, little discrimination between the soft tissues is obtained from the CT, since tumour and the organs at risk (OARs) have similar attenuation coefficients. Magnetic resonance imaging (MRI) is therefore increasingly combined with CT for a better delineation of the tumour and OARs. MRI has proven beneficial for multiple treatment sites such as head-and-neck (HN), prostate, pelvic and brain[1–4].
Transferring the MRI delineated structures to the CT scan requires a registration between the two modalities. This introduces a systematic registration error arising from deformable anatomical changes and inconsistent patient setup at the MRI and CT scan. For example, different rectal and bladder filling could cause a relative deformation of the prostate, and an imperfect setup of HN patients could result in different bending of the neck. The net result is that (ideally non-distorted) MRI delineations with one relative relationship between structures are transferred to the anatomy of the CT scan with a different relative relationship between the structures on which the treatment is planned. Systematic registration errors have been quantified for prostate approx. 2 mm (average displacement), rectum approx. 2 mm (average displacement) and head 1.8 ± 2.2 mm (average displacement and one standard deviation).
An alternative to the CT based RT is so-called MRI-only RT, where MRI is the only modality in all steps of the treatment workflow. It has been demonstrated that an MRI-only simulation can reduce spatial systematic uncertainties by 2 mm compared to a CT based workflow for prostate patients. MRI-only RT can potentially lead to a simplified workflow reducing workload and cost, while easing patient discomfort related to multiple scans especially in palliative cases. Additionally, MRI-only based RT enables the use of co-registered functional MRI for assessment of treatment response and adaptive RT, and recently, it introduces the possibility of real-time MRI-guided RT. The main challenges for MRI-only RT are geometrical distortion and the lack of the electron density information needed for dose calculation and setup verification on bone.
The dosimetric impact of performing dose calculations on MRI as compared to CT has previously been reported, although these differences are somewhat inconsistently described. Kristensen et al. observed a 2% difference for brain tumours in the prescribed dose region while larger differences were observed for volumes enclosed by lower iso-dose levels. Jonsson et al. reported mean differences in monitor units (MUs) of 0.2% with a standard deviation of 0.5 for the prescription point for different MRI bulk density corrected geometries of prostate and thorax patients. Lambert et al. reached a mean difference of approximately 2% for the prostate using a similar approach. These studies used conventional (non-IMRT) techniques and did not present any quantitative measure to determine whether the differences were acceptable in a more general setting.
With MRI-only based RT, no comparison of the obtained dose distribution with CT is available, and hence, the order of the dosimetric uncertainty for a patient is unknown and might be unacceptable, e.g. for a serial OAR. One way to evaluate the reliability of MRI-only based RT would be to state a 95%-probability that a patient has an acceptable maximum deviation from the CT calculated dose distribution for clinical relevant DVH points or uniform equivalent dose conversions.
Here, we present a statistical framework as a possible tool to decide whether deviations based on MRI-only RT as compared to standard CT are of an acceptable order. To illustrate the method, data from dose calculations on MRI scans as compared to CT are investigated for three DVH points of the PTV as recommended by ICRU Report 83. The results are reported for multiple treatment sites and techniques.
Retrospective data from 18 HN patients (with oro- and hypo-pharyngeal cancer) treated with a static intensity-modulated radiotherapy (IMRT) setup, 21 prostate, 10 vesica and 8 pelvic (not prostate or vesica) patients treated with volumetric modulated arc therapy (VMAT) were included. Each patient had a CT-scan (Philips Big Bore CT), a T2 weighted MRI (Philips 1 T Panorama), an MRI delineated structure set including targets and OARs and a clinically approved treatment plan (Eclipse version 10.0, Varian Medical Systems).
A homogeneous density assigned MRI (MRIu) where the entire body was assigned the electron density of water (Hounsfield unit (HU) =0).
A bulk heterogeneous density assigned MRI (MRIb) where, in addition, the CT delineated bone was transferred to the MRI and assigned a bone specific electron density.
For the vesica-, prostate- and pelvic patients the assigned electron density was based on the average age of each diagnostic group. For the HN patients, an age-independent electron density was assigned according to skeleton cranium. Table 1 displays an overview of the electron densities and the corresponding calculated HUs.
The calculated Hounsfield units for bone
Electron density [g/cm3]
Calculated HU (age)
1.33 (30 years)
349 (66.8 years*)
1.22 (90 years)
1.33 (30 years)
356 (64.7 years*)
1.22 (90 years)
1.33 (30 years)
309 (78.8 years*)
1.22 (90 years)
An additional bulk heterogeneous density assigned MRI (MRIb,c) for the HN patients where the effect of air cavities was investigated. Here, CT-delineated air cavities were transferred to the MRI and assigned the electron density equal to that of air. The effect of gas pockets in the bowels was not investigated, since these pockets have an irreproducible interfractional variation in amount and location.
Prior to the MRI-only based dose calculation, the CT based treatment plan and structures (bone and cavity) were registered to the corresponding MRI. The body was outlined separately on the MRI to include possible effects of geometrical distortion. The dose distribution was re-calculated on the density corrected MRI keeping the CT planning parameters, i.e. photon fluence, beam energies, angles, MLC control points, monitor units (MUs) etc. unchanged.
To test if such a CT based re-calculation accurately represents the dose differences present when only MR images are available, 10 prostate treatment plans were optimized on MRI according to our clinical standards and re-calculated on CT with unchanged MRI planning parameters. Subsequently, no significant change in the observed dose difference between the CT- and MRI-planning approaches could be identified for the investigated DVH points. Therefore, the re-calculation on MRI with unchanged CT planning parameters was found to be appropriate.
DVH point investigation
The CT and MRI based dose distributions were compared in a number of relevant DVH points including the PTV coverage recommended in ICRU Report 83. The PTV DVH points were the near-maximum absorbed dose (D2%), the near minimum absorbed dose (D98%) and the median absorbed dose (Dmedian). Elective target volumes such as lymph nodes were not included in the investigation. In addition, DVH points for different OARs were compared. The investigated OARs and corresponding DVH points were based on the clinical guidelines used in our department.
We adopt a statistical approach on the dose coverage of a patient cohort as a measure of the reliability of MRI-only RT. We state that 95% of the patients receiving MRI-only RT should have an uncertainty on the dose calculation of the PTV coverage within 2% with respect to that of CT. Such statistical considerations have previously been accepted and incorporated into clinical practice on PTV margin calculations. Ahnesjö et al. presented a table for the individual components contributing to the overall uncertainty of delivering RT. If the uncertainty is 1% and 2% on the CT and MRI dose calculation, respectively, the uncertainty contribution in dose calculation will be 2.2% adding the errors in quadrature. The resulting overall uncertainty would then be 3.3%, which we here considered acceptable.
A statistical analysis comparing the density corrected MRI versus CT based dose calculations for each diagnostic group was carried out to establish 95% confidence intervals (CI) and test the influence of significance (significance level α = 0.05). An univariate analysis of variance (ANOVA) was performed with an one-way two-tailed ANOVA using the statistical software R version 2.11.0. If the ANOVA displayed a significant difference, a pair wise means comparison was performed with a TukeyHSD analysis for further investigation. The TukeyHSD decreases the risk of detecting a false positive result, i.e. significance, as compared to a t-test. Prior to the statistical analysis, the data were tested and found to be independent and approximately normally distributed with constant variances thereby fulfilling the assumptions of the chosen statistical methods. The dosimetric differences will be presented as percentage difference with respect to CT ± two standard deviations thereby describing the approximate 95% CI of the differences. The CIs were based on the standard deviations calculated from the individual diagnostic groups and not from the pooled data of the ANOVA.
The investigation of rectum for the prostate patients was split into two groups since the prescribed dose showed to be significant in a two-way two-tailed ANOVA.
This is a non-interventional retrospective study assuring the quality of treatment using different image modalities. According to the National Health Research Ethics Committee in Denmark (DNVK), such a study does not require an ethics approval.
Statistical results of the DVH point analysis
Volume (# of patients)
1.3 ± 1.4
−0.0002 ± 1.1
−0.03 ± 1.7
1.4 ± 1.3
−0.02 ± 1.0
A(1.0 · 10-7), B(0.0)
2.0 ± 1.7
0.6 ± 1.6
1.9 ± 2.5
0.8 ± 2.5
0.2 ± 3.0
2.2 ± 3.0
0.8 ± 1.5
−0.008 ± 1.3
−0.06 ± 1.6
−0.7 ± 1.6
1.0 ± 1.9
−0.6 ± 1.2
−0.02 ± 3.0
1.6 ± 2.3
−0.006 ± 1.8
−1.0 ± 2.8
1.2 ± 2.0
−0.4 ± 0.9
0.4 ± 4.3
0.8 ± 3.6
−1.4 ± 3.8
−1.4 ± 3.6
Parotid sin (16)
−1.2 ± 6.8
2.5 ± 7.4
−1.3 ± 8.7
Parotid dxt (16)
1.7 ± 7.4
0.2 ± 6.8
1.4 ± 8.9
0.4 ± 1.2
−0.3 ± 1.3
−0.9 ± 4.8
−1.4 ± 4.2
1.1 ± 1.5
0.2 ± 1.8
0.7 ± 1.6
−0.1 ± 1.7
−0.1 ± 3.0
−0.6 ± 3.2
Femur sin (10)
1.4 ± 3.4
−0.6 ± 1.9
Femur dxt (10)
0.7 ± 2.0
−0.7 ± 2.0
−0.5 ± 3.5
−0.9 ± 3.4
4.2 ± 8.0
3.5 ± 7.6
−0.2 ± 1.3
−0.3 ± 1.2
−0.9 ± 2.0
−1.5 ± 2.0
0.8 ± 1.4
0.03 ± 1.6
Femur sin (8)
0.3 ± 1.1
−0.5 ± 1.3
0.07 ± 1.5
−0.2 ± 3.3
Femur dxt (8)
0.1 ± 3.1
−0.2 ± 2.2
0.3 ± 2.5
−0.4 ± 1.4
The results for the vesica patients showed that MRIu differed significantly from CT and MRIb in Dmedian and D2% for PTV, while no statistical differences were seen for the remaining investigated DVH points within this diagnosis.
The density corrected MRIs for the PTV in HN patients showed a significant difference between MRIu and CT, and, MRIb and MRIb,c respectively. The bulk density corrected MRIs with and without air cavities did not differ remarkable. No significant differences were seen for the OARs but a large standard deviation for the parotid glands could be observed (see Discussion).
Neither the PTV nor the OARs showed significance when comparing the density corrected MRIs with the CT for the pelvic patients. Large standard deviations, however, were seen for the intestines (see Discussion).
MRI-only RT receives increasingly attention as MRI becomes more integrated into modern RT[22, 23]. Further, evidence of including MRI into new treatment sites is emerging. Sparing of the hippocampus in cranial irradiation could be one such example that would require an MRI for palliative treatment otherwise traditionally carried out on CT alone. This raises the question of when it is clinically acceptable to perform RT based on MRI as the only modality.
Here, we suggested a method to evaluate the reliability for introducing MRI-only RT. A population based study comparing CT based dose calculations with those carried out on the segmented MRI of choice (e.g.) should be initiated. Acceptable deviations in clinically relevant DVH points should be established. In this study we accepted a 2% deviation in PTV coverage for 95% of the patients and assumed an uncertainty of 1% on the CT based calculations. If a more conservative uncertainty estimate from CT based calculations is assumed, e.g. about 3% as estimated by ICRP more than 10 years ago, a more restricted uncertainty requirement should be imposed on the MRI deviation ultimately to reach an overall uncertainty of about 4%. We obtained a conservative CI for the MRIb medulla of 4.2%, which could be critical near tolerance doses for myelopathy. Tighter acceptance criteria, e.g. a 1% deviation for a 99% CI could be considered on normal tissue complication probabilities that receive a high priority.
The size of the patient cohort will influence the statistical power and therefore the credibility by which a segmentation strategy fulfill a criterion for MRI-only RT. The size of the patient cohort included in the study can influence the ability to estimate a dosimetric difference with acceptable statistical power. A power analysis on the available CT dose statistics in the relevant DVH points should therefore be carried out. For example, the D98% mean and standard deviation for the prostate patients CT based dose calculation were 1.94 and 0.03 Gy, respectively. With 21 patients, the ability to detect a mean difference of at least 2% could be made with a power of 96%. As seen in the middle panel of Figure 3, such a calculation did not indicate whether a segmentation criteria for MRI-only RT was fulfilled for 95% of the patients (bulk passes and unit does not) but only support the credibility of the statistical analysis. A limitation in the proposed approach was that only a measure of similarity in the DVH points was described without any spatial information of where the differences were located.
Large CIs of 8 and 7% were found for the intestines and parotid glands, respectively. The former could be due to irreproducible gas pockets, which were not transferred to the MRI from CT, and were considered less critical. For the latter, inspections of the CT and MRI showed that the main differences in the body-outline were in the area around the ears, which potentially influenced the calculated dose to the parotid glands. Therefore, this area needs special attention during MRI acquisition and should possibly be corrected for geometrical distortion. Generally, OARs are in lower dose regions and small differences in absolute dose would appear as larger percentage differences. Therefore, the standard deviations in general were found to be larger for the OARs than the PTV DVH points.
The investigation of the HN patients showed similar results when comparing the bulk density corrected MRIs with and without air cavities, this may be related to the fact that no patients with nasopharyngeal cancer, where the dose distributions are expected to be more influenced by air cavities, were included. Hence, the effect of air cavities in nasopharyngeal cancer patients should be investigated further.
For our bulk density correction the majority of the patients showed an uncertainty of 2% or less, on the PTV coverage. The patients that did not fulfil the criteria and their corresponding dose distribution were investigated. For the D98%, 1 prostate (1/21), 1 HN (1/18), 2 vesica (2/10) and 2 pelvic (2/8) patients had deviations above 2% which is consistent with the diagnosis passing the criteria in the middle panel of Figure 3. For the abdominal area, these patients had multiple gas pockets leading to dose valleys in the corresponding MRI dose distributions. The position of the tumor bed of the (post-surgery) HN patient was close to the trachea. This was uncorrected for in the MRIb segmentation with a percentage deviation of 2.3% whereas the deviation was 0.6% in the MRIb,c segmentation. These findings could speak in favor of including air segmentation for all diagnostic groups although this was not the strategy chosen in this study. One vesica patient did not meet the criteria D2%, which could be explained by an abnormally large change in body outline on MRI as compared to CT due to an anatomic deformation (interfractional variation).
This study separated the results for the different diagnostic groups why the impact of the different delivery techniques could not be estimated, as each diagnostic group was treated with the same delivery technique. Since our main results are consistent with those reported previously obtained with 3D conformal RT, the impact is expected to be minimal.
We have suggested a method for establishing a reliable use of MRI-only radiotherapy. A population-based study comparing CT based dose calculations with those obtained on a suggested segmentation of MRI should be initiated and acceptable deviations in clinically relevant DVH points should be established. Such a population-based approach could form a part of the clinical commissioning of MRI-only radiotherapy.
This work was supported by a research grant from Varian Medical Systems, Inc.
- Rasch C, Barillot I, Remeijer P, Touw A, van Herk M, Lebesque JV: Definition of the prostate in CT and MRI: a multi-observer study. Int J Radiat Oncol Biol Phys 1999, 43: 57-66. 10.1016/S0360-3016(98)00351-4View ArticlePubMedGoogle Scholar
- Barillot I, Reynaud-Bougnoux A: The use of MRI in planning radiotherapy for gynaecological tumours. Cancer Imaging 2006, 6: 100-106. 10.1102/1470-7330.2006.0016PubMed CentralView ArticlePubMedGoogle Scholar
- Thiagarajan A, Caria N, Schoder H, Iyer NG, Wolden S, Wong RJ, Sherman E, Fury MG, Lee N: Target volume delineation in oropharyngeal cancer: impact of PET, MRI, and physical examination. Int J Radiat Oncol Biol Phys 2010, 83: 220-227.View ArticleGoogle Scholar
- Aoyama H, Shirato H, Nishioka T, Hashimoto S, Tsuchiya K, Kagei K, Onimaru R, Watanabe Y, Miyasaka K: Magnetic resonance imaging system for three-dimensional conformal radiotherapy and its impact on gross tumor volume delineation of central nervous system tumors. Int J Radiat Oncol Biol Phys 2001, 50: 821-827. 10.1016/S0360-3016(01)01598-XView ArticlePubMedGoogle Scholar
- Roberson PL, McLaughlin PW, Narayana V, Troyer S, Hixson GV, Kessler ML: Use and uncertainties of mutual information for computed tomography/magnetic resonance (CT/MR) registration post permanent implant of the prostate. Med Phys 2005, 32: 473-482. 10.1118/1.1851920View ArticlePubMedGoogle Scholar
- Dean CJ, Sykes JR, Cooper RA, Hatfield P, Carey B, Swift S, Bacon SE, Thwaites D, Sebag-Montefiore D, Morgan AM: An evaluation of four CT-MRI co-registration techniques for radiotherapy treatment planning of prone rectal cancer patients. Br J Radiol 2012, 85: 61-68. 10.1259/bjr/11855927PubMed CentralView ArticlePubMedGoogle Scholar
- Ulin K, Urie MM, Cherlow JM: Results of a multi-institutional benchmark test for cranial CT/MR image registration. Int J Radiat Oncol Biol Phys 2010, 77: 1584-1589. 10.1016/j.ijrobp.2009.10.017PubMed CentralView ArticlePubMedGoogle Scholar
- Nyholm T, Nyberg M, Karlsson MG, Karlsson M: Systematisation of spatial uncertainties for comparison between a MR and a CT-based radiotherapy workflow for prostate treatments. Radiat Oncol 2009, 4: 54. 10.1186/1748-717X-4-54PubMed CentralView ArticlePubMedGoogle Scholar
- 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: 644-651. 10.1016/j.ijrobp.2009.01.065View ArticlePubMedGoogle Scholar
- Nielsen T, Mouridsen K, Maxwell RJ, Stødkilde-Jørgensen H, Østergaard L, Horsman MR: Segmentation of dynamic contrast enhanced magnetic resonance imaging data. Acta Oncol 2008, 47: 1265-1270. 10.1080/02841860802277489View ArticlePubMedGoogle Scholar
- Hu Y, Green OP, Parikh P, Olsen J, Mutic S: TH-E-BRA-07: Initial Experience with the ViewRay System-Quality Assurance Testing of the Imaging Component. Med Phys 2012, 39: 4013.View ArticleGoogle Scholar
- Jonsson JH, Karlsson MG, Karlsson M, Nyholm T: Treatment planning using MRI data: an analysis of the dose calculation accuracy for different treatment regions. Radiat Oncol 2010, 5: 62. 10.1186/1748-717X-5-62PubMed CentralView ArticlePubMedGoogle Scholar
- Kristensen BH, Laursen FJ, Logager V, Geertsen PF, Krarup-Hansen A: Dosimetric and geometric evaluation of an open low-field magnetic resonance simulator for radiotherapy treatment planning of brain tumours. Radiother Oncol 2008, 87: 100-109. 10.1016/j.radonc.2008.01.014View ArticlePubMedGoogle Scholar
- Lambert J, Greer PB, Menk F, Patterson J, Parker J, Dahl K, Gupta S, Capp A, Wratten C, Tang C, et al.: MRI-guided prostate radiation therapy planning: Investigation of dosimetric accuracy of MRI-based dose planning. Radiother Oncol 2011, 98: 330-334. 10.1016/j.radonc.2011.01.012View ArticlePubMedGoogle Scholar
- International Commission on Radiation Units and Measurements (ICRU): Prescribing, Recording, and Reporting Photon-Beam Intensity-Modulated Radiation Therapy (IMRT). ICRU Report 83. Bethesda, MD: International Commission on Radiation Units and Measurements; 2011.Google Scholar
- International Commission on Radiation Units and Measurements (ICRU): Photon, Electron, Proton and Neutron Interaction Data for Body Tissues. ICRU Report 46. Bethesda, MD: International Commission on Radiation Units and Measurements; 1968.Google Scholar
- International Commission on Radiation Units and Measurements (ICRU): Use of Computers in External Beam Radiotherapy Procedures with High-Energy Photons and Electrons. ICRU Report 42. Bethesda, MD: International Commission on Radiation Units and Measurements; 1949.Google Scholar
- 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: 1121-1135. 10.1016/S0360-3016(00)00518-6View ArticlePubMedGoogle Scholar
- Ahnesjo A, Aspradakis MM: Dose calculations for external photon beams in radiotherapy. Phys Med Biol 1999, 44: R99-R155. 10.1088/0031-9155/44/11/201View ArticlePubMedGoogle Scholar
- R Development Core Team: R: A lnguage and environment for statistical computing. R Foundation for Statistical Computing 2009. http://www.R-project.orgGoogle Scholar
- Crawley MJ: The R Book. Chichester: John Wiley & Sons; 2007.View ArticleGoogle Scholar
- Klein EE, Das IJ, Dong L, Rivard MJ: Oncology scan—improvements in dose calculation, deformable registration, and MR-guided radiation delivery. Int J Radiat Oncol Biol Phys 2013, 86: 395-397. 10.1016/j.ijrobp.2013.02.017View ArticleGoogle Scholar
- Klein EE, Chen Z, Chetty IJ, Dogan N: Oncology scan—physics. Int J Radiat Oncol Biol Phys 2012, 84: 871-873. 10.1016/j.ijrobp.2012.08.024View ArticleGoogle Scholar
- Gondi V, Hermann BP, Mehta MP, Tome WA: Hippocampal dosimetry predicts neurocognitive function impairment after fractionated stereotactic radiotherapy for benign or low-grade adult brain tumors. Int J Radiat Oncol Biol Phys 2011, 85: 348-354.View ArticleGoogle Scholar
- Johansson A, Karlsson M, Nyholm T: CT substitute derived from MRI sequences with ultrashort echo time. Med Phys 2011, 38: 2708-2714. 10.1118/1.3578928View ArticlePubMedGoogle Scholar
- Valentin J: Prevention of accidental exposures to patients undergoing radiation therapy. Ann ICRP 2000, 30: 57-61. 10.1016/S0146-6453(01)00043-4View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.