13.07.2015 Aufrufe

Bildverarbeitung 8 – Matching und Segmentierung - IGP - ETH Zürich

Bildverarbeitung 8 – Matching und Segmentierung - IGP - ETH Zürich

Bildverarbeitung 8 – Matching und Segmentierung - IGP - ETH Zürich

MEHR ANZEIGEN
WENIGER ANZEIGEN

Erfolgreiche ePaper selbst erstellen

Machen Sie aus Ihren PDF Publikationen ein blätterbares Flipbook mit unserer einzigartigen Google optimierten e-Paper Software.

<strong>Bildverarbeitung</strong>8 <strong>–</strong> <strong>Matching</strong> <strong>und</strong> <strong>Segmentierung</strong>Frühjahrssemester 2013 | <strong>ETH</strong> <strong>Zürich</strong>Jan Dirk Wegner, Stefan Walk, Michal Havlena <strong>und</strong> Konrad Schindler<strong>ETH</strong> <strong>Zürich</strong>, Institut für Geodäsie <strong>und</strong> Photogrammetrie


Themen heute§ <strong>Matching</strong>o Messen von Ähnlichkeit <strong>und</strong> Least Squares <strong>Matching</strong>o <strong>Matching</strong> von Punktwolken (Gr<strong>und</strong>züge)o Globale Ähnlichkeit: GIST-Merkmal§ <strong>Segmentierung</strong>o Region growing, Watershed, Split-and-Mergeo k-means clusteringo mean-shift<strong>Matching</strong> & <strong>Segmentierung</strong> 2


Template matchingWo ist Waldo?[©Kristen Grauman]Korrelation<strong>Matching</strong> & <strong>Segmentierung</strong>3


Template matchingWo ist Waldo?Template[©Kristen Grauman]Korrelation<strong>Matching</strong> & <strong>Segmentierung</strong>4


ZuordnungSuche nach homologen Merkmalen§Verknüpfungspunkte für dieautomatische Aerotriangulation§homologe Kanten für dieLinienauswertung§Dense matching: Generierung vonTiefenkarten z.B. zur Erstellung vonHöhenmodellen è für jedes Pixel imersten Bild wird das korrespondierendePixel im zweiten Bild gesucht§Suche nach bekannten Merkmalen(z.B. signalisierte Punkte) (siehe auchletzte Vorlesung zum Thema templatematching für Punktextration)<strong>Matching</strong> & <strong>Segmentierung</strong> 5


ZuordnungAufgabe: Merkmal mit grösster Ähnlichkeit zu einem Referenzmerkmal findenEinzelner Grauwert reicht nicht aus, Beschreibung eines Merkmals immer durcheine (M x N) Nachbarschaft è innerhalb eines Fensters, das über das Bild geschoben wird<strong>Matching</strong> & <strong>Segmentierung</strong> 6


ZuordnungAufgabe: Merkmal mit grösster Ähnlichkeit zu einem Referenzmerkmal findenEinzelner Grauwert reicht nicht aus, Beschreibung eines Merkmals immer durcheine (M x N) Nachbarschaft è innerhalb eines Fensters, das über das Bild geschoben wirdIm Luftbildfall: genäherterNormalfall, geringeProjektionsverzerrung zwischenbenachbarten Bildern<strong>Matching</strong> & <strong>Segmentierung</strong> 7


Themen heute§ <strong>Matching</strong>o Messen von Ähnlichkeit <strong>und</strong> Least Squares <strong>Matching</strong>o <strong>Matching</strong> von Punktwolken (Gr<strong>und</strong>züge)o Globale Ähnlichkeit: GIST-Merkmal§ <strong>Segmentierung</strong>o Region growing, Watershed, Split-and-Mergeo k-means clusteringo mean-shift<strong>Matching</strong> & <strong>Segmentierung</strong> 2


Summe der quadratischen Differenzen (SSD)SSD (sum of squared differences) zwischen Grauwerten einer NachbarschaftxxSSDyIJ( I(x,y)J ( x,y))= ∑∑ −x yI2yWertebereich [0..∞](0 → Bildregionen sind identisch)J<strong>Matching</strong> & <strong>Segmentierung</strong> 9


Summe der quadratischen Differenzen (SSD)SSD (sum of squared differences) zwischen Grauwerten einer NachbarschaftxxSSDyIJ( I(x,y)J ( x,y))= ∑∑ −x yI2yWertebereich [0..∞](0 → Bildregionen sind identisch)JVoraussetzungen:Bilder unterscheiden sich radiometrisch nicht <strong>und</strong>geometrisch nur durch eine Translation<strong>Matching</strong> & <strong>Segmentierung</strong> 10


KreuzkorrelationKorrelationskoeffizient (normalized cross-correlation, NCC) zwischen (geordneten,zentrierten) Intensitätswerten (vgl. Vorlesung zu Bildstatistik)xxyIyJWertebereich [-1..1](1/-1→ linear abhängig,0 → unkorreliert)<strong>Matching</strong> & <strong>Segmentierung</strong> 11


KreuzkorrelationKorrelationskoeffizient (normalized cross-correlation, NCC) zwischen (geordneten,zentrierten) Intensitätswerten (vgl. Vorlesung zu Bildstatistik)xxyIyJWertebereich [-1..1](1/-1→ linear abhängig,0 → unkorreliert)Voraussetzungen:Bilder unterscheiden sich geometrisch durch eine Translation <strong>und</strong> radiometrischdurch eine konstante Helligkeitsänderung <strong>und</strong> eine lineare Kontrastskalierung<strong>Matching</strong> & <strong>Segmentierung</strong> 12


Kreuzkorrelation: Beispiele1.0<strong>Matching</strong> & <strong>Segmentierung</strong> 13


Kreuzkorrelation: Beispiele1.0 0.7<strong>Matching</strong> & <strong>Segmentierung</strong> 14


Brendan Allen 15Monika’s group beats us back to the lodge by a long shot, and they’ve already assembleda line for tuna sandwiches, crackers, and fresh slices of mango. There’s not many slices left, so Igrab the pit and use my teeth to scrape off some of the bits of fruit and pulp still attached,dripping juices aro<strong>und</strong> my mouth and into my beard. I keep gnawing. I’ve already learned not tolet good Costa Rican fruit go to waste here <strong>–</strong> their sweetness intoxicates in a way I’ve neverexperienced at home.After lunch, the entire group hikes to set up one last FIT, this time only a short walk fromthe lodge. While Andrew, Devin, Luke and Michael work on tying the strings, the rest of us staybusy by collecting leaf litter from the gro<strong>und</strong> <strong>–</strong> bits of dirt and vegetation crawling with insects.Unlike last time, this one goes off without a hitch. We return to the lodge with our bags of leaflitter, which we use for another kind of trap called a “Berlese funnel”. Again, like the light trapand the FIT, the Berlese funnel’s ingenuity works hand in hand with its simplicity: First, youdump the leaf litter (filled with all kinds of crawling insects) into cloth funnels hanging from asuspended cord. The leaf litter rests on top of a wide metal screen, <strong>und</strong>erneath which the funnelcontinues into an attached baggie of alcohol. You then complete the trap by sealing a light bulbin the top of the funnel, creating something that looks a bit like a Chinese paper lantern <strong>–</strong> onethat happens to be filled with bugs and dirt. As the light heats and dries the leaf litter, the insectscrawl lower and lower into what little moist soil remains, eventually dropping through the screenand into the alcohol baggie. That’s it <strong>–</strong> after a day or two, a thick mass of preserved insectsdangles, waiting for collection.* * *Over the next few days, we fall into a regular routine: breakfast, a morning session ofchecking and replacing traps, lunch, an afternoon session of hand-collecting insects with nets and


Kreuzkorrelation: Beispiele1.0 0.7 0.1 -1.0<strong>Matching</strong> & <strong>Segmentierung</strong> 16


Least-Squares-<strong>Matching</strong>Suche nach der besten Zuordnung über z.B. SSD oder NCC:§ Ausprobieren aller Möglichkeiten in einer Nachbarschaft§ Genauigkeit ≈1 pixel<strong>Matching</strong> & <strong>Segmentierung</strong> 17


Least-Squares-<strong>Matching</strong>Suche nach der besten Zuordnung über z.B. SSD oder NCC:§ Ausprobieren aller Möglichkeiten in einer Nachbarschaft§ Genauigkeit ≈1 pixel§ Nur Translation berücksichtigtGenauigkeitssteigerung§ Annahme: Pixelwerte sind verrauschte Samples desselbenkontinuierlichen Signals§ Gesucht: Transformation für die die Summe derFehlerquadrate minimal wird§ Verschiedene Transformationsmodelle möglichProfile des zu matchendenTemplates verschiedenerIterationenGrauwertprofil im Bild<strong>Matching</strong> & <strong>Segmentierung</strong> 18


Subpixelgenaue VerschiebungModell: Bild J = Bild I mit Verschiebung + Rauschenu,v klein (Grobregistrierung bereits erfolgt) <strong>–</strong> Taylor-Entwicklung möglichElimination der Rauscheinflüsse <strong>und</strong> robuste Bestimmung der Parameter u,v durch lineareAusgleichsrechnungBestimmung von u,v durch Lösen der Normalengleichung<strong>Matching</strong> & <strong>Segmentierung</strong> 19


Subpixelgenaue VerschiebungGegebenenfalls ist es sinnvoll, die Pixel zu gewichten durch (diagonale) Gewichtsmatrix P<strong>Matching</strong> & <strong>Segmentierung</strong> 20


Subpixelgenaue VerschiebungGegebenenfalls ist es sinnvoll, die Pixel zu gewichten durch (diagonale) Gewichtsmatrix P<strong>Matching</strong> & <strong>Segmentierung</strong> 21


Subpixelgenaue VerschiebungGegebenenfalls ist es sinnvoll, die Pixel zu gewichten durch (diagonale) Gewichtsmatrix PWahl der Elemente w(i,j) der Gewichtsmatrix P:1 im Fenster, 0 ausserhalbGauss-Gewichte<strong>Matching</strong> & <strong>Segmentierung</strong> 22


Least-Squares-<strong>Matching</strong>Annahme bisher:Grauwertstruktur ist bis auf Translation gleichè In Wirklicheit: unterschiedliche Projektionsverzerrung <strong>und</strong> RadiometrieGeneralisierung§ allgemeinere Transformation zwischen den beiden Bildern§ Auch hier: Bestimmung der Transformationsparameter nach der Methode derkleinsten Quadrate<strong>Matching</strong> & <strong>Segmentierung</strong> 23


Least-Squares-<strong>Matching</strong>“Zuordnung nach kleinsten Quadraten”§ prinzipiell mit jeder parametrischen Transformation vonGeometrie <strong>und</strong> Bildintensität möglich<strong>Matching</strong> & <strong>Segmentierung</strong> 24


Least-Squares-<strong>Matching</strong>“Zuordnung nach kleinsten Quadraten”§ prinzipiell mit jeder parametrischen Transformation vonGeometrie <strong>und</strong> Bildintensität möglichHäufigstes Modell§ Geometrie: (eingeschränkte) Affintransformation als lokaleNäherung für die projektive Abbildung§Radiometrie: lineare Transformation (wie NCC)<strong>Matching</strong> & <strong>Segmentierung</strong> 25


Least-squares <strong>Matching</strong>Geometrie:TranslationRotation, Skalierung <strong>und</strong> Scherung<strong>Matching</strong> & <strong>Segmentierung</strong> 26


Least-squares <strong>Matching</strong>Geometrie:TranslationRotation, Skalierung <strong>und</strong> ScherungRadiometrie:Offsetlineare Skalierung<strong>Matching</strong> & <strong>Segmentierung</strong> 27


Least-squares <strong>Matching</strong>Geometrie:TranslationRotation, Skalierung <strong>und</strong> ScherungRadiometrie:Offsetlineare SkalierungM x N Beobachtungsgleichungen:Näherungswerte:<strong>Matching</strong> & <strong>Segmentierung</strong> 28


LSQ-<strong>Matching</strong>: BeispielAusgangstemplate1. IterationReferenzbild2. Iteration 3. Iteration<strong>Matching</strong> & <strong>Segmentierung</strong> 29


LSQ-<strong>Matching</strong>: Eigenschaften§ Die zwei radiometrischen Korrekturen der patch-Grauwerte werden nach jederIteration angebracht§ In jeder Iteration werden die patch-Grauwerte nach der Affin-Transformation durchbilineare Interpolation neu bestimmt<strong>Matching</strong> & <strong>Segmentierung</strong> 30


LSQ-<strong>Matching</strong>: Eigenschaften§ Die zwei radiometrischen Korrekturen der patch-Grauwerte werden nach jederIteration angebracht§ In jeder Iteration werden die patch-Grauwerte nach der Affin-Transformation durchbilineare Interpolation neu bestimmt§ Iterative Ausgleichung wird gestoppt, wenn die Änderungen der zu schätzendenParameter unterhalb zuvor definierter Schwellwerte liegen oder eine maximaleAnzahl an Iterationen erreicht ist<strong>Matching</strong> & <strong>Segmentierung</strong> 31


LSQ-<strong>Matching</strong>: Eigenschaften§ Die zwei radiometrischen Korrekturen der patch-Grauwerte werden nach jederIteration angebracht§ In jeder Iteration werden die patch-Grauwerte nach der Affin-Transformation durchbilineare Interpolation neu bestimmt§ Iterative Ausgleichung wird gestoppt, wenn die Änderungen der zu schätzendenParameter unterhalb zuvor definierter Schwellwerte liegen oder eine maximaleAnzahl an Iterationen erreicht ist§ Standardabweichungen aller Parameter werden (häufig zu optimistisch)mitgeschätzt:Residuen (Differenz zwischen transformiertem Bild u. Template)§ Genauigkeiten: 0,01 - 0,3 Pixel für signalisierte Punkte0.5 - 1.5 Pixel für natürliche Oberflächenpunkte<strong>Matching</strong> & <strong>Segmentierung</strong> 32


LSQ-<strong>Matching</strong>Nachteile LSQ-matching:§ gute Näherungswerte für zu schätzende Parameter notwendig§ Probleme bei abrupten Höhenänderungen, z.B. an Gebäuderändernè LSQ-matching kann zu ungewollter Kantenglättung führenè für hoch aufgelöste Oberflächenmodelle in urbanen Gebieten kaum geeignet<strong>Matching</strong> & <strong>Segmentierung</strong> 33


LSQ-<strong>Matching</strong>Nachteile LSQ-matching:§ gute Näherungswerte für zu schätzende Parameter notwendig§ Probleme bei abrupten Höhenänderungen, z.B. an Gebäuderändernè LSQ-matching kann zu ungewollter Kantenglättung führenè für hoch aufgelöste Oberflächenmodelle in urbanen Gebieten kaum geeignetFür komplexe Szenen z.B. Semi-global matching (SGM) <strong>und</strong> Markov Random FieldsWichtiges Thema aktueller Forschung! mehr Details in Photogrammetric Computer Vision bei Konrad Schindler<strong>Matching</strong> & <strong>Segmentierung</strong> 34


Themen heute§ <strong>Matching</strong>o Messen von Ähnlichkeit <strong>und</strong> Least Squares <strong>Matching</strong>o <strong>Matching</strong> von Punktwolken (Gr<strong>und</strong>züge)o Globale Ähnlichkeit: GIST-Merkmal§ <strong>Segmentierung</strong>o Region growing, Watershed, Split-and-Mergeo k-means clusteringo mean-shift<strong>Matching</strong> & <strong>Segmentierung</strong> 2


<strong>Matching</strong> von PunktwolkenÜbergang von 2D-matching von Bildern zu 3D-matching von Punktwolken z.B.von Laserscannern (terrestrisch oder luftgestützt), aus dense matching vonBildern, aus range cameras etc.<strong>Matching</strong> & <strong>Segmentierung</strong>36


<strong>Matching</strong> von PunktwolkenÜbergang von 2D-matching von Bildern zu 3D-matching von Punktwolken z.B.von Laserscannern (terrestrisch oder luftgestützt), aus dense matching vonBildern, aus range cameras etc.Schwierigkeiten:§ Mit der Entfernung abnehmende Punktdichte bei Laserscannern aufgr<strong>und</strong>des polaren Messprinzips matching korrespondierender Bereicheaufgenommen mit unterschiedlicher Punktdichte sehr schwierig§ Sequentielle Abtastung von Scanprofilen führt bei sich bewegendenObjekten zu Verschiebungseffekten§ Sehr grosse Datenmenge etc.è Aktueller Forschungsaspekt am <strong>IGP</strong>, hier nur ein kurzer Einblick<strong>Matching</strong> & <strong>Segmentierung</strong>37


<strong>Matching</strong> von Punktwolken: Nahbereichsphotogrammetrie<strong>Matching</strong> von Punktwolkenaufgenommen mit SLRKameraNikon D80[David Novak]<strong>Matching</strong> & <strong>Segmentierung</strong>38


<strong>Matching</strong> von Punktwolken: Laserscanning[Pascal Theiler]Grobregistrierung mittelskorrespondierenderEbenenschnittpunkteFeinregistrierung mit Iterative Closest Point (ICP) Methode<strong>Matching</strong> & <strong>Segmentierung</strong>39


<strong>Matching</strong> von Punktwolken: Range images[Tobias Kohoutek]GrobregistrierungFeinregistrierungMESA SR4000Time-of-flight range Kamera<strong>Matching</strong> & <strong>Segmentierung</strong>40


Themen heute§ <strong>Matching</strong>o Messen von Ähnlichkeit <strong>und</strong> Least Squares <strong>Matching</strong>o <strong>Matching</strong> von Punktwolken (Gr<strong>und</strong>züge)o Globale Ähnlichkeit: GIST-Merkmal§ <strong>Segmentierung</strong>o Region growing, Watershed, Split-and-Mergeo k-means clusteringo mean-shift<strong>Matching</strong> & <strong>Segmentierung</strong> 2


Globale Ähnlichkeit von BildernStatistische Methoden wie SSD, NCC oder Parameter von Verteilungsfunktionen derGrauwerte geben nicht die räumliche Anordnung wieder.GIST-Merkmal [Antonio Torralba & Aude Oliva, 2001] zur Beschreibung ganzer Szenen[©Torralba & Oliva]<strong>Matching</strong> & <strong>Segmentierung</strong>42


Globale Ähnlichkeit von BildernStatistische Methoden wie SSD, NCC oder Parameter von Verteilungsfunktionen derGrauwerte geben nicht die räumliche Anordnung wieder.GIST-Merkmal [Antonio Torralba & Aude Oliva, 2001] zur Beschreibung ganzer SzenenIdee des spatial envelope:[©Torralba & Oliva]3D-Grauwertoberflächen unterschiedlicher Szenen haben verschiedene dominante räumlicheStrukturen è Erfassung im Frequenzraum + Information über räumliche Anordnung<strong>Matching</strong> & <strong>Segmentierung</strong>43


Globale Ähnlichkeit von Bildern: GISTGlobale Fourier-Spektren verschiedener Bilder natürlicher Szenen[©Torralba & Oliva]<strong>Matching</strong> & <strong>Segmentierung</strong>44


Globale Ähnlichkeit von Bildern: GISTGlobale Fourier-Spektren verschiedener Bilder natürlicher Szenen[©Torralba & Oliva]Fourier-Spektren (mit windowed Fourier transform) von 16 Bildregionen verschiedener Bilder<strong>Matching</strong> & <strong>Segmentierung</strong>45


Globale Ähnlichkeit von Bildern: GISTGlobale Fourier-Spektren verschiedener Bilder natürlicher Szenen[©Torralba & Oliva]Fourier-Spektren (mit windowed Fourier transform) von 16 Bildregionen verschiedener BilderGIST: Berechnung der Fourier-Spektren (global <strong>und</strong> lokal), dann Bestimmung der unkorreliertenInformation <strong>und</strong> Reduktion im Merkmalsraum mittels Principal Component Analysisè Bestimmung repräsentativer Merkmale mit Regression basierend auf Trainingsdaten<strong>Matching</strong> & <strong>Segmentierung</strong>46


Themen heute§ <strong>Matching</strong>o Messen von Ähnlichkeit <strong>und</strong> Least Squares <strong>Matching</strong>o <strong>Matching</strong> von Punktwolken (Gr<strong>und</strong>züge)o Globale Ähnlichkeit: GIST-Merkmal§ <strong>Segmentierung</strong>o Region growing, Watershed, Split-and-Mergeo k-means clusteringo mean-shift<strong>Matching</strong> & <strong>Segmentierung</strong> 2


Wozu <strong>Segmentierung</strong> von Regionen?Optical scene flow è Bewegungsdetektion <strong>und</strong> 3D-Szenenanalyse[©Christoph Vogel]<strong>Matching</strong> & <strong>Segmentierung</strong>48


Wozu <strong>Segmentierung</strong> von Regionen?Optical scene flow è Bewegungsdetektion <strong>und</strong> 3D-Szenenanalyse[©Christoph Vogel]Semantic segmentation è Automatische Kartierung <strong>und</strong> Objektdetektion [©J. Montoya & J.D.Wegner]<strong>Matching</strong> & <strong>Segmentierung</strong>49


<strong>Segmentierung</strong><strong>Segmentierung</strong> bezeichnet allgemein die Extraktion von Primitiven im Bildbzw. die Zusammenfassung von Pixeln aufgr<strong>und</strong> bestimmterZusammengehörigkeitskriterien. Primitive sind in der Regel Punkte, Linienoder Regionen. Hier: Betrachtung von Regionen<strong>Matching</strong> & <strong>Segmentierung</strong>50


<strong>Segmentierung</strong><strong>Segmentierung</strong> bezeichnet allgemein die Extraktion von Primitiven im Bildbzw. die Zusammenfassung von Pixeln aufgr<strong>und</strong> bestimmterZusammengehörigkeitskriterien. Primitive sind in der Regel Punkte, Linienoder Regionen. Hier: Betrachtung von Regionenguthomogene Texturmit shape prior KonturschlechtEigenschaften von Regionen:§ Einheitlichkeit (Homogenität der Grauwerte, Farbe, Textur,...)§§Nachbarschaft (Verb<strong>und</strong>enheit)Regionen sind von gemeinsamer Kontur begrenzt[©Bodo Rosenhahn]Die Kriterien für eine „gute“ <strong>Segmentierung</strong> hängen immer von der Anwendung ab!<strong>Matching</strong> & <strong>Segmentierung</strong>51


Global Thresholding (Schwellwertverfahren)<strong>Segmentierung</strong> basierend auf Bildhistogramm des ganzen BildesUnterteilung in Klassen mittels verschiedener Schwellwerte (thresholds)Zuordnung der Pixel zu Klasse, wenn Wert innerhalb des zugeordneten BereichsSchwellwertHistogramm Glättungsfilter sinnvoll als Vorverarbeitungsschritt, um Rauschen zu minimieren!


Adaptive ThresholdingEingangsbild Globales Histogramm Globaler Schwellwert


Adaptive ThresholdingEingangsbild Globales Histogramm Globaler SchwellwertLokale AnpassungLokale Histogramme pro Bildregion Adaptive Thresholding


Thresholding-MethodenVorteile:einfache (sehr schnelle) Berechnung <strong>und</strong> gute <strong>Segmentierung</strong>sergebnisse beihomogenen Objektenè z.B. gut geeignet bei simpler Vordergr<strong>und</strong>-/Hintergr<strong>und</strong>trennung


Thresholding-MethodenVorteile:einfache (sehr schnelle) Berechnung <strong>und</strong> gute <strong>Segmentierung</strong>sergebnisse beihomogenen Objektenè z.B. gut geeignet bei simpler Vordergr<strong>und</strong>-/Hintergr<strong>und</strong>trennungNachteile:Ø Erfordert evtl. Anpassung der Schwellwerte für neue BilderØ Wahl der Thresholds <strong>und</strong> KlassenanzahlØ Unvollständige <strong>Segmentierung</strong>: Löcher in den RegionenØ Sehr unterschiedliche Größen <strong>und</strong> Formen von Regionen


Themen heute§ <strong>Matching</strong>o Messen von Ähnlichkeit <strong>und</strong> Least Squares <strong>Matching</strong>o <strong>Matching</strong> von Punktwolken (Gr<strong>und</strong>züge)o Globale Ähnlichkeit: GIST-Merkmal§ <strong>Segmentierung</strong>o Region growing, Watershed, Split-and-Mergeo k-means clusteringo mean-shift<strong>Matching</strong> & <strong>Segmentierung</strong> 2


Watershed-<strong>Segmentierung</strong>Annahme:Regionen getrenntdurch hohe Gradienten1. Gradientenbetrag 2. „Grauwertgebirge“<strong>Matching</strong> & <strong>Segmentierung</strong>58


Watershed-<strong>Segmentierung</strong>Annahme:Regionen getrenntdurch hohe Gradienten1. Gradientenbetrag 2. „Grauwertgebirge“3. Fluten der Oberflächebeginnend bei den Minima<strong>Matching</strong> & <strong>Segmentierung</strong>59


Watershed-<strong>Segmentierung</strong>Annahme:Regionen getrenntdurch hohe Gradienten1. Gradientenbetrag 2. „Grauwertgebirge“3. Fluten der Oberflächebeginnend bei den Minima4. Unterbinden der Verschmelzung des Wassers, das ausunterschiedlichen Quellen stammt<strong>Matching</strong> & <strong>Segmentierung</strong>60


Watershed-<strong>Segmentierung</strong>Annahme:Regionen getrenntdurch hohe Gradienten1. Gradientenbetrag 2. „Grauwertgebirge“3. Fluten der Oberflächebeginnend bei den Minima5. <strong>Segmentierung</strong> è„Wasserscheiden“ (rot) sindKonturen der Regionen4. Unterbinden der Verschmelzung des Wassers, das ausunterschiedlichen Quellen stammt<strong>Matching</strong> & <strong>Segmentierung</strong>61


Flutungsalgorithmus§Jedes neu überflutete Pixel ist in Isolation:Es ist nicht zu anderen überfluteten Pixeln der Höhen h


Flutungsalgorithmus§Jedes neu überflutete Pixel ist in Isolation:Es ist nicht zu anderen überfluteten Pixeln der Höhen h


Flutungsalgorithmus§Jedes neu überflutete Pixel ist in Isolation:Es ist nicht zu anderen überfluteten Pixeln der Höhen h


Region GrowingAblauf:1. Selektion der Seed Points (z.B. zufällig, in Raster etc.)2. Festlegung des Homogenitätskriteriums (Farbe, Textur, ...)3. Wachstum (Growing) è Analyse benachbarter Pixel <strong>und</strong> Integration, wenn Distanzentsprechend Homogenitätskriterium kleiner als Schwellwert ist4. Prozedur beendet, wenn keine weiteren Pixel zu Regionen mehr hinzugefügt werden können<strong>Matching</strong> & <strong>Segmentierung</strong>65


Region GrowingNachteile:§Seed Points müssen festgelegt werden§Je nach Wahl des Schwellwerts <strong>und</strong> derVerteilung der Seed Points können Löcherentstehen (Pixel die zu keiner Regiongehören)§Ineinander Laufen benachbarter Regionen§Regionen sehr unterschiedlicher Größe <strong>und</strong> Form<strong>Matching</strong> & <strong>Segmentierung</strong>66


Split & MergeEingangsbild1. Split per Quadtree (vgl.Kompression)Festlegung einesHomogenitätskriteriums èAufteilen jedes Quadranten, fallsnicht homogen2. MergeVergleich benachbarterRegionen <strong>und</strong> Verschmelzung,falls homogenSplitMergeErgebnis<strong>Matching</strong> & <strong>Segmentierung</strong>68


Split & MergeEingangsbild1. Split per Quadtree (vgl.Kompression)Festlegung einesHomogenitätskriteriums èAufteilen jedes Quadranten, fallsnicht homogen2. MergeSplitMergeErgebnisVergleich benachbarterRegionen <strong>und</strong> Verschmelzung,falls homogenVorteil gegenüber Region Growing: keine LöcherNachteil: Ergebnis abhängig von Merge-Reihenfolge<strong>Matching</strong> & <strong>Segmentierung</strong>69


Themen heute§ <strong>Matching</strong>o Messen von Ähnlichkeit <strong>und</strong> Least Squares <strong>Matching</strong>o <strong>Matching</strong> von Punktwolken (Gr<strong>und</strong>züge)o Globale Ähnlichkeit: GIST-Merkmal§ <strong>Segmentierung</strong>o Region growing, Watershed, Split-and-Mergeo k-means clusteringo mean-shift<strong>Matching</strong> & <strong>Segmentierung</strong> 2


K-Mean ClusteringK-Means ist eine unüberwachte Klassifizierungmethode <strong>und</strong> sucht nach KClusterzentren im Merkmalsraum.Ablauf (für normalverteilte Daten):1. Initialisiere die K Clusterzentren z.B. auf Basis des Zufallsprinzips oderregelmäßig entlang der Hauptdiagonale im Merkmalsraum2. Ordne jeden Datenpunkt dem nächsten Cluster zu3. Bestimme µ n <strong>und</strong> Σ k für jedes Cluster4. Wiederhole 2-3 bis eine stabile Situation eintritt <strong>und</strong> sich dieClusterzentren nicht mehr ändern<strong>Matching</strong> & <strong>Segmentierung</strong>71


K-Mean ClusteringK-Means ist eine unüberwachte Klassifizierungmethode <strong>und</strong> sucht nach KClusterzentren im Merkmalsraum.Ablauf (für normalverteilte Daten):1. Initialisiere die K Clusterzentren z.B. auf Basis des Zufallsprinzips oderregelmäßig entlang der Hauptdiagonale im Merkmalsraum2. Ordne jeden Datenpunkt dem nächsten Cluster zu3. Bestimme µ n <strong>und</strong> Σ k für jedes Cluster4. Wiederhole 2-3 bis eine stabile Situation eintritt <strong>und</strong> sich dieClusterzentren nicht mehr ändernVorteile: schnelle Berechnung, keine Löcher, auf beliebige Merkmale anwendbarNachteile: Anzahl K der Clusterzentren muss vorgegeben werden, sehr unterschiedlicheGrößen <strong>und</strong> Formen der Regionen, Pixel eines Clusters im Merkmalsraum müssen imräumlich nicht zwingend benachbart sein è „Ausfransen“ von Regionen<strong>Matching</strong> & <strong>Segmentierung</strong>72


K-means Clustering[©Bishop, 2006]<strong>Matching</strong> & <strong>Segmentierung</strong>73


Themen heute§ <strong>Matching</strong>o Messen von Ähnlichkeit <strong>und</strong> Least Squares <strong>Matching</strong>o <strong>Matching</strong> von Punktwolken (Gr<strong>und</strong>züge)o Globale Ähnlichkeit: GIST-Merkmal§ <strong>Segmentierung</strong>o Region growing, Watershed, Split-and-Mergeo k-means clusteringo mean-shift<strong>Matching</strong> & <strong>Segmentierung</strong> 2


Mean shiftMean shift ist eine Methode zur Bestimmung von Clustern im Merkmalsraum basierendauf der Merkmalsdichte [Comaniciu & Meer, 2002]Pixel werden dafür mit Merkmalen beschrieben (z.B. Farbe, Textur in Nachbarschaft des Pixels,Räumliche Position im Bild), die dann im Merkmalsraum abgebildet werdenè eine Koordinatenachse pro Merkmal<strong>Matching</strong> & <strong>Segmentierung</strong>75


Mean shiftMean shift schätzt lokale Dichten der Merkmale im Merkmalsraumè Beginnend bei zufällig ausgewählten lokalen Startpunkten werden entlang von Suchpfaden(schwarz) lokale Dichtemaxima innerhalb einer Nachbarschaft (kernel width) gesuchtIm gleichen Punkt (rot) endende Suchpfade werden dem gleichen Cluster zugeordnetVorteil: schnelles Verfahren <strong>und</strong> die Anzahl der Cluster muss nicht vorgegeben werdenAber: Anzahl der Cluster bzw. Segmente im Bild hängt ab von der kernel width<strong>Matching</strong> & <strong>Segmentierung</strong>76


<strong>Segmentierung</strong> komplexer SzenenGezeigte Verfahren oft nicht erfolgreich bei Bildern mit komplexer Szenenstrukturè Typischer Ansatz: erst Übersegmentierung in kleine aber regelmäßige Regionen <strong>und</strong>dann Kombination mit Klassifizierungsverfahren z.B. Maximum likelihood, NeuronaleNetze, Random Forests, Boosting, Support Vector Machines für Semantische<strong>Segmentierung</strong>


<strong>Segmentierung</strong> komplexer SzenenGezeigte Verfahren oft nicht erfolgreich bei Bildern mit komplexer Szenenstrukturè Typischer Ansatz: erst Übersegmentierung in kleine aber regelmäßige Regionen <strong>und</strong>dann Kombination mit Klassifizierungsverfahren z.B. Maximum likelihood, NeuronaleNetze, Random Forests, Boosting, Support Vector Machines für Semantische<strong>Segmentierung</strong>§Hilfreich bei komplexen Szenen:Lokale Nachbarschaften <strong>und</strong> globale Zusammenhänge berücksichtigen wie z.B. „Himmelist über dem Haus“, „Autos sind auf der Strasse“, „Türen befinden sich im unteren Teilvon Häuserfassaden“ è z.B. mit graph-basierten Methoden wie Conditional RandomFields (Erforschung in der Gruppe „Photogrammetrie & Fernerk<strong>und</strong>ung“)è in Vorlesung Image Interpretation bei Konrad Schindler

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!