Segmentation of Stochastic Images using ... - Jacobs University
Segmentation of Stochastic Images using ... - Jacobs University
Segmentation of Stochastic Images using ... - Jacobs University
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
2.2 Random Walker <strong>Segmentation</strong><br />
Definition The space <strong>of</strong> functions <strong>of</strong> bounded variation is<br />
{<br />
∫<br />
}<br />
BV(D) = u ∈ L 1 (D) : |Du|dx < ∞<br />
D<br />
. (2.2)<br />
Following [17] and <strong>using</strong> the Lebesgue decomposition theorem (see [32]) the derivative <strong>of</strong> a BVfunction<br />
decomposes into three parts, the absolutely continuous part ∇udx, the jump part D j u and<br />
the Cantor part D c u. This leads us to the definition <strong>of</strong> the class <strong>of</strong> special BV-functions:<br />
Definition The class <strong>of</strong> BV-function for which the Cantor part vanishes, i.e. D c u = 0, is called the<br />
space <strong>of</strong> special functions <strong>of</strong> bounded variation (SBV).<br />
Based on the space SBV we define the space <strong>of</strong> generalized functions <strong>of</strong> bounded variation (GSBV):<br />
Definition The space GSBV consists <strong>of</strong> all functions u ∈ L 1 (D) satisfying<br />
i.e. the truncated function belongs to SBV for all T .<br />
∀T > 0 : u T = sign(u)max(T,|u|) ∈ SBV , (2.3)<br />
In particular, the spaces SBV and GSBV are useful, when we introduce Mumford-Shah segmentation<br />
and the related Ambrosio-Tortorelli approximation.<br />
2.2 Random Walker <strong>Segmentation</strong><br />
Random walker segmentation performs on a single image u : D → IR defined on the image domain D.<br />
Since random walker segmentation is based on pixel values only, we need no additional assumptions<br />
about the smoothness <strong>of</strong> the images. The main idea <strong>of</strong> random walker segmentation is that the user<br />
prescribes a set <strong>of</strong> seed points for the object and the background. From the remaining unseeded<br />
points, random walks start and the percentage <strong>of</strong> random walks reaching the object seeds is the<br />
probability <strong>of</strong> the pixel to belong to the object.<br />
Before we begin to introduce random walker segmentation, we have to define notation for the<br />
graph representation <strong>of</strong> the image. A graph G is a pair G = (V,E) containing vertices or nodes<br />
v ∈ V and edges e ∈ E ⊂ V ×V. We denote an edge connecting the vertices v i and v j by e i j and<br />
identify e i j with e ji , because we are interested in nondirectional graphs only. Every edge has a<br />
weight w(e i j ) =: w i j that describes the costs for <strong>using</strong> the edge. A graph containing edge weights is<br />
a weighted graph. The summation <strong>of</strong> all edge weights for a node i,<br />
d i =<br />
∑<br />
{ j∈V :e i j ∈E}<br />
w(e i j ) , (2.4)<br />
is the degree <strong>of</strong> the node i.<br />
For random walker segmentation, we identify the image u with a graph G. The pixels <strong>of</strong> the digital<br />
image are the graph nodes and every pixel (respectively node) is connected to the neighboring nodes<br />
by a weighted edge. Fig. 2.2 shows a graph corresponding to a 3 × 3 image. For random walker<br />
segmentation the graph weights are<br />
w(e i j ) = exp ( −β(g i − g j ) 2) , (2.5)<br />
where (g i − g j ) 2 is the normalized difference between the image intensities at position i and j:<br />
(g i − g j ) 2 =<br />
The parameter β is the only free parameter that the user chooses.<br />
(u i − u j ) 2<br />
max {k,l∈V :ekl ∈E}(u k − u l ) 2 . (2.6)<br />
9