Die (digitale) Bildverarbeitung nutzt die Mittel der Signalverarbeitung zur Aufbereitung und Speicherung von visuellen Informationen. Im Gegensatz zur Bildbearbeitung, welche sich mit der Manipulation von Bildern zur anschließenden Darstellung beschäftigt, dient die Bildverarbeitung als eine Zwischenstufe zu einer weitergehenden maschinellen Bearbeitung (Bildverstehen, Mustererkennung).

Objekte der Bildverarbeitung

---

Die in der Bildverarbeitung manipulierten Objekte lassen sich nach verschiedenen Gesichtspunkten klassifizieren:

Typ

• Reflexionsbilder (z. B. Kameraaufnahmen, Radaraufnahmen, Multispektralaufnahmen)
• Projektionsbilder (z. B. Röntgenaufnahmen, Ultraschall, Elektrophoretogramme)
• Schematisierte Bilder (z. B. Karten, Pläne, Dokumente)

Manchmal findet sich auch die Unterscheidung zwischen Reflexions- (s. o.), Absorptions- (Röntgen), und Emissionsbildern (astronomische Aufnahmen, Nuklearaufnahmen) in der Fachliteratur.

Codierung

• Rasterbilder (z. B. konventionell auf quadratischem Gitter)
• Kettencodierung (durch Richtungsangaben beschriebene Bildlinien)
• Transformationscodierung (Ergebnis einer Reihenentwicklung)
• Fraktale Codierung (Ausnutzung von Selbstähnlichkeit)

Operationen der Bildverarbeitung

---

Die Operationen in der Bildverarbeitung lassen sich in drei Klassen einteilen: Punktoperatoren, lokale Operatoren und globale Operatoren. Diese Einteilung basiert darauf, wie viele Bildpunkte bei dem jeweiligen Verfahren für die Berechnung eines neuen Grau- oder Farbwertes f(x,y) eines Bildpunktes (x,y) betrachtet werden. Punktoperatoren transformieren jeden Punkt eines Bildes einzeln nur in Abhängigkeit vom Grau- oder Farbwert und ggf. der Position im Bild. Beispiele hierfür sind Histogrammtransformationen wie die Histogrammspreizung oder die Histogrammäqualisation. Lokale Operatoren berechnen einen neuen Farb- oder Grauwert eines Bildpunktes immer auf Basis einer Nachbarschaft oder einer örtlich begrenzten Region um den Punkt. Hier seien als Beispiele Rangordnungsoperatoren oder morphologische Operatoren genannt. Globale Operatoren betrachten für die Transformation eines jeden Pixels immer das gesamte Bild, was beispielsweise bei der Fouriertransformation der Fall ist.

Punktoperatoren

Hauptartikel: Punktoperator (Bildverarbeitung)

Ein Punktoperator T ordnet einem Eingabebild f durch Transformation der Grauwerte der einzelnen Pixel ein Ergebnisbild f* zu. Der Grauwert f(x,y) eines Pixels (x,y) wird dabei nur in Abhängigkeit vom Grauwert selbst und eventuell von der Position des Pixels im Bild modifiziert:

$f^*(x,y)=T\_\{xy\}(f(x,y))$

Ist die Transformation von der Position des Pixels im Bild abhängig, so heißt sie inhomogen. Die Indizes x und y von T sollen diese Abhängigkeit verdeutlichen. In der Mehrheit der Fälle kommen jedoch homogene Transformationen zum Einsatz, bei denen diese Abhängigkeit nicht gegeben ist. Die Indizes werden dann überflüssig:

$f^*(x,y)=T(f(x,y))$

Histogramm

In der digitalen Bildverarbeitung versteht man unter einem Histogramm die statistische Häufigkeit der einzelnen Grau- bzw. Farbwerte in einem Bild. Das Histogramm eines Bildes erlaubt eine Aussage über die vorkommenden Grau- bzw. Farbwerte und über Kontrastumfang und Helligkeit des Bildes.

Typische Punktoperatoren bewirken eine sichtbare Änderung im Histogramm. So wird es bei einer Histogrammspreizung zwecks Erhöhung des Kontrasts auseinandergezogen oder bei einer Histogrammverschiebung zwecks Helligkeitskorrektur durch Addition einer Konstante verschoben. Eine klassische Anwendung von Histogrammen in der Bildverarbeitung liegt in der Äqualisation (auch Egalisierung oder Einebnung genannt), bei der eine Gleichverteilung der Grau- bzw. Farbwerte berechnet wird. Dadurch kann eine bessere Verteilung der Farbgebung erreicht werden, die über eine bloße Kontrastverstärkung hinausgeht.

Schwellwertoperationen

Hauptartikel: Schwellwertverfahren

Sollen Bilder in kleinere Farbräume konvertiert werden, ergibt sich das Problem der Bestimmung adäquater Schwellwerte. Eine klassische Aufgabe ist die Binarisierung eines Grauwertbildes. Typische Schwellwerte können aus dem Grauwerthistogramm eines Bildes gewonnen werden:

• Mittelmäßiger Grauwert
• Mittlerer oder erwarteter Grauwert
• Grauwertmedian
• Hauptsenken (zentralste relative Extremstellen des Grauwerthistogramms)
• Maximum-Entropie-Schwellwert

Andere Methoden bieten ein „Grauwertmischverteilungsmodell“ oder der Intermeans-Algorithmus.

Lokale Operatoren

Filter

Die Operationen, die ein Eingangsbild mit Hilfe einer mathematischen Abbildung in ein Ausgabebild überführen, heißen Filter. Die folgenden Filter beziehen sich im Wesentlichen auf Rasterbilder.

Ziel einer Filterung ist im Allgemeinen eine Verbesserung eines Musters. Im Einzelnen kann dies eine Reduktion störender Anteile, eine Hervorhebung informativer Anteile bzw. die Restaurierung eines idealen Musters sein. Die Abbildung welche ein gegebenes Bild auf ein Ausgabebild abbildet, heißt Transformation.

Die Funktionsfamilie $g\_\backslash mu$ mit $g\_\backslash mu\; =\; \backslash Im\backslash \{\backslash delta\_\backslash mu\backslash \}\; (\backslash mu\; \backslash in\; \backslash mathbb\{R\})$ heißt Impulsantwort des linearen Systems $\backslash Im$. $\backslash delta\_\backslash mu:\backslash mathbb\{R\}\backslash to\backslash mathbb\{R\}$ ist die Diracsche Deltafunktion mit Impuls bei $x=\backslash mu$.

Lineare Filter

Ein wichtiges Unterscheidungsmerkmal von Transformationen ist ihre mathematische Beschaffenheit. Eine Transformation $\backslash Im$ heißt Lineares System wenn für alle Funktionen $f,\; g:\; \backslash mathbb\{R\}\; \backslash to\; \backslash mathbb\{R\}$ und alle Skalare $a,b\; \backslash in\; \backslash mathbb\{R\}$ gilt:

$\backslash Im\; (af+bg)\; =\; a\backslash Im\; (f)\; +\; b\backslash Im\; (g)$

Lineare Systeme (=LSI-Systeme) und damit lineare Filter sind praktisch hinreichend und mathematisch handhabbar. Für gewöhnlich konstruiert man LSI-Systeme mit den zusätzlichen Eigenschaften:

• Kausalität: Für die Impulsantwort $g=g\_j$ gilt für : $g\_j=0$ und damit $h\_j=\backslash sum\_\{k=0\}^\backslash infty\; f\_\{j-k\}g\_k$.
• Stabilität:

Während man symbolisch das Filtern einer Folge ) vor. Man berechnet dann h, indem f und g diskret gefaltet werden:

• 1-dimensional: $h\_k=\backslash sum\_\{j=-\backslash infty\}^\{+\backslash infty\}f\_jg\_\{k-j\}$
• 2-dimensional:$h\_\{kl\}=\backslash sum\_\{i=-\backslash infty\}^\{+\backslash infty\}\backslash sum\_\{j=-\backslash infty\}^\{+\backslash infty\}f\_\{ij\}g\_\{k-i,l-j\}$

= Klassifikation von linearen Filtern

= Lässt sich ein LSI-System $\backslash Im$ eines Filters $h=\backslash Im\backslash \{f\backslash \}$ mit der Impulsantwort $b$ durch

$h\_j=\backslash sum\_\{m=1\}^Mb\_mf\_\{j-m\}$

oder durch seine Frequenzantwort

$G(z)=B(z^\{-1\})=\backslash sum\_\{m=0\}^Mb\_mz^\{-m\}$

beschreiben, so nennt man den Filter FIR-Filter (finite impulse response). Jeder transformierte Wert hängt damit allein von der konstanten Impulsantwort und damit endlich vielen Nachbarn ab.

Es lässt sich die Operation $h=\backslash Im\backslash \{f\backslash \}$ eines LSI-Systems durch eine Differenzengleichung

$\backslash sum\_\{n=0\}^Na\_nh\_\{j-n\}=\backslash sum\_\{m=0\}^Nb\_mf\_\{j-m\}$

beschreiben. Man nennt einen Filter AR-Filter (autoregressive), wenn er sich durch

$h\_j=\backslash sigma\; f\_j-\backslash sum\_\{n=0\}^Na\_nh\_\{j-v\}$

darstellen lässt, oder die Frequenzantwort

$G(z)=\backslash frac\{\backslash sigma\}\{A(z^\{-1\})\}=\backslash sigma/\backslash sum\_\{n=0\}^Na\_nz^\{-m\}$

besitzt. Man nennt einen Filter ARMA-Filter (autoregressive moving average), wenn er sich durch $h\_j=\backslash sum\_\{m=0\}^Mb\_mf\_\{j-m\}-\backslash sum\_\{n=1\}^Na\_nh\_\{j-v\}$

darstellen lässt, oder die Frequenzantwort

$G(z)=\backslash frac\{B(z^\{-1\})\}\{A(z^\{-1\})\}=\backslash sum\_\{m=0\}^Mb\_mz^\{-m\}/\backslash sum\_\{n=0\}^Na\_nz^\{-n\}$

besitzt. Filter der Art AR oder ARMA mit einer rekursiven Berechnungsformel für $h$ nennt man IIR-Filter (infinite impulse response).

= Filtermatrix

= Um in der Praxis die Filter anwenden zu können nutzt man in der Regel eine m x n Filtermatrix. Sie stellt die Impulsantwort des 2-dimensionalen FIR - Filters dar. Jeder Eintrag dieser Matrix entspricht einem Gewicht, mit welchem die Nachbarpunkte des zu transfomierenden Punktes multipliziert werden. Ein Punkt g(x,y) im Ergebnisbild erhält also den Wert, der sich dadurch ergibt, dass man die Punkte seiner Nachbarschaft mit den Werten der Filtermatrix gewichtet aufaddiert und anschließend durch die Gesamtsumme des Betrags der Gewichte teilt. Beispiel: 3x 3 Filtermatrix mit dem

Ableitungsoperator

1	0	-1
1	0	-1
1	0	-1

dem Normierungsfaktor 1/6

und dem Bildausschnitt

0	0	0	0	0
0	1	2	3	0
0	1	1	1	0
0	0	0	0	0

führt zu dem Ergebnis

0	0	0	0	0
0	-3	-2	3	0
0	-3	-2	3	0
0	0	0	0	0

und jeder Pixelwert noch durch 6 geteilt. Man sieht, dass der Rand mit Nullen aufgefüllt wurde, es wäre auch möglich, den Rand ins Unendliche zu erweitern. In der Praxis würde in der Regel noch eine Betragsbildung durchgeführt.

= Separierbarkeit

= In der Praxis nutzt man gerne separierbare Filter. Separierbar ist ein Filter immer dann, wenn man die Impulsantwort $g\_\{ij\}$ folgendermaßen ausdrücken kann: $g\_\{i,j\}\; =\; g1\_\{i\}\; g2\_\{j\}$

Es gilt dann:$h\_\{kl\}=\backslash sum\_\{i=-\backslash infty\}^\{+\backslash infty\}\backslash sum\_\{j=-\backslash infty\}^\{+\backslash infty\}f\_\{ij\}g\_\{k-i,l-j\}\; =\; \backslash sum\_\{i\; =\; -\backslash infty\}^\{+\backslash infty\}g1\_\{k-i\}\; \backslash sum\_\{j=\; -\backslash infty\}^\{+\backslash infty\}\; f\_\{ij\}g2\_\{l-j\}$ Man kann dann die 2D-Faltung aufspalten in zwei 1D-Faltungen. Dies verringert den Aufwand einer Faltung bei einem Bildpunkt von M*N auf M+N Multiplikationen. Wobei M die Breite und N die Höhe der FIR - Filtermatrix ist.

Nichtlineare Filter

= Morphologische Operationen

= Während Grauwertrasterbilder für gewöhnlich als Abbildung von Koordinatenpaaren auf Farbwerte dargestellt werden, bietet sich für Binärbilder eine Darstellung als Menge der gesetzten Pixel an. Operationen können dann als Mengenoperation zwischen den Rasterpunktmengen des Binärbildes und eines Strukturelementes beschrieben werden. Aus den Basisoperationen der Morphologischen Bildverarbeitung, der Erosion und Dilatation, lassen sich die Operationen Öffnung, Schließung und schließlich die Morphologische Glättung definieren.

Filter nach Anwendungsgebieten

= Glättung

= Durch Glättung kann das Bildrauschen vermindert werden, grobere Strukturen bleiben dagegen erhalten. Hier ein Dead_tree_salt_and_pepper.png für ein (künstlich) verrauschtes Bild. Auf das Frequenzspektrum eines Bildes bezogen kommt eine Glättung einem Tiefpassfilter gleich. Typische FIR-Glättungsfilter sind

• Mittelwertfilter (box filter): Ein Mittelwertfilter der Größe $n\backslash times\; m$ wird durch eine $n\backslash times\; m$-Impulsantwortmatrix $g\_\{nm\}=(nm)^\{-1\}$ beschrieben. Die Bildpunkte des transformierten Bildes sind somit die Mittelwerte ihrer $n\backslash times\; m$ Nachbarn. Box-Filter sind nicht isotrop und nicht abklingend und stellen entgegen einer naiven Einschätzung keinen eigentlichen Tiefpassfilter dar.
• Gaußfilter: Gaußfilter der Größe $n\backslash times\; m$ und einer Varianz $\backslash sigma^2$ werden durch eine Impulsantwortmatrix $g\_\{nm\}=\backslash frac\{1\}\{\backslash sigma\backslash sqrt\{2\backslash pi\}\}e^\{-\backslash frac\{n^2+m^2\}\{2\backslash sigma^2\}\}$ beschrieben. Die Nachbarpunkte des Ausgangsbildes gehen damit nicht wie beim box-filter gleichwertig ein, sondern werden entsprechend einer zweidimensionalen Gaußglocke gewichtet. Gaußfilter sind damit isotrop und abklingend. Sie werden als reales optisches Bauteil als Gauß-Filterung der Fourier-Transformierten in der Apodisation eingesetzt. (Dead_tree_salt_and_pepper_gauss_3x3.png).

• Nichtlineare Glättungsfilter, der Grauwert des aktuellen Pixels wird dabei ersetzt durch
  • Medianfilter
    • den Median der Grauwerte der aktuellen Umgebung (Dead_tree_salt_and_pepper_median_3x3.png)
  • k-zentriertes Mittel
    • den Durchschnitt der k-mittleren Grauwerte
  • k-NN-Filter
    • der Median der k-nächsten Grauwerte der Umgebung
  • $\backslash delta$-Nachbarschaftsfilter
    • den Durchschnitt aller Grauwerte deren Abstand kleiner als $\backslash delta$ zum aktuellen Grauwert ist
  • Lee's Filter
    • $h\_\{nm\}\; =\; \backslash mu\_\{nm\}\; +\; (f\_\{nm\}-\backslash mu\_\{nm\})\backslash cdot\backslash frac\{S\_\{nm\}^2-\backslash sigma^2\}\{S\_\{nm\}^2\}$
    • $\backslash mu\_\{nm\}$ ist Grauwertmittel, $S\_\{nm\}^2$ Grauwertvarianz und $\backslash sigma^2$ Rauschenergie der aktuellen Umgebung
  • Minimum-Varianz-Filter
    • es werden die Grauwertmittlel $\backslash mu\_i$ und Varianzen $S\_i^2$ einiger Teilfenster berechnet
    • $h\_\{nm\}\; =\; \backslash mu\_j$ wird gewählt mit $S\_j^2=min\_i\backslash \{S\_i^2\backslash \}$

= Kantenhervorhebung

= Werden anstatt der absoluten Farbwerte die Abweichung zu den benachbarten Punkten signalisiert, so werden die Objektkonturen, d.h. harte Farbübergänge, eines Bildes hervorgehoben, während weiche Übergänge abgeschwächt werden. Dieser Vorgang wird als Kantenhervorhebung bezeichnet.

Häufig finden folgende FIR-Filter Anwendung:

• Ableitungsfilter
• Laplacefilter

Nichtlineare Kantendetektoren:

• Varianzfilter
• Extremalspannenfilter
• Roberts-Kreuz
• Kirsch-Filter
• Gradientenfilter
  • Prewitt-Filter
  • Sobel-Filter
  • Canny-Filter

Globale Operationen

Fouriertransformation

Werden Bilder als Folge von Farbwerten dargestellt, spricht man von einer Darstellung im Zeit-, Impuls- oder auch Orts-Raum. Aus diesem lassen sich Bilder injektiv in den sogenannten Frequenz-Raum überführen, in welchem nicht mehr die Farbwerte selbst, sondern die Frequenz- und Phasenanteile der zugrundeliegenden Punktfolge gespeichert werden. Große Sprünge innerhalb der Farbwertsequenz stehen dann für hohe Frequenzen, weiche Farbwertübergänge für niedrige Frequenzen. Die Fouriertransformierte gibt schließlich an, wie durch alleinige Überlagerung von Sinusfunktionen verschiedener Frequenz und Phasenwinkel, der ursprüngliche Farbwertverlauf rekonstruiert werden kann.

Der Vorteil dieser Darstellung liegt in ihrer höheren Effizienz bei der Anwendung linearer Filter. Im Impuls-Raum bedeutet die Faltung der das Bild beschreibenden Farbwertfolge $f\_i$ mit der Impuls-Antwort $g\_\backslash mu$(der Transformationsabbildung) des Filters eine Summation (der Gewichtung der Impulse mit der Impulsantwort) über die gesamte Filterlänge ein. Im Frequenzraum hingegen berechnet sich die Faltung als Produkt der Fouriertransformierten mit der Frequenzantwort der Transformation $\backslash Im$. Dieser Vorteil bezüglich der Rechenzeit ist so groß, das selbst der Aufwand der Hin- und Rücktransformation in den Frequenzraum in Kauf genommen werden kann. Dieser Aufwand lässt sich mit Hilfe der schnellen Fourier-Transformation (FFT) deutlich senken.

Von Vorteil ist des weiteren die Vereinfachung, Filter direkt im Frequenzraum designen zu können, anstatt Impulsantworten bauen und ihre Wirkung auf das Spektrum berechnen zu müssen. Damit erübrigt sich gleichsam die Transformation der Impulsantwort $g\_k$ in die zugehörige Frequenzantwort $G(z)$.

Radontransformation

Siehe Hauptartikel: Radontransformation

Transformation

Transformationen von diskreten Bildern gestalten sich gegenüber stetigen Bildern ungemein schwieriger, da Transformationen im Allgemeinen nicht gitterkonform sind. Neben der eigentlichen Transformation muss daher geeignet interpoliert werden. Typische Transformationen sind:

• Skalierung
• Rotation
• Neigung
• (Translation)

Normierung

Eine Meta-Manipulation von Bildern stellt das Feld der Normierung dar. Sie dient dem Abstrahieren von irrelevanten Eigenschaften und der Reduktion von Mustervariabilität. Unproduktive Parameter werden dabei auf einen Standardwert gesetzt. Sie ist damit in der Mustererkennung als Vorverarbeitung von großer Bedeutung. Normierbare Parameter sind unter anderen Intensität, Größe, Lage oder Dicke/Stärke.

Klassifikation

Mit Hilfe der Klassifikation (auch Klassifizierung) soll das breite Informationsspektrum eines Originalbildes durch Einteilung von Farb- oder Grauwerten in Klassen vereinfacht werden. Ziel ist es, Flächen mit gleichen bzw. ähnlichen Eigenschaften zu ermitteln und diese verschiedenen Klassen zuzuordnen. Klassifikatoren werden nicht programmiert, sondern anhand von Beispieldaten trainiert. Es gibt zwei Arten, um Klassifikatoren zu trainieren, das unüberwachte und das überwachte Training.

Unüberwachtes Training

Beim unüberwachten Training werden die Klassen über eine Clusteranalyse gebildet. Die Folge ist, dass die resultierenden Klassen nicht beeinflussbar sind, das Ergebnis also nicht vorhersehbar ist. Wenn eine funktionelle Klassifikation gewünscht ist, d.h. die Klassen bestimmten realen Klassen entsprechen sollen (wichtig in der Fernerkundung, um bspw. Nutzungsklassen auf einem Satellitenbild zu identifizieren) müssen die Ergebnis-Klassen anschließend durch Vergleiche mit den Real-Klassen aufwendig ermittelt werden.

Grundsätzlich gibt es divisive Verfahren, bei denen eine anfängliche Gesamtklasse sukzessive in Unterklassen aufgeteilt wird, und agglomerative Verfahren, bei denen die Klassen (zu beginn so viele, wie es Elemente/Pixel gibt) sukzessive bis zu einer Gesamtklasse zusammengefasst werden.

Überwachtes Training

Das überwachte Training besteht aus zwei Schritten:

1. Zuerst werden die Zielklassen mit den jeweiligen (Farb-)Merkmalen durch eine Auswahl so genannter "Trainingsflächen" festgelegt. Dabei sollte möglichst die gesamte Varianz einer Klasse erfasst, die Auswahl von "Mischpixeln" aber vermieden werden.
2. Anschließend erfolgt die Zuweisung aller Pixel zu den entsprechenden Klassen. Dafür gibt es verschiedenen Methoden, z.B.:

• die Maximum Likelihood-Methode: Zuweisung an Hand der größten Zugehörigkeitswahrscheinlichkeit

• die Minimum Distance-Methode: Zuweisung an Hand des nächstliegenden Klassenmittelpunktes
• die Quader-Methode: Zuweisung an Hand von Ober- und Untergrenzen der einzelnen Klassen
• Künstliche neuronale Netze: iterative Optimierung der Verbindung zwischen den einzelnen Einheiten (Pixeln)

Beispiele

---

Bild:Dead_tree_salt_and_pepper.png|Bild mit künstlich hinzugefügtem "Salt-and-Pepper"- Rauschen. Bild:Dead_tree_salt_and_pepper_gauss_3x3.png|Anwendung eines 3x3 Gaußfilters auf das verrauschte Bild. Bild:Dead_tree_salt_and_pepper_median_3x3.png|Anwendung eines 3x3 Medianfilters auf das verrauschte Bild.

Programme und Libraries zur Bildverarbeitung

---

• ImageJ - freies Bildbe- und Bildverarbeitungsprogramm (Java), siehe auch ImageJ
• DIAS - Interaktives Softwaresystem zur digitalen Bildverarbeitung der Friedrich-Schiller-Universität Jena
• ICE - Bildverarbeitungsbibliothek (C++) der Friedrich-Schiller-Universität Jena
• OpenCV - Open-Source-Bildverarbeitungsbibliothek (C++) von Intel, siehe auch OpenCV
• Orasis3D - freies 3D-Bildverarbeitungsprogramm
• PUMA - Bildverarbeitungsbibliothek (C++) der Universität Koblenz-Landau
• VIGRA - Bildverarbeitungsbibliothek (C++) der Universität Hamburg
• HALCON - Bildverarbeitungsbibliothek der Firma MVTec

Siehe auch

---

• Bildauflösung
• Bildsegmentierung
• Snakes

Weblinks

---

• Normierte Auflösung, Equivalente Auflösung - Maximal notwendige Auflösung für Rasterbilder
• Software-Bibliothek

Bildverarbeitung

Image processing | پردازش تصویر | Traitement d'image | การประมวลผลภาพ | Görüntü işleme | 图像处理