- Open Access
Dosimetric evaluation of synthetic CT relative to bulk density assignment-based magnetic resonance-only approaches for prostate radiotherapy
Radiation Oncology volume 10, Article number: 239 (2015)
Magnetic resonance imaging (MRI) has been incorporated as an adjunct to CT to take advantage of its excellent soft tissue contrast for contouring. MR-only treatment planning approaches have been developed to avoid errors introduced during the MR-CT registration process. The purpose of this study is to evaluate calculated dose distributions after incorporating a novel synthetic CT (synCT) derived from magnetic resonance simulation images into prostate cancer treatment planning and to compare dose distributions calculated using three previously published MR-only treatment planning methodologies.
An IRB-approved retrospective study evaluated 15 prostate cancer patients that underwent IMRT (n = 11) or arc therapy (n = 4) to a total dose of 70.2-79.2 Gy. Original treatment plans were derived from CT simulation images (CT-SIM). T1-weighted, T2-weighted, and balanced turbo field echo images were acquired on a 1.0 T high field open MR simulator with patients immobilized in treatment position. Four MR-derived images were studied: bulk density assignment (10 HU) to water (MRW), bulk density assignments to water and bone with pelvic bone values derived either from literature (491 HU, MRW+B491) or from CT-SIM population average bone values (300 HU, MRW+B300), and synCTs. Plans were recalculated using fixed monitor units, plan dosimetry was evaluated, and local dose differences were characterized using gamma analysis (1 %/1 mm dose difference/distance to agreement).
While synCT provided closest agreement to CT-SIM for D95, D99, and mean dose (<0.7 Gy (1 %)) compared to MRW, MRW+B491, and MRW+B300, pairwise comparisons showed differences were not significant (p < 0.05). Significant improvements were observed for synCT in the bladder, but not for rectum or penile bulb. SynCT gamma analysis pass rates (97.2 %) evaluated at 1 %/1 mm exceeded those from MRW (94.7 %), MRW+B300 (94.0 %), or MRW+B491 (90.4 %). One subject’s synCT gamma (1 %/1 mm) results (89.9 %) were lower than MRW (98.7 %) and MRW+B300 (96.7 %) due to increased rectal gas during MR-simulation that did not affect bulk density assignment-based calculations but was reflected in higher rectal doses for synCT.
SynCT values provided closest dosimetric and gamma analysis agreement to CT-SIM compared to bulk density assignment-based CT surrogates. SynCTs may provide additional clinical value in treatment sites with greater air-to-soft tissue ratio.
Radiation therapy treatment planning was developed using computed tomography (CT) as its base imaging modality due to accurate geometric fidelity and the straightforward conversion from linear attenuation coefficients to electron density values. However, a major drawback in CT is its poor soft tissue contrast that makes it difficult to accurately identify and contour soft tissue structures. Conversely, magnetic resonance imaging (MRI) provides excellent soft tissue contrast. Thus, efforts have been made to incorporate MRI into the treatment planning process. MR images would first be registered to the CT simulation image (CT-SIM), and MR-contoured structures would then be transferred onto the CT-SIM for treatment planning. Unfortunately, the registration process introduces additional systematic uncertainties (~1-2 mm for pelvis)  that would propagate throughout the treatment planning workflow. Moreover, having two simulation modalities can be cost-prohibitive while introducing additional burden on the clinical workflow. Therefore, interest has grown in developing an MR-only workflow for radiotherapy treatment planning [2, 3].
While the benefits of MR imaging (i.e. excellent soft tissue contrast and lack of ionizing radiation) are substantial, several logistical difficulties must be overcome before MR-only simulation for radiotherapy would be practical clinically. The first major issue is geometric distortion, which is typically categorized as either machine-specific or patient-specific. The main machine-specific sources of distortion are inhomogeneity of the main magnetic field and non-linearity in the magnetic field produced by gradient coils, while magnetic susceptibility-induced distortion is the dominant form of patient-specific errors. Many papers have been published characterizing the magnitude of these distortions and the post-processing of the image needed to correct for them . Through various correction schemes and the use of MR sequences such as 3D turbo spin echo sequences, it has been reported that distortions can be reduced to within acceptable tolerance for radiotherapy purposes for many treatment sites .
Another significant limitation of MRI is that there is no fundamental relationship between MRI image intensity values and the electron density information required by treatment planning systems to accurately account for heterogeneity within a patient for accurate dose calculation. The most straightforward method to overcome this obstacle is to assume no significant heterogeneities exist and assign a bulk electron density value equivalent to water within the entire body contour [5, 6]. However, this approach introduces uncertainties in the prostate of up to 2.5 % . Variations of this method have been introduced where one or more structures that introduce large heterogeneities (e.g. bone or air) are segmented and assigned appropriate relative electron density values [3, 6, 7], which reduced the dose uncertainty to less than 1.5 %. Other methods that have been introduced are based on the use of atlases  or statistical models [9, 10]. Atlas-based methods may struggle with atypical anatomy, such as cases with local recurrence after radical prostatectomy. Presently, studies incorporating statistical methods have focused mainly in the brain region and rely on sequences that may have difficulty with larger field of view (FOV) body sites such as the pelvis.
We recently introduced a novel, voxel-based, weighted-summation method for generating synthetic CTs (synCTs) from MRI images for male pelvis anatomy . While preliminary dosimetric comparisons were made between synCT and CT-SIM, we build on our initial investigation by evaluating synCT performance relative to three other established approaches for MR-only treatment planning of prostate cancer. In this paper, the use of synCT for dose calculation was compared to the use of bulk density assignment methods in order to elucidate situations where implementation of the synCT algorithm may be beneficial.
Retrospective analysis was performed for 15 early stage (T1 or T2, N0, M0) prostate cancer patients with a median patient age of 71 years (range: 53–96) that were enrolled in an institutional review board-approved study wherein MR-simulation was performed as an adjunct to the CT-simulation process. The planning target volume (PTV) was composed of the prostate and proximal seminal vesicles for 13 patients, while the two remaining patients were treated for local recurrence in the prostate/seminal vesicle bed after radical prostatectomy. All patients underwent radiotherapy delivered using either intensity-modulated radiation therapy (IMRT, n = 11) or volume-modulated arc therapy (VMAT, n = 4) to a total dose of 70.2-79.2 Gy. All patients received the same preparation instructions for both imaging sessions: full bladder and empty rectum, which is consistent with our clinical practice .
A Brilliance Big Bore (Philips Health Care, Cleveland, OH) CT scanner was used to acquire pelvis images using the following parameters: 140 kVp, 500 mAs, 512x512 in-plane image dimensions, 1.28x1.28 mm2 in-plane spatial resolution, and 3 mm slice thickness.
A 1.0 T Panorama High Field Open (Philips Medical Systems, Best, Netherlands) MR simulator was used with rigid, solenoid-based body coils and a flat tabletop insert (Civco, Orange City, IA). Patients were aligned and leveled to their tattoos from CT simulation, and the same immobilization devices (banded feet and shaped foam pad for legs) were used. T1-weighted fast field echo, T2-weighted turbo spin echo, and balanced turbo field echo images were acquired for each patient.
Sequences were collected using the parameters for repetition time (TR), echo time (TE), flip angle (α), FOV, and image resolution provided in Table 1. Total scan time to collect these sequences was approximately 18 min. An inverse T1 image was also generated by subtracting the T1 image from the intensity value below which fell 95 % of the area under the curve of the intensity value histogram for the T1 image. Inverse T1 images provided high intensity values in bone regions for the synCT algorithm.
CT surrogate generation
Four different established approaches of defining MR-derived CT surrogates were generated for comparison with the CT-SIM. The first was a homogeneous electron density relative to water assignment of 1.0 to the entire external contour (MRW) . The second approach, defined as MR (water + bone = 491 HU or MRW+B491) included two segments: relative density assignment of 1.0 for the external contour and relative electron density assignment of 1.27  (equivalent to 491 HU) to bone that was manually contoured on the T2-weighted MR image. While assigning 491 HU to bone has been used in the literature, it has been observed that this value may be too high for the femoral heads  since the higher proportion of cancellous (i.e. spongy) bone within the femoral head would reduce its relative electron density. Consequently, it has been found that calculating the average CT value within the femoral heads for the patient population provides more accurate relative electron density values that produce dose distributions closer to original CT-based treatment plans . We followed this workflow by contouring the femoral heads for each patient and calculating the average CT value over the full set of 15 patients, which yielded a value of ~300 HU. Therefore, the third approach included a relative density assignment of 1.0 to the external contour as well as an assignment of 300 HU (relative electron density ≈ 1.14) to the femoral heads (water + bone = 300 HU or MRW+B300). The final approach was to generate synCTs via the method described in detail in a previous publication . Briefly, manual segmentation of bone was followed by automatic segmentation of four other predetermined classes (air, soft tissue, fat, and fluid) using k-means clustering and morphological operations. Each synCT voxel value was then calculated as a summation of intensity values from the acquired MR images that had been weighted by sequence- and class-specific factors. The various MR-derived surrogates are shown in Fig. 1.
All patients were treated using CT-SIM based plans created in the Eclipse® treatment planning system (Varian Medical Systems, Palo Alto, CA) with physician-contoured target and organ-at-risk (OAR) volumes. Treatments were delivered using either IMRT with seven or nine fields delivering 1.8-2.0 Gy/fraction (dose range: 75.6-79.2 Gy) or VMAT using a single arc to deliver 1.8-2.0 Gy/fraction (dose range: 70–78 Gy). Once generated, each CT surrogate was interpolated onto the original CT-SIM image grid and imported into Eclipse®. The synCT was rigidly registered to the CT-SIM image, and all other CT surrogates were registered by applying those same transformation coordinates derived from the registration of the synCT image. Original CT-SIM treatment plans were copied onto each CT surrogate, and dose was recalculated using original plan parameters. Target and OAR contours were then transferred from the CT-SIM to each CT surrogate for analysis.
Standard dose-volume histogram (DVH) metrics were evaluated for the target and OARs. PTVs were compared using D99, D95, and mean dose to prostate where Dx is equal to the dose delivered to x% of the structure volume. Bladders, prostates, and penile bulbs were analyzed using dosimetric indices defined in QUANTEC guidelines [15–17]. Namely, D15, D25, and D35 were calculated for the bladder  with observed clinical symptoms as an endpoint. D15, D25, and D35 were calculated for the rectum based on their use in evaluating the risk for late rectal toxicity . Finally, D90 for the penile bulb, which is used to assess the risk of radiotherapy-induced erectile dysfunction , was calculated. Absolute differences between dosimetric values calculated using CT-SIM, taken as ground truth, and those calculated using each of the MR-derived CT surrogates were determined. Relative agreement with CT-SIM-derived dose values was evaluated through a paired comparison between each set of pseudo-CT images using nonparametric signed rank tests with the significance level set to p = 0.05. Moreover, 2D gamma analysis (evaluated at 1 %/1 mm and 2 %/2 mm dose difference/distance to agreement) was performed to compare axial dose distributions at isocenter for each of the MR-derived treatment plans to those from the CT-SIM-based plans.
Table 2 compares the population average values along with standard deviation and 95 % confidence intervals for selected DVH metrics calculated using CT-SIM-based treatment plans with those of each CT surrogate-based plans. While synCT yielded the smallest dosimetric differences from CT-SIM for target D95, D99, and mean dose (<0.7 Gy (1 %)) compared to MRW, MRW+B491, and MRW+B300, pairwise comparisons showed that dose reduction was only statistically significant relative to MRW+B491. Figure 2 shows that the synCT tended to have closer agreement overall with CT-SIM values, but that patient-specific results were variable. In addition, using population average values for bone assignment (i.e. MRW+B300) provided significant improvements relative to MRW+B491 for all metrics. For bladder, SynCT-derived dose values displayed closest average agreement with CT-SIM dose values for all metrics, with small (<1 %) but significant improvements relative to MRW, MRW+B491, and MRW+B300. For the rectum and penile bulb, all MR-derived images showed good agreement (difference <1 %) with original CT-SIM calculated values, and no significant differences were observed between synCT, MRW, and MRW+B300. MRW+B491 performed significantly worse for all metrics. Figure 3 displays DVHs for Patient 5, whose OAR results showed average agreement with CT-SIM values, and Patient 8, who displayed the largest deviation from those values (up to 3 %, 2 %, and 3.5 % for bladder, rectum, and penile bulb, respectively). For Patient 8, the bulk density assignment-based approaches produced dose values lower than those produced for CT-SIM because of a higher than normal amount of fatty tissue, where CT-SIM values were lower than bulk assigned HU values. Conversely, while synCT tissue HU values were fairly close to those in CT-SIM, synCT bone HU values tended to be lower for this patient, leading to higher calculated dose values.
Overall, synCT gamma analysis pass rates at 1 %/1 mm (97.2 %) exceeded those for MRW (94.7 %), MRW+B491 (90.4 %), and MRW+B300 (94.0 %) with individual patient results at both 2 %/2 mm and 1 %/1 mm for each CT surrogate shown in Table 3. Generally, gamma analysis conducted at 2 %/2 mm was not sensitive enough to reveal significant differences between CT surrogates (pass rates ranging from 97.8 % (MRW+B491) to 99.8 % (synCT)). Figure 4 displays a typical gamma analysis case where the synCT showed slightly better agreement than alternative CT surrogate methods. One exception to this trend occurred for Patient 11, whose gamma analysis results at 1 %/1 mm were lower for synCT (89.9 %) than for MRW (98.7 %) or MRW+B300 (96.7 %). As Fig. 5a illustrates, there was a large increase in rectal gas during MR simulation relative to that present during CT simulation. While this did not affect the bulk density assignment methodologies tested, Fig. 5 shows that the change in air volume was propagated to the generated synCT (Fig. 5e). The increased air volume led to higher calculated doses for the synCT-derived treatment plan while the bulk density assignment-based methods remained unaffected (Fig. 5g). For this case, synCT rectal doses were systematically higher than all other approaches: MRW by 1.6-2.6 Gy (4.3-5.6 %), MRW+B300 by 1.8-3.0 Gy (4.6-5.8 %), and MRW+B491 by 2.0-3.5 Gy (5.0-6.2 %).
This work sought to incorporate synCTs and previously published bulk density assignment-based CT surrogates into the development of radiotherapy treatment plans and to evaluate the resulting dose distributions in order to determine the impact of integrating synCT into the treatment planning workflow. Overall, target and OAR dose metrics and gamma analysis used for comparison revealed that synCT-derived treatment plans tended to provide closer agreement to CT-SIM than bulk density assignment-based approaches. While 3D gamma analysis has recently emerged as a tool for evaluating the entire irradiated volume of both target and OARs, 2D gamma analysis has been reported to yield more stringent results  and is consistent with our clinical practice. The most significant OAR differences were observed for the bladder. Consistent with the literature , the MRW+B491 tended to perform significantly worse than other MR-derived images due to its high bone relative electron density values that yielded slight underestimations of dose to target and OARs relative to the CT-SIM plan. Use of population average-derived relative electron density values for MRW+B300 provided better agreement with original CT-SIM plan dose distributions than MRW+B491. However, the improvement, though statistically significant, was minimal (~0.5 %), which is consistent with previous reports (~0.6 %) .
Common external contours were used that may mask dosimetric differences between dose distributions calculated using CT-SIM and those using MR-derived CT substitutes due to changes in treatment beam path length. These were added to correct for changes in patient geometry between different patient setups, particularly for cases having a time lapse of up to two weeks between CT and MR simulation. This effect was assumed to be small. Consistent with that assumption, the number of added voxels amounted to <5 % of the whole, and these were mainly concentrated around the abdomen and out of the path of the treatment beams.
With the exception of MRW, which does not take into account tissue heterogeneities, all CT replacement methods that were evaluated required manual segmentation of bone. Manual bone segmentation is time-consuming, making its clinical implementation impractical. Currently, different avenues for automatically segmenting bone in MR images of the prostate are under investigation. Both atlas-based  and statistical model-based approaches [9, 10] have been proposed. In atlas-based methods, probabilistic methods are employed to automatically segment structures of interest. Atlas-based methods rely on the accuracy of the implemented multimodal deformable image registration algorithm and, therefore, appear to be prone to the same MR-CT registration errors that implementation of an MR-only treatment planning workflow attempts to avoid. Statistical methods utilize magnitude information from several MR sequences to identify clusters associated with different segment types in the joint histogram and to determine the probability that a voxel belongs to each of the clusters. Like air, bone appears as a signal void in conventional MR images because of its low proton density and extremely short T2 signal lifetime . Therefore, bone tends to be segmented with air. While ultrashort echo time (UTE) sequences may help to obtain a discernible signal from bone in order to separate bone from air, UTE has mainly been focused on the head region because of difficulties in attaining acceptable image quality at larger FOVs.
While reducing the time required for image post-processing and synCT generation is important, it is also important to reduce the scan time for the patient. The stated time to acquire the requisite images for our synCT workflow is ~18 min for a 1.0 T scanner. By contrast, a CT-SIM scan can be obtained in <1 min. Long scan times may cause degradation of image quality due to patient motion or variation in patient internal anatomy (e.g. changes in bladder and rectal filling). Scan times would be reduced at higher field strengths where increases in signal-to-noise ratio would also enable more aggressive acceleration factors, thereby reducing acquisition time. Future work will explore reducing the number of image sets needed to generate synCTs without introducing significant errors.
Variation in the level of rectal and bladder filling has been shown to deform the shape and to shift the position of the prostate and seminal vesicles by as much as 0.7 cm , leading to a reduction in local control for prostate cancer radiotherapy . Therefore, the level of bladder and rectal filling must be reproducible between simulation and delivery of treatment. There is an increased risk that the level of rectal and bladder filling may change during an MR session simply due to the extended scan time for MR compared to that of CT. To mitigate this risk, it was recommended that steps be taken to ensure the patient has a full bladder  and empty rectum  prior to simulation such as regular patient instruction before and during the course of treatment. In extreme cases, administration of an enema prior to simulation may be warranted .
Results from Patient 11 indicate that synCT may provide greater value in regions with a higher air-tissue ratio, such as the brain or head and neck where airways are more prevalent. Manual segmentation in these regions is more difficult and time-intensive than for the prostate and would only further strain clinical workflow. For future implementations of the synCT workflow in the brain to be viable clinically, some method to automatically and accurately differentiate air from bone would be necessary. A promising method for accurately segmenting air separately from bone using UTE phase information is currently under development in our group.
Overall, synCT-based treatment plans provided dosimetric and gamma analysis values in close agreement to original CT-based treatments plans consistent with previously published bulk density assignment-based substitute CT images, with significant improvement relative to those other substitute CT images observed for bladder.
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.
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–51.
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.
Liney GP, Moerland MA. Magnetic resonance imaging acquisition techniques for radiotherapy planning. Semin Radiat Oncol. 2014;24:160–8.
Chen L, Price Jr RA, Nguyen TB, Wang L, Li JS, Qin L, et al. Dosimetric evaluation of MRI-based treatment planning for prostate cancer. Phys Med Biol. 2004;49:5157–70.
Lambert J, Greer PB, Menk F, Patterson J, Parker J, Dahl K, et al. MRI-guided prostate radiation therapy planning: investigation of dosimetric accuracy of MRI-based dose planning. Radiother Oncol. 2011;98:330–4.
Chin AL, Lin A, Anamalayil S, Teo BK. Feaasibility and limitations of bulk density assignment in MRI for head and neck IMRT treatment planning. J Appl Clin Med Phys. 2014;15:4851.
Dowling JA, Lambert J, Parker J, Salvado O, Fripp J, Capp A, et al. An atlas-based electron density mapping method for magnetic resonance imaging (MRI)-alone treatment planning and adaptive MRI-based prostate radiation therapy. Int J Radiat Oncol Biol Phys. 2012;83:e5–e11.
Johansson A, Karlsson M, Nyholm T. CT substitute derived from MRI sequences with ultrashort echo time. Med Phys. 2011;38:2708–14.
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 radiation therapy. Phys Med Biol. 2013;58:8419–35.
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 data sets for prostate cancer radiation therapy. Int J Radiat Oncol Biol Phys. 2015;91:39–47.
Wen N, Kumarasiri A, Nurushev T, Burmeister J, Xing L, Lui D, et al. An assessment of PTV margin based on actual accumulated dose for prostate cancer radiotherapy. Phys Med Biol. 2013;58:7733–44.
Korhonen J, Kapanen M, Keyrilainen J, Seppala T, Tuomikoski L, Tenhunen M. Absorbed doses behind bones with MR image-based dose calculations for radiotherapy treatment planning. Med Phys. 2013;40:011701.
Hu Y, Zhao W, Du D, Wooten HO, Olsen JR, Gay HA, et al. Magnetic resonance imaging-based treatment planning for prostate cancer: Use of population average tissue densities within the irradiated volume to improve plan accuracy. Pract Radiat Oncol. 2015; epub.
Viswanathan AN, Yorke ED, Marks LB, Eifel PJ, Shipley WU. Radiation dose-volume effects of the urinary bladder. Int J Radiat Oncol Biol Phys. 2010;76:S116–22.
Michalski JM, Gay H, Jackson A, Tucker SL, Deasy JO. Radiation dose-volume effects in radiation-induced rectal injury. Int J Radiat Oncol Biol Phys. 2010;76:S123–9.
Roach III M, Nam J, Gagliardi G, El Naga I, Deasy JO, Marks LB. Radiation dose-volume effects and the penile bulb. Int J Radiat Oncol Biol Phys. 2010;76:S130–4.
Pulliam KB, Huang JY, Howell RM, Followill D, Bosca R, O’Daniel J, et al. Comparison of 2D and 3D gamma analyses. Med Phys. 2014;41:021710.
Delso G, Carl M, Wiesinger F, Sacolick L, Porto M, Hüllner M, et al. Anatomic evaluation of 3-dimensional ultrashort-echo-time bone maps for PET/MR attenuation correction. J Nucl Med. 2014;55:780–5.
Roeske JC, Forman JD, Mesina CF, He T, Pelizzari CA, Fontenla E, et al. Evaluation of changes in the size and location of the prostate, seminal vesicles, bladder, and rectum during a course of external beam radiation therapy. Int J Radiat Oncol Biol Phys. 1995;33:1321–9.
de Crevoisier R, Tucker SL, Dong L, Mohan R, Cheung R, Cox JD, et al. Increased risk of biochemical and local failure in patients with distended rectum on the planning CT for prostate cancer radiotherapy. Int J Radiat Oncol Biol Phys. 2005;62:965–73.
Jain S, Loblaw DA, Morton GC, Danjoux C, Szumacher E, Chu W, et al. The effect of radiation technique and bladder filling on the acute toxicity of pelvic radiotherapy for localized high risk prostate cancer. Radiother Oncol. 2012;105:193–7.
Dr. Glide-Hurst is partially funded by an internal mentored grant and has received funding unrelated to this work from Philips Healthcare. Dr. Chetty and Dr. Movsas hold research agreements with Philips Healthcare and Varian Medical Systems for research outside the submitted work.
Henry Ford Health System holds research agreements with Philips Healthcare (Best, Netherlands) and Varian Medical Systems (Palo Alto, CA).
JK, KG, IJC, and CGH were involved in the conceptual design of the project. JK, KG, and CGH developed treatment plans and collected data. KL, BM, and MSS were involved in patient care, OAR and target contouring, and plan approval. JK and CGH analyzed the datasets, and LS performed statistical analysis on the resulting data. All authors contributed to writing, editing, and giving final approval to the manuscript.
About this article
Cite this article
Kim, J., Garbarino, K., Schultz, L. et al. Dosimetric evaluation of synthetic CT relative to bulk density assignment-based magnetic resonance-only approaches for prostate radiotherapy. Radiat Oncol 10, 239 (2015) doi:10.1186/s13014-015-0549-7
- MR simulation
- Synthetic CT
- Radiotherapy treatment planning
- Radiation oncology