Babies born prematurely are at increased risk of adverse neurodevelopmental outcomes. Recent advances suggest that measurement of brain volumes can help in defining biomarkers for neurodevelopmental outcome. These techniques rely on an accurate segmentation of the MRI data. However, due to lack of contrast, partial volume (PV) effect, the existence of both hypo- and hyper-intensities and significant natural and pathological anatomical variability, the segmentation of neonatal brain MRI is challenging. We propose a pipeline for image segmentation that uses a novel multi-model Maximum a posteriori Expectation Maximisation (MAP-EM) segmentation algorithm with a prior over both intensities and the tissue proportions, a B0 inhomogeneity correction, and a spatial homogeneity term through the use of a Markov Random Field. This robust and adaptive technique enables the segmentation of images with high anatomical disparity from a normal population. Furthermore, the proposed method implicitly models Partial Volume, mitigating the problem of neonatal white/grey matter intensity inversion. Experiments performed on a clinical cohort show expected statistically significant correlations with gestational age at birth and birthweight. Furthermore, the proposed method obtains statistically significant improvements in Dice scores when compared to the a Maximum Likelihood EM algorithm.