Search

Using MRI for radiotherapy treatment planning and image guidance is appealing as it provides superior soft tissue information over CT scans and avoids possible systematic errors introduced by aligning MR to CT images. This study presents a method that generates Synthetic CT (MRCT) volumes by performing probabilistic tissue classification of voxels from MRI data using a single imaging sequence (T1 Dixon). The intensity overlap between different tissues on MR images, a major challenge for voxel-based MRCT generation methods, is addressed by adding bone shape information to an intensity-based classification scheme. A simple pelvic bone shape model, built from principal component analysis of pelvis shape from 30 CT image volumes, is fitted to the MR volumes. The shape model generates a rough bone mask that excludes air and covers bone along with some surrounding soft tissues. Air regions are identified and masked out from the tissue classification process by intensity thresholding outside the bone mask. A regularization term is added to the fuzzy c-means classification scheme that constrains voxels outside the bone mask from being assigned memberships in the bone class. MRCT image volumes are generated by multiplying the probability of each voxel being represented in each class with assigned attenuation values of the corresponding class and summing the result across all classes. The MRCT images presented intensity distributions similar to CT images with a mean absolute error of 13.7 HU for muscle, 15.9 HU for fat, 49.1 HU for intra-pelvic soft tissues, 129.1 HU for marrow and 274.4 HU for bony tissues across 9 patients. Volumetric modulated arc therapy (VMAT) plans were optimized using MRCT-derived electron densities, and doses were recalculated using corresponding CT-derived density grids. Dose differences to planning target volumes were small with mean/standard deviation of 0.21/0.42 Gy for D0.5cc and 0.29/0.33 Gy for D99%. The results demonstrate the accuracy of the method and its potential in supporting MRI only radiotherapy treatment planning. [ABSTRACT FROM AUTHOR]

Separating bone from air in MR data is one of the major challenges in using MR images to derive synthetic CT. The problem is further complicated when the anatomic regions filled with air are altered across scans due to air mobility, for instance, in pelvic regions, thereby the air regions estimated using an ultrashort echo time (UTE) sequence are invalid in other image series acquired for multispectral classification. This study aims to develop and investigate a female pelvic bone shape model to identify low intensity regions in MRI where air is unlikely to be present in support of synthetic CT generation without UTE imaging. CT scans of 30 patients were collected for the study, 17 of them also have corresponding MR scans. The shape model was built from the CT dataset, where the reference image was aligned to each of the training images using B-spline deformable registration. Principal component analysis was performed on B-spline coefficients for a compact model where shape variance was described by linear combination of principal modes. The model was applied to identify pelvic bone in MR images by deforming the corresponding MR data of the reference image to target MR images, where the search space of the deformation process was constrained within the subspace spanned by principal modes. The local minima in the search space were removed effectively by the shape model, thus supporting an efficient binary search for the optimal solution. We evaluated the model by its efficacy in identifying bone voxels and excluding air regions. The model was tested across the 17 patients that have corresponding MR scans using a leave-one-out cross validation. A simple model using the first leading principal mode only was found to achieve reasonable accuracy, where an averaged 87% of bone voxels were correctly identified. Finally dilation of the optimally fit bone mask by 5 mm was found to cover 96% of bone voxels while minimally impacting the overlap with air (below 0.4%). [ABSTRACT FROM AUTHOR]