实践一理想低通滤波器、Butterworth低通滤波器、高斯低通
滤波器
111理想低通滤波器实践代码：
IimreadcouplebmpIrgb2grayIsubplot221imshowItitle原图像sfftshiftfft2Isubplot223imshowabsstitle图像傅里叶变换所得频谱subplot224imshowlogabsstitle图像傅里叶变换取对数所得频谱
absizesa0rou
da2b0rou
db2d10fori1aforj1b
dista
cesqrtia02jb02ifdista
cedh1elseh0e
dsijhsije
de
dsui
t8realifft2ifftshiftssubplot222imshowstitle低通滤波所得图像
IimreadcouplebmpHdo
essizeIHdr020figuresurfHdFacecolori
terpEdgecolor
o
eFacelighti
gpho
g画三维曲面（色）图
112理想低通滤波器实践结果截图：
f121Butterworth低通滤波器实践代码：
clearallI1imreadgirlbmpsubplot221imshowI1xlabela原始图像fdoubleI1数据类型转换gfft2f图像傅里叶转换gfftshiftg傅里叶变换平移
fF2logabsg对傅里叶变换结果取绝对值，然后取对数
subplot222imshowF2I
itialMag
ificatio
fit将计算后的矩阵用图像表
示
colormapjet设置色彩索引图
colorbar
显示色彩索引条
xlabelb原始图像的傅里叶变换图像
N1N2sizeg傅里叶变换图像尺寸

2参数赋初始值
d05

1fixN12数据圆整

2fixN22数据圆整
fori1N1遍历图像像素
forj1N2
dsqrti
12j
22
ifd0
h0
else
h11dd02

e
d
resultijhgij图像矩阵计算处理
e
d
e
d
F3logabsresult对傅里叶变换结果取绝对值，然后取对数
subplot223imshowF3I
itialMag
ificatio
fit
colormapjet设置色彩索引图
colorbar
显示色彩索引条
xlabelc滤波后的傅里叶变换图像
resultifftshiftresult
X2ifft2result
X3ui
t8realX2
subplot224imshowX3
xlabeldButterworth低通滤波图像
I1imreadcouplebmpf1f2freqspacesizeI1meshgridD03rf12f22
4fori1sizeI11
forj1sizeI12trijDDHdij1t
1
e
de
d
fsurfHdFacecolori
terpEdgecolor
o
eFacelighti
gpho
g画三维曲面（色）图
122Butterworth低通滤波器实践结果截图：
131高斯低通滤波器实践代码：
IAimreadgirlbmpf1f2freqspacesizeIAmeshgridD100sizeIA1rf12f22Hdo
essizeIA
ffori1sizeIA1forj1sizeIA2trijDDHdijexpte
d
e
dYfft2doubleIAYfftshiftYYaYHdYaifftshiftYaIarealifft2Yafiguresubplot221imshowui
t8IAtitle原图像subplot222imshowui
t8Iatitle高斯低通滤波处理figuresurfHdFacecolori
terpEdgecolor
o
eFacelighti
gpho
g画三维曲r