Image segmentation using voronoi polygons and MCMC, with application to muscle fibre images

Abstract

We investigate a Bayesian method for the segmentation of muscle fibre images. The images are reasonably well approximated by a Dirichlet tessellation, and so we use a deformable template model based on Voronoi polygons to represent the segmented image. We consider various prior distributions for the parameters and suggest an appropriate likelihood. Following the Bayesian paradigm, the mathematical form for the posterior distribution is obtained (up to an integrating constant). We introduce a Metropolis-Hastings algorithm and a reversible jump Markov chain Monte Carlo algorithm (RJMCMC) for simulation from the posterior when the number of polygons is fixed or unknown. The particular moves in the RJMCMC algorithm are birth, death and position/colour changes of the point process which determines the location of the polygons. Segmentation of the true image was carried out using the estimated posterior mode and posterior mean. A simulation study is presented which is helpful for tuning the hyperparameters and to assess the accuracy. The algorithms work well on a real image of a muscle fibre cross-section image, and an additional parameter, which models the boundaries of the muscle fibres, is included in the final model.Coloured tessellation, Markov chain Monte Carlo, point pattern, regularity, reversible jump, Strauss process,

A neural network approach to analyze cross-sections of muscle fibers in pathological images

“…Combinations of image processing pipelines have also been proposed to segment muscle fiber boundaries. Voronoi and reversible jump Markov chain Monte Carlo has also been used for muscle fiber segmentation [13]. Miazaki et al proposed an image-processing pipeline that consists of enhancement, noise reduction and binarization, and shape analysis of contours generated in the binarization process to detect muscle fiber boundaries [14].…”

A Stochastic Polygons Model for Glandular Structures in Colon Histology Images

“…On the right-hand side of (20) and (21), the first term corresponds to (A1) and the second term corresponds to (A2). Both terms are weighted by the relative area of the objects.…”

Charisma: An integrated approach to automatic H&E-stained skeletal muscle cell segmentation using supervised learning and novel robust clump splitting

“…Both works evaluate their method on histological breast tissue. A similar approach for skeletal muscle tissue is presented by Dryden et al (2006). The segmentation problem is expressed as the maximization of the posterior probability based on the tessellation of the muscle tissue using Markov Chain Monte Carlo methods.…”