收藏 分享(赏)

数字图像处理 降噪滤波 大作业.doc

上传人:精品资料 文档编号:10673789 上传时间:2019-12-19 格式:DOC 页数:27 大小:2.42MB
下载 相关 举报
数字图像处理 降噪滤波 大作业.doc_第1页
第1页 / 共27页
数字图像处理 降噪滤波 大作业.doc_第2页
第2页 / 共27页
数字图像处理 降噪滤波 大作业.doc_第3页
第3页 / 共27页
数字图像处理 降噪滤波 大作业.doc_第4页
第4页 / 共27页
数字图像处理 降噪滤波 大作业.doc_第5页
第5页 / 共27页
点击查看更多>>
资源描述

1、-1-昆明理工大学信息工程与自动化学院学生实验报告( 20152016 学年 第 一 学期 )课程名称:图形图像基础程序设计 开课实验室: 2015 年 12 月 1 日年级、专业、班物联网 131 学号 201310410133 姓名 李哲 成绩实验项目名称 图像综合处理 指导教师 毛存礼教师评语教师签名:年 月 日一、实验目的及内容目的:掌握和熟悉 Matlab 编程环境及语言;掌握图像降噪算法和用途。内容:在课程教学和查阅相关文献资料的基础上,选择下面一个数字图像处理技术专题,实现相应算法进行仿真实验,并完成大作业报告。专题如下:(1)图像增强处理技术;(2)图像降噪处理技术。2、题目分

2、析 利用 matlab 的 GUI 程序设计一个简单实用的图像处理程序。该程序应具备图像处理的常用功能,以满足用户的使用。现设计程序有以下基本功能:1)图像的读取和保存。2)通过自己输入数值,实现图像的旋转。3)图像直方图统计和直方图均衡,要求显示直方图统计,比较直方图均衡后的效果。4)能对图像加入各种噪声,5)并通过几种滤波算法实现去噪并显示结果。6)将图像转化成灰度图像。-2-3.总体设计软件的总体设计界面布局如上图所示分为显示区域与操作区域。上边为显示区域:显示载入原图,以及通过处理后的图像。操作区域:通过功能键实现对图像的各种处理。设计完成后运行的软件界面如下:-3-4、具体设计现介绍

3、各个功能模块的功能与实现。4.1 图像的读取和保存:(1)利用matlab中 “ uigetfile”、“imread” “imshow”实现图像文件的读取与显示:实现代码:function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)-4-% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user da

4、ta (see GUIDATA)filename,pathname=uigetfile(*.jpg;*.bmp;*.tif;*.*,载入图像 );if isequal(filename,0)|isequal(pathname,0)errordlg(没有选中文件 ,出错);return;else file=pathname,filename;global S %设置一个全局变量 S,保存初始图像路径,以便之后的还原操作S=file;x=imread(file);set(handles.axes1,HandleVisibility,ON);axes(handles.axes1);imshow(x)

5、;set(handles.axes1,HandleVisibility,OFF);axes(handles.axes2);imshow(x);handles.img=x;guidata(hObject,handles);end% - Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4 (see GCBO)% eventdata reserved - to be defined in a f

6、uture version of MATLAB% handles structure with handles and user data (see GUIDATA)sfilename ,sfilepath=uiputfile(*.jpg;*.bmp;*.tif;*.*,保存图像文件,untitled.jpg);if isequal(sfilename,sfilepath,0,0)sfilefullname=sfilepath ,sfilename;imwrite(handles.img,sfilefullname);elsemsgbox(你按了取消键, 保存失败 );end(2)图像保存。利

7、用uiputfile和imwrite函数实现图像文件的保存。-5-实现代码:function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)sfilename ,sfilepath=uiputfile(*.jpg

8、;*.bmp;*.tif;*.*,保存图像文件,untitled.jpg);if isequal(sfilename,sfilepath,0,0)sfilefullname=sfilepath ,sfilename;imwrite(handles.img,sfilefullname);elsemsgbox(你按了取消键, 保存失败 );end3)程序的退出。-6-实现代码:function pushbutton5_Callback(hObject, eventdata, handles)% hObject handle to pushbutton5 (see GCBO)% eventdata

9、reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)clc;close all;close(gcf);clear;4.2图像转化为灰度图像因为matlab中较多的图像处理函数支持对灰度图像进行处理,故对图像进行灰度转化十分必要。利用rgb2gray(X)函数对其他图像进行灰度图像的转化。-7-实现代码:function pushbutton11_Callback(hObject, eventdata, handles)%

10、 hObject handle to pushbutton11 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global Taxes(handles.axes2);T=getimage;x=rgb2gray(handles.img); %RGB?imshow(x);handles.img=x;guidata(hObject,handles);4.3图像直方图统计和直方图

11、均衡(1)通过histeq(X)函数实现直方图均衡。此函数只能对灰度图像进行直方图均衡所以要先将彩图转为灰度图像。实现代码:function pushbutton12_Callback(hObject, eventdata, handles)% hObject handle to pushbutton12 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global T-8

12、-axes(handles.axes2);T=getimage;h=histeq(handles.img); imshow(h);handles.img=h;guidata(hObject,handles); (2)直方图统计。通过利用 imhist(X)函数来实现直方图统计。实现代码:function pushbutton13_Callback(hObject, eventdata, handles)% hObject handle to pushbutton13 (see GCBO)% eventdata reserved - to be defined in a future versi

13、on of MATLAB% handles structure with handles and user data (see GUIDATA)axes(handles.axes2);x=imhist(handles.img); %直方图统计x1=x(1:10:256);horz=1:10:256;bar(horz,x1);axis(0 255 0 15000);set(handles.axes2,xtick,0:50:255);set(handles.axes2,ytick,0:2000:15000);4.4加入各种噪声,并通过几种滤波算法实现去噪(1)加入噪声。通过imnoise(I,ty

14、pe,parameters)来加入各种噪声。-9-加入椒盐噪声实现代码:function pushbutton6_Callback(hObject, eventdata, handles)% hObject handle to pushbutton6 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB-10-% handles structure with handles and user data (see GUIDATA)global Taxes(handles.axes2);T=geti

15、mage; prompt=数日椒盐噪声参数 1:;defans=0.02;p=inputdlg(prompt,input,1,defans);p1=str2num(p1);f=imnoise(handles.img,salt imshow(f);handles.img=f;guidata(hObject,handles);加入高斯噪声:-11-实现代码:function pushbutton10_Callback(hObject, eventdata, handles)% hObject handle to pushbutton10 (see GCBO)% eventdata reserved

16、 - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global Taxes(handles.axes2);T=getimage;prompt=输入高斯噪声 1:,输入高斯噪声 2;defans=0,0.02;p=inputdlg(prompt,input,1,defans);p1=str2num(p1);p2=str2num(p2);f=imnoise(handles.img,gaussian,p1,p2);imshow(f);han

17、dles.img=f;guidata(hObject,handles);加入乘性噪声:-12-实现代码:function pushbutton8_Callback(hObject, eventdata, handles)% hObject handle to pushbutton8 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global Taxes(handles.a

18、xes2);-13-T=getimage;prompt=输入乘性噪声 1:;defans=0.02;p=inputdlg(prompt,input,1,defans);p1=str2num(p1);f=imnoise(handles.img,speckle,p1);imshow(f);handles.img=f;guidata(hObject,handles); (2)滤除噪声(椒盐噪声)滤波前中值滤波后-14-实现代码:function pushbutton14_Callback(hObject, eventdata, handles)% hObject handle to pushbutt

19、on14 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global Taxes(handles.axes2);T=getimage;k=medfilt2(handles.img);imshow(k);handles.img=k;guidata(hObject,handles);线性滤波后-15-实现代码:function pushbutton16_Callback(hO

20、bject, eventdata, handles)% hObject handle to pushbutton16 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global Taxes(handles.axes2);T=getimage;h=1 1 1;1 1 1;1 1 1;H=h/9;i=double(handles.img);k=convn(i,h);imsho

21、w(k,);handles.img=k;guidata(hObject,handles);自适应滤波后实现代码:function pushbutton18_Callback(hObject, eventdata, handles)% hObject handle to pushbutton18 (see GCBO)-16-% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global Taxes

22、(handles.axes2);T=getimage;k=wiener2(handles.img,5,5);imshow(k);handles.img=k;guidata(hObject,handles);低通滤波器滤波后实现代码:function pushbutton19_Callback(hObject, eventdata, handles)% hObject handle to pushbutton19 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles struct

23、ure with handles and user data (see GUIDATA)axes(handles.axes2);y1=handles.img; f=double(y1); % 数据类型转换,matlab 不支持图像的无符号整型的计算g=fft2(f); % 傅里叶变换g=fftshift(g); % 转换数据矩阵M,N=size(g);nn=2; %二阶巴特沃斯低通滤波器d0=50; %截止频率 50m=fix(M/2); n=fix(N/2);for i=1:Mfor j=1:N-17-d=sqrt(i-m)2+(j-n)2);h=1/(1+0.414*(d/d0)(2*nn

24、); % 计算低通滤波器传递函数result(i,j)=h*g(i,j);endendresult=ifftshift(result);y2=ifft2(result);y3=uint8(real(y2);imshow(y3); % 显示处理后的图像高通滤波器滤波后实现代码:function pushbutton20_Callback(hObject, eventdata, handles)% hObject handle to pushbutton20 (see GCBO)% eventdata reserved - to be defined in a future version of

25、MATLAB% handles structure with handles and user data (see GUIDATA)axes(handles.axes2);x=(handles.img);f=double(x); % 数据类型转换k=fft2(f); % 傅里叶变换g=fftshift(k); % 转换数据矩阵M,N=size(g);nn=2;d0=25; %截止频率 25-18-m=fix(M/2); n=fix(N/2);for i=1:Mfor j=1:Nd=sqrt(i-m)2+(j-n)2); % 计算高通滤波器传递函数if d=d0h=0;else h=1;endr

26、esult(i,j)=h*g(i,j);endendresult=ifftshift(result);y2=ifft2(result);y3=uint8(real(y2);imshow(y3); % 显示滤波处理后的图像4.5 还原通过一个全局变量保存原始图像路径,在需要还原至原始图像时,重新读取该全局变量即可。实现代码:function pushbutton21_Callback(hObject, eventdata, handles)% hObject handle to pushbutton21 (see GCBO)% eventdata reserved - to be defined

27、 in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global S %还原axes(handles.axes2);y=imread(S);f=imshow(y);handles.img=y;guidata(hObject,handles);5、结果分析软件测试基本成功,课题所要求的功能均能较好实现。但一些功能只支持灰度图像的处理。其中值得一提的是在滤波处理中的低通滤波与高通滤波的效果。由于一般图像中含有较多的低频信息成分高频成分较少,通过低通滤波后,噪声以及高频成分

28、被滤除,图像虽有少量失真,略显模糊,但尚可辨识。但若是通过高通滤波后,大量的有效低频信息被滤除,图像严重失真,不可辨识。-19-当我第一次拿到此次的课题时,感到有些无所适从。虽然,曾经学习过 matlab 的课程,在课程的考核中也取得了较好的成绩,但由于对 matlab 的学习更多的只是停留在理论上的学习,在课时内的试验也只是简单的基础性试验, 所以对 matlab 实际运用不是很熟练。为此,在实践正式开始前,我利用课余时间,重新复习了 matlab 教材,专门借阅了利用matlab 进行图像处理的相关教程,通过索引网络上的相关资料,为课设做了较为充分的准备。在参考了相关材料及源程序,我对自己

29、要做的课设内容有了进一步的了解,并对 matlab 的使用有了更深的体会。当然,在课设的进行过程中,我还是遇到了不少问题。例如,起初由于我对句柄使用以及一些函数使用的不恰当,使得在对图像文件的保存上就遇到了问题,不过最后还是在老师的提示下解决了。随着课设的进行,对 matlab 的的熟悉度逐步加深。总体来说,此次的课程设计,还是较为满意的。它不但鞭策着我去巩固 matlab 的基础理论知识,还提高了我对 matlab 的实际操作运用,使得理论与实践相结合,为进一步学习 matlab 打下坚实的基础;同时,在实践的工程中,也让我体会到一种努力付出并得到回报的满足感觉。参考书目:(五号,宋体加粗

30、)1 数字图像处理(MATLAB)(第二版) 美Rafael C.Gonzalez 电子工业出版社附录:(五号,宋体加粗)function varargout = faded(varargin)% FADED MATLAB code for faded.fig% FADED, by itself, creates a new FADED or raises the existing% singleton*.% H = FADED returns the handle to a new FADED or the handle to% the existing singleton*.% FADED

31、(CALLBACK,hObject,eventData,handles,.) calls the local% function named CALLBACK in FADED.M with the given input arguments.% FADED(Property,Value,.) creates a new FADED or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before faded_OpeningFcn get

32、s called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to faded_OpeningFcn via varargin.% *See GUI Options on GUIDEs Tools menu. Choose “GUI allows only one% instance to run (singleton)“.-20-% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the a

33、bove text to modify the response to help faded% Last Modified by GUIDE v2.5 29-Dec-2015 22:05:27% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct(gui_Name, mfilename, .gui_Singleton, gui_Singleton, .gui_OpeningFcn, faded_OpeningFcn, .gui_OutputFcn, faded_OutputFcn, .gui_L

34、ayoutFcn, , .gui_Callback, );if nargin endif nargoutvarargout1:nargout = gui_mainfcn(gui_State, varargin:);elsegui_mainfcn(gui_State, varargin:);end% End initialization code - DO NOT EDIT% - Executes just before faded is made visible.function faded_OpeningFcn(hObject, eventdata, handles, varargin)%

35、This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to faded (see VARARGIN)% Choose default command line output for fa

36、dedhandles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes faded wait for user response (see UIRESUME)% uiwait(handles.figure1);% - Outputs from this function are returned to the command line.function varargout = faded_OutputFcn(hObject, eventdata, handles) % vara

37、rgout cell array for returning output args (see VARARGOUT);% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)-21-% Get default command line output from handles structurevarargout1 = handles.output;%

38、 - Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject handle to pushbutton2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)filename,pathname=uigetfile(*.j

39、pg;*.bmp;*.tif;*.*,载入图像 );if isequal(filename,0)|isequal(pathname,0)errordlg(没有选中文件 ,出错);return;else file=pathname,filename;global S %设置一个全局变量 S,保存初始图像路径,以便之后的还原操作S=file;x=imread(file);set(handles.axes1,HandleVisibility,ON);axes(handles.axes1);imshow(x);set(handles.axes1,HandleVisibility,OFF);axes(h

40、andles.axes2);imshow(x);handles.img=x;guidata(hObject,handles);end% - Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with

41、 handles and user data (see GUIDATA)sfilename ,sfilepath=uiputfile(*.jpg;*.bmp;*.tif;*.*,保存图像文件,untitled.jpg);if isequal(sfilename,sfilepath,0,0)sfilefullname=sfilepath ,sfilename;imwrite(handles.img,sfilefullname);elsemsgbox(你按了取消键, 保存失败 );end% - Executes on button press in pushbutton5.function pus

42、hbutton5_Callback(hObject, eventdata, handles)% hObject handle to pushbutton5 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)clc;close all;close(gcf);-22-clear;% - Executes on button press in pushbutton6.functio

43、n pushbutton6_Callback(hObject, eventdata, handles)% hObject handle to pushbutton6 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global Taxes(handles.axes2);T=getimage; prompt=数日椒盐噪声参数 1:;defans=0.02;p=inputdlg

44、(prompt,input,1,defans);p1=str2num(p1);f=imnoise(handles.img,salt imshow(f);handles.img=f;guidata(hObject,handles);% - Executes on button press in pushbutton7.function pushbutton7_Callback(hObject, eventdata, handles)% hObject handle to pushbutton7 (see GCBO)% eventdata reserved - to be defined in a

45、 future version of MATLAB% handles structure with handles and user data (see GUIDATA)% - Executes on button press in pushbutton8.function pushbutton8_Callback(hObject, eventdata, handles)% hObject handle to pushbutton8 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% han

46、dles structure with handles and user data (see GUIDATA)global Taxes(handles.axes2);T=getimage;prompt=输入乘性噪声 1:;defans=0.02;p=inputdlg(prompt,input,1,defans);p1=str2num(p1);f=imnoise(handles.img,speckle,p1);imshow(f);handles.img=f;guidata(hObject,handles); % - Executes on button press in pushbutton9.

47、function pushbutton9_Callback(hObject, eventdata, handles)% hObject handle to pushbutton9 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% - Executes on button press in pushbutton10.function pushbutton10_Callbac

48、k(hObject, eventdata, handles)-23-% hObject handle to pushbutton10 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global Taxes(handles.axes2);T=getimage;prompt=输入高斯噪声 1:,输入高斯噪声 2;defans=0,0.02;p=inputdlg(prompt,input,1,defans);p1=str2nu

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > 企业管理 > 管理学资料

本站链接:文库   一言   我酷   合作


客服QQ:2549714901微博号:道客多多官方知乎号:道客多多

经营许可证编号: 粤ICP备2021046453号世界地图

道客多多©版权所有2020-2025营业执照举报