Ventricle Boundary in CT: Partial Volume Effect and Local Thresholds

We present a mathematical frame to carry out segmentation of cerebrospinal fluid (CSF) of ventricular region in computed tomography (CT) images in the presence of partial volume effect (PVE). First, the image histogram is fitted using the Gaussian mixture model (GMM). Analyzing the GMM, we find global threshold based on parameters of distributions for CSF, and for the combined white and grey matter (WGM). The parameters of distribution of PVE pixels on the boundary of ventricles are estimated by using a convolution operator. These parameters are used to calculate local thresholds for boundary pixels by the analysis of contribution of the neighbor pixels intensities into a PVE pixel. The method works even in the case of an almost unimodal histogram; it can be useful to analyze the parameters of PVE in the ground truth provided by the expert.