低通滤波.rar_harderj9a_quite8is_低通滤波_图像低通滤波

close all;clear all;clc; %利用理想低通滤波器对图像进行滤波 I=imread('coins.png'); I=im2double(I); M=2*size(I,1); %滤波器的行数 N=2*size(I,2); %滤波器的列数 u=-M/2:(M/2-1); v=-N/2:(N/2-1); [U,V]=meshgrid(u,v); D=sqrt(U.^2+V.^2); D0=80; %截止频率设置为80 H=double(D<=D0); %理想低通滤波器 J=fftshift(fft2(I,size(H,1),size(H,2))); %时域图像转换到频域 K=J.*H; %滤波处理 L=ifft2(ifftshift(K)); %傅立叶反变换 L=L(1:size(I,1),1:size(I,2)); figure; subplot(121);imshow(I); %显示原图像 subplot(122),imshow(L); %显示滤波后的图像 --------------------- 作者：Ibelievesunshine 来源：CSDN 原文：https://blog.csdn.net/ibelievesunshine/article/details/79965821 版权声明：本文为博主原创文章，转载请附上博文链接！