数字图像处理——陷波滤波算法

陷波滤波

1.什么是陷波滤波？

滤波器种类：低通、高通、带通、带阻、全通。其中的带阻滤波器就是指能通过大多数频率分量、但将某些范围的频率分量衰减到极低水平的滤波器，而陷波滤波器就是带阻滤波器的特殊情况，即不再是“某些范围的频率分量”而是指“阻带范围极小的频率分量”。

2.陷波带阻滤波算法

n阶布特沃斯陷波带阻滤波器表达式：

零相移滤波器必须是关于原点对称的，因此，一个中心位于(u0,v0)的陷波在位置(-u0,-v0)必须有一个对应的陷波

3.陷波带通滤波算法

n阶布特沃斯陷波带通滤波器表达式:

4.陷波滤波算法Matlab实现

镇楼图：

Matlab代码：

close all;
clear all;
clc;
%------------------- Input ---------------------------%
%读取图像
I = imread('.\Images\test1.tif');

%------------------- Image Processing ----------------%
%图像填充
[M,N] = size(I);
P = 2*M;
Q = 2*N;
fillimage = uint8(zeros(P,Q));
fillimage(1:M,1:N) = I;

%频域中心化
for x=1:P
    for y=1:Q
        h(x,y)=(-1)^(x+y);
    end
end
fillimagecenter = h.*double