MATLAB 边缘检测，图像滤波 工具

使用App Designer 写了一款可以实现如下功能的小程序：

• 卷积运算
• 中值滤波
• 高斯滤波
• 边缘检测
• 噪声添加
• 灰度化、二值化、反色

使用方式

1，导入图片 2，图像预处理 点击左上角将图片灰度化，黑白化或者，或添加噪声 注：该步骤为预处理只对原图生效 注：二值化及噪声添加的各种数据都可设置：

3，设置模板大小、横纵梯度算子、方差、卷积次数 注：

• 模板卷积及中值滤波只需设置横向梯度算子（即模板）模板卷积会自动除以权重和
• 中值滤波在模板区域输入1
• 边缘检测需同时设置横向及纵向梯度算子
• 高斯滤波只需设置space-Sigma及range-Sigma
• 输入数值可以为公式形式，例如sqrt(5)
• 空缺位置默认为0

4，开始卷积 举边缘检测例子 5，存储图像 存储图像按钮存储的为当前图像 批量存储按钮存储的是每次卷积的图像及预处理图像

输入批量存储图片名称及格式： 批量存储结果： 其中_0的图片为预处理图片

完整代码

function Template_Convolution_App %该程序用于实现图像模板卷积实验 %功能包括： %图像平滑（中值滤波，均值滤波，高斯滤波，双边滤波等） %图像边缘检测 %图像噪声添加 %全局变量>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> global TCA_Fig TCA_Axes
global TCA_BG_SE TCA_BG_35 TCA_Spinner PGdisplay
global Hx Hy SigmaS SigmaR T33 times
global Hx_ Hy_ Gx_ Gy_ Gxy_ SigmaS_ SigmaR_ Gs_ Gr_ Gsr_ Wx_ Wy_ coe_
global oriPic subOriPic noisePic convolutionPic
global PicSet
times=0; %界面制作>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> TCA_Fig=uifigure('units','pixels'); TCA_Fig.Position=[10,50,930,600]; TCA_Fig.NumberTitle='off'; TCA_Fig.MenuBar='none'; TCA_Fig.Name='模板卷积App'; TCA_Fig.Color=[1,1,1].*0.98; TCA_Axes=uiaxes('Parent',TCA_Fig); TCA_Axes.Position=[15,15,570,570]; TCA_Axes.XTick=[]; TCA_Axes.YTick=[]; TCA_Axes.XColor='none'; TCA_Axes.YColor='none'; TCA_Axes.Toolbar.Visible='off'; hold(TCA_Axes,'on'); %按钮 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - uibutton(TCA_Fig,'Text','导 入 图 片','BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],... 'FontWeight','bold','Position',[600,550,120,35],'FontSize',16,'ButtonPushedFcn',@loadPicFcn); uibutton(TCA_Fig,'Text','返 回 原 图','BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],... 'FontWeight','bold','Position',[600,505,120,35],'FontSize',16,'ButtonPushedFcn',@selectOriFcn); uibutton(TCA_Fig,'Text','清 空 数 据','BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],... 'FontWeight','bold','Position',[600,460,120,35],'FontSize',16,'ButtonPushedFcn',@clearAllFcn); uibutton(TCA_Fig,'Text','开 始 卷 积','BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],... 'FontWeight','bold','Position',[600,120,120,35],'FontSize',16,'ButtonPushedFcn',@convolutionFcn); uibutton(TCA_Fig,'Text','存 储 图 片','BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],... 'FontWeight','bold','Position',[600,75,120,35],'FontSize',16,'ButtonPushedFcn',@savePicFcn); uibutton(TCA_Fig,'Text','批 量 存 储','BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],... 'FontWeight','bold','Position',[600,30,120,35],'FontSize',16,'ButtonPushedFcn',@saveAllPicFcn); %按钮组- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - %TCA:Template Convolution App %BG :Button Group %SE :Smooth or Edge Detection
TCA_BG_SE=uibuttongroup(TCA_Fig); TCA_BG_SE.Title=' 模板卷积方法'; TCA_BG_SE.Position=[600,310,120,130]; TCA_BG_SE.FontSize=14; set(TCA_BG_SE,'SelectionChangedFcn',@changeCTFcn);%CT:Convolution Type % BG_S=uiradiobutton(TCA_BG_SE); BG_S.Position=[15 80 100 20]; BG_S.FontSize=13; BG_S.Text='模板运算'; BG_E=uiradiobutton(TCA_BG_SE); BG_E.Position=[15 55 100 20]; BG_E.FontSize=13; BG_E.Text='边缘检测'; BG_M=uiradiobutton(TCA_BG_SE); BG_M.Position=[15 30 100 20]; BG_M.FontSize=13; BG_M.Text='中值滤波'; BG_B=uiradiobutton(TCA_BG_SE); BG_B.Position=[15 5 100 20]; BG_B.FontSize=13; BG_B.Text='双边滤波'; %35:3*3 size Template or 5*5 size Template
TCA_BG_35=uibuttongroup(TCA_Fig); TCA_BG_35.Title=' 模板大小'; TCA_BG_35.Position=[600,210,120,80]; TCA_BG_35.FontSize=14; set(TCA_BG_35,'SelectionChangedFcn',@changeTSFcn);%CT:Template Size % BG_3=uiradiobutton(TCA_BG_35); BG_3.Position=[15 30 100 20]; BG_3.FontSize=13; BG_3.Text='3*3 模板'; BG_5=uiradiobutton(TCA_BG_35); BG_5.Position=[15 5 100 20]; BG_5.FontSize=13; BG_5.Text='5*5 模板'; %模板组- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - for i=1:5 for j=1:5 Hx(i,j)=uieditfield(TCA_Fig,'text'); set(Hx(i,j),'Position',[705+35*j,550-35*i,30,30]); set(Hx(i,j),'HorizontalAlignment','center'); set(Hx(i,j),'FontSize',15); set(Hx(i,j),'FontColor',[1 1 1].*0.45); Hy(i,j)=uieditfield(TCA_Fig,'text'); set(Hy(i,j),'Position',[705+35*j,330-35*i,30,30]); set(Hy(i,j),'HorizontalAlignment','center'); set(Hy(i,j),'FontSize',15); set(Hy(i,j),'FontColor',[1 1 1].*0.45); set(Hy(i,j),'Editable','off'); end
end % T33=[1 1;1 2;1 3;1 4;1 5;5 1;5 2;5 3;5 4;5 5;2 1;3 1;4 1;2 5;3 5;4 5]; for i=1:size(T33,1),set(Hx(T33(i,1),T33(i,2)),'Editable','off');end % uilabel(TCA_Fig,'Text','     横 向 梯 度 算 子 (Hx)','HorizontalAlignment','left',... 'BackgroundColor',[0 0.2510 0.4510],'FontColor',[1 1 1],'FontWeight','bold',... 'Position',[740,555,170,30],'FontSize',13); uilabel(TCA_Fig,'Text','     纵 向 梯 度 算 子 (Hy)','HorizontalAlignment','left',... 'BackgroundColor',[0 0.2510 0.4510],'FontColor',[1 1 1],'FontWeight','bold',... 'Position',[740,335,170,30],'FontSize',13); %- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - uilabel(TCA_Fig,'Text',' 卷 积 次 数','HorizontalAlignment','left',... 'BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],'FontWeight','bold',... 'Position',[600,170,80,30],'FontSize',13); TCA_Spinner=uispinner(TCA_Fig); TCA_Spinner.Position=[675,170,50,30]; TCA_Spinner.Limits = [0 200]; TCA_Spinner.RoundFractionalValues='on'; %- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - uilabel(TCA_Fig,'Text','  space-Sigma','HorizontalAlignment','left',... 'BackgroundColor',[0 0.2510 0.4510],'FontColor',[1 1 1],'FontWeight','bold',... 'Position',[740,110,120,30],'FontSize',13); uilabel(TCA_Fig,'Text','  range-Sigma','HorizontalAlignment','left',... 'BackgroundColor',[0 0.2510 0.4510],'FontColor',[1 1 1],'FontWeight','bold',... 'Position',[740,70,120,30],'FontSize',13); SigmaS=uieditfield(TCA_Fig,'text'); SigmaS.Position=[860,110,50,30]; SigmaS.HorizontalAlignment='center'; SigmaS.FontSize=15; SigmaS.Editable='off'; SigmaR=uieditfield(TCA_Fig,'text'); SigmaR.Position=[860,70,50,30]; SigmaR.HorizontalAlignment='center'; SigmaR.FontSize=15; SigmaR.Editable='off'; %噪声与通道- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Menu=uimenu(TCA_Fig); Menu.Text='图像预处理'; % Menu_1=uimenu(Menu); Menu_1.Text='色彩操作'; % Menu_1_1=uimenu(Menu_1); Menu_1_1.Text='灰度图像(多通道)'; set(Menu_1_1,'MenuSelectedFcn',@C_1_1_Fcn); Menu_1_2=uimenu(Menu_1); Menu_1_2.Text='灰度图像(单通道)'; set(Menu_1_2,'MenuSelectedFcn',@C_1_2_Fcn); Menu_1_3=uimenu(Menu_1); Menu_1_3.Text='二值图像(多通道)'; set(Menu_1_3,'MenuSelectedFcn',@C_1_3_Fcn); Menu_1_4=uimenu(Menu_1); Menu_1_4.Text='二值图像(单通道)'; set(Menu_1_4,'MenuSelectedFcn',@C_1_4_Fcn); Menu_1_5=uimenu(Menu_1); Menu_1_5.Text='反色图像'; set(Menu_1_5,'MenuSelectedFcn',@C_1_5_Fcn); % Menu_2=uimenu(Menu); Menu_2.Text='噪声操作'; % Menu_2_1=uimenu(Menu_2); Menu_2_1.Text='添加椒盐噪声'; set(Menu_2_1,'MenuSelectedFcn',@N_2_1_Fcn); Menu_2_2=uimenu(Menu_2); Menu_2_2.Text='添加高斯噪声'; set(Menu_2_2,'MenuSelectedFcn',@N_2_2_Fcn); Menu_2_3=uimenu(Menu_2); Menu_2_3.Text='添加泊松噪声'; set(Menu_2_3,'MenuSelectedFcn',@N_2_3_Fcn); %进程指示器- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - %'  正在进行第i次卷积 . . .' PGdisplay=uilabel(TCA_Fig,'Text','','HorizontalAlignment','left',... 'BackgroundColor',[0.8912 0.9727 1],'FontColor',[0.1 0.1 0.3],'FontWeight','bold',... 'Position',[740,30,170,35],'FontSize',15); %回调函数>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> %图片导入- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function loadPicFcn(~,~) try [filename,pathname]=uigetfile({'*.*';'*.jpg';'*.tif';'*.bmp';'*.png'},'File Selector'); oriPic=imread([pathname,filename]); subOriPic=oriPic; noisePic=[]; convolutionPic=[]; [XLim,YLim,~]=size(oriPic); hold(TCA_Axes,'off'); imshow(oriPic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); TCA_Axes.XLim=[0,max([XLim,YLim])]; TCA_Axes.YLim=[0,max([XLim,YLim])]; catch end
end %返回原图- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function selectOriFcn(~,~) subOriPic=oriPic;noisePic=[];convolutionPic=[]; hold(TCA_Axes,'off'); imshow(oriPic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); PGdisplay.Text=''; end %清空数据- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function clearAllFcn(~,~) oriPic=[]; subOriPic=[]; noisePic=[]; convolutionPic=[]; for ii=1:5 for jj=1:5 set(Hx(ii,jj),'Value',''); set(Hy(ii,jj),'Value',''); end
    end
    PGdisplay.Text=''; SigmaS.Value=''; SigmaR.Value=''; hold(TCA_Axes,'off'); imshow(uint8(ones(100,100).*255),'Parent',TCA_Axes); TCA_Axes.XLim=[0 100]; TCA_Axes.YLim=[0 100]; hold(TCA_Axes,'on'); end %开始卷积- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function convolutionFcn(~,~) selectedButton=TCA_BG_SE.SelectedObject; Type35=TCA_BG_35.SelectedObject.Text(1); if Type35=='3' Hx_=zeros(3,3);Hy_=zeros(3,3); for ii=1:3 for jj=1:3 if isempty(get(Hx(ii+1,jj+1),'Value')) Hx_(ii,jj)=0; else Hx_(ii,jj)=eval(get(Hx(ii+1,jj+1),'Value')); end
            end
        end for ii=1:3 for jj=1:3 if isempty(get(Hy(ii+1,jj+1),'Value')) Hy_(ii,jj)=0; else Hy_(ii,jj)=eval(get(Hy(ii+1,jj+1),'Value')); end
            end
        end else Hx_=zeros(5,5);Hy_=zeros(5,5); for ii=1:5 for jj=1:5 if isempty(get(Hx(ii,jj),'Value')) Hx_(ii,jj)=0; else Hx_(ii,jj)=eval(get(Hx(ii,jj),'Value')); end
            end
        end for ii=1:5 for jj=1:5 if isempty(get(Hy(ii,jj),'Value')) Hy_(ii,jj)=0; else Hy_(ii,jj)=eval(get(Hy(ii,jj),'Value')); end
            end
        end
    end if isempty(SigmaS.Value) SigmaS_=0; else SigmaS_=eval(SigmaS.Value); end if isempty(SigmaR.Value) SigmaR_=0; else SigmaR_=eval(SigmaR.Value); end if ~isempty(noisePic) convolutionPic=noisePic; else convolutionPic=subOriPic; end switch selectedButton.Text case '模板运算' for ii=1:TCA_Spinner.Value
                PGdisplay.Text=['  正在进行第',num2str(ii),'次卷积 . . .'];pause(0.05); convolutionPic=edgeDetectionOpt(convolutionPic); PicSet.(['p',num2str(ii)])=convolutionPic; end case '边缘检测' for ii=1:TCA_Spinner.Value
                PGdisplay.Text=['  正在进行第',num2str(ii),'次卷积 . . .'];pause(0.05); convolutionPic=edgeDetectionOpt(convolutionPic); PicSet.(['p',num2str(ii)])=convolutionPic; end case '中值滤波' for ii=1:TCA_Spinner.Value
                PGdisplay.Text=['  正在进行第',num2str(ii),'次卷积 . . .'];pause(0.05); convolutionPic=midFilterImage(convolutionPic); PicSet.(['p',num2str(ii)])=convolutionPic; end case '双边滤波' for ii=1:TCA_Spinner.Value
                PGdisplay.Text=['  正在进行第',num2str(ii),'次卷积 . . .'];pause(0.05); convolutionPic=bilateralFilterImage(convolutionPic); PicSet.(['p',num2str(ii)])=convolutionPic; end   
    end if TCA_Spinner.Value~=0 hold(TCA_Axes,'off'); imshow(convolutionPic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); end
    times=TCA_Spinner.Value; PGdisplay.Text='  卷积试验结束'; end %存储图像- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function savePicFcn(~,~) try [filename,pathname]=uiputfile({'*.*';'*.jpg';'*.tif';'*.bmp';'*.png'}); if ~isempty(convolutionPic) imwrite(convolutionPic,[pathname,filename]); elseif ~isempty(noisePic) imwrite(noisePic,[pathname,filename]); else imwrite(subOriPic,[pathname,filename]); end catch end
end %批量存储- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function saveAllPicFcn(~,~) if ~isempty(convolutionPic) try [filename,pathname]=uiputfile({'*.*';'*.jpg';'*.tif';'*.bmp';'*.png'}); num=abs(filename); pos=find(num==46); for ii=1:times imwrite(PicSet.(['p',num2str(ii)]),[pathname,filename(1:pos-1),'_',num2str(ii),filename(pos:end)]); end if ~isempty(noisePic) imwrite(noisePic,[pathname,filename(1:pos-1),'_0',filename(pos:end)]); else imwrite(subOriPic,[pathname,filename(1:pos-1),'_0',filename(pos:end)]); end catch end
    end
end %TCA_BG_SE- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function changeCTFcn(~,~) selectedButton=TCA_BG_SE.SelectedObject; switch selectedButton.Text case '模板运算',onOff1='on';onOff2='off';onOff3='off';clear1=0;clear2=1;clear3=1; case '边缘检测',onOff1='on';onOff2='on';onOff3='off';clear1=0;clear2=0;clear3=1; case '中值滤波',onOff1='on';onOff2='off';onOff3='off';clear1=0;clear2=1;clear3=1; case '双边滤波',onOff1='off';onOff2='off';onOff3='on';clear1=1;clear2=1;clear3=0; end for ii=1:5 for jj=1:5 set(Hx(ii,jj),'Editable',onOff1); set(Hy(ii,jj),'Editable',onOff2); if clear1==1,set(Hx(ii,jj),'Value','');end if clear2==1,set(Hy(ii,jj),'Value','');end
        end
    end
    SigmaS.Editable=onOff3; SigmaR.Editable=onOff3; if clear3==1,SigmaS.Value='';SigmaR.Value='';end changeTSFcn(); end %TCA_BG_35- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function changeTSFcn(~,~) selectedButton=TCA_BG_SE.SelectedObject; Type35=TCA_BG_35.SelectedObject.Text(1); switch selectedButton.Text case '模板运算',onOff1=1;onOff2=0; case '边缘检测',onOff1=1;onOff2=1; case '中值滤波',onOff1=1;onOff2=0; case '双边滤波',onOff1=0;onOff2=0; end if onOff1==1 if Type35=='3' for ii=1:size(T33,1) set(Hx(T33(ii,1),T33(ii,2)),'Editable','off'); set(Hx(T33(ii,1),T33(ii,2)),'Value','') end else for ii=1:size(T33,1) set(Hx(T33(ii,1),T33(ii,2)),'Editable','on'); end
        end
    end if onOff2==1 if Type35=='3' for ii=1:size(T33,1) set(Hy(T33(ii,1),T33(ii,2)),'Editable','off'); set(Hy(T33(ii,1),T33(ii,2)),'Value','') end else for ii=1:size(T33,1) set(Hy(T33(ii,1),T33(ii,2)),'Editable','on'); end
        end
    end
    PGdisplay.Text=''; end %噪声与通道- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function C_1_1_Fcn(~,~) NChannels=size(subOriPic,3); NChannel=uint8(zeros([size(oriPic,1),size(oriPic,2)])); for ii=1:NChannels
        NChannel=NChannel+subOriPic(:,:,ii)./NChannels; end for ii=1:3 subOriPic(:,:,ii)=NChannel; end hold(TCA_Axes,'off'); imshow(subOriPic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); noisePic=[]; convolutionPic=[]; PGdisplay.Text=''; end % function C_1_2_Fcn(~,~) NChannels=size(subOriPic,3); NChannel=uint8(zeros([size(oriPic,1),size(oriPic,2)])); for ii=1:NChannels
        NChannel=NChannel+subOriPic(:,:,ii)./NChannels; end
    subOriPic=NChannel; hold(TCA_Axes,'off'); imshow(subOriPic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); noisePic=[]; convolutionPic=[]; PGdisplay.Text=''; end % function C_1_3_Fcn(~,~) createDataExtractObj('二值化阈值','150',@C_1_3_Fcn_part2) end
function C_1_3_Fcn_part2(~,~) NChannels=size(subOriPic,3); NChannel=uint8(zeros([size(oriPic,1),size(oriPic,2)])); for ii=1:NChannels
        NChannel=NChannel+subOriPic(:,:,ii)./NChannels; end NChannel(NChannel>=coe_)=255; NChannel(NChannel<coe_)=0; for ii=1:3 subOriPic(:,:,ii)=NChannel; end hold(TCA_Axes,'off'); imshow(subOriPic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); noisePic=[]; convolutionPic=[]; PGdisplay.Text=''; end % function C_1_4_Fcn(~,~) createDataExtractObj('二值化阈值','150',@C_1_4_Fcn_part2) end
function C_1_4_Fcn_part2(~,~) NChannels=size(subOriPic,3); NChannel=uint8(zeros([size(oriPic,1),size(oriPic,2)])); for ii=1:NChannels
        NChannel=NChannel+subOriPic(:,:,ii)./NChannels; end NChannel(NChannel>=coe_)=255; NChannel(NChannel<coe_)=0; subOriPic=NChannel; hold(TCA_Axes,'off'); imshow(subOriPic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); noisePic=[]; convolutionPic=[]; PGdisplay.Text=''; end % function C_1_5_Fcn(~,~) NChannels=size(subOriPic,3); for ii=1:NChannels subOriPic(:,:,ii)=255-subOriPic(:,:,ii); end hold(TCA_Axes,'off'); imshow(subOriPic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); noisePic=[]; convolutionPic=[]; PGdisplay.Text=''; end % function N_2_1_Fcn(~,~) createDataExtractObj('椒盐噪声密度','0.1',@N_2_1_Fcn_part2) end
function N_2_1_Fcn_part2(~,~) noisePic=imnoise(subOriPic,'salt & pepper',coe_); hold(TCA_Axes,'off'); imshow(noisePic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); convolutionPic=[]; PGdisplay.Text=''; end % function N_2_2_Fcn(~,~) createDataExtractObj('高斯噪声标准差','0.1',@N_2_2_Fcn_part2) end
function N_2_2_Fcn_part2(~,~) noisePic=imnoise(subOriPic,'gaussian',0,coe_); hold(TCA_Axes,'off'); imshow(noisePic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); convolutionPic=[]; PGdisplay.Text=''; end % function N_2_3_Fcn(~,~) noisePic=imnoise(subOriPic,'poisson'); hold(TCA_Axes,'off'); imshow(noisePic,'Parent',TCA_Axes); hold(TCA_Axes,'on'); convolutionPic=[]; PGdisplay.Text=''; end %功能函数>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> %图像扩张- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function expandImage=expandImageFcn(oriImage) NChannels=size(oriImage,3); rows=size(oriImage,1); cols=size(oriImage,2); expandImage=uint8(zeros([rows+2,cols+2,NChannels])); for ii=1:NChannels
        tempMat=uint8(zeros([rows+2,cols+2])); tempMat(2:rows+1,2:cols+1)=oriImage(:,:,ii); tempMat(2:rows+1,1)=oriImage(:,1,ii); tempMat(2:rows+1,cols+2)=oriImage(:,cols,ii); tempMat(1,2:cols+1)=oriImage(1,:,ii); tempMat(rows+2,2:cols+1)=oriImage(rows,:,ii); tempMat(1,1)=oriImage(1,1,ii); tempMat(rows+2,1)=oriImage(rows,1,ii); tempMat(1,cols+2)=oriImage(1,cols,ii); tempMat(rows+2,cols+2)=oriImage(rows,cols,ii); expandImage(:,:,ii)=tempMat; end
end %梯度算子、拉普拉斯算子、模板运算- - - - - - - - - - - - - - - - - - - - - function newPic=edgeDetectionOpt(oriImage) NChannels=size(oriImage,3); rows=size(oriImage,1); cols=size(oriImage,2); sep=(size(Hx_,1)-1)/2; Wx_=sum(sum(Hx_)); Wy_=sum(sum(Hy_)); Wx_(Wx_==0)=1; Wy_(Wy_==0)=1; newPic=uint8(zeros([rows,cols,NChannels])); expandImage=expandImageFcn(oriImage); for ii=1:sep-1 expandImage=expandImageFcn(expandImage); end for ii=1+sep:rows+sep for jj=1+sep:cols+sep for k=1:NChannels
                tempMat=expandImage(ii-sep:ii+sep,jj-sep:jj+sep,k); Gx_=sum(sum(double(tempMat).*Hx_)); Gy_=sum(sum(double(tempMat).*Hy_)); Gxy_=sqrt(Gx_^2+Gy_^2)/Wx_/Wy_; newPic(ii-sep,jj-sep,k)=uint8(Gxy_); end
        end
    end
end %中值滤波- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function newPic=midFilterImage(oriImage) NChannels=size(oriImage,3); rows=size(oriImage,1); cols=size(oriImage,2); sep=(size(Hx_,1)-1)/2; Hx_(Hx_~=0)=1; Hx_(Hx_==0)=nan; newPic=uint8(zeros([rows,cols,NChannels])); expandImage=expandImageFcn(oriImage); for ii=1:sep-1 expandImage=expandImageFcn(expandImage); end for ii=1+sep:rows+sep for jj=1+sep:cols+sep for k=1:NChannels
                tempMat=expandImage(ii-sep:ii+sep,jj-sep:jj+sep,k); tempMat2=double(tempMat).*Hx_; newPic(ii-sep,jj-sep,k)=uint8(median(tempMat2(~isnan(tempMat2)))); end
        end
    end
end %双边滤波- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - function newPic=bilateralFilterImage(oriImage) NChannels=size(oriImage,3); rows=size(oriImage,1); cols=size(oriImage,2); sep=(size(Hx_,1)-1)/2; [xmesh,ymesh]=meshgrid(1:size(Hx_,1)); Ds2=(xmesh-sep-1).^2+(ymesh-sep-1).^2; newPic=uint8(zeros([rows,cols,NChannels])); expandImage=expandImageFcn(oriImage); for ii=1:sep-1 expandImage=expandImageFcn(expandImage); end for ii=1+sep:rows+sep for jj=1+sep:cols+sep for k=1:NChannels
                tempMat=expandImage(ii-sep:ii+sep,jj-sep:jj+sep,k); tempMat2=double(tempMat); Dr2=tempMat2-tempMat2(sep+1,sep+1); Dr2=Dr2.^2; Gs_=exp(-(Ds2)./(2*SigmaS_*SigmaS_)); Gr_=exp(-(Dr2)./(2*SigmaR_*SigmaR_)); Gsr_=Gs_.*Gr_; newPic(ii-sep,jj-sep,k)=uint8(sum(sum(tempMat2.*Gsr_))/sum(sum(Gsr_))); end
        end
    end
end %数据返回器>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> function createDataExtractObj(str,presetValue,Fcn) DE_Fig=uifigure('units','pixels'); DE_Fig.Position=[TCA_Fig.Position(1),TCA_Fig.Position(2)+TCA_Fig.Position(4)-100,200,100]; DE_Fig.NumberTitle='off'; DE_Fig.MenuBar='none'; DE_Fig.Name='参数设置'; DE_Fig.Color=[1,1,1].*0.98; uilabel(DE_Fig,'Text',[' ',str],'HorizontalAlignment','left',... 'BackgroundColor',[0 0.2510 0.4510],'FontColor',[1 1 1],'FontWeight','bold',... 'Position',[10,52,140,35],'FontSize',15); coe=uieditfield(DE_Fig,'text'); coe.Position=[130,52,60,35]; coe.HorizontalAlignment='center'; coe.FontSize=15; coe.Value=presetValue; uibutton(DE_Fig,'Text','确 认 设 置','BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],... 'FontWeight','bold','Position',[10,10,85,35],'FontSize',14,'ButtonPushedFcn',@ensureDEOFcn); uibutton(DE_Fig,'Text','取 消 设 置','BackgroundColor',[0.31 0.58 0.80],'FontColor',[1 1 1],... 'FontWeight','bold','Position',[105,10,85,35],'FontSize',14,'ButtonPushedFcn',@closeDEOFcn); function ensureDEOFcn(~,~) coe_=eval(coe.Value); Fcn(); delete(DE_Fig) end
    function closeDEOFcn(~,~) delete(DE_Fig) end
end
end