Low Pass Filter in 2D

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

pospiech

Forum-Century


	Beiträge: 130

	Anmeldedatum: 03.09.07

	Wohnort: ---

	Version: ---

Verfasst am: 11.01.2010, 14:31 Titel: Low Pass Filter in 2D

Ich habe mal ein Beispiel zusammengebastelt, das zeigt was ich filtern möchte.

Da ich mit allen Dokumentationen zu Filtern von Matlab nicht weiterkommen, da diese alle nur für 1D Daten ausgelegt sind und praktisch keine Beispiele enthalten (die ich verstehe) habe ich mal probiert einen Filter selber zu basteln (komplett naiv ausgedacht).

Ist das ein brauchbarer Ansatz oder geht es grundsätzlich anders und besser?

Code:

clear; clc;

N = 100;
xAxis = linspace(-5,5,N);
yAxis = linspace(-5,5,N);

data = zeros([length(xAxis) length(yAxis)]);

%% generate gaussian with random noise
w = 2;
for i=1:length(xAxis)
for j=1:length(yAxis)
x = xAxis(i);
y = yAxis(j);
r = sqrt(x^2+y^2);
data(i,j) = exp(-(r^2/w^2));
data(i,j) = data(i,j) + rand(1)*0.1;
end
end

%% generate super-gaussian filter function
Nf = 256;
xAxisf = linspace(-5,5,Nf);
yAxisf = linspace(-5,5,Nf);
filter = zeros([Nf Nf]);
wfilter = 2;
for i=1:Nf
for j=1:Nf
x = xAxisf(i);
y = yAxisf(j);
r = sqrt(x^2+y^2);
filter(i,j) = exp(-(r^2/wfilter^2)^4);
end
end

%% Apply filter

% FFT
Xf = fftshift(data);
Xf = fft2(data,Nf,Nf);
Xf = fftshift(Xf);
% filter
for i=1:Nf
for j=1:Nf
Xf(i,j) = Xf(i,j) * filter(i,j)^2;
end
end

figure(2);
surf(abs(Xf)); shading interp;

% iFFt
Xf = ifftshift(Xf);
Xif = ifft2(Xf);
Xif = Xif(1:N,1:N);
Xif = abs(Xif);

figure(5); surf(filter); shading interp;
figure(1); surf(data); shading interp;
figure(3); surf(Xif); shading interp;
figure(4); surf(Xif-data); shading interp;

Funktion ohne Link?

Das hauptproblem hierbei ist die Anpassung der Filterbreite


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.