数字图像处理学（第四版）MATLAB函数源代码

PAGE1-MATLAB函数源代码functionvarargout=DIP(varargin)%DIPMATLABcodeforDIP.fig%DIP,byitself,createsanewDIPorraisestheexisting%singleton*.%%EdittheabovetexttomodifytheresponsetohelpDIP%LastModifiedbyGUIDEv2.527-May-201211:43:05gui_Singleton=1;gui_State=struct('gui_Name',mfilename,...'gui_Singleton',gui_Singleton,...'gui_OpeningFcn',@DIP_OpeningFcn,...'gui_OutputFcn',@DIP_OutputFcn,...'gui_LayoutFcn',[],...'gui_Callback',[]);ifnargin&&ischar(varargin{1})gui_State.gui_Callback=str2func(varargin{1});endifnargout[varargout{1:nargout}]=gui_mainfcn(gui_State,varargin{:});elsegui_mainfcn(gui_State,varargin{:});end%Endinitializationcode-DONOTEDIT%ExecutesjustbeforeDIPismadevisible.functionDIP_OpeningFcn(hObject,eventdata,handles,varargin)%Thisfunctionhasnooutputargs,seeOutputFcn.%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%varargincommandlineargumentstoDIP(seeVARARGIN)%ChoosedefaultcommandlineoutputforDIPhandles.output=hObject;%Updatehandlesstructureguidata(hObject,handles);%UIWAITmakesDIPwaitforuserresponse(seeUIRESUME)%uiwait(handles.DIP);set(handles.Menu2,'Enable','off');%在打开图片前其他下拉菜单不可用set(handles.Menu3,'Enable','off');set(handles.Menu4,'Enable','off');set(handles.Menu5,'Enable','off');set(handles.Menu6,'Enable','off');set(handles.Menu7,'Enable','off');set(handles.Menu8,'Enable','off');set(handles.Menu9,'Enable','off');setappdata(handles.DIP,'IM',0);%Outputsfromthisfunctionarereturnedtothecommandline.functionvarargout=DIP_OutputFcn(hObject,eventdata,handles)%varargoutcellarrayforreturningoutputargs(seeVARARGOUT);%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%Getdefaultcommandlineoutputfromhandlesstructurevarargout{1}=handles.output;%Executesonbuttonpressinfile_open.functionfile_open_Callback(hObject,eventdata,handles)%hObjecthandletofile_open(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM%定义一个全局变量im[filename,pathname]=uigetfile({'*.bmp';'*.tif';'*.png';'*.jpg'},'selectpicture');%选择图片路径ifisequal(filename,0)||isequal(pathname,0),return;endstr=[pathnamefilename];%合成路径+文件名IM=imread(str);%读取图片axes(handles.axes1);%使用第一个axesimshow(IM);%显示图片%axes(handles.axes2);%clareset;set(handles.Menu2,'Enable','on');%在打开图片后其他下拉菜单可用set(handles.Menu3,'Enable','on');set(handles.Menu4,'Enable','on');set(handles.Menu5,'Enable','on');set(handles.Menu6,'Enable','on');set(handles.Menu7,'Enable','on');set(handles.Menu8,'Enable','on');set(handles.Menu9,'Enable','on');setappdata(handles.DIP,'IM',IM);%functionMenu1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionMenu1_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu1_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM%定义一个全局变量IM[filename,pathname]=uigetfile({'*.bmp';'*.tif';'*.png';'*.jpg'},'selectpicture');%选择图片路径ifisequal(filename,0)||isequal(pathname,0),return;endstr=[pathnamefilename];%合成路径+文件名IM=imread(str);%读取图片axes(handles.axes1);%使用第一个axesimshow(IM);%显示图片%axes(handles.axes2);%clareset;set(handles.Menu2,'Enable','on');%在打开图片后其他下拉菜单可用set(handles.Menu3,'Enable','on');set(handles.Menu4,'Enable','on');set(handles.Menu5,'Enable','on');set(handles.Menu6,'Enable','on');set(handles.Menu7,'Enable','on');set(handles.Menu8,'Enable','on');set(handles.Menu9,'Enable','on');setappdata(handles.DIP,'IM',IM);%functionMenu1_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu1_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalBW%定义处理后的图片BW这个全局变量[filename,pathname,filterindex]=...uiputfile({'*.bmp';'*.tif';'*.png';'*.jpg'},'savepicture');%存储图片路径iffilterindex==0return%如果取消操作，返回elsestr=[pathnamefilename];%合成路径+文件名axes(handles.axes2);%使用第二个axesimwrite(BW,str);%写入图片信息，即保存图片end%functionMenu1_3_Callback(hObject,eventdata,handles)%hObjecthandletoMenu1_3(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)close(gcf)%关闭当前Figure窗口句柄%functionMenu1_4_Callback(hObject,eventdata,handles)%hObjecthandletoMenu1_4(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;axes(handles.axes2);imshow(IM);[filename,pathname,filterindex]=...uiputfile({'*.bmp';'*.tif';'*.png';'*.jpg'},'savepicture');iffilterindex==0returnelsestr=[pathnamefilename];axes(handles.axes2);imwrite(IM,str);%写入图片信息，即保存图片end%functionMenu2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionMenu2_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionMenu2_1_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_1_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=imresize(IM,0.3,'nearest');axes(handles.axes2);imshow(BW,[]);%functionMenu2_1_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_1_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=imresize(IM,0.1,'bilinear');axes(handles.axes2);imshow(BW,[]);%functionMenu2_1_3_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_1_3(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=imresize(IM,0.1,'bicubic');axes(handles.axes2);imshow(BW,[]);%functionMenu2_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;H=rotate;%functionMenu2_3_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_3(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;set(handles.main_txt,'FontName','宋体');set(handles.main_txt,'String','请选择要裁剪的局域，并双击选定区域以显示！');waitforbuttonpress;BW=imcrop(IM);axes(handles.axes2);imshow(BW);set(handles.main_txt,'FontName','Calibri');set(handles.main_txt,'String','Welcome!');%functionMenu2_4_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_4(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionMenu2_4_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_4_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;snoise=0.15*randn(size(IM));BW=imadd(IM,im2uint8(snoise));axes(handles.axes2);imshow(BW);%functionMenu2_4_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_4_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=imnoise(IM,'salt&pepper',0.05);axes(handles.axes2);imshow(BW);%functionMenu2_5_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_5(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;len=25;theta=35;PSF=fspecial('motion',len,theta);BW=imfilter(IM,PSF,'circular','conv');axes(handles.axes2);imshow(BW);%functionMenu2_6_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_6(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)set(handles.main_txt,'FontName','宋体');set(handles.main_txt,'String','请用鼠标选择任意几个像素点后按回车以显示所选像素点的数据值！');vals=impixel;set(handles.main_txt,'FontName','Calibri');set(handles.main_txt,'String','Welcome!');%functionMenu2_7_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_7(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)set(handles.main_txt,'FontName','宋体');set(handles.main_txt,'String','请用鼠标选择一线段后按回车以显示轨迹强度图！');axes(handles.axes1);improfile;%c=improfile;%axes(handles.axes2);%imshow(c);set(handles.main_txt,'FontName','Calibri');set(handles.main_txt,'String','Welcome!');%functionMenu2_8_Callback(hObject,eventdata,handles)%hObjecthandletoMenu2_8(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;[counts,x]=imhist(IM,64);axes(handles.axes2);stem(x,counts,'k-','k.');%functionMenu4_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionMenu4_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=imadjust(IM);axes(handles.axes2);imshow(BW);%functionMenu4_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionMenu4_2_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_2_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=histeq(IM);axes(handles.axes2);imshow(BW);figure(1);subplot(1,2,1);imhist(IM);title('原始图像直方图');subplot(1,2,2);imhist(BW);title('直方图均衡化后图像直方图');%functionMenu4_2_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_2_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;hgram=ones(1,256);BW=histeq(IM,hgram);axes(handles.axes2);imshow(BW);figure(2);subplot(1,2,1);imhist(IM);title('原始图像直方图');subplot(1,2,2);imhist(BW);title('直方图规定化后图像直方图');%functionMenu4_3_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_3(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionMenu4_3_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_3_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;h=[111;101;111];%8邻域平均滤波h=h/8;BW=conv2(IM,h);axes(handles.axes2);imshow(BW);%functionMenu4_3_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_3_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=medfilt2(IM);%用于3*3滤波窗口对图像进行中值滤波处理%若用[m,n],K=medfilt2(IM,[m,n])axes(handles.axes2);imshow(BW);%functionMenu4_3_3_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_3_3(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;[f1,f2]=freqspace(size(IM),'meshgrid');生成频率序列矩阵Hd=ones(size(IM));%构造滤波器大小r=sqrt(f1.^2+f2.^2);%构造滤波器决策函数Hd(r>0.5)=0;%构造滤波器，半径r>0.5Y=fft2(double(IM));%对I进行傅里叶变换Y=fftshift(Y);%频谱平移Ya=Y.*Hd;%滤波Ya=ifftshift(Ya);%反变换BW=ifft2(Ya);axes(handles.axes2);imshow(uint8(BW));%functionMenu4_4_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_4(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionMenu4_4_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_4_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;h=[-1,-1,-1;-1,9,-1;-1,-1,-1];%八领域BW=imfilter(IM,h);axes(handles.axes2);imshow(uint8(BW));%functionMenu4_4_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu4_4_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;[f1,f2]=freqspace(size(IM),'meshgrid');%生成频率序列矩阵Hd=ones(size(IM));%构造滤波器大小r=sqrt(f1.^2+f2.^2);%构造滤波器决策函数Hd(r<0.1)=0;%构造滤波器，半径r<0.1Y=fft2(double(IM));%对I进行傅里叶变换Y=fftshift(Y);%频谱平移Ya=Y.*Hd;%滤波Ya=ifftshift(Ya);%反变换BW=ifft2(Ya);axes(handles.axes2);imshow(uint8(BW));%functionMenu5_3_1_Callback(hObject,eventdata,handles)%hObjecthandletoMenu5_3_1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=edge(IM,'roberts');axes(handles.axes2);imshow(BW);%functionMenu5_3_2_Callback(hObject,eventdata,handles)%hObjecthandletoMenu5_3_2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=edge(IM,'sobel');axes(handles.axes2);imshow(BW);%functionMenu5_3_3_Callback(hObject,eventdata,handles)%hObjecthandletoMenu5_3_3(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=edge(IM,'prewitt');axes(handles.axes2);imshow(BW);%functionMenu5_3_4_Callback(hObject,eventdata,handles)%hObjecthandletoMenu5_3_4(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=edge(IM,'canny');axes(handles.axes2);imshow(BW);%functionMenu5_3_5_Callback(hObject,eventdata,handles)%hObjecthandletoMenu5_3_5(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalIM;globalBW;BW=edge(IM,'zerocross');axes(handles.axes2);imshow(BW);%functionMenu5_4_Callback(hObject,eventdata,handles)%hObjecthandletoMenu5_4(seeGCBO)%eventdatareserved-tobe