首页 基于MATLAB的图像处理的课程设计

基于MATLAB的图像处理的课程设计

举报
开通vip

基于MATLAB的图像处理的课程设计目录 1、​ 课程设计目的…………………..………………….…………3 2、​ 课程设计要求………………….………………………….….3 3、​ 课程设计的内容………………..……………………….……3 4、​ 题目分析…………………………………………….………..3 5、​ 总体设计……………….…………………………….……….4 6、​ 具体设计……………………………………………..…….…5 6.1、​ 文件……………….………………………………….……….5 6.1.1、打开………..……………………………...

基于MATLAB的图像处理的课程设计
目录 1、​ 课程设计目的…………………..………………….…………3 2、​ 课程设计 要求 对教师党员的评价套管和固井爆破片与爆破装置仓库管理基本要求三甲医院都需要复审吗 ………………….………………………….….3 3、​ 课程设计的内容………………..……………………….……3 4、​ 题目分析…………………………………………….………..3 5、​ 总体设计……………….…………………………….……….4 6、​ 具体设计……………………………………………..…….…5 6.1、​ 文件……………….………………………………….……….5 6.1.1、打开………..…………………………………….…………5 6.1.2、保存…………………………..………………….…………5 6.1.3、退出……………………..……………………………….…5 6.2、编辑…………………………..…………………………………5 6.2.1、灰度………………………..………………….……………5 6.2.2、亮度………………………..………………….……………6 6.2.3、截图…………………………..……………….……………7 6.2.4、缩放…………………………..…………….………………7 6.3、旋转…………………………………..…………………………9 6.3.1、上下翻转…………………….…………….………………9 6.3.2、左右翻转………………….…………….…………………9 6.3.3任意角度翻转…………………………….…………………9 6.4、噪声………………………………..…………..………………10 6.5、滤波……………………………………..………..……………10 6.6、直方图统计……………………………………………………11 6.7、频谱分析………………………………………..…..…………12 6.7.1、频谱图………………………… …… …… ……….……12 6.7.2、通过高通滤波器………………… …… .……………….12 6.7.3、通过低通滤波器……………………………….. ……….13 6.8、灰度图像处理………………………………………… . …….14 6.8.1、二值图像……………………………………………….. .14 6.8.2、创建索引图像……………………………………… .…..14 6.9、颜色模型转换……………………………..… ..……………14 6.10、操作界面设计……………………………………..…………15 7、​ 程序调试及结果分析………………….……………………15 8、​ 心得体会…………………………….………………………16 9、​ 参考文献…………………………….………………………17 10、​ 附录………………………….………………………………18 基于MATLAB的图像处理的课程设计 摘要: 数字图像处理技术是20世纪60年代发展起来的一门新兴学科,随着图像处理理论和方法的进一步完善,使得数字图像处理技术在各个领域得到了广泛应用,并显示出广阔的应用前景。MATLAB既是一种直观、高效的计算机语言,同时又是一个科学计算平台。它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具。根据它提供的500多个数学和 工程 路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理 函数,工程技术人员和科学工作者可以在它的集成环境中交互或编程以完成各自的计算。 MATLAB中集成了功能强大的图像处理工具箱。由于MATLAB语言的语法特征与C语言极为相似,而且更加简单,更加符合科技人员对数学表达式的书写格式,而且这种语言可移植性好、可扩展性强,再加上其中有丰富的图像处理函数,所以MATLAB在图像处理的应用中具有很大的优势。 关键词:MATLAB,数字图像处理 1、​ 课程设计目的 综合运用MATLAB工具箱实现图像处理的GUI程序设计。 2、​ 课程设计要求 1)熟悉和掌握MATLAB 程序设计方法 2)掌握MATLAB GUI 程序设计 3)学习和熟悉MATLAB图像处理工具箱 4)学会运用MATLAB工具箱对图像进行处理和分析 3、​ 课程设计的内容 学习MATLAB GUI程序设计,利用MATLAB图像处理工具箱,设计和实现自己的Photoshop 。要求:按照软件工程方法,根据需求进行程序的功能分析和界面设计,给出设计详细说明。然后按照自己拟定的功能要求进行程序设计和调试。 以下几点是程序必须实现的功能。 1)图像的读取和保存。 2)设计图形用户界面,让用户能够对图像进行任意的亮度和对比度变化调整,显示和对比变换前后的图像。 3)设计图形用户界面,让用户能够用鼠标选取图像感兴趣区域,显示和保存该选择区域。 4)编写程序通过最近邻插值和双线性插值等算法将用户所选取的图像区域进行放大和缩小整数倍的操作,并保存,比较几种插值的效果。 5)图像直方图统计和直方图均衡,要求显示直方图统计,比较直方图均衡后的效果。 6)能对图像加入各种噪声,并通过几种滤波算法实现去噪并显示结果。比较去噪效果。 4、​ 题目分析 信息化社会中,计算机在各种信息处理中发挥着重要的作用。我们可以借助计算机,对数字图像进行处理,以达到不同的效果。根据题目的要求,除了实现要求的功能外,还有很多的功能需要用到。 (1)、将一个RGB图像转换为灰度图像。 (2)、可以对图像做各种变换,如旋转等。 (3)、有时并不需要图像显示其细节部分,只要其轮廓,这时候不要很高的灰度级。可以把图像转换为二值图像,进行图像腐蚀,或是创建索引图像等。 (4)、分析一个图像的频谱特征,利用傅里叶变换,将图像从空间域变换到频域,然后进行各种处理,经过高通滤波器或是低通滤波器。 (5)、为了科学地定量描述和使用颜色,人们提出了各种颜色模型,按用途可分为三类:计算颜色模型,视觉颜色模型和工业颜色模型。有时为了不同的需要,要对颜色模型进行转换。 5、​ 总体设计 由于要实现的功能并不是很多,所以在排版的过程中,把各个功能都安排在目录栏上,整体安排如下图所示: 同时在调节亮度时,虽然可以同对话框的形式输入调节的比例系数,但是这样效果不好了,不容易调节,因此这里考虑用滚动条来调节。因此,总体的设计界面如下图所示: 6、​ 具体设计 6.1、​ 文件 6.1.1、打开 为了让使用者更方便的使用,所以在设计的时候,通过对话框的形式来选择文件,选择uigetfile函数来实现,uigetfile函数显示一个打开文件对话框,该对话框自动列出当前路径下的目录和文件,由于这个GUI程序的操作对象是图像文件,所以设置这里的缺省后缀名为“.bmp”。 Uigetfile函数的调用格式为[name,path]=yigetfile(…), 在按下对话框中的执行按钮“打开”后,返回选择的文件名和路径,分别保存到“name”和“path”中。如果按下取消按钮或是发生错误,则返回值是0。 根据返回值的情况,如果是0,则弹出提示错误的对话框,否则,通过imread函数读出图像数据,把图像数据赋值给全局变量handles.img。 6.1.2、保存 同样也通过对话框的形式来保存图像数据,通过uigetfile函数选择文件名和路径,用getimage(gca)取出坐标2变换后的图像数据保存到变量i,最后用imwrite函数,把数据i存到指定的文件。 6.1.3、退出 退出比较简单,程序如下所示: clc; close all; close(gcf); 6.2、​ 编辑 6.2.1、灰度 由于RGB图像是三维图像,所以图像数据是一个三维数组,为了显示灰度图像,把三维图像降为二维,可以只取其中的二维数据,实现方法程序为: y=(handles.img(:,:,1)); %当然也可以选择(:,:,2) 或(:,:,3) imshow(y); 但是这样的话,根据程序所选的不同,图像数据也不同,显示也就不一样。 另一种方法就是,运用rgb2gray函数实现彩色图像到灰度图像的转换。程序为: y=rgb2gray(handles.img);  imshow(y); 这个程序只能用于RGB图像转换灰度图像,当原始图像本来就是灰度图像时,运行该程序时就会出错,但是使用者在使用时有时并不知道这些,为了使该程序更加完善,应该在使用者原先图像时灰度图像时使用该功能时,应该要显示提示类信息。所以在开始时应该要有一个RGB图像或是灰度图像的判断过程。完整的程序如下: if isrgb(handles.img) y=rgb2gray(handles.img);  imshow(y); else msgbox('这已经是灰度图像','转换失败'); end 如果原图是RGB,执行该操作的结果如下图: 如果原图本身已经是灰度图像了,执行该操作弹出如 右图所示的提示对话框 6.2.2、亮度 用imadjust函数,其调用格式如下: g=imadust(f,[low_in high_in],[low_out high_out]),gamma) gamma 表示映射性质,默认值是1 表示线性映射。 由于该函数有五个参数需要输入,为了方便用户改变,所以这里设计一个输入对话框,用户通过对话框把五个参数赋值给[low_in high_in],[low_out high_out],gamma这五个参数,如下一组命令建立了如图所示的输入对话框: prompt={'输入参数1','输入参数2','输入gamma'}; defans={'[0 0.7]','[0 1]','1'}; p=inputdlg(prompt,'输入参数',1,defans); 但是,这种方法并不能很好的让用户能够对图像进行任意的亮度和对比度变化调整,有时并不事先知道参数的值要多少,也不关心,而是任意调节的,直到满意为止。所以应该用滑动条来调节图像的亮度和对比度,这样更适合用户的使用习惯。 由于imadjust函数有五个参数,所以原则上需要设计五个滑动条来调节对比度,这对用户来说显然比较麻烦,因此在设计的时候固定其中的三个参数,通过调节两个参数的值来改变亮度和对比度。 [0 handles.beta],[0 1],handles.gm,这里的变量handles.beta和handles.gm就通过滑动条得到,滑动条设计如下图: 亮度调整的tag名为ld,取值范围0~1,gamma值的tag名为gamma,取值范围为0~5。获取滑动条参数的程序如下: handles.beta=get(handles.ld,'value'); handles.gm=get(handles.gamma,'value'); 执行该操作,调节滑动条到上图所示位置,结果如下图: 6.2.3、截图 在MATLAB中,用函数imcrop实现对图像的剪切操作。该操作剪切的是图像中的一个矩形子图,用户可以通过参数指定这个矩形四个顶点的坐标,也可以交互地用鼠标选取这个矩形。 Imcrop函数的调用格式如下: y=imcrop(handles.img); 不管handles.img是三维的还是二维数据,该函数都能进行操作。下图就是对三维图像的截图: 6.2.4、缩放 在MATLAB中,用函数imresize来实现对图像的放大或缩小。插值方法可选用三种方法,最近邻插值,双线性插值,双三次插值。 该函数的调用格式如下: B=imresize(A,m,method) 其中: 参数method用于指定插值的方法,可选的值为“nearest”(最近邻法),“bilinear”(双线性插值)、“bicubic”(双三次插值),缺省值为“nearest”。 B=imresizee(A.m,method)表示返回原图A的m倍放大图像(m小于1时实际上是缩小); 下图就是采用邻近插值法的放大和缩小图像,参数值保持默认设置: 虽然处理后看不出放大的效果,这是由于坐标轴限制的原因,如果把处理后的图片保存起来,再把处理后的文件打开,就可以看到比较明显的放大效果。 缩小后的结果如下: 6.3、​ 旋转 6.3.1、上下翻转 函数flipud是实现一个二维矩阵的上下翻转,如 a=[1 2;3 4],经过该函数处理后,原矩阵变为[3 4;1 2];所以利用该函数也可以对图像进行上下翻转处理,但由于该函数针对二维数据的处理,所以在写程序时,要对RGB图像和灰度图像分开处理,这就要用到isrgb函数来判断,如果是灰度图像,则可以直接用这个函数进行处理,否则就要对RGB图像进行降维处理, for k=1:3 y(:,:,k)=flipud(x(:,:,k)); end 处理结果如图: 6.3.2、左右翻转 对图像的左右翻转也可以用fliplr函数来处理,同样的,也要对灰度和彩色图像分开处理,处理结果如图: 6.3.3任意角度翻转 用函数imrotate来实现对图像的插值旋转。 该函数的调用格式如下: B=imrotate(A,angle,method,’crop’) 其中,参数method用于指定插值的方法,可选的值可以有三种,分别为邻近插值,双线性插值,双三次插值,缺省时为邻近插值,参数angle代表旋转的角度。 一般来说,旋转后的图像会比原图大,用户可以指定“crop”参数对旋转后的图像进行剪切(取图像的中间部分),使返回的图像与原图大小相同。执行结果为: 6.4、​ 噪声 经常用到的噪声有三种,高斯噪声,椒盐噪声,乘性噪声,可以通过以下三个函数来实现: y=imnoise(handles.img,'gaussian',p1,p2);%高斯噪声 y=imnoise(x,'salt & pepper',p1); %椒盐噪声 y=imnoise(handles.img,'speckle',p1); %乘性噪声 p1,p2的参数也通过输入对话框的形式得到,原图加入高斯噪声后结果如下所示: 6.5、​ 滤波 这里选择三种滤波方法, k=medfilt2(handles.noise_img);%中值滤波 k=wiener2(handles.noise_img,[5,5]);%自适应滤波 k=filter2(fspecial('average',3),handles.noise_img)/255;%平滑滤波 同样的,这些函数也是针对二维数据,所以要先判断是彩色图像还是灰度图像,然后分别进行处理,下面是对彩色图像的自适应滤波处理: i=handles.noise_img; if isrgb(i) a=handles.noise_img(:,:,1); b=handles.noise_img(:,:,2); c=handles.noise_img(:,:,3); k(:,:,1)=wiener2(a,[5,5]); k(:,:,2)=wiener2(b,[5,5]); k(:,:,3)=wiener2(c,[5,5]); imshow(k); 执行结果如图: 对其他方法的滤波程序也类似,由于把各个滤波方法放在一个选择框里,所以程序要用以下的选择语句: switch str case '中值滤波' case '自适应滤波' case '平滑滤波' end 6.6、​ 直方图统计 用imhist函数对图像数据进行直方图统计, x=imhist(handles.img(:,:,1)); bar(horz,x); 其中,x矩阵的数据是0~255灰度值的统计个数,如果直接对x矩阵数据进行图形图显示,由于有256个数据,在坐标系中就会很密集,为了更清楚的显示条形图,所以在程序设计时,把x数据进行部分提取, x1=x(1:10:256); horz=1:10:256; bar(horz,x1); 除了显示数据的直方图统计外,还可以对图像进行均衡处理,所用到的函数是histeq,这两个函数同样只使用于二维数据,所以也要对二维和三维数据分开处理。 直方图显示和均衡后的图像分别如下图所示 6.7、​ 频谱分析 6.7.1、频谱图 为了得到图像的频谱图,先要对数据进行傅里叶变换,用fft2函数对二维数据进行快速傅里叶变换,同时为了更好的观察频谱图,需要把fft2变换后的数据进行平移,利用fftshift函数,把快速傅里叶变换的DC 组件移到光谱中心。这样图像能量的低频成分将集中到频谱中心,图像上的边缘、线条细节信息等高频成分将分散在图像频谱的边缘。如下图所示: 6.7.2、通过高通滤波器 axes(handles.axes2); x=(handles.img); if isrgb(x) msgbox('这是彩色图像,不能通过高通滤波器','失败'); else y1=imnoise(x,'gaussian'); %加高斯噪声 f=double(y1); % 数据类型转换 k=fft2(f); % 傅立叶变换 g=fftshift(k); % 转换数据矩阵 [M,N]=size(g); nn=2; d0=3; %截止频率为3 m=fix(M/2); n=fix(N/2); for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2); % 计算高通滤波器传递函数 if d<=d0 h=0; else h=1; end result(i,j)=h*g(i,j); end end result=ifftshift(result); y2=ifft2(result); y3=uint8(real(y2)); imshow(y3); end 6.7.3、通过低通滤波器 axes(handles.axes2); x=(handles.img); if isrgb(x) msgbox('这是彩色图像,不能通过低通滤波器','失败'); else y1=imnoise(x,'salt & pepper'); % 叠加椒盐噪声 f=double(y1); % 数据类型转换,不支持图像的无符号整型的计算 g=fft2(f); % 傅立叶变换 g=fftshift(g); % 转换数据矩阵 [M,N]=size(g); nn=2; % 二阶巴特沃斯(Butterworth)低通滤波器 d0=10; %截止频率为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+0.414*(d/d0)^(2*nn));% 计算低通滤波器传递函数 result(i,j)=h*g(i,j); end end result=ifftshift(result); y2=ifft2(result); y3=uint8(real(y2)); imshow(y3); % 显示滤波处理后的图像 end 6.8、​ 灰度图像处理 6.8.1、二值图像 用j=im2bw(x);来对灰度图像到二值图像的转换。转换结果为: 6.8.2、创建索引图像 用X = grayslice(I,n)函数来实现,转换后得到的图像为: 6.9、​ 颜色模型转换 下面是RGB颜色模型到HSV模型的转换程序: axes(handles.axes2); x=(handles.img); if isrgb(x) HSV=rgb2hsv(x); imshow(HSV); else msgbox('这是灰度图像,不能转换','转换失败'); end 转换后的结果为: 同样的,转化为ntsc和ycbcr模型,只要改为相应的函数即可。函数如下: rgb2ntsc(x);% rgb模型转ntsc模型 rgb2ycbcr(x);% rgb模型转ycbcr模型 6.10、​ 操作界面设计 为了使整个操作界面更加的美观,对背景颜色和字体颜色进行设置。 7、​ 程序调试及结果分析 在程序设计过程中,碰到很多的问题。 1、​ 一个函数只能对灰度图像处理,不能对RGB图像处理,那么如何才能对RGB图像处理呢? 这主要是对MATLAB函数的不够清楚,用到的很多函数是针对二维数据的,而RGB图像的数据是一个三维矩阵,所以处理要与灰度图像不同,在开始的时候,我认为应该找一个能够应用于三维矩阵的函数,结果却没找到,后来想到可以把三维数据进行降维处理,同样使用二维的函数,只要是同样处理三次。比如,彩色图像的滤波处理,直方图均衡等。 2、​ 同一个操作对灰度图像可以使用,当用户选择的是彩色图像时,该操作就会出错? 这是由于设计程序时,本身程序是有针对性的,有些程序只能对灰度图像有效,有些对彩色图像有效,但是用户并不清楚这些,所以在设计的过程中就要考虑全面,要分开设计。在本次课程设计中,我均对每个程序的开始时,用if isrgb(x)进行判断。 3、​ 在没有加入噪声的情况下,点击“中值滤波”或是其它滤波,会提示错误? 这是由于在滤波程序设计的过程中,开始用到的变量是handles.noise_img,而这个变量是在加入噪声时候才定义的,所以在没有加入噪声的情况下,点击各个滤波就会弹出变量没有定义的错误,解决方法就是在文件打开的时候就给定义handles.noise_img=x。 4、​ 本次设计存在一个比较大的问题,就是每次操作都是独立的,比如:要对图像加入噪声,然后在此基础上进行亮度调整,截图,频谱分析等等就不行。 我认为应该是整个程序设计过程中对变量的设置没有做好,应该每次操作后,把处理后的数据保存在一个全局变量,这样还要对处理后的数据进行在处理时,只要把这个全局变量作为原始数据带入就可。 8、​ 心得体会 1、在这次课程设计过程中,感触很深,由于对MATLAB图像处理的函数不熟悉,导致自己走了很多的弯路,比如在设计图像左转90度和右转90度时,由于开始并不知道imrotate函数,只知道上下翻转flipud和左右翻转fliplr函数,想着要怎么用这两个函数来实现左转和右转呢,如果当当只用矩阵转置的话,并不能达到要求,后来想可以结合flipud和fliplr函数来实现,下面程序是实现左转90度: axes(handles.axes2); x=(handles.img); if isrgb(handles.img) a=x(:,:,1); b=x(:,:,2); c=x(:,:,3); e=a'; g=b'; f=c'; y(:,:,1)=e; y(:,:,2)=g; y(:,:,3)=f; axes(handles.axes2); for k=1:3 f(:,:,k)=flipud(y(:,:,k)); end imshow(f); else m=x'; y=flipud(m); imshow(y); end 同理,也可以通过转置和fliplr函数实现右转90度。 但是后来发现其实只要用imrotate函数就可以解决问题。通过这次经历后,后来在设计其它程序时,尽量找MATLAB自带的现成函数,而不是一碰到问题就自己想算法,这样可以节省较多的时间。 2、在理工科的专业应用背景下,用matlab 进行相关计算与仿真编程的优势非常突出。特定的问题处理算法,我们通常都以M文件的文本形式给定最终的解决 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 ,自己设计的程序是在MATLAB环境下,用MATLAB语言编写的,这对于有安装MATLAB软件的计算机上运行并不存在什么问题,关键是一般的计算机很少有去安装MATLAB软件的,那么要在这些机子上运行该程序要怎么办呢,难道要先安装MATLAB软件,这显然太麻烦了,对于一个通用的,比较成熟的解决方案,我们当然期望它能应用到更多的场合,而.exe (可执行)文件可运行于所有的通用WINDOWS操作系统,为此,将M文件转换成.exe文件倒是个不错的想法。 通过查找资料,了解到在生成可执行性文件之前,需要进行一系列编译环境方面的配置: 要将用Matlab语言编写的函数文件编译成可独立执行的*.exe文件(即可脱离Matalab环境的执行程序),首先要安装和配置好Matlab Compiler,一般来说,在安装Matlab时就已经安装了相应版本的Matlab Compiler。只是不同版本的Matlab,其编译器的使用方法有一定的差异,这一点要引起一定的注意。 在确定安装好Matlab Compiler后,还需要对Compiler进行适当的配置,方法是在Matlab命令窗口输入: Mbuild –setup 然后根据提示执行相应的操作,使用者可根据自己计算机中现有编译器的情况选择合适的编译器,如VC++ 6.0、VC++7.0、Bland C的编译器等,目前Matlab好象还不支持VC++8.0(我计算机安装的就是VC++2005,Matlab就无法识别)。当然,如果你的计算机里根本就没有安装其他任何语言的编译器,也可选择Matlab自带的Lcc编译器,其实这个编译器对大多数用户已经够用了(我就是选择的Matlab自带的Lcc编译器)。 配置好编译器后,自然就是对自己编写的M文件进行编译了。 将M文件编译为独立可执行文件的语法是: >>mcc –m fun1.m fun2.m….. 其中fun1就是最后的可执行文件的名称。 3、通过本次课程设计,使自己对MATLAB GUI设计流程有了比较深刻的体会,同时也了解了一般软件设计的过程。在设计过程中碰到了很多的问题,通过这些问题,使自己分析问题,解决问题的能力得到了较大的提高。 9、​ 参考文献 【1】​ 郑阿奇,曹戈,赵阳.MATLAB实用教程[M].北京:电子工业出版社 【2】​ 程卫国,冯峰,姚东,徐听.MATLAB5.3应用指南[M].北京:人民邮电出版社 【3】​ 陈杨.MATLAB 6.X图像编程与图像处理[M].西安:西安电子科技大学出版社 10、​ 附录 function varargout = two(varargin) % Begin initialization code - DO NOT EDIT gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @two_OpeningFcn, ... 'gui_OutputFcn', @two_OutputFcn, ... 'gui_LayoutFcn', [] , ... 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); end % End initialization code - DO NOT EDIT % --- Executes just before erzhi is made visible. function two_OpeningFcn(hObject, eventdata, handles, varargin) % This function has no output args, see OutputFcn. --------------------------------------------------------- %接下来是菜单的程序 --------------------------------------------------------- function file_Callback(hObject, eventdata, handles) -------------------------------------------------------- %文件打开程序 function open_Callback(hObject, eventdata, handles) [name,path]=uigetfile({'*.bmp'},'载入图像'); if isequal(name,0)|isequal(path,0) errordlg('没有选中文件','出错'); return; else x=imread([path,name]); axes(handles.axes1); imshow(x); handles.img=x; handles.noise_img=x; guidata(hObject,handles) end --------------------------------------------------- function save_Callback(hObject, eventdata, handles) %文件保存 [filename,pathname] = uiputfile('*.bmp','图片保存为'); if isequal([filename,pathname],[0,0]) errordlg('没有保存','出错'); return; else file=strcat(pathname,filename); (handles.axes2); i=getimage(gca); imwrite(i,file); end ----------------------------------------------------- function exit_Callback(hObject, eventdata, handles) clc; close all; close(gcf); --------------------------------------------------------- %编辑菜单的程序 ------------------------------------------------------ function edit_Callback(hObject, eventdata, handles) ------------------------------------------------- function huidu_Callback(hObject, eventdata, handles) %灰度处理 axes(handles.axes2); if isrgb(handles.img) y=rgb2gray(handles.img); %RGB imshow(y); else msgbox('这已经是灰度图像','转换失败'); end ------------------------------------------------------ function liangdu_Callback(hObject, eventdata, handles) %亮度处理 prompt={'输入参数1','输入参数2','输入gamma'}; defans={'[0 0.7]','[0 1]','1'}; p=inputdlg(prompt,'输入参数',1,defans); p1=str2num(p{1}); p2=str2num(p{2}); p3=str2num(p{3}); gamma=p3; x=(handles.img); y=imadjust(x,p1,p2,gamma); axes(handles.axes2); imshow(y); ------------------------------------------------------ function jietu_Callback(hObject, eventdata, handles) %截图 set(handles.axes2,'HandleVisibility','ON'); axes(handles.axes2); y=imcrop(handles.img); imshow(y); handles.Timage=y; ------------------------------------------------------ function fangda_Callback(hObject, eventdata, handles) ------------------------------------------------------ function lingjz_Callback(hObject, eventdata, handles) %邻近插值放大 axes(handles.axes2); prompt={'输入放大倍数:'}; defans={'2'}; p=inputdlg(prompt,'输入放大倍数',1,defans); p1=str2num(p{1}); y=imresize(handles.img,p1,'nearest'); %最近邻插值法放大 imshow(y); ------------------------------------------------------ function shuangxianfa_Callback(hObject, eventdata, handles) %双线放大 axes(handles.axes2); prompt={'输入放大倍数:'}; defans={'2'}; p=inputdlg(prompt,'输入放大倍数',1,defans); p1=str2num(p{1}); y=imresize(handles.img,p1,'bilinear'); %最近邻插值法放大 imshow(y); ------------------------------------------------ function shuangsanci_Callback(hObject, eventdata, handles) %双三次放大 axes(handles.axes2); prompt={'输入放大倍数:'}; defans={'2'}; p=inputdlg(prompt,'输入放大倍数',1,defans); p1=str2num(p{1}); y=imresize(handles.img,p1,'bicubic'); imshow(y); ------------------------------------------------------- function suoxiao_Callback(hObject, eventdata, handles) %邻近插值缩小 axes(handles.axes2); prompt={'输入放大倍数:'}; defans={'0.2'}; p=inputdlg(prompt,'输入放大倍数',1,defans); p1=str2num(p{1}); y=imresize(handles.img,p1,'nearest'); %最近邻插值法缩小 imshow(y); ------------------------------------------------------ function shuangxian_Callback(hObject, eventdata, handles) %双线性缩小 axes(handles.axes2); prompt={'输入放大倍数:'}; defans={'0.2'}; p=inputdlg(prompt,'输入放大倍数',1,defans); p1=str2num(p{1}); y=imresize(handles.img,p1,'bilinear'); %最近邻插值法缩小 imshow(y); -------------------------------------------------------- function shuangsancisuox_Callback(hObject, eventdata, handles) %双三次缩小 axes(handles.axes2); prompt={'输入放大倍数:'}; defans={'0.2'}; p=inputdlg(prompt,'输入放大倍数',1,defans); p1=str2num(p{1}); y=imresize(handles.img,p1,'bicubic'); imshow(y); ------------------------------------------------------ function xuanzhuan_Callback(hObject, eventdata, handles) ------------------------------------------------------- function updown_Callback(hObject, eventdata, handles) %上下翻转 axes(handles.axes2); x=(handles.img); if isrgb(handles.img) for k=1:3 y(:,:,k)=flipud(x(:,:,k));%上下翻转函数 end imshow(y); else x=(handles.img); y=flipud(x); imshow(y); end --------------------------------------------------------- function leftright_Callback(hObject, eventdata, handles) %左右翻转 axes(handles.axes2); if isrgb(handles.img) x=(handles.img); for k=1:3 y(:,:,k)=fliplr(x(:,:,k));%左右翻转函数 end imshow(y); else x=(handles.img); y=fliplr(x); imshow(y); end ------------------------------------------------------ function zeft90_Callback(hObject, eventdata, handles) %左转90度 axes(handles.axes2); x=(handles.img); y=imrotate(x,90); imshow(y); ------------------------------------------------------- function right90_Callback(hObject, eventdata, handles) %右转90度 axes(handles.axes2); x=(handles.img); y=imrotate(x,-90); imshow(y); ------------------------------------------------------ function other_Callback(hObject, eventdata, handles) %任意角度旋转 axes(handles.axes2); prompt={'输入参数1:'}; defans={'30'}; p=inputdlg(prompt,'输入参数',1,defans); p1=str2num(p{1}); y=imrotate(handles.img,p1); imshow(y); ----------------------------------------------------- function gs_Callback(hObject, eventdata, handles) %加入高斯噪声 axes(handles.axes2); prompt={'输入参数1:','输入参数2'}; defans={'0','0.02'}; p=inputdlg(prompt,'输入参数',1,defans); p1=str2num(p{1}); p2=str2num(p{2}); y=imnoise(handles.img,'gaussian',p1,p2); imshow(y); handles.noise_img=y; guidata(hObject,handles); ------------------------------------------------------- function jy_Callback(hObject, eventdata, handles) %加入椒盐噪声 prompt={'输入参数1:'}; %对话框的设置,用户输入的是字符串 defans={'0.02'}; %缺省值 p=inputdlg(prompt,'输入参数',1,defans); p1=str2num(p{1}); %字符串转化为数值 axes(handles.axes2); x=(handles.img); y=imnoise(x,'salt & pepper',p1); imshow(y); handles.noise_img=y; guidata(hObject,handles); -------------------------------------------------------- function cx_Callback(hObject, eventdata, handles) %加入乘性噪声 axes(handles.axes2); prompt={'输入参数1:'}; defans={'0.02'}; p=inputdlg(prompt,'输入参数',1,defans); p1=str2num(p{1}); y=imnoise(handles.img,'speckle',p1); imshow(y); handles.noise_img=y; guidata(hObject,handles); ----------------------------------------------------- function zhifangtutongji_Callback(hObject, eventdata, handles) --------------------------------------------------------- function red_Callback(hObject, eventdata, handles) %R直方图 set(handles.axes2,'HandleVisibility','ON'); axes(handles.axes2); x=imhist(handles.img(:,:,1)); %直方图统计 x1=x(1:10:256); horz=1:10:256; bar(horz,x1); set(handles.axes2,'xtick',0:50:255); --------------------------------------------------------- function gray_Callback(hObject, eventdata, handles) %G直方图 set(handles.axes2,'HandleVisibility','ON'); axes(handles.axes2); if isrgb(handles.img) x=imhist(handles.img(:,:,2)); %直方图统计 x1=x(1:10:256); horz=1:10:256; bar(horz,x1); set(handles.axes2,'xtick',0:50:255); else msgbox('这是灰度图像','旋转失败'); end ------------------------------------------------------ function blue_Callback(hObject, eventdata, handles) %B直方图 set(handles.axes2,'HandleVisibility','ON'); axes(handles.axes2); if isrgb(handles.img) x=imhist(handles.img(:,:,3)); %直方图统计 x1=x(1:10:256); horz=1:10:256; bar(horz,x1); %axis([0 255 0 150000]); set(handles.axes2,'xtick',0:50:255); %set(handles.axes2,'ytick',0:2000:15000); else msgbox('这是灰度图像','旋转失败'); end --------------------------------------------------------- function junheng_Callback(hObject, eventdata, handles) %直方图均衡 set(handles.axes2,'HandleVisibility','ON'); axes(handles.axes2); if isrgb(handles.img) a=histeq(handles.img(:,:,1)); b=histeq(handles.img(:,:,2)); c=histeq(handles.img(:,:,3)); k(:,:,1)=a; k(:,:,2)=b; k(:,:,3)=c; imshow(k); else h=histeq(handles.img); %直方图均衡 imshow(h); end -------------------------------------------------------- %频谱分析 ----------------------------------------------------- function pinpu_Callback(hObject, eventdata, handles) -------------------------------------------------------- function pinputu_Callback(hObject, eventdata, handles) %显示频谱图 axes(handles.axes2); x=(handles.img); if isrgb(x) m=fft2(x(:,:,1)); y=fftshift(m); imshow(log(abs(y)),[]); else m=fft2(x); y=fftshift(m); imshow(log(abs(y)),[]); end -------------------------------------------------------- function frequency_Callback(hObject, eventdata, handles) %低通滤波器 axes(handles.axes2); x=(handles.img); if isrgb(x) msgbox('这是彩色图像,不能通过低通滤波器','失败'); else y1=imnoise(x,'salt & pepper'); % 叠加椒盐噪声 f=double(y1); % 数据类型转换,MATLAB不支持图像的无符号整型的计算 g=fft2(f); % 傅立叶变换 g=fftshift(g); % 转换数据矩阵 [M,N]=size(g); nn=2; % 二阶巴特沃斯(Butterworth)低通滤波器 d0=10; %截止频率为10 m=fix(M/2); n=fix(N/2)
本文档为【基于MATLAB的图像处理的课程设计】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
下载需要: 免费 已有0 人下载
最新资料
资料动态
专题动态
is_347086
暂无简介~
格式:doc
大小:1MB
软件:Word
页数:25
分类:
上传时间:2010-10-23
浏览量:34