This monograph gives an introduction to the fundamentals of Markovian modeling in image segmentation as well as a brief overview of recent advances in the field. Segmentation is considered in a common framework, called image labeling, where the problem is reduced to assigning labels to pixels. In a probabilistic approach, label dependencies are modeled by Markov random fields (MRF) and an optimal labeling is determined by Bayesian estimation, in particular maximum a posteriori (MAP) estimation. The main advantage of MRF models is that prior information can be imposed locally through clique potentials. The primary goal is to demonstrate the basic steps to construct an easily applicable MRF segmentation model and further develop its multiscale and hierarchical implementations as well as their combination in a multilayer model. MRF models usually yield a non-convex energy function. The minimization of this function is crucial in order to find the most likely segmentation according to the MRF model. Besides classical optimization algorithms, like simulated annealing or deterministic relaxation, we also present recently introduced graph cutbased algorithms. We briefly discuss the possible parallelization techniques of simulated annealing, which allows efficient implementation on, e.g., GPU hardware without compromising convergence properties of the algorithms. While the main focus of this monograph is on generic model construction and related energy minimization methods, many sample applications are also presented to demonstrate the applicability of these models in real life problems such as remote sensing, biomedical imaging, change detection, and color-and motion-based segmentation. In real-life applications, parameter estimation is an important issue when implementing completely data-driven algorithms. Therefore some basic procedures, such as expectation-maximization, are also presented in the context of color image segmentation. Note: A sample implementation of the most important segmentation algorithms is available in grey scale at http://dx.doi.org/ 10.1561/2000000035 demogray and in color at http://dx.doi.org/ 10.1561/2000000035 democolor. ix 4 Graph Cut 4.1 Exact MAP of Binary MRFs via Standard Maxflow/Mincut 4.2 Solving Multilabel and Higher Order MRFs via GraphCut 4.3 An Example: Interactive Segmentation of Fluorescent Microscopic Images 5 Parameter Estimation and Sample Applications 5.1 Unsupervised Image Segmentation 5.2 Classification of Synthetic Aperture Radar Images 5.3 Multilayer MRF Models 6 Conclusion Acknowledgments References Dedication "To the memory of my mother" Zoltan Kato "To the memory of my beloved sister Elise who passed away in August 2012" Josiane Zerubia xi 1.1 Image Segmentation 3Ω is finite, although huge. A widely accepted standard, also motivated by the human visual system [121,162], is to construct this probability measure in a Bayesian framework [37,161,214]: We shall assume that we have a set of observed (Y ) and hidden (X) random variables. In our context, any observed value y ∈...