图像高通与低通滤波处理.docVIP

实验五 图像高通和低通滤波处理 实验目的 掌握常用频域高通滤波器和低通滤波器的设计。进一步加深理解和掌握图像频谱的特点和频域高通滤波及频域低通滤波的原理。理解图象高通滤波和低通滤波的处理过程和特点。 实验要求 1.实验之前要预习 2.独立完成程序的编写 3.写出实验报告 4.实验每组1人 实验设备：每组计算机一台 实验原理 低通滤波器：容许低频信号通过，但减弱（或减少）频率高於截止频率的信号的通过。用在绘制长期走势或均化。 高通滤波器：容许高频信号通过、但减弱（或减少）频率低於截止频率信号通过的滤波器。强调细节。 D0是从频率矩形中点测得的截止频率长度，它将以D0为半径的圆周内的所有频率置零，而毫不衰减地通过圆周外的任何频率。但其物理上是不可实现的。 巴特沃斯高通滤波器的传递函数为： 式中D0为截止频率距远点距离。与低通滤波器的情况一样，可认为巴特沃斯高通型滤波器比IHPF更平滑。 高斯高通滤波器传递函数为： 高通滤波器能够用高斯型低通滤波器的差构成。这些不同的滤波器有更多的参数，因此能够对滤波器的形状进行更多的控制。 实验内容 1.对数字图象进行低通滤波处理 2.对数字图象进行高通滤波处理 3.比较和分析所得到的结果。 实验过程 流程图 试验程序 %----------------------------选择图片路径及显示---------------------------% global f %选择图片路径 [filename,pathname]= ... %...表示与下行连接 uigetfile({*.bmp;*.jpg;*gif;*tif},选择图片); %打开文件类型，对话框名称 f=imread([pathname filename]); axes(handles.axes1); %在axes1显示原图像 imshow(f) ; title(原始图像); % --- Executes on selection change in popupmenu1. function popupmenu1_Callback(hObject, eventdata, handles) % hObject handle to popupmenu1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: contents = get(hObject,String) returns popupmenu1 contents as cell array % contents{get(hObject,Value)} returns selected item from popupmenu1 global f va=get(handles.popupmenu1,Value); val=get(hObject,Value); switch val %用switch语句设置选项 case 1 %巴特沃斯高通滤波器 g=fft2(f); % 傅立叶变换 g=fftshift(g); % 转换数据矩阵 [M,N]=size(g); nn=2; % 二阶 Butterworth高通滤波器 d0=20; %截止频率为10 m=fix(M/2); n=fix(N/2); for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2); h=1/(1+d0/(0.414*d)^(2*nn));% 计算高通滤波器传递函数 result(i,j)=h*g(i,j); end end result=ifftshift(result); y2=ifft2(result); y3=uint8(real(y2)); axes(handles.axes2); imshow(y3); % 显示滤波处理后的图像 title