基于某某matlab对图像进行高通、低通、带通滤波
- 格式:doc
- 大小:899.78 KB
- 文档页数:27
matlab自带的滤波器函数
Matlab自带的滤波器函数可以用于对信号进行滤波处理,常用的函数有:
1. fir1函数:设计一阶低通、高通、带通、带阻滤波器的FIR 数字滤波器,可自定义通带和阻带的截止频率。
2. cheby1函数:设计ChebyshevI型低通、高通、带通、带阻数字滤波器,可自定义通带和阻带的截止频率和最大通带波纹。
3. butter函数:设计Butterworth型低通、高通、带通、带阻数字滤波器,可自定义通带和阻带的截止频率和滤波器阶数。
4. filtfilt函数:对信号进行双向滤波处理,可避免滤波后信号的相位畸变和滞后。
这些函数可以在Matlab的Signal Processing Toolbox中找到,可根据需要选择合适的函数进行滤波处理。
- 1 -。
在MATLAB 中,你可以使用内置的butter函数来设计带通滤波器。
以下是一个示例,展示了如何创建一个带通滤波器,并将其应用于一个信号:生成一个带噪声的信号
Fs = 1000; 采样频率
t = 0:1/Fs:1-1/Fs; 时间向量
x = sin(2*pi*50*t) + sin(2*pi*120*t); 合成信号,包含50Hz和120Hz 的成分
xn = x + 2*randn(size(t)); 添加高斯白噪声
设计带通滤波器
[b,a] = butter(4, [50 120]/Fs, 'bandpass'); 'bandpass' 表示带通滤波器应用滤波器
y = filter(b, a, xn);
绘制原始信号和滤波后的信号
subplot(2,1,1);
plot(t, xn);
title('原始信号');
subplot(2,1,2);
plot(t, y);
title('滤波后的信号');
在这个示例中,我们首先生成了一个包含两个频率成分(50Hz和120Hz)的合成信号,并添加了高斯白噪声。
然后,我们使用butter函数设计了一个带通滤波器,该滤波器的通带范围为50Hz到120Hz。
最后,我们使用filter函数将滤波器应用于信号,并绘制了原始信号和滤波后的信号。
利用MATLAB设计巴特沃斯低通数字滤波器引言数字滤波器是数字信号处理中的重要组成部分,可以用于去除信号中的噪音和不需要的频率成分。
巴特沃斯滤波器是一种常见的数字滤波器,被广泛应用于信号处理领域。
本文将介绍如何利用MATLAB设计巴特沃斯低通数字滤波器,并给出详细的步骤和示例代码。
设计步骤利用MATLAB设计巴特沃斯低通数字滤波器主要包括以下步骤:1.设计滤波器的参数2.计算滤波器的传递函数3.绘制滤波器的幅频响应曲线4.通过频域图像观察滤波器的性能下面将分别介绍每个步骤的详细操作。
设计滤波器的参数巴特沃斯低通数字滤波器的参数包括截止频率和阶数。
截止频率决定了滤波器的通频带,阶数决定了滤波器的陡峭程度。
通过MATLAB的butter()函数可以方便地设计巴特沃斯低通数字滤波器。
该函数的参数为滤波器的阶数和截止频率。
示例代码如下:order = 4; % 阶数cutoff_freq = 0.4; % 截止频率[b, a] = butter(order, cutoff_freq);计算滤波器的传递函数通过设计参数计算得到滤波器的传递函数。
传递函数是一个复数,包括了滤波器的频率响应信息。
使用MATLAB的freqz()函数可以计算滤波器的传递函数。
该函数的参数为滤波器的系数b和a,以及频率取样点的数量。
示例代码如下:freq_points = 512; % 频率取样点数量[h, w] = freqz(b, a, freq_points);绘制滤波器的幅频响应曲线经过计算得到的传递函数能够提供滤波器的幅频响应信息。
通过绘制幅频响应曲线,可以直观地观察滤波器的频率特性。
使用MATLAB的plot()函数可以绘制滤波器的幅频响应曲线。
该函数的参数为频率点和传递函数的幅值。
示例代码如下:magnitude = abs(h); % 幅值plot(w/pi, magnitude);xlabel('归一化频率');ylabel('幅值');title('巴特沃斯低通数字滤波器幅频响应');通过频域图像观察滤波器的性能通过绘制滤波器的频域图像,可以直观地观察滤波器对不同频率的信号的响应情况。
利用MATLAB对彩色图像进行低通滤波
王潇;李晓光;李金成
【期刊名称】《电脑知识与技术》
【年(卷),期】2006(000)003
【摘要】彩色图像的数据量比较大,在进行压缩与传输时会增加工作负担.通过对图像的频谱特性分析可以知道,一幅图像中的低频成分要远远高于高频成分,大部分图像信息集中在低频部分.因此,在不过于苛刻地计较图像的再现质量的场合下,如果将代表图像边缘,跳越部分,以及颗粒噪声的高频部分通过低通滤波器滤除的话就会减少图像数据量,进而减少下一步的工作量.本文主要以真彩色RGB图像为例,介绍一种对其进行低通滤波的方法.
【总页数】2页(P165-165,184)
【作者】王潇;李晓光;李金成
【作者单位】北京交通大学电子信息工程学院,北京,100044;北京交通大学电子信息工程学院,北京,100044;北京交通大学电子信息工程学院,北京,100044
【正文语种】中文
【中图分类】TP391
【相关文献】
1.利用Matlab软件对医学彩色图像的空间转换 [J], 张兰凤;郭丹
2.基于低通滤波的彩色图像多模式自适应系统 [J], 杨博;周孝宽
3.利用MATLAB对彩色图像进行低通滤波 [J], 王潇
4.利用小波和分层聚类进行彩色图像分割 [J], 叶吉祥;陈香华;谭冠政
5.基于同态高低通滤波与多尺度Retinex的低照度彩色图像增强 [J], 张江鑫;杨惠因版权原因,仅展示原文概要,查看原文内容请购买。
计算机视觉(⼆)-matlab之理想低通滤波器,布特沃斯低通、⾼斯低通,理想⾼通、布特沃斯。
未整理完!在滤波器之前,先讲解傅⾥叶变换理想低通滤波器f = imread('Fig0441.tif');f = im2double(f);% 计算填充图像⼤⼩[M,N] = size(f);M2 = 2*M;N2 = 2*N;% 傅⾥叶变换F = fftshift(fft2(f,M2,N2));figure;imshow(mat2gray(log(1+abs(F))));title('傅⾥叶频谱');% 设计滤波器% ⽣成⽹格坐标u = -N:N-1;v = -M:M-1;[U,V] = meshgrid(u,v);% 设计滤波器D = hypot(U,V);D0 = 30; % 截⽌频率[10 30 60 160 460]H = mat2gray( D <= D0 );% 理想低通滤波器figure;imshow(H);title('理想低通滤波器');% 频域滤波G = F.*H;figure;imshow(mat2gray(log(1+abs(G))));title('频域滤波'); 1-sum(sum(abs(G).^2))/sum(sum(abs(F).^2))g0 = ifft2(fftshift(G));g = g0(1:M,1:N);g = real(g);figure;imshow(g);title('滤波后的图像');figure;subplot(121);imshow(f);title('原图');subplot(122);imshow(g);title('滤波后的图像');结果:原图:布特沃斯低通滤波器f = imread('Fig0441.tif');f = im2double(f);% 计算填充图像⼤⼩[M,N] = size(f);M2 = 2*M;N2 = 2*N;% 傅⾥叶变换F = fftshift(fft2(f,M2,N2));figure;imshow(mat2gray(log(1+abs(F))));title('傅⾥叶频谱'); % 设计滤波器% ⽣成⽹格坐标u = -N:N-1;v = -M:M-1;[U,V] = meshgrid(u,v);% 设计滤波器D = hypot(U,V);D0 = 30; % 截⽌频率[10 30 60 160 460]H = mat2gray(1./(1+((D./D0).^4)));figure;imshow(H);title('布特沃斯低通滤波器(n=2)');% 频域滤波G = F.*H;figure;imshow(mat2gray(log(1+abs(G))));title('频域滤波'); 1-sum(sum(abs(G).^2))/sum(sum(abs(F).^2))g0 = ifft2(fftshift(G));g = g0(1:M,1:N);g = real(g);figure;imshow(g);title('滤波后的图像');figure;subplot(121);imshow(f);title('原图');subplot(122);imshow(g);title('滤波后的图像'); matlab结果:⾼斯低通滤波器⾼斯低通滤波器代码:H = exp((-D.^2)/(2*(D0).^2));figure;imshow(H);title('滤波器(D0=100)');结果:理想⾼通滤波器理想⾼通滤波器代码:H = mat2gray( D >= D0 );%理想⾼通滤波器figure;imshow(H);title('滤波器()');结果:布特沃斯⾼通滤波器布特沃斯⾼通滤波器代码:H = 1./(1+((D0./D).^4));%布特沃斯⾼通滤波器figure;imshow(H);title('滤波器(n=2)');结果:⾼斯⾼通滤波器⾼斯⾼通滤波器代码:H = 1-exp((-D.^2)/(2*(D0).^2));%⾼斯⾼通滤波器figure;imshow(H);title('滤波器');matlab结果:混合⾼频强调滤波[file,path] = uigetfile({'*.png';'*.jpg';},'选择图⽚');f = imread([path,file]);f = im2double(f);%计算填充图像⼤⼩[M,N] = size(f);M2 = 2*M;N2 = 2*N;% 傅⾥叶变换F = fftshift(fft2(f,M2,N2));figure;imshow(mat2gray(log(1+abs(F))));title('傅⾥叶频谱'); % 设计滤波器% ⽣成⽹格坐标u = -N:N-1;v = -M:M-1;[U,V] = meshgrid(u,v);% 设计滤波器D = hypot(U,V);D0 = 40; %截⾄频率[10 30 60 160 460]H = 1-exp((-D.^2)/(2*(D0).^2));%¸⾼斯⾼通滤波G = (F.^(-1)).*((0.5+0.75.*H).*F);%figure;imshow(H);title('滤波器');figure;imshow(mat2gray(log(1+abs(G))));title('频率滤波'); % 傅⾥叶逆变换g0 = ifft2(fftshift(G));g = g0(1:M,1:N);g = real(g);i=histeq(g);figure;imshow(g);title('滤波后的图像');figure;imshow(i);title('直⽅图均衡');结果:原图:reference:李卫军,肖宛昂,董肖莉,覃鸿⽼师《视觉信息处理及FPGA 实现》课程等越是憧憬,越要风⾬兼程。
基于MATLAB的数字图像处理的设计与实现摘要数字图像处理是一门新兴技术,随着计算机硬件的发展,数字图像的实时处理已经成可能,由于数字图像处理的各种算法的出现,使得其处理速度越来越快,能更好的为人们服务。
数字图像处理是一种通过计算机采用一定的算法对图形图像进行处理的技术。
目的:改善医学图像质量,使图像得到增强。
方法:利用Matlab工具箱函数,采用灰度直方图均衡化和高通滤波的方法对一幅X线图像进行增强处理。
结果:用直方图均衡化的算法,将原始图像密集的灰度分布变得比较稀疏,处理后的图像视觉效果得以改善。
高通滤波对于局部细节增强显著,高通滤波后使不易观察到的细节变得清晰。
结论:使用Matlab工具箱大大简化了编程工作,为医学图像处理提供了一种技术平台。
经过直方图均衡化和高通滤波处理后的医学图像,视觉效果得到改善。
关键词:MATLAB;直方图均衡化;高通滤波;图像增强AbstractDigital image processing is an emerging technology, with the development of computer hardware, real—time digital image processing has become possible due to digital image processing algorithms to appear,making it faster and faster processing speed,better for people services .Digital image processing is used by some algorithms computer graphics image pro cessing technology. Objective:To improve the quality of medical image by enhancing the details。
matlab频域低通滤波Matlab频域低通滤波的过程包括以下几个步骤:准备工作、载入信号、转换信号、频域滤波、逆傅里叶变换、结果分析和可视化。
1. 准备工作:在进行频域低通滤波之前,需要先了解滤波的基本概念和频域变换的原理。
频域滤波是一种在频域中操作信号的方法,它能够消除或削弱不需要的频率成分,实现对信号的有选择性地处理。
常见的频域变换方法包括傅里叶变换和离散傅里叶变换。
2. 载入信号:使用Matlab自带的`audioread`函数或者其他适用的载入函数,将待处理的音频文件读入到工作环境中。
例如,我们可以读取一个.wav格式的音频文件并将其存储为一个向量。
matlab[x, Fs] = audioread('audio.wav');3. 转换信号:由于频域滤波需要将信号转换为频域表示,因此需要对信号进行频域变换。
在Matlab中,可以使用`fft`函数对信号进行快速傅里叶变换(FFT)。
matlabX = fft(x);4. 频域滤波:根据滤波的需求,选择合适的滤波器类型进行滤波。
常用的低通滤波器有理想低通滤波器、巴特沃斯低通滤波器和高斯低通滤波器。
这里以理想低通滤波器为例。
理想低通滤波器的特点是在截止频率之前保留信号的全部频谱成分,并在截止频率之后完全削弱信号的高频成分。
使用以下代码可以生成一个理想低通滤波器的频域响应:matlabN = length(X);cutoff_freq = 1000; 设定截止频率H = zeros(N, 1);H(1:cutoff_freq) = 1;H(N-cutoff_freq+2:N) = 1;5. 滤波:将信号的频谱与滤波器的频谱相乘,可以实现滤波效果。
matlabY = X .* H;6. 逆傅里叶变换:将频域滤波后的信号进行逆变换,转换回时域表示。
在Matlab中,可以使用`ifft`函数对信号进行逆傅里叶变换。
matlaby = ifft(Y);7. 结果分析和可视化:经过频域低通滤波后得到的信号,可以通过时域分析和频谱分析来验证和评估滤波效果。
低通滤波matlab
低通滤波是指只允许低于某一频率范围的信号通过的滤波器,而阻止高于该频率范围的信号通过。
在MATLAB中,可以使用以下函数来实现低通滤波:
1. butter()函数:用于设计巴特沃斯滤波器,可设置滤波器的阶数、截止频率等参数。
2. fir1()函数:用于设计FIR(有限脉冲响应)滤波器,可设置滤波器的截止频率、窗函数等参数。
3. filter()函数:用于应用滤波器,输入原始信号和滤波器的系数,输出滤波后的信号。
例如,以下代码可以实现对信号x进行低通滤波,保留频率在0-1000Hz范围内的信号:
matlab
fs = 10000; %采样频率
fc = 1000; %截止频率
Wn = fc/(fs/2); %归一化截止频率
[b,a] = butter(6,Wn,'low'); %设计6阶巴特沃斯滤波器
y = filter(b,a,x); %应用滤波器
其中,b和a分别代表巴特沃斯滤波器的Numerator和Denominator系数,6表示滤波器的阶数,x是原始信号,y是滤波后的信号。
需要注意的是,滤波器的阶数越高,滤波效果通常越好,但也会增加计算量和延迟时间。
matlab巴特沃斯低通滤波函数巴特沃斯低通滤波是一种常用的数字滤波方法,用于滤除信号中高频成分,保留低频信号。
在MATLAB中,可以通过调用内置函数`butter`来实现巴特沃斯低通滤波。
本文将介绍这个函数的使用方法,并给出一个简单的示例。
## 巴特沃斯低通滤波概述巴特沃斯低通滤波器是一种无限脉冲响应(IIR)滤波器,它的频率响应曲线为-3dB截止频率。
在信号处理中,巴特沃斯低通滤波器可以用于去除信号中的高频噪声,使得信号更加平滑。
## MATLAB中的`butter`函数MATLAB提供了`butter`函数来设计巴特沃斯滤波器。
该函数的基本语法格式如下:```[b, a] = butter(n, Wn, 'type')```其中,`n`是滤波器阶数,`Wn`为归一化的截止频率,范围为[0, 1],具体取值要根据实际需求设定。
参数`'type'`用来指定滤波器类型,可选的值有'low'、'high'、'bandpass'、'stop'。
`butter`函数将返回滤波器的分子系数`b`和分母系数`a`,用于滤波器的差分方程表达式。
## 示例假设我们有一个包含噪声的信号`x`,现在希望设计一个巴特沃斯低通滤波器来去除信号中的高频成分。
下面是具体的步骤:1. 生成一个包含噪声的信号`x`,作为示例输入信号。
```matlabFs = 1000; % 采样率t = 0:1/Fs:1; % 时间长度x = sin(2*pi*50*t) + sin(2*pi*120*t) + randn(size(t));```2. 设计巴特沃斯低通滤波器。
```matlabfc = 150; % 截止频率Wn = fc/(Fs/2); % 归一化的截止频率n = 6; % 滤波器阶数[b, a] = butter(n, Wn, 'low');```3. 使用设计好的巴特沃斯低通滤波器对信号进行滤波。
【导言】1. 序Matlab 是一种强大的数学软件,广泛应用于科学计算、数据分析和图像处理等领域。
其丰富的函数库和灵活的编程环境使得它成为许多研究人员和工程师的首选工具之一。
2. 研究背景数据低通滤波是一种常见的信号处理技术,用于去除信号中的高频噪声,平滑信号曲线,提取信号的潜在趋势。
在实际工程和科学研究中,低通滤波常常被用于处理传感器数据、音频信号、图像等各种类型的信号。
3. 目的本文旨在介绍基于 Matlab 的数据低通滤波算法,包括算法原理、实现步骤和应用范例,帮助读者了解该算法的基本原理和实际应用,同时通过具体的代码示例和实验结果来验证算法的有效性。
【算法原理】1. 信号与频率信号可以分解为不同频率的分量,高频分量对应着信号的快速变化部分,而低频分量对应着信号的缓慢变化部分。
低通滤波就是通过滤波器去除信号中的高频分量,保留低频分量,从而平滑信号。
2. 离散时间信号的滤波在数字信号处理中,通常采用差分方程表示滤波器的行为。
对于离散时间信号,可以使用差分方程描述数字滤波器的输入输出关系,其中包括滤波器的系数和延迟项。
3. Matlab 中的滤波器设计工具Matlab 提供了丰富的滤波器设计工具,包括基于频率响应的滤波器设计、基于窗函数的滤波器设计、基于优化算法的滤波器设计等多种方法。
用户可以根据具体的需求选择合适的滤波器设计方法。
4. 低通滤波器的设计低通滤波器通常具有截止频率,截止频率之上的信号被滤除,而截止频率之下的信号被保留。
在 Matlab 中,可以通过设计滤波器的频率响应来实现低通滤波器的设计。
【算法实现】1. 滤波器设计使用 Matlab 提供的滤波器设计工具,根据具体的要求设计低通滤波器。
用户可以设定滤波器的截止频率、滤波器类型(巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器等)和滤波器阶数等参数。
2. 滤波器与信号处理将设计好的低通滤波器应用到信号处理中,可以使用 Matlab 中的滤波函数对信号进行滤波处理。
标准实用文案 文档 数字图像处理三级项目
—高通、低通、带通滤波器 标准实用文案
文档 摘要 在图像处理的过程中,消除图像的噪声干扰是一个非常重要的问题。利用matlab软件,采用频域滤波的方式,对图像进行低通和高通滤波处理。低通滤波是要保留图像中的低频分量而除去高频分量,由于图像中的边缘和噪声都对应图像傅里叶频谱中的高频部分,所以低通滤波可以除去或消弱噪声的影响并模糊边缘轮廓;高通滤波是要保留图像中的高频分量而除去低频分量,所以高通滤波可以保留较多的边缘轮廓信息。低通滤波器有巴特沃斯滤波器和高斯滤波器等等,本次设计使用的低通滤波器为****。高通滤波器有巴特沃斯滤波器、高斯滤波器、Laplacian高通滤波器以及Unmask高通滤波器等等,本次设计使用巴特沃斯高通滤波器。
1、 频域低通滤波器:设计低通滤波器包括 butterworth and Gaussian (选择标准实用文案
文档 合适的半径,计算功率谱比),平滑测试图像test1和2。
实验原理分析 根据卷积定理,两个空间函数的卷积可以通过计算两个傅立叶变换函数的乘积的逆变换得到,如果f(x, y)和h(x, y)分别代表图像与空间滤波器,F(u, v)和H(u, v)分别为响应的傅立叶变换(H(u, v)又称为传递函数),那么我们可以利用卷积定理来进行频域滤波。
在频域空间,图像的信息表现为不同频率分量的组合。如果能让某个范围内的分量或某些频率的分量受到抑制,而让其他分量不受影响,就可以改变输出图的频率分布,达到不同的增强目的。
频域空间的增强方法的步骤: (1)将图像从图像空间转换到频域空间; (2)在频域空间对图像进行增强; (3)将增强后的图像再从频域空间转换到图像空间。 低通滤波是要保留图像中的低频分量而除去高频分量。图像中的边缘和噪声都对应图像傅里叶频谱中的高频部分,所以低通滤波可以除去或消弱噪声的影响并模糊边缘轮廓。理想低通滤波器具有传递函数: 标准实用文案
文档 其中D0为制定的非负数,D(u,v)为点(u,v)到滤波器中心的距离。
功率谱比的定义: g(,)(,)fPuvLPuv
其中,(,)
fPuv为滤波前图像的功率谱,g(,)Puv为滤波后图像的功率谱。
频率计算公式为:2(,)(,)fPuvFuv,2(,)(,)gPuvGuv。
① Butterworth滤波器设计: 理想低通滤波器在数学上定义得很清楚,在计算机模拟中也可实现,但在截断频率处直上直下的理想低通滤波器是不能用实际的电子器件实现的。 n阶Butterworth低通滤波器(BLPF)的传递函数(截止频率距原点的
距离为 )定义如下:
(1) 其中,22(uv)=+v-N/2D,(u-M/2)()。 (2)
不同于ILPF,BLPF变换函数在通带与被滤除的频率之间没有明显的截断。对于有平滑传递函数的滤波器,定义一个截止频率的位置并使H(u,v)幅度降到其最大值的一部分。在式(1)中,当D(u,v)=D0时,H(u,v)=0.5(从最大值降到它的50%)。 一阶的巴特沃斯滤波器没有振铃,在二阶中振铃通常很微小,这是因为与理想低通滤波器相比,它的通带与阻带之间没有明显的跳跃,在高低频率间的过渡
nDvuDvuH2
0/),(11),(
0D标准实用文案
文档 比较光滑。巴特沃斯低通滤波器的处理结果比理想滤波器的要好,但阶数增高时振铃便成为一个重要因素。本次实验中设计实现了二阶巴特沃斯滤波器。 根据以上原理设计Butterworth低通滤波器,其处理结果如下图示: 源图像test1.pgmButterworth低通滤波图test1.pgm
源图像test2.tifButterworth低通滤波图test2.tif 理想低通滤波器有明显的振铃现象,而巴特沃斯滤波器的效果较好。 计算得test1的功率谱比L=0.9939。test2的功率谱比为0.9902。 ②Gaussian滤波器设计:
二维高斯低通滤波器,其传递函数的形式为: 22(,)/2(uv)=1-eDuvH, (3)
其中,22(uv)=+v-N/2D,(u-M/2)()。表示高斯曲线扩展的程度。使=D0,
可以将滤波器表示为: 标准实用文案 文档 220(,)/2(uv)=eDuvDH, (4)
其中,D0是截止频率。当D(u,v)=D0时,滤波器下降到它最大值的0.607倍处。 由于高斯低通滤波器的傅里叶反变换也是高斯的,这就是说通过式(3)或式(4)的傅里叶反变换而得到的空间高斯滤波器将没有振铃。
根据以上分析,设计Gaussian低通滤波器,处理结果如下: 源图像test1.pgmGaussian低通滤波后的test1.pgm(r=5)
Gaussian低通滤波后的test1.pgm(r=15)Gaussian低通滤波后的test1.pgm(r=80)标准实用文案
文档 Gaussian低通滤波后的test1.pgm(r=80)Gaussian低通滤波后的test1.pgm(r=230)
源图像test2.tifGaussian低通滤波后的test2.tif(r=5) Gaussian低通滤波后的test2.tif(r=15)Gaussian低通滤波后的test2.tif(r=30) 标准实用文案
文档 Gaussian低通滤波后的test2.tif(r=80)Gaussian低通滤波后的test2.tif(r=230)
可见,当滤波器的半径不同时,对应的滤波效果也不同。半径越小,平滑效果越明显,但半径过小,会使得图像变得模糊不清。 计算得test1(r=5)的功率谱比L= 0.4674。test2(r=5)的功率谱比为L=0.2930。
2、 频域高通滤波器:设计高通滤波器包括butterworth and Gaussian,在频域增强边缘。选择半径和计算功率谱比,测试图像test3,4:
实验原理分析 高通滤波是要保留图像中的高频分量而除去低频分量。理想高通滤波器传递函数表示为:
00),(1),(0),(DvuDDvuDvuH如如标准实用文案 文档 ① Butterworth滤波器设计: n阶Butterworth高通滤波器(BLPF)的传递函数(截止频率距原点的
距离为 )定义如下:
(5) 其中, 22(uv)=+v-N/2D,(u-M/2)()。 (6)
不同于ILPF,BLPF变换函数在通带与被滤除的频率之间没有明显的截断。对于有平滑传递函数的滤波器,定义一个截止频率的位置并使H(u,v)幅度降到其最大值的一部分。在式(1)中,当D(u,v)=D0时,H(u,v)=0.5(从最大值降到它的50%)。 根据以上原理设计Butterworth高通滤波器,其处理结果如下图示: 源图像test3corrupt.pgmButterworth高通滤波图test3corrupt.pgm
0DnvuDDvuH2
0),(11),(
标准实用文案
文档 源图像test4 copy.bmpButterworth高通滤波图test4 copy.bmp
计算得test3的功率谱比为0.0851。test4的功率谱比为0.0547。 ②Gaussian滤波器设计: 二维高斯高通滤波器,其传递函数的形式为: 22(,)/2(uv)=1-eDuvH, (7)
其中,22(uv)=+v-N/2D,(u-M/2)()。表示高斯曲线扩展的程度。使=D0,
可以将滤波器表示为: 220(,)/2(uv)=1-eDuvDH, (8)
其中,D0是截止频率。当D(u,v)=D0时,滤波器下降到它最大值的0.607倍处。 由于高斯低通滤波器的傅里叶反变换也是高斯的,这就是说通过式(7)或式(8)的傅里叶反变换而得到的空间高斯滤波器将没有振铃。
根据以上分析,设计Gaussian高通滤波器,处理结果如下: 标准实用文案
文档 源图像test3corrupt.pgmGaussian高通滤波后的test3corrupt.pgm(r=5)
Gaussian高通滤波后的test3corrupt.pgm(r=15)Gaussian高通滤波后的test3corrupt.pgm(r=30) 源图像test4 copy.bmpGaussian高通滤波后的test4 copy.bmp(r=5)