Pseudo-CT generation by conditional inference random forest for MRI-based radiotherapy treatment planning

Dose calculation from MRI is a topical issue. New treatment systems combining a linear accelerator with a MRI have been recently being developed. MRI has good soft tissue contrast without ionizing radiation exposure. However, unlike CT, MRI does not provide electron density information necessary for dose calculation. We propose in this paper a machine learning method to simulate a CT from a target MRI and co-registered CT-MRI training set. Ten prostate MR and CT images have been considered. Firstly, a reference image was randomly selected in the training set. A common space has been built thanks to affine registrations between the training set and the reference image. Multiscale image descriptors such as spatial information, gradients and texture features were extracted from MRI patches at dilïerent levels of a Gaussian pyramid and used as voxel-wise characteristics in the learning scheme. A Conditional Inference Random Forest (CIRF) modelled the relation between MRI descriptors and CT patches. For validation, test images were spatially normalized and the same descriptors were computed to generate a new pCT. Leave-one out experiments were performed. We obtained a MAE = 45.79 (pCT vs CT). Dose volume histograms inside PTV and organs at risk are in close agreement. The D98% was 0.45 % (inside PTV) and the 3D gamma pass rate (1mm, 1%) was 99,2%. Our method has better results than direct bulk assignment. And the results suggest that the method may be used for dose calculations in an MR based planning system.

[1]  Jurgen Fripp,et al.  Automatic Substitute Computed Tomography Generation and Contouring for Magnetic Resonance Imaging (MRI)-Alone External Beam Radiation Therapy From Standard MRI Sequences. , 2015, International journal of radiation oncology, biology, physics.

[2]  Ninon Burgos,et al.  Robust CT Synthesis for Radiotherapy Planning: Application to the Head and Neck Region , 2015, MICCAI.

[3]  Yaozong Gao,et al.  Estimating CT Image From MRI Data Using Structured Random Forest and Auto-Context Model , 2016, IEEE Transactions on Medical Imaging.

[4]  Jan J W Lagendijk,et al.  The feasibility of utilizing pseudo CT-data for online MRI based treatment plan adaptation for a stereotactic radiotherapy treatment of spinal bone metastases , 2014, Physics in medicine and biology.

[5]  Romain Viard,et al.  Monte Carlo calculation based on hydrogen composition of the tissue for MV photon radiotherapy. , 2015, Journal of applied clinical medical physics.

[6]  D. Low,et al.  A technique for the quantitative evaluation of dose distributions. , 1998, Medical physics.

[7]  Juha Korhonen,et al.  Dosimetric characterization of MRI-only treatment planning for brain tumors in atlas-based pseudo-CT images generated from standard T1-weighted MR images. , 2016, Medical physics.

[8]  Tiina Seppälä,et al.  A dual model HU conversion from MRI intensity values within and outside of bone segment for MRI-based radiotherapy treatment planning of prostate cancer. , 2013, Medical physics.

[9]  Jan J W Lagendijk,et al.  MRI/linac integration. , 2008, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[10]  Jan-Jakob Sonke,et al.  A fast algorithm for gamma evaluation in 3D. , 2007, Medical physics.

[11]  Adam Johansson,et al.  CT substitute derived from MRI sequences with ultrashort echo time. , 2011, Medical physics.

[12]  Peter B. Greer,et al.  Structure-Guided Nonrigid Registration of CT-MR Pelvis Scans with Large Deformations in MR-Based Image Guided Radiation Therapy , 2013, CLIP.

[13]  Olivier Salvado,et al.  An atlas-based electron density mapping method for magnetic resonance imaging (MRI)-alone treatment planning and adaptive MRI-based prostate radiation therapy. , 2012, International journal of radiation oncology, biology, physics.

[14]  Mika Kapanen,et al.  T1/T2*-weighted MRI provides clinically relevant pseudo-CT density data for the pelvic bones in MRI-only based radiotherapy treatment planning , 2013, Acta oncologica.

[15]  Adam Johansson,et al.  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 , 2013, Acta oncologica.

[16]  Daniel Rueckert,et al.  Patch-Based Segmentation without Registration: Application to Knee MRI , 2013, MLMI.

[17]  Steve Webb,et al.  Radiotherapy treatment planning of prostate cancer using magnetic resonance imaging alone. , 2003, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[18]  Jurgen Fripp,et al.  Robust inverse-consistent affine CT-MR registration in MRI-assisted and MRI-alone prostate radiation therapy , 2015, Medical Image Anal..

[19]  Olivier Salvado,et al.  MRI-guided prostate radiation therapy planning: Investigation of dosimetric accuracy of MRI-based dose planning. , 2011, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.