Classification of bones from MR images in torso PET-MR imaging using a statistical shape model

There have been exclusive features for hybrid PET/MRI systems in comparison with its PET/CT counterpart in terms of reduction of radiation exposure, improved soft-tissue contrast and truly simultaneous and multi-parametric imaging capabilities. However, quantitative imaging on PET/MR is challenged by attenuation of annihilation photons through their pathway. The correction for photon attenuation requires the availability of patient-specific attenuation map, which accounts for the spatial distribution of attenuation coefficients of biological tissues. However, the lack of information on electron density in the MR signal poses an inherent difficulty to the derivation of the attenuation map from MR images. In other words, the MR signal correlates with proton densities and tissue relaxation properties, rather than with electron density and, as such, it is not directly related to attenuation coefficients. In order to derive the attenuation map from MR images at 511keV, various strategies have been proposed and implemented on prototype and commercial PET/MR systems. Segmentation-based methods generate an attenuation map by classification of T1-weighted or high resolution Dixon MR sequences followed by assignment of predefined attenuation coefficients to various tissue types. Intensity-based segmentation approaches fail to include bones in the attenuation map since the segmentation of bones from conventional MR sequences is a difficult task. Most MR-guided attenuation correction techniques ignore bones owing to the inherent difficulties associated with bone segmentation unless specialized MR sequences such as ultra-short echo (UTE) sequence are utilized. In this work, we introduce a new technique based on statistical shape modeling to segment bones and generate a four-class attenuation map. Our segmentation approach requires a torso bone shape model based on principle component analysis (PCA). A CT-based training set including clearly segmented bones of the torso region of 20 clinical studies was designed. Using this training set, a bone atlas was trained taking advantage of PCA analysis. Our active shape segmentation technique uses the trained shape model to segment bones from user defined initial seed points. The segmentation algorithm was evaluated using 10 clinical datasets (aligned MR and CT pairs). The resulting attenuation maps were compared to corresponding attenuation maps derived from CT resulting in a mean relative difference less than 7%.