数字图像滤波器平滑

最新推荐文章于 2024-10-07 23:32:36 发布

今天要来电代码吗

最新推荐文章于 2024-10-07 23:32:36 发布

阅读量180

点赞数

分类专栏：数字图像处理文章标签：计算机视觉 opencv 人工智能

本文链接：https://blog.csdn.net/wzy1414/article/details/124989558

版权

数字图像处理专栏收录该内容

4 篇文章 0 订阅

订阅专栏


clear
close all
clc

I=imread('Fig5.jpg'); 
figure, imshow(I);

%实现中心平移
%f(x,y)*(-1)^(x+y)
%TODO：
%A(i,j) = A(i,j) * power(-1, i+j);
[m, n] = size(I);
for i=1:m
    for j=1:n
        I(i,j) = I(i,j) * power(-1, i+j);
    end
end
    
%将空间域图像转换到频率域
%二维离散傅里叶变换，并显示频谱图像
%TODO：
A = fft2(I);

figure,imshow(log(abs(A)),[]);

% 实现低通滤波器H
%TODO：
% D0=?
% (floor(m/2),floor(n/2))
% D: (y,x)->(floor(m/2),floor(n/2))
H = zeros(m,n);
D0=25;
for i=1:m
    for j=1:n
        D = (i-m/2).^2+(j-n/2).^2;
        D = sqrt(D);
        if(D<=D0)
            H(i,j) = 1;
        end            
    end
end

%频率域图像平滑处理
%TODO：
A = A.*H;

%将处理后的图像在空间域进行显示
A = ifft2(A);

for i=1:m
    for j=1:n
        A(i,j) = A(i,j) * power(-1, i+j);
    end
end

%对s进行二维反离散的Fourier变换后，取复数的实部转化为无符号8位整数
figure,imshow(uint8(real(A)));