Wie ist die Beziehung zwischen Sigma und Radius? Ich habe gelesen, dass Sigma dem Radius entspricht. Ich verstehe nicht, wie Sigma in Pixeln ausgedrückt wird. Oder ist "radius" nur ein Name für Sigma, der nicht mit Pixeln zusammenhängt?
Hier spielen drei Dinge eine Rolle. Die Varianz ( ), der Radius und die Anzahl der Pixel. Da es sich um eine zweidimensionale Gaußfunktion handelt, ist es sinnvoll, stattdessen von der Kovarianzmatrix zu sprechen . Wie dem auch sei, diese drei Konzepte sind schwach miteinander verwandt.σ2Σ
Zunächst wird der 2-D-Gauß durch die Gleichung gegeben:
G( z ) = 1( 2 π)2| Σ |-------√e- 12( z - μ )TΣ- 1 ( z - μ )
Wobei ein Spaltenvektor ist, der die und Koordinate in Ihrem Bild enthält. Also ist , und ist ein Spaltenvektor, der den Mittelwert Ihrer Gaußschen Funktion in und Richtung . x y z = [ x y ] μ x y μ = [ μ x μ y ]zxyz=[xy]μxyμ=[μxμy]
Beispiel:
Nehmen wir nun an, wir setzen die Kovarianzmatrix und . Ich werde auch die Anzahl der Pixel auf x einstellen . Außerdem wird mein 'Gitter', in dem ich dieses PDF auswerte , sowohl in als auch in von bis gehen . Dies bedeutet, dass ich eine Rasterauflösung von . Das ist aber völlig willkürlich. Mit diesen Einstellungen erhalte ich das Bild der Wahrscheinlichkeitsdichtefunktion auf der linken Seite. Nun, wenn ich die "Varianz" (wirklich die Kovarianz) so ändere, dass μ = [ 0 0 ] 100 100 - 10 10 x y 10 - ( - 10 )Σ=[1001]μ=[00]100100−1010xyΣ=[ 9 0 0 9 ]10−(−10)100=0.2Σ=[9009] und alles andere gleich, ich bekomme das Bild rechts.
Die Anzahl der Pixel ist für beide gleich, x , aber wir haben die Varianz geändert. Angenommen, wir machen stattdessen dasselbe Experiment, verwenden jedoch stattdessen x Pixel, aber ich habe immer noch einen Wert zwischen und . Dann hat mein Gitter eine Auflösung von . Wenn ich die gleichen Covarianzen wie zuvor verwende, erhalte ich Folgendes:100 20 20 - 10 10 10 - ( - 10 )1001002020−101010−(−10)20=1
So müssen Sie das Zusammenspiel dieser Variablen verstehen. Wenn Sie den Code mögen, kann ich das auch hier posten.
Wie wähle ich Sigma?
Die Wahl der Varianz / Kovarianz-Matrix Ihres Gauß-Filters ist stark anwendungsabhängig. Es gibt keine "richtige" Antwort. Das ist so, als würde man fragen, welche Bandbreite man für einen Filter wählen soll. Auch dies hängt von Ihrer Anwendung ab. In der Regel möchten Sie einen Gauß-Filter so auswählen, dass Sie eine beträchtliche Anzahl von Hochfrequenzkomponenten in Ihrem Bild auf Null setzen. Eine Sache, die Sie tun können, um ein gutes Maß zu erhalten, ist, die 2D-DFT Ihres Bildes zu berechnen und ihre Koeffizienten mit Ihrem 2D-Gauß-Bild zu überlagern. Hier erfahren Sie, welche Koeffizienten stark bestraft werden.
Wenn Ihr Gauß-Bild beispielsweise eine so breite Kovarianz aufweist, dass es viele Hochfrequenzkoeffizienten Ihres Bildes umfasst, müssen Sie seine Kovarianzelemente verkleinern.
Der Parameter Sigma reicht aus, um die Gaußsche Unschärfe aus einer kontinuierlichen Sicht zu definieren. In der Praxis sind Bilder und Faltungskerne jedoch diskret. Wie wählt man eine optimale diskrete Approximation des kontinuierlichen Gaußschen Kernels?
Die diskrete Approximation ist näher am kontinuierlichen Gaußschen Kernel, wenn ein größerer Radius verwendet wird. Dies kann jedoch zu Lasten der zusätzlichen Berechnungsdauer gehen.
Idealerweise würde man einen Wert für Sigma auswählen und dann einen Radius berechnen, der es ermöglicht, den entsprechenden kontinuierlichen Gaußschen Kern genau darzustellen. Für einen gegebenen Approximationsfehler muss der Radius umso größer sein, je größer das Sigma ist.
Interessanterweise kann dies sehr kompliziert werden, um es richtig zu machen. Wenn Sie die Gaußsche Matrix konstruieren, ist dies die beste Lösung, um den kontinuierlichen Kernel abzutasten, oder gibt es bessere Näherungen? Wie normalisiere ich den berechneten diskreten Kernel, um das Abschneiden zu berücksichtigen? etc.
Als Referenz bietet die Funktion Gaußsche Matrix in Mathematica verschiedene Möglichkeiten zur Berechnung einer diskreten Gaußschen Matrix, z. B. unter Verwendung der diskreten Bessel-Approximation. Standardmäßig ist der Radius 2 * Sigma, was bedeutet, dass bei Sigma 1 die Matrix 5x5 ist.
quelle
Es stellt sich heraus, dass die Zeilen des Pascalschen Dreiecks sich recht gut einem Gaußschen annähern und den praktischen Vorteil haben, dass sie ganzzahlige Werte haben, deren Summe eine Potenz von 2 ist (wir können diese Werte genau als Ganzzahlen, Festpunktwerte oder Gleitkommazahlen speichern ). Nehmen wir zum Beispiel an, wir möchten einen 7x7-Gauß-Kernel konstruieren. Dazu können wir die 7. Zeile des Pascal-Dreiecks wie folgt verwenden:
Beachten Sie, dass dieser Filter den geringsten Einfluss auf die Ecken hat, während der Wert für die ganze Zahl erhalten bleibt. Sie können den Mittelwert 20/64 verwenden, um das entsprechende Standardabweichungs-Sigma zu bestimmen, das in diesem Fall 64 / (20 * sqrt (2 * pi)) = 1,276 für den angenäherten Gaußschen Wert ist. Sie können den Gaußschen Wert grafisch darstellen, um zu sehen, dass dies eine hervorragende Anpassung ist.
Ein guter Ausgangspunkt für die Bestimmung einer vernünftigen Standardabweichung für einen Gaußschen Kern ist also das Pascalsche Dreieck (auch bekannt als Binomialkoeffizienten ) - für ein (N + 1) x (N + 1) -Filter, das der obigen Konstruktionsverwendung entspricht
Wolfram Alphas Gaußsche Matrix [3] verwendet nur r / 2 = 1,5. Seltsamerweise liefert GaussianMatrix [{3,1.276}] nicht den gleichen 2D-Filter wie meine und ist für x, y zwischen -3 und 3 nicht der folgende:
Ich bin mir nicht sicher warum nicht? Mein 2D-Filter passt hervorragend.
quelle