频率域滤波之带阻滤波器matlab代码

代码主要实现了两点功能

1、增加频率域噪声

2、用带阻滤波器去掉频率域噪声

效果：

原图像和加噪声后的图像（多了很多波纹）

频率域图像

带阻滤波器

最后的图像

上代码：

clear;
close all;
clc;I=imread('m.tif');
I=im2double(I);
[m,n]=size(I);r = 0:m-1;%Row coordinates.
c = 0:n-1;%Column coordinates.
[C, R] = meshgrid(c, r);
A=0.01;
%半径2
r=2;
%划区域10个
nn=10;
f = zeros(m, n);
for i=1:nntemp=(i-0.5)*pi/nn;u0=r*sin(temp);v0=r*cos(temp);f = f+A*sin(u0*R + v0*C);
end
J=I+f;
figure;
imshow(J);%傅里叶变换
IF=fftshift(fft2(I));
JF=fftshift(fft2(J));
%对数变换将高频细节扩展出来
IF=log(1+abs(IF));
JF=log(1+abs(JF));
figure;
subplot(1,2,1), imshow(IF,[]), title('原图像');
subplot(1,2,2), imshow(JF,[]), title('周期噪声图像');%% %%%%%%%%%%%%%%%%%%%%%%%%%%
%处理J
k=fft2(J);
g=fftshift(k);
imshow(log(abs(g)+1),[]);[N1,N2]=size(g);W=10;
n=2;
d0=148;
n1=round(N1/2);
n2=round(N2/2);
Homo = zeros(N1, N2);
for i=1:N1 for j=1:N2d=sqrt(((i-n1)*N2/N1)^2+((j-n2))^2);Home(i,j)=1/(1+((d*W)/(d^2-d0^2))^(2*n));end
end
figure;
mesh(abs(Home(1:10:N1,1:10:N2)));
figure;
imshow(Home,[]);%y=g.*Home;
y=g.*Home;
figure;
imshow(log(abs(y)+1),[]);y=ifftshift(y);
E1=ifft2(y);
E2=real(E1);subplot(1,3,1), imshow(I,[]), title('原图像');
subplot(1,3,2), imshow(J,[]), title('周期噪声图像');
subplot(1,3,3), imshow(E2,[]), title('处理后图像');

问题：

为什么加噪声以后频率域会出现十字形状，为什么对周围造成了影响？

频率域滤波之带阻滤波器matlab代码相关推荐

频率域滤波matlab函数,频率域滤波（2） - osc_t0xmr011的个人空间 - OSCHINA - 中文开源技术交流社区...
一.频率域滤波的基本步骤: 1)使用函数tofloat把输入图像转换为浮点图像(im2double函数也可以) [f,revertclass] = tofloat(f) 2)使用函数paddedsiz ...
MATLAB中实现图像的空间域滤波和频率域滤波
1. 空间域滤波空间域滤波是指在图像空间中借助模板对图像领域进行操作,处理图像每一个像素值.主要分为线性滤波和非线性滤波两类,根据功能可分为平滑滤波器和锐化滤波器.平滑可通过低通来实现,平滑的目的有 ...
matlab频率域滤波器,频率域滤波的MATLAB设计与实现_课程设计
频率域滤波的MATLAB设计与实现_课程设计综合课程设计设计题目频率域滤波的 MATLAB 设计与实现专业名称班级学号学生姓名指导教师设计时间目录摘要 .- 3 -1. 数字图像处理 - 1 ...
频率域滤波matlab函数,6、频率域滤波
1.基础卷积定理折叠误差补零当处理DFT时,图像及其变换是周期的.在周期接近函数非零部分的持续周期时,对周期函数进行卷积会引起相邻周期的串扰. 这种称为折叠误差的串扰可通过补零方法来避免. ...
图像处理中的傅里叶变换和频率域滤波概念
写在前面的话作者是一名在读的硕士研究僧,方向是机器视觉.由于视觉是一门相对复杂的学科,作者在课堂上学到的东西只是非常浅显的内容,我们老师说是,领我们进了个门.现在打算利用图书馆和网络上的资源进行自学 ...
数字图像处理-频率域滤波原理
from:https://blog.csdn.net/forrest02/article/details/55510711?locationNum=15&fps=1 写在前面的话作者是一名在 ...
OpenCV —— 频率域滤波（傅里叶变换，低通和高通滤波，带通和带阻滤波，同态滤波）
频率域滤波基本概念傅里叶变换二维离散的傅里叶变换快速傅里叶变换傅里叶幅度谱与相位谱谱残差显著性检测卷积与傅里叶变换的频率域滤波低通滤波和高通滤波带通和带阻滤波同态滤波基本概念 ...
OpenCV实现频率域滤波——以高斯低通滤波去噪为例
最近由于作业原因,试着用OpenCV实现频率域滤波,但是OpenCV中并没有像MATLAB中fftshift这样的中心化操作,所以我写了一个频率域滤波的函数,以后用频率域滤波的时候在主函数中调用即可. ...
数字图像处理-频率域滤波
一.实验原理频率域滤波是对图像进行傅里叶变换,将图像由图像空间转换到频域空间,然后在频率域中对图像的频谱作分析处理,以改变图像的频率特征,原理是用傅里叶变换表示的函数特征完全可以通过傅里叶反变换来重 ...

频率域滤波之带阻滤波器matlab代码

频率域滤波之带阻滤波器matlab代码相关推荐

最新文章

热门文章