Both Markov Random Field (MRF) and Marked Point Process (MPP) models have their limitations in image analysis. While MRF is a pixel-based representation, it is useful for imposing local constraints, but global constraints are not easily modeled. On the contrary, it is convenient to model global constraints, such as geometric shape and object interactions, within the MPP framework. However, such object-based MPP model has limited capability for imposing local constraints such as pixel-wise interactions. We have proposed a combined MRF/MPP model to incorporate local and global constraints within one unified framework. Unlike MRF or MPP, which only utilizes either local or global constraint for image analysis, our hybrid model addresses both constraints at pixel and object level simultaneously.
Let denote the segmented image,
the object field and
the observed image.
and
are defined on a lattice
;
represents a pixel site.
and
denote the space of all possible realizations of
and
.
is an object field realization, with each
being an object. Let
be the pixel site region projected by
onto the image lattice. Then we can define the segmentation of object
as
, i.e., the segmentation
restricted to the region
.
A non-homogeneous Gibbs process is introduced to represent the hybrid model. For this process, the neighborhood system is defined as: and
are neighbors, denoted
, if
intersects
. The energy function for the Gibbs distribution is
where is an energy potential describing how well object
fits the observed image,
is an energy potential based on the likelihood of the MAP segmentation for object
, and
is a geometric prior, describing the interactions between objects.
We use the Candy model [1] as the object model, which describes each line segment as an object. Usually three types of line segments are considered, which consist of free segment, single segment and double segment. Free segment has no ending point connected, single segment has one ending point connected while double segment has both ending points connected (illustrated as red segment in Fig. 1). Then a prior is proposed to describe the interactions between line segments. A detailed explanation of such prior can be found in [1].
Fig. 1. three types of line segments
We integrated such Candy model into our hybrid model and applied it to a silicate image. The target is to extract the silicate structure and obtain the corresponding segmentation. We perform optimization with our model using reversible jump Markov chain Monte Carlo (RJMCMC). We compared our hybrid model with original MPP model for object identification part and MRF based Graph Cuts method for segmentation. The results (Fig.2 ) show our proposed algorithm performs better than basic MPP at object level as well as MRF based Graph Cuts method at segmentation level.
[1] Point Processes for Unsupervised Line Network Extraction in Remote Sensing (Caroline Lacoste, Xavier Descombes, and JosianeZerubia, Transaction on PAMI, 2005)