如何编程带阻滤波器

如何编程带阻滤波器

本文介绍了如何编程带阻滤波器的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!

问题描述

你好

我有一个图像处理功课,需要帮助

我有一个吵闹的相机男人照片

当我使用fft变换来确切知道噪声的位置我发现它像一个圆圈所以我必须使用带阻滤波器来过滤图像

我在matlab中的代码是:

i = imread('D:/ MY UNIVERSITY / FIFTH / IMAGE PROCESSING / practical / homework / 9.jpg');

i = im2double(i);

i1 = fft2(i);

i2 = fftshift(i1);

figure,imshow(log(1 + abs(i2)),[]);

[x,y] = meshgrid(-128:127,-128:127);

z = sqrt(x。^ 2 + y。^ 2);

c =(z< 30&z> 40);

c1 = fft2(c);

c2 = fftshift(c1);

figure,imshow(log(1 + abs(c2)),[]);

cm = c2。* i2;

cm1 = ifft2(cm) ;

图,imshow(log(1 + abs(cm1)),[]);



但它不起作用,当我展示过滤器照片时,它是完全黑色的,而不是我需要过滤图像

是不是我做了什么

任何帮助都会受到赞赏

hello
I've got an Image processing homework and need help
I've got a noisy camera man photo
when I used fft transform to know exactly where the noise is located I found it like a circle so I have to filter the image using a band-stop filter
my code in matlab is:
i = imread('D:/MY UNIVERSITY/FIFTH/IMAGE PROCESSING/practical/homework/9.jpg');
i=im2double(i);
i1 = fft2(i);
i2 = fftshift(i1);
figure,imshow(log(1+abs(i2)),[]);
[x,y]=meshgrid(-128:127,-128:127);
z=sqrt(x.^2+y.^2);
c = (z<30 & z>40);
c1 = fft2(c);
c2 = fftshift(c1);
figure,imshow(log(1+abs(c2)),[]);
cm = c2.*i2;
cm1=ifft2(cm);
figure,imshow(log(1+abs(cm1)),[]);

but it didn't work, when I showed the filter photo it was totally black and it isn't what I need to filter the image
Is it wrong what I did
any help will be appreciated

推荐答案


这篇关于如何编程带阻滤波器的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持!

08-20 00:19