图像高通和低通滤波处理.doc

《图像高通和低通滤波处理.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验五图像高通和低通滤波处理一、实验目的掌握常用频域高通滤波器和低通滤波器的设计。进一步加深理解和掌握图像频谱的特点和频域高通滤波及频域低通滤波的原理。理解图象高通滤波和低通滤波的处理过程和特点。二、实验要求1.实验之前要预习2.独立完成程序的编写3.写出实验报告4.实验每组1人实验设备：每组计算机一台三、实验原理低通滤波器：容许低频信号通过，但减弱（或减少）频率高於截止频率的信号的通过。用在绘制长期走势或均化。高通滤波器：容许高频信号通过、但减弱（或减少）频率低於截止频率信号通过的滤波器。强调细节。二维理想高通滤波器的传递

2、函数为：D0是从频率矩形中点测得的截止频率长度，它将以D0为半径的圆周内的所有频率置零，而毫不衰减地通过圆周外的任何频率。但其物理上是不可实现的。巴特沃斯高通滤波器的传递函数为：式中D0为截止频率距远点距离。与低通滤波器的情况一样，可认为巴特沃斯高通型滤波器比IHPF更平滑。高斯高通滤波器传递函数为：高通滤波器能够用高斯型低通滤波器的差构成。这些不同的滤波器有更多的参数，因此能够对滤波器的形状进行更多的控制。一、实验内容1.对数字图象进行低通滤波处理2.对数字图象进行高通滤波处理3.比较和分析所得到的结果。二、实验过程1.流

3、程图2.试验程序%----------------------------选择图片路径及显示---------------------------%globalf%选择图片路径[filename,pathname]=...%...表示与下行连接uigetfile({'*.bmp';'*.jpg';'*gif';'*tif'},'选择图片');%打开文件类型，对话框名称f=imread([pathnamefilename]);axes(handles.axes1);%在axes1显示原图像imshow(f);title('原始

4、图像');%---Executesonselectionchangeinpopupmenu1.functionpopupmenu1_Callback(hObject,eventdata,handles)%hObjecthandletopopupmenu1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%Hints:contents=get(hObj

5、ect,'String')returnspopupmenu1contentsascellarray%contents{get(hObject,'Value')}returnsselecteditemfrompopupmenu1globalfva=get(handles.popupmenu1,'Value');val=get(hObject,'Value');switchval%用switch语句设置选项case1%巴特沃斯高通滤波器g=fft2(f);%傅立叶变换g=fftshift(g);%转换数据矩阵[M,N]=size

6、(g);nn=2;%二阶Butterworth高通滤波器d0=20;%截止频率为10m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt((i-m)^2+(j-n)^2);h=1/(1+d0/(0.414*d)^(2*nn));%计算高通滤波器传递函数result(i,j)=h*g(i,j);endendresult=ifftshift(result);y2=ifft2(result);y3=uint8(real(y2));axes(handles.axes2);imshow(y3);%显

7、示滤波处理后的图像title('巴特沃斯高通滤波');case2%巴特沃斯低通滤波g=fft2(f);%傅立叶变换g=fftshift(g);%转换数据矩阵[M,N]=size(g);nn=2;%二阶巴特沃斯(Butterworth)低通滤波器d0=30;%截止频率为10m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt((i-m)^2+(j-n)^2);h=1/(1+0.414*(d/d0)^(2*nn));%计算低通滤波器传递函数result(i,j)=h*g(i,j);endend

8、result=ifftshift(result);y2=ifft2(result);y3=uint8(real(y2));axes(handles.axes3);imshow(y3);%显示滤波处理后的图像title('巴特沃斯低通滤波')case3%理想高通滤波k=fft2(f);%傅立