Fourier transform - filtering, Electrical Engineering

Fourier transform  (filtering)

(i) Perform low pass filtering in the frequency domain. Write and m-file lowfft.m which does this operation.

lowfft.m

function lowfft
im=imread('class_f.png');
imd=double(im);FI=fft2(imd);
phase=angle(FI);
amplitude=abs(FI);
figure;imshow(log(amplitude),[]);
figure;imshow(phase,[]);

%Where is the energy concentrated? where is low frequency in that spectrum?

help fftshift
figure;imshow(log(fftshift(amplitude)),[]);
figure;imshow(fftshift(phase),[]);

%Where is the energy concentrated? Where is the low frequency in that spectrum?

FIc=fftshift(FI);
amplitude=abs(FIc);
figure;imshow(log(amplitude),[]);
[h w]=size(FIc);
w2=uint8(w/2)
h2=uint8(h/2)
s=50;
mask=zeros(h,w);
mask(h2-s:h2+s,w2-s:w2+s)=1;
figure;imshow(mask,[]);

%What is the role of the box filter here?

FFIc=FIc.*mask;
phase=angle(FFIc);
amplitude=abs(FFIc);
figure;imshow(log(amplitude),[]);
figure;imshow(phase,[]);
FFI=ifftshift(FFIc);
RI=ifft2((FFI));
figure;
imshow(real(RI),[]);

%What happens when you apply the mask filer?
%Where do the artifacts come from?

(ii) Perform high pass filtering in the frequency domain. Write and m-file highfft.m which does this operation.
replace the mask by

mask=ones(h,w);
mask(h2-s:h2+s,w2-s:w2+s)=0;

(iii) Filter the image for different values of parameter s.
What happens when you vary the size of the box filter?

Posted Date: 2/26/2013 12:05:39 AM | Location : United States







Related Discussions:- Fourier transform - filtering, Assignment Help, Ask Question on Fourier transform - filtering, Get Answer, Expert's Help, Fourier transform - filtering Discussions

Write discussion on Fourier transform - filtering
Your posts are moderated
Related Questions

Can someone please guide me on how to compute a nash equilibrium in matlab

Switching Characteristics During Turn off The techniques  used for turning off a thyristor is known  as  commutation techniques. A  thyristor can be can be  turn off by reducin

what are the functions of interpolar winding and compensating winding?

in DC series generator we consider Eg=V+Ia(Ra+Rse).But in DC shunt generator why it is only Eg=V+IaRa and not Eg=V+Ia(Ra*Rsh/Ra+Rsh)

application and demerits

Q. A negative impedance converter is used, as shown in Figure. Show that the load current i L is given by v in /R, which is independent of Z L . Note that since the load sees a cu

Equipment parameter data: The schematic diagrams for existing substations have to be prepared along with information of power transformer rating and numbers, impedance values,