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

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Filterung im Frequenzraum; Hoch-/Tiefpassfilter

 

lorenz1980
Forum-Anfänger

Forum-Anfänger


Beiträge: 20
Anmeldedatum: 05.07.08
Wohnort: Berlin
Version: 2010b
     Beitrag Verfasst am: 05.07.2008, 09:37     Titel: Filterung im Frequenzraum; Hoch-/Tiefpassfilter
  Antworten mit Zitat      
Hi, ich bekam folgende Aufgabe in einer meiner Vorlesungen und bin ein wenig (sehr Sad ) unschlüssig über meine Arbeit dazu:

Hier die Aufgabenstellung.

Ich habe dazu folgendes zusammengebastelt:

Code:

function [ output_args ] = prog( pic_in, k0, howto )
%PROG Summary of this function goes here
%   Detailed explanation goes here

pic=imread(pic_in);
[m,n] = size(pic); %Höhe x Breite
%500x500 pic = 1..501x1.501
pic_dbl=double(pic);

for i=1:m
    for j=1:n
        signum(i,j)=(-1)^(i+j);
    end
end

temppic=signum.*pic_dbl;

pic_ft=fft2(temppic);

for u=1:m
    for v=1:n
        pic_ft_c(u+m/2,v+n/2)=pic_ft(u, v);
    end
end

switch howto
    case 'itp'
        %idealer Tiefpassfilter
        h=zeros(u+m/2,v+n/2);
        for i=1:u+m/2
            for j=1:v+n/2
                if sqrt(i^2+j^2)<=k0
                    h(i,j)=1;
                end
            end
        end
    case 'ihp'
        %idealer Hochpassfilter
        h=zeros(u+m/2,v+n/2);
        for i=1:u+m/2
            for j=1:v+n/2
                if sqrt(i^2+j^2)>k0
                    h(i,j)=1;
                end
            end
        end
    case 'gtp'
        %Gauß-TP
        h=exp(-sqrt((u+m/2)^2+(v+n/2)^2)/(2*k0^2));
    case 'ghp'
        %Gauß-HP
        h=1-exp(-sqrt((u+m/2)^2+(v+n/2)^2)/(2*k0^2));
end

pic_temp=ifft2(h.*pic_ft_c);
imshow(real(pic_temp))
pause
%pic_new=signum.*pic_temp;

%imshow(real(pic_new))
%pause
figure(1); close;
 


Mein(e) Problem(e):
Wozu verschiebe ich das Bild um m/2, n/2?
Das führt dazu das das Ausgangsbild 50% größer ist als der Input.
Ebenso führt das dazu, dass man das neue Bild nicht mit signum multiplizieren kann, da die Matrizen natürlich unterschiedlich groß sind.


Hoffe, dass jemand meine Unklarheiten beseitigen und mir helfen kann, das Ding auf die Beine zu bekommen :-/



PS: Gibt es eine Matlab-Funktion, die das macht, was ich machen soll (zum Vergleichen der Ergebnisse)?[/code]
Private Nachricht senden Benutzer-Profile anzeigen


Neues Thema eröffnen Neue Antwort erstellen



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 | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2024 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.