Morphologische Filter
Morphologische Filter
Morphologische Filter
Sie wollen auch ein ePaper? Erhöhen Sie die Reichweite Ihrer Titel.
YUMPU macht aus Druck-PDFs automatisch weboptimierte ePaper, die Google liebt.
<strong>Morphologische</strong> <strong>Filter</strong><br />
Vorlesung FH-Hagenberg<br />
DSB<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Mathematische Morphologie<br />
• Binäre Morphologie<br />
– Strukturelement<br />
– Grundlegende Operatoren<br />
• Erosion, Dilation<br />
– Abgeleitete Operatoren<br />
• Open, close<br />
• Grauwert-Morphologie<br />
Foliensatz W. Burger, W. Backfrieder<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Ziel: Strukturabhängige Segmentierung<br />
artefacts<br />
(a) (b) (c)<br />
(a) fingerprint grey-level, (b) thresholded image, (c)<br />
morphological filter, artefacts removed<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Binärbilder<br />
nur 2 Pixelwerte:<br />
1<br />
0<br />
p th<br />
BLACK WHITE<br />
0 1 (255)<br />
Binarisierung durch Schwellwertbildung:<br />
q<br />
0 255<br />
P th ...Schwellwert<br />
„Threshold“<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
p
MEDIAN-<strong>Filter</strong><br />
Original gefiltert<br />
• dünne Elemente werden eliminiert<br />
• dicke Elemente bleiben erhalten<br />
Ergebnis ist abhängig von lokaler Bildstruktur!<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Kann man <strong>Filter</strong> bauen, die auf<br />
bestimmte Strukturen<br />
gezielt reagieren?<br />
Bsp.: kleine Flecken eliminieren<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
IDEE:<br />
1. alle Regionen “schrumpfen” (kleine verschwinden)<br />
2. übrige Regionen wieder wachsen lassen<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
“Schrumpfen” einer Region (1)<br />
Original<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
“Schrumpfen” einer Region (2)<br />
Abschälen<br />
einer Schicht<br />
“Erosion”<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
“Schrumpfen” einer Region (3)<br />
Ergebnis<br />
“Erosion”<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Bsp. “Erosion”<br />
[Weeks]<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
“Wachsen” einer Region (1)<br />
Original<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
“Wachsen” einer Region (2)<br />
Hinzufügen<br />
einer Schicht<br />
“Dilation”<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
“Wachsen” einer Region (3)<br />
Ergebnis<br />
“Dilation”<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Bsp. “Dilation”<br />
[Weeks]<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
• EROSION, DILATION: wie definiert?<br />
• Gibt es weitere solche Operationen?<br />
<strong>Morphologische</strong> <strong>Filter</strong><br />
Strukturelement<br />
(analog zur<br />
Impulsantwort):<br />
H =<br />
•<br />
• • •<br />
•<br />
„Hot Spot“<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Darstellung als Mengen:<br />
Bsp.:<br />
0<br />
1<br />
2<br />
3<br />
-1<br />
0<br />
1<br />
0 1 2 3<br />
• • •<br />
-1 0 1<br />
• •<br />
A = {(1,1), (2,1), (2,2)}<br />
(Bild)<br />
H = {(0,0), (1,0)}<br />
(Strukturelement)<br />
(x a ,y a )<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Dilation<br />
B = A ⊕ H<br />
Ergebnis Bild Strukturelement<br />
A ⊕ H = {(x,y) | (x,y) = (x a ,y a ) + (x h ,y h )}<br />
für alle möglichen Paare<br />
(x a ,y a ) ∈ A, (x h ,y h ) ∈ H<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
0<br />
1<br />
2<br />
3<br />
Bsp.: Dilation<br />
0 1 2 3<br />
• • •<br />
A ⊕ H<br />
-1 0 1<br />
⊕<br />
-1<br />
0<br />
=<br />
1<br />
• •<br />
A ⊕ H = {<br />
(1,1)+(0,0) , (1,1)+(1,0) ,<br />
(2,1)+(0,0) , (2,1)+(1,0) ,<br />
(2,2)+(0,0) , (2,2)+(1,0) }<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
0<br />
1<br />
2<br />
3<br />
0 1 2 3<br />
• • •<br />
•<br />
•
0<br />
1<br />
2<br />
3<br />
Bsp.: Dilation<br />
0 1 2 3<br />
• • •<br />
A ⊕ H = B<br />
-1 0 1<br />
⊕<br />
-1<br />
0<br />
=<br />
1<br />
• •<br />
Wenn am Hotspot das Element<br />
in A gesetzt ist, dann werden in B<br />
alle Elemente von H markiert!<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
0<br />
1<br />
2<br />
3<br />
0 1 2 3<br />
• • •<br />
•<br />
•
0<br />
1<br />
2<br />
3<br />
Erosion<br />
0 1 2 3<br />
• • •<br />
A H = B<br />
-1<br />
0<br />
1<br />
-1 0 1<br />
• •<br />
Wenn alle Elemente von H<br />
in A Überlappung finden, nur dann<br />
wird in B der Hot-Spot markiert!<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
=<br />
0<br />
1<br />
2<br />
3<br />
0 1 2 3<br />
•
Erosion<br />
Ergebnis<br />
B = A H<br />
Bild Strukturelement<br />
A H= {(x,y) | (x,y) + (x h ,y h ) ∈ A}<br />
für alle möglichen Punkte<br />
(x h ,y h ) ∈ H<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
0<br />
1<br />
2<br />
3<br />
Bsp.: Erosion<br />
0 1 2 3<br />
• • •<br />
A H= {(1,1)}<br />
A H = B<br />
-1<br />
0<br />
1<br />
-1 0 1<br />
• •<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
=<br />
0<br />
1<br />
2<br />
3<br />
0 1 2 3<br />
•<br />
weil (1,1)+(0,0)=(1,1) ∈A<br />
und (1,1)+(1,0)=(2,1) ∈A
Dilation<br />
Erosion<br />
A ⊕ H = H ⊕ A kommutativ<br />
(A ⊕ B) ⊕ C = A ⊕ (B ⊕ C) Kettenregel<br />
(assoziativ)<br />
A H = H A<br />
NICHT kommutativ<br />
(A B) C) = A (B ⊕ C) Kettenregel<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Typische Strukturelemente<br />
•<br />
•<br />
• •<br />
•<br />
4-Neighborhood<br />
• • •<br />
• • •<br />
• • •<br />
8-Neighborhood<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
•<br />
•<br />
•<br />
• • •<br />
•<br />
•<br />
•<br />
•<br />
• •<br />
•<br />
•<br />
• •<br />
• •<br />
Small Disk<br />
•<br />
•<br />
•
Paarweise Anwendung (Erosion + Dilation)<br />
Erosion<br />
(H)<br />
“Opening”<br />
(A H) ⊕ H<br />
Dilation<br />
(H)<br />
= A H<br />
Kleine Strukturen werden eliminiert<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Bsp. “Opening”<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Bsp. “Opening”<br />
[Weeks]<br />
(A H) ⊕ H<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Paarweise Anwendung (Dilation + Erosion)<br />
Dilation<br />
(H)<br />
(A ⊕ H) H<br />
“Closing”<br />
Erosion<br />
(H)<br />
= A H<br />
Kleine Abstände werden geschlossen<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Bsp. “Closing”<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Bsp. “Closing”<br />
[Weeks]<br />
(A ⊕ H) H<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
<strong>Morphologische</strong>r <strong>Filter</strong><br />
• Mengen aller<br />
Translationen<br />
• Vereinigung aus<br />
Objekt und<br />
Strukturelement ist<br />
nicht leer<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Segmentierung: Threshold<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Segmentierung:<br />
Erosion<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Segmentierung:<br />
Connect<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Segmentie-rung:<br />
Dilation<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Binäre <strong>Morphologische</strong> <strong>Filter</strong><br />
(Beispiele)<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Binäre Dilation<br />
(Disk r=2.5, 5.0, 20)<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Binäre Erosion<br />
(Disk r=1.0, 2.5, 5.0)<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Binäres Opening<br />
= Erosion + Dilation<br />
(Disk r=1.0, 2.5, 5.0)<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Binäres Closing<br />
= Dilation + Erosion<br />
(Disk r=1.0, 2.5, 5.0)<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
<strong>Morphologische</strong> Konturextrakion<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Outline<br />
Schnittmenge<br />
Erosion<br />
H: 4-neighborhood H: 8-neighborhood<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
4-neighborhood filter -><br />
8-neighborhood contour<br />
Outline<br />
8-neighborhood filter -><br />
4-neighborhood contour<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Skeletonize<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
"Intelligent"<br />
Morphology
Zusammenfassung:<br />
• <strong>Morphologische</strong> <strong>Filter</strong> (MF) reagieren selektiv<br />
auf Bildstrukturen<br />
• Erosion/Dilation sind grundlegende Operationen<br />
• Opening/Closing = Hintereinanderausführen<br />
von Erosion/Dilation<br />
• Eindeutige Spezifikation durch<br />
“Strukturelement”<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Graustufen Morphologie<br />
Morphologie II<br />
W. Burger, W. Backfrieder<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Grayscale Morphology<br />
GS Dilation:<br />
( I<br />
GS Erosion:<br />
( I<br />
⊕ H ) ( u,<br />
v)<br />
H ) ( u,<br />
v)<br />
=<br />
=<br />
I( u,<br />
v),<br />
H ( i,<br />
j)<br />
∈<br />
max<br />
( i,<br />
j)<br />
∈H<br />
min<br />
( i,<br />
j)<br />
∈H<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
R<br />
{ I(<br />
u + i,<br />
v + j)<br />
+ H ( i,<br />
j)<br />
}<br />
{ I(<br />
u + i,<br />
v + j)<br />
− H ( i,<br />
j)<br />
}
Grayscale Dilation<br />
6 7 3 4<br />
5 6 6 8<br />
6 4 5 2<br />
6 4 2 3<br />
A +<br />
A<br />
H<br />
7 8 4<br />
6 8 7<br />
7 5 6<br />
⊕<br />
max<br />
H<br />
1 1 1<br />
1 2 1<br />
1 1 1<br />
8 9<br />
7 9<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
=
Grayscale Dilation<br />
(Disk)<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
r=2.5<br />
r=10
Grayscale Erosion<br />
6 7 3 4<br />
5 6 6 8<br />
6 4 5 2<br />
6 4 2 3<br />
A −<br />
A<br />
H<br />
5 6 2<br />
4 4 6<br />
5 3 4<br />
min<br />
H<br />
1 1 1<br />
1 2 1<br />
1 1 1<br />
2 1<br />
1 1<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
=
Grayscale Erosion<br />
(Disk)<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
r=2.5<br />
r=10
Grayscale Opening<br />
(Disk)<br />
(A H) ⊕ H<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
r=2.5<br />
r=10
Grayscale Closing<br />
(Disk)<br />
(A ⊕ H) H<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder<br />
r=2.5<br />
r=10
Geometrische Interpretation<br />
Opening<br />
Closing<br />
[Parker]<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Color Morphology • kein spezieller Ansatz<br />
• separate Anwendung auf R/G/B<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Erosion<br />
Dilation<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Opening<br />
Closing<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder
Zusammenfassung:<br />
• <strong>Morphologische</strong> <strong>Filter</strong> (MF) reagieren selektiv auf<br />
Bildstrukturen<br />
• Erosion/Dilation sind grundlegende Operationen<br />
• Opening/Closing = Hintereinanderausführen von<br />
Erosion/Dilation<br />
• Eindeutige Spezifikation durch “Strukturelement”<br />
• MF werden vorwiegend auf Binärbilder angewandt<br />
(auch für Graubilder definiert)<br />
Digitale Signal & Bildverarbeitung Werner Backfrieder