Tabla de contenido
Principio de filtro de paso alto
3. Filtro de paso alto Butterworth
Código (filtro de paso alto de Butterworth)
4. Filtro de paso alto gaussiano
1. Introducción
Principio de filtro de paso alto:
Atenúa o rechaza los componentes de baja frecuencia, permitiendo que pasen los componentes de alta frecuencia
efecto:
Haga que la imagen sea más nítida y resalte el borde de la imagen.
(Nota: en términos generales, el filtrado de paso alto no tiene efecto en la supresión de ruido . Si simplemente se utiliza el filtrado de paso alto, es posible que sea difícil lograr una mejora satisfactoria en la calidad de la imagen debido al ruido grave. Para mejorar los detalles de la imagen y suprimir el ruido , filtrado de mejora de alta frecuencia , que en realidad se compone de un filtro de paso alto y un filtro de paso total, de modo que la información de baja frecuencia se puede retener sobre la base del filtrado de paso alto)
2. Filtro de paso alto ideal
Generar fórmula:
(donde D0 es la frecuencia de corte del filtro de paso alto ideal)
3. Filtro de paso alto Butterworth
Generar fórmula:
(donde D0 es la frecuencia de corte del filtro de paso alto Butterworth, n es el orden del filtro Butterworth, que se utiliza para controlar la inclinación del filtro)
Código (filtro de paso alto Butterworth):
clear all;clc;
I = imread('D:\resource_photo\1(1).png');
figure(1);
subplot(221),imshow(I);
title('原始图像');
k=fft2(im2double(I));
s=fftshift(fft2(im2double(I)));
subplot(222),imshow(k);
title('二维傅里叶变换后图像');
subplot(223),imshow(log(abs(s)));
title('中心频谱图像')
[a,b]=size(s);
a0=round(a/2);
b0=round(b/2);
d0=50; % 将理想高通滤波器的截止频率D0设置为50
for i=1:a %双重for循环计算频率点(i,j)与频域中心的距离D(i,j)=sqrt((i-round(a/2)^2+(j-round(b/2)^2))
for j=1:b
distance=sqrt((i-a0)^2+(j-b0)^2);
if distance<=d0 % 根据理想高通滤波器产生公式,当D(i,j)<=D0,置为0
h=0;
else
h=1; % 根据理想高通滤波器产生公式,当D(i,j)>D0,置为1
end
s(i,j)=h*s(i,j);% 频域图像乘以滤波器的系数
end
end
% real函数取元素的实部
s=real(ifft2(ifftshift(s)));% 最后进行二维傅里叶反变换转换为时域图像
subplot(224),imshow(s,[]);
title('理想高通滤波所得图像');
resultado de la operación:
4. Filtro de paso alto gaussiano
Generar fórmula:
(donde D0 es la frecuencia de corte del filtro de paso alto gaussiano)
