Lineare Filter

Ableitungen

Man versuche nicht die Fourier-Reihen zu differenzieren !

Formel für Ableitunge gewinnt man durch die Konstruktion eines Interpolationspolynoms, das Bilden der Ableitung des
Polynoms und dem Auswerten an einer bestimmte Position.

1. Ordnung

LinearFilters_1.gif

LinearFilters_2.gif

1. Ableitung LinearFilters_3.gif LinearFilters_4.gif

2. Ordnung

LinearFilters_5.gif

LinearFilters_6.gif

1. Ableitung LinearFilters_7.gif LinearFilters_8.gif LinearFilters_9.gif
2. Ableitung LinearFilters_10.gif LinearFilters_11.gif LinearFilters_12.gif

3. Ordnung

LinearFilters_13.gif

LinearFilters_14.gif

1. Ableitung LinearFilters_15.gif LinearFilters_16.gif LinearFilters_17.gif LinearFilters_18.gif LinearFilters_19.gif
2. Ableitung LinearFilters_20.gif LinearFilters_21.gif LinearFilters_22.gif LinearFilters_23.gif LinearFilters_24.gif

Rauschen

Ableitungen sind extrem empfindlich gegen Rauschen.

In[12]:=

LinearFilters_25.gif

Out[12]=

LinearFilters_26.gif

Gauß Filter

Einsatz für die Rauschunterdrückung.

LinearFilters_27.gif

LinearFilters_28.gif

In[16]:=

LinearFilters_29.gif

Out[16]=

LinearFilters_30.gif

Ableitungen

Erste Ableitungen

In[18]:=

LinearFilters_31.gif

Out[18]=

LinearFilters_32.gif

In[20]:=

LinearFilters_33.gif

Out[20]=

LinearFilters_34.gif

Laplace Filter

In[15]:=

LinearFilters_35.gif

Out[15]=

LinearFilters_36.gif

Nullinie (Kanten) für Gaussfilterweiten 1 (Schwarz),2 (Rot), 4 (Grün), 8 (Blau)

LinearFilters_37.gif

Das kontinuierliche Variiren der Breite des Gauß-Kerns wird als "scale space" bezeichnet. Es verbindet den lokalen Charakter den Ableitung mit dem globalen Charakter der Faltung mit einem Gauß-Filter.

LinearFilters_38.gif

Laplace Gauss Filter werden auch mit mehr als einer Filterweite verwendet, und das Maximum über die FIlterweiten verwendet.

LinearFilters_39.gif

LinearFilters_40.gif

Diffusion & Gauß Filter

Eine Diffusions- oder Wärmeleitungsgleichung hat im einfachsten (zweidimensionalen) Fall die Form:

LinearFilters_41.gif

oder mit LinearFilters_42.gif und LinearFilters_43.gif

LinearFilters_44.gif

Dabei sei u(x,y) das Grauwertbild. Die Lösung dieser Gleichung im zweidimensionalen lautet

LinearFilters_45.gif

Für einen Anfangswert LinearFilters_46.gif erhält man

LinearFilters_47.gif

Dies ist aber gerade die kontinuierliche Version der Faltung eines Bildes LinearFilters_48.gif  mit einem Gauß-Kern der Breite LinearFilters_49.gif.

Isotrope Nichtlineare Diffusion

Diffusion/Gauß-Filter habe unterdrücken das Rauschen in Bildern sehr gut, zerstören aber die Bildkannten.  Perona und Malik führten dazu eine orts- und bildabhänginge Diffusionskonstante g ein, die vom Betrag des Gradienten abhängt

LinearFilters_50.gif

Die Funktion g(z) hat dabei die Form

LinearFilters_51.gif

LinearFilters_52.gif

LinearFilters_53.gif

Lösung der nichtlinearen Diffusionsgleichung für   ν=0.1,0.5,1,2,4,8 bei t=0.8 und w=7

LinearFilters_54.gif

Veränderung der Lösung mit der Zeit t=1,3,5,7,9,11,13,15 mit ν=0.1

LinearFilters_55.gif

Literatur

P.Perona, J.Malik, Scale space and edge detection using anisotropic diffusion, Proc.IEEE Comp.Soc.Workshop on Computer Vision (Miami Beach, Nov.30 – Dec.2, 1987), IEEE Computer Society Press, Washington, 16 –22, 1987.

P.Perona, J.Malik, Scale space and edge detection using anisotropic diffusion, IEEE Trans.Pattern Anal.Mach.Intell., Vol.12, 629 –639, 1990

Vorbereitung

In[3]:=

LinearFilters_56.gif

In[4]:=

LinearFilters_57.gif

Out[4]=

LinearFilters_58.gif

In[5]:=

LinearFilters_59.gif

In[6]:=

LinearFilters_60.gif

In[7]:=

LinearFilters_61.gif

In[8]:=

LinearFilters_62.gif

Out[8]=

LinearFilters_63.gif

In[9]:=

LinearFilters_64.gif

In[10]:=

LinearFilters_65.gif

In[11]:=

LinearFilters_66.gif

Out[11]=

LinearFilters_67.gif

LinearFilters_68.gif

LinearFilters_69.gif

LinearFilters_70.gif

LinearFilters_71.gif

LinearFilters_72.gif

LinearFilters_73.gif

LinearFilters_74.gif

LinearFilters_75.gif

LinearFilters_76.gif

LinearFilters_77.gif

LinearFilters_78.gif

LinearFilters_79.gif

Spikey Created with Wolfram Mathematica 7.0