Stereotactic radiosurgery (SRS) plans created using synthetic computed tomography (CT) images derived from magnetic resonance imaging (MRI) data may offer the advantage of inhomogeneity correction by convolution algorithms, as is done for CT-based plans. We sought to determine and validate the clinical significance and accuracy of synthetic CT images for inhomogeneity correction in MRI-only stereotactic radiosurgery plans for treatment of brain tumors.
In this retrospective study, data from two patients with brain metastases and one with meningioma who underwent imaging with multiple modalities and received frameless SRS treatment were analyzed. The SRS plans were generated using a convolution algorithm to account for brain inhomogeneity using CT and synthetic CT images and compared with the original clinical TMR10 plans created using MRI images.
Synthetic CT-derived SRS plans are comparable with CT-based plans using convolution algorithm, and for some targets, based on location, they provided better coverage and a lower maximum dose.
The results suggest similar dose delivery results for CT and synthetic CT-based treatment plans. Synthetic CT plans offered a noticeable improvement in target dose coverage and a more gradual dose fall-off relative to TMR10 MRI-based plans. The major disadvantage is a slightly increased dose (by 0.37%) to nearby healthy tissue (brainstem) for synthetic CT-based plans relative to those created using clinical MRI images, which may be a problem for patients undergoing high-dose treatment.
Magnetic resonance imaging (MRI) is rapidly becoming one of the most important imaging modalities used in stereotactic radiosurgery (SRS) treatment planning for intracranial tumors. SRS is used to correct functional abnormalities and address small tumors of the brain. It can deliver precisely targeted radiation in fewer high-dose treatments than traditional therapy, an approach which can help preserve healthy tissue . During the planning process, MRI and computed tomography (CT) images, and recently, cone beam computed tomography (CBCT) images, are acquired for delineation of tumor and organs at risk (OAR).
The convolution algorithm has recently been introduced as an optimized radiation planning tool in the Gamma Knife planning system. Its use allows correction for brain tissue inhomogeneity and provides more accurate dose calculation. CT images can be used with these algorithms; they provide electron density maps for convolution-based dose calculation and can be co-registered with MR images as a primary dataset to correct for geometric distortion . The use of CT images alone for planning is not ideal due to their lack of quality soft tissue contrast; however, the process of registration with MRI images is associated with some significant uncertainties. CBCT reference images can be used for geometrical distortion correction and pre-treatment positioning on a Gamma Knife Icon machine (Leksell Gamma Knife Icon, Stockholm, Sweden) [2, 3].
Precise target definition is essential for successful SRS treatment. MRI provides superior soft-tissue contrast, physiological information, and multiplanar imaging capabilities, making it an ideal modality for target localization. Therefore, in general, the contours from MRI images are transferred to CT or CBCT images for treatment planning after fusion [4, 5]. However, multi-modality imaging is time-consuming for patients and can introduce systematic positioning errors of up to 2 mm during scanning and fusion, both reducing treatment accuracy and increasing patients’ exposure to imaging radiation .
The use of MRI images alone for radiation treatment planning (RTP), and especially SRS planning, has some limitations, including a lack of electron density information for accurate dose calculation, geometrical distortion of MRI images which affect the accuracy of target definition, and finally, patient localization as a result of lack of bony structure verification during patient setup on the machine . Thus, CT is still the standard imaging modality of choice for treatment planning, and in SRS planning, if convolution is used for correction of tissue inhomogeneity rather than a tissue maximum ratio (TMR) algorithm [7-9].
Many groups are currently working on generating synthetic CT images derived from MRI images to address the first two limitations of RTP based on MRI images alone. This new process will potentially reduce workload, cost, geometrical mismatch associated with the co-registration process, and finally, the patient’s radiation exposure [7, 10-12].
There are different approaches to generate synthetic CT images. One of the simplest and most straightforward is a segmentation-based approach, which involves manual segmentation of each region-of-interest (ROI) over bony anatomy, soft tissue (fat and water), and air on MR images, with assignment of a bulk density to each region [13-15]. The accuracy of this technique has been evaluated for simple 3D conformal, intensity-modulated radiation therapy (IMRT) and SRS plans and has shown good dosimetric agreement for simple 3D conformal plans. However, dose deviations of up to 5% have been reported for IMRT head and neck cancer plans, as well as unacceptable values for SRS brain metastasis plans [16, 17].
The second approach uses an MRI atlas created from paired, co-registered MR-CT datasets to generate a conjugate electron density atlas. One common drawback in this technique is errors arising from mis-registration of atlas data with patient data; this issue is especially common in patients with surgical implants and/or other anatomical abnormalities, which are very common in post-operation radiotherapy patients [18, 19].
Finally, the third approach is the classification-based approach used in this study. This technique aims to classify or cluster different tissue types (i.e., muscle, adipose tissue, and bone) based on pixel intensities from multiple MR pulse sequences. One common problem with this technique is that in order to allow absolute comparison between two patients undergoing the same treatment, the same scanner and parameters (i.e. magnet strength, radiofrequency (RF) coil setup, and MRI pulse sequence) must be used because of the dependency of voxel signal intensity on these factors. Nonetheless, promising results using this technique have been reported for head, neck, and brain tumors [20, 21].
In current practice, during SRS planning using MRI images as a primary dataset, any tissue heterogeneity within the skull is ignored when using the TMR10 algorithm, even though the convolution algorithm is available for such correction. One of the drawbacks of using this algorithm is its speed and potential delay in patient treatment time. While the change in dose distribution when comparing the convolution and TMR10 algorithms is known to an extent, the specific changes that may lead to clinically significant differences are not well understood [22, 23].
Therefore, in this paper, we propose a method for MRI-only SRS planning using synthetic CT images using a classification-based approach, derived directly from MRI images. We investigated the accuracy and dosimetric differences between synthetic CT and CT-based plans created using a convolution-based algorithm, chosen to account for brain tissue inhomogeneity. The SRS treatment plans created in this way were compared with clinical MRI-based plans created using the TMR10 algorithm.
Materials & Methods
The study was approved by the University of Mississippi Medical Center institutional review board. In this retrospective study, we analyzed data from four patients (two with brain metastases and one with meningioma), who had undergone MRI, CT and CBCT scans and received frameless SRS treatment. As part of our clinical workflow for each patient we planned on MRI images co-registered with CT, and in some cases CBCT, for geometric distortion correction. Each plan may include multiple target brain metastases, and some patients had multiple plans with different isocenters. We used all plans in this analysis: eight plans from four patients, for a total of 11 brain metastasis, brainstem, and meningioma target tumors. The specific lesion sites varied from patient to patient.
Synthetic CT images were generated from MR images using syngo.via RT Image Suite (Siemens Healthineers, Erlangen, Germany) using a fuzzy c-means method. Note that the patient setup in MRI was different from that for CT and treatment. Experienced physicians rigidly registered the CT and MRI image datasets and defined the target volume on MRI images, and a medical physicist created SRS treatment plans using the TMR10 algorithm on MRI images. SRS doses were prescribed in the range of 13 to 20 Gray (Gy) to 50% to 75% isodose lines according to the lesion sizes. SRS plans were prepared using Leksell GammaPlan Elekta. In the next step, the volumes were transferred to synthetic CT and CT images. This information was used for planning on both synthetic CT and CT datasets using the convolution algorithm, with the same dose prescription as on primary MRI images. Dose maps were calculated on the synthetic CT and on the original CT, using the same plan parameters.
CT images were acquired using a Siemens SOMATOM Definition AS CT scanner (Siemens Healthineers, Erlangen, Germany) at 120 kVp and 220 mA with the patient in treatment position; each 0.6-mm-thick slice had a resolution of 512×512 pixels, with pixel spacing of 0.586 mm.
MR images were acquired using a Siemens Aera 1.5T diagnostic magnet (Siemens Healthineers, Erlangen, Germany). High-resolution magnitude and phase images were acquired for field mapping after automated shimming over the entire head volume, distortion correction filter on for system inherited gradient non-linearity correction and (Time of Echo 1 (TE1) /Time of Echo 2 (TE2)/Time of Repetition (TR) = 2.46/7.38/12 ms, 550 Hz/pixels, 1 mm3 isotropic sagittal 3D acquisition, standard 20-channel head coil).
Phase images were subtracted and unwrapped to produce field maps (in-house software, IDL 8.2, Boulder, CO, USA and Mathworks, Natick, MA, USA). We derived the final displacement map to correct for residual geometrical distortion using a machine-related distortion map created using a QUASAR MRID3D phantom (Modus Medical Devices Inc., London, Ontario, Canada) (Figure 1) and a B0-derived map. The final displacement map was applied to all MRI images after rescaling using the open-source Analysis of Functional Neuroimaging (AFNI) software package for brain imaging.
All patients were scanned, post-contrast, with a three-dimensional, T1-weighted, gradient-echo sequence (T1-MPRAGE), 1 mm3 isotropic axial 3D acquisition, bandwidth = 160 Hz/pixel, TR= 2200 ms, TE = 3.74 ms, T1 Vibe Dixon, voxel size: 1.0 x 1.3 x 1.3 mm (to generate images of fat and water), TR = 4.10 ms, TE1/TE2 = 1.23 ms/2.46 ms, Flip Angel (FA) = 9 degrees, bandwidth = 1090 Hz/pixel. T2 SPACE, voxel size 1 mm isotropic (anatomical scan appropriate for contouring), TR = 3200 ms, TE = 410 ms, bandwidth = 751 Hz/pixel PETRA, isotropic 1 mm, TR =3.32 ms, TE = 0.07 ms, FA = 6 degree, sagittal, bandwidth = 403 Hz/pixel (identification of air for the purpose of defining an air mask to exclude such voxels from calcification), voxel size 1 mm isotropic, vascular imaging (T2-weighted gradient echo, appropriate to create a threshold-based intensity mask to separate flowing blood), voxel size 3 mm, matrix: 256 x 230, TR = 8.6 ms, TE = 4 ms, FA = 20 degree, bandwidth = 320 Hz/pixel.
To study the accuracy and differences of synthetic CT plans created using the convolution algorithm, we initially visually checked the similarity of dose distribution calculated from both CT and synthetic CT images and compared it qualitatively with dose distribution for TMR10-based MRI plans. Then, by quantitively examining clinically relevant metrics, such as D100 (Gy), D95 (Gy), maximum, minimum, and mean dose were extracted for comparison of dose to target (brain metastasis or meningioma tumor) and OAR (brainstem).
We demonstrate the 2D dose distribution for one of the brain metastases on three planes as an example (Figure 2). Dose distributions appeared similar for both CT and synthetic CT plans and were more conformal for both compared with TMR10-derived plans. Overall, the quantitative data (Table 1) show a 1.84% decrease in maximum point dose inside the targets for synthetic CT and a 1.77% decrease in CT plans compared with TMR10-MRI. The average D100 and D95 for synthetic CT showed a 0.11% increase and 0.19% decrease compared with a 0.67% increase and 0.71% increase in CT-based compared to MRI plans. Overall, synthetic CT plans provided better coverage and a lower maximum dose.
For MRI-only radiotherapy, especially SRS planning, the use of synthetic CT images could potentially involve the application of automatic segmentation, the use of functional MR images, and even of computational texture analysis techniques such as radiomics. Very soon, larger MRI vendors, e.g. GE and Siemens, will offer FDA-approved MRI-only treatment planning platforms . Generation of synthetic CT images and patient-specific distortion correction will be more accessible and streamlined, compatible with busy clinical workflow.
Thus, it is necessary to validate the accuracy of synthetic CT images and to ensure these images’ geometric and dosimetric accuracy is comparable with that of current CT images. Of necessity, this will involve the use of these images for RTP at different sites and both qualitative and quantitative evaluation of the accuracy of dose distribution. This procedure will be more accessible as more MRI-Linac machines become available. In this project, we looked at the application of synthetic CT images during SRS planning, a procedure in which MRI images are frequently used as a primary planning dataset.
Current SRS planning procedures involve the use of MR images without heterogeneity correction. Dose calculation on standard CT scans requires separate registration to MR for delineation of the target and normal tissue . Synthetic CT images generated directly from MR have the potential to take advantage of the improved delineation clarity of MR while maintaining heterogeneity-corrected dose calculation within the treatment planning system. The use of synthetic CT has immediate applications for SRS planning, as well as future implications for developing real-time MR image guidance and plan adaptation for external beam radiotherapy. The incorporation of a convolution algorithm in newer SRS planning systems will improve dose calculation even in the face of tissue inhomogeneity, and especially in regions with air-tissue interfaces and sharp tissue changes.
The results of this study show similarity between CT and synthetic CT-based SRS plans. Synthetic CT offered a noticeable improvement in target dose coverage and a more gradual dose fall off. The disadvantage is an increased dose (by 0.37%) to brainstem, an OAR, relative to plans created using clinical MRI images, which could be significant at higher doses. Future directions for this work will involve 1) further modification of the synthetic CT algorithm to improve its performance; 2) generation of a custom electron density curve to integrate synthetic CT images into the treatment planning system; and 3) study of the use of synthetic CT as the primary calculation image set for different types of treatment plan.
- Phillips MH, Stelzer KJ, Griffin TW, Mayberg MR, Winn HR: Stereotactic radiosurgery: a review and comparison of methods. J Clin Oncol. 1994, 12:1085-1099. 10.1200/JCO.1922.214.171.1245
- Fallows P, Wright G, Harrold N, Bownes P: A comparison of the convolution and TMR10 treatment planning algorithms for Gamma Knife® radiosurgery. J Radiosurgery SBRT. 2018, 5:157-167. Accessed: March 10, 2019: https://www.oldcitypublishing.com/journals/jrsbrt-home/jrsbrt-issue-contents/jrsbrt-volume-5-number-2-2018/jrsbrt-5-2....
- Zhang Q, Chan MF, Burman C, Song Y, Zhang M: Three independent one-dimensional margins for single-fraction frameless stereotactic radiosurgery brain cases using CBCT. Med Phys. 2013, 40:121715. Accessed: March 10, 2019: 10.1118/1.4829517
- Tzikas A, Karaiskos P, Papanikolaou N, et al.: Investigating the clinical aspects of using CT vs. CT-MRI images during organ delineation and treatment planning in prostate cancer radiotherapy. Technol Cancer Res Treat. 2011, 10:231-242. 10.7785/tcrt.2012.500198
- Tan J, Lim Joon D, Fitt G, et al.: The utility of multimodality imaging with CT and MRI in defining rectal tumour volumes for radiotherapy treatment planning: a pilot study. J Med Imaging Radiat Oncol. 2010, 54:562-568. 10.1111/j.1754-9485.2010.02212.x
- Andreasen D, Leemput KV, Hansen RH, Andersen JAL, Edmund JM: Patch-based generation of a pseudo CT from conventional MRI sequences for MRI-only radiotherapy of the brain. Med Phys. 2015, 42:1596-1605. 10.1118/1.4914158
- 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.065
- Chen L, Nguyen T-B, Jones E, et al.: Magnetic resonance-based treatment planning for prostate intensity-modulated radiotherapy: creation of digitally reconstructed radiographs. Int J Radiat Oncol Biol Phys. 2007, 68:903-911. 10.1016/j.ijrobp.2007.02.033
- 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-54
- Baldwin LN, Wachowicz K, Thomas SD, Rivest R, Fallone BG: Characterization, prediction, and correction of geometric distortion in 3 T MR images. Med Phys. 2007, 34:388-399. 10.1118/1.2402331
- Baldwin LN, Wachowicz K, Fallone BG: A two-step scheme for distortion rectification of magnetic resonance images. Med Phys. 2009, 36:3917-3926. 10.1118/1.3180107
- Eilertsen K, Vestad L, Geier O, Skretting A: A simulation of MRI based dose calculations on the basis of radiotherapy planning CT images. Acta Oncol. 2008, 47:1294-1302. 10.1080/02841860802256426
- 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 Lond Engl. 2010, 5:62. 10.1186/1748-717X-5-62
- Lambert J, Greer PB, Menk F, 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.012
- Kapanen M, Collan J, Beule A, Seppälä T, Saarilahti K, Tenhunen M: Commissioning of MRI-only based treatment planning procedure for external beam radiotherapy of prostate. Magn Reson Med. 2013, 70:127-135. 10.1002/mrm.24459
- 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:97-104. 10.1120/jacmp.v12i4.3522
- Chin AL, Lin A, Anamalayil S, Teo BK: Feasibility and limitations of bulk density assignment in MRI for head and neck IMRT treatment planning. J Appl Clin Med Phys. 2014, 15:100-111. 10.1120/jacmp.v15i5.4851
- Dowling JA, Lambert J, Parker J, 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:5-11. 10.1016/j.ijrobp.2011.11.056
- Stanescu T, Jans H, Pervez N, Stavrev P, Fallone B: A study on the magnetic resonance imaging (MRI)-based radiation treatment planning of intracranial lesions. Phys Med Biol. 2008, 53:3579-3593. 10.1088/0031-9155/53/13/013
- Rank CM, Tremmel C, Hünemohr N, Nagel AM, Jäkel O, Greilich S: MRI-based treatment plan simulation and adaptation for ion radiotherapy using a classification-based approach. Radiat Oncol Lond Engl. 2013, 8:51. 10.1186/1748-717X-8-51
- Hsu S, Cao Y, Huang K, Feng M, Balter J: 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-8435. 10.1088/0031-9155/58/23/8419
- Zhou J, Cao Y, Balter J, Chaudhary N, Pandey AS, Thompson G, Tsien CI: Imaging techniques in stereotactic radiosurgery. Principles and Practice of Stereotactic Radiosurgery. Chin L, Regine W (ed): Springer, New York, NY; 2015. 11-23. 10.1007/978-1-4614-8363-2_2
- Alexander E, Kooy HM, van Herk M, et al.: Magnetic resonance image-directed stereotactic neurosurgery: use of image fusion with computerized tomography to enhance spatial accuracy. J Neurosurg. 1995, 83:271-276. 10.3171/jns.1995.83.2.0271
- Kim MM, Parmar HA, Aryal MP, Mayo CS, Balter JM, Lawrence TS, Cao Y: Developing a pipeline for multiparametric MRI-guided radiation therapy: initial results from a phase II clinical trial in newly diagnosed glioblastoma. Tomography. 2019, 5:118-126. 10.18383/j.tom.2018.00035
Evaluation of the Geometric and Dosimetric Accuracy of Synthetic Computed Tomography Images for Magnetic Resonance Imaging-only Stereotactic Radiosurgery
Ethics Statement and Conflict of Interest Disclosures
Human subjects: Consent was obtained by all participants in this study. University of Mississippi Medical Center Institutional Review Board issued approval #2017-0132. The study was conducted under UMMC IRB approval. Animal subjects: All authors have confirmed that this study did not involve animal subjects or tissue. Conflicts of interest: In compliance with the ICMJE uniform disclosure form, all authors declare the following: Payment/services info: All authors have declared that no financial support was received from any organization for the submitted work. Financial relationships: All authors have declared that they have no financial relationships at present or within the previous three years with any organizations that might have an interest in the submitted work. Other relationships: All authors have declared that there are no other relationships or activities that could appear to have influenced the submitted work.
We would like to thank Siemens Medical for providing us with support, and especially the Syngo.via frontier synthetic CT toolbox used to generate images.
Cite this article as:
Fatemi A, Kanakamedala M R, Yang C, et al. (April 08, 2019) Evaluation of the Geometric and Dosimetric Accuracy of Synthetic Computed Tomography Images for Magnetic Resonance Imaging-only Stereotactic Radiosurgery. Cureus 11(4): e4404. doi:10.7759/cureus.4404
Received by Cureus: March 13, 2019
Peer review began: March 27, 2019
Peer review concluded: April 04, 2019
Published: April 08, 2019
© Copyright 2019
Fatemi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License CC-BY 3.0., which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.