Skip to main content

A review of substitute CT generation for MRI-only radiation therapy


Radiotherapy based on magnetic resonance imaging as the sole modality (MRI-only RT) is an area of growing scientific interest due to the increasing use of MRI for both target and normal tissue delineation and the development of MR based delivery systems. One major issue in MRI-only RT is the assignment of electron densities (ED) to MRI scans for dose calculation and a similar need for attenuation correction can be found for hybrid PET/MR systems. The ED assigned MRI scan is here named a substitute CT (sCT). In this review, we report on a collection of typical performance values for a number of main approaches encountered in the literature for sCT generation as compared to CT. A literature search in the Scopus database resulted in 254 papers which were included in this investigation. A final number of 50 contributions which fulfilled all inclusion criteria were categorized according to applied method, MRI sequence/contrast involved, number of subjects included and anatomical site investigated. The latter included brain, torso, prostate and phantoms. The contributions geometric and/or dosimetric performance metrics were also noted. The majority of studies are carried out on the brain for 5–10 patients with PET/MR applications in mind using a voxel based method. T1 weighted images are most commonly applied. The overall dosimetric agreement is in the order of 0.3–2.5%. A strict gamma criterion of 1% and 1mm has a range of passing rates from 68 to 94% while less strict criteria show pass rates > 98%. The mean absolute error (MAE) is between 80 and 200 HU for the brain and around 40 HU for the prostate. The Dice score for bone is between 0.5 and 0.95. The specificity and sensitivity is reported in the upper 80s% for both quantities and correctly classified voxels average around 84%. The review shows that a variety of promising approaches exist that seem clinical acceptable even with standard clinical MRI sequences. A consistent reference frame for method benchmarking is probably necessary to move the field further towards a widespread clinical implementation.


Dose calculations performed on scans from magnetic resonance imaging (MRI) were first reported around the millennium when MRI emerged as a complimentary modality to computed tomography (CT) in the delineation step of the radiotherapy (RT) chain [1, 2]. As MRI provides superior soft tissue contrast and delineation precision as compared to CT [38], the concept of carrying out all steps of the RT chain on MRI as the sole modality, so-called MRI-only RT, could provide a favorable workflow. MRI-only RT would further remove a systematic registration error when transferring MRI delineated structures to the CT which has been reported to be in the order of 2–5 mm for various treatment sites [913]. As CT is used for positioning of the patient at treatment, registration errors introduce a spatial systematic uncertainty. The dosimetric impact of a systematic error will increase when the radiation is aimed at small structures or when the target is close to sensitive organs. This could be the case for small tumors or the hippocampus in the brain [14] with a structure radius in the order of a possible registration error or when a standard PTV or PRV margin has to be compromised to maintain an acceptable therapeutic ratio. An example of the MRI to CT registration variability for a prostate and nasopharynx case is illustrated in Fig. 1.

Fig. 1

Variability of multiple registrations between MRI and the corresponding CT for prostate (top) and nasopharynx (bottom). a: One marker (of three) indicated by a white circle on the axial MRI. b: Two markers shown by the white dots on the sagittal CT. The multiple thin white lines are the MRI delineated clinical target volume (CTV) transferred to the CT based on the marker registration following department protocol from 7 different observers. The protocol is based on a rigid automatic (mutual information) registration for a limited FOV around the prostate followed by a manual adjustment to match the markers in the three planes. The outermost white line was the planned target volume (PTV) applied. The data are taken from reference [59]. c: The gross target volume (inner) and CTV (outer) as delineated on the axial MRI. d: The multiple thin white lines are the MRI based CTV transferred to the CT by 6 different observers (sagital CT slice shown). The registration is based on a rigid automatic (mutual information) registration followed by a manual fine adjustment. The outermost white line was again the applied PTV. The data are taken from reference [60]

In addition, MRI-only will decrease the number of scans and associated patient discomfort, and, reduce the planning related costs [15]. The benefits of MRI-only RT would further increase in a workflow with repeated imaging, e.g. weekly scans, for response assessment and/or treatment adaptation.

A number of concerns related to MRI-only RT exist. One major challenge of performing dose calculations on MRI is the lack of correspondence between the voxel intensity and the associated attenuation property of the tissue. Unlike CT images where the voxel intensity directly reflects the radiological characteristics of the tissue, MRI intensities rather correlate with tissue proton density and the magnetic relaxation, i.e. the inertia of the dipole moment [16]. This leads to voxel ambiguity for tissues such as bone and air which both appear dark on the MRI although they have very different attenuation coefficients. The focus of this review is strategies for dealing with this ambiguity. Further challenges constitute scanner induced geometrical distortions arising from gradient non-linearity and magnet inhomogeneities and patient induced artifacts such as susceptibility and chemical shifts [13]. Specific problems for algorithms converting the MRI signal into a CT number further constitute normalization of absolute signal intensities and data correction strategies such as bias field correction. These topics are considered out of scope of this review.

The increased use of MRI for target and normal tissue delineation in RT in general and two device-driven events have facilitated scientific activity for assigning electron densities to MRI imagesFootnote 1. The first event is the commercial availability of clinical integrated hybrid PET/MRI systems around 2010-2011 [17]. Unlike traditional PET/CT systems where the CT scan is used for attenuation correction of the PET signal needed in quantitative PET volume estimates such as the standard uptake volume (SUV), attenuation coefficients need to be assigned to the MRI scan in hybrid PET/MRI systems to make a similar attenuation correction. The second event is the commercial availability of integrated MRI guided systems in external beam RT around 2014 [18]. These systems can provide MRI scans for patient setup based on soft tissues and monitor the tumor movement during treatment delivery. The systems would be able to calculate and adapt the dose distribution at each given fraction if a dose calculation can be performed on the MRI scan which requires electron densities to be assigned. The increased focus on MRI in RT and the introduction of the imaging and treatment devices is reflected in the number of publications as illustrated in Fig. 2.

Fig. 2

The number of articles versus time after applying the search string and exclusion criterion 1 provided in the text. The publications are grouped into proposed applications of the described method and sorted according to publication year. A boost in the amount of publications can be inspected around 2010 (PET/MRI) and 2013-14 (MRI guided EBRT)

In this review, we will report on a collection of typical performance values for a number of main approaches encountered in the literature for conversion of MR data to electron density or HU maps relevant for RT. The generated map is called “synthetic CT”, “substitute CT”,”pseudo CT” or similar, i.e. no common terminology is currently established. In the following, we will use the term substitute CT (sCT) since the acronym for pseudo CT (pCT) is often used for the planning CT in adaptive RT studies [19, 20]. This field of research covers both diagnostic and therapeutic radiology as well as MRI. Further, the research field expands into automatic organ segmentation and image analysis in general. As a result, a diverse amount of approaches with different scientific traditions, terminology and endpoints in mind have been reported in the literature. Consequently, we have had to make simplifications and compromise details in order to preserve an overview and to further categorize and score the methods with the aim of providing results relevant for RT. Therefore, a direct comparison between the reported performance metrics presented here is not valid. Rather, the idea is to provide an overview of the multiple strategies investigated in the literature along with a general order of accuracy that can be reached.

We have categorized the investigated methods into three main approaches. These are termed voxel [2123], atlas [24, 25] and hybrid [26, 27]. The latter use a combination of the voxel and atlas based approach. Similar approaches are typically categories as segmentation, sequence/image contrast, template and atlas in PET/MRI terminology [16, 28]. The voxel based approach primarily uses information about voxel intensities (contrasts) in the MR images to assign electron densities. No or limited information about the location of voxels is included in this category. The voxel based methods are dominated by the concept of machine learning in which part of the data is used to train (optimize) a model which is then applied on the remaining MRI data to predict the CT numbers. If a study contains n patients, this is usually done by training the model on n-1 patients and then predict the sCT on the remaining patient in a rotating scheme known as leave-one-out cross validation. The electron density assignment (CT number) can be made on the basis of generic values (e.g. from ICRU report 46 [29]) to bulk groups of voxels. Alternatively, the assignment can happen on a continuous scale by including patient specific CT numbers in a training phase. In contrast, the atlas based approaches focus on aligning the location of a patients MRI voxels to the corresponding location of a MRI voxels in an atlas through registration. The atlas can either be a single or average (template) patient or contain a number of patients (often termed multi-atlas). The atlas contains a pre-known correlation between the MRI voxels and the value of interest, e.g. CT number or organ label. Once the alignment has taken place, the atlas CT number can be assigned to the patients’ MRI scan and hence converting it into a sCT scan.

A large number of MRI sequences/contrasts for electron density assignment have been reported in the literature. We have chosen to divide the MR input images into four main contrasts / sequences categories. The first two categories are simply termed T1 weighted (T1w) and T2 weighted (T2w). They are based on common clinical MRI sequences which rely on either the longitudinal (T1) or the transverse (T2 or T2*) tissue relaxation to produce image contrast. The T1 and T2 relaxation is determined from multiple refocusing pluses during the repetition time while T2* describes the relaxation of the free induction decay (FID) produced in the receiver antenna coil. Two main pulse sequences exist for MR image acquisition: spin echo (SE) and gradient echo (GE). The SE MR signal intensity is roughly proportional to ρ[1-exp(TR/T1)]exp(-TE/T2) where ρ (proton density), T1 and T2 are tissue properties and TR (repetition time) and TE (echo time) are sequence parameters. The equation is only valid if TR > > TE which is usually the case, and, in general T1 > T2 > T2* relaxation [30]. T1w images (short TR, short TE) are preferred for visualizing anatomy while T2w images (long TR, long TE) are usually the choice for visualizing pathology. The third category comprises the Dixon family of fat-water separating sequences and is collectively termed Dixon [31]. It is based on the chemical shift between the resonance frequencies of fat and water and can be weighted towards T1, T2 or ρ as it is (typically) a SE sequence. The fourth category of MRI sequences is based on dual ultrashort echo time (dUTE) to visualize solid structures with a very short T2 relaxation time such as the bone [32, 33]. In dUTE image acquisition, a first signal is collected right after the excitation, and a second using the GE technique at a longer nominal echo-time. The first image is ρ weighted or T1w depending on the flip-angle and the second will have a T2*w or T1w contrast depending on the echo-time and the flip-angle. T2* is only possible to realize with gradient echo sequences.

Material and methods

Literature search

This review reports on a collection of typical performance values for substitute CT generation rather than giving a detailed theoretical background of the methods used to predict substitute CT. To give a fair representation on the scientific activity within this field, we performed a literature search in the Scopus database November 2015 [34]. Index terms such as Medical Subject Headings (MeSH) terms were not used to define a search due to the wide diversity of sciences involved in the research field and the consequent lack of common terminology. Instead, a collection of common keywords found in a number of MRI-only RT and PET/MRI articles were organized in a logical search string defining the inclusion criteria:

  • TITLE-ABS-KEY [(“PET MRI” OR “MR PET”) AND NOT (functional OR diffusion OR fdg-spect)]


  • TITLE-ABS-KEY [(radiotherapy OR “radiation therapy”) AND (“magnetic resonance imaging” OR “magnetic resonance” OR mri OR mr) AND NOT chemotherapy]


  • TITLE-ABS-KEY [“Attenuation correction” OR “computed tomography substitute” OR “substitute CT” OR “pseudo CT” OR “MRI only” OR “MRI alone”]

where TITLE-ABS-KEY indicate either the title, abstract or keywords of the paper. This resulted in 254 papers and we further added 7 papers/abstracts that for various reasons were not found in the structured literature search (e.g. strange keywords, conference abstracts etc.). Three exclusion steps were then introduced. Exclusion 1 was defined as papers having TITLE-ABS-KEY on the following:

  • Diagnostics and delineations based on MRI only.

  • Brachytherapy

  • CT to MRI registration error/IGRT studies

  • Subjects specific for PET/MRI: field-of-view (FOV) truncation, effects of headphone and coils, etc.

  • PET/MRI specific corrections: time-of-flight, line source, maximum likelihood for attenuation and activity (MLAA).

The application of exclusion 1 reduced the number of papers to 117. These included relevant investigations for assigning electron densities to MRI scans for the purpose of PET attenuation correction, MRI-only RT or both (see Fig. 2). Exclusion 2 intended to only include studies which presented novel methods for electron density assignment. Abstracts and manuscripts including the following were excluded:

  • Review articles

  • Book series or only insufficient abstracts available.

  • RT feasibility/comparative studies: Dose calculations incl. and excl. CT transferred structures such as bone and air cavities to the MRI. CT bulky assigned density vs. normal CT etc.

  • PET/MRI feasibility/comparative studies: Difference in SUV or similar by applying CT based vs. MRI corrected attenuation maps incl. or excl. CT transferred structures.

  • Focus on MRI artifacts and distortion quantifications and corrections.

  • MRI-only based workflow descriptions.

After exclusion 2, 73 papers presenting novel correction methods remained (see Fig. 3). Whenever multiple method approaches and /or MRI sequences were used these were collectively categorized as “hybrid” and “multiple”, respectively.

Fig. 3

Categorization of contributions after applying exclusion criterion 2. The papers were sorted according to their main method (left), used MRI sequences/contrasts (middle) and number of subjects, i.e. phantom, patients or volunteers (right). In the above categories, the following simplifications were made: head = brain, whole body = torso, cervix = prostate (only 1 study), UTE = dUTE, and water/fat separating MRI sequences = Dixon. Volunteers and phantoms were categorized as patients. Some papers included description of multiple methods which were included in the histograms as separate studies, hence the term “published studies” for the ordinate

Exclusion 3 intended to only consider studies which included a quantitative performance metric of the resulting sCT scan which would be relevant or applicable for RT purposes. The following papers were excluded:

  • No reported quantitative performance metric.

  • Reported performance metric not relevant to MRI-only RT, e.g. differences in SUV, linear correlation coefficients of activity estimates, etc.

The final 50 papers are shown in Table 1 and the selection process is summarized in Fig. 4. The final papers were arranged in the main categories as described in the introduction and further subcategorized within each main approach and applied MRI sequences when possible. The latter was limited to the four overall sequence categories.

Table 1 Table summarizing performance over different substitute CT approaches
Fig. 4

The adopted strategy for inclusion of papers with reported metrics in this review

Method categories

The methods were organized with subcategories according to a main voxel, atlas or hybrid approach in Table 1. Voxel based sCT generation utilizes the contrast in the MR image independently of the voxels spatial location. This makes it a potential computational attractive approach. The voxels can be converted into CT numbers in multiple ways which we have subcategorized below.

Semi-automatic refers to some kind of manual intervention from the user to make the method work, e.g. delineation of the bone or manually established intensity thresholds below/above which a voxel is categorized into a certain tissue category. Threshold covers methods which use the tissue relaxation constant to differentiate between different tissues. For the dUTE sequence/contrast, the T2 relaxation time of the bone (0.5–2 ms depending on magnetic field strength) can as an example be used to categorize voxel which decays > 1/(0.5–2 ms) into bone and voxels which remain constant (or decay slowly) over the short acquisition time into soft tissue [35]. Probabilistic refers to methods which can assign a probability of a voxel to belong to different tissue classes with e.g. a corresponding bulk electron density (sCT) or organ label (auto contouring). This could be done by assuming that the MR intensities come from a mixture of K normal distributions (tissue classes) with a corresponding mean and standard deviation. The initial mixture can be estimated with an expectation maximization algorithm through unsupervised training, i.e. an electron density (CT number) is assigned to each tissue class subsequently. For each voxel, a probability can then be calculated for all tissue classes and the voxel can be assigned to the tissue class for which the highest probability was calculated. This is known as Bayesian statistics. Markow Random Fields include the tissue classes of neighbor voxels in the probability calculation of a given voxel [22]. Fuzzy c-means clustering use a similar strategy of dividing the MRI voxels into distinct clusters (tissue classes). Cluster similarity coefficients are then calculated for each voxel which is then assigned to the tissue which it resembles the most. A number of different similarity measures exist. Regression collects methods which correlate MRI intensities to CT numbers through (statistical) regression on a continuous scale. This can be done in a fashion similar to the Bayesian approach here termed Gaussian by including co-registered CT intensities in an initial supervised training phase of the data that establish the prior mixture of the K tissue classes [23]. Discriminant analysis, Principal Component Analysis and Random Forest are other strategies of performing such a regression between MR and corresponding CT data. Sinogram use a forward projection CT like approach to transform the MRI scan into raw MR data where the different tissues are subsequently identified. Neural network describes supervised training of a correlation model in a hidden layer with an MRI input layer and a CT output layer. Pattern recognition in a voxel setting compares an MRI pattern, e.g. a cluster of 3x3x3 size MR voxels known as a patch, with a pre-established correlation between MR patterns and CT numbers obtained through supervised training. Different measures for pattern similarity can be used such as the (normalized intensity) Euclidean distance. Hybrid voxel methods combine a voxel based method with somewhat loose information of the voxels location, e.g. distance from center of the brain.

Atlas based sCT generation use the location of an MRI voxel to establish the corresponding CT number by aligning the voxel to an atlas with a pre-known correlation between the MRI voxel location and corresponding CT number. This is potentially more computational challenging as each patient's MRI has to be aligned with an atlas with no possibility of exploiting a pre-training model (except for the atlas building). In an atlas setting, Pattern Recognition compares similarities of patient MRI patches (see above) with atlas MRI patches within a limited search volume after alignment. Deformable refers to methods which use deformable registration, i.e. non-rigid registration, to assign CT numbers from an atlas to a given MRI scan. The patient's MRI is first registered non-linearly with the atlas MRI. This could be one registration if the atlas consists of one average or template patient. Otherwise, the MRI has to be individually registered to all MRIs in a multi-atlas which is computationally less attractive. Alternative, one can settle with one (entrance) registration if the multi-atlas is internally registered, i.e. a registration map between the entrance atlas and the other atlases has been pre-established. The deformation map between the patient and atlas MRI is then applied to the corresponding atlas CT and the sCT produced. If multiple atlases are used, a fused CT number can be applied.

Hybrid atlas methods combine multiple methods within the atlas based category. This could be a method which combines deformable registration with a (patch) pattern recognition approach to minimize the influence of registrations which resulted in a poor alignment.

The final Hybrid approach combines categories of the voxel and atlas based approaches. This could for example be a calculation of two probability density functions (PDFs) for each voxel; one based on deformable registration (atlas) and the other based on Bayesian statistics (voxel). The two probabilities are then combined into a unified posterior PDF which determines the final assignment of CT number to the voxel [26]. An example of atlas and voxel based sCT generation for the pelvis and brain can be seen in Fig. 5.

Fig. 5

Examples of sCT generation for the pelvic (top) and brain (bottom). a: An axial CT slice of the pelvic from a prostate patient. b: The corresponding sCT slice created with an atlas patch based approach [49]. c: An axial CT slice of a brain patient. d: The corresponding sCT slice created with a voxel Gaussian mixture regression based approach [75]

Performance metrics

Three common metrics reported in the literature to score the performance of a given sCT generation method were chosen. The first metric, ΔDose, describes the dosimetric agreement when performing dose calculations on the sCT as compared to the standard CT. This is usually quantified as the percentage difference in either single characteristic points, e.g. iso-center or dose prescription point, or in dose volume histogram (DVH) points. The general equivalent uniform dose (gEUD) has also been used to describe the biologically relevant differences of the entire DVH. Another commonly reported metric for quantifying dosimetric differences is the gamma index [36]. This metric covers spatially correlated dose deviation in both the high and low dose regions. Whenever differences in multiple dose metrics were reported, e.g. multiply DVH points, a collectively representative value, e.g. the mean of all the deviations, was chosen. The second metric is the mean absolute error (MAE) and describes the absolute voxel-wise difference in HU defined as

$$ M A E=\frac{1}{N}{\displaystyle \sum_{i=1}^N}\left| C{T}_i- sC{T}_i\right| $$

where N is the number of voxels, CT is the standard CT and sCT is the substitute CT. This metric is typically lowered as the number of voxel similar to water or air increases, e.g. moving from the brain to the pelvis or including air outside the body outline (see Fig. 6a). The MAE shows a great variation over the different tissue regions. For the data presented in Fig. 6a, the shown median MAE of 87 HU for all tissue in the skull covers a median MAE of 216, 36 and 261 HU for air, soft tissue and bone, respectively.

Fig. 6

Performance metrics dependence on the region of interest and CT threshold number for the bone. All metrics are calculated on substitute CTs from 3D T1w MR images of 6 brain patients [47] using an atlas patch based method [48]. a: The boxplot shows the MAE within the body outline (left) and the whole field of view (FOV, right). The medians were 87 and 49 HU for the body and FOV MAE, respectively, and were significantly different (p < 0.002). b: The Dice similarity coefficient (DSC) metric for bone as a function of threshold CT number. c: Receiver operating characteristic (ROC) curve of the sCT bone as a function of CT threshold number (thres). True positive (TP) = sCT > thres & CT > thres, false positive (FP) = sCT > thres & CT < thres, true negative (TN) = sCT < thres & CT < thres and false negative (FN) = sCT < thres & CT > thres. Sensitivity = TP/(TP + FN) and specificity = TN/(TN + FP). The threshold was varied from 100 (right) to 3000 (left) HU in steps of 100 HU. Only voxels > 100 HU on the CT, i.e. the bone region, was included in the evaluation to keep the TN number (non-bone tissue on the sCT and CT) to a reasonable number

The third metric, the dice similarity coefficient (DSC) for bone [37], is a geometric score describing the overlap between the CT and sCT bone volumes. It is defined as

$$ D S{C}_{bone}=\frac{2\left({V}_{C T}{\displaystyle \cap }{V}_{sCT}\right)}{V_{C T}+{V}_{sCT}} $$

where V is the volume of the bone on the CT and substitute CT (sCT), respectively. A similar measure, the Jaccard coefficient (JAC) can be converted to the DSC through the relation DSC = 2 · JAC/(1 + JAC) [38]. Another commonly reported metrics is the sCT specificity and sensitivity for different tissues such as the bone. The bone DSC, sensitivity and specificity will depend on the threshold value set for the CT number of the bone (see Fig. 6b and c). Other metrics are otherwise noted in the legend of Table 1.

The performance metrics reported in the literature cannot be directly compared due to issues such as patient selection and exclusion criteria which are often underreported and can introduce a bias. Further, the amount of preprocessing included in the algorithm such as data normalization and bias field correction will affect the final result. Still, methods performing equally well in the same body region should produce performance metrics within the same gross interval.


The statistics reported in Fig. 2 and Fig. 3 indicate that most studies are carried out on the brain for 5-10 patients with PET/MR applications in mind using a voxel based method. It is common to use model input data that coincide with the image data used for the delineation of the target volume or organs at risk. The main benefits are avoidance of unnecessary registrations in the workflow to compensate for intra examination patient motion and to keep the examination time as short as possible. Therefore evaluation of T1w input data is common for the brain region, while T2w input data is common for the pelvic region. dUTE has the benefit of enabling separation of cortical bone and air, but has not been reported useful for delineation purposes. Our review shows that at present point dUTE has only been evaluated for intra-cranial conditions or phantoms (Table 1). Dixon based sequences enable separation of water and fat signal, and is currently often used for attenuation correction of PET data in hybrid PET/MR scanners. Dixon sequences also tend to be fast and the in-phase sequence has been reported useful for identification of fiducial makers in the prostate [39]. Table 1 shows a wide diversity in terms of the methods and MRI sequences investigated in the literature. Overall, the dosimetric agreement is in the order of 0.3–2.5%. A strict gamma criterion of 1% and 1mm has a range of passing rates from 68 to 94% while less strict criteria show pass rates > 98%. Given the relatively small order of dosimetric disagreement, the residual distortion (i.e. remaining distortion after applying distortion correction procedures) present in the MR and hence sCT images subject to dose calculations seem to be of minor importance in order to reach an acceptable dosimetric accuracy. Rather, these seem to be more critical for accurate target and OAR delineation [40]. The MAE is between 80 and 200 HU for the brain with a majority of values lying in the 120–140 HU interval. The MAE values are around 40 HU for the prostate (pelvis region). The Dice score for bone is between 0.5 and 0.95 across the different methods, MRI sequences/contrasts and anatomical sites. The specificity and sensitivity range from 75 to 98%. As is apparent from Fig. 6c, increasing the specificity will decrease the sensitivity and vice versa. A compromise seems to be in the upper 80s for both quantities. Correctly classified voxels average around 84% for the different methods. No strict relationship between the dosimetric and geometric agreement as scored by the metrics common in the literature is present. Further, when a large number of patients are present in a study, i.e. more than 20 patients, the dice score seem to be in the lower range (<0.85) of the reported values probably caused by a larger diversity in the patient material.


Table 1 is not a complete list of all methods used for generating a sCT, due to the search strategy, time of search and exclusion criteria applied, some novel methods such as sparse representation [41] and Random Forest with auto-context modelling [42, 43] are not included. With these limitations in mind, our investigation shows that sCT generation for MRI-only based radiotherapy or PET/MRI attenuation correction seems to be a comprehensively tested area of research given the variety of investigated methods which are able to produce a sCT from an MRI scan. This presents an advantage in the sense that a broad material is currently available to further develop on. This strength, however, also presents a challenge for the field. There is no obvious method or MRI contrast(s) which seem to be clearly favorable from the others and hence no clear indication as to which path of promising methods to pursue. The data in Table 1do not indicate that inclusion of more MR contrasts in the generation of the sCT automatically increase the accuracy. This is encouraging as extra sequences both increase the total acquisition time and overall complexity of the method and the workflow. Notably, the dice score of the voxel hybrid approaches using Random Forest (Random F) in Table 1 creates an almost perfect overlap between CT and sCT bone volumes especially considering the high threshold of 600 HU applied in these studies (see Fig. 6b). These reported results should encourage other researchers to reproduce this method on different datasets especially since a pure Random F approach on dUTE images does not produce similar high performance scores. The average dosimetric deviations reported in Table 1 should be used with some caution. Korsholm et al. reported on dose deviations using bulk density corrections for multiple treatment sites [44]. With a bulk density correction for the bone, the average deviation for median PTV dose of the prostate was practically zero but covers a range of -1.1 to 1.1% representing the 95% confidence interval. They further argued that the 95% confidence interval should be within a 2% dosimetric deviation to produce clinical acceptable results.

The focus of the present review is the conversion of image data acquired with MRI to electron density or HU maps to facilitate a so called MR-only treatment planning workflow. In addition and not limited to the MR-only approach, use of MR in radiotherapy requires MR data with a minimum of geometrical distortions. It has been shown that distortions due to non-linearity in the spatial encoding gradients can be successfully corrected using deterministic algorithms [45], and chemical shift artifacts and distortions caused by susceptibility effects can be minimized using a sufficient bandwidth [46]. There is, however, still a need to further confirm these results and develop efficient quality control techniques, but this is out of the scope of the current review.


One could argue for a need to unify the efforts to localize promising candidate for sCT generation. A possible way is to standardize the calculation of the performance metrics reported in the literature. It is clear from Fig. 6 that the MAE, DSCbone, specificity and sensitivity metrics have their limitations and further do not necessarily reflect the corresponding dosimetric performance of the method. Quantitative metrics that more unambiguously reflect a correlation between the geometrical and dosimetrical agreement are needed and these should further display an independence of parameters such as selected CT number threshold and field-of-view. An example could be to score the sCT-CT difference in bins covering the HU scale independent of the number of voxel present in each bin [47]. Another example could be to use differences in radiologic (water equivalent) path lengths which represent both a geometric and dosimetric property [48, 49]. The accuracy of the performance metric should further be related to the application in question, e.g. RT or PET/MRI. To benchmark results, another possible way would be to make datasets consisting of MRI scans with a variety of contrasts and corresponding CT scans for different anatomical sites public available. Such a dataset currently exists for quantitative imaging of biomarkers [50] and a similar dataset for sCT generation could serve as a mandatory step for method benchmarking before publication. Issues related to the algorithms processing of MRI signal normalization and correction would further become apparent on such a dataset. Currently, many studies only benchmark the sCT against the gold standard CT. In addition, a dosimetric comparison between a proposed method and the investigated MRI scan(s) set to water should be included in a study to give a perspective of the reported quantities. Another issue is a number of relevant items which are often not addressed. These cover inclusion or exclusion criteria for the investigated patients which could create a bias, computation time needed for a given method to work and restrictions in the possible clinical implementation of the method.

Clinical implementation

So far, the authors are aware of two institutes which currently have implemented in-house developed MRI-only methods clinically. The first institute use a voxel based dual regression approach on the prostate and have currently treated around 150Footnote 2 patients [21, 51]. The second institute use a voxel based probabilistic approach with fuzzy c-means on the brain and have currently treated around 30 whole brain and 15Footnote 3 focal brain cases, respectively [52]. Further, commercial solutions for sCT generation are becoming available [27, 53]. Stereotactic radiosurgery treatment planning of brain tumors have been carried out on MRI scans set to water for decades [54]. Brain radiosurgery was outside the scope of this review but one could question this practice given published literature which demonstrates an improvement in dosimetric accuracy as compared to a pure water based dose calculation for the brain, see e.g. [55]. For the brain and other treatment sites, it has further been shown that a bulk density assignment is probably sufficient for RT treatment planning [44, 56].

Treatment delivery and quality assurance (QA) of an MRI-only based RT workflow are other important issues which need to be addressed for clinical implementation of MRI-only RT. In terms of commissioning an MRI-only workflow, the guidelines as to what tolerances which are acceptable are limited. As mentioned earlier, a 2% dosimetric agreement between sCT and CT based dose calculation seems to be of an acceptable order [44] and most of the investigated methods demonstrate an agreement better than this. For kV X-ray based image-guided RT (IGRT) delivery systems, the cone beam CT (CBCT) seems to provide an acceptable solution for patient setup of both brain and pelvic patients (marker match is probably sufficient for prostate) [47, 57, 58]. The CBCT can further be used for patient specific QA of the generated sCT as it provides an independent estimate of the CT numbers [47]. sCT QA verification efforts using the electronic portal imaging device (EPID) has also been proposed [53]. All of the above elements should be considered when formulating an MRI-only based RT protocol.


In summary, a variety of promising approaches for substitute CT generation exist which seem to provide results acceptable for clinical implementation. This also includes methods based on clinical simple standard MRI sequences/contrasts. However, the field suffers from a current lack of an established benchmarking method and reporting consistency, which challenge the commitment from interested vendors and hence risk a delay for a broad clinical implementation.


  1. 1.

    In photon-based RT, the main interest is to convert the CT number (HU) into an electron density (relative to water) for dose calculation purposes. Therefore, CT number and electron density is used interchangeable throughout the text.

  2. 2.

    Korhonen J. Personal correspondance. 2015.

  3. 3.

    Balter JM. Personal correspondance. 2015.


  1. 1.

    Beavis AW, Gibbs P, Dealey RA, Whitton VJ. Radiotherapy treatment planning of brain tumours using MRI alone. Br J Radiol. 1998;71:544–8.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Chen LL, Price RA, Wang L, Li JS, Qin LH, McNeeley S, Ma CMC, Freedman GM, Pollack A. MRI-based treatment planning for radiotherapy: Dosimetric verification for prostate IMRT. Int J Radiat Oncol Biol Phys. 2004;60:636–47.

    Article  PubMed  Google Scholar 

  3. 3.

    Rasch C, Steenbakkers R, van Herk M. Target definition in prostate, head, and neck. Semin Radiat Oncol. 2005;15:136–45.

    Article  PubMed  Google Scholar 

  4. 4.

    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:121.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Fiorentino A, Caivano R, Pedicini P, Fusco V. Clinical target volume definition for glioblastoma radiotherapy planning: magnetic resonance imaging and computed tomography. Clin Transl Oncol. 2013;15:754–8.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Barillot I, Reynaud-Bougnoux A. The use of MRI in planning radiotherapy for gynaecological tumours. Cancer Imaging. 2006;6:100–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Thiagarajan A, Caria N, Schoder H, Iyer N, Wolden S, Wong RJ, Kraus DH, 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;78:S428.

    Article  Google Scholar 

  8. 8.

    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–7.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    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–82.

    Article  PubMed  Google Scholar 

  10. 10.

    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–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    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–9.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Daisne JF, Sibomana M, Bol A, Cosnard G, Lonneux M, Gregoire V. Evaluation of a multimodality image (CT, MRI and PET) coregistration procedure on phantom and head and neck cancer patients: accuracy, reproducibility and consistency. Radiother Oncol. 2003;69:237–45.

    Article  PubMed  Google Scholar 

  13. 13.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Soltys SG, Kirkpatrick JP, Laack NN, Kavanagh BD, Breneman JC, Shih HA. Is Less, More? The Evolving Role of Radiation Therapy for Brain Metastases. Int J Radiat Oncol Biol Phys. 2015;92:963–6.

    Article  PubMed  Google Scholar 

  15. 15.

    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.

    Article  PubMed  Google Scholar 

  16. 16.

    Hofmann M, Pichler B, Scholkopf B, Beyer T. Towards quantitative PET/MRI: a review of MR-based attenuation correction techniques. Eur J Nucl Med Mol Imaging. 2009;36:93–104.

    Article  Google Scholar 

  17. 17.

    Delso G, Fürst S, Jakoby B, Ladebeck R, Ganter C, Nekolla SG, Schwaiger M, Ziegler SI. Performance measurements of the Siemens mMR integrated whole-body PET/MR scanner. J Nucl Med. 2011;52:1914–22.

    Article  PubMed  Google Scholar 

  18. 18.

    Lagendijk JJ, Raaymakers BW, Van den Berg CA, Moerland MA, Philippens ME, van Vulpen M. MR guidance in radiotherapy. Phys Med Biol. 2014;59:R349–69.

    Article  PubMed  Google Scholar 

  19. 19.

    Behrens CF, Eiland RB, Sjostrom D, Maare C, Paulsen RR, Samsoe E. Adaptive RT for Head-and-Neck Cancer: The Usefulness of Deformable Image Registration. Int J Radiat Oncol Biol Phys. 2012;84:S775.

    Article  Google Scholar 

  20. 20.

    Nijkamp J, Marijnen C, van Herk M, van Triest B, Sonke J-J. Adaptive radiotherapy for long course neo-adjuvant treatment of rectal cancer. Radiother Oncol. 2012;103:353–9.

    Article  PubMed  Google Scholar 

  21. 21.

    Korhonen J, Kapanen M, Keyrilainen J, Seppala T, Tenhunen M. A dual model HU conversion from MRI intensity values within and outside of bone segment for MRI-based radiotherapy treatment planning of prostate cancer. Med Phys. 2014;41(1):011704.

    Article  PubMed  Google Scholar 

  22. 22.

    Edmund JM, Kjer HM, Van Leemput K, Hansen RH, Andersen JAL, Andreasen D. A voxel-based investigation for MRI-only radiotherapy of the brain using ultra short echo times. Phys Med Biol. 2014;59:7501–19.

    Article  PubMed  Google Scholar 

  23. 23.

    Johansson A, Karlsson M, Nyholm T. CT substitute derived from MRI sequences with ultrashort echo time. Med Phys. 2011;38:2708–14.

    Article  PubMed  Google Scholar 

  24. 24.

    Sjolund J, Forsberg D, Andersson M, Knutsson H. Generating patient specific pseudo-CT of the head from MR using atlas-based regression. Phys Med Biol. 2015;60:825–39.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Dowling JA, Lambert J, Parker J, Salvado O, Fripp J, Capp A, Wratten C, Denham JW, Greer PB. 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.

    Article  PubMed  Google Scholar 

  26. 26.

    Gudur MSR, Hara W, Le QT, Wang L, Xing L, Li RJ. A unifying probabilistic Bayesian approach to derive electron density from MRI for radiation therapy treatment planning. Phys Med Biol. 2014;59:6595–606.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Siversson C, Nordstrom F, Nilsson T, Nyholm T, Jonsson J, Gunnlaugsson A, Olsson LE. Technical Note: MRI only prostate radiotherapy planning using the statistical decomposition algorithm. Med Phys. 2015;42:6090–7.

    Article  PubMed  Google Scholar 

  28. 28.

    Wagenknecht G, Kaiser H-J, Mottaghy FM, Herzog H. MRI for attenuation correction in PET: methods and challenges. MAGMA. 2013;26:99–113.

    Article  PubMed  Google Scholar 

  29. 29.

    International Commission on Radiation Units and Measurements (ICRU). Photon, Electron, Proton and Neutron Interaction Data for Body Tissues. ICRU Report 46. Bethesda: International Commission on Radiation Units and Measurements; 1992.

  30. 30.

    Bushberg JT, Boone JM. The essential physics of medical imaging. 3rd edition. Philadelphia: Lippincott Williams & Wilkins; 2011.

  31. 31.

    Dixon WT. Simple proton spectroscopic imaging. Radiology. 1984;153:189–94.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Robson MD, Gatehouse PD, Bydder M, Bydder GM. Magnetic resonance: An introduction to ultrashort TE (UTE) imaging. J Comput Assist Tomogr. 2003;27:825–46.

    Article  PubMed  Google Scholar 

  33. 33.

    Rahmer J, Blume U, Bornert P. Selective 3D ultrashort TE imaging: comparison of "dual-echo" acquisition and magnetization preparation for improving short-T 2 contrast. MAGMA. 2007;20:83–92.

    Article  PubMed  Google Scholar 

  34. 34.

    Falagas ME, Pitsouni EI, Malietzis GA, Pappas G. Comparison of pubmed, scopus, web of science, and google scholar: strengths and weaknesses. FASEB J. 2008;22:338–42.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Keereman V, Fierens Y, Broux T, De Deene Y, Lonneux M, Vandenberghe S. MRI-Based attenuation correction for PET/MRI using ultrashort echo time sequences. J Nucl Med. 2010;51:812–8.

    Article  PubMed  Google Scholar 

  36. 36.

    Low DA, Harms WB, Mutic S, Purdy JA. A technique for the quantitative evaluation of dose distributions. Med Phys. 1998;25:656–61.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Dice LR. Measures of the Amount of Ecologic Association between Species. Ecology. 1945;26:297–302.

    Article  Google Scholar 

  38. 38.

    Kosman E, Leonard KJ. Similarity coefficients for molecular markers in studies of genetic relationships between individuals for haploid, diploid, and polyploid species. Mol Ecol. 2005;14:415–24.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    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–35.

    Article  PubMed  Google Scholar 

  40. 40.

    Seibert TM, White NS, Kim G-Y, Moiseenko V, McDonald CR, Farid N, Bartsch H, Kuperman J, Karunamuni R, Marshall D. Distortion inherent to magnetic resonance imaging (MRI) can lead to geometric miss in radiosurgery planning. Pract Radiat Oncol. 2016;6(6):e319–28.

    Article  PubMed  Google Scholar 

  41. 41.

    Chen Y, Juttukonda M, Lee YZ, Su Y, Espinoza F, Lin W, Shen D, Lulash D, An H. MRI based attenuation correction for PET/MRI via MRF segmentation and sparse regression estimated CT. IEEE 11th International Symposium on Biomedical Imaging (ISBI). 2014;1364–67.

  42. 42.

    Huynh T, Gao Y, Kang J, Wang L, Zhang P, Lian J, Shen D. Estimating CT Image from MRI Data Using Structured Random Forest and Auto-context Model. IEEE Trans Med Imaging. 2016;35.1:174–83.

  43. 43.

    Andreasen D, Edmund JM, Zografos V, Menze BH, Van Leemput K. Computed tomography synthesis from magnetic resonance images in the pelvis using multiple random forests and auto-context features. SPIE Medical Imaging, Proceedings. 2016;9784.

  44. 44.

    Korsholm ME, Waring LW, Edmund JM. A criterion for the reliable use of MRI-only radiotherapy. Radiat Oncol. 2014;9:16.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Walker A, Liney G, Metcalfe P, Holloway L. MRI distortion: considerations for MRI based radiotherapy treatment planning. Australas Phys Eng Sci Med. 2014;37:103–13.

    Article  PubMed  Google Scholar 

  46. 46.

    Liney GP, Moerland MA. Magnetic resonance imaging acquisition techniques for radiotherapy planning. Semin Radiat Oncol. 2014;24:160–68.

  47. 47.

    Edmund JM, Andreasen D, Mahmood F, Van Leemput K. Cone beam CT guided treatment delivery and planning verification for magnetic resonance imaging only radiotherapy of the brain. Acta Oncol. 2015;54:1496–500.

    Article  PubMed  Google Scholar 

  48. 48.

    Andreasen D, Van Leemput K, Hansen RH, Andersen JA, 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–605.

    Article  PubMed  Google Scholar 

  49. 49.

    Andreasen D, Van Leemput K, Edmund JM. A patch-based pseudo-CT approach for MRI-only radiotherapy in the pelvis. Med Phys. 2016;43:4742–52.

    Article  PubMed  Google Scholar 

  50. 50.

    Korhonen J. Magnetic resonance imaging-based radiation therapy-Methods enabling the radiation therapy treatment planning workflow for prostate cancer patients by relying solely on MRI-based images throughout the process. Radiological Society of North America. Quantitative Imaging Data Warehouse (QIDW). Available at: 2015.

  51. 51.

    Korhonen J. Magnetic resonance imaging-based radiation therapy-Methods enabling the radiation therapy treatment planning workflow for prostate cancer patients by relying solely on MRI-based images throughout the process. Aalto University publication series DOCTORAL DISSERTATIONS 35/2015, ISBN 978-952-60-6123-8 (printed), ISBN 978-952-60-6124-5 (pdf), Helsinki, Finland.

  52. 52.

    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.

    Article  PubMed  Google Scholar 

  53. 53.

    Frantzen-Steneker M. SP-0022: Routine QA of an MR-only workflow. Radiother Oncol. 2015;115:S12.

    Article  Google Scholar 

  54. 54.

    Schad LR, Blüml S, Hawighorst H, Wenz F, Lorenz WJ. Radiosurgical treatment planning of brain metastases based on a fast, three-dimensional MR imaging technique. Magn Reson Imaging. 1994;12:811–9.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    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–9.

    Article  PubMed  Google Scholar 

  56. 56.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    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:1085–91.

    Article  PubMed  Google Scholar 

  58. 58.

    Korhonen J, Kapanen M, Sonke JJ, Wee L, Salli E, Keyrilainen J, Seppala T, Tenhunen M. Feasibility of MRI-based reference images for image-guided radiotherapy of the pelvis with either cone-beam computed tomography or planar localization images. Acta Oncol. 2015;54:889–95.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Toft-Jensen P, Geertsen P, Hansen RH, Kahlen A, Lindberg H, Kristensen BH. PO-0863 precision of nitinol prostate fiducial marker definition on T2 MRI in clinical practice. Radiother Oncol. 2012;103:S337–8.

    Article  Google Scholar 

  60. 60.

    Edmund J, Andreasen D, Kjer H, Van Leemput K. SP-0510: Dose planning based on MRI as the sole modality: Why, how and when? Radiother Oncol. 2015;115:S248–9.

    Article  Google Scholar 

  61. 61.

    Yu H, Caldwell C, Balogh J, Mah K. Toward magnetic resonance-only simulation: segmentation of bone in MR for radiation therapy verification of the head. Int J Radiat Oncol Biol Phys. 2014;89:649–57.

    Article  PubMed  Google Scholar 

  62. 62.

    Aitken A, Giese D, Tsoumpas C, Schleyer P, Kozerke S, Prieto C, Schaeffter T. Improved UTE-based attenuation correction for cranial PET-MR using dynamic magnetic field monitoring. Med Phys. 2014;41:012302.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Buerger C, Tsoumpas C, Aitken A, King AP, Schleyer P, Schulz V, Marsden PK, Schaeffter T. Investigation of MR-based attenuation correction and motion compensation for hybrid PET/MR. Nucl Sci IEEE Trans. 2012;59:1967–76.

    Article  Google Scholar 

  64. 64.

    Edmund JM, Kjer HM, Hansen RH. Auto-segmentation of bone in MRI-only based radiotherapy using ultra short echo time. Radiother Oncol. 2012;103 Suppl 1:S75.

    Article  Google Scholar 

  65. 65.

    Aasheim LB, Karlberg A, Goa PE, Haberg A, Sorhaug S, Fagerli UM, Eikenes L. PET/MR brain imaging: evaluation of clinical UTE-based attenuation correction. Eur J Nucl Med Mol Imaging. 2015;42:1439–46.

    Article  PubMed  Google Scholar 

  66. 66.

    Khateri P, Rad HS, Fathi A, Ay MR. Generation of attenuation map for MR-based attenuation correction of PET data in the head area employing 3D short echo time MR imaging. Nucl Instrum Methods Phys Res Section A: Accelerators Spectrometers Detectors Assoc Equip. 2013;702:133–6.

    CAS  Article  Google Scholar 

  67. 67.

    Cabello J, Lukas M, Förster S, Pyka T, Nekolla SG, Ziegler SI. MR-based attenuation correction using ultrashort-echo-time pulse sequences in dementia patients. J Nucl Med. 2015;56:423–9.

    Article  PubMed  Google Scholar 

  68. 68.

    Berker Y, Franke J, Salomon A, Palmowski M, Donker HCW, Temur Y, Mottaghy FM, Kuhl C, Izquierdo-Garcia D, Fayad ZA, et al. MRI-based attenuation correction for hybrid PET/MRI systems: A 4-class tissue segmentation technique using a combined ultrashort-echo-time/dixon MRI sequence. J Nucl Med. 2012;53:796–804.

    Article  PubMed  Google Scholar 

  69. 69.

    Juttukonda MR, Mersereau BG, Chen Y, Su Y, Rubin BG, Benzinger TL, Lalush DS, An H. MR-based attenuation correction for PET/MRI neurological studies with continuous-valued attenuation coefficients for bone through a conversion from R2* to CT-Hounsfield units. Neuroimage. 2015;112:160–8.

    Article  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Khateri P, Rad HS, Jafari AH, Ay MR. A novel segmentation approach for implementation of MRAC in head PET/MRI employing Short-TE MRI and 2-point Dixon method in a fuzzy C-means framework. Nucl Instrum Methods Phys Res Section A: Accelerators Spectrometers Detectors Assoc Equip. 2014;734:171–4.

    CAS  Article  Google Scholar 

  71. 71.

    Khateri P, Rad HS, Jafari AH, Kazerooni AF, Akbarzadeh A, Moghadam MS, Aryan A, Ghafarian P, Ay MR. Generation of a Four-Class Attenuation Map for MRI-Based Attenuation Correction of PET Data in the Head Area Using a Novel Combination of STE/Dixon-MRI and FCM Clustering. Mol Imaging Biol. 2015;17:1–9.

  72. 72.

    Su KH, Hu LZ, Stehning C, Helle M, Qian PJ, Thompson CL, Pereira GC, Jordan DW, Herrmann KA, Traughber M, et al. Generation of brain pseudo-CTs using an undersampled, single-acquisition UTE-mDixon pulse sequence and unsupervised clustering. Med Phys. 2015;42:4974–86.

    Article  PubMed  PubMed Central  Google Scholar 

  73. 73.

    Rank CM, Tremmel C, Hunemohr N, Nagel AM, Jakel O, Greilich S. MRI-based treatment plan simulation and adaptation for ion radiotherapy using a classification-based approach. Radiat Oncol. 2013;8:51.

    Article  PubMed  PubMed Central  Google Scholar 

  74. 74.

    Rank CM, Hunemohr N, Nagel AM, Rothke MC, Jakel O, Greilich S. MRI-based simulation of treatment plans for ion radiotherapy in the brain region. Radiother Oncol. 2013;109:414–8.

    Article  PubMed  Google Scholar 

  75. 75.

    Johansson A, Garpebring A, Karlsson M, Asklund T, Nyholm T. Improved quality of computed tomography substitute derived from magnetic resonance (MR) data by incorporation of spatial information--potential application for MR-only radiotherapy and attenuation correction in positron emission tomography. Acta Oncol. 2013;52:1369–73.

    Article  PubMed  Google Scholar 

  76. 76.

    Jonsson JH, Johansson A, Soderstrom K, Asklund T, Nyholm T. Treatment planning of intracranial targets on MRI derived substitute CT data. Radiother Oncol. 2013;108:118–22.

    Article  PubMed  Google Scholar 

  77. 77.

    Jonsson JH, Akhtari MM, Karlsson MG, Johansson A, Asklund T, Nyholm T. Accuracy of inverse treatment planning on substitute CT images derived from MR data for brain lesions. Radiat Oncol. 2015;10:1.

    Article  Google Scholar 

  78. 78.

    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.

    Article  PubMed  Google Scholar 

  79. 79.

    Kapanen M, Tenhunen M. T1/T2*-weighted MRI provides clinically relevant pseudo-CT density data for the pelvic bones in MRI-only based radiotherapy treatment planning. Acta Oncol. 2013;52:612–8.

    Article  PubMed  Google Scholar 

  80. 80.

    Korhonen J, Kapanen M, Keyriläinen J, Seppälä T, Tuomikoski L, Tenhunen M. Influence of MRI-based bone outline definition errors on external radiotherapy dose calculation accuracy in heterogeneous pseudo-CT images of prostate cancer patients. Acta Oncol. 2014;53:1100–6.

    Article  PubMed  Google Scholar 

  81. 81.

    Ay MR, Akbarzadeh A, Ahmadian A, Zaidi H. Classification of bones from MR images in torso PET-MR imaging using a statistical shape model. Nucl Instrum Methods Phys Res Section A: Accelerators Spectrometers Detectors Assoc Equip. 2014;734:196–200.

    Article  Google Scholar 

  82. 82.

    Yang X, Fei B. A skull segmentation method for brain MR images based on multiscale bilateral filtering scheme. SPIE Medical Imaging, Proceedings. 2010;7623.

  83. 83.

    Wagenknecht G, Kops ER, Kaffanke J, Tellmann L, Mottaghy F, Piroth MD, Herzog H. CT-based evaluation of segmented head regions for attenuation correction in MR-PET systems. IEEE Nuclear Science Symposuim & Medical Imaging Conference. 2010;2793–97.

  84. 84.

    Ribeiro AS, Kops ER, Herzog H, Almeida P. Skull segmentation of UTE MR images by probabilistic neural network for attenuation correction in PET/MR. Nucl Instrum Methods Phys Res Section A: Accelerators Spectrometers Detectors Assoc Equip. 2013;702:114–6.

    CAS  Article  Google Scholar 

  85. 85.

    Navalpakkam BK, Braun H, Kuwert T, Quick HH. Magnetic resonance–based attenuation correction for PET/MR hybrid imaging using continuous valued attenuation maps. Invest Radiol. 2013;48:323–32.

    Article  PubMed  Google Scholar 

  86. 86.

    Ribeiro AS, Kops ER, Herzog H, Almeida P. Hybrid approach for attenuation correction in PET/MR scanners. Nucl Instrum Methods Phys Res Section A: Accelerators Spectrometers Detectors Assoc Equip. 2014;734:166–70.

    Article  Google Scholar 

  87. 87.

    Yang X, Fei B. Multiscale segmentation of the skull in MR images for MRI-based attenuation correction of combined MR/PET. J Am Med Inform Assoc. 2013;20:1037–45.

    Article  PubMed  PubMed Central  Google Scholar 

  88. 88.

    Chan S-L, Gal Y, Jeffree RL, Fay M, Thomas P, Crozier S, Yang Z. Automated Classification of Bone and Air Volumes for Hybrid PET-MRI Brain Imaging. In Digital Image Computing: Techniques and Applications (DICTA), 2013 International Conference on. IEEE; 2013:1-8

  89. 89.

    Sjolund J, Jarlideni AE, Andersson M, Knutsson H, Nordstrom H. Skull Segmentation in MRI by a Support Vector Machine Combining Local and Global Features. In Pattern Recognition (ICPR), 2014 22nd International Conference on. IEEE; 2014:3274-3279

  90. 90.

    Akbarzadeh A, Gutierrez D, Baskin A, Ay M, Ahmadian A, Alam NR, Lövblad K-O, Zaidi H. Evaluation of whole-body MR to CT deformable image registration. J Appl Clin Med Phys. 2013;14(4):4163.

    CAS  PubMed  Google Scholar 

  91. 91.

    Uh J, Merchant TE, Li Y, Li X, Hua C. MRI-based treatment planning with pseudo CT generated through atlas registration. Med Phys. 2014;41:051711.

    Article  PubMed  PubMed Central  Google Scholar 

  92. 92.

    Stanescu T, Jans HS, Pervez N, Stavrev P, Fallone BG. A study on the magnetic resonance imaging (MRI)-based radiation treatment planning of intracranial lesions. Phys Med Biol. 2008;53:3579–93.

    CAS  Article  PubMed  Google Scholar 

  93. 93.

    Schreibmann E, Nye JA, Schuster DM, Martin DR, Votaw J, Fox T. MR-based attenuation correction for hybrid PET-MR brain imaging systems using deformable image registration. Med Phys. 2010;37:2101–9.

    Article  PubMed  Google Scholar 

  94. 94.

    Kops ER, Hautzel H, Herzog H, Antoch G, Shah NJ. Comparison of template-based versus CT-based attenuation correction for hybrid MR/PET scanners. IEEE Trans Nucl Sci. 2015;62:2115–21.

    Article  Google Scholar 

  95. 95.

    Whelan B, Kumar S, Dowling J, Begg J, Lambert J, Lim K, Vinod SK, Greer PB, Holloway L. Utilising pseudo-CT data for dose calculation and plan optimization in adaptive radiotherapy. Australas Phys Eng Sci Med. 2015;38(4):561–8.

    Article  PubMed  Google Scholar 

  96. 96.

    Hofmann M, Steinke F, Scheel V, Charpiat G, Farquhar J, Aschoff P, Brady M, Scholkopf B, Pichler BJ. MRI-Based Attenuation Correction for PET/MRI: A Novel Approach Combining Pattern Recognition and Atlas Registration. J Nucl Med. 2008;49:1875–83.

    Article  PubMed  Google Scholar 

  97. 97.

    Dowling JA, Sun J, Pichler P, Rivest-Hénault D, Ghose S, Richardson H, Wratten C, Martin J, Arm J, Best L. Automatic substitute computed tomography generation and contouring for magnetic resonance imaging (MRI)-alone external beam radiation therapy from standard MRI sequences. Int J Radiat Oncol Biol Phys. 2015;93:1144–53.

    Article  PubMed  Google Scholar 

  98. 98.

    Merida I, Costes N, Heckemann RA, Drzezga A, Forster S, Hammers A. Evaluation of several multi-atlas methods for PSEUDO-CT generation in brain MRI-PET attenuation correction. In IEEE 12th International Symposium on Biomedical Imaging (ISBI). 2015;1431–34.

  99. 99.

    Chowdhury N, Toth R, Chappelow J, Kim S, Motwani S, Punekar S, Lin H, Both S, Vapiwala N, Hahn S. Concurrent segmentation of the prostate on MRI and CT via linked statistical shape models for radiotherapy planning. Med Phys. 2012;39:2214–28.

    Article  PubMed  PubMed Central  Google Scholar 

  100. 100.

    Ladefoged CN, Benoit D, Law I, Holm S, Kjær A, Højgaard L, Hansen AE, Andersen FL. Region specific optimization of continuous linear attenuation coefficients based on UTE (RESOLUTE): application to PET/MR brain imaging. Phys Med Biol. 2015;60:8047.

    CAS  Article  PubMed  Google Scholar 

Download references


This work was initiated as a result of a symposium on MR-guided Radiotherapy in Heidelberg, Germany, 2015 organized by Prof. Dr. Oliver Jäkel.


This work was partly supported by a research grant from Varian Medical Systems, Inc.

Availability of data and materials

Not applicable.

Authors’ contributions

JME carried out the database search, collected and analyzed the data, and drafted the manuscript. TN contributed to the interpretation of the results and critically revised the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study reviews the results of published scientific literature and hence is not subject to ethnic approval. The data used for illustrative purposes were retrospective collected from studies and presentations previously published by the authors as described in the figure legends.

Author information



Corresponding author

Correspondence to Jens M. Edmund.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Edmund, J.M., Nyholm, T. A review of substitute CT generation for MRI-only radiation therapy. Radiat Oncol 12, 28 (2017).

Download citation


  • Magnetic Resonance Imaging
  • Magnetic Resonance Imaging Sequence
  • Mean Absolute Error
  • Compute Tomography Number
  • Electronic Portal Imaging Device