FFT auf einem nicht rechteckigen Teil des Bildes

7

Ich muss eine FFT für ein Bild durchführen, um das Rauschen zu reduzieren, aber das Problem ist, dass ich nicht das vollständige Bild benötige, sondern nur einen Kreis in der Mitte. Die Ränder sind ein festes Rig, daher interessiert mich nicht, was es anzeigt, aber es hat einen Einfluss auf das Ergebnis der FFT.

Gibt es eine Möglichkeit, einfach ein kreisförmiges Teil auszuschneiden und dieses für die FFT zu verwenden? Beachten Sie, dass bei Verwendung von schwarzem Hintergrund die Kante zwischen Hintergrund- und Bilddaten erhebliche Auswirkungen hat.

fft

— SinisterMJ
quelle

7

Anstatt eine harte Kante zwischen den interessierenden Bilddaten und dem Hintergrund zu haben, könnten Sie eine zweidimensionale verjüngte Fensterfunktion verwenden , wie dies häufig bei der Spektralanalyse der Fall ist. Sie könnten zunächst ein Gaußsches Fenster ausprobieren , das für einen zweidimensionalen Fall ungefähr so aussehen würde:

w [x, y] = e^{- \frac{{(x - \frac{N_{x} - 1}{2})}^{2}}{2 {(σ_{x} \frac{N_{x} - 1}{2})}^{2}}} e^{- \frac{{(y - \frac{N_{y} - 1}{2})}^{2}}{2 {(σ_{y} \frac{N_{y} - 1}{2})}^{2}}}

$w[x,y] = e^{-\frac{\left(x-\frac{N_x-1}{2}\right)^2}{2\left(\sigma_x \frac{N_x-1}{2}\right)^2}} e^{-\frac{\left(y-\frac{N_y-1}{2}\right)^2}{2\left(\sigma_y \frac{N_y-1}{2}\right)^2}}$

$N_x$ und sind die Abmessungen der gewünschten Transformation in bzw. Richtung, und und sind Parameter, mit denen Sie die Form des Fensters steuern können. Bei kleinen Werten konzentriert sich der größte Teil der Energie in der Fensterfunktion auf die Mitte, wobei dieser Effekt mit zunehmendem abnimmt . $N_y$ $x$ $y$ $\sigma_x$ $\sigma_y$ $\sigma$ $\sigma$

— Jason R.
quelle

Wäre ein flacheres Fenster nicht besser, da er das Zeug in der Mitte nicht verzerren möchte?

— Jim Clay

Sicher; Das Obige wurde nur als Beispiel für ein Fenster gegeben. Es gibt viele Fensterfunktionen, die verwendet werden könnten, obwohl sie normalerweise nur für den eindimensionalen Fall beschrieben werden.

— Jason R

Okay danke. Ich habe bereits über die Verwendung einer Fensterfunktion nachgedacht, wusste aber nicht, ob es etwas Besseres / völlig anderes gibt.

— SinisterMJ

3

Verwenden Sie eine ~~konische Flat-Top-~~ Fensterfunktion mit einer flachen Top- und Kantenverjüngung, z. B. ein sich verjüngendes Cosinusfenster oder ein Tukey-Fenster (drehen Sie es um eine Polarachse, um es zu einer kreisförmigen 2D-Schablone zu machen) und ein Null-Pad, wie es für jedes Fasten erforderlich ist Faltungsfilterung. Das Fenster, das an den Kanten auf Null geht, verringert den Effekt der Kantendiskontinuitäten vom Kreis auf die FFT. Mischen Sie nach dem Filtern oder einer anderen Verarbeitung das Ergebnis wieder, indem Sie die subtraktive Umkehrung Ihrer Fensterfunktion für die vorhandenen Bilddaten verwenden.

— hotpaw2
quelle

0

Eine einfache Lösung könnte darin bestehen, die Randbedingungen Ihres nicht rechteckigen interessierenden Bereichs in das Rechteck um ihn herum zu verteilen.

Sie können dies mit einem Algorithmus für den nächsten Nachbarn durchführen.

— Danny Varod
quelle