MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn moreOpportunities for recent engineering grads.

Apply Today**New to MATLAB?**

Asked by lina
on 27 Jul 2013

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

*No products are associated with this question.*

Answer by Youssef Khmou
on 27 Jul 2013

Edited by Youssef Khmou
on 27 Jul 2013

Accepted answer

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')));

elfrida
on 13 Aug 2013

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

Image Analyst
on 13 Aug 2013

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.

## 1 Comment

## islam habeeb (view profile)

Direct link to this comment:http://www.mathworks.com/matlabcentral/answers/83339#comment_265460

Implement non Gaussian for IPCA in Matlab | Crunch Modo