Gibt es analytische Ansätze für den Entwurf von quadratischen Polynomfiltern?

19

Das folgende Papier beschreibt eine Anwendung des Teager-Kaiser-Energieoperators zur Röntgenbildverbesserung:

Reinhard Bernstein, Michael S. Moore und Sanjit K. Mitra, "Einstellbare quadratische Filter zur Bildverbesserung" Proc. Internationale IEEE-Konferenz für Bildverarbeitung (ICIP), Santa Barbara, CA, vol. 1, S. 287-290, Okt. 1997. http://vision.ece.ucsb.edu/publications/view_abstract.cgi?52

Die Autoren entwickeln Intuition für das Verhalten des Filters durch Analogie mit einem ähnlichen linearen Operator (dh "Somit entspricht die Ausgabe eines Teager-Filters ungefähr einer mit dem lokalen Mittelwert gewichteten Hochpassfilterantwort" ). Präzisionshalber meine ich mit quadratischen Polynomfiltern nichtlineare, nichtrekursive Filter, die durch eine verkürzte Volterra-Reihe wie folgt vollständig charakterisiert werden können (für den 1D-Fall):

y(n)=∑m1=0N1-1h1(m1)x(n-m1)+∑m1=0N2-1.∑m2=0N2-1h2(m1,m2)x(n-m1)x(n-m2)

Es scheint, dass die meisten Ansätze für den Entwurf von Polynomfiltern niedriger Ordnung Systemidentifikations-Frameworks beinhalten, aber ohne tiefes Verständnis, warum die geschätzten Filter funktionieren. Sind analytische Ansätze, die auf linearen Analogien basieren, derzeit Stand der Technik, oder gibt es bekannte mathematische Werkzeuge, die verwendet werden können?

filters Datageist
quelle
Ich habe vorher den Teager-Kaiser-Energieoperator verwendet. Ich weiß, dass es sich hervorragend dazu eignet, kurze Impulse aus verrauschten Daten herauszuholen (das Gegenteil eines Medianfilters). Es kann auch rosa Rauschen weiß machen. Für reine Töne ist die Ausgabe eine Konstante (die Energie des Tons).
Spacey
@ Mohammad Interessant. Ich habe es nur als Beispiel für die Frage verwendet, aber die rosa-zu-weiß-Eigenschaft war mir nicht bewusst. Dank dafür!
Datageist
Hier ist meine Matlab-Implementierung, wenn Sie damit herumspielen möchten. ( dl.dropbox.com/u/4724281/TKEO.m )
Spacey
@datageist Hast du jemals eine Antwort darauf gefunden? Kannst du posten, was du hast? Dies ist die am häufigsten gestellte unbeantwortete Frage! =)
Phonon
@Phonon Ich habe irgendwann eine Art Teilantwort gefunden, aber es war nicht sehr befriedigend. Ich werde versuchen, es bald aufzuschreiben (ish).
Datageist

Antworten:

1

Keine wirkliche Antwort (daher ist dies Community-Wiki), aber ich dachte, wir sollten den Code von @ Mohammad erfassen:

%Mohammad Z

%Teager-Kaiser Non-Linear Energy Operator. 
function [out] = TKEO(x)
    N = length(x);
    x = x(:).';
    temp = x(2:N-1).^2 - x(3:N).*x(1:N-2);
    out = [temp(1) temp temp(end)];    
end
Peter K.
quelle