Anisotrope Diffusion

Die ursprüngliche Variante der Diffusionsgleichung

09NonlinearFilters_1.gif

mit 09NonlinearFilters_2.gif und 09NonlinearFilters_3.gif betrachtete D als Skalar. Ist 09NonlinearFilters_4.gif ein Tensor (eine Matrix)

09NonlinearFilters_5.gif

so handelt es sich um eine anisotrope Diffusion. Macht man den Diffusionstensor 09NonlinearFilters_6.gif von den Bilddaten abhänging, so kann man lineare Struktuen verstärken. Die beiden Eigenvektoren des Diffusionstenors seien

09NonlinearFilters_7.gif

09NonlinearFilters_8.gif mit i=1,2

Zur Konstruktion eines geeigneten Diffusionstensors startet man mit einem geglätteten Bild 09NonlinearFilters_9.gif, also ein mit einem Gauß Filter der Weite σ gefalteten Bild. Aus dem Gradienten 09NonlinearFilters_10.gif lässt sich eine Matrix 09NonlinearFilters_11.gif konstruieren. Diese Matrix ist positive, und symmetrisch und die Eigenvektoren sind parallel und senkrecht zum Gradienten. Die Eigenwerte sind 09NonlinearFilters_12.gifund 0. Diese matrix wird Komponentenweise mit einem zweiten Gauß-Kern geglättet und man erhält

09NonlinearFilters_13.gif

Die Eigenwerte sind

09NonlinearFilters_14.gif

Der erste Eigenvektor dazu lautet

09NonlinearFilters_15.gif

Für die Koherenzverbesserung wählt man eine Diffusionstensor, der im Koordiantensystem der beiden Eigenvektoren die  beiden Eigenwerte

09NonlinearFilters_16.gif

hat.

09NonlinearFilters_17.gif

09NonlinearFilters_18.gif

Sollen die Kanten verstärt werden wählt man als Eigenwerte des Diffusionstensors und (ρ=0)

09NonlinearFilters_19.gif

09NonlinearFilters_20.gif

09NonlinearFilters_21.gif

Shock Filter

Shock filter sind eine weitere Methode zur Bildverbesserung, wie der Name schon sagt erzeugen Shock-Filter homogene Gebiete im Bild, die von scharfen Kanten getrennt sind. Die Kanten entstehen an den Kanten des Bildes, dazwischen wird das Bild geglättet.

09NonlinearFilters_22.gif

09NonlinearFilters_23.gif

Für shock Filter wird die Gleichung

09NonlinearFilters_24.gif

mit dem Gradienten 09NonlinearFilters_25.gif und dem Laplace Operator 09NonlinearFilters_26.gif. Der Laplace-Operator Δ u  dient dabei als Kantendetektor. Die Kanten-Detektion is noch effizienter in Richtung des Gradienten η|| ∇ u und mit einer Glättung von u mit einem Gauß-Kern der Weite σ. Mit 09NonlinearFilters_27.gif und der zweifachen Richtungsableitung 09NonlinearFilters_28.giflässt sich der Filter als  

09NonlinearFilters_29.gif

schreiben.

09NonlinearFilters_30.gif

09NonlinearFilters_31.gif

09NonlinearFilters_32.gif

09NonlinearFilters_33.gif

Im Falle, das lineare Strukturen für die Koherenzverbesserung verstärkt werden sollen, wird nicht die Richtung des geglätteten Gradienten η sonder der dominate Eigenvektror des geglätteten Strukturtensors 09NonlinearFilters_34.gif verwendet.

09NonlinearFilters_35.gif

09NonlinearFilters_36.gif

09NonlinearFilters_37.gif

09NonlinearFilters_38.gif

Total Variation Filter

Diese Filter basieren auf einem Variationsprinzip, die Lösung soll das Funktional

09NonlinearFilters_39.gif

bezüglich u extremal werden lassen. 09NonlinearFilters_40.gif ist dabe das Ausgangsbild und λ ein freier Lagrange Parameter, der bestimmt wie weit sich die Lösung vom Ausgangsbild entfernen kann. Dies führt zu der Gleichung

09NonlinearFilters_41.gif

Eine Möglichkeit diese Gleichung zu lösen besteht darin, eine "Zeitabhängikeit" einzu führen und im Grenzwert t→∞  die Lösung zu erhalten

09NonlinearFilters_42.gif

Der Lagrange Parameter λ kann aus der Standardabweichung des Rauschens im Bild abgeschätzt werden.

09NonlinearFilters_43.gif

09NonlinearFilters_44.gif

Literatur

J.Weickert : Anisotropic Diffusion in Image Processing.ECMI Series, Teubner, Stuttgart, 1998.

J.Weickert : Coherence - enhancing shock filters.In B.Michaelis, G.Krell (Eds.) : Pattern Recognition.Lecture Notes in Computer Science, Vol.2781, Springer, Berlin, 1 - 8, 2003.

Tony F.Chan, Stanley Osher,Jianhong Shen, The Digital TV Filter and Nonlinear Denoising, IEEE  Transactions on image processing, Vol. 10, No. 2, Feb. 2001

Vorbereitung

09NonlinearFilters_45.gif

09NonlinearFilters_46.gif

09NonlinearFilters_47.gif

Spikey Created with Wolfram Mathematica 7.0