Segmentation by discrete watersheds Part 1: Watershed cuts Jean - - PowerPoint PPT Presentation

segmentation by discrete watersheds part 1 watershed cuts
SMART_READER_LITE
LIVE PREVIEW

Segmentation by discrete watersheds Part 1: Watershed cuts Jean - - PowerPoint PPT Presentation

Segmentation by discrete watersheds Part 1: Watershed cuts Jean Cousty Four-Day Course on Mathematical Morphology in image analysis Bangalore 19-22 October 2010 J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology 1/36 An


slide-1
SLIDE 1

Segmentation by discrete watersheds Part 1: Watershed cuts

Jean Cousty Four-Day Course

  • n

Mathematical Morphology in image analysis

Bangalore 19-22 October 2010

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

1/36

slide-2
SLIDE 2

An applicative introduction to segmentation in medicine

Magnetic Resonance Imagery (MRI) is more and more used for cardiac diagnosis

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

2/36

slide-3
SLIDE 3

An applicative introduction to segmentation in medicine

A cardiac MRI examination includes three steps: Spatio-temporal acquisition (cin´ e MRI) im

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

2/36

slide-4
SLIDE 4

An applicative introduction to segmentation in medicine

A cardiac MRI examination includes three steps: Spatio-temporal acquisition (cin´ e MRI) Spatio-temporal acquisition during contrast agent injection (Perfusion) im

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

2/36

slide-5
SLIDE 5

An applicative introduction to segmentation in medicine

A cardiac MRI examination includes three steps: Spatio-temporal acquisition (cin´ e MRI) Spatio-temporal acquisition during contrast agent injection (Perfusion) Volumic acquisition after the evacuation of the contrast agent (delayed enhanced MRI) im

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

2/36

slide-6
SLIDE 6

Medical problem #1

Problem Visualizing objects of interests in 3D or 4D images rendu2

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

3/36

slide-7
SLIDE 7

Medical problem #2

Problem Determining measures useful for cardiac diagnosis

Infarcted volumes, ventricular volumes, ejection fraction, myocardial mass, movement . . .

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

4/36

slide-8
SLIDE 8

Technical problem

Problem Segmentation of object of interest

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

5/36

slide-9
SLIDE 9

Technical problem

Problem Segmentation of object of interest A morphological solution Watershed

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

5/36

slide-10
SLIDE 10

Watershed: introduction

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

6/36

slide-11
SLIDE 11

Watershed: introduction

For topographic purposes, the watershed has been studied since the 19th century (Maxwell, Jordan, . . . )

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

6/36

slide-12
SLIDE 12

Watershed: introduction

One hundred years later (1978), it was introduced by Digabel and Lantu´ ejoul for image segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

6/36

slide-13
SLIDE 13

Watershed: introduction

One hundred years later (1978), it was introduced by Digabel and Lantu´ ejoul for image segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

6/36

slide-14
SLIDE 14

Watershed: introduction

One hundred years later (1978), it was introduced by Digabel and Lantu´ ejoul for image segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

6/36

slide-15
SLIDE 15

Watershed: introduction

One hundred years later (1978), it was introduced by Digabel and Lantu´ ejoul for image segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

6/36

slide-16
SLIDE 16

Watershed: introduction

One hundred years later (1978), it was introduced by Digabel and Lantu´ ejoul for image segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

6/36

slide-17
SLIDE 17

Watershed: introduction

One hundred years later (1978), it was introduced by Digabel and Lantu´ ejoul for image segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

6/36

slide-18
SLIDE 18

Watershed: problem #1

Problem How to define the watershed of digital image?

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

7/36

slide-19
SLIDE 19

Watershed: problem #1

Problem How to define the watershed of digital image? Which mathematical framework(s)?

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

7/36

slide-20
SLIDE 20

Watershed: problem #1

Problem How to define the watershed of digital image? Which mathematical framework(s)? Which properties?

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

7/36

slide-21
SLIDE 21

Watershed: problem #1

Problem How to define the watershed of digital image? Which mathematical framework(s)? Which properties? Which algorithms ?

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

7/36

slide-22
SLIDE 22

Watershed: problem #2

Problem In practice: over-segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

8/36

slide-23
SLIDE 23

Over-segmentation and region merging

Solution 1 Region merging methods consist of improving an initial segmentation by progressively merging pairs of neighboring regions

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

9/36

slide-24
SLIDE 24

Over-segmentation and region merging

Solution 1 Region merging methods consist of improving an initial segmentation by progressively merging pairs of neighboring regions Example : delayed enhanced cardiac MRI [DOUBLIER03]

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

9/36

slide-25
SLIDE 25

Over-segmentation and region merging

Solution 1 Region merging methods consist of improving an initial segmentation by progressively merging pairs of neighboring regions Example : delayed enhanced cardiac MRI [DOUBLIER03]

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

9/36

slide-26
SLIDE 26

Over-segmentation and region merging

Solution 1 Region merging methods consist of improving an initial segmentation by progressively merging pairs of neighboring regions Example : delayed enhanced cardiac MRI [DOUBLIER03]

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

9/36

slide-27
SLIDE 27

Over-segmentation and region merging

Solution 1 Region merging methods consist of improving an initial segmentation by progressively merging pairs of neighboring regions

infarctus myocarde cavité sanguine

Example : delayed enhanced cardiac MRI [DOUBLIER03]

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

9/36

slide-28
SLIDE 28

Over-segmentation

Solution 2 Seeded watershed (or marker based watershed)

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

10/36

slide-29
SLIDE 29

Over-segmentation

Solution 2 Seeded watershed (or marker based watershed) Methodology proposed by Beucher and Meyer (1993)

1 Recognition 2 Delineation (generally done by watershed) 3 Smoothing

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

10/36

slide-30
SLIDE 30

Over-segmentation

Solution 2 Seeded watershed (or marker based watershed) Methodology proposed by Beucher and Meyer (1993)

1 Recognition 2 Delineation (generally done by watershed) 3 Smoothing

Semantic information taken into account at steps 1 and 3 To kow more about this framework, wait for the second lecture of today

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

10/36

slide-31
SLIDE 31

Outline

1 Defining discrete watersheds is difficult

Grayscale image as vertex weighted graphs Region merging problems

2 Watershed in edge-weighted graphs

Watershed cuts: definition and consistency Minimum spanning forests: watershed optimality

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

11/36

slide-32
SLIDE 32

Defining discrete watersheds is difficult

Can we draw a watershed of this image?

2 2 2 2 2 40 30 30 30 40 40 20 20 20 40 40 40 20 40 40 1 5 20 5 1

Image equipped with the 4-adjacency

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

12/36

slide-33
SLIDE 33

Defining discrete watersheds is difficult

Can we draw a watershed of this image?

A A A A A 40 30 30 30 40 40 20 20 20 40 40 40 20 40 40 B 5 20 5 C

Image equipped with the 4-adjacency Label the pixels according to catchment basins letters A,B and C

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

12/36

slide-34
SLIDE 34

Defining discrete watersheds is difficult

Possible drawings

A A A A A 40 30 30 30 40 40 20 20 20 40 40 40 20 40 40 B 5 20 5 C A A A A A A A A A A 40 20 20 20 40 B B 20 C C B B 20 C C

Topographical watershed

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

13/36

slide-35
SLIDE 35

Defining discrete watersheds is difficult

Possible drawings

A A A A A 40 30 30 30 40 40 20 20 20 40 40 40 20 40 40 B 5 20 5 C A A A A A A A A A A C C C C C B B C C C B B C C C

Flooding from the minima

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

13/36

slide-36
SLIDE 36

Defining discrete watersheds is difficult

Possible drawings

A A A A A 40 30 30 30 40 40 20 20 20 40 40 40 20 40 40 B 5 20 5 C A A A A A A A A A A 40 A A A 40 B 40 A 40 C B B 20 C C

Flooding with divide

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

13/36

slide-37
SLIDE 37

Defining discrete watersheds is difficult

Possible drawings

A A A A A 40 30 30 30 40 40 20 20 20 40 40 40 20 40 40 B 5 20 5 C A A A A A 40 30 30 30 40 B B 20 C C B B 20 C C B B 20 C C

Topological watershed

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

13/36

slide-38
SLIDE 38

Defining discrete watersheds is difficult

Possible drawings

A A A A A 40 30 30 30 40 40 20 20 20 40 40 40 20 40 40 B 5 20 5 C A A A A A 40 30 30 30 40 B B 20 C C B B 20 C C B B 20 C C

Topological watershed Conclusion Not easy to define watersheds of digital images

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

13/36

slide-39
SLIDE 39

Defining discrete watersheds is difficult

Region merging: Problem #1

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

14/36

slide-40
SLIDE 40

Defining discrete watersheds is difficult

Region merging: Problem #1

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

14/36

slide-41
SLIDE 41

Defining discrete watersheds is difficult

Region merging: Problem #1

A B C D

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

14/36

slide-42
SLIDE 42

Defining discrete watersheds is difficult

Region merging: Problem #1

A B C D

x

y

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

14/36

slide-43
SLIDE 43

Defining discrete watersheds is difficult

Region merging: Problem #1

A B C D

Problem : “When 3 regions meet”, [PAVLIDIS-77]

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

14/36

slide-44
SLIDE 44

Defining discrete watersheds is difficult

Region merging: Problem #1

adjacence directe adjacence indirecte ?

Problem : “When 3 regions meet”, [PAVLIDIS-77] Is there some adjacency relations (graphs) for which any pair of neighboring regions can always be merged, while preserving all

  • ther regions?
  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

14/36

slide-45
SLIDE 45

Defining discrete watersheds is difficult

Region merging: Problem # 2

A cleft is a set of vertices from which a point cannot be removed while leaving unchanged the number of connected components of its complement

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

15/36

slide-46
SLIDE 46

Defining discrete watersheds is difficult

Region merging: Problem # 2

A cleft is a set of vertices from which a point cannot be removed while leaving unchanged the number of connected components of its complement A cleft is thin if all its vertices are adjacent to its complement

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

15/36

slide-47
SLIDE 47

Defining discrete watersheds is difficult

Region merging: Problem # 2

A cleft is a set of vertices from which a point cannot be removed while leaving unchanged the number of connected components of its complement A cleft is thin if all its vertices are adjacent to its complement

Problem : Thick cleft (or binary watershed)

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

15/36

slide-48
SLIDE 48

Defining discrete watersheds is difficult

Region merging: Problem # 2

A cleft is a set of vertices from which a point cannot be removed while leaving unchanged the number of connected components of its complement A cleft is thin if all its vertices are adjacent to its complement

Problem : Thick cleft (or binary watershed) Is there some graphs in which any cleft is thin?

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

15/36

slide-49
SLIDE 49

Defining discrete watersheds is difficult

The familly of watersheds

Topological watersheds Only watersheds that preserve the altitudes of the passes On pixels Fusion graphs Link between thinness, region merging, and watersheds On edges Watershed cuts Optimality, drop of water principle Power watersheds Framework for seeded image segmentation (graph cuts, random walker, . . .) Energy minimization q = 2 = ⇒ uniqueness Ultrametric watersheds Hierarchical segmentation On complexes Simplicial stacks Link between collapse, watersheds and optimal spanning forests

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

16/36

slide-50
SLIDE 50

Watershed in edge-weighted graphs

Watershed in edge-weighted graphs

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

17/36

slide-51
SLIDE 51

Watershed in edge-weighted graphs

Watershed in edge-weighted graphs

Let G = (V , E) be a graph. Let F be a map from E to R

5 5 5 5 5 2 4 4 3 2 4 5 2 5 3 4 4 3

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

17/36

slide-52
SLIDE 52

Watershed in edge-weighted graphs

Image and edge-weighted graph

For applications to image analysis V is the set of pixels E corresponds to an adjacency relation on V , (e.g., 4- or 8-adjacency in 2D) F is a “gradient” of I: The altitude of u, an edge between two pixels x and y, represents the dissimilarity between x and y

F(u) = |I(x) − I(y)|.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

18/36

slide-53
SLIDE 53

Watershed in edge-weighted graphs

Regional minima

00 1 2 2 2 3 3 4 4 5 5 5 7 7 7 6 7 2

Definition A subgraph X of G is a minimum of F (at altitude k) if: X is connected; and k is the altitude of any edge of X; and the altitude of any edge adjacent to X is strictly greater than k

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

19/36

slide-54
SLIDE 54

Watershed in edge-weighted graphs

Extension

a subgraph X Definition (from Def. 12, (Ber05)) Let X and Y be two non-empty subgraphs of G We say that Y is an extension of X (in G) if X ⊆ Y and if any component of Y contains exactly one component of X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

20/36

slide-55
SLIDE 55

Watershed in edge-weighted graphs

Extension

a subgraph X an extension Y of X Definition (from Def. 12, (Ber05)) Let X and Y be two non-empty subgraphs of G We say that Y is an extension of X (in G) if X ⊆ Y and if any component of Y contains exactly one component of X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

20/36

slide-56
SLIDE 56

Watershed in edge-weighted graphs

Extension

a subgraph X an extension Y of X Definition (from Def. 12, (Ber05)) Let X and Y be two non-empty subgraphs of G We say that Y is an extension of X (in G) if X ⊆ Y and if any component of Y contains exactly one component of X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

20/36

slide-57
SLIDE 57

Watershed in edge-weighted graphs

Extension

a subgraph X an extension Y of X Definition (from Def. 12, (Ber05)) Let X and Y be two non-empty subgraphs of G We say that Y is an extension of X (in G) if X ⊆ Y and if any component of Y contains exactly one component of X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

20/36

slide-58
SLIDE 58

Watershed in edge-weighted graphs

Extension

a subgraph X an extension Y of X Definition (from Def. 12, (Ber05)) Let X and Y be two non-empty subgraphs of G We say that Y is an extension of X (in G) if X ⊆ Y and if any component of Y contains exactly one component of X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

20/36

slide-59
SLIDE 59

Watershed in edge-weighted graphs

Extension

a subgraph X an extension Y of X Definition (from Def. 12, (Ber05)) Let X and Y be two non-empty subgraphs of G We say that Y is an extension of X (in G) if X ⊆ Y and if any component of Y contains exactly one component of X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

20/36

slide-60
SLIDE 60

Watershed in edge-weighted graphs

Extension

a subgraph X an extension Y of X Definition (from Def. 12, (Ber05)) Let X and Y be two non-empty subgraphs of G We say that Y is an extension of X (in G) if X ⊆ Y and if any component of Y contains exactly one component of X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

20/36

slide-61
SLIDE 61

Watershed in edge-weighted graphs

Extension

a subgraph X an extension Y of X Definition (from Def. 12, (Ber05)) Let X and Y be two non-empty subgraphs of G We say that Y is an extension of X (in G) if X ⊆ Y and if any component of Y contains exactly one component of X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

20/36

slide-62
SLIDE 62

Watershed in edge-weighted graphs

Graph cut

a subgraph X a (graph) cut S for X Definition (Graph cut) Let X be a subgraph of G and S ⊆ E, a set of edges. We say that S is a (graph) cut for X if S is an extension of X and if S is minimal for this property

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

21/36

slide-63
SLIDE 63

Watershed in edge-weighted graphs

Watershed cut

The church of Sorbier (a topographic intuition) Definition (drop of water principle) The set S ⊆ E is a watershed cut of F if S is an extension of M(F) and if for any u = {x0, y0} ∈ S, there exist x0, . . . , xn and y0, . . . , ym, two descending paths in S such that:

1 xn and ym are vertices of two distinct minima of F; and 2 F(u) ≥ F({x0, x1}) if n > 0 and F(u) ≥ F({y0, y1}) if m > 0

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

22/36

slide-64
SLIDE 64

Watershed in edge-weighted graphs

Watershed cut: example

2 6 6 2 6 5 5 7 3 6 5 3 1 2 5 6 7 5 7 3 4 5 3 8

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

23/36

slide-65
SLIDE 65

Watershed in edge-weighted graphs

Watershed cut: example

2 6 6 2 6 5 5 7 3 6 5 3 1 2 8 5 6 7 5 7 3 4 5 3

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

23/36

slide-66
SLIDE 66

Watershed in edge-weighted graphs

Watershed cut: example

2 6 6 2 6 5 5 7 3 6 5 3 1 2 8 5 6 7 5 7 3 4 5 3

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

23/36

slide-67
SLIDE 67

Watershed in edge-weighted graphs

Watershed cut: example

2 6 6 2 6 5 5 7 3 6 5 3 1 2 8 5 6 7 5 7 3 4 5 3

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

23/36

slide-68
SLIDE 68

Watershed in edge-weighted graphs

Watershed cut: example

2 6 6 2 6 5 5 7 3 6 5 3 1 2 8 5 6 7 5 7 3 4 5 3

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

23/36

slide-69
SLIDE 69

Watershed in edge-weighted graphs

Watershed cut: example

2 6 6 2 6 5 5 7 3 6 5 3 1 2 8 5 6 7 5 7 3 4 5 3

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

23/36

slide-70
SLIDE 70

Watershed in edge-weighted graphs

Steepest descent

2 6 6 2 6 5 5 7 3 6 5 3 1 2 8 5 6 7 5 7 3 4 5 3 Definition Let π = x0, . . . , xl be a path in G. The path π is a path with steepest descent for F if: ∀i ∈ [1, l], F({xi−1, xi}) = min{xi−1,y}∈E F({xi−1, y})

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

24/36

slide-71
SLIDE 71

Watershed in edge-weighted graphs

Catchment basins by a steepest descent property

Definition Let S be a cut for M(F), the minima of F We say that S is a basin cut of F if, from each point of V to M(F), there exists, in the graph induced by S, a path with steepest descent for F

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

25/36

slide-72
SLIDE 72

Watershed in edge-weighted graphs

Catchment basins by a steepest descent property

Theorem (consistency) An edge-set S ⊆ E is a basin cut of F if and only if S is a watershed cut of F

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

26/36

slide-73
SLIDE 73

Watershed in edge-weighted graphs

Illustration to grayscale image segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

27/36

slide-74
SLIDE 74

Watershed in edge-weighted graphs

Illustration to grayscale image segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

27/36

slide-75
SLIDE 75

Watershed in edge-weighted graphs

Illustration to mesh segmentation

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

28/36

slide-76
SLIDE 76

Watershed in edge-weighted graphs

Watershed optimality?

Problem Are watersheds optimal segmentations? Which combinatorial optimization problem do they solve?

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

29/36

slide-77
SLIDE 77

Watershed in edge-weighted graphs

Relative minimum spanning forest: an image intuitition

cut forest spanning regions

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

30/36

slide-78
SLIDE 78

Watershed in edge-weighted graphs

Relative forest: a botanical intuition

A tree (Lal Bagh)

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

31/36

slide-79
SLIDE 79

Watershed in edge-weighted graphs

Relative forest: a botanical intuition

Cuting the roots yield a forest of several trees

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

31/36

slide-80
SLIDE 80

Watershed in edge-weighted graphs

Relative forest: a botanical intuition

Roots may contain cycles

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

31/36

slide-81
SLIDE 81

Watershed in edge-weighted graphs

Relative forest

a subgraph X a forest Y relative to X Definition Let X and Y be two non-empty subgraphs of G. We say that Y is a forest relative to X if:

1 Y is an extension of X; and 2 any cycle of Y is also a cycle of X

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

32/36

slide-82
SLIDE 82

Watershed in edge-weighted graphs

Minimum spanning forest

The weight of a forest Y is the sum of its edge weights i.e.

  • u∈E(Y ) F(u).
  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

33/36

slide-83
SLIDE 83

Watershed in edge-weighted graphs

Minimum spanning forest

The weight of a forest Y is the sum of its edge weights i.e.

  • u∈E(Y ) F(u).

Definition We say that Y is a minimum spanning forest (MSF) relative to X

if Y is a spanning forest relative to X and if the weight of Y is less than or equal to the weight of any other spanning forest relative to X

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

33/36

slide-84
SLIDE 84

Watershed in edge-weighted graphs

Minimum spanning forest: example

5 5 5 5 5 2 4 4 3 2 4 5 2 5 3 4 4 3

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

34/36

slide-85
SLIDE 85

Watershed in edge-weighted graphs

Minimum spanning forest: example

5 5 5 5 5 2 4 4 3 2 4 5 2 5 3 4 3 4

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

34/36

slide-86
SLIDE 86

Watershed in edge-weighted graphs

Minimum spanning forest: example

5 5 5 5 5 2 4 4 3 2 4 5 2 5 3 4 3 4

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

34/36

slide-87
SLIDE 87

Watershed in edge-weighted graphs

Minimum spanning forest: example

5 5 5 5 5 2 4 4 3 2 4 5 2 5 3 4 3 4

If Y is a MSF relative to X, there exists a unique cut S for Y and this cut is also a cut for X;

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

34/36

slide-88
SLIDE 88

Watershed in edge-weighted graphs

Minimum spanning forest: example

5 5 5 5 5 2 4 4 3 2 4 5 2 5 3 4 3 4

If Y is a MSF relative to X, there exists a unique cut S for Y and this cut is also a cut for X; In this case, we say that S is a MSF cut for X.

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

34/36

slide-89
SLIDE 89

Watershed in edge-weighted graphs

Watershed optimality

Theorem An edge-set S ⊆ E is a MSF cut for the minima of F if and only if S is a watershed cut of F

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

35/36

slide-90
SLIDE 90

Watershed in edge-weighted graphs

Minimum spanning tree

Computing a MSF ⇔ computing a minimum spanning tree

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

36/36

slide-91
SLIDE 91

Watershed in edge-weighted graphs

Minimum spanning tree

Computing a MSF ⇔ computing a minimum spanning tree Best algorithm [CHAZEL00]: quasi-linear time

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

36/36

slide-92
SLIDE 92

Watershed in edge-weighted graphs

Minimum spanning tree

Computing a MSF ⇔ computing a minimum spanning tree Best algorithm [CHAZEL00]: quasi-linear time Problem Can we reach a better complexity for computing watershed cuts?

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

36/36

slide-93
SLIDE 93

Watershed in edge-weighted graphs

Minimum spanning tree

Computing a MSF ⇔ computing a minimum spanning tree Best algorithm [CHAZEL00]: quasi-linear time Problem Can we reach a better complexity for computing watershed cuts? A morphological solution To know the answer, come back after the coffee break

  • J. Serra, J. Cousty, B.S. Daya Sagar : Course on Math. Morphology

36/36