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

1 个评论

Implement non Gaussian for IPCA in Matlab | Crunch Modo

请先登录,再进行评论。

 采纳的回答

Youssef  Khmou
Youssef Khmou 2013-7-27
编辑:Youssef Khmou 2013-7-27
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')));

3 个评论

Well done thank you brother :)
if format of my picture is .jpg, can I use this program? or I must convert to .tif??? thankyou
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.

请先登录,再进行评论。

更多回答(0 个)

类别

帮助中心File Exchange 中查找有关 Images 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by