11.03.2014 Views

Segmentation of Stochastic Images using ... - Jacobs University

Segmentation of Stochastic Images using ... - Jacobs University

Segmentation of Stochastic Images using ... - Jacobs University

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!