作业帮 > 综合 > 作业

用matlab编写使用二维高斯低通滤波器的函数,使用大小为3*3,标准差分别为0.5,1.5,2.5

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/07/03 12:31:33
用matlab编写使用二维高斯低通滤波器的函数,使用大小为3*3,标准差分别为0.5,1.5,2.5
对灰度图像,真彩色图像,伪彩色图像进行滤波
高斯低通滤波对图像进行低通滤波代码,这个是灰度图像的
彩色图像就是取出RGB三个分量进行滤波再整合
I1=imread('flower.jpg'); %读取图像
I=rgb2gray(I1); %将图像变为灰度图象
%将灰度图像的二维不连续Frourier变换的零频率成分移到频谱的中心
s=fftshift(fft2(I));
[M,N]=size(s); %分别返回s的行数到M中,列数到N中
n=2; %对n赋初值
%GLPF滤波,d0=5,15,30(程序中以d0=30为例)
d0=30; %初始化d0
n1=floor(M/2); %对M/2进行取整
n2=floor(N/2); %对N/2进行取整
for i=1:M
for j=1:N
d=sqrt((i-n1)^2+(j-n2)^2); %点(i,j)到傅立叶变换中心的距离
h=1*exp(-1/2*(d^2/d0^2)); %GLPF滤波函数
s(i,j)=h*s(i,j); %GLPF滤波后的频域表示
end
end
s=ifftshift(s); %对s进行反FFT移动
%对s进行二维反离散的Fourier变换后,取复数的实部转化为无符号8位整数
s=uint8(real(ifft2(s)));
figure; %创建图形图像对象
imshow(s); %显示GLPF滤波处理后的图像
title('GLPF滤波(d0=30)'); %为经GLPF滤波后的图像添加标题