(精品word)高通滤波器matlab程序代码-豆柴文库

您所在位置：网站首页 / (精品word)高通滤波器matlab程序代码.doc / 文档详情

在线预览结束，喜欢就下载吧，查找使用更方便

5 金币

下载文档

/ 2

下载提示

如果您无法下载资料，请参考说明：

1、部分资料下载需要金币，请确保您的账户上有足够的金币

2、已购买过的文档，再次下载不重复扣费

3、资料包下载后请先用软件解压，在使用对应软件打开

文本预览

(精品word)高通滤波器matlab程序代码
(精品word)高通滤波器matlab程序代码
(精品word)高通滤波器matlab程序代码
%高斯低通滤波器
RGB=imread（’132。jpg'）;
I0=rgb2gray（RGB)；
subplot(2,3,1），imshow（I0）;title（'原图’）；
I1=imnoise(I0，’gaussian’）；％对原图像加噪声
subplot（2，3,2)，imshow（I1);title（’加入噪声后’)
%将灰度图像的二维不连续Fourier变换的零频率成分移到频谱的中心
s=fftshift(fft2（I1)）;
subplot（2，3,3)，imshow（log(1+abs(s）），［］)；title（’fftshift')；
[M，N]=size（s）；%分别返回s的行数到M中，列数到N中
％GLPF滤波
d0=50；%初始化d0
n1=floor(M/2）;%对M/2进行取整
n2=floor(N/2);％对N/2进行取整
fori=1：M
forj=1：N
d=sqrt（(i-n1)^2+(j-n2）^2）;%点（i，j）到傅立叶变换中心的距离
h（i，j）=1*exp(—1/2＊(d^2/d0^2）);%GLPF滤波函数
s(i，j）=h（i，j)＊s(i，j）；％GLPF滤波后的频域表示
end
end
s=ifftshift(s）;%对s进行反FFT移动
％对s进行二维反离散的Fourier变换后，取复数的实部转化为无符号8位整数
s=uint8(real(ifft2(s))）;
subplot(2,3，4）,imshow(h)；title(’传递函数')；％显示GHPF滤波器的传递函数
subplot（2，3，5）,imshow(s）;title(’GLPF滤波（d0=50）')；%显示GLPF滤波处理后的图像