Automatic Segmentation of Bone Selective MR Images for Visualization and Craniometry of the Cranial Vault

Solid-state MRI has been shown to provide a radiation-free alternative imaging strategy to CT. However, manual image segmentation to produce bone-selective MR-based 3D renderings is time and labor intensive, thereby acting as a bottleneck in clinical practice. The objective of this study was to evaluate an automatic multi-atlas segmentation pipeline for use on cranial vault images entirely circumventing prior manual intervention, and to assess concordance of craniometric measurements between pipeline produced MRI and CT-based 3D skull renderings.

Dual-RF, dual-echo, 3D UTE pulse sequence MR data were obtained at 3T on 30 healthy subjects along with low-dose CT images between December 2018 to January 2020 for this prospective study. The four-point MRI datasets (two RF pulse widths and two echo times) were combined to produce bone-specific images. CT images were thresholded and manually corrected to segment the cranial vault. CT images were then rigidly registered to MRI using mutual information. The corresponding cranial vault segmentations were then transformed to MRI. The "ground truth" segmentations served as reference for the MR images. Subsequently, an automated multi-atlas pipeline was used to segment the bone-selective images. To compare manually and automatically segmented MR images, the Dice similarity coefficient (DSC) and Hausdorff distance (HD) were computed, and craniometric measurements between CT and automated-pipeline MRI-based segmentations were examined via Lin's concordance coefficient (LCC).

Automated segmentation reduced the need for an expert to obtain segmentation. Average DSC was 90.86 ± 1.94%, and average 95th percentile HD was 1.65 ± 0.44 mm between ground truth and automated segmentations. MR-based measurements differed from CT-based measurements by 0.73-1.2 mm on key craniometric measurements. LCC for distances between CT and MR-based landmarks were vertex-basion: 0.906, left-right frontozygomatic suture: 0.780, and glabella-opisthocranium: 0.956 for the three measurements.

Good agreement between CT and automated MR-based 3D cranial vault renderings has been achieved, thereby eliminating the laborious manual segmentation process. Target applications comprise craniofacial surgery as well as imaging of traumatic injuries and masses involving both bone and soft tissue.