Gaussian High Pass Filter

Hello, I try to apply Gaussian High Pass Filter function

a=imread('img'); figure(1) imshow(a) [m n]=size(a); f_transform=fft2(a); f_shift=fftshift(f_transform); p=m/2; q=n/2; d0=70; for i=1:m for j=1:n distance=sqrt((i-p)^2+(j-q)^2); low_filter(i,j)=exp(-(distance)^2/(2*(d0^2))); end end filter_apply=f_shift.*low_filter; image_orignal=ifftshift(filter_apply); image_filter_apply=abs(ifft2(image_orignal)); figure(2) imshow(image_filter_apply,[])

but I got this error

??? Error using ==> times Matrix dimensions must agree.

can someone help me

[islam habeeb]

Implement non Gaussian for IPCA in Matlab | Crunch Modo

Answer 1

[Youssef Khmou]

hi lina,

I tried your program and it works fine, replace the first command with :

 a=im2double(imread('circuit.tif'));

I think the error you get is due to the nature of the sample you use, maybe its RGB image, if so you need to change it to grayscale image or process each channel individually, try this as first attempt :

 a=im2double(rgb2gray(imread('img')));

[lina]

Well done thank you brother :)

[elfrida]

if format of my picture is .jpg, can I use this program? or I must convert to .tif??? thankyou

[Image Analyst]

No, that has nothing to do with it. No matter what format the image is stored in on disk, once imread gets done with it, it's just a regular array of numbers.